地震屬性分析在深水濁積水道刻畫中的應用

樂 靖， 田 楠， 蔡文濤， 高云峰

(中海油研究總院 開發研究院，北京 100028)

?

地震屬性分析在深水濁積水道刻畫中的應用

樂 靖， 田 楠， 蔡文濤， 高云峰

(中海油研究總院 開發研究院，北京 100028)

開發初期, 由于缺乏鉆測井資料，難以刻畫深水濁積水道的發育期次，砂體內部疊置關系認識不清，導致油田開發效果不佳。為了表征深水濁積水道的發育期次和內部充填結構，對深水濁積復合水道期次劃分和邊界識別進行了研究，在傳統的地震屬性分析基礎上，形成了應用地震屬性分析技術識別深水濁積復合水道邊界和期次劃分的新思路。首先進行屬性優選，采用譜分解和多屬性融合的技術手段，結合地層切片分析，綜合測井資料來刻畫深水濁積復合水道邊界和期次，這對復合水道砂體疊置關系和連通性的認識有很大的幫助，指導了油氣田開發方案的設計以及優化。

深水濁積水道； 多屬性分析； 譜分解； 地層切片

0 引言

深水油田開發初期, 由于鉆井成本高，井數量有限，難以通過測井資料認清儲層特征，對井間的儲層非均質性和砂體展布特征的認識更多地依賴地震數據。然而深水濁積水道特有的形成條件和沉積環境使其內部伴生復雜的水道砂體，巖性和地層厚度在縱、橫向上變化大，給水道發育期次解釋及水道內部充填結構和演化特征、沉積相展布、高分辨率層序地層研究帶來極大困難。在開發初期，希望盡可能精細地表征砂體的空間展布、砂體之間的疊置關系，以便設計出最合理、高效的開發方案，提高經濟效益。針對油田開發生產的迫切需求，基于傳統的地震屬性分析技術，探索研究了應用地震多屬性分析技術識別深水濁積水道邊界及期次劃分的新思路，采用解釋性處理提高地震資料品質，進行屬性優選，采用譜分解和多屬性融合的技術手段，結合地層切片分析，綜合測井資料來刻畫深水濁積復合水道邊界和期次，該技術思路在實際應用中取得了較好的效果。

1 應用思路及關鍵技術

地震屬性是經過加工之后的地震數據，能夠從不同角度形象地反映地下的構造特征和儲層特征[]。因此，地震屬性分析是地震解釋的重要技術手段，傳統的地震屬性分析主要是優選某一種屬性來研究地下的儲層特征，由于地震資料本身的多解性，這種做法也具有較大的不確定性。不同的敏感屬性在反映儲層特征上各有優缺點，對這些敏感屬性綜合分析，可以有效綜合各種敏感屬性的優點反映儲層信息，降低地震資料的多解性，更有效地識別深水濁積水道的邊界和展布規律。將地震多屬性分析技術應用到深水濁積水道邊界識別及期次劃分的新思路如圖1所示，主要步驟如下：①分析地震資料品質，明確濁積水道的地震響應特征，對地震資料做解釋性處理，為下一步工作提供高信噪比的地震數據；②以地震資料同相軸為基本等時地層單元, 結合沉積地層認識，建立等時地層格架；③提取地震屬性，分析屬性對復合水道的響應，優選若干復合水道的敏感屬性；④譜分解地震沿層屬性，并采用RGB融合技術融合譜分解的地震屬性，初步識別復合水道的邊界和內部情況；⑤提取地層切片，分析復合水道內部次一級水道的沉積演化規律，認清各支水道之間的疊置關系；⑥綜合前面屬性分析及地層切片分析的認識，結合水道在地震剖面的響應特征和測井、沉積相的認識，精細解釋復合水道內部的各支次一級水道。

圖1 地震多屬性分析在深水濁積水道期次劃分中的應用思路Fig.1 Application of seismic multi attribute analysis in the sub division of deep water turbidity channel

這一技術思路中的關鍵技術為：

1)地震資料解釋性處理技術。地震資料中的噪音降低了資料的信噪比，掩蓋或降低了儲層在地震屬性中的響應特征，影響地震屬性分析的效果。因此在提取地震屬性之前，先對地震資料做解釋性處理，進一步壓制地震資料中的噪音，提高資料的信噪比。采用的技術方法主要有中值濾波及擴散濾波，①中值濾波可以有效去除地震資料的隨機噪音；②擴散濾波能有效增強異常體的顯示(增強斷層、沉積體、特殊巖性體的顯示)。采用地震資料解釋性處理技術，能有效去除隨機噪音，增強斷層和沉積體的顯示，提高地震資料的信噪比，為地震屬性分析提供一個良好的數據基礎。

2)等時地層格架建立技術。等時地層格架能夠更好地避免常規地震切片中存在的穿時問題[2]，在地震資料同相軸為基本等時地層單元的假設下, 結合沉積地層認識，通過井震標定賦予地震同相軸對應的地質層序界面的意義，根據地震資料同相軸產狀變化進行層位解釋。為提高層位解釋精度，采用層位自動追蹤技術，為下一步分析提供精細的等時地層格架。

3)譜分解技術。深水濁積水道的厚度變化很快，在水道邊界砂體變薄，而地震數據的分辨率有限，難以精確識別水道邊界的薄砂體并刻畫水道邊界。不同厚度的砂體其地震振幅在橫向上是有差異的，譜分解技術基于振幅的差異，尋找對薄砂體有最好反映的頻率切片，在平面上結合沉積相和測井信息分析水道砂體的展布和邊界[3]。因此，采用譜分解技術，可以有效地識別深水濁積水道，并刻畫深水濁積水道的邊界。

4)RGB(紅、綠、藍三原色融合顯示)融合技術。地震資料本身存在的多解性導致進行單一屬性分析時不確定性較大，為了降低屬性分析的不確定性，不能僅依靠一種屬性進行分析。由于不同的屬性能夠分別從不同的角度反映地下的儲層特征，把這些對儲層響應較敏感的屬性進行綜合顯示分析，能夠綜合各種屬性的優勢，突出儲層的響應特征，降低單一屬性分析的不確定性，RGB融合技術就是一種能將多種屬性融合起來進行綜合分析，能夠有效降低單一屬性顯示的不確定性的技術方法[4]。將多種屬性以RGB模式融合顯示之后，相比單一屬性顯示，能更好地反映砂體內部的細節變化，結合沉積模式能更好地分析砂體的疊置關系和展布范圍[5]。對于深水濁積水道，僅僅依靠單一屬性的分辨率很難識別，而基于多種地震屬性，通過RGB顏色融合技術，可以大大增強復合水道的視覺異常，減少單一地震屬性的多解性，可以大幅度提高深水濁積水道的識別精度。

5)地層切片技術。在屬性分析過程當中，僅局限于地震屬性本身，不結合沉積學的觀點和知識進行屬性分析是難以真正認識儲層的變化的，由此提出了地層切片技術[6-7]。地層切片技術要求以沉積學的觀點去認識不同時期的屬性變化，有助于更有效地認清沉積體的沉積演化規律，以及刻畫沉積體的期次和邊界。地層切片的參考界面是不隨頻率而變化的等時地震同相軸，因此地層切片最接近于等時，適用于非席狀、非水平的不規則地層，充分考慮了地層厚度的橫向變化，克服了穿時帶來的不足，可以有效地識別深水濁積水道，刻畫深水濁積水道內部的非均質性。

2 實際應用效果

研究區位于西非的剛果盆地，盆地內發育了世界第二大規模的深水扇-剛果扇。研究區位于中新世剛果扇北部，目的層為上中新統。中新統目的層段主要發育深水濁積水道、遠端朵葉體和半深海泥巖3種沉積相，且以侵蝕沖刷面為層序界面，大面積發育的復合水道為層序底部，半深海泥巖密集段為層序頂部[9]。儲層段是大面積發育的復合水道，水道數量繁多，期次復雜，相互切割，內部關系難以理清；水道內儲層非均質性強，橫向變化大，同一水道內地震反射特征相似，但井上的含砂率相差較大。開發初期井少，多期、多支水道間砂體疊置關系和連通性需要有效刻畫，以便于設計合理、高效的開發方案。針對生產需求，在研究區探索研究應用地震多屬性分析技術識別深水濁積水道邊界及期次劃分的新思路。

2.1 地震資料解釋性處理

研究區的地震資料信噪比較低，隨機噪音模糊了斷層、復合水道的邊界，影響屬性分析的精度。需采用地震資料解釋性處理技術對地震資料進行處理，從圖2可以看出，處理之后，去除了隨機噪音，提高了資料的信噪比，增強了斷層、復合水道邊界，降低了屬性分析的多解性，提高了屬性識別復合水道的精度。

圖2 解釋性處理前后相干切片對比圖Fig.2 Comparison of coherent slices before and after interpretation processing(a)處理前相干切片；(b)處理后相干切片

2.2 對深水濁積水道邊界進行識別

與周圍的泥巖的弱反射特征相比，深水濁積水道主要表現為強振幅、中低頻的雜亂發射，形成在泥巖中高頻弱反射背景上的強的丘狀或透鏡狀反射結構。因此提取振幅、頻率、相干等屬性能有效地反映水道的變化，識別水道砂體與泥巖的邊界。圖3分別展示了目的層段沿層相干切片，以及20 Hz、30 Hz、40 Hz、50 Hz的目的層段沿層相干切片。相干切片能在一定程度上刻畫較大級次的水道復合體，但是難以刻畫水道主體，而單頻相干切片能更好地刻畫水道主體。由圖3可以看出，當頻率為20 Hz、30 Hz時，水道的主體及邊界非常清楚，斷層比較清晰，頻率增大到40 Hz、50 Hz時，水道邊界逐漸變模糊。因此頻率為30 Hz的沿層相干切片能較好的識別各期水道的主體及邊界。不同頻率的相干切片

圖3 20 Hz、30 Hz、40 Hz、50 Hz的沿層相干屬性Fig.3 Coherence properties along the horizon of 20 Hz、30 Hz、40 Hz、50 Hz(a) 20 Hz；(b) 30 Hz；(c)40 Hz；(d) 50 Hz

在水道識別上各有優缺點，如果僅通過其中某一頻率的相干切片識別水道邊界，不確定性較大且對水道內部的變化表征得不夠精細，因此優選不同頻率體融合顯示[8]。如圖4所示，將譜分解得到的10 Hz、40 Hz、70 Hz不同頻段的數據分別用為R(紅色)、G(綠色)、B(藍色)三原色進行融合顯示，清楚地刻畫了水道的邊界, 并對水道內部的非均質性有一定程度地反映, 較常規屬性更精細地刻畫了濁積水道的邊界以及內部的細節變化。

2.3 水道發育期次的劃分。

水道復合體通常由多期次一級水道構成，采用譜分解、RGB融合技術之后，可以識別水道復合體的邊界，并認識到次一級水道之間具有復雜的相互疊置、切割關系，為了進一步精細刻畫水道期次，采用了地層切片技術[9]。

將復合水道的頂底面作為參考層，按照8等分內插得到8張地層切片，并沿切片提取均方根屬性。從圖5可以看出，按照從老至新排列后觀察到的復合水道內部期次的展布、演化特征為：

圖5 水道發育情況在地層切片上的反映Fig.5 Channels development reflected in the strata slices(a)地層切片1；(b)地層切片2；(c)地層切片3；(d)地層切片4；(e)地層切5；(f)地層切片6；(g)地層切片7；(h)地層切片8

1)在地層切片1上顯示的紅黃色及黃綠色強振幅反映了兩期水道的存在，主要發育了兩支北西—南東向水道B和C的雛形，水道內部顯示為藍色背景下的紅黃色及黃綠色斑點分布，即富泥質背景下的砂巖零星分布。

2)在地層切片2上，水道B、C的輪廓開始變得清晰，水道邊界比較清楚，內部振幅增強，紅黃色的強振幅增多，兩支水道在南東方向相互疊置和切割，在北西方向兩支水道是分開的。

3)在地層切片3上，水道B、C在北西向的振幅變弱，南東向的振幅較強，紅黃色斑點較多，兩支水道在南東方向依然相互疊置和切割，兩支水道在北西方向是分離的，邊界清晰。

4)在地層切片4上，水道B、C在北西向的振幅稍微變強，紅黃色斑點增多，而且在兩支水道之間出現了第三支水道A的雛形。水道A在北西向的振幅較弱，它是在水道B、C之后形成的，在南東向它與水道B、C的頂部相互疊置、切割。在北西向水道A、B、C的邊界比較清楚，相互之間是分離的。南東向的振幅最強，水道A疊置在水道B的上方，從平面上看兩支水道在南東向幾乎重疊。

5)在地層切片5上，水道C基本消失了，水道B在北西向的振幅變弱，振幅顏色偏向于藍色，水道B也在逐漸消失；水道A變得更清晰，振幅也在逐漸增強，振幅顏色偏向于紅黃色。在南東向水道A的振幅最強，以紅黃色斑點為主，水道A疊置在水道B的上方，從平面上看兩支水道在南東向幾乎重疊。

6)在地層切片6上，水道B基本消失了，水道A的邊界比較清晰，振幅比較強，振幅顏色以紅黃色為主。

7)在地層切片7上，僅有水道A存在，水道邊界仍然清晰，水道內部是弱振幅背景下的強振幅斑點分布，往南北兩個方向振幅變弱。

8)在地層切片8上，水道A在南北兩個方向的振幅變得更弱，以藍色的弱振幅為主，僅有水道中部的振幅較強，以紅黃色斑點為主，這表明水道A也在逐漸消失。

從這8個地層切片可以看出，各次一級水道的發育演化規律為：水道B和C幾乎同時出現，接著水道A發育；水道復合體內部由多期次一級水道構成，后期水道對已有水道存在一定程度的切割和疊置；將切片按時間順序排列，能識別出各個次一級水道的發育期次，明確各個水道之間的切割、疊置關系。

2.4 井震結合解釋水道期次

根據前面從地層切片分析得到的水道發育演化規律，結合從巖心分析、測井資料分析得到的沉積模式的認識，在地震資料上開展水道期次解釋，解釋出水道的期次及展布特征。如圖6所示，從圖6可以看出，水道的地震響應是雜亂無章的中強反射，外形呈透鏡狀，水道內部反射頻率低，連續性差。該復合水道是侵蝕-加積混合型水道，垂向疊加，側向疊置。結合地層切片得到的水道發育演化規律，將復合水道解釋為3個次一級水道，利用多屬性的分析結果相互驗證，增強了結果的可靠性。

圖6 地震剖面上的水道解釋結果Fig.6 Interpretation results of the channel on seismic profile

3 結論

通過對傳統的屬性分析技術進行優化組合，相位、相干、振幅等屬性有效地識別和刻畫了深水濁積復合水道，創新的地震多屬性分析思路在實際應用中取得了良好的效果。研究的主要認識如下：

1)深水濁積水道由于特有的形成條件和沉積環境使其內部伴生復雜的水道砂體，巖性和地層厚度在縱、橫向上變化大，對應的地震反射通常是雜亂無章的中強反射，僅依靠單一的地震屬性難以解釋水道發育期次及內部結構，采用地震多屬性分析的新思路，從提高資料的信噪比入手，結合各種多屬性分析的技術手段能夠有效地刻畫水道的發育期次及內部結構。

2)相干屬性經過分頻、RGB融合處理后，對復合水道的邊界識別更清楚，對復合水道的內部結構也有了更清晰地反映。

3)地層切片是在不隨頻率變化的等時地震同相軸之間線性內插形成的，所提供的地質信息更接近于等時，克服了時間切片和沿層切片的不等時帶來的不足，使地質體和沉積體的識別更加精確、合理，清晰地刻畫了水道的發育期次和演化規律。

4)深水濁積復合水道內部通常由多期次一級水道構成，后期水道對已有水道存在一定程度的切割和疊置。將地層切片按時間順序由老到新排列，能識別出水道的發育期次和演化規律，再結合測井相，可以精細解釋次一級水道，明確次一級水道間的疊置關系及連通性，為開發方案的設計提供精細的儲層描述成果。

[1] CHOPRA S, MARFURT KJ. Seismic Attributes for Prospect Identification and Reservoir Characterization[C]．Tulsa: Society of Exploration Geophysicists, 2007:1-7.

[2] 朱筱敏，董艷蕾，胡廷惠，等．精細層序地層格架與地震沉積學研究——以泌陽凹陷核桃園組為例[J]．石油與天然氣地質，2011，32(4)：615-624.

ZHU X M ，DONG Y L，HU T H，et a1．Seismic Sedimentology Study of Fine Sequence Stratigraphic Framework：a Case Study of the Hetaoyuan Formation in the Biyang Sag[J]．Oil&Gas Geology，2011，32(4)：615- 624.(In Chinese)

[3] 劉喜武，寧俊瑞，劉培體，等．地震時頻分析與分頻解釋及頻譜分解技術在地震沉積學與儲層成像中的應用[J]．地球物理學進展，2009，24(5)：1679-1688．

LIU X W，NING J R，LIU P T，et al．Seismic Time Frequency Analysis for Frequency Decomposition With Applications to Seismic Sedimentology and Reservior Imaging[J]．Progress in Geophysics，2009，24(5)：1679-1688. (In Chinese)

[4] 曹鑒華．RGB混頻顯示技術及其在河道識別中的應用[J]．勘探地球物理進展，2010，33(5)：355-358．

CAO J H．RGB Color Blending and Its Application In Channel Recognition[J]．Progress in Exploration Geophysics，2010，33(5)：355-358．(In Chinese)

[5] 姜秀娣，翁斌，劉亞茹，等．分頻混色技術在高精度地震解釋中的應用[J]．地球物理學進展，2013，28(2)：0882-0888.

JIANG X D，WENG B，LIU Y R，et al．Application Of Spectral Decomposition RGB Plotting Technique For Spectral Components in High Accuracy Seismic Interpretation[J]．Progress in Geophys．2013，28(2)：0882-0888. (In Chinese)

[6] ZENG H L， BACKUS M M，BARROW KT ,et al．Office Stratal Slicing, Part I: Realistic 3-D Seismic Model[J]．Geophysics， 1998，63( 2) : 502-513.

[7] ZENG H L，HENRY S C AND RIOLA J P．Stratal Slicing, Part II: Real 3-D Seismic Data[J]．Geophysics，1998，63( 2) : 514-522.

[8] 郎曉玲，彭仕宓，康洪全，等．利用多屬性體分類技術預測扇三角洲砂體[J]．西南石油大學學報：自然科學版，2010，32(1)：57-62.

LANG X L， PENG S M, KANG H Q， et al．Predict Fan Delta And Body by Using Multi-attributes Volume Classification Technology[J]．Journal of Southwest Petroleum University ： Science & Technology Edition, 2010，32( 1 ) ：57-62. (In Chinese)

[9] 謝清惠，鄧宏文．地層切片技術在中新統剛果扇中的應用[J]．大慶石油地質與開發，2013，32(6)：135-140.

XIE Q H,DENG H W. Application Of Strata Slicing Technique In Miocene Congo Fan [J]．Petroleum Geology and Oilfield Development in Daqing, 2013，32(6)：135-140. (In Chinese)

Application of seismic attribute analysis in the sub division of deep water turbidity channel

LE Jing, TIAN Nan, CAI Wen-tao, GAO Yun-feng

(CNOOC Research Institute Development Research Institute, Beijing 100028,China)

In the early stage of development, it is difficult to describe the development period of the deep water turbidity channel and the internal structure of the sand body because of the lack of logging data, which leads to the poor development of the oilfield. In order to characterize the developmental period and internal structure of the deep water turbidity channel, it is studied that how to identify the boundary and internal structure of the deep water turbidity composite channel. Based on the traditional seismic attributes analysis, a new method for identifying the boundary and stage division of the deep water turbidity composite channel is formed by using seismic attribute analysis. Attribute optimization, adopting spectral decomposition and multi attribute fusion, strata slice analysis, and combining logging data to characterize the boundary and internal structure of the deep water turbidity composite channel. This is help to know superimposed relationship and connectivity of deep water turbidity composite channel, which is to guide the design and optimization of oil and gas field development plan.

deep water turbidity channel； multi attribute analysis； spectral decomposition； strata slice

2015-08-30 改回日期：2015-10-08

樂靖(1984-)，男，碩士，從事地震綜合解釋，E-mail：lejing@cnooc.com.cn。

1001-1749(2016)05-0686-06

P 631.4

A

10.3969/j.issn.1001-1749.2016.05.19