基于Alford旋轉的轉換波各向異性校正技術

劉紅愛 馬昭軍 甘其剛

1.成都理工大學地球物理學院 2.中國石化西南油氣分公司勘探開發研究院德陽分院

基于Alford旋轉的轉換波各向異性校正技術

劉紅愛1馬昭軍2甘其剛2

1.成都理工大學地球物理學院 2.中國石化西南油氣分公司勘探開發研究院德陽分院

橫波穿過裂縫時會分裂為沿著裂縫的快橫波和垂直裂縫的慢橫波,在轉換波處理時為了提高徑向分量的成像質量,通常需要進行橫波分裂方位各向異性校正處理。為此,通過構建互相垂直的徑向和橫向分量對,實現了A lfo rd旋轉方法對轉換波數據的快慢橫波分離,并完成了轉換橫波分裂方位各向異性校正處理。首先,在淺層時窗內根據快、慢橫波分離方法得到裂縫方向以及快、慢橫波,并求取快、慢橫波之間的時延;然后,將時延應用于慢橫波作時延補償;最后,將快橫波和時延補償后的慢橫波旋轉回原來的徑向和橫向方向,得到新的徑向和橫向分量。依照該方法對川西地區地震勘探資料中深層或深層時窗內的轉換波數據進行處理,大大改善了轉換波徑向分量的成像質量。

轉換波 各向異性 A lfo rd旋轉 快波 慢波 四川盆地 西 地震資料 成像質量

橫波不平行也不垂直穿過裂縫時會發生分裂,質點振動沿裂縫走向時,傳播速度快,稱為快橫波,質點振動垂直裂縫走向時傳播速度慢,稱為慢橫波[1]。橫波與裂縫方位呈一定的角度時才會發生分裂,分裂的快慢橫波的強弱與裂縫的強度密切相關。當進行縱波源激發的轉換波勘探時,下行縱波遇到分界面時會轉換成橫波,橫波在向上傳播過程中遇到裂縫時,同樣會發生分裂。因此,兩個橫波分量檢波器接收到的信號中既有快橫波也有慢橫波信息。在三維轉換波處理時,需要將兩個橫波分量旋轉到炮點和檢波點方向(徑向分量)和垂直于炮點和檢波點方向(橫向分量)。在寬方位轉換波數據上,受到快慢橫波的影響,徑向分量的方位道集同相軸呈現“正弦”或者“余弦”狀并表現出較大的方位時差,而橫向分量方位道集每間隔90°會發生極性反轉[2]。在寬方位轉換波數據的成像處理時,需要將橫向分量的能量校正到徑向分量上以消除徑向分量的方位時差并提高成像質量[3]。

筆者基于A lford旋轉方法,實現了寬方位轉換波的快慢橫波分離,并形成了方位各向異性校正技術。根據四川盆地西部地區的寬方位三維三分量資料,完成了轉換橫波分裂方位各向異性校正處理。

1 方法原理

1.1 快慢橫波分離技術

構造接收信號矩陣:

式中S1為快橫波;S2為慢橫波。

對式(3)正交旋轉變換后,可得S=RTVR。

在實際數據處理中,分別記S的對角元素為S1、S12、S21和 S2,在轉換橫波分裂分析處理中,根據 A lford旋轉理論S的對角元素展開式為:

反旋轉后得到 R1、T1和 R2、T2分別是:

根據式(4)可以得到分離快慢橫波的兩種判別方法:第一種方法是在忽略噪聲干擾與快慢橫波不同的衰減影響條件下,快慢橫波應該是相似的,求取和的互相關便可分離快慢橫波[4-5],這種方法可以得到角度—時延譜;第二種方法是最小化斜對角線元素來分離快慢橫波。實際數據處理中,常常采用第一種方法來分離快慢橫波。

1.2 方位各向異性校正技術

轉換橫波分裂分析方位各向異性校正技術,主要實現以下幾個步驟。

1)根據徑向和橫向方位數據確定一個淺層分析時窗,在分析時窗內按照前面所述的快慢橫波分離技術求取快慢橫波。如果有多個方位角數據,則需要分別求取各個方位的快慢橫波,然后進行累加。

2)根據快慢橫波相似原則,求取快慢橫波互相關值,解釋角度—時延譜得到各向異性方向和快慢橫波時延。

3)將快慢橫波時延用于慢橫波上,作時延補償。

4)根據各向異性方位角,利用反旋轉技術將快橫波和時延補償后的慢橫波反旋轉回原來的徑向和橫向方向上,得到各向異性校正后的徑向和橫向分量。

5)在完成淺層時窗橫波分裂方位各向異性校正后的數據上,進行深層數據的分析和處理。

2 實際數據處理

將寬方位轉換波的徑向和橫向分量數據分成不同的方位角扇區,扇區劃分一般為36個或者12個。在每個方位角扇區內限制最大偏移距,以滿足各個方位數據均衡,并按照方位各向同性處理思路進行處理[6-8]。根據川西地區寬方位轉換波資料的特點設置淺層分析時窗時間為2.3～2.5 s。在該分析時窗內,首先確定各向異性方向、快慢橫波時差,再進行慢橫波時延補償,最后旋轉回原來的徑向和橫向分量。

圖1 淺層分析時窗內某個位置的角度—時延譜圖

圖2 淺層方位各向異性校正前、后徑向分量疊加剖面圖

圖1為淺層分析時窗內某個CDP位置的角度—時延圖,從圖1中可以看出,該位置的各向異性方向約90°(即東西方向),快慢橫波時差超過40 m s,表明淺層各向異性特征明顯。圖2是淺層方位各向異性校正前后的徑向分量疊加剖面對比,從圖2中可以看出,經過方位各向異性校正后的剖面的連續性得到了較大的提高,有效反射得到了加強,成像質量得到了明顯的改善。與此同時,淺層分析時窗以下的徑向分量數據的能量也發生了變化,說明各向異性具有“傳遞”作用,上層的各向異性會影響下層。

校正了淺層方位各向異性后,在此數據基礎上,設置中深層分析時窗為3.5～3.7 s。圖3為中深層分析時窗內與圖1相同CDP位置的角度—時延譜。對比圖3和圖1可知,中深層快慢橫波時差明顯小于淺層。這主要是因為校正淺層的方位各向異性后,中深層各向異性主要由地下裂縫引起,其強度要遠弱于上覆地層引起的方位各向異性。同理,利用橫波分裂方位各向異性校正方法對中深層數據進行處理,圖4為中深層校正處理后的徑向分量疊加剖面。對比圖4和圖2-b可知,中深層的成像質量得到了進一步的提高。

圖4 中深層方位各向異性校正后的徑向分量疊加剖面圖

3 結論

1)A lford旋轉是橫波分裂分析的較好方法,在構建正交的徑向和橫向分量對的情況下可以用于地面轉換波處理。

2)受近地表非均質性等因素的影響,轉換橫波方位各向異性特征主要體現在淺層。由此,為了提高轉換波的成像質量,至少需要消除淺層的方位各向異性。實際數據成像處理時,較好的實現方式是分為淺、中和深層3個時窗進行處理。

3)通過方位各向異性校正處理后,可以大大提高了徑向分量剖面的成像質量,這樣的數據更有利于縱、橫波聯合解釋和縱、橫波疊后聯合反演等研究。

[1]CRAM PIN S.A review of wave motion in anisotropic and cracked ealastic-media[J].Wave Motion,1981,3(4): 343-391.

[2]L IXiangyang,CRAMPIN S.Line-transfo rm techniques for p rocessing shear-wave anisotropy in four-component seismic data[J].Geophysics,1993,58(2):240-256.

[3]馬昭軍,唐建明,徐天吉.多波多分量地震勘探技術研究進展[J].勘探地球物理進展,2010,33(4):247-253.

[4]唐建侯,賀振華,黃德濟.角度譜的原理及計算方法[J].石油物探,2003,42(2):196-199.

[5]唐建侯,楊慧珠,賀振華,等.P-SV快慢波旋轉分離公式研究[J].天然氣工業,2004,24(12):44-47.

[6]黃東山,謝芳,李忠,等.多波多分量技術在廣安構造中的初步應用[J],天然氣工業,2006,26(8):44-46.

[7]馬昭軍,唐建明.疊前時間偏移技術在三維轉換波處理中的應用[J].石油物探,2007,46(2):174-180.

[8]馬昭軍,唐建明,楊振武,等.三維轉換波處理方法研究與應用[J].天然氣工業,2007,27(增刊1):213-215.

Alford Rotation-based an isotropic correction for the converted waves

Liu Hongai1,Ma Zhaojun2,Gan Qigang2

(1.School of Geophysics,Chengdu University of Technology,Chengdu,Sichuan 610059,China;2.Deyang B ranch of Exp loration and Development Research Institute,Sinopec Southw est B ranch,Deyang,Sichuan 618000,China)

NATUR.GAS IND.VOLUM E 31,ISSUE 3,pp.41243,3/25/2011.(ISSN 100020976;In Chinese)

When crossing a fracture,a shear wave may sp lit into a rapid shear wave p ropagating along the fracture and a slow shear wave p ropagating perpendicular to the fracture.In order to imp rove the imaging quality of radial component,azimuthal anisotropic correction for shear wave sp litting should be performed during the p rocessing of the converted waves.After creating a component pair w ith the radial and lateral components perpendicular to each other,the Alfo rd Rotation Method is used to split the converted waves into a rapid and a slow shear wave,and to perfo rm azimuthal anisotropic co rrection for shearwave sp litting.First,fracture o2 rientation and the rapid and slow shear waves are obtained in the shallow layer,and the time lag between the rapid and slow shear waves is also calculated.Then,the time lag is used to compensate the slow shear wave.Finally,the rapid shear wave and the time2 lag2compensated slow shear wave are rotated to their o riginal radial and lateral orientation respectively,obtaining new radial and lat2 eral components.With thismethod adop ted to p rocess the converted wave data in medium2to2deep o r deep layers in the western Si2 chuan Basin,the imaging quality of the radial component of the converted waves has been significantly imp roved.

converted wave,anisotropy,A lford Rotation,rapid wave,slow wave,Sichuan Basin,west,seismic data,imaging qual2 ity

國家科技重大專項專題“碎屑巖層系大中型油氣田富集規律與勘探關鍵技術”(編號:2008ZX05002-004)、中國石化重點科技攻關項目“HTI介質地震資料疊前偏移技術開發研究”(編號:P09076-T01)。

劉紅愛,女,1982生,碩士研究生;主要從事地震資料處理研究工作。地址:(610059)四川省成都市成都理工大學榕樹園8-2-2。電話:15183860725。E-mail:mazj_mazj@yahoo.com.cn

劉紅愛等.基于A lfo rd旋轉的轉換波各向異性校正技術.天然氣工業,2011,31(3):41-43.

10.3787/j.issn.1000-0976.2011.03.010

(修改回稿日期 2010-01-12 編輯 韓曉渝)

DO I:10.3787/j.issn.100020976.2011.03.010

L iu Hongai,bo rn in 1982,is studying for an M.Sc.degree,focusing on her research on seismic p rocessing.

Add:Mail Box 822 of Rongshuyuan,Chengdu University of Technology,Chengdu,Sichuan 610059,P.R.China

Mobile:+86215183860725 E2mail:mazj_mazj@yahoo.com.cn