Stokastisk poromekanisk analys förutspår en betydande överskridande sannolikhet för jordbävningen i Pohang, Sydkorea 2017, Mw 5.5

Varför en människoskapad jordbävning spelar roll

År 2017 skakade en magnitud 5,5-jordbävning staden Pohang i Sydkorea, orsakade byggnadsskador och överraskade forskare eftersom händelsen kopplades till ett geotermiskt projekt snarare än en naturlig förskjutning längs en ytförbindelse. Att förstå hur mänskliga aktiviteter som djup vätskeinjektion kan utlösa så kraftiga skalv är avgörande om vi vill kunna utöka lågkolenergiproduktion utan att utsätta närboende samhällen för risk. I denna studie används ett fysikbaserat, sannolikhetsdrivet angreppssätt för att ställa en enkel men långtgående fråga: givet vad vi (inte) vet om berggrunden och spänningarna under marken, hur sannolikt var det att en jordbävning av Pohangs storlek skulle inträffa?

Hur energiprojekt kan väcka dolda förkastningar

Förbättrade geotermiska system utvinner värme genom att injicera högtrycksvatten djupt under marken för att öppna befintliga sprickor och förbättra vattencirkulationen. I Pohang skedde injektionen i granitisk berggrund på ungefär 4,2 kilometers djup, nära en redan existerande förkastning som inte bröt igenom till ytan. När trycksatt vatten tränger in i bergarterna höjs trycket i de små porerna mellan mineralpartiklarna och bergmassans sätt att bära upp spänningar ändras något. Dessa subtila förändringar kan minska friktionen på en förkastning och låta den glida. I Pohang visar flera bevislinjer att huvudskalvet inträffade nära injektionsbrunnen, längs en mogen förkastningsyta vars exakta orientering och spänningstillstånd var, och fortfarande är, dåligt begränsade.

Att omvandla osäkerhet till en sannolikhetsprognos

De flesta tidigare analyser av Pohang försökte rekonstruera en enda ”bästa” modell av underjorden, och antog att förkastningen redan var extremt nära brott så att även små förändringar i spänning kunde utlösa ett skalv. Fältmätningar och seismiska data pekar emellertid på att denna förkastning var mer stabil än den enkla bilden tillåter. Istället för att satsa på ett scenario behandlar författarna centrala underjordiska egenskaper — såsom riktning och storlek på spänningarna, förkastningens lutning och friktionen på dess yta — som osäkra storheter. Därefter använder de en teknik kallad Monte Carlo-simulering: tusentals något olika, men fysiskt plausibla, underjordiska världar genereras och för varje sådan räknar de ut hur vätsketrycket sprider sig, hur bergarterna reagerar mekaniskt och om förkastningen skulle börja glida och i så fall hur mycket.

Simulera hur förkastningar svarar på injektion

För att hålla det stora antalet försök beräkningsmässigt hanterbart använder teamet analytiska formler istället för tunga numeriska modeller för att beskriva hur injektionen höjer portrycket kring brunnen och hur den förändringen påverkar det omgivande spänningsfältet. De undersöker två realistiska sätt som förkastningen kan röra sig på, båda involverande en sned blandning av vertikalt och sidledes glid. I deras baskasus, med medelvärden för berg- och spänningsegenskaper, förblir förkastningen faktiskt stabil trots injektionen — vilket tydligt står i kontrast till den verkliga jordbävningen. När de låter de osäkra parametrarna variera inom intervall som stöds av mätningar och laboratorietester producerar vissa realisationer bara obetydliga, oupptäckbara skalv, medan andra genererar mycket större händelser. Genom att omvandla storleken på det glidande förkastningsområdet i varje realisation till en jordbävningsmagnitud bygger de upp en fullständig sannolikhetsfördelning över möjliga utfall.

Hur sannolik var Pohang‑jordbävningen?

Simuleringarna visar att, under förhållanden relevanta för Pohang, skulle den största möjliga inducerade jordbävningen i princip kunna närma sig magnitud 7, men sådana händelser är mycket osannolika. Mer talande är den skattade chansen att överskrida vissa magnituder. För skalv lika stora som den faktiska 2017-händelsen (Mw 5,5) förutspår modellen en överskridandesannolikhet på mellan ungefär 7 % och 15 %, beroende på vilken glidrörelsemönster som antas. Detta intervall överensstämmer nära med sannolikheten som oberoende härleds från den observerade sekvensen av mindre skalv vid platsen. Analysen visar också en tydlig koppling mellan hur nära en förkastning ligger att kollapsa innan injektionen och storleken på det följande skalvet. I Pohang, när det ursprungliga ”säkerhetsmarginalet” på förkastningen sjunker under ungefär 0,1–0,2 megapascal, kan även måttliga poromekaniska störningar få den att brista i ett skadligt brott.

Vad detta betyder för framtida geoenergiprojekt

För en lekman är huvudbudskapet att Pohang‑jordbävningen inte var en slumpmässig olycka, inte heller en oundviklig följd av geotermisk utveckling, utan en kvantifierbar risk som beror på hur kritiskt belastade närliggande förkastningar är och hur mycket vi vet om dem. Denna studie visar att genom att kombinera fysikbaserade modeller med systematisk osäkerhetsanalys är det möjligt att i förväg uppskatta sannolikheten att injektion utlöser skalv av en viss storlek. Den varnar för att förkastningar som redan ligger nära brott kan ge upphov till skadliga jordbävningar från relativt små tryckförändringar, och antyder att traditionella "trafikljussystem" baserade enbart på övervakning av små händelser kanske inte räcker. I stället kommer noggrann platskarakterisering och adaptiva, modellinformerade riskbedömningar — av den typ som demonstreras här — att vara avgörande om vi ska nyttja djupa underjordiska resurser säkert och ansvarsfullt.

Citering: 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

Nyckelord: inducerad seismisk aktivitet, geotermisk energi, sprickstabilitet, vätskeinjektion, jordbävningsrisk