(第九講)地震層析成象-資料下載頁
2024-11-19 01:34本頁面
【正文】 地震走時層析成像算法,當然，不一定非取 的L2模極小不可，也可以取任意階的模，如L2p模極小來求慢度的修改增量。其中P＝1，2，……。但是當P1時涉及開方運算，速度太慢，一般很少采。只有在令 時，可導出ART迭代的最簡單修正公式 （8）,第二十六頁，共五十頁。,三、地震走時層析成像算法,上式說明走時差平均地分配給每一條射線j通過地單元，而不考慮像元內(nèi)射線地長短。由式（6）和式（8）可得ART方法迭代公式,式中： 為系數(shù)矩陣 的分量；k=0,1,…為迭代次數(shù)。當 時，取 ，其中ε為給定允許誤差。,第二十七頁，共五十頁。,三、地震走時層析成像算法,3.2.3 SIRT算法（simultaneous iterative reconstruction technique） SIRT算法就是聯(lián)合迭代重建技術，用一個或多個優(yōu)化準則使得的解估計唯一。按照二次最優(yōu)化準則，得到SIRT典型的迭代修正公式 （11） 式中： ， ， 為非負實數(shù)；u為松弛系數(shù)；W為對稱矩陣；A為系數(shù)矩陣；B為非負定矩陣；C為對稱正定矩陣；B和C又稱為平滑矩陣；I為單位矩陣； 為與 同維的向量，它與圖像數(shù)值的某種先驗知識有關。,第二十八頁，共五十頁。,三、地震走時層析成像算法,矩陣W,B,C和向量 要根據(jù)實際問題來選擇，如對最 小二乘法和最小范數(shù)解估計，選 ，B=I, 及 , 可作為初始猜測（ 為數(shù)據(jù)協(xié)方差矩陣 的逆、 為圖像向量的協(xié)方差矩陣）。,第二十九頁，共五十頁。,三、地震走時層析成像算法,可以看出，SIRT重建算法的特點是在某一輪迭代中，所有像元的值（即圖像向量 ）都用前一輪的迭代結果 來修正。為了得到最簡單情況下的SIRT迭代修正公式，可令式（11）中 ＝ ＝0，適當選取矩陣W，C和重新定義不隨迭代變化的松弛系數(shù)u，則式（11）可寫為,,第三十頁，共五十頁。,四、實例,第三十一頁，共五十頁。,四、實例,第三十二頁，共五十頁。,四、實例,第三十三頁，共五十頁。,四、實例,第三十四頁，共五十頁。,四、實例,第三十五頁，共五十頁。,四、實例,第三十六頁，共五十頁。,四、實例,第三十七頁，共五十頁。,四、實例,第三十八頁，共五十頁。,四、實例,第三十九頁，共五十頁。,四、實例,第四十頁，共五十頁。,四、實例,第四十一頁，共五十頁。,四、實例,第四十二頁，共五十頁。,四、實例,第四十三頁，共五十頁。,四、實例,第四十四頁，共五十頁。,四、實例,第四十五頁，共五十頁。,四、實例,第四十六頁，共五十頁。,四、實例,第四十七頁，共五十頁。,四、實例,第四十八頁，共五十頁。,四、實例,,第四十九頁，共五十頁。,內(nèi)容總結,地震勘探新方法技術。四、實例。此后，地學界借助醫(yī)學CT思想，利用地震波的傳播對地殼乃至上地幔結構開始進行半定量研究。求解方程組（3）的最簡單和粗糙的離散圖像重構方法之一是反投影法。代數(shù)重建是按射線依次修改有關像元的圖像向量的一類迭代算法。上式說明走時差平均地分配給每一條射線j通過地單元，而不考慮像元內(nèi)射線地長短,第五十頁，共五
化學相關推薦
鄂ICP備17016276號-1