Chiyambi cha Algorithm ya Full Waveform Inversion (FWI)
1. Pendauluan
Mu geophysics yamakono, kufunika "kowona" kapangidwe ka pansi pa Dziko Lapansi kukuchulukirachulukira. Makampani opanga mphamvu, kuchepetsa masoka, kufufuza za geothermal, komanso kafukufuku wa tectonic amafuna zitsanzo zatsatanetsatane za pansi pa nthaka: momwe liwiro la mafunde a seismic limasinthira ndi kuya, komwe malire a miyala ali, komanso momwe kusiyana kochepa kungakhudzire kufalikira kwa mafunde. Njira imodzi yamphamvu kwambiri pachifukwa ichi ndi Full Waveform Inversion (FWI), njira yosinthira yomwe imagwiritsa ntchito chidziwitso chonse kuchokera ku mafunde a seismic kuti iyerekezere magawo enieni a pansi pa nthaka.
FWI nthawi zambiri imatchedwa "muyezo wagolide" wa kusinthasintha kwa seismic chifukwa imatha kupanga zithunzi zapansi panthaka zapamwamba kwambiri, kuposa njira zachikhalidwe zomwe zimagwiritsa ntchito nthawi yoyenda kapena ma amplitude osavuta. Komabe, mphamvu iyi imabwera ndi mtengo wake: FWI imafuna kupanga mafunde molondola, zida zowerengera zambiri, komanso njira zowongolera mosamala kuti tipewe kugwera mu mayankho olakwika.
Nkhaniyi ikufotokoza mfundo zoyambira za FWI, zigawo zazikulu za njira yake yogwiritsira ntchito, komanso zovuta ndi njira zomwe zimagwiritsidwa ntchito pokhazikitsa.
-
2. Kodi Full Waveform Inversion ndi chiyani?
Mwachidule, FWI ndi njira yopezera chitsanzo cha pansi pa nthaka chomwe chikugwirizana bwino ndi deta yoyeserera ya seismic ndi deta yowonedwa ya seismic. "Kuyenerera kwabwino kwambiri" kumatanthauzidwa ndi ntchito yolunjika (yosayenera), monga kusiyana kwa masikweya pakati pa deta yowonedwa ndi yopangidwa nthawi iliyonse komanso pa wolandila aliyense.
Kusiyana kwakukulu pakati pa FWI ndi njira zina zosinthira ndi:
- FWI imagwiritsa ntchito mawonekedwe onse a mafunde (gawo, kuchuluka, kusokoneza, kuchulukitsa, kusinthasintha), osati kungosankha nthawi yofika.
- FWI imadalira kuthetsa mafunde (acoustic, elastic, kapena anelastic) kuti iyerekezere deta yopangidwa.
– FWI ndi vuto lalikulu la kukonza zinthu mopanda mzere, chifukwa magawo a chitsanzo (monga kuthamanga kwa P- kapena S-wave) amatha kuwerengera mamiliyoni a maselo pa gridi ya 2D/3D.
-
3. Zigawo Zazikulu mu FWI
3.1 Deta Yoyang'anitsitsa ndi Deta Yopangidwa
FWI imafuna:
- Deta yowonera: zojambula za zivomerezi za m'munda (kusonkhanitsa zipolopolo) kuchokera ku magwero osiyanasiyana ndi olandira.
- Deta yopangidwa: zotsatira za kuyerekezera kwa manambala kwa kufalikira kwa mafunde pa chitsanzo cha kanthawi.
Deta yopangidwa imawerengedwa pothetsa equation ya mafunde (monga equation ya mafunde a acoustic):
\[
\frac{1}{v^2(\mathbf{x})}\frac{\partial^2 p}{\partial t^2} – \nabla^2 p = s(\mathbf{x},t)
\]
kumene \( v(\mathbf{x}) \) ndi liwiro la mafunde, \( p \) ndi kuthamanga/mafunde a scalar, ndipo \( s \) ndiye gwero.
3.2 Chitsanzo Choyamba
FWI imakhudzidwa kwambiri ndi chitsanzo choyamba. Ngati chitsanzo choyamba chili kutali kwambiri ndi momwe zinthu zilili, njira yogwiritsira ntchito njira imeneyi imatha kusinthasintha nthawi ndi nthawi, zomwe zimapangitsa kuti kusiyana kwa gawo pakati pa deta yopangidwa ndi yowonedwa kupitirire theka la nthawi, zomwe zimapangitsa kuti kusintha kwa njirayo kukankhire njirayo mbali yolakwika.
Chitsanzo choyamba nthawi zambiri chimapezeka kuchokera ku:
- kusintha kwa nthawi yoyenda (tomography),
- zitsanzo za liwiro lalikulu kuchokera ku kusanthula kwa geological,
- kapena njira yogwiritsira ntchito ma scale ambiri (kuyambira pa ma frequency otsika).
3.3 Ntchito Yoyenera (Yosayenera)
Ntchito zofanana:
\[
J(m)=\frac{1}{2}\sum_{s}\sum_{r}\int (d_{\text{syn}}(t; m)-d_{\text{obs}}(t))^2 \, dt
\]
kumene \( m \) ndi chitsanzo cha chitsanzo (monga liwiro), \( s \) ndi chiyambi cha chizindikiro, ndipo \( r \) ndi chizindikiro cha wolandila.
Kupatula pa L2 misfit yakale, palinso njira zina zochepetsera kulumpha kwa cycle, mwachitsanzo:
- zosakwanira pogwiritsa ntchito envelopu,
- zosakwanira gawo lokha,
- zoyendera bwino (Wasserstein),
– kapena zosefera zofanana.
3.4 Kuwerengera kwa Magiredi: Njira Yogwirizana ndi Boma
Chimodzi mwa zinthu zofunika kwambiri pa FWI ndi momwe mungawerengere bwino ma gradients. Chifukwa chiwerengero cha ma parameter a chitsanzo ndi chachikulu kwambiri, kuchotsa manambala mwachindunji n'kosatheka. Yankho lake ndi njira ya adjoint-state.
Chithunzi chosavuta kumva:
1. Kupangira chitsanzo cha kutsogolo: werengerani gawo la mafunde amtsogolo kuchokera ku gwero lomwe lili mu chitsanzo chamakono.
2. Werengani deta yotsalira \( \Delta d = d_{\text{syn}} – d_{\text{obs}} \).
3. Kupanga chitsanzo chogwirizana: ikani chotsaliracho ngati "gwero lobwerera" kuchokera pamalo olandirira ndikuchifalitsa mmbuyo munthawi yake.
4. Kugwirizana pakati pa minda ya mafunde a kutsogolo ndi adjoint kumapanga kusintha kwa ma parameter a chitsanzo.
Ndi njira iyi, mtengo wa makompyuta a gradient ndi wofanana ndi kawiri kuposa ma wave modeling pa gwero lililonse (forward + adjoint), kotero ndi waukulu koma wotheka pa HPC/GPU.
3.5 Ndondomeko Yosinthira Chitsanzo (Kukonza)
Mukapeza gradient \( \nabla J \) , chitsanzocho chimasinthidwa pogwiritsa ntchito njira yokonzera, mwachitsanzo:
- Kutsika Kwambiri (kosavuta),
- Conjugate Gradient,
– L-BFGS (yogwiritsidwa ntchito kwambiri chifukwa ndi yothandiza pamavuto akuluakulu),
– kapena njira ya Newton/Quasi-Newton.
Zosintha zoyambira:
\[
m_{k+1} = m_k – \alpha_k \, H_k^{-1}\nabla J(m_k)
\]
kumene \( \alpha_k \) ndi kutalika kwa sitepe , ndipo \( H_k^{-1} \) ndi kuyandikira kwa Hessian kotsutsana (monga mu L-BFGS).
-
4. Kayendedwe ka Ntchito ka FWI Algorithm (Kachidule)
Kawirikawiri, FWI imachitika mobwerezabwereza:
1. Sankhani chitsanzo choyambirira \( m_0 \).
2. Pa gwero lililonse:
– chitani chitsanzo chamtsogolo → deta yopangidwa,
- werengerani zotsalira poyerekeza ndi deta yowonera,
- chitani chitsanzo chogwirizana,
- kusonkhanitsa ma gradient.
3. Ikani chowongolera (monga chowongolera kuwala kapena kusalala).
4. Sinthani chitsanzocho ndi njira zowongolera.
5. Bwerezani mpaka kufika pamlingo wokwanira kapena wofanana.
Kawirikawiri FWI imayendetsedwa ndi ma scale ambiri, kuyambira pa ma frequency otsika (kukonza zigawo zazikulu za chitsanzo) kenako nkukwera kupita ku ma frequency apamwamba (kuwonjezera tsatanetsatane).
-
5. Mavuto Aakulu mu FWI
5.1 Kudumpha Njinga
Ili ndiye vuto lodziwika bwino. Ngati deta yopangidwa ndi yomwe yawonedwa siili "mu gawo," kusakwanira kwa L2 kungayambitse kukonzedwanso kwa minima yabodza. Yankho lalikulu:
- kuyambira pa ma frequency otsika kwambiri,
- kukonza chitsanzo choyambirira (tomography),
- kugwiritsa ntchito njira zina zosayenerera,
- gwiritsani ntchito njira zowerengera nthawi ndi njira zosankhira deta.
5.2 Mtengo wa Makompyuta
Kujambula zithunzi za 3D FWI zokhala ndi magwero osiyanasiyana kungafunike mafunde ambirimbiri kapena mamiliyoni ambiri. Izi zimafuna:
- makompyuta ofanana (gulu, GPU),
- kusunga kukumbukira (kufufuza kuti musunge mafunde),
- komanso njira zosankhira ma source subset (source encoding kapena mini-batch ala machine learning).
5.3 Kusagwirizana kwa Fiziki ndi Phokoso
Deta ya m'munda ili ndi phokoso, zotsatira za zida, anisotropy, attenuation (Q), complex topography, ndi kusatsimikizika kwa gwero. Ngati forward modeling ndi yosavuta kwambiri (monga, acoustic pamene medium ndi elastic), zotsatira za inversion zitha kukhala zosagwirizana.
5.4 Ma Parameter a Chitsanzo
Kusankha magawo oti musinthe (monga, \( v_p \), kuchuluka, \( v_s \), anisotropy, Q) kumakhudza kukhudzika ndi kukhazikika. Magawo ambiri angayambitse kusinthana ndikuwonjezera kulumikizana.
-
6. Kugwiritsa Ntchito FWI
FWI imagwiritsidwa ntchito kwambiri pa:
- kufufuza mafuta ndi gasi: kukonza njira yodziwira liwiro la magalimoto, kuthandiza kusamuka ndi kutanthauzira kapangidwe kake,
– kutentha kwa dziko: kujambula madera osweka ndi kusintha kwa lithological,
- seismology yapadziko lonse/yachigawo: kujambula zithunzi za mantle ndi crust,
- uinjiniya wa geotechnical ndi pafupi ndi pamwamba: kupanga mapu osaya a maziko, ngalande, ndi kuchepetsa kuwononga nthaka.
Ubwino wa FWI umaonekera bwino pamene njirayo ndi yovuta ndipo deta ili ndi zambiri zokhudza mafunde: ma diffraction ang'onoang'ono, njira zambiri, ndi ma multiples omwe nthawi zambiri amaonedwa kuti ndi osokoneza akhoza kukhala magwero a chidziwitso.
-
7. Mapeto
Full Waveform Inversion ndi njira yosinthira mafunde yomwe imagwiritsa ntchito chidziwitso chochuluka chomwe chili mu ma waveforms kuti ipange ma high-resolution subsurface models. Chinsinsi chake chopambana chili mu ma wave equation modeling olondola, ma gradient calculation ogwira ntchito pogwiritsa ntchito njira ya adjoint-state, ndi njira zowongolera zomwe zimapewa minima yapafupi monga cycle skipping. Ngakhale kuti imafuna makompyuta ambiri komanso kapangidwe kake kosamala, FWI yatsimikizira kuti ndi njira yamphamvu kwambiri pa geoimaging ndi kufufuza zinthu.
Ngati mukufuna, nditha kupitiliza ndi nkhani zina zaukadaulo—monga, kukambirana za adjoint gradient derivation mwalamulo, zitsanzo za FWI pseudocode, kapena njira zambiri zogwiritsira ntchito band-pass ndi windowing zogwiritsa ntchito ma multiscale kuti mupeze deta yeniyeni.
