基于分傾角掃描的相干屬性在斷層識別中的應用

(1.中國石油 大慶油田有限責任公司 第六采油廠 第四油礦， 大慶 163114；2.中國石油 東方地球物理公司研究院， 涿州 072751；3.中國石油 西南油氣田分公司 西北氣礦， 綿陽 621000)

?

基于分傾角掃描的相干屬性在斷層識別中的應用

秦 瑞1， 白 軍2， 鄭 超3， 徐 爽2， 鄒 雯2， 馬 聰2， 曹 淼2

(1.中國石油 大慶油田有限責任公司 第六采油廠 第四油礦， 大慶 163114；2.中國石油 東方地球物理公司研究院， 涿州 072751；3.中國石油 西南油氣田分公司 西北氣礦， 綿陽 621000)

現代體屬性的發展使斷層識別能力得到很大提高，尤其以相干屬性為主的一系列斷層識別方法得到了廣泛應用，并取得較好的效果。然而直接從地震數據提取相干屬性對于資料的信噪比要求較高，且由于未考慮地層傾角的影響往往難以準確全面地識別斷層。為解決這一問題，對研究區的地震資料進行了分角度傾角掃描和構造導向濾波，并在此基礎上提取相干屬性，從而提高了對斷層識別的精度。

分傾角掃描； 構造導向濾波； 相干屬性； 斷層識別

0 引言

隨著屬性技術的不斷發展，對于斷層識別精度的要求也越來越高。通常在構造相對簡單的區域，地震資料品質一般較好，斷層相對容易識別。對于構造復雜的區域，斷層的規模不同、地層傾角變化劇烈，地震資料的信噪比低，斷層的識別難度較大。趙明章等[1]認為基于原始地震數據體直接進行相干體分析的方式簡單快速但受地震資料品質的限制，對于斷層的解釋往往不太準確，尤其是小斷層的響應很可能被噪聲覆蓋掉；同時Chopra等[2]認為屬性的發展正朝著多屬性綜合應用的方向發展。為此，本文聯合多種屬性，有針對性的對斷層進行預測，預測效果相對于常規方法有很大改善。具體來說是通過對地震數據進行分角度傾角掃描，再將不同的傾角體作為導向依據對地震數據進行構造導向濾波，繼而再進行相干分析對斷層進行識別。

1 方法原理

1.1 傾角掃描

傾角、方位角是確定任一平面在空間中展布特征的屬性。傾角(Dip，θ)是偏離水平面的角度，在地質上又可分解為在x和y方向的視傾角；方位角(Azi，Φ)指反射軸與水平面的交線偏離正北方向的角度[2](圖1)。

圖1 傾角、方位角示意圖[1]Fig.1 The schematic diagram of dip angle and azimuth angle[1]

傾角和方位角的屬性體現在兩個方面：

1)聯合傾角、方位角屬性可以指示地層的產狀，解釋確定的局部反射層面，直接用于估算反射層的不連續性。

2)傾角、方位角信息可以為數據濾波提供導向依據，突出特定方向上反射特征的不連續性。通過估算地震反射界面的傾角和方位角為構造導向濾波提供導向依據，濾波后可在一定程度上提高地震資料的信噪比，降低白噪聲和相干噪聲，同時保留斷層、裂縫或其他有價值構造的反射信息, 保證了其不連續性, 提高了對斷層的分辨率[1]。

目前，常用的傾角估算法包括基于復數道分析的相位對齊法、離散掃描最大相關法與梯度構造張量法[3]。離散傾角掃描最大相關法就是選取了分析點周圍多個時窗進行掃描，求取相似度最大的窗口作為分析點處的傾角、方位角估算窗口，從而大大提高了分析點的傾角與方位角精度，對于識別小斷層、小構造具有重要的意義[4]。

1.2 相干屬性

相干是表征地震道之間相似程度的量。主要分析以目標點為中心的時窗內相鄰地震道波形的相似性，以判斷地層的連續性，即利用波形之間的相似性來反映地層的不連續性特征。根據地震勘探理論，剖面上連續性好的同相軸往往對應于地下連續性較好的地質體(如水平地層或傾斜地層等)；連續性差的同相軸對應于連續性較差的地質體(如斷層等)[5～8]。計算地震相干體主要目的是突出地震數據中不連續的部分，強調的是地震道之間的不相關性。通過地震數據體來計算局部地震波形的相似性，可以得到三維地震相關性的估計值；地震道之間的波形特征發生變化時，會導致局部的道與道之間相關性的突變，因此可以利用相干來預測斷層、地層巖性突變和特殊地質體的發育區[3]。

地震相干體的算法很多，至今主要經歷了三次算法的變革，這里主要采用第三代算法——特征值分析算法。用特征值法進行相干估算是一種計算輸入地震道波形相似性的方法，特征值法也叫做主分量分析法，該算法將數學中的矩陣特征結構引入到相干分析中,利用特征結構來計算相干，通過計算協方差矩陣的特征值來得到相干屬性[2]。

Gersztenkorn等[5]把2道互相關的相干體算法推廣到任意J道，這種方法需要初始估算的傾角、方位角，而特征值結構的相干體算法本身就可以實現按一定間隔進行傾角或方位角掃描。首先從地震數據體中以一定的時窗范圍(t±k)提取樣點矢量，在時間t的視傾角p和q決定的平面上，通過對數據樣點進行插值，生成新數據矩陣的行和列，并對其進行自相關和互相關運算，生成協方差矩陣C。此時，協方差矩陣的元素Cij可以表示為：

qyi)-u(t,p,q)]*[uj(t+

kΔt-pxi-qyi)-u(t,p,q)]

(1)

其中樣本均值u(t,p,q)為：

(2)

通過以上公式定義式(3)為第三代相干體的C3相干值：

(3)

式中：分母是矩陣的跡，代表了協方差矩陣的能量；分子λ1是最大特征值，代表了占優的能量。根據協方差矩陣，計算其本征值及其對應的本征向量[9]。

1.3 基于分傾角掃描的相干屬性

在實際應用中往往會遇到這樣的情況，同一工區內分布有多個構造區帶，不同的構造區帶地層傾角往往變化較大，通常會直接由地震數據或使用單傾角進行掃描后提取相干屬性，但這兩種方法都不能準確全面的識別斷層信息；由地震資料直接提取相干屬性往往受地震資料的品質影響而無法清晰識別斷層；通過單傾角掃描處理之后再提取相干屬性雖然效果相對前者較好，但由于傾角掃描范圍和步長單一會使得一些有用的斷層信息被忽視掉。為此，作者提出了分傾角掃描的相干技術來解決這一問題。基于分傾角掃描的相干技術實現步驟如下：①針對同一工區使用不同的傾角范圍和掃描步長進行傾角掃描得到多個不同的傾角體；②在此基礎上進行適當構造導向濾波；③提取相應的相干屬性，綜合多個相干屬性來識別不同地層傾角區域的斷層，最終能夠對不同級別斷層有一個準確全面的認識。具體來說就是首先對于地層傾角有一個初步估算，然后劃分出傾角掃描時的傾角范圍和步長，小的傾角范圍使用小步長可以突出地層平緩區小尺度的斷層，大的傾角范圍使用大的步長可以突出高陡區大尺度的斷層；為了使顯示效果更好可適當使用構造導向濾波技術，最后利用不同傾角掃描的數據體分別提取相干屬性，針對不同區域使用不同的相干屬性來輔助斷層解釋，以達到全面準確認識斷層的目的。[10-13]

2 應用分析

研究的區域主要分布兩個構造區帶①高陡區；②平緩區(圖2)。高陡區地層傾角較大，約有15°，主要分布有高角度的逆斷層，斷距相對較大；平緩區地層近似水平，總體不超過1°，主要是由灰巖溶蝕形成的陷落柱(類似于小斷層)以及小斷距的逆斷層。兩個區帶地層傾角不同，斷層尺度也不同，直接提取相干屬性或是單一傾角掃描再提相干屬性都不能全面識別斷層。因此，根據研究區地地層傾角和構造的分布特征，首先分別對地震數據進行小傾角范圍(0°～1°)和大傾角范圍(5°～15°)的傾角掃描得到目的層段兩個不同的傾角體,在此基礎上再進行相干分析。通過對比多種方式得到的相干屬性對斷層的識別效果，我們發現基于分傾角掃描的相干屬性對于斷層的識別更清晰更全面(圖3)。經過傾角掃描和構造導向濾波后提取的相干屬性(圖3(b)和3(c))相對與直接從成果數據中提取的相干屬性(圖3(a))相比，前者對斷層的分辨能力明顯增強，平面上噪聲干擾明顯降低，剖面上斷點相對清晰。

圖2 區域構造剖面Fig.2 The structural section of research area

圖3 基于原始地震資料和分傾角掃描 的相干屬性對比圖Fig.3 The coherence attribute contrast figure of the original and separated dip scanning seismic data(a)原始地震資料;(b)0°～1°傾角掃描;(c)5°～15°傾角掃描;(d)斷裂特征

同時不同的傾角掃描范圍對不同尺度斷層的識別能力也不同：傾角掃描范圍為0°～1°的相干屬性(圖3(b))明顯對于平緩區小斷距的斷層(陷落柱)的識別能力較好，但同時對高陡區的斷層識別能力變差，噪聲干擾明顯增強；相反，傾角掃描范圍為5°～15°的相干屬性(圖3(c))對于高陡區斷距相對較大的逆斷層識別效果較好，反而是平緩區的小斷距斷層(陷落柱)基本無法識別。

基于以上認識，結合大小兩種不同的傾角掃描得到的相干屬性來全面的識別全區斷層，在用原始地震數據提取的相干屬性和剖面斷層識別結果進行檢驗，最終準確全面的認識了全區斷層，并繪制了目的層斷裂特征圖(圖3(d))。

3 結論

通過對比以上應用效果，發現將分傾角掃描技術與三維相干屬性結合，能夠更全面、更有效地識別斷裂。相比于常規地震相干解釋技術，該方法以地層傾角作為導向因素，對于識別特定走向的斷層針對性更強，效果更明顯。同時分傾角掃描的方法又可以根據地層傾角的不同全面地識別出不同區域不同尺度的斷層。總的來說，該方法適用性強、效果明顯，值得在其他工區推廣。

需要指出的是，本方法在應用過程中有以下幾方面需要注意：

1)在傾角估算的過程中，最大傾角和掃描步長的選取對于識別斷層的效果有很大影響，傾角范圍過小，雖然對于小斷層的刻畫更為精細但地震同向軸的連續性變差，可能會有一些假象出現；傾角范圍過大，干擾可能減少但對于細節刻畫的精細程度會降低。

2)作者進行傾角掃描后，由于本區構造走向較為單一，所以選擇作為導向的傾角體是垂直于構造走向的。對于構造更為復雜的區域，除分傾角掃描外，在濾波時應考慮方位信息，從多個方位進行濾波再做相干，這樣對于不同走向的斷層識別效果可能更好。

3)分傾角掃描實質上是對原始資料做了某種意義上的濾波用以凸顯其某些不連續的地震反射特征，濾波難免會損失一些信息或造成一些假象，這些都應當在后期做好質控工作，減少這些因素的影響。

[1] 趙明章，范雪輝，劉春芳．利用構造導向濾波技術識別復雜斷塊圈閉[J]．石油地球物理勘探，2011，46(增刊1)：128-133.

ZHAO M Z,FAN X H,LIU C F.Fault-blocked traps identification based on structure-oriented filtering[J].Oil Geophysical Prospecting,2011,46(suppl 1):128-133.(In Chinese)

[2] CHOPRA S,MARFURT K J.Seismic attribute expression of differential compaction[J].The Leading Edge,2012,31(12)：1418-1421.

[3] 王霞，汪關妹,劉東琴.地震體屬性分析技術及應用[J].石油地球物理勘探，2012,47(增刊1)：382-389.

WANG X,WANG G M,LIU D Q.Seismic volume attributes analysis technology and application[J]. Oil Geophysical Prospecting,2012,47(suppl 1):382-389. (In Chinese)

[4] 蔡涵鵬，賀振華，李亞林．基于多窗口相干性的傾角導向主分量濾波[J]．石油地球物理勘探，2014，49(3)：486-494.

CAI H P,HE Z H,LI Y L.Structure-oriented filtering based on multi-window dip scanning[J]. Oil Geophysical Prospecting,2014,49(3):486-494.(In Chinese)

[5] 王大偉,劉震,陳小宏.地震相干技術的進展及其在油氣勘探中的應用[J].地質科技情報，2005，24(2)：56-58.

WANG D W,LIU Z,CHEN X H.The progress of the seismic coherence technology and its application in oil and gas exploration[J].Geological Science and Technology Information,2005,24(2):56-58. (In Chinese)

[6] 徐中清,周海明.復雜斷塊油田精細勘探開發技術[M].北京:石油工業出版社, 2002.

XU Z Q,ZHOU H M.Technology of complex reservoirs exploration and development[M].Beijing:Petroleum industry press,2002.(In Chinese)

[7] 楊威，賀振華，陳學華.結構方位濾波在體曲率屬性中的應用[J].石油物探，2011,50(1)：27-32.

YANG W,HE Z H,CHEN X H.The application of structure bearing filter in curvature volume attributes.Geophysical Prospecting for Petroleum,2011,50 (1):27-32.(In Chinese)

[8] 王彥輝,王歡,李綱.小斷層的識別方法[J].大慶石油地質與開發，2000,10：31-32.

WANG Y H,WANG H,LI G.Identification methods of small fault[J].Petroleum Geology and Oilfield Development in Daqing,2000,10:31-32.(In Chinese)

[9] MARFURT KJ,KIRLIN RL.3D seismic attributes using a semblance-based coherency algorithm[J].Geophysics,1998,63(3)：1024-1035.

[10]張馳，張軍華，黨玉芳.基于多窗口傾角掃描的曲率屬性識別小斷裂方法研究[R].中國地球科學聯合學術年會,2014：1076.

ZHANG C,Z J H,DANG Y F.Fault identify- cation of curvature attributes based on multi-window dip scanning[R].Chinese earth science joint academic annual meeting,2014:1076.(In Chinese)

[11]崔世凌，張軍華，王偉.地震紋理屬性在JJD工區斷層識別中的應用[J].物探化探計算技術，2010,32(3)：304-309

CUI S L, ZHANG J H,WANG W.Fault identify- cation based on seismic texture attributes in the work area JJD[J].Computing Techniques for Geophysical and Geochemical Exploration,2010,32(3):304-309.(In Chinese)

[12]雷英成.傾角掃描地震紋理對比度屬性研究及應用[J].斷塊油氣田,2011，18(4)：465-468.

LEI Y C.Research and application of seismic texture contrast with dip sweeping[J].Fault-Block Oil and Gas Field,2011,18(4):465-468.(In Chinese)

[13]GIJS C FEHMERS,CHRISTIAN F W HOCKER.Fsat structural interpretation with structure-oriented filtering[J].Geophysics,2003,68(4)：1286-1293.

The application of coherence attribute based on separated dip scanning in fault identification

QIN Rui1， BAI Jun2，ZHENG Chao3，XU Shuang2，ZOU Wen2，MA Cong2，CAO Miao2

(1. The Fourth Oilfield of Sixth Oil Production Plant,Daqing Oilfield,CNPC,Daqing 163114,China； 2.Geophysical Reaserch Institute,BGP,CNPC,Zhuozhou 072751,China; 3. Northwest Gas Field,Southwest Oil and Gasfield Company,CNPC,Mianyang 621000,China)

The development of modern volume attributes has improved the fault identification ability, especially the coherent attributes have been widely used and acquired great effect. However, extracting the coherence attributes from seismic data requires higher S/N ratio, and it is difficult to identify the fault comprehensively without considering the influence of stratigraphic dip. To solve this problem, analysis of separated dip scanning and structure-oriented filtering has been made and then the coherence attributes have been extracted to improve the accuracy of fault identification in this study.

separated dip scanning; structural orient fliter; coherence attribute; fault identification

2015-09-11 改回日期：2015-11-20

秦瑞(1990-)，男，碩士，主要從事地震解釋與儲層預測方面的工作，E-mail：849887410@qq.com。

1001-1749(2016)05-0656-04

P 631.4

A

10.3969/j.issn.1001-1749.2016.05.13