Entwodiksyon sou Algorit Envèsyon Fòm Vag Konplè (FWI) la
1. Pendahuluan
Nan jeofizik modèn, bezwen pou "wè" estrikti souteren Latè a ap ogmante. Endistri enèji a, mezi pou diminye dezas , eksplorasyon jeotèmal, e menm rechèch tektonik mande pou modèl souteren detaye: kijan vitès ond sismik yo chanje ak pwofondè, ki kote limit kouch wòch yo ye, ak kijan ti eterogeneite ka afekte pwopagasyon ond yo. Youn nan apwòch ki pi pwisan pou objektif sa a se Envèsyon Fòm Ond Konplè (FWI), yon algorithm envèsyon ki itilize enfòmasyon konplè ki soti nan fòm ond sismik yo pou estime paramèt fizik souteren yo.
Yo souvan refere a FWI kòm "etalon an lò" nan envèsyon sismik paske li ka pwodui imaj souteren ki gen gwo rezolisyon, depase apwòch konvansyonèl ki sèlman itilize tan vwayaj oswa anplitid senp. Sepandan, pouvwa sa a vini ak yon pri: FWI mande yon modèl vag egzat, resous enfòmatik sibstansyèl , ak estrateji optimize atansyon pou evite tonbe nan solisyon ki pa kòrèk.
Atik sa a prezante konsèp debaz FWI yo, prensipal konpozan algorithm li a, ansanm ak defi ak estrateji jeneral nan aplikasyon li.
-
2. Kisa Envèsyon Fòm Vag Konplè ye?
Senpleman, FWI se pwosesis pou jwenn yon modèl souteren ki pi byen matche ak done sismik simile yo ak done sismik obsève yo. "Pi bon ajisteman" defini pa yon fonksyon objektif (maladat), tankou diferans kare ant done obsève yo ak done sentetik yo nan chak pwen tan ak nan chak reseptè.
Diferans fondamantal ant FWI ak lòt metòd envèsyon yo se:
– FWI itilize tout fòm ond lan (faz, anplitid, entèferans, miltip, difraksyon), pa sèlman chwa lè arive a.
– FWI apiye sou rezoud ekwasyon ond lan (akoustik, elastik, oswa anelastik) pou simile done sentetik.
– FWI se yon pwoblèm optimize non lineyè sou gwo echèl, paske paramèt modèl yo (pa egzanp, vitès ond P oswa S) ka konte plizyè milyon selil sou yon griy 2D/3D.
-
3. Konpozan prensipal nan FWI
3.1 Done Obsèvasyonèl ak Done Sentetik
FWI egzije:
– Done obsèvasyon: anrejistreman sismik sou teren (rasanbleman imaj) ki soti nan plizyè sous ak reseptè.
– Done sentetik: rezilta simulation nimerik pwopagasyon ond sou yon modèl tanporè.
Done sentetik yo kalkile lè yo rezoud ekwasyon ond lan (pa egzanp, ekwasyon ond akoustik la):
\[
\frac{1}{v^2(\mathbf{x})}\frac{\partial^2 p}{\partial t^2} – \nabla^2 p = s(\mathbf{x},t)
\]
kote \(v(\mathbf{x}) \) se vitès vag la, \(p \) se presyon/vag eskalè a, epi \(s \) se sous la.
3.2 Modèl Inisyal
FWI trè sansib a modèl inisyal la. Si modèl inisyal la twò lwen kondisyon reyèl yo, algorithm nan ka fè eksperyans sote sik, ki se lè diferans faz ant done sentetik yo ak done obsève yo plis pase mwatye yon peryòd, sa ki lakòz gradyan an pouse solisyon an nan move direksyon.
Modèl inisyal la anjeneral jwenn nan:
– envèsyon tan vwayaj (tomografi),
– modèl vitès makro ki soti nan analiz jewolojik,
– oubyen yon apwòch milti-echèl (kòmanse ak frekans ki ba).
3.3 Fonksyon Objektif (Enkonpatibilite)
Fonksyon objektif komen:
\[
J(m)=\frac{1}{2}\sum_{s}\sum_{r}\int (d_{\text{syn}}(t; m)-d_{\text{obs}}(t))^2 \, dt
\]
kote \(m\) se yon paramèt modèl (pa egzanp, vitès), \(s\) se endèks sous la, epi \(r\) se endèks reseptè a.
Anplis defo klasik L2 a, genyen tou lòt defo pou diminye sote sik, pa egzanp:
– move anfòm ki baze sou anvlòp,
– enkonpatibilite faz sèlman,
- transpò pi bon (Wasserstein),
– oubyen filtè ki koresponn yo.
3.4 Kalkil Gradyan: Metòd Eta Adjwen
Youn nan karakteristik kle FWI se kijan pou kalkile gradyan yo avèk efikasite. Piske kantite paramèt modèl la tèlman gwo, derivasyon nimerik dirèk enposib. Solisyon an se metòd eta adjwen an.
Imaj entwitif la:
1. Modèl pou pi devan: kalkile chan vag pou pi devan an apati sous la nan modèl aktyèl la.
2. Kalkile rezidyèl done a \( \Delta d = d_{\text{syn}} – d_{\text{obs}} \).
3. Modèl adjwen: enjekte rezidyèl la kòm yon "sous retou" apati pozisyon reseptè a epi pwopaje li bak nan tan.
4. Korelasyon ki genyen ant chan vag dirèk ak adjwen yo pwodui yon gradyan nan paramèt modèl la.
Avèk teknik sa a, pri kalkil gradyan an apeprè ekivalan a de fwa modèl vag pou chak sous (pou pi devan + adjwen), kidonk li toujou gwo men posib sou HPC/GPU.
3.5 Plan Mizajou Modèl (Optimizasyon)
Yon fwa yo jwenn gradyan \( \nabla J \) a, yo mete modèl la ajou avèk yon metòd optimize, pa egzanp:
– Desant ki pi apik la (ki pi senp lan),
– Gradyan konjige,
– L-BFGS (souvan itilize paske li efikas pou gwo pwoblèm),
– oubyen metòd Newton/Quasi-Newton.
Mizajou debaz yo:
\[
m_{k+1} = m_k – \alpha_k \, H_k^{-1}\nabla J(m_k)
\]
kote \( \alpha_k \) se longè etap la, epi \( H_k^{-1} \) se apwoksimasyon envès Hessian an (pa egzanp nan L-BFGS).
-
4. Workflow Algorit FWI a (Konsiz)
An jeneral, FWI fonksyone iteratifman:
1. Chwazi modèl inisyal la \( m_0 \).
2. Pou chak sous:
– fè modèlizasyon pou pi devan → done sentetik,
– kalkile rezidyèl yo parapò ak done obsèvasyon yo,
– fè modèl adjwen,
– akimilasyon gradyan.
3. Aplike prekondisyone (pa egzanp, konpansasyon ekleraj oswa lisaj).
4. Mete modèl la ajou ak metòd optimize yo.
5. Repete jiskaske konvèjans lan oswa limit iterasyon an rive.
Tipikman, FWI fèt sou plizyè echèl, kòmanse nan frekans ki ba (korije konpozan modèl gwo echèl) epi answit monte nan frekans ki pi wo (ajoute detay).
-
5. Prensipal defi nan FWI
5.1 Sote Sik
Sa a se pwoblèm ki pi byen koni an. Lè done sentetik ak done obsève yo pa "an faz", move anfòm L2 ka mennen optimize a fo minimòm lokal yo. Solisyon jeneral la:
– kòmanse ak frekans ki ba anpil,
– amelyore modèl inisyal la (tomografi),
– lè l sèvi avèk lòt bagay ki pa anfòm,
– aplike estrateji pou divize an fenèt tan ak seleksyon done.
5.2 Pri Konpitasyonèl
3D FWI ak plizyè sous ka mande plizyè milye rive plizyè milyon similasyon fòm ond. Sa mande:
– informatique paralèl (gwoup, GPU),
– ekonomize memwa (pwen kontwòl pou konsève chan vag yo),
– ansanm ak estrateji seleksyon sou-ansanm sous (kodaj sous oswa mini-batch ala aprantisaj otomatik).
5.3 Fizik ak Bri Enkonpatibilite
Done teren yo gen ladan bri, efè enstriman, anizotropi, atenyasyon (Q), topografi konplèks, ak ensètitid sous. Si modèl pou pi devan an twò senp (pa egzanp, akoustik lè milye a elastik), rezilta envèsyon yo ka gen yon patipri.
5.4 Paramètizasyon Modèl
Chwazi paramèt pou envèse yo (pa egzanp, \(v_p \), dansite, \(v_s \), anizotropi, Q) afekte sansiblite ak estabilite. Twòp paramèt ka lakòz konpwomi epi vin pi mal konvèjans lan.
-
6. Aplikasyon FWI
FWI lajman itilize pou:
– eksplorasyon lwil oliv ak gaz: amelyore rezolisyon modèl vitès, ede migrasyon ak entèpretasyon estrikti yo,
– jeotèmal: katografi zòn frakti ak chanjman litolojik,
– sismoloji mondyal/rejyonal: imaj manto ak kwout,
– jeni jeoteknik ak jeni toupre sifas: katografi sifas pou fondasyon, tinèl, ak mitigasyon.
Avantaj FWI yo evidan lè mwayen an konplèks epi done yo rich nan enfòmasyon sou ond: ti difraksyon, miltiptrajektwa, ak miltip ke yo anjeneral konsidere kòm entèferans kapab aktyèlman sous enfòmasyon.
-
7. Konklizyon
Envèsyon Fòm Ond Konplè a se yon algorithm envèsyon sismik ki itilize enfòmasyon rich ki nan fòm ond yo pou konstwi modèl souteren ki gen gwo rezolisyon. Kle siksè li se nan modèl ekwasyon ond ki egzat, kalkil gradyan efikas lè l sèvi avèk metòd eta adjwen an, ak estrateji optimize ki evite minima lokal tankou sote sik. Malgre gwo demand enfòmatik li yo ak konsepsyon workflow atansyon, FWI te pwouve li se yon apwòch trè pwisan nan imaj Latè ak eksplorasyon resous.
Si ou vle, mwen ka kontinye ak plis atik teknik—pa egzanp, diskite sou derivasyon gradyan adjwen an pi fòmèlman, egzanp pseudokòd FWI, oswa estrateji miltiskal ki baze sou band-pass ak windowing pou done reyèl.
