孫 超

（山東省煤田地質(zhì)規(guī)劃勘察研究院，濟南，250104）

0 引言

P-SV轉(zhuǎn)換波動校正是轉(zhuǎn)換波處理的關(guān)鍵步驟，實現(xiàn)轉(zhuǎn)換波動校正直接影響到P-SV轉(zhuǎn)換波地震數(shù)據(jù)處理過程中的速度分析、轉(zhuǎn)換波疊加和轉(zhuǎn)換波偏移的效果[1]。P-SV轉(zhuǎn)換波時距方程是決定轉(zhuǎn)換波動校正質(zhì)量的重要部分，構(gòu)建高精度轉(zhuǎn)換波動校正方程引起了許多學(xué)者的關(guān)注[2-4]。

常用的3種對轉(zhuǎn)換波旅行時方程的擬合即Taylor級數(shù)展開法、雙平方根方程（DSR方程）和時移雙曲方程法。在Hake（1984）的基礎(chǔ)上，Tsvankin和Thomsen給出了橫向各向同性介質(zhì)轉(zhuǎn)換波的動校正方程，該動校正方程存在遠偏移距動校正精度低的問題[5]。Thomsen推導(dǎo)出了轉(zhuǎn)換波高階Taylor級數(shù)逼近的轉(zhuǎn)換波動校正方程，但是炮點與檢波點之間的距離和目的層深度之比被限制在1.0以內(nèi)[3]。在雙平方根方程（DSR方程）基礎(chǔ)上，Li和Yuan推導(dǎo)了轉(zhuǎn)換波四參數(shù)動校正方程，同時給出了轉(zhuǎn)換波兩參數(shù)動校正方程，可以有效提高轉(zhuǎn)換波動校正精度。劉洋推導(dǎo)了能反映傳播介質(zhì)垂向非均質(zhì)性的轉(zhuǎn)換波三參數(shù)動校正方程[4]。Ursin和Stovas總結(jié)了轉(zhuǎn)換波動校正方程，指出時移雙曲轉(zhuǎn)換波動校正方程的各參數(shù)物理意義不明確，動校正方程比較復(fù)雜，難以用于實際生產(chǎn)。目前，雙平方根方程（DSR方程）仍是轉(zhuǎn)換波動校正處理中的主要方向[6]。

雙平方根方程（DSR方程）在轉(zhuǎn)換波動校正過程中涉及到轉(zhuǎn)換波轉(zhuǎn)換點的計算。在水平層狀均勻介質(zhì)中，轉(zhuǎn)換波的射線路徑是不對稱的，即轉(zhuǎn)換點到炮點的距離不等于轉(zhuǎn)換點到檢波點之間的距離。因此，轉(zhuǎn)換波動校正處理不再針對共中心點道集，而是共轉(zhuǎn)換點道集[7]。確定P-SV轉(zhuǎn)換波轉(zhuǎn)換點位置一直是人們關(guān)注的問題，F(xiàn)romm采用垂直漸近線逼近得到轉(zhuǎn)換點位置的近似值，這一轉(zhuǎn)換點位置的確定只是與炮檢距和縱橫波速度比有關(guān)[8]。Taylor等推導(dǎo)了水平界面轉(zhuǎn)換點的精確公式[9]。傅旦丹也提出一種傾斜界面轉(zhuǎn)換點位置的迭代算法[10]。苑春方給出了水平介質(zhì)以及傾斜介質(zhì)的解析解[11-13]。馮晅等提出了一種針對多層傾斜介質(zhì)的最佳角度搜索法，該方法同時能給出最佳P-SV傳播路徑[14]。

目前，轉(zhuǎn)換波動校正方程假設(shè)炮、檢點位于同一水平基準(zhǔn)面，即滿足地表一致性假設(shè)。但是，在起伏劇烈且橫向速度變化劇烈的山區(qū)，難以滿足地表一致性假設(shè)，往往造成轉(zhuǎn)換波動校正的誤差較大。為解決這一問題，本文給出一種起伏地表的轉(zhuǎn)換波動校正方程，很好地解決了適應(yīng)地形起伏劇烈的勘探區(qū)的具體問題。

1 起伏地表條件下P-SV轉(zhuǎn)換波動校正方程的構(gòu)建

1.1 起伏地表條件下P-SV轉(zhuǎn)換波轉(zhuǎn)換點的確定

由公式（1）可得：

圖1 起伏地表下P-SV轉(zhuǎn)換波射線路徑

圖2 P-SV轉(zhuǎn)換波動校正傳播示意圖

將式（5）兩邊平方，整理后，得到關(guān)于炮點到轉(zhuǎn)換點水平距離的四次方程：

式中：

對方程（5）兩邊平方，分類移項整理后可將其寫為以下迭代形式：

1.2 起伏地表條件下P-SV轉(zhuǎn)換波動校正方程

如圖2所示，轉(zhuǎn)換波射線路徑中炮點為S，檢波點為R，縱波速度為，轉(zhuǎn)換波速度為，炮點到轉(zhuǎn)換點的水平距離是，檢波點到轉(zhuǎn)換點的水平距離是， 炮點到反射界面的高程為，檢波點的高程為，M點的高程為,反射界面的高程為。那么準(zhǔn)確的轉(zhuǎn)換波時距曲線可以表示為：

M點處，縱波的垂直旅行時為，轉(zhuǎn)換橫波旅行時為。

方程（11）為轉(zhuǎn)換波動校正方程。

2 數(shù)值模擬測試

為了檢驗本文構(gòu)建的起伏地表下P-SV轉(zhuǎn)換波動校正方程的有效性，設(shè)計了起伏地表條件地質(zhì)模型。圖3為地表起伏條件下的地質(zhì)模型，該模型存在4個波阻抗界面，且各反射界面均為水平狀況，模型參數(shù)如表1。

圖3 地層模型

表1 模型二參數(shù)表

采用射線追蹤方法計算了該模型條件下的地震記錄（圖4a），圖4b和圖4c分別為DSR方程和本文方程對地震記錄進行動校正處理的結(jié)果?？梢钥闯?，對于起伏地表條件，采用本文動校正方法可以取得較好的動校正效果，然而DSR方程動校正效果較差且存在較大的誤差。在地表一致性假設(shè)條件下，地震波在近地表為垂直入射和出射，相較于縱波出射，轉(zhuǎn)換橫波速度較小，尤其是在一些地表起伏較大且低速層速度較高的地區(qū)，因不再滿足地表一致性假設(shè)條件，導(dǎo)致誤差較大。此外，轉(zhuǎn)換點的計算是建立在水平地表這一條件下的，因此真實轉(zhuǎn)換點位置與計算獲取的轉(zhuǎn)換點位置存在誤差；然而采用新方程進行校正時，盡可能獲得準(zhǔn)確的轉(zhuǎn)換點位置，可以提高轉(zhuǎn)換波動校正的精度。

3 資料處理

圖4 轉(zhuǎn)換波動校正效果對比圖

利用河南YC煤礦三維三分量地震資料，測試轉(zhuǎn)換波動校正方法的實際應(yīng)用效果。針對轉(zhuǎn)換波資料進行預(yù)處理：靜校正、面波壓制和反褶積等，獲得較高品質(zhì)的預(yù)處理質(zhì)料，將其應(yīng)用于轉(zhuǎn)換波動校正分析。圖5a為原始共轉(zhuǎn)換點道集，反射波同相軸表現(xiàn)出類雙曲型特征。圖5b為采用常規(guī)動校正方法對共轉(zhuǎn)換點道集進行動校正的效果，結(jié)果表明動校正可以將近、中偏移距的反射波同相軸校平，具有一定的動校正效果。圖5c為采用本文提出的起伏地表下的轉(zhuǎn)換波動校正方法處理后的效果，結(jié)果表明本文提出的方法能夠有效校平共轉(zhuǎn)換點道集中反射波同相軸的連續(xù)性。通過對比圖5a～5c表明，本文動校正方程可以實現(xiàn)轉(zhuǎn)換波資料的動校正處理，其在近偏移距和遠偏移距范圍內(nèi)都具有較好的動校正效果。圖6是采用不同的動校正方程進行轉(zhuǎn)換波動校正的疊加剖面，圖6a為采用常規(guī)轉(zhuǎn)換波動校正方法處理后的疊加剖面，圖6b為本文方法處理后的疊加剖面，通過對比圖6a和6b表明，本文提出的轉(zhuǎn)換波動校正方法能夠?qū)D(zhuǎn)換波數(shù)據(jù)進行成像，增強疊加剖面的連續(xù)性，提高分析時窗內(nèi)的成像效果。

圖5 共轉(zhuǎn)換點道集動校正效果對比圖

圖6 疊加剖面效果對比圖

4 結(jié)論

轉(zhuǎn)換波動校正是轉(zhuǎn)換波資料處理中的關(guān)鍵步驟之一?；诘乇硪恢滦约僭O(shè)條件下的轉(zhuǎn)換波動校正方法，將地表假設(shè)為水平層狀。但是，在地表起伏劇烈且橫向速度變化較大的山區(qū)，難以滿足地表一致性假設(shè)，傳統(tǒng)的轉(zhuǎn)換波動校正方法往往造成轉(zhuǎn)換波動校正的誤差增加，影響了P-SV轉(zhuǎn)換波動校正的準(zhǔn)確性。本文推導(dǎo)了起伏地表條件下轉(zhuǎn)換波動校正公式，并給出起伏地表條件下轉(zhuǎn)換點求解公式。通過起伏地表地質(zhì)模型測試表明：本文P-SV轉(zhuǎn)換波動校正方法，能夠取得較好的動校正效果。實際地震數(shù)據(jù)測試結(jié)果表明，起伏地表條件下轉(zhuǎn)換波動校正方法能夠提高轉(zhuǎn)換波動校正效果。