Ynlieding ta it Full Waveform Inversion (FWI) algoritme
1. Pendahuluan
Yn 'e moderne geofysika nimt de needsaak om de ûndergrûnske struktuer fan 'e ierde te "sjen" ta. De enerzjysektor, rampbestriding , geotermyske eksploraasje en sels tektonysk ûndersyk fereaskje detaillearre ûndergrûnske modellen: hoe't de snelheid fan seismyske weagen feroaret mei de djipte, wêr't de grinzen fan 'e rotslaach lizze, en hoe't lytse heterogeniteiten ynfloed kinne hawwe op weachfersprieding. Ien fan 'e machtichste oanpakken foar dit doel is Full Waveform Inversion (FWI), in ynverzje-algoritme dat folsleine ynformaasje fan seismyske golffoarmen brûkt om fysike parameters fan 'e ûndergrûn te skatten.
FWI wurdt faak oantsjutten as de "gouden standert" fan seismyske ynverzje, om't it ûndergrûnske ôfbyldings mei hege resolúsje kin produsearje, en dêrmei konvinsjonele oanpakken dy't allinich reistiid of ienfâldige amplitude brûke, oertreffe. Dizze krêft komt lykwols tsjin in priis: FWI fereasket krekte weachmodellering, substansjele kompjûterboarnen en soarchfâldige optimalisaasjestrategyen om te foarkommen dat men yn ferkearde oplossingen telâne komt.
Dit artikel yntrodusearret de basisbegripen fan FWI, de wichtichste komponinten fan syn algoritme, en algemiene útdagings en strategyen yn 'e ymplemintaasje dêrfan.
-
2. Wat is folsleine golffoarmynverzje?
Simpelwei sein, FWI is it proses fan it finen fan in ûndergrûnsmodel dat simulearre seismyske gegevens it bêste oerienkomt mei waarnommen seismyske gegevens. "Bêste fit" wurdt definiearre troch in objektive funksje (misfit), lykas it kwadraatferskil tusken waarnommen en synthetyske gegevens op elk tiidpunt en by elke ûntfanger.
De fûnemintele ferskillen tusken FWI en oare ynverzjemetoaden binne:
– FWI brûkt de hiele golffoarm (faze, amplitude, ynterferinsje, mearfâld, diffraksje), net allinich it selektearjen fan oankomsttiid.
– FWI is ôfhinklik fan it oplossen fan 'e golffergeliking (akoestysk, elastysk of anelastysk) om synthetyske gegevens te simulearjen.
– FWI is in grutskalich net-lineêr optimalisaasjeprobleem, om't de modelparameters (bygelyks P- of S-weachsnelheden) miljoenen sellen op in 2D/3D-raster telle kinne.
-
3. Haadkomponinten yn FWI
3.1 Waarnimmingsgegevens en synthetyske gegevens
FWI fereasket:
– Waarnimmingsgegevens: seismyske fjildopnamen (shot collect) fan ferskate boarnen en ûntfangers.
– Syntetyske gegevens: resultaten fan numerike simulaasje fan weachfoarútgong op in tydlik model.
Syntetyske gegevens wurde berekkene troch de weachfergeliking op te lossen (bygelyks de akoestyske weachfergeliking):
\[
\frac{1}{v^2(\mathbf{x})}\frac{\partial^2 p}{\partial t^2} – \nabla^2 p = s(\mathbf{x},t)
\]
wêrby't \(v(\mathbf{x}) \) de golfsnelheid is, \(p \) de skalêre druk/golf is, en \(s \) de boarne is.
3.2 Inisjele model
FWI is tige gefoelich foar it earste model. As it earste model te fier fan 'e werklike omstannichheden ôf is, kin it algoritme syklusferskowing ûnderfine, wat is as it fazeferskil tusken de synthetyske en waarnommen gegevens mear as in heale perioade is, wêrtroch't de gradiënt de oplossing yn 'e ferkearde rjochting triuwt.
It earste model wurdt meastentiids krigen fan:
– reistiidynverzje (tomografy),
– makrosnelheidsmodellen út geologyske analyze,
– of in mearskalige oanpak (begjinnende fan lege frekwinsjes).
3.3 Objektyffunksje (Misfit)
Algemiene doelfunksjes:
\[
J(m) = \frac{1}{2}\sum_{s}\sum_{r}\int (d_{syn}}(t; m) - d_{obs}}(t))^2 \, dt
\]
wêrby't \(m \) in modelparameter is (bygelyks snelheid), \(s \) de boarne-yndeks is, en \(r \) de ûntfanger-yndeks is.
Neist de klassike L2-misfit binne der ek alternative misfits om syklusferskowing te ferminderjen, bygelyks:
– envelop-basearre misfit,
– allinnich-faze misfit,
- optimaal ferfier (Wasserstein),
– of oerienkommende filters.
3.4 Gradiëntberekkening: Adjoint-State-metoade
Ien fan 'e wichtichste skaaimerken fan FWI is hoe't gradiënten effisjint berekkenje kinne wurde. Omdat it oantal modelparameters sa grut is, is direkte numerike ôflieding ûnmooglik. De oplossing is de adjoint-state-metoade.
It yntuïtive byld:
1. Foarútmodellering: berekkenje it foarútgolffjild fan 'e boarne yn it hjoeddeiske model.
2. Bereken it dataresidu (Δd = d_syn – d_obs).
3. Adjoint-modellering: ynjeksjearje it residu as in "weromkearboarne" fan 'e ûntfangerposysje en ferspriede it werom yn 'e tiid.
4. De korrelaasje tusken de foarútgeande en oanswettende weachfjilden produseart in gradiënt foar de modelparameters.
Mei dizze technyk binne de berekkeningskosten fan 'e gradiënt sawat lykweardich oan twa kear de wellemodellering per boarne (foarút + oanswettend), dus it is noch altyd grut, mar útfierber op HPC/GPU.
3.5 Modelupdateskema (Optimalisaasje)
Sadree't de gradiënt \( \nabla J \) krigen is, wurdt it model bywurke mei in optimalisaasjemetoade, bygelyks:
– Steilste ôfdaling (ienfâldichste),
– Konjugearre gradiënt,
– L-BFGS (faak brûkt om't it effisjint is foar grutte problemen),
– of Newton/Quasi-Newton-metoade.
Basis updates:
\[
m_{k+1} = m_k – \α_k \, H_k^{-1} \nabla J(m_k)
\]
wêrby't \( \alpha_k \) de staplingte is, en \(H_k^{-1} \) de inverse Hessian-benadering is (bygelyks yn L-BFGS).
-
4. FWI-algoritme-workflow (bondige)
Yn 't algemien rint FWI iteratyf:
1. Selektearje it earste model \(m_0 \).
2. Foar elke boarne:
– foarútmodellering dwaan → synthetyske gegevens,
– berekkenje de residuen tsjin de observaasjegegevens,
- oanslutende modellering dwaan,
- gradiëntakkumulaasje.
3. Tapasse foarkondisjonearring (bygelyks ljochtkompensaasje of glêdmeitsjen).
4. Aktualisearje it model mei optimalisaasjemetoaden.
5. Werhelje oant de konverginsje of de iteraasjelimyt berikt is.
Typysk wurdt FWI op meardere skaal útfierd, begjinnend by lege frekwinsjes (korrizjearjen fan grutskalige modelkomponinten) en dan omheech gien nei hegere frekwinsjes (tafoegjen fan detail).
-
5. Wichtichste útdagings yn FWI
5.1 Syklus oerslaan
Dit is it bekendste probleem. As syntetyske en waarnommen gegevens net "yn faze" binne, kin L2-misfit liede ta optimalisaasje ta falske lokale minima. De algemiene oplossing:
- begjinnend fan tige lege frekwinsjes,
– it earste model (tomografy) ferbetterje,
- it brûken fan alternative misfits,
– tiidfinster- en gegevensseleksjestrategyen tapasse.
5.2 Berekkeningskosten
3D FWI mei meardere boarnen kin tûzenen oant miljoenen golffoarmsimulaasjes fereaskje. Dit fereasket:
– parallelle kompjûtergebrûk (kluster, GPU),
– ûnthâldbesparring (kontrôlepunten om golffjilden te bewarjen),
– lykas ek strategyen foar it selektearjen fan boarnesubsets (boarnekodearring of mini-batch ala masinelearen).
5.3 Mismatch tusken natuerkunde en lûd
Fjildgegevens befetsje rûs, ynstruminteffekten, anisotropie, ferswakking (Q), komplekse topografy en boarne-ûnwissichheid. As de foarútmodellering te ienfâldich is (bygelyks akoestysk as it medium elastysk is), kinne de ynverzjeresultaten fertekene wêze.
5.4 Modelparameterisaasje
It kiezen fan hokker parameters om te ynvertearjen (bygelyks, \(v_p \), tichtens, \(v_s \), anisotropie, Q) beynfloedet gefoelichheid en stabiliteit. Tefolle parameters kinne ôfwagings feroarsaakje en de konverginsje fergrutsje.
-
6. FWI-oanfraach
FWI wurdt breed brûkt foar:
– oalje- en gaseksploraasje: ferbetterjen fan de resolúsje fan snelheidsmodellen, helpe by migraasje en ynterpretaasje fan struktueren,
– geotermysk: it yn kaart bringen fan breuksônes en litologyske feroarings,
– wrâldwide/regionale seismology: ôfbylding fan mantel en korst,
– geotechnyske en tichtby-oerflak technyk: ûndjippe kartering foar fundearrings, tunnels en mitigaasje.
De foardielen fan FWI binne dúdlik as it medium kompleks is en de gegevens ryk binne oan weachynformaasje: lytse diffraksjes, multipaden en mearfâlden dy't normaal as ynterferinsje beskôge wurde, kinne eins boarnen fan ynformaasje wêze.
-
7. Konklúzje
Folsleine golffoarmynverzje is in seismysk ynverzje-algoritme dat de rike ynformaasje yn golffoarmen brûkt om ûndergrûnmodellen mei hege resolúsje te bouwen. De kaai ta sukses leit yn krekte golffergelikingmodellering, effisjinte gradiëntberekkeningen mei de adjoint-state-metoade, en optimalisaasjestrategyen dy't lokale minima lykas syklusferskowing foarkomme. Nettsjinsteande de hege berekkeningseasken en soarchfâldige workflowûntwerp hat FWI bewiisd in heul krêftige oanpak te wêzen yn ierdôfbylding en boarnenûndersyk.
As jo wolle, kin ik trochgean mei mear technyske ferfolchartikels - bygelyks it formeler besprekken fan 'e adjunkte gradiëntôflieding, foarbylden fan FWI-pseudokode, of band-pass- en windowing-basearre multiskaalstrategyen foar echte gegevens.
