Stochastyczna analiza poromechaniczna prognozuje istotne prawdopodobieństwo przekroczenia dla trzęsienia ziemi Mw 5.5 w Pohang (Korea Południowa) z 2017 roku

Dlaczego trzęsienie wywołane działalnością człowieka ma znaczenie

W 2017 roku trzęsienie o magnitudzie 5,5 wstrząsnęło miastem Pohang w Korei Południowej, uszkadzając budynki i zaskakując naukowców tym, że zostało powiązane z projektem geotermalnym, a nie z naturalnym przesunięciem uskoku. Zrozumienie, w jaki sposób działania ludzkie, takie jak głęboki wtrysk płynów, mogą wywołać tak silne wstrząsy, jest kluczowe, jeśli chcemy rozwijać niskoemisyjne źródła energii bez narażania pobliskich społeczności. Niniejsze badanie wykorzystuje podejście oparte na fizyce i prawdopodobieństwie, aby postawić proste pytanie o duże konsekwencje: biorąc pod uwagę to, co (nie) wiemy o skałach i naprężeniach pod ziemią, jak prawdopodobne było wystąpienie trzęsienia o wielkości takiej jak w Pohang?

Jak projekty energetyczne mogą obudzić ukryte uskoki

Zaawansowane systemy geotermalne pozyskują ciepło przez wtrysk wysokociśnieniowej wody na dużą głębokość, aby otworzyć istniejące spękania i poprawić cyrkulację wody. W Pohang wtrysk miał miejsce do granitoidalnych skał na głębokości około 4,2 km, w pobliżu wcześniej istniejącego uskoku, który nie sięgał powierzchni. Gdy w skały wnika woda pod ciśnieniem, podnosi ciśnienie w drobnych porach między ziarnami minerałów i nieznacznie zmienia sposób, w jaki masa skalna przenosi naprężenia. Te subtelne zmiany mogą zmniejszyć tarcie na uskokach, pozwalając im się przesunąć. W Pohang kilka linii dowodów wskazuje, że główne trzęsienie zdarzyło się blisko studni wtryskowej, wzdłuż dojrzałej płaszczyzny uskoku, której dokładna orientacja i stan naprężeń były i nadal są słabo określone.

Przekształcanie niepewności w prognozę prawdopodobieństwa

Większość wcześniejszych analiz trzęsienia w Pohang próbowała odtworzyć pojedynczy „najlepszy” model podziemny, zakładając, że uskok był już ekstremalnie bliski awarii, tak że nawet drobne zmiany naprężeń mogły go wywołać. Jednak pomiary terenowe i dane sejsmiczne sugerują, że ten uskok był bardziej stabilny niż pozwala na to taki prosty obraz. Zamiast polegać na jednym scenariuszu, autorzy traktują kluczowe właściwości podziemne — takie jak kierunek i wielkość naprężeń, kąt uskoku oraz tarcie na powierzchni uskoku — jako wielkości niepewne. Następnie stosują technikę zwaną symulacją Monte Carlo: generują tysiące nieco różnych, lecz fizycznie wiarygodnych wersji podziemnego środowiska i dla każdej obliczają, jak rozprzestrzenia się ciśnienie płynu, jak skały reagują mechanicznie oraz czy uskok by się przesunął, i w jakim stopniu.

Symulowanie reakcji uskoków na wtrysk

Aby utrzymać dużą liczbę prób w granicach wykonalności obliczeniowej, zespół używa formuł analitycznych zamiast ciężkich modeli numerycznych do opisu, jak wtrysk podnosi ciśnienie porowe wokół studni i jak ta zmiana wpływa na otaczające pole naprężeń. Badają dwa realistyczne sposoby, w jakie uskok mógłby się poruszyć, oba obejmujące skośne połączenie ruchu pionowego i bocznego. W wersji podstawowej, z przeciętnymi właściwościami skał i naprężeń, uskok pozostaje stabilny mimo wtrysku — co wyraźnie stoi w sprzeczności z rzeczywistym trzęsieniem. Gdy pozwalają, by parametry niepewne zmieniały się w ramach zakresów popartych pomiarami i testami laboratoryjnymi, niektóre realizacje dają tylko maleńkie, niewykrywalne wstrząsy, podczas gdy inne generują znacznie większe zdarzenia. Przekształcając wielkość obszaru przemieszczenia uskoku w każdej realizacji na magnitudę trzęsienia, budują pełny rozkład prawdopodobieństwa możliwych rezultatów.

Jak prawdopodobne było trzęsienie w Pohang?

Symulacje pokazują, że w warunkach istotnych dla Pohang największe możliwe trzęsienie indukowane teoretycznie mogłoby zbliżyć się do magnitudy 7, lecz takie zdarzenia są bardzo mało prawdopodobne. Znacznie bardziej informujący jest szacowany odsetek przekroczeń określonych magnitud. Dla trzęsień tak silnych jak rzeczywiste zdarzenie z 2017 roku (Mw 5,5) model prognozuje prawdopodobieństwo przekroczenia w przybliżeniu między około 7% a 15%, w zależności od przyjętego wzoru poślizgu. Ten zakres dobrze zgadza się z prawdopodobieństwem wywnioskowanym niezależnie z zaobserwowanej sekwencji mniejszych wstrząsów na miejscu. Analiza ujawnia także wyraźny związek między tym, jak blisko awarii znajdował się uskok przed wtryskiem, a wielkością następującego trzęsienia. W Pohang, gdy początkowy „margines bezpieczeństwa” na uskok spada poniżej około 0,1–0,2 megapascalów, nawet umiarkowane poromechaniczne zaburzenia mogą doprowadzić do uszkadzającego pęknięcia.

Co to oznacza dla przyszłych projektów geoenergetycznych

Dla laika kluczowy wniosek jest taki, że trzęsienie w Pohang nie było anomalią ani nieuniknionym skutkiem rozwoju geotermii, lecz ryzykiem możliwym do wyliczenia, które zależy od tego, jak krytycznie naprężone są pobliskie uskoki i ile o nich wiemy. Badanie pokazuje, że łącząc modele oparte na fizyce z systematyczną analizą niepewności, można z wyprzedzeniem oszacować prawdopodobieństwo, że wtrysk wywoła trzęsienia o danej wielkości. Ostrzega, że uskoki już bliskie awarii mogą wygenerować niszczące trzęsienia przy stosunkowo niewielkich zmianach ciśnienia, i sugeruje, iż tradycyjne systemy „światła drogowego” oparte wyłącznie na monitorowaniu małych zdarzeń mogą być niewystarczające. Zamiast tego konieczne będą staranna charakterystyka miejsc oraz adaptacyjne, oparte na modelach oceny ryzyka — takie jak zaprezentowane w tym badaniu — jeśli chcemy bezpiecznie i odpowiedzialnie wykorzystywać zasoby głębokiego podłoża.

Cytowanie: Wu, H., Vilarrasa, V., Parisio, F. et al. Stochastic poromechanical analysis forecasts a notable exceedance probability for the 2017 Pohang, South Korea, M_w 5.5 earthquake. Commun Earth Environ 7, 236 (2026). https://doi.org/10.1038/s43247-026-03268-7

Słowa kluczowe: sejsmiczność indukowana, energia geotermalna, stabilność uskoku, wtrysk płynów, ryzyko trzęsień ziemi