Tam Dalğa Formasının İnversiyası (FWI) Alqoritminə Giriş
1. Pendahuluan
Müasir geofizikada Yer kürəsinin yeraltı strukturunu "görməyə" ehtiyac artır. Enerji sənayesi, fəlakətlərin azaldılması, geotermal kəşfiyyat və hətta tektonik tədqiqatlar ətraflı yeraltı modellər tələb edir: seysmik dalğa sürətinin dərinliklə necə dəyişdiyi, süxur təbəqəsinin sərhədlərinin harada yerləşdiyi və kiçik heterojenliklərin dalğaların yayılmasına necə təsir göstərə biləcəyi. Bu məqsəd üçün ən güclü yanaşmalardan biri, yeraltı fiziki parametrləri qiymətləndirmək üçün seysmik dalğa formalarından tam məlumatlardan istifadə edən inversiya alqoritmi olan Tam Dalğa Formasının İnversiyasıdır (FWI).
FWI tez-tez seysmik inversiyanın "qızıl standartı" adlandırılır, çünki o, yalnız səyahət vaxtı və ya sadə amplituda istifadə edən ənənəvi yanaşmaları üstələyir və yüksək qətnaməli yeraltı görüntülər yarada bilir. Lakin bu güc öz qiymətini ödəyir: FWI səhv həllərə düşməmək üçün dəqiq dalğa modelləşdirməsi, əhəmiyyətli hesablama resursları və diqqətli optimallaşdırma strategiyaları tələb edir.
Bu məqalədə FWI-nin əsas anlayışları, alqoritminin əsas komponentləri və onun tətbiqindəki ümumi çətinliklər və strategiyalar təqdim olunur.
-
2. Tam Dalğa Formasının İnversiyası nədir?
Sadə dillə desək, FWI, simulyasiya edilmiş seysmik məlumatları müşahidə olunmuş seysmik məlumatlarla ən yaxşı şəkildə uyğunlaşdıran yeraltı modelin tapılması prosesidir. "Ən yaxşı uyğunluq" obyektiv funksiya (uyğunsuzluq) ilə müəyyən edilir, məsələn, hər zaman nöqtəsində və hər qəbuledicidə müşahidə olunmuş və sintetik məlumatlar arasındakı kvadrat fərq.
FWI və digər inversiya metodları arasındakı əsas fərqlər aşağıdakılardır:
– FWI yalnız gəliş vaxtının seçilməsindən deyil, bütün dalğa formasından (faza, amplituda, müdaxilə, çoxluq, difraksiya) istifadə edir.
– FWI sintetik məlumatları simulyasiya etmək üçün dalğa tənliyinin (akustik, elastik və ya anelastik) həllinə əsaslanır.
– FWI genişmiqyaslı qeyri-xətti optimallaşdırma problemidir, çünki model parametrləri (məsələn, P- və ya S-dalğa sürətləri) 2D/3D şəbəkəsində milyonlarla hüceyrəni saya bilər.
-
3. FWI-də əsas komponentlər
3.1 Müşahidə Məlumatları və Sintetik Məlumatlar
FWI tələb edir:
– Müşahidə məlumatları: müxtəlif mənbələrdən və qəbuledicilərdən sahə seysmik qeydləri (atışların toplanması).
– Sintetik məlumatlar: müvəqqəti model üzərində dalğa yayılmasının ədədi simulyasiyasının nəticələri.
Sintetik məlumatlar dalğa tənliyini (məsələn, akustik dalğa tənliyini) həll etməklə hesablanır:
\[
\frac{1}{v^2(\mathbf{x})}\frac{\partial^2 p}{\qismən t^2} – \nabla^2 p = s(\mathbf{x},t)
\]
burada \(v(\mathbf{x}) \) dalğa sürəti, \(p \) skalyar təzyiq/dalğa və \(s \) mənbədir.
3.2 İlkin Model
FWI ilkin modelə çox həssasdır. İlkin model faktiki şərtlərdən çox uzaqdırsa, alqoritm dövr atlaması ilə qarşılaşa bilər, bu da sintetik və müşahidə edilən məlumatlar arasındakı faza fərqinin yarım dövrdən çox olması və qradiyentin həlli səhv istiqamətə itələməsinə səbəb olur.
İlkin model adətən aşağıdakılardan əldə edilir:
– səyahət vaxtı inversiyası (tomoqrafiya),
– geoloji analizdən əldə edilən makro sürət modelləri,
– və ya çoxmiqyaslı yanaşma (aşağı tezliklərdən başlayaraq).
3.3 Məqsəd Funksiyası (Uyğunsuzluq)
Ümumi məqsəd funksiyaları:
\[
J(m)=\frac{1}{2}\sum_{s}\sum_{r}\int (d_{\text{syn}}(t; m)-d_{\text{obs}}(t))^2 \, dt
\]
burada \(m\) model parametridir (məsələn, sürət), \(s\) mənbə indeksidir və \(r\) qəbuledici indeksidir.
Klassik L2 uyğunsuzluğundan əlavə, velosiped atlamalarını azaltmaq üçün alternativ uyğunsuzluqlar da mövcuddur, məsələn:
– zərf əsaslı uyğunsuzluq,
– yalnız faza uyğunsuzluğu,
- optimal nəqliyyat (Wasserstein),
– və ya uyğun filtrlər.
3.4 Qradiyent Hesablanması: Qoşulma Vəziyyəti Metodu
FWI-nin əsas xüsusiyyətlərindən biri qradiyentləri necə səmərəli hesablamaqdır. Model parametrlərinin sayı çox olduğundan, birbaşa ədədi çıxarış mümkün deyil. Həll yolu qoşma vəziyyət metodudur.
İntuitiv şəkil:
1. İrəli modelləşdirmə: cari modeldə mənbədən irəli dalğa sahəsini hesablayın.
2. Məlumat qalığını hesablayın \( \Delta d = d_{\text{syn}} – d_{\text{obs}} \).
3. Qoşma modelləşdirmə: qalığı qəbuledici mövqeyindən “qayıdış mənbəyi” kimi yeridin və zamanla geriyə doğru yayın.
4. İrəli və bitişik dalğa sahələri arasındakı korrelyasiya model parametrlərinə bir qradiyent yaradır.
Bu texnika ilə qradiyentin hesablama dəyəri təxminən hər mənbəyə düşən dalğa modelləşdirməsinin ikiqatına bərabərdir (irəli + qoşma), buna görə də hələ də böyükdür, lakin HPC/GPU-da mümkündür.
3.5 Model Yeniləmə Sxemi (Optimallaşdırma)
Qradiyent \( \nabla J \) əldə edildikdən sonra, model optimallaşdırma metodu ilə yenilənir, məsələn:
– Ən dik eniş (ən sadə),
– Qoşulma Qradiyenti,
– L-BFGS (böyük problemlər üçün səmərəli olduğuna görə geniş istifadə olunur),
– və ya Nyuton/Kvazi-Nyuton metodu.
Əsas yeniləmələr:
\[
m_{k+1} = m_k – \alpha_k \, H_k^{-1}\nabla J(m_k)
\]
burada \( \alpha_k \) addım uzunluğu, \( H_k^{-1} \) isə tərs Hess yaxınlaşmasıdır (məsələn, L-BFGS-də).
-
4. FWI Alqoritmi İş Axını (Qısa)
Ümumiyyətlə, FWI təkrarlanan şəkildə işləyir:
1. İlkin modeli seçin \( m_0 \).
2. Hər bir mənbə üçün:
– irəli modelləşdirmə → sintetik məlumatlar edin,
– qalıqları müşahidə məlumatlarına əsasən hesablayın,
– əlavə modelləşdirmə aparmaq,
– qradiyent yığılması.
3. İlkin şərtlər tətbiq edin (məsələn, işıqlandırma kompensasiyası və ya hamarlaşdırma).
4. Modeli optimallaşdırma metodları ilə yeniləyin.
5. Konvergensiya və ya iterasiya həddinə çatana qədər təkrarlayın.
Tipik olaraq, FWI aşağı tezliklərdən başlayaraq (böyük miqyaslı model komponentlərini korreksiya etməklə) və daha yüksək tezliklərə keçərək (detallar əlavə etməklə) çoxmiqyaslı şəkildə idarə olunur.
-
5. FWI-də əsas çətinliklər
5.1 Velosipeddən Atlama
Bu, ən çox bilinən problemdir. Sintetik və müşahidə edilən məlumatlar "fazada" olmadıqda, L2 uyğunsuzluğu optimallaşdırmanın yalançı lokal minimuma çatmasına səbəb ola bilər. Ümumi həll yolu:
– çox aşağı tezliklərdən başlayaraq,
– ilkin modeli (tomoqrafiya) təkmilləşdirmək,
– alternativ uyğunsuzluqlardan istifadə etməklə,
– zaman pəncərəsi və məlumat seçimi strategiyalarını tətbiq edin.
5.2 Hesablama Xərci
Çoxsaylı mənbələrə malik 3D FWI minlərlə ilə milyonlarla dalğa forması simulyasiyasını tələb edə bilər. Bunun üçün aşağıdakılar tələb olunur:
– paralel hesablama (klaster, GPU),
– yaddaş qənaəti (dalğa sahələrini saxlamaq üçün yoxlama nöqtəsi),
– həmçinin mənbə alt çoxluğu seçim strategiyaları (mənbə kodlaşdırması və ya mini-toplu maşın öyrənməsi).
5.3 Fizika və səs-küyün uyğunsuzluğu
Sahə məlumatları səs-küy, alət effektləri, anizotropiya, zəifləmə (Q), mürəkkəb topoqrafiya və mənbə qeyri-müəyyənliyini ehtiva edir. İrəli modelləşdirmə çox sadədirsə (məsələn, mühit elastik olduqda akustik), inversiya nəticələri qərəzli ola bilər.
5.4 Model Parametrləşdirməsi
Hansı parametrlərin tərsinə çevriləcəyini seçmək (məsələn, \( v_p \), sıxlıq, \( v_s \), anizotropiya, Q) həssaslığa və sabitliyə təsir göstərir. Həddindən artıq çox parametr kompromislərə səbəb ola və konvergensiyanı pisləşdirə bilər.
-
6. FWI Ərizəsi
FWI aşağıdakılar üçün geniş istifadə olunur:
– neft və qaz kəşfiyyatı: sürət modellərinin həllini təkmilləşdirmək, strukturların miqrasiyasına və interpretasiyasına kömək etmək,
– geotermal: çatlama zonalarının və litoloji dəyişikliklərin xəritələşdirilməsi,
– qlobal/regional seysmologiya: mantiya və yer qabığının görüntülənməsi,
– geotexniki və səthə yaxın mühəndislik: bünövrələr, tunellər və təsirlərin azaldılması üçün dayaz xəritələşdirmə.
FWI-nin üstünlükləri mühit mürəkkəb olduqda və məlumatlar dalğa məlumatlarına zəngin olduqda aydın olur: adətən müdaxilə hesab edilən kiçik difraksiyalar, çoxyolluluqlar və çoxluqlar əslində məlumat mənbəyi ola bilər.
-
7. Nəticə
Tam Dalğa Formasının İnversiyası, yüksək qətnaməli yeraltı modellər qurmaq üçün dalğa formalarında olan zəngin məlumatlardan istifadə edən seysmik inversiya alqoritmidir. Onun uğurunun açarı dəqiq dalğa tənlik modelləşdirməsində, qoşma vəziyyət metodundan istifadə edərək səmərəli qradiyent hesablamalarında və dövr atlamaları kimi lokal minimumların qarşısını alan optimallaşdırma strategiyalarındadır. Yüksək hesablama tələblərinə və diqqətli iş axını dizaynına baxmayaraq, FWI geogörüntüləmə və resurs kəşfiyyatında güclü bir yanaşma olduğunu sübut etmişdir.
İstəsəniz, daha çox texniki məqalə ilə davam edə bilərəm - məsələn, qoşma qradiyent törəməsini daha rəsmi şəkildə müzakirə etmək, FWI psevdokod nümunələri və ya real məlumatlar üçün zolaq keçirmə və pəncərə əsaslı çoxmiqyaslı strategiyalar.
