一种机器学习辅助的两相流油藏随机建模的网格粗化方法

技术领域

本发明属于油藏数值模拟领域，具体涉及一种基于机器学习的油藏数值模拟网格粗化方法。

背景技术

建立储层地质模型，是油藏描述的核心与关键。储层建模方法分为确定性建模和随机建模。确定性建模是以井下测得的确定性资料为基础，推测井间确定的、唯一的储层参数。然而地下储层是许多复杂地质活动综合作用的结果，以非常有限的资料无法描述储层井间预测的不确定性。储层随机建模是地质统计学在石油工程中的广泛应用中发展起来的，是利用计算机科学，地质学、数学等多学科领域的技术再现储层的离散型、连续型和非均质性的方法。随机建模是以已知的信息为基础，以随机函数为理论，应用随机模拟方法产生多个可选的、等概率的储层模型(又称模型实现，realizations)的方法。虽然预测结果可能有任意多个，但其结果可以真实地反应储层属性的非均质性，多种结果也表现出了储层属性的不确定性，其中有好的结果，也有坏的结果，可以提供不同的参考。

随机建模技术产生了大量等概率的地质模型实现，每个模型实现都要进行油藏数值模拟，这使得模拟工作量增加数十倍乃至更多。网格粗化是降低模拟工作量的有效方法。网格粗化就是在尽可能保留细尺度模型的储层参数情况下，把细网格系统的数据信息转换到粗网格系统的过程。此方法可有效降低网格数量，大大降低计算成本，提高计算速度。传统的网格粗化算法需要对每个粗尺度网格系统进行大量“重复性”计算，而机器学习方法正适合处理这种“重复性”工作。

发明内容

本发明的目的是针对上述不足，提出了一种在保证计算精度的前提下，能显著提高网格粗化速度的机器学习辅助的两相流油藏随机建模的网格粗化方法。

本发明具体采用如下技术方案：

一种机器学习辅助的两相流油藏随机建模的网格粗化方法，包括以下步骤：

(1)通过随机建模方法建立若干地质模型实现，确定模型实现将要粗化成粗网格的尺寸和数量；

(2)选取部分地质模型实现记为模型集P1，其余模型记为模型集P2，对模型集P1中模型实现进行相对渗透率粗化计算，得到P1中粗网格的粗尺度相对渗透率数据；

(3)对所有地质模型实现的渗透率进行数据预处理；

(4)提取P1和P2中每个粗网格的渗透率分布特征，用机器学习算法对所有粗网格的渗透率分布特征进行归类，得到P1和P2中粗网格的类别；

(5)将P1中粗网格的粗尺度相对渗透率数据赋值给相似度最高的P2中的网格，得到P2近似的粗尺度相对渗透率数据；

(6)使用P1和P2中粗网格粗尺度的相对渗透率数据进行油藏数值模拟计算。

优选地，对所有地质模型实现的渗透率进行数据预处理的处理方式为：先对所有细网格的渗透率做ln对数变换，然后对每个粗网格进行特征放缩处理，特征放缩处理式为(1)所示：

其中，x

优选地，步骤(2)中选取地质模型实现划分模型集P1，在所有模型实现中按比例随机选取。

优选地，步骤(2)中选取地质模型实现划分模型集P1，通过手动选取部分有明显地质特征的模型实现。

优选地，步骤(4)中提取粗网格的渗透率分布特征，直接采用数据预处理后粗网格中每个细网格的渗透率数据。

优选地，步骤(4)中提取粗网格的渗透率分布特征，通过把数据预处理后的粗网格中细网格的渗透率场画成图，以图片识别的方式对粗网格进行归类。

优选地，步骤(4)中提取粗网格的渗透率分布特征，对渗透率分布特征进行降维处理以加快计算机处理速度。

优选地，步骤(4)中用机器学习算法对所有粗网格的渗透率分布特征进行归类，归类方法为分类或聚类算法。

本发明具有如下有益效果：

该方法可以避免随机建模后对每个地质模型实现进行网格粗化计算，仅需对少数地质模型实现进行网格粗化计算，通过机器学习算法识别相似网格，对其余模型实现中相似网格快速进行结果赋值。

该方法可显著提高随机建模网格粗化速度，对提高随机建模模拟效率，加快油藏数值模拟智能化，具有很高的实用价值。

附图说明

图1为一种机器学习辅助的两相流油藏随机建模的网格粗化方法的流程框图；

图2为实施例二的流程图；

图3为实施例二中的地质模型之一的渗透率分布及其特征缩放的示例图；

图4为实施例二、实施例三和实施例四的计算时间对比图；

图5为实施例二、实施例三和实施例四产量不确定性分析结果对比图；

图6为实施例一和实施例四在不同时间下所有地质模型实现的总流量散点拟合直线对比。

具体实施方式

下面结合附图和具体实施例对本发明的具体实施方式做进一步说明：

实施例一

结合图1机器学习辅助的两相流油藏随机建模的网格粗化方法，包括以下步骤：

(1)通过随机建模方法建立若干地质模型实现，确定模型实现将要粗化成粗网格的尺寸和数量。

(2)选取部分地质模型实现记为模型集P1，其余模型记为模型集P2，对模型集P1中模型实现进行相对渗透率粗化计算，得到P1中粗网格的粗尺度相对渗透率数据。选取地质模型实现划分模型集P1，在所有模型实现中按比例随机选取，或手动选取部分有明显地质特征的模型实现。

(3)对所有地质模型实现的渗透率进行数据预处理，处理的处理方式为：先对所有细网格的渗透率做ln对数变换，然后对每个粗网格进行特征放缩处理，特征放缩处理式为(1)所示：

其中，x

(4)提取P1和P2中每个粗网格的渗透率分布特征，用机器学习算法对所有粗网格的渗透率分布特征进行归类，得到P1和P2中粗网格的类别。

提取粗网格的渗透率分布特征，可以直接采用数据预处理后粗网格中每个细网格的渗透率数据，也可把数据预处理后的粗网格中细网格的渗透率场画成图，以图片识别的方式对粗网格进行归类，达到更直观的可视化效果；

提取粗网格的渗透率分布特征，可以对渗透率分布特征进行降维处理以加快计算机处理速度，降维处理方法包括但不限于PCA、T-SNE、LDA、MDS等。

用机器学习算法对所有粗网格的渗透率分布特征进行归类，归类方法可以是分类或聚类算法，具体地，包括但不限于KNN、K-Means、SVM、随机森林、朴素贝叶斯等。

(5)将P1中粗网格的粗尺度相对渗透率数据赋值给与之同类(或相似度最高)的P2中的网格，得到P2近似的粗尺度相对渗透率数据。

(6)使用P1和P2中粗网格粗尺度的相对渗透率数据进行油藏数值模拟计算。

实施例二

基于上述方法，通过机器学习KNN算法，对数据预处理后的地质模型实现进行网格分类，整体方法流程图如图2，结合图3，具体步骤如下：

第一步：生成110个地质模型实现，每个实现由200×200个细网格组成，目标为将细尺度模型粗化成由尺寸为10×20的粗网格系统组成的粗尺度模型。考虑水驱油的两相渗流问题，设定边界条件为上下边界封闭，左侧为注水端，右侧为采油端。

第二步：随机选取10个模型实现(记为P1)进行相对渗透率粗化计算，得到P1中粗网格的粗尺度相对渗透率数据，其余100个模型实现记为P2。

第三步：为更好地提取每个粗网格的非均质性，以便于机器学习算法识别，对所有地质模型渗透率进行数据预处理，处理方法为：对所有细网格的渗透率值做ln对数变换，并对所有粗网格做特征放缩处理。数据处理后的网格示意图如图3所示，数据处理后粗网格的非均质性更加突出。

特征放缩处理公式为：

其中，x

第四步：以粗网格中细网格的渗透率数据为特征，用KNN算法对粗网格进行分类，得到P2和P1中粗网格的分类关系。

第五步：将P1中粗网格的粗尺度相对渗透率数据赋值给P2中与之同类的粗网格。

第六步：使用P1和P2中粗网格的粗尺度相对渗透率，对110个地质模型实现进行油藏数值模拟计算，得到产量等数据，记录上述步骤花费时间合计为7029秒。

实施例三：

对所有地质模型实现不进行网格粗化，直接以细尺度模型进行油藏数值模拟计算，得到产量等数据，记录数值模拟计算花费时间合计为445658秒。

实施例四：

对所有地质模型实现进行传统的两相流网格粗化计算(不以机器学习方法辅助)，然后进行油藏数值模拟计算，得到产量等数据，记录网格粗化数值模拟花费时间合计为66555秒。

图4为细尺度数值模拟、网格粗化数值模拟和机器学习辅助网格粗化的数值模拟计算时间对比，机器学习辅助网格粗化的数值模拟效率提高明显，效率相比细尺度数值模拟提高63.4倍，相比网格粗化数值模拟提高9.5倍。

图5为对随机建模总产量的不确定性分析结果对比。P90、P50和P10分别代表所有模型实现数据中有90％、50％和10％小于该数据，代表总产量概率的预估。PVI为无量纲时间，代表不同的注水阶段，计算公式为

图6为PVI＝0.5和PVI＝1，在不同模拟方法下，110个地质模型实现的总流量q

当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。