灰质背景下基于敏感曲线的储层预测方法——以孤北洼陷沙三中亚段为例

乌洪翠

（中国石化胜利油田分公司勘探开发研究院）

灰质背景下基于敏感曲线的储层预测方法
——以孤北洼陷沙三中亚段为例

乌洪翠

（中国石化胜利油田分公司勘探开发研究院）

针对灰质泥岩影响储层预测的问题，以孤北洼陷沙三中亚段为例，提出一种有效剔除灰质泥岩影响的储层预测方法。首先，通过对比常规测井曲线，认为声波时差、密度、自然伽马、自然电位以及波阻抗曲线对区分灰质泥岩和砂岩均不敏感；其次，结合岩石物理学特征及测井数据分析，表明泥质含量曲线对灰质泥岩和砂岩区分效果明显，可以作为区分二者的敏感曲线；最终，提出基于优选的敏感曲线，通过建立的等时地层格架，利用地质统计学反演技术提高储层预测的精度。预测结果表明，该方法能够有效剔除灰质泥岩的影响，准确刻画扇体形态，且与实钻结果具有较高的吻合度，达到了勘探部署的要求。

灰质泥岩；储层预测；敏感曲线；孤北洼陷；泥质含量

孤北洼陷经过近 50年的勘探，目前整体处于高勘探程度阶段，中深层岩性砂体成为最具潜力的勘探目标之一[1-3]。在以浊积砂体、扇三角洲前缘扇体为代表的岩性油藏勘探中，如何精确描述储层是勘探成败的关键[4-6]。孤北洼陷西南部沙三中亚段广泛发育的灰质泥岩与砂岩在测井响应及地震反射特征方面的相似性，影响了储层预测的精度。

储层预测中去除灰质泥岩影响的方法目前主要集中在两个方面：一是利用地震属性，通过提取灰质泥岩集中发育段相应频率的数据体或能量半衰时实现对灰质泥岩的剔除[7-9]；二是采用拓频的方法，通过子波分解技术重构和提高原始地震资料的品质，一定程度上达到去除灰质泥岩影响的目的[10-11]。由此可见，剔除灰质泥岩的影响多集中于对地震数据的处理，没有做到将地质、地震、测井、录井等信息有效地综合利用。本文旨在通过优选一条能够有效区分灰质泥岩和砂岩的测井敏感曲线，借助精细等时地层格架[12]和地质统计学反演技术，形成一套融合测井、地震、地质认识的储层预测技术，为灰质背景下的岩性油藏勘探提供技术支撑。

1 基本思想

孤北洼陷位于济阳坳陷沾化凹陷东北部，西与埕东凸起相邻，南与孤岛凸起相接，东、北分别毗邻长堤潜山、桩西潜山。受古地貌及新生代断陷作用的影响，孤北洼陷形成了五号桩洼陷、孤北西次洼和孤北低隆起“两洼一隆”的构造格局 (图1)。孤北洼陷作为沾化凹陷的油气富集区,发育多套有利储盖组合，目前已发现了沙河街组、东营组和馆陶组等多套含油层系，其中沙三段为主力含油层系。沙三段油藏类型以构造油藏和岩性油藏为主，累计探明石油地质储量5000×104t以上。

图1 孤北洼陷构造纲要图

岩性油藏勘探中对厚度小、横向变化快的岩性体进行描述和预测时，受分辨率的限制，仅凭地震资料无法获得理想的预测结果。特别是在受灰质泥岩干扰的地区，更需要井震信息结合，充分利用测井资料井点处储层与围岩的某种测井响应特征差异[13]，以及地震资料空间上的波形属性响应差异，综合预测有利储层的分布。首先，测井曲线具有很高的纵向分辨能力和精度，优选区分灰质泥岩和砂岩敏感的曲线参与储层预测，是识别有利储层的基础；其次，在宏观地质规律的指导下，按各砂组占沙三中亚段总厚度比例，建立各砂组精细的等时地层格架，从而保证储层预测横向的等时性和纵向的针对性；最终，在等时地层格架约束下，通过地质统计学反演将稀疏脉冲波阻抗反演结果和敏感测井曲线有机融合，用测井数据弥补波阻抗数据纵向分辨率的不足，提高储层预测精度（图2）。

图2 储层预测流程图

2 敏感测井曲线优选

2.1 测井曲线敏感性分析

以研究区z394井为例，自然伽马、密度、声波时差、电阻率曲线区分灰质泥岩和砂岩的效果不明显，灰质泥岩和砂岩测井响应特征相似，但二者与泥岩不同（图3）。统计表明：砂岩波阻抗值整体位于灰质泥岩区间内；泥岩的密度、声波时差曲线幅值与灰质泥岩、砂岩虽有部分重叠，但相较灰质泥岩和砂岩，区分效果较为明显（表1）。自然电位曲线通常可以区分泥岩和砂岩，但容易受地层水矿化度、钻井液电阻率、钻井液侵入、含油性等因素的影响。研究区位于海陆交互地带，钻井液电阻率变化大，对自然电位曲线造成影响。

图3 z394 井不同岩性测井响应特征

表1 不同岩性的测井曲线幅值

在排除掉声波时差、密度、自然伽马、电阻率和自然电位曲线后，急需优选一种能表征灰质泥岩和砂岩差别的曲线。从岩石物理学特征考虑，灰质泥岩作为泥岩的一种，必然在泥质含量方面与砂岩存在区别。表征岩性差异的曲线除了常用的自然伽马外，还有泥质含量曲线。

通过单井泥质含量曲线的对比发现，泥质含 量 曲 线 能 够 有 效 区 分 3117～3124m 段 砂 岩 和3133～3135m、3136～3138m、3140～3148m段灰质泥岩、灰质油泥岩（图4）。通过对多口井泥质含量曲线统计发现，砂岩的泥质含量低于 30%，灰质泥岩为 66%～92%，泥岩高于 79%；泥岩和灰质泥岩存在重叠区，但二者与砂岩分界线明显。初步分析认为，泥质含量曲线能够有效区分研究区内灰质泥岩和砂岩,可以作为区分灰质泥岩和砂岩的敏感曲线。

图4 泥质含量曲线区分灰质泥岩和砂岩效果图

2.2 测井曲线计算

泥质含量曲线通常由自然伽马曲线通过公式计算得到。首先，用自然伽马相对幅度的变化计算泥质含量指数IGR：

式中 GR目的——目的层的自然伽马值；

GRmax——纯泥岩层的自然伽马值；

GRmin——纯砂岩层的自然伽马值。

通常IGR值介于 0～1，用下式将IGR转化为泥质含量Vsh：

式中 GCUR——希尔奇指数，在古近系—新近系取值 3.7，前古近系取值 2.0[14]。

受母源成分的影响，不同层段的自然伽马曲线幅度不同，在实际计算过程中应根据沉积旋回，分段计算泥质含量，最大程度提高泥质含量计算精度。

2.3 测井曲线可靠性分析

泥质含量曲线作为一种非实测得到的曲线，需要验证其区分储层与围岩的适用性和可靠性。首先，泥质含量曲线是通过其他曲线计算得到，其计算数据源是实测的，且数据经过分析优选，真实可靠；其次，泥质含量曲线是通过经典公式计算，计算过程中还会考虑不同层段的泥岩基线值，且受录井资料约束；最终，除了前述数据源和计算过程外，最重要的是计算结果要经得起实钻资料的检验，能够准确刻画灰质泥岩与砂岩间的差异。

为此，沿物源方向分别选取分属扇三角洲不同亚相的 3 口井（z392 井、z397 井、z26 井）进行对比。z392 井位于扇三角洲前缘，靠近物源方向，3014.0～3023.5m、3036.0～3045.0m 段发育两套细砂岩，与 3054.0～3080.0m 段发育的多套灰质泥岩区分效果明显；z397 井位于扇三角洲前缘，远离物源，泥质含量曲线存在一定程度齿化，但灰质油泥岩、泥岩与砂岩的差异性明显；z26 井仅发育 2.0m/1 层细砂岩，从区分效果来看，z26 井不如 z392 井、z397井区分效果明显，主要因为 z26 井位于前扇三角洲，岩性以砂泥岩互层为主，储层岩性多为灰质砂岩、泥质砂岩，随着砂岩中泥质、灰质成分的增多[15]，与泥岩的差异性变弱，但根据泥质含量曲线依然能够识别有利储层（图5）。总之，泥质含量曲线能够有效区分砂岩与灰质泥岩，只是在扇体的不同部位区分效果存在一定差异性。扇三角洲前缘区分效果明显，前扇三角洲则相对较差。

3 实现过程

3.1 精细等时地层格架

切片方法通常有3种方式：时间切片、沿层切片、地层切片[16-17]，针对不同的地层结构采用不同的建立方式。沙三中亚段自上而下分为5个砂组，各砂组厚度较薄（20m 左右）。分析认为研究区沙三中亚段沉积较薄、连续沉积、无超剥，但地层厚度差异大、倾角大，因此时间切片、沿层切片不能保证储层预测的等时性，而依据地震反射特征追踪各砂组底界具有较大的不确定性且精度不够。沙三中亚段各砂组纵向地层厚度变化均匀，为此，采用地层切片技术基本原理，统计各砂组占沙三中亚段总厚度比例 (表2)，纵向上按比例求取各砂组底界，横向上在地震数据的约束下内插外推，建立各砂组等时地层格架。通过该方法建立的地层格架，能够最大限度避免沿层切片和时间切片所造成的穿时现象，克服大尺度地层切片造成的预测无针对性等问题，使地层格架更接近真实的地质结构，且与预测砂体厚度更匹配。

图5 泥质含量曲线区分灰质泥岩与砂岩效果验证图

表2 沙三中亚段各砂组厚度所占比例统计表

3.2 地质统计学反演

地质统计学反演是提供一个在某种概率条件下，既满足数据的地质统计学特征又满足地质、测井和地震信息的三维储层参数概率模型，其适用于井距小于一个砂体长度的高勘探程度地区。目前研究区钻井密度为 0.9 口 /km2左右，适用地质统计学反演的方法。

地质统计学反演的主要步骤是求取空间变差函数[18]，包括求取横向变差函数和纵向变差函数。横向变差函数主要是从波阻抗数据体中计算得出，纵向变差函数通过对上述敏感测井曲线的计算得到。

地质统计学反演的技术关键是各种统计学参数的选取。变差函数中主要涉及变程和基台值两个参数；计算变差函数后还需要拟合一个理论模型来表征储层的空间结构，研究区选取的是指数模型；反演还会涉及砂体的宽度、长度、厚度等参数，这些参数的选取主要是通过地质认识、实钻井和统计得到。

研究区以波阻抗数据和优选的敏感测井曲线作为输入，利用地质统计学反演技术，得到反映储层变化的泥质含量数据体。

4 应用效果

利用上述方法对研究区沙三中亚段进行了储层预测。以沙三中亚段4砂组为例，从预测结果平面图上看扇体特征刻画更加清晰，较之前均方根振幅属性具有较大改善（图6）。孤岛凸起物源的扇体发育在孤北低隆起西翼，发育 gb342、gb35、gb107、z393 等4个扇体，以z393扇体规模最大，扇体呈北西向展布，可延伸至z26井到z898 井之间；z26井前端发育浊积扇体，浊积扇体与 z331 块发育的埕东凸起物源的近岸水下扇界线明显。

剖面分析认为，未参与储层预测井与实钻结果吻合度较高，单砂体及砂体间预测准确，预测精度在85% 以上，能够达到勘探部署的要求。

总体而言，运用该方法所获得的储层预测结果与先前单一的均方根属性预测结果相比较，在反映扇体形态方面具有很大改进，扇体展布、延伸距离、扇体间相互关系均得到了很好刻画。

图6 沙三中亚段 4 砂组储层预测平面图对比

5 结论

（1）由自然伽马曲线出发通过公式计算得到泥质含量曲线，通过校验计算过程和结果认为，该曲线能够有效区分研究区灰质泥岩和砂岩，且效果理想。

（2）在孤北洼陷应用基于敏感曲线的储层预测方法，达到了在灰质背景下精细描述储层并完整、准确反映扇体形态的目的，储层的展布规律与地质认识、实钻结果吻合度较高，完善了储层预测方法，实现了多尺度地质资料的综合应用。

（3）在该方法描述储层的基础上，结合构造、烃源岩等地质认识和已钻井分析，共描述有利圈闭22 个，其中在 z26 井北部、z48 井西部为探索浊积扇体的含油性部署探井两口，均取得良好勘探效果。

参考文献

[1] 秦永霞，尚延安，方旭庆，范云 .孤北洼陷西斜坡沙三段中亚段储集体特征 [J].油气地质及采收率，2003,10( 增刊 1):10-11. Qin Yongxia,Shang Yanan,Fang Xuqing,Fan Yun.Reservoir characteristics of Sha-3 central Member in Gubei subsag west slope [J]. Petroleum Geology and Recovery Efficiency, 2003,10(Supp.1):10-11.

[2] 王永诗，庄文山 .孤北洼陷沙三、四段沉积体系分析 [J].油气地质与采收率，2003,10(3):7-10. Wang Yongshi, Zhuang Wenshan. Analysis on sedimentary system in Es3and Es4of Gubei subsag [J]. Petroleum Geology and Recovery Efficiency, 2003,10(3):7-10.

[3] 杨剑萍，赵卫卫，姜在兴 .沾化凹陷孤北油田古近系沙三段扇三角洲 沉 积 特 征 及 油 气 储 层 意 义 [J]. 石 油 与 天 然 气 地 质，2003,24(2): 157-161. Yang Jianping, Zhao Weiwei, Jiang Zaixing. Sedimentary characteristics and reservoir significance of Eogene Es3fan delta in Gubei oilfield , Zhanhua depression [J]. Oil & Gas Geology, 2003,24(2):157-161.

[4] 王永诗，王勇，朱德顺，丁桔红，尚冰，朱家俊 .东营凹陷北部陡坡带砂砾岩优质储层成因 [J].中国石油勘探，2016,21(2):28-36. Wang Yongshi, Wang Yong, Zhu Deshun, Ding Juhong, Shang Bing, Zhu Jiajun. Genetic mechanism of high-quality glutenite reservoirs at the steep slope in northern Dongying sag [J]. China Petroleum Exploration, 2016,21(2):28-36.

[5] 张志虎，孙玉红，金力钻，李松林，黄晓 .致密砂岩气层测井评价新方法研究 [J].特种油气藏，2015,22(5):46-49. Zhang Zhihu, Sun Yuhong, Jin Lizuan, Li Songlin, Huang Xiao.A new method of logging interpretation for tight sandstone gas zones [J]. Special Oil & Gas Reservoirs, 2015,22(5):46-49.

[6] 成大伟，袁选俊，周川闽，谭聪，汪梦诗 . 测井岩性识别方法及应用——以鄂尔多斯盆地中西部长 7 油层组为例 [J].中国石油勘探，2016,21(5):117-126. Cheng Dawei, Yuan Xuanjun, Zhou Chuanmin, Tan Cong, Wang Mengshi. Logging-lithology identification methods and their application: a case study on Chang 7 Member in central-western Ordos Basin, NW China [J]. China Petroleum Exploration, 2016,21(5):117-126.

[7] 林中凯，姜蕾 .临南洼陷沙三段浊积砂体地震特征与精细描述 [J].石油天然气学报，2012,34(4):73-77. Lin Zhongkai, Jiang Lei. Seismic characteristics and fine description of turbidite sandbody in the Third Member of Shahejie Formation in Linnan sag [J]. Journal of Oil and Technology Gas, 2012,34(4):73-77.

[8] 张营革 .能量半衰时属性在浊积岩储层预测中的应用 [J].石油物探，2013,52(6):662-668. Zhang Yingge. The applicaion of energy half-time attribute to turbidite reservoir prediction [J]. Geophysical Prospecting for Petroleum, 2013,52(6):662-668.

[9] 刘文辉，潘和平，李健伟，赵景 .鄂尔多斯盆地大牛地气田铝土质泥岩储层的测井评价 [J].天然气工业，2015,35(5):24-30. Liu Wenhui, Pan Heping, Li Jianwei, Zhao Jing. Well logging evaluation on bauxitic mudstone reservoirs in the Daniudi Gasfield, Ordos Basin [J]. Natural Gas Industry, 2015,35(5):24-30.

[10]于正军 .灰质背景下浊积岩储层地震响应特征及识别方法——以东营凹陷董集洼陷为例 [J].油气地质与采收率，2014,21(2):95-97. Yu Zhengjun. Seismic response characteristics and recognition method of turbidity under carbonate depositional environment-a case in Dongji sag of Dongying sag [J]. Petroleum Geology and Recovery Efficiency, 2014,21(2):95-97.

[11]王惠勇，陈世悦，张云银，谭明友，于景强 .东营凹陷浊积岩优质储层预测技术 [J].石油地球物理勘探，2014,49(4):776-782. Wang Huiyong, Chen Shiyue, Zhang Yunyin, Tan Mingyou, Yu Jingqiang. Turbidite high-quality reservoir prediction in Dongying depression [J]. Oil Geophysical Prospecting, 2014,49(4):776-782.

[12]谷志猛，王德英，王明臣，杨海风，张藜 .海上岩性油气藏勘探的关键思路与技术方法——以渤海海域石臼坨凸起中段为例 [J].中国石油勘探，2015,20(4):54-62. Gu Zhimeng, Wang Deying, Wang Mingchen, Yang Haifeng, Zhang Li. Key exploration approach and techniques for offshore lithologic reservoirs: a case study on the middle section of Shijiutuo salient in Bohai sea area [J]. China Petroleum Exploration, 2015,20(4):54-62.

[13]方雪芹，刘德菊，田军，王建民，高君 .测井资料在松辽盆地岩性储层预测中的应用 [J].石油物探，2007,46(4):496-500. Fang Xueqin, Liu Deju, Tian Jun,Wang Jianmin, Gao Jun. Application of logging data in synthetically predicting lithologic reservoir of northern Songliao Basin [J]. Geophysical Prospecting for Petroleum,2007,46(4):496-500.

[14]王贵文，郭荣坤 .测井地质学 [M].北京 :石油工业出版社，2000. Wang Guiwen, Guo Rongkun. Logging geology [M]. Beijing: Petroleum Industry Press, 2000.

[15]廖 飞 燕， 马 庆 林， 孙 新 销， 程 刚， 曾 文 光， 卢 政 环， 等 .北 部 湾盆地福山凹陷古近系流沙港组岩性油气藏勘探 [J].中国石油勘探，2015,20(2):43-50. Liao Feiyan, Ma Qinglin, Sun Xinxiao, Cheng Gang,Zeng Wenguang,Lu Zhenghuan, et al. Study of lithologic reservoir of Paleogene Liushagang Formation in Fushan depression of Beibu Bay Basin [J]. China Petroleum Exploration, 2015,20(2): 43-50.

[16]张义娜，朱筱敏，刘长利 .地震沉积学及其在中亚南部地区的应用 [J].石油勘探与开发，2009,36(1):74-79. Zhang Yina, Zhu Xiaomin, Liu Changli. Seismic sedimentology and its application in south of Central Asia [J]. Petroleum Exploration and Development, 2009,36(1):74-79.

[17]王江，林东成，李鹏 ,杨微 .地层切片技术在复杂勘探区岩相识别中的应用 [J].石油天然气学报，2010,32(6):62-65. Wang Jiang, Lin Dongcheng, Li Peng, Yang Wei. Application of stratal slicing in lithofacies identification of complex exploration area [J]. Journal of Oil and Gas Technology, 2010,32(6):62-65.

[18]王雅春，王璐 .地质统计学反演在杏北西斜坡区储层预测中的应用 [J].地球物理进展，2013,28(5):2554-2560. Wang Yachun, Wang Lu. Application of geostatistical inversion to reservoir prediction in western slope of the northern Xingshugang oilfield [J]. Progress in Geophysics, 2013,28(5): 2554-2560.

A reservoir prediction method based on sensitivity curve in the background of lime mudstone: a case study on middle Es3in Gubei subsag

Wu Hongcui
(Exploration and Development Research Institute, Sinopec Shengli Oilfi eld Company)

In this paper, the middle section of the third member of Shahejie Formation (middle Es3) in the Gubei subsag was taken as an example to investigate the effects of lime mudstone on reservoir prediction. A reservoir prediction method which can eliminate such effects was proposed. Firstly, comparisons on conventional logs indicate that acoustic (AC), density (DEN), gamma-ray (GR), spontaneous potential (SP) and wave impedance (WI) logs are not sensitive in distinguishing lime mudstone from sandstone. Secondly, petrophysical characteristics and logging data were analyzed, revealing that shale content curve play a remarkable role in distinguishing lime mudstone from sandstone, so it can be used as the sensitivity curve for lime mudstone and sandstone distinguishing. Finally, it was proposed to improve the reservoir prediction precision utilizing geostatistical inversion technology based on the preferred sensitivity curve and the established isochronous stratigraphic framework. The prediction result of this method is in higher coincidence with the actual drilling result. It is demonstrated that this method can eliminate the effects of lime mudstone effectively, characterize the morphology of fans accurately and meet the requirements of exploration deployment.

t: lime mudstone, reservoir prediction, sensitivity curve, Gubei subsag, shale content

P631.8

A

10.3969/j.issn.1672-7703.2017.03.007

国家科技重大专项“渤海湾盆地精细勘探关键技术”(2011ZX0506-003)。

乌洪翠 (1982-)，女，山东聊城人，硕士，2008 年毕业于中国石油大学（华东），工程师，现主要从事特殊岩性油气藏储层评价工作。地址 :山东省东营市聊城路 2 号胜利油田勘探开发研究院，邮政编码：257015。E-mail:wuhongcui.slyt@sinopec.com

2015-11-17；修改日期：2017-03-16