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

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

1方法原理

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

正變換 ：

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

2算法實現(xiàn)

TAUP變換算法可以在時間－空間域和頻率－波數(shù)域?qū)崿F(xiàn)，且正反變換過程可逆。

時間－空間域：

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

3快速TAUP變換模塊使用說明

3.1基本功能

TAUP模塊將疊前數(shù)據(jù)從T-X域變換到T-P域，或從T-P域反變換到T-X域，與SLNFOW和SLNINV的功能一樣，但是這里采用了快速算法，計算速度比SLNFOW和SLNINV快，而且增加了MASK功能，可以有效地壓制變換過程中高P值部分引起的假頻。TAUP模塊提供兩種算法，一種在T-X域進(jìn)行，一種在F-K域進(jìn)行。T-X域變換的計算時間只有SLNFOW和SLNINV的70%左右；F-K域變換的計算速度更快，僅有SLNFOW和SLNINV的30%左右。

3.2 模塊使用說明

3.2.1 TAUP模塊參數(shù)表

3.2.2 TAUP模塊參數(shù)說明

CALL卡是指調(diào)用卡。其中TAUP是指模塊名；TYPE是指定變換所采用的算法。

FORWARD卡指定義正變換的參數(shù)。其中FORWARD指正變換卡片標(biāo)識；PMIN是指定義變換的最小P值。 DP指定義射線參數(shù)P的增 量。NP指定義變換的總射線數(shù)，如果該值小于記錄道數(shù)，則取記錄道數(shù)；FMIN指定義頻率的最小值。

MASK 卡定義TAPER道數(shù)與方式。其中MASK指MASK卡片標(biāo)識；LTAP指左邊濾波道數(shù)； LPOW指左邊濾波指數(shù)； RTAP指右邊濾波道數(shù)；RPOW指右邊濾波指數(shù)。

INVERSE卡指定義反變換的參數(shù)。其中INVERSE指反變換卡片標(biāo)識；NTR指最大記錄道數(shù)；FMIN指頻率最小值，該卡只在TYPE為FK時有效；NRHO指RHO濾波點數(shù)，該卡只在TYPE為TX時有效。

4 TAUP變換應(yīng)用效果及分析

4.1模型數(shù)據(jù)

為了驗證TAUP正反變換的正確性，用理論合成記錄做TAUP正反變換處理，從模型數(shù)據(jù)來看，變換結(jié)果是正確的，如圖3所示。

4.2實際野外數(shù)據(jù)

實際數(shù)據(jù)選取某地區(qū)86W285測線，該測線950炮，5s，2ms采樣，120道，2020個CDP。選取該測線的某個處理流程作為測試流程，該流程包含F(xiàn)OCUS版本的TAUP變換模塊SLNFOW和 SLNINV。分別用新編寫的TX域正反變換、FK域正反變換和FK域正變換TX域反變換來替換這兩個模塊，變換參數(shù)一致。在SUN工作站上用一個CPU分別運行這幾個作業(yè)，結(jié)果如表2及圖4、圖5所示。（ 單位：小時）

結(jié)束語

從上面的模型數(shù)據(jù)和實際數(shù)據(jù)應(yīng)用效果和效率我們可以得出以下結(jié)論：

１ 新編寫的快速TAUP變換模塊計算結(jié)果正確。

２ 新模塊的TX域算法比FOCUS模塊快1.3倍左右，F(xiàn)K域算法比FOCUS模塊快3倍左右。

參考文獻(xiàn)

[1]渥·伊爾瑪茲.地震數(shù)據(jù)處理[M].北京:石油出版社，1994.

[2]Ken H M，Dick P，Arthur B，et al.A comparison of three mul2

tiple─attenuation methods applied to a hard water─bottom

data set[J]，The Leading Edge，1999，18(1):120─126.

[3]牟永光.地震勘探資料數(shù)字處理方法[M].北京:石油工出版社，1981年.