張正峰，郭 建，，韓立國，徐善輝，鞏向博

1吉林大學 地球探測科學與技術(shù)學院，長春 130026 2中國科學院地質(zhì)與地球物理研究所，北京 100029

1 引 言

P－SV地震勘探技術(shù)在氣云下成像、PP反射差地區(qū)目標成像、碎屑巖巖性構(gòu)圖及提高淺層分辨率（小于1000m深度）等方面已證實有效，并且在各向異性參數(shù)求取、彈性反演方面有著重要前景［1－4］，尤其是在裂縫檢測［5－6］、儲層預測［7］和探測隱伏活動斷裂［8］方面已取得較好成果.隨著三維地震勘探［9－10］在很多地區(qū)的成功實踐應用，P－SV靜校正將成為一個很有挑戰(zhàn)性的課題，它遠比P波復雜，是因為近地表低速帶的厚度、速度和密度的橫向變化對P與S來說是不同的，在同一位置上S波靜校正的數(shù)值通常比P波的大2～10倍［11］.因此，P－SV波靜校正的難點在于接收點的S波靜校正.P－SV波靜校正可以分為反射波法和折射波法兩大類，對于反射波法，比較有代表性的有：P.W.Cary［12］提出的共檢波點疊加道相關(guān)法，主要用來求取短波長靜校正量，郭桂紅［13］用曲線擬合圓滑消除構(gòu)造影響、潘樹林［14］用構(gòu)造校正分別對其進行了改進，隨后潘樹林［15］用初至疊加法拾取初至，根據(jù)縱波和轉(zhuǎn)換波的初至時差計算了轉(zhuǎn)換波的靜校正量，并在宣漢地區(qū)得到很好效果；對于折射波法，比較有代表性的有：Armin W Schafter［16］提出的轉(zhuǎn)換波折射靜校正法，主要用來求取長波長靜校正量，李彥鵬［17］提出的轉(zhuǎn)換波延遲時靜校正法，它根據(jù)P波和S波的初至時差來求取S波靜校正，Heelan［18］很好的計算了折射波的振幅強度，為轉(zhuǎn)換波靜校正提供了便利，LIU［19］對折射波法進行了深入研究，得出大波阻抗差是產(chǎn)生折射橫波的條件，李彥鵬［20］將轉(zhuǎn)換波延遲時和共檢波點疊加相關(guān)結(jié)合，得出一個針對厚風化層覆蓋區(qū)的有效轉(zhuǎn)換波靜校正方法.這些方法在一定條件下均取得了不錯的效果，為研究轉(zhuǎn)換波靜校正提供了理論基礎(chǔ).

頻率波數(shù)域偏移是一種效率高、效益好的方法，Dan Kosloff［21］、Bengt Fornberg［22］等人對其做了很多改進，John R.Berryhill［23］最早指出可以將其應用到P波靜校正中，隨后很多學者在這方面做了工作，Reshef［24］用相移法以逐步累加向下波場延拓的方式完成了起伏地表的基準面校正，Bevc D［25］用克希霍夫算子以波場向上外推的方式完成了起伏地表的基準面校正，王祥春［26］通過對起伏地表進行變換，用雙平方根法將地震波場向下延拓完成了起伏地表的靜校正，楊鍇［27］用有限差分法以逐步累加向下延拓實現(xiàn)了非水平觀測面的基準面校正，隨后又將其應用到三維波動方程的基準面校正中［28］，李錄明［29］做了基于復雜表層模型的三維波場延拓研究，盧回憶［30］用快速Fourier變換相移、分裂步和一階退化三種不同的算子完成了波動方程的基準面校正，王紅落［31］研究分析了波動方程基準面延拓的進展，陳樹文［32］研究了復雜地質(zhì)體的波動方程基準面技術(shù)，王建民［33］研究了波場延拓的速度分析，李翔［34］用變換坐標系的方法推導出一種新的延拓算子，用近地表速度向上延拓完成了復雜地表的靜校正，可見，頻率波數(shù)域波場延拓在P波靜校正中已取得良好效果，Li［35］提出可以將其應用到轉(zhuǎn)換波靜校正中.另外，許多學者在轉(zhuǎn)換波的其他方面也做了深入研究，崔杰［36］提出了高分辨率的P－SV波速度分析方法，馮晅［37－38］分別研究了P－SV波轉(zhuǎn)換點的最佳角度搜索方法和P－SV波的AVO反演，李天成［39］研究了P－SV波轉(zhuǎn)換點和各向異性參數(shù)的關(guān)系，蔡曉剛［40］研究了P－SV波的反投射系數(shù)，陳學良［41］研究了場地鉆井臺陣觀測中P－SV波的分離方法，蘆俊［42］設(shè)計了P－SV波的采集觀測系統(tǒng)，這都為轉(zhuǎn)換波的靜校正提供了理論幫助.本文在此基礎(chǔ)上，提出基于波場延拓的變換坐標轉(zhuǎn)換波靜校正，首先對起伏地表函數(shù)進行坐標變換，將炮點和檢波點映射到同一水平面上，使地震數(shù)據(jù)的頻率波數(shù)域和時間空間域能夠相對應；接著波場分離，提取出P－SV地震數(shù)據(jù)，將P－SV地震數(shù)據(jù)分成共炮點道集和共檢波點道集；然后根據(jù)新坐標系下的波動方程延拓公式向上延拓P－SV地震數(shù)據(jù)到基準面，共炮點道集用近地表S波速度，共檢波點用P波速度，恢復起伏地表到基準面之間的真實波場；最后還原到原始坐標系，取出基準面P－SV地震數(shù)據(jù)并完成校正.該方法不但適用于起伏平緩地區(qū)，而且可以在地形復雜地區(qū)應用.

2 方法要點

2.1 變換坐標轉(zhuǎn)換波波場延拓靜校正

變換坐標轉(zhuǎn)換波波長延拓靜校正方法可簡單分為四步：

1）坐標變換，將炮點和檢波點映射到同一水平面上；

2）波場分離，提取P－SV地震數(shù)據(jù)，P－SV地震數(shù)據(jù)分成共炮點道集和共檢波點道集；

3）根據(jù)新坐標系下的波動方程延拓公式向上延拓P－SV地震數(shù)據(jù)到基準面，共炮點道集用近地表S波速度，共檢波點用P波速度；

4）將延拓到基準面的波場還原到原始坐標系，即是校正到基準面的P－SV地震數(shù)據(jù).

圖1 基于波場延拓的變換坐標轉(zhuǎn)換波靜校正原理圖（a）原始坐標下炮點和檢波點空間位置；（b）變換坐標下炮點和檢波點空間位置；（c）變換坐標系下延拓炮點和檢波點；（d）原始坐標系下延拓炮點和檢波點；（e）在原始坐標系下選取浮動水平基準面；（f）完成轉(zhuǎn)換波靜校正.Fig.1 The schematic diagram of transformating coordination converted wave static correction based on wavefiled continuation

2.2 公式推導

公式推導的目的是為了將炮點和檢波點映射到水平面上，解決延拓靜校正頻率波數(shù)域和時間空間域相對應的問題.假設(shè)地表函數(shù)為f（x），f（x）可以通過拉格朗日插值函數(shù)求取，由二維波動方程出發(fā)，有

式（1）中，x和z為二維坐標，v是速度，u是波場.

新舊坐標轉(zhuǎn)換關(guān)系式為

式（2）中，x′和z′為新二維坐標，將（2）式根據(jù)波動方程（1）進行轉(zhuǎn)變，可得

公式（4）稱為變換坐標系下的波動方程，其中

對公式（4）兩邊做二維傅里葉變換，得到變換坐標系下的頻率波數(shù)域波動方程

為解決波場延拓問題，需對公式（5）進行求解.邊界條件是

式（6）中u0（x′，z′＝0，t′）為已知變換坐標系下的地震數(shù)據(jù).

對波場函數(shù)U（kx′，z′，ω）進行分離變量，有

式（8）為針對z的二階線性常微分方程，通過求解，可以得到最后的延拓公式

公式（10）即為轉(zhuǎn)換波頻率波數(shù)域波動方程延拓靜校正的公式，它根據(jù)地表起伏函數(shù)，通過變換坐標將炮點和檢波點映射到水平面上，解決了地震數(shù)據(jù)的頻率波數(shù)域和時間空間域的相對應問題.在延拓過程中不僅易于選擇基準面，而且還能真實的還原波場，能夠得到準確的靜校正結(jié)果.

3 模擬數(shù)據(jù)處理

圖2是1000m×1000m大小的模型數(shù)據(jù)，在400m和700m處各有一反射層，每層縱波和橫波速度分別是：1500m／s和900m／s、2000m／s和1150m／s、2500m／s和1450m／s.圖3是地表高程，最大高程和最小高程相差60m.

圖4是模型正演結(jié)果.通過Tesseral軟件模擬得到，中間放炮，99道接收，10m道間距，0.3ms采樣率，1.2s總采樣時間，共4000個采樣點.從圖4中可以看出，由于地形起伏較劇烈，造成有些地方直達波也接收不到，能量斷斷續(xù)續(xù)，但地下500ms和750ms處的兩條P－P波看的還是很清楚，700ms和1000ms處的兩條P－S波隱隱可見，接近900ms處的是一條多次波.

圖2 正演模型Fig.2 The forward model

圖5是通過波場分離得到的P－SV地震數(shù)據(jù)，采用的方法是τ－p變換法［43］.圖6是地震數(shù)據(jù)的延拓過程，根據(jù)高程最大差60m，可以選擇延拓60步，每一步1m.圖7是本方法靜校正結(jié)果，根據(jù)最大高程－70m，把最佳基準面選定在－69m水平位置，從圖中可以看出校正結(jié)果呈現(xiàn)偽雙曲線的形狀，符合轉(zhuǎn)換波理論特征.

為了驗證本方法的可靠性和精確性，我們在延拓過程中隨機抽取了第40道的單道數(shù)據(jù)，藍色和紅0.12，這對正演子波的平均振幅來說是非常低的，達到了可以忽略的程度.這表明了本方法的可行性.

4 實際數(shù)據(jù)處理

為了驗證本方法的實用性，本文對東北某油田的一炮P－S波原始數(shù)據(jù)做了簡單處理，該地區(qū)地表有起伏，低降速帶S波變化劇烈，并且有較強的面波干擾.可以利用初至波自動拾取技術(shù)，通過共炮檢距校正統(tǒng)計方法，獲取近地表P波和SV波的速度.觀測系統(tǒng)道間距為10m，共450道接收，圖10為其中S1測線上的第20炮記錄，從原始記錄中的S波反射不夠光滑以及地下S波反射曲線連續(xù)不明顯可以看出地形起伏對地震記錄是有一定影響的.

圖11為通過本文方法靜校正的結(jié)果，延拓填充速度為真實地表速度.從圖中可以看出，本方法校正的S波同相軸很明顯，線條圓潤光滑，連續(xù)性很好，信噪比也很高，并且把350道附近有明顯凹陷的P波直達波也較平.圖12為隨機抽取的第200道的原始單道數(shù)據(jù)（藍色）和延拓20步（紅色）后的記錄比較，從圖中可以得知在延拓過程中地震子波的振幅（A）、頻率、相位等屬性基本沒有改變，能量也沒有損失.圖13為轉(zhuǎn)換波靜校正前的疊加數(shù)據(jù)，圖14為轉(zhuǎn)換波靜校正后的疊加數(shù)據(jù)，通過對比我們可以發(fā)現(xiàn)，校正后的反射層更清晰，連續(xù)性更好，同相軸更光滑，尤其是100道左右中間層的位置，同相軸明顯的顯現(xiàn)了出來.此結(jié)果充分驗證了變換坐標系下頻率波數(shù)域波動方程延拓靜校正的真實性和可行性.

5 認識與結(jié)論

頻率波數(shù)域波場延拓靜校正已在縱波靜校正中取得了一定成果，主要針對地形起伏劇烈，地表條件復雜的地區(qū)，這恰好為轉(zhuǎn)換波靜校正提供了思路.本文根據(jù)此思路研制新方法，通過坐標變換解決了起伏地表地震數(shù)據(jù)的頻率波數(shù)域和時間空間域的相對應問題，在延拓過程中使用真實的近地表P波和S波速度，恢復了真實的波場，保證了轉(zhuǎn)換波靜校正的精度，同時解決了基于地表一致性假設(shè)帶來的校正畸變問題.雖然將頻率波數(shù)域波場延拓應用到轉(zhuǎn)換波靜校正中，還存在不少問題尚待進一步解決，但仍可為多波多分量后續(xù)處理研究提供思路和方法.

致 謝 感謝北京吉奧菲斯科技有限責任公司的大力支持和李翔博士富有意義的建議.

（

）

［1］ 何樵登.地震勘探原理和方法.北京：地質(zhì)出版社，1986.

He Q D.The Principles and Methods of the Seismic Exploration.Beijing：Geological Publishing House（in Chinese），1986.

［2］ 黃緒德，楊文霞.轉(zhuǎn)換波地震勘探.北京：石油工業(yè)出版社，2008.

Huang X D，Yang W X.The Converted Wave Seismic Exploration.Beijing：The Petroleum Industry Press（in Chinese），2008.

［3］ Luynn H，Stewart R R，Garotta R，et al.4C－in the future－a word from the“gurus”.TheLeadingEdge，2001，20（9）：978.

［4］ Stewart R R，Gaiser J E，Brown R J，et al.Converted－wave seismic exploration：applications.Geophysics，2003，68（1）：40－57.

［5］ 程冰潔，徐天吉.轉(zhuǎn)換波方位各向異性裂縫檢測技術(shù)研究及應用.地球物理學進展，2012，27（2）：575－581.

Cheng B J，Xu T J.Research and application of fracture detection technique using azimuthal anisotropy of C－wave.ProgressinGeophys.（in Chinese），2012，27（2）：575－581.

［6］ 楊曉，王真理，喻岳鈺.裂縫型儲層地震檢測方法綜述.地球物理學進展，2010，25（5）：1785－1794.

Yang X，Wang Z L，Yu Y Y.The overview of seismic techniques in prediction of fracture reservoir.Progressin Geophys.（in Chinese），2010，25（5）：1785－1794.

［7］ 單剛義，韓立國，張麗華等.轉(zhuǎn)換波在儲層預測中的研究分析.地球物理學進展，2010，25（1）：282－287.Shang G Y，Han L G，Zhang L H，et al.Research and analysis of converted wave in reservoir prediction.Progress inGeophys.（in Chinese），2010，25（1）：282－287.

［8］ 劉保金，趙成彬，酆少英等.應用三分量淺層地震反射方法探測隱伏活動斷裂.地球物理學報，2012，55（8）：2676－2686.

Liu B J，Zhao C B，F(xiàn)eng S Y，et al.Application of the threecomponent shallow seismic reflection method to probing buried active faults.ChineseJ.Geophys.（in Chinese），2012，55（8）：2676－2686.

［9］ 付雷，劉財，王建民等.大慶地區(qū)P－SV波地震資料處理技術(shù).地球物理學進展，2004，19（2）：414－419.

Fu L，Liu C，Wang J M，et al.Seismic processing of P－SV wave in Daqing area.ProgressinGeophys.（in Chinese），2004，19（2）：414－419.

［10］ 耿麗娟.新疆哈密戈壁灘區(qū)煤層三維地震勘探實踐.地球物理學進展，2005，20（2）：393－398.

Geng L J.3Dseismic wave survey for coal layer in Hami desert area of Xinjiang.ProgressinGeophys.（in Chinese），2005，20（2）：393－398.

［11］ Tatham R H，McCormack M D，Neitzel E B，et al.Multicomponent Seismology in Petroleum Exploration.Tulsa：Soc.Expl.Geophys.，1991.

［12］ Cary P W，Eaton D W S.A simple method for resolving large converted－wave（P－SV）static.Geophysics，1993，58（3）：429－433.

［13］ 郭桂紅，王德利，何樵登等.轉(zhuǎn)換波的靜校正.吉林大學學報（地球科學版），2003，33（4）：542－544，554.

Guo G H，Wang D L，He Q D，et al.Statics of converted wave.TheJournalofJilinUniversity（EarthSciences Edition）（in Chinese），2003，33（4）：542－544，554.

［14］ 潘樹林.P－SV轉(zhuǎn)換波靜校正方法研究［博士論文］.成都：成都理工大學，2008.

Pan S L.The research for P－SV converted－wave static correction［Ph.D.thesis］（in Chinese）.Chengdu：Chengdu University of Technology，2008.

［15］ 潘樹林，高 磊，吳 波等.初至疊加法在宣漢地區(qū)轉(zhuǎn)換波靜校正中的應用.地球物理學進展，2010，25（3）：932－938.

Pal S L，Gao L，Wu B，et al.Application on the superposition method of the first arrival in P－SV static correction in the Xuanhan area.ProgressinGeophys.（in Chinese），2010，25（3）：932－938.

［16］ Schafer A W.The determination of converted－wave statics using P refractions together with SV refractions.CREWS Research Report，1991，13：51－66.

［17］ Li Y P.A new method for converted wave statics correction.72nd Annual International Meeting，SEG，Expanded Abstract，2002：979－981.

［18］ Hall D.Converted waves in refraction surveys over markers with variable depth.Geophysics，1964，29（5）：733－744.

［19］ Liu Y，Wei X C.Propagation characteristics of converted refracted wave and its application in static correction of converted wave.ScienceinChinaSeriesD：EarthSciences，2008，51（S2）：226－232.

［20］ 李彥鵬，馬在田，孫鵬遠等.厚風化層覆蓋區(qū)轉(zhuǎn)換波靜校正方法.地球物理學報，2012，55（2）：614－621.

Li Y P，Ma Z T，Sun P Y，et al.Converted－wave static correction method for thick weathering area.ChineseJ.Geophys.（in Chinese），2012，55（2）：614－621.

［21］ Dan K，John S，Zvi K，et al.Velocity and interface depth determination by tomography of depth migrated gathers.Geophysics，1996，61（5）：1511－1523.

［22］ Bengt F.The pseudospectral method：Accurate representation of interfaces in elastic wave calculations.Geophysics，1998，53（5）：625－637.

［23］ Berryhill J R.Wave－equation datuming.Geophysics，1979，44（8）：1329－1344.

［24］ Reshef M.Depth migration from irregular surfaces with depth extrapolation methods.Geophysics，1991，56（1）：119－122.

［25］ Bevc D.Flooding the topography：Wave－equation datuming of land data with rugged acquisition topography.Geophysics，1997，62（5）：1558－1569.

［26］ 王祥春，劉學偉.變換坐標系下相移法起伏地表地震波場延拓.地球物理學進展，2005，20（3）：677－680.

Wang X C，Liu X W.Downward continuing the seismic record of topography using coordination transformated method.ProgressinGeophys.（in Chinese），2005，20（3）：677－680.

［27］ 楊鍇，王華忠，程玖兵等.非水平觀測面有限差分法疊前波動方程基準面校正.石油地球物理勘探，2002，37（2）：154－162.

Yang K，Wang H Z，Cheng J B，et al.Prestack wave equation datum correction from non－flat surface with finitedifference scheme.OilGeophysicalProspecting（in Chinese），2002，37（2）：154－162.

［28］ 楊鍇，程玖兵，劉玉柱等.三維波動方程基準面校正方法的應用研究.地球物理學報，2007，50（4）：1232－1240.

Yang K，Cheng J B，Liu Y Z，et al.A study on the application of the 3－D wave－equation－datuming.ChineseJ.Geophys.（in Chinese），2007，50（4）：1232－1240.

［29］ 李錄明，羅省賢.三維波場延拓復雜表層模型校正方法及應用.地球物理學進展，2005，20（4）：1027－1034.

Li L M，Luo S X.The Method and application of 3Dwave field continuation correction based on complicated surface layer model.ProgressinGeophys.（in Chinese），2005，20（4）：1027－1034.

［30］ 盧回憶，符力耘，蔣韜.快速Fourier變換波動方程基準面校正方法研究.地球物理學進展，2010，25（4）：1313－1322.

Lu H Y，F(xiàn)u L Y，Jiang T.Wave－equation datuming based on fast fourier transform.ProgressinGeophys.（in Chinese），2010，25（4）：1313－1322.

［31］ 王紅落，陳中州，常旭.波動方程基準面延拓研究進展.地球物理學進展，2004，19（2）：230－234.

Wang H L，Cheng Z Z，Chang X.The recent progress in wave－equation datuming.ProgressinGeophys.（in Chinese），2004，19（2）：230－234.

［32］ 陳樹文，劉洪.可用于復雜地質(zhì)體的波動方程基準面技術(shù).地球物理學進展，2002，17（3）：365－369.

Cheng S W，Liu H.The usability of the wave－equation datuming for complicated geological body.Progressin Geophys.（in Chinese），2002，17（3）：365－369.

［33］ 王建民，劉財，楊寶俊等.地質(zhì)雷達技術(shù)在巴彥查干地區(qū)地震勘探中的應用.地球物理學進展，2002，17（4）：745－752.

Wang J M，Liu C，Yang B J，et al.The application of ground penetrating radar for seismic exploration in Bayanchagan area.ProgressinGeophys.（in Chinese），2002，17（4）：745－752.

［34］ 李翔.基于復雜地表波場延拓的靜校正方法研究［碩士論文］.長春：吉林大學，2009.

Li X.Static correction based on wavefield continuation in complex area［Master′s thesis］（in Chinese）.Changchun：Jilin University，2009.

［35］ Li G F，Peng S P.Static corrections for low S／N ratio converted－wave seismic data.AppliedGeophysics，2008，5（1）：44－49.

［36］ 崔杰，韓立國，廉玉廣等.高分辨率P－SV波速度分析.地球物理學進展，2009，24（6）：2036－2043.

Cui J，Han L G，Lian Y G，et al.High resolution velocity analysis of P－SV wave.ProgressinGeophys.（in Chinese），2009，24（6）：2036－2043.

［37］ 馮晅，劉財，郭繼茹等.多層傾斜層狀介質(zhì)P－SV轉(zhuǎn)換波轉(zhuǎn)換點計算的最佳角度搜索法.地球物理學報，2011，54（2）：320－328.

Feng X，Liu C，Guo J R，et al.Estimation of P－SV conversion point by searching optimal angle in a dipping or horizontal layered medium.ChineseJ.Geophys.（in Chinese），2011，54（2）：320－328.

［38］ 馮晅，舒夢珵，王兆國等.P－P波和P－SV波聯(lián)合AVO反演研究.地球物理學進展，2008，23（5）：1556－1559.

Feng X，Shu M C，Wang Z G，et al.Joint AVO inversion using P－P and P－SV seismic data.ProgressinGeophys.（in Chinese），2008，23（5）：1556－1559.

［39］ 李天成，牛濱華，孫春巖等.VTI介質(zhì)中P－SV波轉(zhuǎn)換點與各向異性參數(shù)關(guān)系.地球物理學進展，2007，22（5）：1522－1526.

Li T C，Niu B H，Sun C Y，et al.The relationship of VTI media P－SV wave converted Point and anisotropy parameters.ProgressinGeophys.（in Chinese），2007，22（5）：1522－1526.

［40］ 蔡曉剛，陳曉非.傾斜裂隙介質(zhì)反透射系數(shù)研究.地球物理學報，2009，52（5）：1253－1262.

Cai X G，Chen X F.Study on the reflection－transmission coefficients of elastic waves in the media with dipping fractures.ChineseJ.Geophys.（in Chinese），2009，52（5）：1253－1262.

［41］ 陳學良，金星.場地鉆井臺陣觀測中斜入射P波和SV波的近似剝離方法.地球物理學進展，2009，24（5）：1636－1645.

Cheng X L，Jin X.Approximate decomposition method of oblique incidence P，SV wave fields in observation of site borehole array.ProgressinGeophys.（in Chinese），2009，24（5）：1636－1645.

［42］ 蘆俊，王赟，孟召平等.優(yōu)化P－P、P－SV波聯(lián)合采集觀測系統(tǒng)設(shè)計方法.地球物理學進展，2003，18（4）：715－723.

Lu J，Wang Y，Meng Z P，et al.Method of optimizing joint P－P、P－SV acquisition geometry design.ProgressinGeophys.（in Chinese），2003，18（4）：715－723.

［43］ 馮晅，張先武，劉財?shù)?帶有多道相關(guān)的拋物線Radon變換法分離P－P、P－SV波.地球物理學報，2011，54（2）：304－309.

Feng X，Zhang X W，Liu C，et al.Separating P－P and P－SV wave by parabolic Radon transform with multiple coherence.ChineseJ.Geophys.（in Chinese），2011，54（2）：304－309.