-
相渗曲线能够综合反映流体在多孔介质中的渗流特征，是用于油田开发参数计算、动态分析、油藏数值模拟研究中的重要资料[1]。在油田不断开发的过程中，要了解油藏物性和渗流规律动态变化，需要实时获取储层相对渗透率变化[2]。传统的储层相渗数据的获取主要有3种方法：室内实验法、经验公式法和测井资料解释法[3-9]。室内实验法能获得完整的相渗曲线，但过程比较精细、昂贵且耗时；经验公式法是在缺乏实验数据的情况下进行相关估计，精度低、误差大；测井资料解释法只是得出在某一含水饱和度下的相对渗透率值，无法反映油水的相对渗透率变化情况。近年来，随着大数据、人工智能技术的发展，石油行业已开始利用人工智能方法获取储层物性参数[10-14]。为了获取更高精度的渗透率，开始尝试采用机器学习算法对储层渗透率进行预测。但是从调研结果来看，目前人工智能仅限于某一个状态下的储层渗透率预测，未能实现动态渗透率预测，即储层相渗曲线的预测。
-
目前，胜利油田已有大量的水驱储层相渗预测样本资料，尤其是能够反映储层状态的测井曲线，每口井都有多条不同类型的测井曲线，这些资料符合大数据技术特点：重复性、差异性和可预测性。因此，基于油田大数据，优选人工智能算法，开展储层相渗曲线仿真成为可能。
-
笔者提出一种新的水驱储层相渗曲线仿真方法，将测井曲线值作为输入，相渗曲线值作为输出，经过测井曲线敏感参数的优选、训练集的构建、机器学习算法的优选、模型的优化等一系列流程，生成仿真相渗曲线，通过多属性条件约束优化模型，提高相渗曲线预测精度，从而填补了基于机器学习算法预测完整相渗曲线的空白。该方法应用于油藏工程和油藏数值模拟研究中，可有效提升地质研究精度和效率，具有良好的应用前景。
-
1 相渗曲线样本构建
-
1.1 数据准备
-
相渗曲线主要由端点值（即最大油相相渗值、最大水相相渗值、残余油饱和度和束缚水饱和度） 和等渗点等参数来表征，而影响相渗曲线形态的因素主要包括饱和历程、润湿性、岩石孔隙几何形态、大小、分布、原油黏度和温度等[15]。影响因素大部分是岩石流体物性参数，需要通过取样并进行岩心分析才能获得。胜利油田共 6 万多口井，其中做过岩心分析的井仅 4 000 多口，数据来源少。在储层预测研究中，应用测井资料能够间接反映储层的渗透性[16]，从而预测储层相渗曲线。研究收集的测井数据几乎每口井都有，占比为 100%。因此，将测井曲线值作为主要特征参数输入，相渗曲线数据做标签建立模型预测。
-
1.2 测井参数敏感性分析
-
从反映地层特性的全面性以及数据量方面考虑，选取自然伽马（GR）、自然电位（SP）、井径 （CAL）、八侧向电阻率（RFOC）、感应电导率 （COND）、4 m 底部梯度电阻率（R4）、微电位电阻率 （RNML）、微梯度电阻率（RLML）、2.5 m 底部梯度电阻率（R25）、密度（DEN）、中子（CNL）、声波时差 （AC）、总孔隙度（PORT）13个测井参数，分别建立各个测井参数与相渗曲线的相关性，在此基础上，排除高度相关的测井参数，最终分析认为 GR，CAL， SP，AC，R25，COND，DEN等7个测井参数能够全面反映岩性、电性、孔隙度三方面特征。
-
选取孔隙度、渗透率、束缚水饱和度、残余油饱和度、束缚水下油相相对渗透率、残余油下水相相对渗透率 6 个特征参数，利用 K-Means 聚类方法对收集到的水驱储层280多个相渗曲线样品进行编号并自动聚类，通过调参求出 K-Means聚类方法最佳的分类数为5，5类相渗曲线之间的曲线形态差异较大，从中选出具有代表性的样品，对应的样本编号分别为3，105，125，126和228。
-
利用上述5个相渗曲线样品对应的岩石样品顶深与底深，分别截取对应深度的测井曲线段并求平均值（表1），再对测井曲线进行相对变化率计算，得到对应的7个测井参数的相对变化率（表2）。
-
由计算得到的测井参数相对变化率绘制成雷达图（图1），图中5个相渗曲线类型分别用不同颜色代表，雷达图中各测井参数对应轴上的点代表5 个样品各测井参数的相对变化率，同一个样品点相连围成的面积越大，且同一对应轴上的值分布范围越大，表示该测井参数对相渗曲线越敏感。分析认为AC，SP，DEN 这 3 个测井参数不敏感；GR，CAL，R25， COND 这 4 个测井参数对于相渗曲线而言是敏感的，将其作为预测模型输入特征参数。
-
图1 基于测井曲线的储层相渗曲线影响因素敏感性分析
-
Fig.1 Sensitive analysis of factors affecting reservoir relative permeability curves based on logging curves
-
1.3 样本构建
-
选取 GR，CAL，R25，COND 和相渗实验中的储层含水饱和度（Sw）共 5个参数作为预测相渗曲线的特征输入，将油相相对渗透率（Kro）和水相相对渗透率 （Krw）作为预测相渗曲线的特征输出。
-
基于原始的样本数据，进行数据分析、数据量及分布情况统计等数据预处理；采用分箱法处理异常值；采取去重处理、噪音数据处理等，对缺失的测井曲线值采用 XGBoost 算法进行补全；将岩心相渗数据与岩心深度段相对应的测井参数融合构建相渗曲线预测样本集（表3）。
-
2 模型训练与优化
-
2.1 模型训练结果对比
-
采用机器学习算法对胜坨区块的相渗曲线进行预测，使用样本数量为 278 个。样本数据的 90% 作为训练数据，10% 作为测试数据；其中训练数据的 90% 用于模型的训练，10% 用于模型的交叉验证。分别选用随机森林，AdaBoost，GBDT和XGBoost 这4种机器学习算法对油相相对渗透率和水相相对渗透率进行预测，并对预测精度进行分析对比。
-
注：同一个测井位置对应多个岩心样品、多条相渗曲线的情况下，先对相渗曲线做归一化处理。
-
预测结果表明，4 种算法预测精度均在 85% 以上，其中XGBoost算法预测效果最好，油相相对渗透率预测精度达 90.4%，水相相对渗透率预测精度达 93.8%；其次为 GBDT 算法，油相相对渗透率预测精度为 90.0%，水相相对渗透率预测精度为 92.1%；随机森林与 AdaBoost 算法的油相相对渗透率预测精度分别为 86.0% 和 87.2%，水相相对渗透率预测精度分别为 90.2% 和 92.1%。XGBoost 算法优于其他算法的原因是：①应用正则化项防止过拟合。②不仅使用了一阶导数，还使用二阶导数，损失更精确。 ③可并行优化。④考虑了训练数据为稀疏值的情况，大大提升算法效率。⑤支持列抽样，不仅能降低过拟合，还能减少计算时间。
-
2.2 模型优化
-
2.2.1 地质因素约束优化
-
考虑油藏的地质特征，在原有预测模型的基础上添加、融合油藏地质特征参数进行模型约束，主要包括：①坐标数据作为沉积相类型的反映与约束，结合井斜数据，计算出每个井层的坐标数据，间接反映油藏沉积类型。②地质分层数据作为测井数据离散标准化的深度约束，为了消除由于非地质因素造成的误差，需要对测井资料进行标准化处理。③通过查阅岩心综合录井资料，将岩心样品进行归位。
-
水驱储层相渗曲线预测模型地质因素约束优化后，预测精度有所提升，其中XGBoost算法油相相对渗透率预测精度达 93.0%，水相相对渗透率预测精度达96.3%，优化后精度提升超过2%。
-
2.2.2 曲线端点约束优化
-
束缚水饱和度（Sw1）和残余油饱和度（Sw2）端点值决定了相渗曲线的区间范围（图2）。基于相渗曲线的预测样本中输入特征除了测井曲线值外，还包括含水饱和度数据。在模型的应用过程中，测井资料充足，但相渗实验中的含水饱和度资料缺乏。因此必须把含水饱和度的端点值预测出来，才能更好地应用模型。
-
图2 相渗曲线
-
Fig.2 Relative permeability curve
-
在原样本集的基础上，输入特征为GR，CAL，R25 和 COND，输出特征为含水饱和度，采用 XGBoost 算法的结果可多输出的特点同时预测束缚水饱和度和残余油饱和度。其中，束缚水饱和度的预测精度为93%，残余油饱和度的预测精度为94.9%。
-
在模型预测相渗曲线的应用过程中，结合饱和度端点值预测，以预测的2个端点值为边界，等间隔划分相渗曲线的含水饱和度，将含水饱和度和测井曲线值共同作为输入，经过模型预测后，即可得到不同含水饱和度下的油水相对渗透率，绘制完整的相渗曲线。
-
3 模型验证及应用
-
3.1 模型验证
-
利用验证集中的样本进行相渗曲线预测，并与其真实值做比较，由相渗曲线模型的预测结果（图3）可以看出，油相相对渗透率与水相相对渗透率的预测值与真实值曲线几乎重合，计算得到油相相对渗透率与水相相对渗透率的绝对误差均小于 10%。说明该模型预测效果较好。
-
图3 相渗曲线预测值与真实值对比
-
Fig.3 Comparison between predicted values and true values of relative permeability curves
-
3.2 实例应用
-
为便于水驱储层预测模型的应用，研发可视化展示功能，通过确定区块及预测层段，自动提取测井资料，即可实现油水相渗曲线的动态生成。基于本文得到的预测模型，生成 1 543 口井的相渗曲线数据，应用于孤东七区西馆上段 63+4、孤东油田七区西馆上段 52+3、孤岛油田中一区馆陶组 3-6、孤岛油田西区馆陶组3-4、坨21块沙二段9等5个水驱储层区块单元的数值模型建立。相渗数据主要用于模型中剩余油分布和高耗水条带判识方面，采用不同层段相渗曲线赋值代替使用平均值，充分反映层间及层内的非均质性，提升建模准确度。以孤东油田七区西馆上段 52+3 含油饱和度分布为例，常规做法是取该单元或相似单元的残余油饱和度，取值为 0.21，但利用模型预测结果，获得 129个不同井层相渗曲线，残余油饱和度为 0.15～0.37，充分体现油藏模型的非均质性，该区块建模精度提升0.056。该成果实现油藏不同空间位置相渗曲线的赋值，广泛应用于区块研究中，大大提升油藏区块或单元相渗曲线的获取效率。
-
4 结论
-
基于机器学习算法的水驱储层相渗曲线仿真预测优势在于能够实现每口井每个层段的相渗曲线预测，利用储层相渗预测模型，预测精度大于90 %，能够达到矿场应用精度。该方法不受实验条件、油田开发阶段等限制，充分利用油田测井资料，实时生成井层相渗曲线，准确反映油藏渗流特征和储层渗透率变化规律，为油藏工程研究中油田、区块、井组等不同层级提供相渗数据，并解决油藏工程数值模拟时目标块都用同一条相渗曲线的问题，使数值模拟更准确。
-
参考文献
-
[1] 秦同洛，李璗，陈元千.实用油藏工程方法［M］.北京：石油工业出版社，1989.QIN Tongluo，LI Dang，CHEN Yuanqian.Practical methods of pe⁃ troleum reservoir engineering［M］.Beijing：Petroleum Industry Press，1989.
-
[2] 姜瑞忠，乔欣，滕文超，等.储层物性时变对油藏水驱开发的影响［J］.断块油气藏，2016，23（6）：768-771.JIANG Ruizhong，QIAO Xin，TENG Wenchao，et al.Impact of physical properties time variation on waterflooding reservoir de⁃ velopment［J］.Fault-Block Oil and Gas Field，2016，23（6）：768-771.
-
[3] 王东琪，殷代印.水驱油藏相对渗透率曲线经验公式研究［J］.岩性油气藏，2017，29（3）：159-164.WANG Dongqi，YIN Daiyin.Empirical formulas of relative perme⁃ ability curve of water drive reservoirs［J］.Lithologic Reservoirs，2017，29（3）：159-164.
-
[4] 谢伟伟，王晓冬，李晨，等.基于流线应用特征线方法求解双重孔隙介质两维两相渗流问题［J］.油气地质与采收率，2018，25（1）：94-99.XIE Weiwei，WANG Xiaodong，LI Chen，et al.Application of char⁃ acteristics method for two-dimensional，two-phase seepage prob⁃ lem in dual-porosity medium based on streamline［J］.Petroleum Geology and Recovery Efficiency，2018，25（1）：94-99.
-
[5] 吕伟峰，秦积舜，吴康云，等.低渗岩石孔渗及相对渗透率测试方法综述［J］.特种油气藏，2011，18（3）：1-6.LÜ Weifeng，QIN Jishun，WU Kangyun，et al.An overview of test⁃ ing methods of poroperm and relative permeability for low permea⁃ bility rocks［J］.Special Oil & Gas Reservoirs，2011，18（3）：1-6.
-
[6] 王继强，岳圣杰，朱孟高，等.计算油田特高含水期相渗曲线新方法［J］.石油化工应用，2017，36（8）：62-64，88.WANG Jiqiang，YUE Shengjie，ZHU Menggao，et al.A new meth⁃ od of calculating relative permeability curve in ultra-high watercut oilfield［J］.Petrochemical Industry Application，2017，36（8）：62-64，88.
-
[7] 崔传智，郑文乾，李立峰，等.基于动态数据反演的相渗曲线及应用效果［J］.石油钻采工艺，2019，41（4）：516-520.CUI Chuanzhi，ZHENG Wenqian，LI Lifeng，et al.The relative permeability curve based on dynamic data inversion and its appli⁃ cation effect［J］.Oil Drilling & Production Technology，2019，41（4）：516-520.
-
[8] 杨宇，周文，邱坤泰，等.计算相对渗透率曲线的新方法［J］.油气地质与采收率，2010，17（2）：105-107，118.YANG Yu，ZHOU Wen，QIU Kuntai，et al.A new method of calcu⁃ lating relative permeability curve［J］.Petroleum Geology and Re⁃ covery Efficiency，2010，17（2）：105-107，118.
-
[9] 马东，吴华，曾鸣.从电阻率数据中得到相对渗透率的新方法 ［J］.石油与天然气地质，2015，36（4）：695-700.MA Dong，WU Hua，ZENG Ming.A new method for estimating rel⁃ ative permeability with resistivity data［J］.Oil & Gas Geology，2015，36（4）：695-700.
-
[10] 董满仓，柴汝宽，辛晶，等.致密储层动态渗透率预测模型［J］.陕西科技大学学报，2018，36（1）：90-95.DONG Mancang，CHAI Rukuan，XIN Jing，et al.Study on predic⁃ tion of dynamic permeability based on improved BP neural net⁃ work［J］.Journal of Shaanxi University of Science & Technology，2018，36（1）：90-95.
-
[11] 张雁，柳成志，秦秋寒，等.利用人工神经网络预测砂岩储层渗透率［J］.大庆石油学院学报，2005，29（4）：10-11，32.ZHANG Yan，LIU Chengzhi，QIN Qiuhan，et al.Predicting reser⁃ voir permeability of sandstone by means of artificial neural net⁃ work［J］.Journal of Daqing Petroleum Institute，2005，29（4）：10-11，32.
-
[12] 王猛，刘志杰，杨玉卿，等.基于区域测井大数据和实验资料的储层流动单元渗透率建模方法［J］.地球物理学进展，2021，36（1）：274-280.WANG Meng，LIU Zhijie，YANG Yuqing，et al.Permeability cal⁃ culation in reservoir flow unit based on regional logging big data and experimental data［J］.Process in Geophysics，2021，36（1）：274-280.
-
[13] 夏宏泉，张贤辉，范翔宇，等.基于神经网络法的逐点渗透率测井解释研究［J］.西南石油学院学报，2001，23（1）：11-13.XIA Hongquan，ZHANG Xianhui，FAN Xiangyu，et al.The study on predicting each point’s permeability based on neural network and log data［J］.Journal of Southwest Petroleum Institute，2001，23（1）：11-13.
-
[14] 杨建，杨程博，张岩，等.基于改进神经网络的渗透率预测方法 ［J］.岩性油气藏，2011，23（1）：98-102.YANG Jian，YANG Chengbo，ZHANG Yan，et al.Permeability prediction method based on improved BP neural network［J］.Lith⁃ ologic Reservoirs，2011，23（1）：98-102.
-
[15] 孙雷，闫成海，潘毅，等.长岩芯驱替中不同位置的相渗曲线计算［J］.西安石油大学学报：自然科学版，2014，36（2）：139-144.SUN Lei，YAN Chenghai，PAN Yi，et al.Calculation of relative permeability curve at different position in long core displacement experiment［J］.Journal of Southwest Petroleum University：Sci⁃ ence & Technology Edition，2014，36（2）：139-144.
-
[16] 王庚阳，钟祖兰.利用测井资料确定储层油水相渗透率方法 ［J］.石油学报，1991，12（4）：75-80.WANG Gengyang，ZHONG Zulan.A method for determining the oil/water phase permeability using well logging data［J］.Acta Petrolei Sinica，1991，12（4）：75-80.
-
摘要
相渗曲线是油气田开发研究中的一项重要基础资料。采用常规室内实验法获得相渗曲线费用昂贵且耗时，测试样品少，难以代表整个油藏的特征；经验公式法估算获得的结果精度低且误差大。为了实时、准确获得水驱储层相渗曲线，采用机器学习算法进行仿真预测。通过测井参数敏感性分析，融合相渗曲线数据，构建水驱储层相渗曲线仿真样本集。在此基础上，优选机器学习算法进行地质因素约束优化以及曲线端点约束优化，实现相渗曲线智能可视化生成。研究结果表明：该方法能实现每口井每个层段的相渗曲线预测，预测精度大于90%，能准确反映油藏渗流特征和储层渗透率变化规律，具有较高的实际应用价值和良好的推广应用前景。
Abstract
The relative permeability curve is an important basic datum in oil and gas field development research. The tradi- tional experimental method is expensive and time-consuming，and the relative permeability curve obtained by a few test samples is difficult to represent the characteristics of the whole reservoir. The results from the empirical formula method have low precision and large error. To obtain the relative permeability curves of the water-flooding reservoir in real time and accurately，this paper utilizes machine learning algorithms for simulation and prediction. The simulation sample set of the water-flooding reservoir relative permeability curves is constructed by analyzing the sensitivity of logging parameters and integrating the data of the relative permeability curves. On this basis，the machine learning algorithm is selected to opti- mize the geological factor constraint and the curve endpoint constraint，and the intelligent visualization generation of the rel- ative permeability curves is realized. The results show that this method can realize the prediction of the relative permeabili- ty curves for each well and each section，and the prediction accuracy is more than 90%. It can accurately reflect the flow characteristics of reservoir and the variation law of reservoir permeability，possessing high practical application value and a good prospect of popularization and application.
