渤海錦州A構造古近系地震儲層預測研究

于 茜， 邊立恩， 沈洪濤

(中海石油(中國)有限公司 天津分公司渤海石油研究院，天津 300459)

渤海錦州A構造古近系地震儲層預測研究

于 茜， 邊立恩， 沈洪濤

(中海石油(中國)有限公司 天津分公司渤海石油研究院，天津 300459)

渤海錦州A構造古近系東營組二段(E3d2)從上至下分為東二上段(E3d2u)和東二下段(E3d2L)，其中東二下段(E3d2L)一油組為辮狀河三角洲前緣沉積，地震剖面上廣泛發育疊瓦狀前積反射特征。首先通過地震相分析、古地貌恢復、分頻技術及正演模擬等分析了前積反射的成因, 認為該前積反射的形成主要反映了間歇型辮狀河三角洲前緣席狀砂順斜坡方向的前積作用,是由多個砂體沿著物源方向前積、疊置和遷移所形成的。在此基礎上，基于新三維地震資料，結合地震多屬性綜合分析、砂地比擬合等開展儲層預測，對前積砂體的平面展布特征進行刻畫，最終對錦州A構造東二下段儲層有利發育區做出預測，為該構造的油氣勘探評價提供了有利的借鑒。

疊瓦狀前積反射； 古地貌； 正演模擬； 分頻； 儲層預測； 地震屬性

0 研究區概況

錦州A構造位于渤海遼西北洼北端，該構造自1989年到1991年基于二維地震資料解釋成果共鉆探4口井，從已鉆井(圖1)油氣顯示情況可以看到，該區主要的成藏層系位于東二下亞段，各井鉆遇的砂體厚度普遍較薄，多數在1 m～10 m之間，儲層不發育，但油氣顯示很活躍。可見，儲層預測是制約該構造油氣勘探取得突破的最關鍵因素。2013年通過對該區新三維地震資料的精細解釋，發現東二下段上部廣泛發育前積型地震反射，該套前積反射有望成為該區尋找優質儲層的突破口。

1 地震相特征分析

1.1 層位標定

通過合成地震記錄制作，對遼西北洼錦州A構造東營組地層進行井震層位標定(圖2)，以建立該巖性段井震特征對應關系。從圖2可以看到，東二下段一油組為兩波谷夾一中-強振幅、連續寬相位波峰復波反射，該復波反射反映了多個薄砂層組合所形成的復波地震反射，復波的寬度和振幅與薄砂層組所在位置砂層的總厚度及單砂層厚度有關。如J5井東二下段一油組單砂層厚度和砂層總厚度較大，對應的反射波振幅相對較弱，復波的寬度較大，頻率明顯較低；相鄰的J1、J2井一油組單砂層厚度和砂層總厚度較小，對應的反射波振幅相對較強，復波特征不明顯，波形寬度較小，頻率明顯變高。上述地震響應特征在砂、泥巖薄互層模型地震正演剖面(圖3)中得到證實，即地震剖面上無法分辨單個薄層砂體，多個薄砂層的疊合體對應一個復合波形，復合波形的振幅和頻率與單砂層的厚度有密切關系。如單層砂層厚度為1 m～4 m時(如圖3(a)),薄砂層組對應一個低頻、強振幅單波反射，而單層薄砂層厚度增大到4 m～8 m時(如圖3(b))，砂層組對應一個振幅較弱、波形寬緩的復波反射。上述特征與研究區井震響應規律基本一致。可見，對于錦州A構造東二下段地層而言，低頻寬相位、弱振幅的復波反射特征代表儲層更加發育。

圖1 錦州A構造東營組油氣顯示連井對比圖Fig.1 The comparison of oil-gas shows in Dongying formation of Jinzhou A structure

圖2 錦州A構造井震標定結果Fig.2 Results of seismic-well calibration for Jinzhou A structure

圖3 不同厚度薄砂層模型及其正演結果對比Fig.3 The comparison of different thickness of thin sand model and its modeling results(a)單層砂層厚度為1 m～4 m;(b)單層薄砂厚度為4 m～8 m

1.2 物源方向識別

物源方向是沉積研究的一項重要內容，它在確定物源位置、物源性質、沉積作用與沉積體系、沉積物的搬運路徑及儲層的空間展布等方面起著重要的作用[1-2]，物源方向的識別是宏觀上分析儲層時空分布規律的重要依據。

1.2.1 古地貌分析

古溝谷是物源區向盆地深凹區輸送沉積物的通道，是沉積物搬運的直接證據，良好的輸砂通道保證前方沉積砂體的持續供應，控制著主力砂體的發育方向[3-5]。圖4是采用層拉平法恢復的東二下段沉積時的古地貌，從圖4可以看到，在靠近北東向邊界斷層的下降盤，規模性的古溝谷比較發育，因此由凸起向凹陷伸展并順著古溝谷的方向是該區主要的物源供給方向(圖4中紅色箭頭所指)。

1.2.2 分頻技術

圖4 東二下段沉積前古地貌圖Fig.4 The ancient landform map of the formation before E3d2L

分頻技術是一項新型的儲層研究技術，它通過將時間域的地震數據轉換到頻率域，對地震數據的頻率域特征進行表征，進而研究儲層的沉積特征或發育規律[6-7]。相似的沉積環境和地層組合，其在頻率域的響應往往相近。圖5是通過對三維地震數據體做頻譜分析后所提取的單頻歸一化振幅切片，根據分頻技術原理，低頻為厚層的響應，高頻為薄層的響應，從圖5可以看到，隨著頻率的降低(20 Hz～50 Hz),單頻切片的振幅譜峰值(紅-黃色)由凸起向凹陷區不斷擴展延伸，指示了沉積物的厚度由凸起向凹陷區逐漸增厚的過程，即古溝谷發育的位置儲層可能更加發育。同時，單頻切片的這種變化特征也反映了物源的方向(圖5中藍色箭頭所指)。

圖5 單頻振幅切片對比Fig.5 Comparison of single frequency slices(a)50 Hz;(b)40 Hz;(c)30 Hz;(d)20 Hz

1.3 前積地震反射特征分析

1.3.1 地震相特征

基于物源方向的認識，通過對錦州A構造精細地震解釋，發現順著物源方向該構造在東二下段一油組發育明顯的疊瓦狀前積地震反射(圖6(a))，各前積朵葉體首、尾具有明顯臺階狀反射特征，且上、下呈部分疊置關系，朵葉體內部呈低頻、弱振幅反射特征。而對于垂直物源的方向，在地震剖面上基本上看不到疊瓦狀的前積地震反射，剖面上主要以反映前積體內部的低頻、弱連續性、弱振幅反射特征為主(圖6(b))。

圖6 前積反射地震剖面Fig.6 Seismic profile of foreset reflection (a)順著物源方向；(b)垂直物源方向

1.3.2 正演模擬分析

為了研究前積反射特征的成因，對J5井的鉆井及測井資料進行統計得到一油組前積砂體的主要巖石物理參數如下：砂巖速度為3 012 m/s、密度為2 800 kg/m3,泥巖速度為2 520 m/s、密度為2 750 kg/m3,參考上述巖石物理參數，并結合圖6(a)所示的反射特征設計了地質模型，并采用聲波方程有限差分數值模擬技術進行地震正演模擬(圖7)。

圖7 地質模型及其地震正演剖面Fig.7 Geological model and its seismic forward modeling (a) 無前積砂體；(b)有前積砂體(傾角較大);(c)有前積砂體(傾角較小)

由圖7正演結果可見：當無前積砂體時，不存在前積反射特征(圖7(a))；當前積砂體有一定傾角且砂體較厚時，出現與圖6(a)相似的疊瓦狀前積反射特征(圖7b)；當前積砂體傾角較小、砂體較薄時，前積反射特征不明顯(圖7(c))，由正演結果可以看到，前積反射的形成需要在一定物源供給、斜坡地形和較強水動力條件下，在順古水流方向上容易形成前積砂體的疊瓦狀前積反射特征，圖6所示的前積反射代表了多個砂體順斜坡方向的前積、疊置和側向遷移作用，是優質儲層發育的指示。因此，對于錦州A構造，儲層預測的重點即在于識別或刻畫圖6所示的前積反射的平面展布特征。

2 儲層預測

2.1 地震多屬性分析

在新三維地震資料基礎上，利用地震屬性對東二下段一油組發育的低頻、弱反射前積砂體進行平面展布范圍預測，屬性提取時的時窗頂、底分別為東二下段一油組的頂、底面。本次屬性分析提取了能量類、頻率類、相位類、幾何類、統計類等五大類的多種屬性，由于每一種地震屬性只能從某一個角度反映地震波的特征，而實際的地震波常常是地下多種地質因素綜合響應的結果[8-13]，因此利用屬性交會的方法可以進一步降低屬性解釋的多解性。結合研究區的儲層地質條件及地震反射特征，本次交會分析通過對初步篩選出的十種屬性進行了交會統計(圖8)。圖8中，橫坐標列出的十種屬性與縱坐標列出的十種屬性及其順序完全相同，以交會圖的對角線呈軸對稱分布。由圖8可以看到，同一種屬性(橫、縱坐標屬性相同)在交會圖的對角線上呈對稱分布，且交會成一條直線，表明兩者完全相關(圖8中黑色箭頭所指)。任何兩種不同的屬性在交會圖上均形成一個交會面，其中反映共性地質信息的點則會在交會面上收斂、聚集，呈朵狀或錐狀分布(如圖8中藍色箭頭所指)，交會面上的點越收斂、聚集，說明反映的共性信息越多，越為有效的屬性。

圖9 錦州A構造東二下段一油組地震屬性對比Fig.9 Comparison of seismic attributes of No.1 oil group in E3d2L of Jinzhou A structure(a)弧長；(b)主頻率；(c)瞬時相位

結合交會分析結果，選取交會圖上屬性點更為聚集、收斂并呈朵狀分布的弧長屬性、主頻率屬性及瞬時相位屬性來從不同角度描述儲層的特征及發育規律：

1)弧長屬性(圖9(a))，由于研究區的前積反射主要表現為低頻寬相位特征，因此其對應的弧長屬性就表現為低值，即屬性的強弱直接反映了該區前積反射即儲層的發育程度。屬性分析結果與已鉆井的鉆遇情況吻合度非常高，即J1、J2、J3井儲層不發育，弧長屬性上表現為藍色的高值區，J5井儲層相對發育，屬性上表現為紅色的低值區。 另外，由該屬性分布結果可以看到，圖9(a)中粉色橢圓圈所示的構造高部位儲層可能會更加發育。

2)主頻率屬性(圖9(b))，紅色代表低頻，即儲層比較發育，厚度比較大；藍色代表高頻，指示儲層不發育,與已鉆井情況也是吻合的。

3)瞬時相位屬性(圖9(c))所示，該屬性對儲層邊界或橫向連通性的刻畫比較明顯，屬性中不連續的地方，基本反映了儲層的邊界范圍。

綜合以上屬性分析結果，認為優選出的3種屬性從不同角度揭示了該區的儲層發育規律，在構造北部邊界斷層下降盤的構造高部位，有望鉆遇到優質的辮狀河三角洲前積砂體。

2.2 砂地比預測

基于上述弧長屬性與已鉆井揭示的儲層發育程度的高吻合性，利用研究區4口已鉆井在目的層東二下段一油組、二油組及三油組的弧長屬性與其相應的砂地比統計結果進行交會分析(圖10)，發現二者呈較好的線性相關性，并以此進行了線性擬合，擬合結果為y=-0.029x+57.71，其中y代表砂地比，x為弧長屬性值。

圖10 弧長屬性與砂地比線性擬合圖Fig.10 The linear fitting chart of Arc length attribute and sand ratio

最終利用弧長屬性的平面分布及該擬合公式，對整個研究區的砂地比分布進行了平面預測(圖11)，該圖也同樣揭示了構造北部儲層的富集特征。

圖11 錦州A構造東二下段一油組砂地比分布Fig.11 Sand ratio distribution of No.1 oil group in E3d2L of Jinzhou A structure

3 結論

1)錦州A構造東營組東二下段(E3d2L)一油組廣泛發育的疊瓦狀前積反射，是由辮狀河三角洲前緣多個砂體前積、疊置和遷移所形成的，前積反射特征發育的位置指示了優質儲層的發育區。

2)儲層預測結果與已鉆井吻合度高，有效揭示了該構造儲層的空間發育規律，認為構造北部邊界斷層下降盤的構造高部位為優質儲層發育區，可能會鉆遇到厚層辮狀河三角洲前積砂體，具有較大的勘探評價潛力。儲層預測結果，為該區的油氣勘探提供了借鑒，也進一步降低了井位部署的風險，目前依據儲層預測結果部署的井位已經獲得通過，有望帶動該區迎來油氣勘探新的發現。

[1] 王冰. 蘇里格氣田盒8儲層沉積物物源方向探討[J].石油化工應用，2010，29(11)：38-41.

WANG B. A discussion on source direction from the He8 reservoir in the Sulige gas field[J].Petrochemical Industry Application, 2010,29(11):38-41.(In Chinese)

[2] 張曉寶，何海清，陳國俊. 吐哈盆地草南1井中侏羅統物源方向綜合判識及其研究意義[J].沉積學報，1997，15(4)：109-114.

ZHANG X B, HE H Q, CHEN G J. A comprehensive judgement of the source direction of the middle jurassic of Chaonan 1 well in the Tuha basin and its research significiance[J].Acta Sedimentologica Sinica, 1997,15(4):109-114. (In Chinese)

[3] 趙俊興，陳洪德，時志強. 古地貌恢復技術方法及其研究意義——以鄂爾多斯盆地侏羅紀沉積前古地貌研究為例[J]. 成都理工學院學報，2001，28(3)：1589-1596.

ZHAO J X, CHEN H D, SHI Z Q. The way and implications of rebuilding palaeogeomorphology-Taking the research of palaeogeomorphology of the Ordos Basin before Jurassic deposition as example[J].Acta Sedimentologica Sinica,2001,28(3):1589-1596. (In Chinese)

[4] 張宇焜，胡曉慶，牛濤，等. 古地貌對渤海石臼坨凸起古近系沉積體系的控制作用[J]. 吉林大學學報(地球科學版)，2015，45(6)：42-47.

ZHANG Y K,HU X Q,NIU T,et al. Controlling of paleogeomorphology to paleogene sedimentary system of Shijiutuo Uplift in Bohai basin[J].Journal of Jilin University(Earth Science Edition),2015, 45(6):42-47. (In Chinese)

[5] 劉峰，裴健翔，汪洋，等. 古地貌對海底扇沉積過程的控制及與油氣富集的關系—以鶯歌海盆地東方區黃流組一段為例[J]. 中國海上油氣，2015,27(4)：37-46.

LIU F ,PEI J X,WANG Y,et al. Palaeogeomorphologic control on sedimentary process of submarine fans and hydrocarbon accumulation:a case study of Member 1 of Huangliu Formation in DF area,Yinggehai basin[J].China Offshore Oil and Gas, 2015,27(4):37-46. (In Chinese)

[6] 朱定，孫素琴，蔣進勇. 分頻解釋技術在塔河油田儲層預測中的應用[J].石油與天然氣地質，2010，31(1)：49-53.

ZHU D,SUN S Q,JIANG J Y. Application of frequency division interpretation technique to reservoir prediction in Tahe oilfield[J]. Oil & Gas Geology,2010,31(1):49-53. (In Chinese)

[7] 邊立恩，于茜，韓自軍，等. 基于譜分解技術的儲層定量地震解釋[J].石油與天然氣地質，2011，32(5)：718-723.

BIAN L E,YU Q,HAN Z J,et al.Quantitative seismic interpretation of reservoirs based on the spectral decomposition technique[J].Oil & Gas Geology, 2011,32(5): 718-723. (In Chinese)

[8] 竇松江，蔡明俊，陳金鳳，等. 地震屬性分析在河道砂體內部構型研究中的應用[J].石油地質與工程，2009，23(5)：46-49.

DOU S J,CAI M J,CHEN J F,et al. The application of the seismic attribution analysis on the channel sand body architecture[J].Petroleum Geology and Engineering, 2009,23(5):46-49. (In Chinese)

[9] 武麗，于文芹，楊立強，等. 地震屬性分析技術在波斯灣盆地中部S區塊儲層預測中的應用[J].石油地球物理勘探，2008，43(5)：557-561.

WU L, YU W Q, YANG L Q,et al. Application of seismic attributes analytic technique to reservoir prediction in some block of central part of Persian Gulf basin[J]. Oil Geophysical Prospecting,2008,43(5): 557-561. (In Chinese)

[10] 趙斌. 地震屬性約束定量預測儲層參數—以渤海D油田4號砂體為例[J].石油地質與工程，2012，26(2)：31-34. (In Chinese)

ZHAO B.With the restraint of seismic attributes to predict reservoir parameters quantificationally—an example from No.4 sand body of D oilfield in Bohai[J]. Petroleum Geology and Engineering, 2012,26(2):31-34. (In Chinese)

[11] 邊立恩，屈勇，于茜，等. 渤海BZ油田儲層定量預測研究[J].石油地質與工程，2016，30(3)：77-80. (In Chinese)

BIAN L E,QU Y,YU Q,et al. Study on reservoir quantitative prediction of BZ oilfield in Bohai sea[J].Petroleum Geology and Engineering, 2016,30(3):77-80. (In Chinese)

[12] 姚威，吳沖龍，史原鵬，等. 利用地震屬性融合技術研究洪浩爾舒特凹陷下白堊統沉積相特征[J].石油地球物理勘探，2013，48(4)：634-642.

YAO W,WU C L,SHI Y P,et al. Sedimentary characteristics research on lower Cretaceous in Honghaoershute Depression using seismic multi-attribute fusion[J].Oil Geophysical Prospecting, 2013,48(4):634-642. (In Chinese)

[13] 李密，李少華，王紅賓，等. 多屬性概率神經網絡技術在ML油田巖性油氣藏預測中的應用[J].石油地質與工程，2013，27(3)：50-52.

LI M,LI S H,WANG H B,et al. Application of multi-attribute probabilistic neural network inversion in lithologic reservoir hydrocarbon prediction[J].Petroleum Geology and Engineering, 2013,27(3):50-52. (In Chinese)

StudyonpaleogenereservoirpredictionofJinzhouAstructureinBohai

YU Qian, BIAN Li'en, SHEN Hongtao

(Bohai Oilfield Research Institute, Tianjin Branch of CNOOC Ltd., Tianjin 300459, China)

Member two of Dongying formation of Paleogene can be divided into two lithologic segments, the upper part(E3d2) and the lower part(E3d2) from top to bottom in Jinzhou A structure. The first oil group of the lower part is of the braided river delta front, where the characteristics of imbricated foreset reflections is very obvious in the seismic profile. Firstly, the origin of the progradational reflection is analyzed by seismic facies analysis, the recovery of the ancient landform, the spectral decomposition technique and the forward modeling. This concludes that the progradational reflection shows progradation action of sheet sand stones of the intermittent braided river delta front along the direction of slope and is composed of progradation, superposition and lateral migration of multiple thin sand bodies along the direction of provenance. On this basis, reservoir prediction is carried out by multiple seismic attributes comprehensive analysis and sand ratio fitting based on the new 3D seismic data and the later distribution characteristics of the foreset sand bodies is described. Finally this paper forecasts the favorable development zone of the reservoir in the lower part of member two of Dongying formation which provides a favorable target for oil-gas exploration and evaluation of the structure.

imbricated foreset reflections; ancient landform; forward modeling; spectral decomposition; reservoir prediction; seismic attribute

2016-11-24 改回日期： 2017-06-17

國家重大科技專項(2011ZX05023-002)

于茜(1983-)，女，碩士，工程師，主要從事地震資料解釋與儲層預測研究工作，E-mail：120776161@qq.com。

1001-1749(2017)06-0808-08

P 631.4

A

10.3969/j.issn.1001-1749.2017.06.14