引用本文
郑卓, 王建立, 夏良冰, 刘广明, 张弘. 基于高斯混合模型聚类分析的解释方法在埕北油田气顶区混合油气层的应用[J]. 录井工程, 2022,33(3): 73-77.
ZHENG Zhuo, WANG Jianli, XIA Liangbing, LIU Guangming, ZHANG Hong. Application of interpretation method based on Gaussian Mixture Model cluster analysis to mixed oil and gas reservoirs of gas cap area in Chengbei Oilfield[J]. Mud Logging Engineering, 2022,33(3): 73-77.
Doi:10.3969/j.issn.1672-9803.2022.03.012  
Permissions
Copyright©2022, 《录井工程》杂志社
《录井工程》杂志社 所有
基于高斯混合模型聚类分析的解释方法在埕北油田气顶区混合油气层的应用
郑卓①, 王建立②, 夏良冰①, 刘广明①, 张弘①
①中海油能源发展股份有限公司工程技术分公司
②中海石油(中国)有限公司天津分公司

作者简介： 郑卓 工程师,1992年生,2014年毕业于中国地质大学(北京)资源勘探工程(能源)专业,主要从事海上油气勘探开发的工程地质及数字化应用工作。通信地址:300452 天津市滨海新区滨海智谛山。电话:(022)25800043。E-mail:zhengzhuo2@qq.com

收稿日期: 2022-07-29
摘要

埕北低凸起油气成藏构造是受燕山、喜山期多次区域性抬升影响而形成,埕北油田东二上亚段构造高部位区域气层发育,经过二十余年生产开采,气顶区已非原始气层状态,气藏规模变小且分布不均,常规录井、测井方法难以判别气油界面。在油田综合调整过程中,为了确保水平井着陆在纯油层中,规避混合油气层对产能的影响,收集埕北油田12口井的气测数据,运用高斯混合模型进行数据关联性分析,得出关键分离系数 α值,将分离系数与经验公式相结合,建立了适用于埕北油田气顶区的混合油气层解释方法,其在钻井过程中应用能实时判别并解释混合油气层,从而指导水平井准确着陆。

关键词: 高斯混合模型; 聚类分析; 混合油气层; 分离系数; 埕北油田
中图分类号:TE132.1 文献标志码:A
Application of interpretation method based on Gaussian Mixture Model cluster analysis to mixed oil and gas reservoirs of gas cap area in Chengbei Oilfield
ZHENG Zhuo①, WANG Jianli②, XIA Liangbing①, LIU Guangming①, ZHANG Hong①
①CNOOC Ener Tec-Drilling & Production Co., Tianjin 300452, China
②CNOOC China Limited Tianjin Company, Tianjin 300452, China
Abstract

The hydrocarbon accumulation structure of Chengbei low uplift was formed by the regional uplift of Yanshan and Himalayan epoch, and the regional gas reservoirs developed in the high position of the upper segment of the second Member for the Dongying Formation in Chengbei Oilfield. After more than 20 years of production, the gas cap area is not the original gas reservoir state, gas pool has the features of smaller scale and uneven distribution, and the conventional mud logging and well logging methods are difficult to judge the gas-oil interface. In the process of comprehensive adjustment of the oil field, in order to ensure that the horizontal well lands in the oil reservoirs and avoid the impact of mixed oil and gas reservoirs on the productivity, the gas logging data of 12 wells in Chengbei Oilfield were collected, the Gaussian Mixture Model was used for data association analysis, and the key separation coefficient ( α) was obtained.Combining the separation coefficient with the empirical formula, an interpretation method suitable for mixed oil and gas reservoirs of gas cap area in Chengbei Oilfield was established, which can be used to discriminate and interpret mixed oil and gas reservoirs in real time during the drilling process, so as to guide the accurate landing of horizontal wells.

Keyword: Gaussian Mixture Model; cluster analysis; mixed oil and gas reservoirs; separation coefficient; Chengbei Oilfield
0 引言

埕北油田位于渤海西部海域, 立足于一套井网, 一套开发层系, 依靠边水能量开发[1], 目前, 采出程度已达43%。长期高强度水驱生产方式影响储层流体性质的稳定性, 使储层气、油、水流体特征差异性降低, 含水率高达90%以上[2]。

为提高埕北油田整体产能, 大规模采用调整作业的模式优化井网, 经过多年实施“ 油气同采” 的开发模式, 致使埕北油田开发后期的气顶区储层压力降低、部分下部流体运移至上部气层中, 以油气混合的状态存在, 最终导致气顶区气层经过多年生产逐渐转变为混合油气层。目前气顶区域在调整井着陆作业过程中无法准确辨识混合油气层, 着陆于储层顶部混合油气层将影响油田调整井网部署, 同时也将加大钻井作业安全风险, 数据挖掘技术能从海量的数据中提取出有效信息, 发挥出每一项信息的长处, 并组合形成综合优势参数[3]。基于高斯混合模型创新的混合油气层解释方法, 对埕北油田调整作业有实质性指导意义。

1 埕北油田区域概况
1.1 地质概况

埕北油田构造位于济阳坳陷与渤中坳陷之间, 油田气顶油藏是渤海黄河口凹陷典型的一类油气藏, 顶部气层圈闭范围及规模由北部的北西向与南部的东西向断层决定。埕北油田东二上亚段油组为北东走向的断背斜, 受古地貌和后期构造作用的影响, 油田的主体部位拱升, 形成较典型的披覆背斜构造, 是油气富集的有利区带[4]。东二上亚段砂体夹于大套暗色湖相泥岩之间, 具有良好的盖层, 断裂两盘岩性以砂泥对接为主, 侧向封堵性好, 有优越的成藏条件[5]。开发储层以中细砂岩为主, 岩石类型主要为岩屑长石砂岩, 泥质胶结、疏松, 泥质夹层分布不稳定。东二上亚段油组储层厚度17~35 m, 储层平均孔隙度为28.6%, 平均有效渗透率为1 817.2 mD, 具有高孔、特高渗特征, 砂体叠置连片分布, 横向展布较为稳定。

埕北油田气顶区油气藏边、底水能量弱, 在早期“ 避气开发” 的开发模式下, 气顶区在开发后期发生气窜、油侵等现象, 由于气顶区不规则分布的油藏特性, 导致气顶指数差异大[6]。埕北油田气顶区油藏埋深为1 110.0~1 990.0 m, 主力单元油层厚度8.0~12.0 m, 气层厚度2.0~6.0 m, 目前油气界面混乱。

1.2 气顶区录测井特征

东二上亚段油组气藏是位于构造高部位的气顶, 气油界面为海拔-1 635 m左右, 气顶区探明含气面积约2.6 km2, 圈闭分布范围较广, 气藏储层延展性好。埕北油田在勘探开发前期, 为保证储层高开采率, 采用岩屑录井、气测录井、随钻伽马测井、随钻电阻率测井、随钻中子密度测井等较齐备的录测井方法; 但是在油田后期开发调整作业中, 为降低油田开发作业成本, 优选录测井作业方案, 仅采用岩屑录井、气测录井、随钻伽马测井、随钻电阻率测井这4项常规录测井方法。

通过对油田气顶区开发初期的录测资料进行分析发现, 气顶区气层在实钻中的气全量较油层高出10%~20%, 气层中子、密度反向交会特征明显, 电阻率呈突刺状, 油层电阻率曲线特征较平缓, 中子、密度交会不明显, 符合常规气层响应特征[7]。如图1所示, CB-5井为开发初期正常气层响应特征, CB-22井为开发后期混合油气层响应特征。在油田开发初期, 气顶区录测井响应特征明显, 常规油气界面识别方法适用于该油田的初期开发阶段。

图1 开发初期(CB-5井)与后期(CB-22井)气顶区地质完井综合图

通过对开发后期作业的12口井的资料统计分析发现, 气顶区的混合油气层与油层在自然伽马、电阻率、气全量等录测特征相似性较高, 着陆储层电阻率为10~20 Ω · m, 自然伽马为60~80 API, 气全量为10%~20%。目前使用的常规录测井解释方法符合率低于65%, 已无法满足实时解释混合油气层的要求, 为此很有必要建立适用于油田现状的混合油气层解释方法。

2 混合油气层解释方法
2.1 解释方法介绍

聚类算法是一种统计分类的方法, 目的是找出所研究参数之间存在的不同程度的相关性, 根据参数集找出各参数间的关联性, 以不同组分气测值为参数原本, 找到相同层位、不同深度条件下气测组分值之间的关联性, 使各参数之间的相似度尽可能小[8]。高斯混合模型依据聚类算法能有效地将无关联性的数据通过非监督训练迭代计算, 建立标准化数据关联项, 最终得到混合油气层组分参数之间的关联系数, 进而可将关联系数应用于后续作业井中。

埕北油田是经过二十余年生产开发的老油田, 气顶区开发调整井数量较多且资料齐全, 优选着陆层位于气顶区的12口井气测数据作为数据因子(包括全烃Tg、C1、C2、C3、iC4、nC4、iC5、nC5共8项气测参数), 将其作为基准项参与混合油气层解释方法的建立。

2.2 图板及解释标准建立

对埕北油田气顶区已钻的12口井进行分析与总结, 发现C1/C4在钻遇混合油气层时, 特征呈现极高值后突降(图2), 具有一定的指导意义。根据油层重组分高, 气层重组分低的固定特征, 采取重组分与轻组分比值法, 给定各组分值固定系数, 可得到归一化范围区间的无量纲值; 根据该无量纲值区间范围划分油气水界限, 且埕北油田气顶区混合油气层的特殊性, 不适用于气油比大于890的产气层工业气流标准[9]; 根据埕北油田已投产的12口井的气油比总结经验, 投产现状认为气油比高于100即可解释为混合油气层。以此建立了适用于埕北油田混合油气层含油指数(Io)公式及含气指数(Ig)经验公式:

Io= 10000×(C4+C5)2C1+C2+C3

Ig=lg C1+C2+C3(C4+C5)2

混合油气层解释方法使用高斯混合模型进行聚类分析, 将收集的气测数据库采用高斯混合模型作为调度函数, 首先使用极大似然估计算法迭代更新估计高斯混合模型中参数[10], 得到各项气测参数的标准值, 通过确定分离系数α 使气测参数能各向分离形成簇, 最终计算得出的分离系数α (表1)即为各项气测参数的计算系数值。

图2 埕北油田CB-32井混合油气层气体综合分析图

表1 埕北油田混合油气层分离系数α 值

将高斯混合模型聚类分析得出的分离系数α 值导入含油、气指数经验公式中, 得出以下气顶区混合油气层优化后含油指数(Io')、含气指数(Ig')解释公式:

Io'= 10000×(0.9C4+0.8C5)22.3C1+1.4C2+1.6C3

Ig'=lg 2.3C1+1.4C2+1.6C3(0.9C4+0.8C5)2

以优化后含气指数为横坐标, 优化后含油指数为纵坐标, 将气顶区气层解释公式计算得出的含油、气指数投入Excel建立含油、气指数解释图板, 形成直观的含油、气指数交会图板(图3)。根据埕北油田混合油气层解释图板判定:含油指数> 2.0, 且含气指数< 1.0解释为油层; 含油指数在1.5~3.0之间, 且含气指数在0.5~2.0之间解释为油气混合层; 含油指数< 1.5, 且含气指数> 1.5解释为气层。

图3 混合油气层解释图板

3 应用效果

将混合油气层解释图板应用于埕北油田气顶区近期3口正钻井中, 后经投产验证, 混合油气层解释图板在水平段储层解释中符合率较高, 现将这3口井的分析数据进行展示, 计算得出的结论表明混合油气层的含气指数驼峰明显, 与实际生产情况(表2)吻合。

表2 埕北油田近期作业3口井解释结论

埕北油田CB-18井是其中一口水平井, 储层着陆靶点位置处于油气混合层发育区块, 储层岩性为细砂岩, 受油气混合层影响, 实钻过程中的气测参数及测井参数无法实时解释着陆储层流体性质, 给着陆作业带来极大的挑战。在该井储层着陆过程中使用混合油气层解释图板解释得出(表3):在1 933.2~1 975.9 m井段含油指数介于2.2~2.5之间, 含气指数介于1.3~1.6之间, 图板解释为混合油气层; 在1 975.9~2 103.7 m井段含油指数介于2.7~3.5之间, 含气指数介于0.5~0.8之间, 图板解释为油层。通过钻后测井曲线解释判定水平段(1 980.7~2 103.7 m)流体性质为气层, 将混合油气层解释图板与测井解释结论对比分析(图4), 油气混合层与油层界面在1 988.5 m, 油气解释图板结论与投产结论基本一致。

表3 埕北油田CB-18井储层流体性质解释结论

图4 CB-18井混合油气层图板解释与测井解释对比

埕北油田CB-18井水平段(1 980.7~2 103.7 m)下入气密筛管后进行投产作业, 初期产能显示产油量为30.30 m3/d, 产气量为382 m3/d, 且气油比仅为18.19。通过产能判定本井水平段为油层, 所产出的气为混合气。生产结果与测井解释为气层的结论不符, 印证了本文混合油气层解释方法的准确性。

4 结论

本文通过深度挖掘气测数据在复杂混合油气层的应用效果, 依托真实的录井、测井、投产等多项参数, 将高斯混合模型的分离系数α 值与录井解释经验公式相结合, 形成适用于埕北油田开发后期的气顶区混合油气层的特定解释方法。

该解释方法具有一定局限性, 文中计算所得的分离系数仅适用于埕北油田开发后期的气顶区混合油气层。对于其他具有类似气层的油田, 尚需通过实际参数对分离系数α 值进行重新标定, 并需重新拟合经验公式后使用。

该解释图板在渤海油田广泛应用前需开展一系列测录井数据标定工作, 以提高录井定性解释油田开发后期的混合油气层准确性, 提升录井参数智能化解释应用效果。

(编辑 卜丽媛)

参考文献
[1] 苏彦春. 渤海湾盆地埕北油田高采收率主控地质因素剖析[J]. 中国海洋大学学报(自然科学版), 2016, 46(8): 87-95.
SUN Yanchun. Analysis on the main controlling geological factors about high recovery in Chengbei Oilfield, Bohai Gulf Basin[J]. Periodical of Ocean University of China, 2016, 46(8): 87-95. [本文引用:1]
[2] 李伟. 海上“特双高”油田剩余油挖潜关键技术研究与实践: 以埕北油田调整挖潜为例[J]. 中国石油和化工标准与质量, 2019, 39(19): 229-230, 232.
LI Wei. Research and practice on key techniques of remaining oil potential exploitation in offshore "super double high" oilfields: Taking Chengbei Oilfield as an example[J]. China Petroleum and Chemical Stand ard and Quality, 2019, 39(19): 229-230, 232. [本文引用:1]
[3] 郭素杰, 姜维寨, 孟庆峰, 等. 数据挖掘技术在复杂油气层评价中的应用探索[J]. 录井工程, 2019, 30(4): 1-7.
GUO Sujie, JIANG Weizhai, MENG Qingfeng, et al. Application exploration of data mining technology in complex oil and gas reservoir evaluation[J]. Mud Logging Engineering, 2019, 30(4): 1-7. [本文引用:1]
[4] 杨文杰. 歧口凹陷埕北断阶带油气运聚成藏机制研究[D]. 青岛: 中国石油大学(华东), 2010.
YANG Wenjie. Research on mechanism of the hydrocarbon migration and accumulation in Chengbei step-fault zone, Qikou Depression[D]. Qingdao: China University of Petroleum (EastChina), 2010. [本文引用:1]
[5] 马魁勇. 埕岛油田东营组油气成藏规律分析[J]. 内蒙古石油化工, 2012, 38(10): 115-118.
MA Kuiyong. Study on hydrocarbon accumulation of Dongying Formation, Chengdao Oilfield[J]. Inner Mongolia Petrochemical Industry, 2012, 38(10): 115-118. [本文引用:1]
[6] 苏彦春. 埕北油田次生底水形成机制及边水推进特征[J]. 石油地质与工程, 2015, 29(2): 90-92.
SU Yanchun. Secondary bottom water formation mechanism and the characteristics of edge water incursion in Chengbei Oilfield[J]. Petroleum Geology and Engineering, 2015, 29(2): 90-92. [本文引用:1]
[7] 孙玉娇. 油气水层解释方法研究[J]. 信息系统工程, 2019(8): 110.
SUN Yujiao. Study on interpretation methods of oil, gas and water layers[J]. China CIO News, 2019(8): 110. [本文引用:1]
[8] 尹楠. 基于高斯混合模型的期望最大化聚类算法[J]. 统计与决策, 2017(4): 87-89.
YIN Nan. The expectation maximization clustering algorithm based on Gaussian mixture model[J]. Statistics & Decision, 2017(4): 87-89. [本文引用:1]
[9] 全国国土资源标准技术委员会. 海上石油天然气储量计算规范: DZ/T 0252-2013[S]. 北京: 中国标准出版社, 2013.
National Technical Committee on Land and Resources Stand ards. Regulation of offshore petroleum reserves estimation: DZ/T 0252-2013[S]. Beijing: Stand ards Press of China, 2013. [本文引用:1]
[10] 王宏伟, 柴秀俊. 基于高斯混合模型聚类的非均匀采样系统的多模型切换辨识[J/OL]. 控制与决策: 1-9[2021-05-26].
WANG Hongwei, CHAI Xiujun. Multi-model switching identification for non-uniformly sampled systems based on Gaussian mixture model clustering[J/OL]. Control and Decision: 1-9[2021-05-26]. [本文引用:1]