基于地震屬性分析技術的棗35斷塊火山巖分布預測

周天琪

(中國地質大學(武漢)資源學院，湖北 武漢 430073)

季嶺

(中國石油勘探開發研究院研究生部，北京 100083)

謝叢姣

(中國地質大學(武漢)資源學院，湖北 武漢 430073)

杜樂樂

(中國地質大學(北京)研究生部，北京 100083)

棗35斷塊位于黃驊坳陷風化店構造北端李天木斷層下降盤的自來屯構造，整體構造為東南高、西北低，構造形態較簡單[1-2]，為沙三段火山巖裂縫型-孔隙型雙重介質稠油油藏。該區火山巖厚度不均，巖性復雜，大部分井缺少巖性描述資料，給勘探開發帶來了很多的困擾。為此，筆者在分析研究前人研究成果的基礎上[3-5]，采用井震聯合分析技術，歸納總結火山巖體地震剖面的反射特征、火山巖體不同部位的振幅譜特征以及火山巖體頂部振幅的變化規律。通過將火山巖地震相特征、振幅類信息、分頻信息3者緊密結合，可以很好識別火山口的發育特征與位置，并繪制出火山巖的橫向展布特征。

1 火山巖地震相反射特征的識別與分析

由于火山巖密度較大，與圍巖間形成較強的波阻抗界面，因此火山巖的地震反射具有較強的獨特性和可識別特征，故分析火山巖地震相的反射特征能為火山巖識別工作建立良好的基礎。在三維地震剖面上，通過對研究區內26口井進行合成地震記錄的制作，對沙三段火山巖頂底進行層位標定，發現研究區內火山巖主要地震反射特征為3類，分別是丘狀地震反射、柱狀地震反射和層狀地震反射。

1.1 丘狀地震反射特征

圖1 丘狀地震反射原始地震剖面

丘狀地震反射是巖漿沿斷層運移的表現特征。研究區內丘形反射(見圖1)主要表現為頂部界面為連續的強軸反射軸；內部多為連續性較差的弱軸反射軸，反射結構呈波狀、珠串狀斷續分布，互為交叉；在部分丘形反射兩翼可見沉積地層出現超覆沉積現象。由于古近紀時期的風化剝蝕，大部分火山體已被風化剝蝕，因此部分剖面可見丘形反射頂部出現磨平、塌陷的現象。

1.2 柱狀地震反射特征

柱狀地震反射是火山口在地震剖面上的表現。其整體形態呈煙囪狀，2條陡直斷裂直通基底；頂面反射軸較為連續，振幅較強；內部同相軸多不連續，呈蚯蚓狀、波狀呈斷續分布，反射較為雜亂。原因為火山頸內部巖性結構變化快，分布范圍小，其中含玄武巖等多種噴發巖，故導致內部反射軸成層性差[6]。可在LINE2141、LINE2154剖面上觀察到柱狀地震反射(見圖2)，其中Z78井(LINE2157)、Z66井鉆遇火山角礫巖。根據巖性資料以及前人所描繪的火山巖巖相模式[7]，推斷柱狀地震反射是研究區內火山口位置。

(a) LINE2141 (b) LINE2154圖2 柱狀地震反射原始地震剖面

圖3 層狀地震反射地震剖面

圖4 棗35斷塊火山巖底(0～20ms)RMS地震屬性圖

1.3 層狀地震反射特征

層狀地震反射為層狀溢流相火山巖(致密玄武巖，氣孔玄武巖)在地震剖面上的反映，頂層同相軸較為連續、平直，多具有強振幅等地震反射；對下伏地層具有屏蔽作用，使得上下地層反射變弱[8]。由于研究區內沙三段火山巖巖性結構多變且厚度不大，故該反射特征不具典型性，無法較好劃分溢流相的分布范圍(見圖3)。

2 火山巖地震屬性分析

利用地震屬性技術可以快速有效的識別地下巖體，火山巖獨特的反射特征為利用地震屬性快速識別提供了基礎。

最大波峰振幅(Maximum Peak Amplitude)、總絕對振幅值(Total Absolute Amplitude)、總振幅(Total Amplitude)等振幅類地震屬性可以較好反映火山巖分布范圍；由于能量類屬性與振幅值平方密切相關，故能量類屬性同樣可以用來進行火山巖的識別，如總能量、平均能量等。

筆者研究選用均方根振幅屬性(以下簡稱 RMS)，從計算方法上看，振幅平均前先平方，所以RMS對振幅異常更加敏感[8-9]。在具體操作上，合適的時窗提取長度不僅有利于火山巖范圍的刻畫，更有利于火山口的識別，研究沿沙三段火山巖底部(0～20ms)提取了RMS信息(見圖4)。

結合對火山巖地震相反射特征的分析，可知火山口處、李天木斷層振幅較小，溢流相振幅較大。通過對RMS屬性圖分析，在J15-21井東可以清楚的識別出小區域范圍內同心圓狀低振幅異常區，即內圓為高振幅，外圓為低振幅屬性；同樣在Z78井、J24-22井處也可識別出近似同心圓狀低振幅異常區，這些區域與地震剖面上柱狀地震反射分布范圍一致，為火山口的平面分布位置。由圖4高振幅區可以預測，棗35斷塊的火山巖多分布于該區北部和東南部，北部火山巖振幅高于東南部火山巖振幅，其原因主要在于火山口處于構造高處(即研究區的東南方)，熔巖流動方向主要有北向、南西向2個方向，由于研究區東南高、西北低，熔巖主要沿著北-北西向溢流，最終導致北部火山巖厚度明顯大于西南部。

3 火山巖地震分頻信息分析

不同厚度的地質體有著不同的頻率特征，采用分頻解釋技術將火山巖地震反射在頻率域進行分頻振幅處理，以便確定哪個頻段的振幅能準確刻畫沙三段火山巖的平面分布規律[3]，可達到識別、劃分火山巖的橫向變化的目的。

結合地震相和地震屬性的分析結果，對研究區內3個柱狀地震反射、12個丘狀地震反射以及20個層狀地震反射進行了振幅譜(見圖5)分析，落實研究區火山口頂部平均頻率為20Hz，火山口內部頻率段為12～16Hz；溢流相頂部頻率段為18～25Hz；丘狀構造頂部為18～25Hz。經過上述統計可得知，火山巖體的頂部頻率多高于內部頻率，這是因為火山巖體頂部多具有高頻屏蔽作用，在火山巖頂界面會出現瞬時高頻，但由于火山巖地層對高頻成分吸收嚴重，使得火山巖內部具有低頻特性[5]。

圖5 LINE2154剖面上火山頸頻譜圖

圖6 棗35塊沙三段火山巖20Hz振幅譜圖

筆者選用目的層調諧體對火山巖表征進行預測。以沙三段火山巖頂部作為目的層段，時窗長度為30ms(1/2相位為15ms左右)，利用SpecDecomp得到分頻地震數據體，在0～80Hz的頻段內以1Hz為間隔觀察頻率切片，根據不同火山巖體的特定頻段找出其相對應的振幅圖(見圖6)。由此觀察，15～30Hz的振幅圖與實鉆火山巖分布范圍相符，而20Hz的振幅譜對于火山口的位置刻畫有一定幫助(該區地震資料的主頻為25Hz)。

4 棗35斷塊沙三段火山巖分布范圍預測以及巖相劃分

圖7 棗35塊沙三段火山巖巖相分布圖

通過地層對比了解該區不同井內火山巖變化特征，通過棗35斷塊油藏不同巖性的測井響應特征[5]，對研究區內缺少錄井資料的井進行火山巖統計，發現火山巖廣泛分布于沙三段，利用人工合成記錄建立時深關系，標定沙三段火山巖所在層位；其次，利用柱狀地震反射特征在平面構造圖上標出火山口；最后，利用振幅屬性(能量屬性)、分頻后的振幅圖對火山巖分布進行詳細刻畫。根據火山巖平面分布的結果，結合鉆井揭示的火山巖巖性特征，編制了棗35塊沙三段火山巖巖相分布圖(見圖7)。

棗35塊沙三段火山巖厚度受到多個火山口多期次噴發的影響，呈片狀分布(見圖7)，火山巖主要發育于研究區北部和東南部。溢流相分布范圍較廣，北部溢流相火山巖較東南部更厚。火山口主要分布于J15-21井東部，ZG53-5井北部以及J25-23井南方向，其中J25-23井南東方向火山口為研究區內主火山口，其周圍火山巖分布范圍較大，在主火山口地震剖面上觀察到明顯的火山頸，且其振幅呈同心圓環狀，故爆發相火山巖主要位于主火山口附近；小火山口分布于李天木斷層兩側，規模較小，火山口呈不規則條帶狀，故推測該噴發類型為裂隙式噴發。研究區內的李天木斷層對火山巖的控制作用十分顯著，是巖漿運移的重要通道；火山沉積相(凝灰巖等)則呈不規則狀分布于整個研究區內，出露范圍較小。試油數據表明，初步認為火山口周圍爆發相具有較好的儲集性能，但是由于面積較小不能作為該區主要儲層；該區溢流相火山巖體儲層物性較好。因此，溢流相區是下一步生產開發的有利區域。

5 結論

1)地震振幅屬性和分頻屬性能較好地預測棗園油田火山巖分布，并且能夠較準確確定火山口位置。

2)棗園油田棗35斷塊沙三段發育了3個火山巖體，其中對于該地區火山巖儲層分布范圍具有控制作用的主火山口分布于棗35斷塊東南部；該區北部火山巖厚度明顯大于西南部。

3)該區內溢流相火山巖是該地區主要的油氣儲集體，是下一步生產開發的有利區域；爆發相火山巖具有良好的儲集性能。

[參考文獻]

[1]王華崇，冉啟全，胡永樂. 大港棗園油田火成巖巖相[J]. 石油勘探與開發，2004， 31(5)：21-24.

[2]鄭亞斌，王延斌，冉啟全. 棗35斷塊火山巖儲層特征研究[J]. 石油學報，2006，27(4)：54-58.

[3]冉啟全，胡永樂，任寶生. 火成巖巖性識別方法及其應用研究——以大港棗園油田棗35塊火成巖油藏為例[J]. 中國海上油氣，2005，17(1)：25-30.

[4]湯浩哲，許永忠，劉永福，等. 地震屬性技術在哈拉哈塘地區火成巖識別中的應用[J]. 能源技術與管理，2012(4)：20-22.

[5]鄭亞斌，王延斌，韓德馨. 地震反演技術在火山巖儲集層預測中的應用——以棗35塊為例[J]. 新疆石油地質，2006，27(6)：746-748.

[6]張殿成，何德友，尚廣弟，等. 松遼盆地汪家屯東火山巖儲層預測[J]. 石油物探，2000，39(2)：36-43.

[7]陸基孟. 地震勘探原理[M].北京：石油大學出版社，1993.

[8]彭成敏，趙磊，吳擁軍，等. 文安斜坡沙河街組火山巖空間展布研究[J]. 科學技術與工程，2013，13(4)：980-984.

[9]王元君，王貴重，王建忠，等.地震屬性在英買力地區火成巖描述中的應用[J]. 石油物探，2007，46(1)：48-52.