Inngangur að reikniritinu fyrir heildarbylgjuform (FWI)
1. Pendahuluan
Í nútíma jarðeðlisfræði er þörfin fyrir að „sjá“ uppbyggingu jarðar undir yfirborði jarðar að aukast. Orkuiðnaðurinn, aðgerðir til að draga úr náttúruhamförum , jarðhitaleit og jafnvel jarðskorpurannsóknir krefjast ítarlegra undir yfirborðslíkana: hvernig hraði jarðskjálftabylgna breytist með dýpi, hvar mörk berglaga eru staðsett og hvernig lítill fjölbreytileiki getur haft áhrif á útbreiðslu bylgna. Ein öflugasta aðferðin í þessu skyni er Full Waveform Inversion (FWI), snúningsreiknirit sem notar allar upplýsingar úr jarðskjálftabylgjuformum til að meta eðlisfræðilegar breytur undir yfirborði.
FWI er oft kallaður „gullstaðallinn“ í jarðskjálftaviðsnúningi þar sem hann getur framleitt neðanjarðarmyndir í hárri upplausn, sem er betri en hefðbundnar aðferðir sem nota aðeins ferðatíma eða einfalda sveifluvídd. Hins vegar kostar þetta afl sitt: FWI krefst nákvæmrar bylgjulíkanagerðar, umtalsverðra reikniauðlinda og vandlegra hagræðingaraðferða til að forðast að falla í rangar lausnir.
Þessi grein kynnir grunnhugtök FWI, helstu þætti reiknirits þess og almennar áskoranir og aðferðir við framkvæmd þess.
-
2. Hvað er full bylgjuformsbreyting?
Einfaldlega sagt er FWI ferlið við að finna jarðskjálftalíkan sem passar best við hermdar jarðskjálftagögn og mæld jarðskjálftagögn. „Besta aðlögun“ er skilgreind með hlutlægu falli (misræmi), svo sem ferningmismuninum á milli mældra og tilbúins gagna á hverjum tímapunkti og við hvern móttakara.
Grundvallarmunurinn á FWI og öðrum aðferðum við andhverfu er:
– FWI notar alla bylgjuformið (fasa, sveifluvídd, truflun, margfeldi, dreifing), ekki bara komutíma.
– FWI byggir á því að leysa bylgjujöfnuna (hljóðbylgju-, teygju- eða anelastíska bylgjujöfnu) til að herma eftir tilbúnum gögnum.
– FWI er stórfellt ólínulegt bestunarvandamál, því líkanbreyturnar (t.d. P- eða S-bylgjuhraði) geta talið milljónir frumna á 2D/3D risti.
-
3. Helstu þættir í FWI
3.1 Athugunargögn og tilbúin gögn
FWI krefst:
– Athugunargögn: jarðskjálftamælingar á vettvangi (skotsöfnun) frá ýmsum aðilum og móttökutækjum.
– Tilbúin gögn: niðurstöður tölulegrar hermunar á bylgjuútbreiðslu á tímabundnu líkani.
Tilbúin gögn eru reiknuð með því að leysa bylgjujöfnuna (t.d. hljóðbylgjujöfnuna):
\[
\frac{1}{v^2(\mathbf{x})}\frac{\partial^2 p}{\partial t^2} – \nabla^2 p = s(\mathbf{x},t)
\]
þar sem \(v(\mathbf{x}) \) er bylgjuhraðinn, \(p \) er þrýstingurinn/bylgjan á skalarstigi og \(s \) er upptök bylgjunnar.
3.2 Upphafleg fyrirmynd
FWI er mjög viðkvæmt fyrir upphafslíkaninu. Ef upphafslíkanið er of langt frá raunverulegum aðstæðum getur reikniritið orðið fyrir hringrásarhoppi, sem er þegar fasamunurinn á milli tilbúnu og mældra gagna er meira en hálfur tímabil, sem veldur því að hallinn ýtir lausninni í ranga átt.
Upphafslíkanið er venjulega fengið úr:
– umsnúningur ferðatíma (sneiðmyndataka),
– stórhraðalíkön úr jarðfræðilegri greiningu,
– eða fjölþátta nálgun (byrjað á lágum tíðnum).
3.3 Markfall (mispassun)
Algengar markmiðsföll:
\[
J(m) = \frac{1}{2}\sum_{s}\sum_{r}\int (d_{\text{syn}}(t; m) - d_{\text{obs}}(t))^2 \, dt
\]
þar sem \(m \) er líkanbreyta (t.d. hraði), \(s \) er upptökuvísitala og \(r \) er móttökuvísitala.
Auk klassíska L2 misræmisins eru einnig til aðrar misræmi til að draga úr hoppi í hringrás, til dæmis:
– umslagsbundin misræmi,
– fasa-eingöngu misræmi,
- ákjósanlegur flutningur (Wasserstein),
– eða samsvarandi síur.
3.4 Útreikningur á halla: Aðlögunarástandsaðferð
Einn af lykilþáttum FWI er hvernig á að reikna út stigul á skilvirkan hátt. Þar sem fjöldi líkanbreytna er svo mikill er bein töluleg útleiðsla ómöguleg. Lausnin er aðlögunarástandsaðferðin.
Innsæismyndin:
1. Framvirk líkangerð: reiknaðu framvirka bylgjusviðið út frá upptökum í núverandi líkani.
2. Reiknið út gagnaleifina (Δd = d_{syn} – d_{obs}).
3. Aðliggjandi líkanagerð: sprautaðu leifinni inn sem „endurkomugjafa“ frá móttökustaðnum og dreift henni aftur í tímann.
4. Fylgni milli frambylgjusviðsins og aðliggjandi bylgjusviðsins framkallar halla á líkanbreyturnar.
Með þessari tækni er reiknikostnaður hallans nokkurn veginn jafngildur tvöföldum bylgjulíkönum á hverja uppsprettu (framvirk + aðliggjandi), þannig að hann er samt mikill en framkvæmanlegur á HPC/GPU.
3.5 Uppfærslukerfi líkans (hagræðing)
Þegar hallatölunni \( \nabla J \) hefur verið náð er líkanið uppfært með bestunaraðferð, til dæmis:
– Brattasta niðurleiðin (einfaldasta),
– Samtengdur halli,
– L-BFGS (algengt notað vegna þess að það er skilvirkt við stór vandamál),
– eða Newton/Quasi-Newton aðferðin.
Grunnuppfærslur:
\[
m_{k+1} = m_k – \alpha_k \, H_k^{-1}\nabla J(m_k)
\]
þar sem \( \alpha_k \) er skreflengdin og \(H_k^{-1} \) er öfug hessísk nálgun (t.d. í L-BFGS).
-
4. FWI reiknirit vinnuflæði (hnykktur)
Almennt séð keyrir FWI endurtekið:
1. Veldu upphafslíkanið \(m_0 \).
2. Fyrir hverja uppsprettu:
– gera framvirka líkön → tilbúin gögn,
– reikna út leifarnar miðað við athugunargögnin,
– gera aðlögunarlíkön,
– uppsöfnun halla.
3. Beita skal forstillingu (t.d. birtujöfnun eða sléttun).
4. Uppfærðu líkanið með bestunaraðferðum.
5. Endurtakið þar til samleitni eða ítrunarmörkum er náð.
Venjulega er FWI keyrt á mörgum skala, byrjað á lágum tíðnum (leiðrétting á íhlutum stórra líkana) og síðan fært upp í hærri tíðni (bætið við smáatriðum).
-
5. Helstu áskoranir í erlendum fjármálum
5.1 Hringrásarslepp
Þetta er þekktasta vandamálið. Þegar tilbúnar og mældar upplýsingar eru ekki „í fasa“ getur misræmi í L2 leitt til bestunarinnar sem leiðir til falskra staðbundinna lámarka. Almenna lausnin:
– byrjar á mjög lágum tíðnum,
– bæta upphaflega líkanið (sneiðmyndatöku),
– að nota aðra mislíka,
– beita tímaglugga og aðferðum til að velja gögn.
5.2 Reiknikostnaður
Þrívíddar FWI með mörgum uppsprettum getur krafist þúsunda til milljóna bylgjuformshermuna. Þetta krefst:
– samsíða útreikningar (klasi, GPU),
– minnissparnaður (eftirlitspunktur til að vista bylgjusvið),
– sem og aðferðir til að velja undirmengi heimilda (heimildarkóðun eða smáhópavinnsla ala vélanám).
5.3 Ósamræmi í eðlisfræði og hávaða
Gögn úr reitnum innihalda hávaða, áhrif mælitækja, ósamhverfu, deyfingu (Q), flókna landslagsfræði og óvissu í upptökum. Ef framvirk líkanagerð er of einföld (t.d. hljóðfræðileg þegar miðillinn er teygjanlegur) geta niðurstöður umsnúningsins verið skekktar.
5.4 Líkanstillingar
Að velja hvaða breytur á að snúa við (t.d. \(v_p \), eðlisþyngd, \(v_s \), anisótrópía, Q) hefur áhrif á næmi og stöðugleika. Of margir breytur geta valdið málamiðlunum og versnað samleitni.
-
6. Umsókn um FWI
FWI er mikið notað fyrir:
– olíu- og gasleit: bæta upplausn hraðalíkana, aðstoða við flutning og túlkun mannvirkja,
– jarðvarmi: kortlagning sprungusvæða og jarðfræðilegra breytinga,
– hnattræn/svæðisbundin jarðskjálftafræði: möttul- og jarðskorpumyndun,
– jarðtækni og verkfræði nálægt yfirborði: grunn kortlagning fyrir grunna, jarðgöng og mótvægisaðgerðir.
Kostir FWI eru augljósir þegar miðillinn er flókinn og gögnin rík af bylgjuupplýsingum: litlar bylgjubrot, fjölslóðir og margfeldi sem venjulega eru talin truflun geta í raun verið upplýsingauppsprettur.
-
7. Niðurstaða
Full Waveform Inversion er jarðskjálftaumsnúningsreiknirit sem notar ríkulegar upplýsingar í bylgjuformum til að smíða neðanjarðarlíkön með mikilli upplausn. Lykillinn að árangri þess liggur í nákvæmri bylgjujöfnulíkönum, skilvirkum hallaútreikningum með aðlögunarástandsaðferðinni og hagræðingaraðferðum sem forðast staðbundin lágmörk eins og hopp í hringrás. Þrátt fyrir miklar reiknikröfur og vandaða vinnuflæðishönnun hefur FWI reynst mjög öflug aðferð í jarðmyndun og auðlindakönnun.
Ef þú vilt get ég haldið áfram með tæknilegri framhaldsgreinar — til dæmis umræðu um aðlögunarhalla með formlegri hætti, dæmi um FWI gervikóða eða bandpass- og gluggatengdar fjölkvarðaaðferðir fyrir raunveruleg gögn.
