快速ＴＡＵＰ變換模塊在ＦＯＣＵＳ下的實現(xiàn)和應用

在目前地震資料處理軟件中，因FOCUS圖形界面好，工具菜單布置合理，相對其它處理軟件所占比例要高，然而有些模塊運算速度相對其它處理軟件運算速度要慢很多，F(xiàn)OCUS下處理人員要求能改進相應的模塊，使之提高處理速度。

近幾年，利用TAUP域處理來壓制多次波和提高信噪比已成為一項非常有特色的處理技術，應用在多個區(qū)塊的資料處理中，取得了非常好的效果。而在實際生產中發(fā)現(xiàn)，ITA系統(tǒng)下的TAUP變換模塊運算速度很快，而FOCUS系統(tǒng)下的TAUP變換模塊運算速度很慢，致使在用FOCUS系統(tǒng)進行批量生產時，TAUP域處理的效率很低，不利于提高生產效率。本文考慮采用一些快速算法，編寫FOCUS系統(tǒng)下的快速TAUP變換模塊。經過實際數據的測試對比，新編快速TAUP變換模塊在相同變換參數情況下的計算速度比原有模塊快1.3-3倍；可以提高TAUP域處理的效率。

1方法原理

對TAUP變換的原理，這里不在敘述，僅列出變換公式。設地震波場為 ， 為變換波場，TAUP正、反變換由下式定義：

正變換 ：

從時間域的直線變換到TAUP域為一點，雙曲線變換到TAUP為橢圓，如圖1、圖2。利用同相軸在TAUP域的特點，可以壓制線性干擾、多次波以及波場分離等處理。

2算法實現(xiàn)

TAUP變換算法可以在時間－空間域和頻率－波數域實現(xiàn)，且正反變換過程可逆。

時間－空間域：

為了提高變換速度，用 C語言編寫了FFT、褶積、插值等子程序，與Fortran語言混合編程，實現(xiàn)FOCUS系統(tǒng)下的快速TAUP變換模塊，運算速度大幅度提高。另外，增加了MASK功能，可以有效地壓制變換過程中高p值部位引起的假頻現(xiàn)象。……