Муқаддима ба алгоритми инверсияи пурраи мавҷ

Муқаддима ба алгоритми инверсияи пурраи мавҷ (FWI)

1. Пендахулуан

Дар геофизикаи муосир, ниёз ба "дидан"-и сохтори зеризаминии Замин меафзояд. Саноати энергетика, коҳиши офатҳои табиӣ, ҷустуҷӯи геотермалӣ ва ҳатто таҳқиқоти тектоникӣ ба моделҳои муфассали зеризаминӣ ниёз доранд: чӣ гуна суръати мавҷҳои сейсмикӣ бо умқ тағйир меёбад, сарҳадҳои қабатҳои сангӣ дар куҷо ҷойгиранд ва чӣ гуна гетерогении хурд метавонад ба паҳншавии мавҷ таъсир расонад. Яке аз пуриқтидортарин равишҳо барои ин мақсад табдили пурраи мавҷ (FWI) мебошад, ки алгоритми табдилест, ки маълумоти пурраро аз шаклҳои мавҷҳои сейсмикӣ барои арзёбии параметрҳои физикии зеризаминӣ истифода мебарад.

FWI аксар вақт ҳамчун "стандарти тиллоӣ"-и инверсияи сейсмикӣ номида мешавад, зеро он метавонад тасвирҳои зеризаминии баландсифатро тавлид кунад ва аз усулҳои анъанавӣ, ки танҳо вақти сафар ё амплитудаи оддиро истифода мебаранд, зиёдтар бошад. Аммо, ин қудрат бо нарх меояд: FWI ба моделсозии дақиқи мавҷҳо, захираҳои назарраси ҳисоббарорӣ ва стратегияҳои бодиққати оптимизатсия ниёз дорад, то аз дучор шудан ба роҳҳои ҳалли хато пешгирӣ карда шавад.

Дар ин мақола мафҳумҳои асосии FWI, ҷузъҳои асосии алгоритми он ва мушкилот ва стратегияҳои умумӣ дар татбиқи он муаррифӣ карда мешаванд.

-

2. Инверсияи пурраи шакли мавҷ чист?

Ба таври оддӣ, FWI раванди ёфтани модели зеризаминӣ аст, ки маълумоти сейсмикии симулятсияшударо бо маълумоти сейсмикии мушоҳидашуда беҳтарин мувофиқат мекунад. "Мувофиқати беҳтарин" бо функсияи объективӣ (нодурустӣ), ба монанди фарқи квадратии байни маълумоти мушоҳидашуда ва синтетикӣ дар ҳар нуқтаи вақт ва дар ҳар як қабулкунанда муайян карда мешавад.

Фарқиятҳои асосии байни FWI ва дигар усулҳои инверсия инҳоянд:
– FWI на танҳо интихоби вақти расиданро, балки тамоми шакли мавҷро (фаза, амплитуда, интерференсия, зарб, дифракция) истифода мебарад.
– FWI барои тақлид кардани маълумоти синтетикӣ ба ҳалли муодилаи мавҷӣ (акустикӣ, эластикӣ ё анеластикӣ) такя мекунад.
– FWI як масъалаи оптимизатсияи ғайрихаттии миқёси калон аст, зеро параметрҳои модел (масалан, суръатҳои мавҷи P ё S) метавонанд миллионҳо ҳуҷайраҳоро дар шабакаи 2D/3D шумуранд.

Хонед  Фаҳмиши асосии анизотропияи сейсмикӣ

-

3. Ҷузъҳои асосӣ дар FWI

3.1 Маълумоти мушоҳидавӣ ва маълумоти синтетикӣ
FWI талаб мекунад:
– Маълумоти мушоҳида: сабтҳои сейсмикии саҳроӣ (ҷамъоварии тир) аз манбаъҳо ва қабулкунакҳои гуногун.
– Маълумоти синтетикӣ: натиҷаҳои моделсозии рақамии паҳншавии мавҷ дар модели муваққатӣ.

Маълумоти синтетикӣ бо роҳи ҳалли муодилаи мавҷ (масалан, муодилаи мавҷи акустикӣ) ҳисоб карда мешавад:

\[
\frac{1}{v^2(\mathbf{x})}\frac{\partial^2 p}{\qisman t^2} – \nabla^2 p = s(\mathbf{x},t)
\]

ки дар он \(v(\mathbf{x}) \) суръати мавҷ, \(p \) фишори скалярӣ/мавҷ ва \(s \) манбаъ аст.

3.2 Модели ибтидоӣ
FWI ба модели ибтидоӣ хеле ҳассос аст. Агар модели ибтидоӣ аз шароити воқеӣ хеле дур бошад, алгоритм метавонад гузариши даврро аз сар гузаронад, ки ин вақте рух медиҳад, ки фарқияти фаза байни маълумоти синтетикӣ ва мушоҳидашуда беш аз ним давра аст, ки боиси тела додани градиенти ҳалли масъала ба самти нодуруст мегардад.

Модели ибтидоӣ одатан аз инҳо гирифта мешавад:
– инверсияи вақти сафар (томография),
– моделҳои суръати макро аз таҳлили геологӣ,
– ё равиши бисёрҷониба (аз басомадҳои паст сар карда).

3.3 Функсияи мақсаднок (Номувофиқӣ)
Функсияҳои мақсадноки маъмул:

\[
J(m)=\frac{1}{2}\sum_{s}\sum_{r}\int (d_{\text{syn}}(t; m)-d_{\text{obs}}(t))^2 \, dt
\]

ки дар он \( m \) параметри модел аст (масалан, суръат), \( s \) индекси манбаъ ва \( r \) индекси қабулкунанда аст.

Илова бар номутобиқатии классикии L2, номутобиқатиҳои алтернативӣ низ барои кам кардани партофтани даврҳо мавҷуданд, масалан:
– номувофиқатии бар асоси лифофа,
– номувофиқатии танҳо марҳила,
- интиқоли оптималӣ (Вассерштейн),
– ё филтрҳои мувофиқ.

3.4 Ҳисобкунии градиент: Усули ҳолати пайваст
Яке аз хусусиятҳои асосии FWI ин ҳисоб кардани самараноки градиентҳо мебошад. Азбаски шумораи параметрҳои модел хеле зиёд аст, ҳосилаи мустақими рақамӣ ғайриимкон аст. Ҳалли масъала усули ҳолати пайвастшавӣ мебошад.

Тасвири интуитивӣ:
1. Моделсозии мустақим: майдони мавҷи мустақимро аз манбаъ дар модели ҷорӣ ҳисоб кунед.
2. Бақияи маълумотро ҳисоб кунед \( \Delta d = d_{\text{syn}} – d_{\text{obs}} \).
3. Моделсозии иловагӣ: боқимондаро ҳамчун "манбаи бозгашт" аз мавқеи қабулкунанда ворид кунед ва онро дар тӯли вақт ба ақиб паҳн кунед.
4. Ҳамбастагии байни майдонҳои мавҷи рост ва ҳамшафат градиентро ба параметрҳои модел ба вуҷуд меорад.

Хонед  Усули поляризатсияи индуксионӣ дар ҷустуҷӯи маъданҳо

Бо ин усул, арзиши ҳисобкунии градиент тақрибан ба ду баробар моделсозии мавҷ барои як манбаъ (пеш + пайваст) баробар аст, аз ин рӯ он то ҳол калон аст, аммо дар HPC/GPU имконпазир аст.

3.5 Нақшаи навсозии модел (беҳсозӣ)
Пас аз ба даст овардани градиенти \( \nabla J \), модел бо истифода аз усули оптимизатсия навсозӣ мешавад, масалан:
– Нишебтарин фуруд (соддатарин),
– Градиенти конъюгатӣ,
– L-BFGS (одатан истифода мешавад, зеро он барои мушкилоти калон самаранок аст),
– ё усули Нютон/Квази-Ньютон.

Навсозиҳои асосӣ:

\[
m_{k+1} = m_k – \alpha_k \, H_k^{-1}\nabla J(m_k)
\]

ки дар он \( \alpha_k \) дарозии қадам аст ва \( H_k^{-1} \) тақриби баръакси Гессиан аст (масалан, дар L-BFGS).

-

4. Ҷараёни кории алгоритми FWI (мухтасар)

Умуман, FWI ба таври такрорӣ иҷро мешавад:

1. Модели ибтидоии \(m_0 \)-ро интихоб кунед.
2. Барои ҳар як манбаъ:
– моделсозии пешакӣ → маълумоти синтетикӣ анҷом диҳед,
– боқимондаҳоро бо маълумоти мушоҳида ҳисоб кунед,
- моделсозии муштаракро анҷом диҳед,
– ҷамъшавии градиент.
3. Пешшарткуниро истифода баред (масалан, ҷуброни равшанӣ ё ҳамворкунӣ).
4. Моделро бо усулҳои оптимизатсия навсозӣ кунед.
5. То расидан ба ҳадди конвергенсия ё лимити такрорӣ такрор кунед.

Одатан, FWI бисёрҷониба иҷро мешавад, ки аз басомадҳои паст оғоз мешавад (ислоҳи ҷузъҳои модели миқёси калон) ва сипас ба басомадҳои баландтар мегузарад (илова кардани тафсилот).

-

5. Мушкилоти асосӣ дар FWI

5.1 Гузариш аз даврзанӣ
Ин мушкили маъруфтарин аст. Вақте ки маълумоти синтетикӣ ва мушоҳидашуда "дар марҳила" нестанд, номувофиқатии L2 метавонад ба беҳсозӣ ба минимумҳои маҳаллии бардурӯғ оварда расонад. Ҳалли умумӣ:
– аз басомадҳои хеле паст сар карда,
– такмил додани модели ибтидоӣ (томография),
– истифодаи номувофиқатҳои алтернативӣ,
– истифодаи стратегияҳои тирезабандии вақт ва интихоби маълумот.

5.2 Арзиши ҳисоббарорӣ
FWI-и сеченака бо манбаъҳои гуногун метавонад аз ҳазорҳо то миллионҳо симулятсияи шакли мавҷро талаб кунад. Ин талаб мекунад:
– ҳисоббарории мувозӣ (кластер, GPU),
– захира кардани хотира (нуқтаи санҷиш барои захира кардани майдонҳои мавҷ),
– инчунин стратегияҳои интихоби зермаҷмӯи манбаъ (рамзгузории манбаъ ё омӯзиши мошинии хурд-пакетӣ).

Хонед  Харитасозии сохторҳои зеризаминӣ бо истифода аз GPR

5.3 Номувофиқатии физика ва садо
Маълумоти саҳроӣ дорои садо, таъсири асбобҳо, анизотропия, сустшавӣ (Q), топографияи мураккаб ва номуайянии манбаъ мебошанд. Агар моделсозии пешрафта хеле содда бошад (масалан, акустикӣ вақте ки муҳит чандир аст), натиҷаҳои инверсия метавонанд ғаразнок бошанд.

5.4 Параметризатсияи модел
Интихоби параметрҳое, ки бояд инверт карда шаванд (масалан, \( v_p \), зичӣ, \( v_s \), анизотропия, Q) ба ҳассосият ва устуворӣ таъсир мерасонад. Параметрҳои аз ҳад зиёд метавонанд боиси мувозинат шаванд ва ҳамгироиро бадтар кунанд.

-

6. Аризаи FWI

FWI ба таври васеъ барои инҳо истифода мешавад:
– ҷустуҷӯи нафт ва газ: беҳтар кардани ҳалли моделҳои суръат, мусоидат ба муҳоҷират ва тафсири сохторҳо,
– геотермалӣ: харитасозии минтақаҳои шикастагӣ ва тағйироти литологӣ,
– сейсмологияи ҷаҳонӣ/минтақавӣ: аксбардории мантия ва қишри замин,
– муҳандисии геотехникӣ ва наздик ба рӯизаминӣ: харитасозии наонқадар чуқур барои таҳкурсӣ, нақбҳо ва чораҳои коҳишдиҳанда.

Бартариҳои FWI вақте аёнанд, ки муҳит мураккаб аст ва маълумот аз маълумоти мавҷӣ бой аст: дифракцияҳои хурд, бисёрпайвастагӣ ва зарбҳо, ки одатан интерференсия ҳисобида мешаванд, дар асл метавонанд манбаи иттилоот бошанд.

-

7. Хулоса

Инверсияи пурраи мавҷ як алгоритми инверсияи сейсмикӣ мебошад, ки аз маълумоти бойи мавҷуд дар шаклҳои мавҷ барои сохтани моделҳои зеризаминии баландсифат истифода мебарад. Калиди муваффақияти он дар моделсозии дақиқи муодилаи мавҷ, ҳисобҳои самараноки градиент бо истифода аз усули ҳолати пайваст ва стратегияҳои оптимизатсия мебошад, ки аз минимумҳои маҳаллӣ, ба монанди гузарондани даврҳо, канорагирӣ мекунанд. Сарфи назар аз талаботи баланди ҳисоббарорӣ ва тарҳрезии бодиққати ҷараёни корӣ, FWI ҳамчун як равиши пуриқтидор дар геотасвиркунӣ ва ҷустуҷӯи захираҳо исбот шудааст.

Агар шумо хоҳед, ман метавонам бо мақолаҳои бештари техникӣ идома диҳам - масалан, баррасии расмӣтар дар бораи ҳосилаи градиенти ҳамроҳ, мисолҳои псевдокоди FWI ё стратегияҳои бисёрмиқёси банд-гузариш ва тирезасозӣ барои маълумоти воқеӣ.

Шарҳ гузоред