佟愷林,李　瑞,潘樹林,吳　波

(1.“油氣藏地質(zhì)及開發(fā)工程”國家重點實驗室,四川成都610059;2.西南石油大學地球科學與技術(shù)學院,四川成都610500;3.中國石油化工股份有限公司西南油氣分公司,四川成都610041;4.中國石油化工股份有限公司多波地震技術(shù)重點實驗室,四川成都610041)

?

多波聯(lián)合的轉(zhuǎn)換波折射靜校正技術(shù)及應用

佟愷林1,李瑞1,潘樹林2,吳波3,4

(1.“油氣藏地質(zhì)及開發(fā)工程”國家重點實驗室,四川成都610059;2.西南石油大學地球科學與技術(shù)學院,四川成都610500;3.中國石油化工股份有限公司西南油氣分公司,四川成都610041;4.中國石油化工股份有限公司多波地震技術(shù)重點實驗室,四川成都610041)

轉(zhuǎn)換波地震勘探是海上和陸上多波多分量地震勘探的主要工作方法。由于轉(zhuǎn)換波具有傳播速度低、吸收衰減大、能量相對較弱等特點,因此其地震資料靜校正問題突出。常規(guī)轉(zhuǎn)換波靜校正方法在表層縱橫波速度比不詳、橫波初至不易識別的情況下,效果都不太理想。針對轉(zhuǎn)換波靜校正難點,提出了多波聯(lián)合的轉(zhuǎn)換波折射靜校正技術(shù),該技術(shù)通過P-P-P折射波初至拾取技術(shù)和P-P-SV折射波共檢波點初至疊加技術(shù)的聯(lián)合,實現(xiàn)了P-P-SV折射波初至的成像及準確識別；以此為基礎,再聯(lián)合利用P-P-P折射波折射靜校正技術(shù)和P-P-SV折射波延遲時差提取技術(shù),求取P-P-SV折射波檢波點延遲時間。最后,根據(jù)P-P-SV折射波檢波點延遲時間,求解出P-SV波檢波點基準面靜校正量。實際數(shù)據(jù)測試結(jié)果表明,多波聯(lián)合的轉(zhuǎn)換波折射靜校正技術(shù)提高了轉(zhuǎn)換波檢波點靜校正量計算精度,能夠明顯提高轉(zhuǎn)換波資料成像質(zhì)量。

轉(zhuǎn)換波;折射;靜校正;延遲時間;初至疊加

地震波在傳播過程中,轉(zhuǎn)換波(P-SV)相對縱波(P-P)具有傳播速度低、吸收衰減大、能量相對較弱、低速帶厚度不一致、傳播路徑不相同等特點,轉(zhuǎn)換波的表層靜校正量更大,變化更加復雜,導致P-SV波靜校正量通常是相同位置P-P波靜校正量的數(shù)倍,P-SV波靜校正量難以準確估計;轉(zhuǎn)換波記錄信噪比通常較低,初至附近P-P波折射、P-SV波折射等信息混雜在一起,使得P-SV波初至不易分辨,難以拾取,這些難點造成了轉(zhuǎn)換波靜校正處理問題突出。

在P-SV波靜校正處理過程中,由于炮點靜校正量是通過P-P波靜校正處理得到,因此難點是求取P-SV波檢波點靜校正量。實際資料處理多利用P-P波靜校正量中的檢波點靜校正量乘以一個系數(shù)來獲取P-SV波檢波點靜校正量。雖然此方法能夠提高轉(zhuǎn)換波資料的成像質(zhì)量,但是,由于求得的轉(zhuǎn)換波檢波點靜校正量精度不夠,會導致轉(zhuǎn)換波資料處理成果中經(jīng)常出現(xiàn)假構(gòu)造現(xiàn)象,且在低信噪比地區(qū)還可能出現(xiàn)無法成像的情況。

針對轉(zhuǎn)換波靜校正的特點,前人進行了大量研究,并提出了很多方法。楊海申等[1]提出利用轉(zhuǎn)換波與縱波初至的互相關(guān)函數(shù)疊加結(jié)果來推測轉(zhuǎn)換波初至位置的算法為轉(zhuǎn)換波初至拾取方法提供了新的思路。這種算法在轉(zhuǎn)換波初至特別清晰的時候,可以大體上確定出轉(zhuǎn)換波初至的位置,在這個大體位置附近使用手工拾取精確的轉(zhuǎn)換波初至。潘樹林等[2]總結(jié)出轉(zhuǎn)換波初至判定的標準,提出使用初至疊加法確定轉(zhuǎn)換波初至,并用τ-p變換自動拾取轉(zhuǎn)換波初至的思路。該算法提出在CRP(共檢波點)道集中使P-P-SV波(P波入射,在折射界面上以P波速度滑行,然后以SV波形式出射)同相折射初至能量通過疊加得到增強,而其它信息能量被衰減,從而分辨出P-P-SV波折射初至。

雖然轉(zhuǎn)換波靜校正方法較多[3-7],但是,目前仍然沒有一套能夠適用于實際生產(chǎn)的轉(zhuǎn)換波表層靜校正技術(shù)。基于此,結(jié)合縱波和轉(zhuǎn)換波的特點,提出了一種新的多波聯(lián)合的轉(zhuǎn)換波折射靜校正技術(shù),以求解轉(zhuǎn)換波檢波點靜校正量。該方法利用成熟的縱波折射波(P-P-P波,即P波入射,在折射界面上以P波速度滑行,然后以P波形式出射)靜校正技術(shù),求取較準確的折射層厚度和縱波檢波點延遲時,并結(jié)合共檢波點初至疊加技術(shù)和P-P-SV波延遲時差法,求解轉(zhuǎn)換波檢波點靜校正量[8-12]。該技術(shù)應用于轉(zhuǎn)換波實際資料處理,能夠提高轉(zhuǎn)換波檢波點靜校正量精度,改善轉(zhuǎn)換波資料成像質(zhì)量,構(gòu)造形態(tài)更加真實[13-19]。

1　多波聯(lián)合的轉(zhuǎn)換波折射靜校正基本方法

分析P-P-P折射波初至延遲時法靜校正技術(shù)、P-P-SV折射波共檢波點初至疊加技術(shù)以及P-P-SV折射波延遲時差法的方法原理,并以這3種方法為基礎,進行多波聯(lián)合的轉(zhuǎn)換波折射靜校正。

1.1P-P-P折射波初至延遲時法靜校正技術(shù)

延遲時是指地震波從激發(fā)點經(jīng)折射界面滑行傳播到接收點所用的時間相對于地震波從激發(fā)點直接以折射層速度傳播到接收點所用的時間差。這個延遲時可以分解為炮點延遲時和檢波點延遲時[3]。

如圖1所示,設A為炮點,B為檢波點,A點和B點的表層厚度分別為Za和Zb,臨界角為α,則A點激發(fā),B點接收到的折射初至時間為：

(1)

檢波點延遲時間為：

(2)

對工區(qū)內(nèi)有效的炮點-檢波點對,建立一個超大規(guī)模的、超定的方程組,解這個方程組,得到檢波點延遲時(dP)、檢波點處的折射層厚度(Zb,即B點的表層厚度)和縱波折射波速度(v1)。

圖1　折射模型

1.2P-P-SV折射波共檢波點初至疊加技術(shù)

對多個水平層或傾斜界面情況下的P波和P-SV波的折射波時距曲線進行推導,發(fā)現(xiàn)P-P-P折射波和P-P-SV折射波的時距曲線在同一構(gòu)造下曲線斜率相同,即在同一構(gòu)造下兩者相互平行。

基于以上發(fā)現(xiàn),利用共檢波點初至疊加技術(shù)確定縱波初至與轉(zhuǎn)換波初至時差(Δt),具體過程為：

1) 首先將轉(zhuǎn)換波記錄所對應縱波分量的P-P-P折射初至波進行準確拾取;

2) 將拾取到的初至,應用于對應的P-SV波記錄上,所有初至水平校正至某個時間點,記錄為T0;

3) 將P-SV波中的每一道使用對應的移動量進行校正,移動量為P-P-P折射初至時間和某個固定時間與T0的差值;

4) 然后將改變位置之后的P-SV波CRP道集進行疊加;

5) 拾取疊加道中能量最強的波峰所對應的時間位置作為P-P-SV折射初至與P-P-P折射初至的時間差。結(jié)果如圖2所示。

圖2　某地區(qū)轉(zhuǎn)換波初至疊加a 轉(zhuǎn)換波資料CRP道集(紅線為縱波初至,藍線為轉(zhuǎn)換波初至); b 轉(zhuǎn)換波資料CRP疊加道

1.3轉(zhuǎn)換波延遲時差法

轉(zhuǎn)換波延遲時差法[1]的原理為：轉(zhuǎn)換波檢波點延遲時等于轉(zhuǎn)換波與縱波初至時差及縱波檢波點延遲時之和。

如圖3所示,從激發(fā)點A到接收點B的縱波與轉(zhuǎn)換波初至時間差為：

圖3　縱波和轉(zhuǎn)換波初至傳播路徑

Δt=[TAC+TCE+TEB]-[TAC+TCD+TDB]

=TDE+TEB-TDB

=TEB-TEF+TEF+TDE-TDB

=(TEB-TEF)-(TDB-TDF)

(3)

dS=Δt+dP

(4)

dP=TDB-TDF

dS=TEB-TEF

(5)

式中：dP為縱波檢波點延遲時;dS為轉(zhuǎn)換波檢波點延遲時;Δt為轉(zhuǎn)換波與縱波的初至時差。

1.4多波聯(lián)合的轉(zhuǎn)換波折射靜校正技術(shù)

1.4.1技術(shù)思路

利用高質(zhì)量的P-P-P折射波初至信息獲得P-P波近地表模型,以此為基礎,獲取高質(zhì)量的P-P-SV折射波初至信息,通過建立檢波點處P-P-P折射波與P-P-SV折射波的參數(shù)轉(zhuǎn)換關(guān)系,獲得檢波點處P-SV波近地表模型,從而實現(xiàn)P-SV波基準面靜校正。

1.4.2基本原理

利用P-P-P折射波初至延遲時法靜校正技術(shù)、P-P-SV折射波共檢波點初至疊加技術(shù)以及P-P-SV折射波延遲時差法分別得到每一折射層檢波點處折射層厚度(Zb)、折射層速度(v1)、轉(zhuǎn)換波檢波點延遲時(dS)?；谶@些信息,采用多波聯(lián)合的轉(zhuǎn)換波折射靜校正技術(shù)綜合反演得到每一折射層相應的橫波速度(vS)。

(6)

(7)

式中：β為轉(zhuǎn)換波折射角。

聯(lián)合求解(6)式和(7)式,可得：

(8)

由公式(8),可知β=arctan[Zb/(dSv1)],代入(6)式即可求得vS。

當各檢波點對應折射層的厚度和橫波速度已知,利用高程靜校正技術(shù)就可以求解出各檢波點轉(zhuǎn)換波靜校正量:

(9)

式中:vc為填充速度，通常根據(jù)工區(qū)特點給定一個常速度;zbi為第i層的檢波點處折射層厚度;vSi為第i層的橫波速度;Eg為物理點高程;Ed為固定基準面;Δt為檢波點轉(zhuǎn)換波靜校正量。

1.4.3實現(xiàn)方法

方法的實現(xiàn)步驟為：

1) 精確拾取P-P波初至信息;

2) 利用P-P波初至信息,采用縱波折射波靜校正技術(shù),獲得P-P-P折射波近地表模型參數(shù),包括檢波點折射層厚度、P-P-P折射波檢波點延遲時間、P-P-P折射波折射層速度;

3) 利用P-P波初至信息,采用P-P-SV波共檢波點初至疊加技術(shù),獲得P-P-SV折射波初至與P-P-P折射波初至的時間差值;

4) 利用P-P-SV折射波初至與P-P-P折射波初至的時間差值和P-P-P折射波檢波點延遲時間,采用P-P-SV折射波延遲時差技術(shù),獲得P-P-SV折射波的檢波點延遲時間;

5) 結(jié)合檢波點折射層厚度、P-P-P折射波檢波點延遲時間、P-P-P折射波折射層速度、P-P-SV折射波的檢波點延遲時間,采用多波聯(lián)合的轉(zhuǎn)換波折射靜校正技術(shù),獲得P-SV波檢波點處的橫波速度;

6) 綜合以上P-SV波的近地表結(jié)構(gòu),利用基于近地表模型的靜校正技術(shù)可以求解出P-SV波檢波點基準面靜校正量。

2　實際資料應用

利用多波聯(lián)合的轉(zhuǎn)換波折射靜校正技術(shù)對某工區(qū)3D3C資料進行轉(zhuǎn)換波靜校正處理,該資料共10398炮,每炮12個排列,共35322個檢波點。資料信噪比較低,低降速帶橫向變化劇烈。具體處理過程如下：

1) 首先采用P-P-P折射波初至延遲時法靜校正技術(shù)對縱波資料進行處理,得到每個檢波點對應折射層的厚度、速度和每個檢波點的縱波延遲時間。

2) 采用P-P-SV折射波共檢波點初至疊加技術(shù)確定每個檢波點的轉(zhuǎn)換波初至時差,圖4為轉(zhuǎn)換波資料CRP疊加剖面轉(zhuǎn)換波初至時差拾取效果,可以看出,初至疊加效果較好,拉平后的縱波初至時間下方明顯出現(xiàn)強能量的轉(zhuǎn)換波初至,有利于初至拾取。

3) 對步驟2)得到的轉(zhuǎn)換波初至時差和縱波檢波點延遲時間,采用P-P-SV折射波延遲時差法求取轉(zhuǎn)換波檢波點延遲時間。

4) 最后采用多波聯(lián)合的轉(zhuǎn)換波折射靜校正技術(shù),結(jié)合每個檢波點的轉(zhuǎn)換波延遲時間,及其對應折射層的厚度、速度,計算出轉(zhuǎn)換波檢波點靜校正量。

圖5對比了縱波和轉(zhuǎn)換波檢波點靜校正量?？梢钥闯?計算出的轉(zhuǎn)換波檢波點靜校正量分布范圍為-19～109ms,遠大于縱波檢波點靜校正量分布范圍(-12～70ms),并且圖5 中靜校正量曲線特征與轉(zhuǎn)換波初至時差大致相當,也與近地表高程特征相符合,因此可以認為計算出的轉(zhuǎn)換波檢波點靜校正量正確。

圖6分別為轉(zhuǎn)換波資料施加轉(zhuǎn)換波初至時差后的遠、中、近偏移距炮集和轉(zhuǎn)換波CRP疊加剖面對比圖?？梢钥闯?當施加本方法計算的轉(zhuǎn)換波靜校正量后,轉(zhuǎn)換波炮集的初至更加平滑,且深層的反射波呈現(xiàn)出雙曲線特征,與P-P-SV折射波基本原理相符,說明計算出的轉(zhuǎn)換波靜校正量是正確的。

圖7a是施加縱波炮點靜校正量的疊加剖面,可看出,剖面上構(gòu)造起伏不大,但信噪比較低,成像質(zhì)量不高;圖7b是施加縱波炮點靜校正量和縱波檢波點靜校正量的疊加剖面,可看出,加入縱波檢波點靜校正量后,疊加剖面成像質(zhì)量有所提高,但改善程度不大;圖7c是施加縱波炮點靜校正量和轉(zhuǎn)換波檢波點靜校正量的疊加剖面。比較圖7c和圖7b可以看出,加入轉(zhuǎn)換波檢波點靜校正量后的疊加剖面,淺、中、深層成像質(zhì)量均得到了改善,構(gòu)造特征更加清晰,沒有出現(xiàn)串相位現(xiàn)象,說明轉(zhuǎn)換波檢波點靜校正量比縱波檢波點靜校正量準確。

圖4　轉(zhuǎn)換波資料CRP疊加剖面轉(zhuǎn)換波初至時差拾取效果(紅線為轉(zhuǎn)換波初至)a 道號14130到21194; b 道號21195到28259;

圖5　縱波和轉(zhuǎn)換波檢波點靜校正量對比a 縱波檢波點靜校正量; b 轉(zhuǎn)換波檢波點靜校正量

圖6　轉(zhuǎn)換波資料施加轉(zhuǎn)換波初至時差后的近(a)、中(b)、遠(c)偏移距炮集和轉(zhuǎn)換波CRP疊加剖面對比

圖7　施加不同靜校正量靜校正后轉(zhuǎn)換波疊加剖面對比a 縱波炮點靜校正量; b 縱波炮點靜校正量和縱波檢波點靜校正量; c 縱波炮點靜校正量和轉(zhuǎn)換波檢波點靜校正量

3　結(jié)束語

針對轉(zhuǎn)換波靜校正處理中轉(zhuǎn)換波資料信噪比低、轉(zhuǎn)換波初至不易分辨和轉(zhuǎn)換波檢波點靜校正量大等難點,充分利用縱波折射延遲時靜校正、疊加法確定轉(zhuǎn)換波初至時差和轉(zhuǎn)換波延遲時差提取三項技術(shù)的優(yōu)勢,建立了一種新的適用于實際資料處理的多波聯(lián)合的轉(zhuǎn)換波折射靜校正技術(shù)。實際地震資料處理效果表明：轉(zhuǎn)換波靜校正問題得到了很好的解決,單炮記錄上轉(zhuǎn)換波同相軸連續(xù)性增強;疊加剖面上消除了串相位現(xiàn)象、信噪比明顯提高、構(gòu)造形態(tài)更加合理、成像效果顯著改善。表明多波聯(lián)合轉(zhuǎn)換波折射靜校正技術(shù)具有很好的應用前景。

[1]楊海申,李彥鵬,陳海青.轉(zhuǎn)換波延遲時靜校正[J].石油地球物理勘探,2006,41(1):13-16

YANG H S,LI Y P,CHEN H Q.P-SV converted wave static corrections by delay times[J].Oil Geophysical Prospecting,2006,41(1):13-16

[2]潘樹林,高磊,周熙襄,等.一種改進的P-SV轉(zhuǎn)換波靜校正方法[J].石油物探,2007,46(2):143-146

PAN S L,GAO L,ZHOU X X,et al.An improved static correction method of converted P-SV wave[J].Geophysical Prospecting for Petroleum,2007,46(2):143-146

[3]MIKE C.反射波地震勘探靜校正技術(shù)[M].李培明,柯本善,譯.北京:石油工業(yè)出版社,2004:1-358

MIKE C.Static corrections for seismic reflection surveys[M].LI P M,KE B S,translators.Beijing:Petroleum Industry Press,2004:1-358

[4]CARY P W,EATON W S.A simple method for resolving large converted wave (P-SV) statics[J].Geophysics,1993,58(3):429-433

[5]馬昭軍,唐建明.基于構(gòu)造時間控制的P-SV波靜校正方法[J].物探化探計算技術(shù),2008,30(5):373-376

MA Z J,TANG J M.A P-SV wave static correction method based on controlling structure time[J].Computing Techniques for Geophysical and Geochemical Exploration,2008,30(5):373-376

[6]RICHARD R,VAN D I.Static correction for PS-wave surface seismic surveys,recent advances in shear wave technology for reservoir characterization[J].SEG/EAGE Summer Research Workshop,2000：233-235

[7]FJELLANGER J P,BOEN F,RONNING K J.Successful use of converted wave data for interpretation and well optimization on Grane[J].Expanded Abstracts of 76thAnnual Internat SEG Mtg,2006:1138-1142

[8]EDME P,KRAGH E.Near-surface S-wave velocity estimation from P-wave polarization analysis[J].Expanded Abstracts of 79thAnnual Internat SEG Mtg,2009:4289-4292

[9]DAVID W E,DON C L.P-SV stacking charts and binning periodicity[J].Geophysics,1992,57(5):745-748

[10]ROTHMAN D H.Nonlinear inversion,statistical mechanics,and residual statics estimation[J].Geophysics,1985,50(12):2784-2796

[11]RONEN J,CLAERBOUT J F.Surface-consistent residual statics estimation by stack power maximization[J].Geophysics,1985,50(12):2759-2767

[12]STEWART R R,GAISER J E,Brown R J,et al.

Converted-wave seismic exploration:application[J].Geophysics,2003,68(1):40-57

[13]桑運云,孫軍曉,焦淑萍,等.起伏地表下基于拋物插值的最短路徑射線追蹤[J].石油物探,2014,53(2):142-148

SANG Y Y,SUN J X,JIAO Y P,et al.Shortest path raytracing based on parabolic traveltime interpolation in irregular topography[J].Geophysical Prospecting for Petroleum,2014,53(2):142-148

[14]陳海峰,李向陽,趙波,等.三維轉(zhuǎn)換波處理技術(shù)研究及應用[J].石油物探,2013,52(5):512-518

CHEN H F,LI X Y,ZHAO B,et al.3D converted wave seismic data processing technology and its application[J].Geophysical Prospecting for Petroleum,2013,52(5):512-518

[15]劉財,楊慶節(jié),耿美霞.多波多分量地震資料速度掃描法轉(zhuǎn)換波識別技術(shù)[J].石油物探,2013,52(4):383-387

LIU C,YANG Q J,GENG M X.P-SV wave identification technique by velocity scanning for multiwave multi-component seismic data[J].Geophysical Prospecting for Petroleum,2013,52(4):383-387

[16]劉玉萍,潘衛(wèi)國,丁飛浩,等.轉(zhuǎn)換波組合靜校正方法應用及效果[J].地球物理學進展,2014,29(2):805-814

LIU Y P,PAN W G,DING F H,et al.Application and effect of combined static correction method of converted wave[J].Progress in Geophysics,2014,29(2):805-814

[17]查文鋒,于小磊.淺表層結(jié)構(gòu)復雜區(qū)靜校正的優(yōu)化選取及應用[J].地球物理學進展,2014,29(2):958-965

CHA W F,YU X L.The optimal selection and comprehensive application of the static correction in the complex shallow surface structure areas[J].Progress in Geophysics,2014,29(2):958-965

[18]陳海峰,李向陽,錢忠平,等.聯(lián)合縱波確定轉(zhuǎn)換波接收點大的靜校正量[J].地球物理學進展,2014,29(1):306-311

CHEN H F,LI X Y,QIAN Z P,et al.Estimating C-wave receiver large statics associated with p-wave imaging[J].Progress in Geophysics,2014,29(1):306-311

[19]張正峰,郭建,韓立國,等.基于波場延拓的變換坐標轉(zhuǎn)換波靜校正[J].地球物理學報,2013,56(5):1676-1685

ZHANG Z F,GUO J,HAN L G,et al.Transformation coordination converted wave static correction based on wave field continuation[J].Chinese Journal of Geophysics,2013,56(5):1676-1685

(編輯：朱文杰)

Refraction static correction technique for multi-wave joint converted wave and its application

TONG Kailin1,LI Rui1,PAN Shulin2,WU Bo3,4

(1.StateKeyLaboratoryofOilandGasReservoirGeologyandExploitation,ChengduUniversityofTechnology,Chengdu610059,China;2.SchoolofGeoscienceandTechnology,SouthwestPetroleumUniversity,Chengdu610500,China;3.SinopecSouthwestCompany,Chengdu610041,China;4.SinopecKeyLaboratoryofMultiComponentsSeismicTechnology,Chengdu610041,China)

The converted wave seismic exploration is the main method of multi-wave and multi-component seismic exploration applied to onshore and offshore exploration.However,the converted wave is characterized by low velocity,serious absorption and attenuation,weak energy and so on,which highlights static correction of converted wave.Research of converted wave static correction undertaken by many researchers is of little help,when surface P-velocity/S-velocity is unknown and first break of S-wave is difficult to identify.We propose a refraction static correction technique for multi-wave combined converted wave aimed at the difficulty.This technique achieved P-P-SV first break imaging and identification by the combination of the first break pickup for P-P-P refraction wave and the first break stacking in CRP gather for P-P-SV refraction wave.On this basis,to solve the delay time of receiver point for P-P-SV refraction wave we combined the refraction static correction of P-P-P refraction wave with the delay time extraction of P-P-SV refraction wave.At last,according to the delay times of receiver point for P-P-SV refraction wave,the datum level statics of P-SV receiver point can be obtained.The experimental result shows that the refraction static correction technique of multi-wave combined converted wave improved the calculation accuracy of converted wave static correction and the imaging quality.

converted wave,refraction,static correction,delay time,first break stacking

2016-04-15;改回日期：2016-06-07。

佟愷林(1987—),男,博士,主要從事地震資料處理及解釋等工作。

國家自然科學基金項目(41204101)和四川省重點實驗室開放基金項目(2015trqdz03)聯(lián)合資助。

P631

A

1000-1441(2016)05-0674-08

10.3969/j.issn.1000-1441.2016.05.006

This research is financially supported by the National Natural Science Foundation of China (Grant No.41204101) and the Opening Foundation for Sichuan Provincial Key Laboratory (Grant No.2015trqdz03).