基于LSTM的致密低渗油藏压裂井网气驱效果预测方法与流程

基于lstm的致密低渗油藏压裂井网气驱效果预测方法
技术领域
1.本发明是关于一种基于lstm的致密低渗油藏压裂井网气驱效果预测方法，涉及油藏工程和深度学习领域。


背景技术：

2.压裂气驱是当前开采致密低渗油藏最有效提高采收率的技术之一。准确模拟预测压裂气驱过程是实现致密低渗油藏注气高效开发的关键理论基础。目前进行生产效果预测的方法主要有传统油藏工程方法和油藏数值模拟方法。
3.传统油藏工程方法仅能建立较简单的预测模型，对油藏类型、开发阶段、工作制度等差异难以考虑，导致预测精度有限。油藏数值模拟方法可以模拟实际生产井条件，预测任意措施类型的效果，但是需要以准确的地质模型和历史拟合为基础，模拟拟合预测周期较长，调参繁琐，不适合短期、局部的快速调整预测。


技术实现要素：

4.针对上述问题，本发明的目的是提供一种能够对油藏生产进行准确且高效预测，准确复制油藏模拟器的行为，进而提高生产效率的基于lstm的致密低渗油藏压裂井网气驱效果预测方法、系统、电子设备及介质。
5.为实现上述目的，本发明采取以下技术方案：
6.第一方面，本发明提供一种基于lstm的致密低渗油藏压裂井网气驱效果预测方法，包括：
7.建立致密低渗油藏压裂井网注气开发双重介质数值模拟基础模型；
8.基于数值模拟基础模型，建立多套致密低渗油藏注气开发双重介质数值模拟模型，形成致密低渗油藏注气开发数据库；
9.筛选用于构建日产油、地层压力、采收率代理模型相关的特征参数；
10.将筛选出的特征参数使用基于均方根反向传播的lstm算法，分别建立并训练日产油、地层压力及采收率的预测代理模型。
11.进一步地，基于数值模拟基础模型，建立多套致密低渗油藏注气开发双重介质数值模拟模型，形成致密低渗油藏注气开发数据库包括：
12.设计正交试验，用于筛选有代表性试验方案；
13.基于数值模拟基础模型data文件和正交试验方案，使用python编程调用eclipse软件生成多套致密低渗油藏注气开发双重介质数值模拟模型data文件，分别模拟不同地质\流体\油藏生产过程；
14.使用python自动提取eclipse软件运行多套注气开发数值模拟模型的结果数据，生成规范统一数据格式，包括日产油、地层压力、采收率。
15.进一步地，设计正交试验使用allpairs软件。
16.进一步地，筛选用于构建日产油、地层压力、采收率代理模型相关的特征参数基于
随机森林算法。
17.第二方面，本发明提供的一种基于lstm的致密低渗油藏压裂井网气驱效果预测系统，该系统包括：
18.基础模型建立单元，被配置为建立致密低渗油藏压裂井网注气开发双重介质数值模拟基础模型；
19.数据库建立单元，被配置为基于数值模拟基础模型，建立多套致密低渗油藏注气开发双重介质数值模拟模型，形成致密低渗油藏注气开发数据库；
20.特征参数筛选单元，被配置为筛选用于构建日产油、地层压力、采收率代理模型相关的特征参数；
21.预测模型建立单元，被配置为将筛选出的特征参数使用基于均方根反向传播的lstm算法，分别建立并训练日产油、地层压力及采收率的预测代理模型。
22.第三方面，本发明还提供的一种电子设备，所述电子设备至少包括处理器和存储器，所述存储器上存储有计算机程序，所述处理器运行所述计算机程序时执行所述的方法。
23.第四方面，本发明还提供一种计算机存储介质，其上存储有计算机可读指令，所述计算机可读指令可被处理器执行以实现所述的方法。
24.本发明由于采取以上技术方案，其具有以下优点：
25.1、本发明将深度神经网络算法与油藏压裂注气开采技术相结合，实现基于数据模型准确且高效预测致密低渗油藏压裂井网气驱生产效果，为油藏模拟提供了新途径，在解决复杂的油藏工程和提高采收率模拟预测问题时，有效兼顾易用性和精确性；
26.2、本发明使用代理模型将深度学习与油藏工程相结合，提高生产效率，使用代理模型实现基于数据模型对油藏生产准确且高效进行预测，准确地复制油藏模拟器的行为，提高生产效率，对致密低渗油藏开发与管理有重大意义；
27.3、本发明基于长短期记忆单元(lstm)的模型是对递归神经网络(rnn)的扩展，它能够清晰解决梯度消失问题，成为建模时间序列的流行选择，lstm模型实际上扩展了递归神经网络内存，使它们能够保存和学习输入的长期依赖关系；
28.4、本发明应用随机森林算法筛选包含地质与工程因素的特征，简单、计算成本小，防止过拟合，油田生产由一系列相互作用的复杂参数控制，包括地质参数、与流体压力/体积/温度及其性质相关参数、与油田开发决策相关参数，所有参数重要性可用数字来描述和估计，用于确定包括日产油、地层压力、采收率的油藏整体性能，为防止样本数据量小产生过拟合风险，本发明首先应用随机森林算法研究不同输入变量的影响，进而确定参与预测的特征参数；
29.5、本发明自动化建立、运行数值模拟模型，提取数据，基于eclipse数值模拟软件，使用python编程，实现自动生成、运行多套致密低渗油藏注气开发数值模拟模型，并自动提取结果数据；
30.6、本发明使用均方根传播算法，抑制梯度的锯齿下降，均方根传播(rmsprop)是geoff hinton提出的一种自适应学习率方法，rmsprop将动量思想引入梯度的累加计算中，使用指数衰减平均以丢弃遥远过去历史，使其能够在找到凸碗状结构后快速收敛，抑制梯度的锯齿下降；
31.综上，本发明为油藏模拟提供了新途径，在解决复杂的油藏工程和提高采收率模
拟预测问题时，可有效兼顾易用性和精确性。
附图说明
32.通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。在整个附图中，用相同的附图标记表示相同的部件。在附图中：
33.图1是本发明实施例的基于lstm的致密低渗油藏压裂井网气驱效果预测方法流程图；
34.图2是本发明实施例的致密低渗压裂井网注气开发数值模拟模型示意图；
35.图3是本发明实施例lstm算法预测代理模型真实值预测值对比图，(a)lstm算法预测日产油代理模型真实值预测值对比；(b)lstm算法预测地层压力代理模型真实值预测值对比；(c)lstm算法预测采收率代理模型真实值预测值对比；
36.图4是本发明实施例的电子设备结构图。
具体实施方式
37.下面将参照附图更详细地描述本发明的示例性实施方式。虽然附图中显示了本发明的示例性实施方式，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了能够更透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
38.应理解的是，文中使用的术语仅出于描述特定示例实施方式的目的，而无意于进行限制。除非上下文另外明确地指出，否则如文中使用的单数形式“一”、“一个”以及“所述”也可以表示包括复数形式。术语“包括”、“包含”、“含有”以及“具有”是包含性的，并且因此指明所陈述的特征、步骤、操作、元件和/或部件的存在，但并不排除存在或者添加一个或多个其它特征、步骤、操作、元件、部件、和/或它们的组合。文中描述的方法步骤、过程、以及操作不解释为必须要求它们以所描述或说明的特定顺序执行，除非明确指出执行顺序。还应当理解，可以使用另外或者替代的步骤。
39.为了便于描述，可以在文中使用空间相对关系术语来描述如图中示出的一个元件或者特征相对于另一元件或者特征的关系，这些相对关系术语例如为“内部”、“外部”、“内侧”、“外侧”、“下面”、“上面”等。这种空间相对关系术语意于包括除图中描绘的方位之外的在使用或者操作中装置的不同方位。
40.近年来，深度学习(dl)算法的快速发展为解决工程领域复杂多变量问题提供了有效途径。递归神经网络(rnn)已经成为建模时间序列的流行选择，基于长短期记忆单元(lstm)的模型是对递归神经网络的扩展，能够清晰解决梯度消失问题。本发明将深度神经网络算法与油藏压裂注气开采技术相结合，实现基于数据模型对油藏生产准确且高效预测，有效提高了生产效率。
41.实施例一
42.如图1所示，本实施例提供的基于lstm的致密低渗油藏压裂井网气驱效果预测方法，包括：
43.s1、为实现使用python编程调用eclipse软件自动化生成、运行多套数值模拟方
案，本实施例使用eclispe软件建立致密低渗油藏压裂直井/水平井混合井网注气开发双重介质数值模拟基础模型。
44.本实施例的数值模拟基础模型用于对致密低渗油藏直井\水平井混合压裂井网气驱生产进行模拟，基础参数包括网格大小、网格数量、介质属性、网格总数、顶深、油藏中心生产井和注气井参数等。
45.具体地，如图2所示，本实施例建立的数值模拟基础模型的基本参数为：网格大小为10m*10m*10m，网格数量为41*76*8，笛卡尔坐标系下油藏z方向上4层为基质，下4层为天然裂缝，网格总数为24928个，顶深为1815m。本实施例的油藏中心一口水平生产井，位于实际油藏的第3层，长度为500m，井径为0.1m，两端有射孔，中间段射孔位置、数量不等，模拟生产20年。采用eclipse软件的e300模拟器局部网格加密技术对中心水平生产井的多条水力裂缝进行模拟，等效裂缝导流能力为0.1
‑
0.5μm2·
m，等效裂缝网格宽度为0.1m。在致密低渗注气开发模型中，对4口垂直注气井进行水力压裂。每口注气井压裂1条沿笛卡尔坐标系下x方向、缝长为60m固定长度的水力裂缝，等效导流能力为0.1
‑
0.5μm2·
m，等效裂缝网格宽度为0.1m。水力裂缝渗透率采用公式(1)进行计算：
46.k
f
×
w
f
＝k
′
f
×
δ
f
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
47.式中，w
f
和δ
f
分别为水力裂缝、等效水力裂缝网格的宽度，单位为m；k
f
和k
′
f
分别为水力裂缝、等效水力裂缝网格的渗透率，单位为μm2。
48.s2、基于数值模拟基础模型，使用python编程调用eclipse软件建立多套致密低渗油藏注气开发双重介质数值模拟模型，形成致密低渗油藏注气开发数据库，包括：
49.s21、使用allpairs软件设计正交试验，涉及包括地质参数、与流体压力/体积/温度及其性质相关参数、与油田开发决策相关的因素等26个：地层厚度、地层初始压力、地层温度、岩石压缩系数、水平井长度、水平井水力裂缝条数、水平井水力裂缝半长总长度、水平井最小水力裂缝半长、射孔位置、基质孔隙度、天然裂缝孔隙度、天然裂缝渗透率、注气井裂缝渗透率、注气量、注气气体组分比例、原油粘度、束缚水粘度、多个相渗参数(油水系统油相指数、油水系统水相指数、油水系统最大水相相对渗透率、油水系统残余油饱和度、束缚水饱和度、油气系统油相指数、油气系统气相指数、油气系统最大气相相对渗透率、油气系统残余油饱和度)，用于筛选有代表性试验方案；
50.s22、基于s1中使用eclipse建立的基础模型data文件和s21正交试验方案，使用python编程调用eclipse生成多套致密低渗油藏注气开发双重介质数值模拟模型data文件，实现自动生成、运行多套致密低渗油藏注气开发数值模拟模型，分别模拟不同地质\流体\油藏生产过程，具体地，本实施例通过eclipse软件产生例如571种数值模拟方案。
51.s23、使用python自动提取eclipse软件运行多套注气开发数值模拟模型的结果数据，生成规范统一数据格式，包括日产油(fopr)、地层压力(fpr)、采收率(foe)等。
52.s3、基于随机森林算法筛选用于构建日产油(fopr)、地层压力(fpr)、采收率(foe)代理模型相关的特征参数。
53.具体地，本实施例应用随机森林算法研究各个特征参数对日产油(fopr)、地层压力(fpr)、采收率(foe)不同输入变量的影响程度，综合筛选包括地质参数、与流体压力/体积/温度及其性质相关参数、与油田开发决策相关的特征，涉及地层初始压力、岩石压缩系数、裂缝条数、裂缝半长总长度、射孔位置、基质孔隙度、天然裂缝孔隙度、天然裂缝渗透率、
注气井裂缝渗透率、注气量、注气气体组分比例、多个相渗参数(油气系统油相指数、油气系统最大气相相对渗透率、油气系统残余油饱和度、束缚水饱和度)共15个特征参数，作为后续深度学习的输入参数，减小样本过拟合风险。
54.s4、将筛选出的特征参数使用基于均方根反向传播(rmsprop)的长短期记忆单元(lstm)算法，分别建立并训练日产油(fopr)、地层压力(fpr)及采收率(foe)的预测代理模型。
55.具体地，将筛选出的特征参数使用基于均方根反向传播(rmsprop)的长短期记忆单元(lstm)算法，经过minmax归一化，并按照设定要求划分为训练集(90％)和测试集(10％)，分别训练日产油(fopr)、地层压力(fpr)、采收率(foe)预测代理模型，在测试集上进行盲测。
56.本实施例采用的评价标准有拟合系数(r2)、均方根误差(rmse)、平均绝对误差(mae)，如公式(2)、(3)、(4)所示。
[0057][0058][0059][0060]
其中，y为真实值，为预测值，为真实值的平均值，n为样本个数。
[0061]
预测代理模型可以根据以往数据预测未知，准确地复制油藏模拟器的行为。使用长短期记忆单元(lstm)算法训练三种预测代理模型，所使用的超参数如表1所示，预测代理模型研究结果如表2所示，基于lstm算法的预测代理模型真实值预测值对比如图3所示，经过特征参数筛选、正确划分训练集测试集，减小样本过拟合风险后，基于lstm算法的预测代理模型能被用来挖掘多变量油藏数据之间复杂非线性关系，在解决复杂的油藏工程和提高采收率模拟预测问题时，可有效兼顾易用性和精确性。
[0062]
表1神经网络超参数
[0063][0064]
表2预测代理模型研究结果
[0065][0066]
实施例二
[0067]
上述实施例一提供了基于lstm的致密低渗油藏压裂井网气驱效果预测的方法，与之相对应地，本实施例提供一种基于lstm的致密低渗油藏压裂井网气驱效果预测的系统。本实施例提供的系统可以实施实施例一的基于lstm的致密低渗油藏压裂井网气驱效果预测的方法，该系统可以通过软件、硬件或软硬结合的方式来实现。为了描述的方便，描述本实施例时以功能分为各种单元分别描述。当然，在实施时可以把各单元的功能在同一个或多个软件和/或硬件中实现。例如，该系统可以包括集成的或分开的功能模块或功能单元来执行实施例一各方法中的对应步骤。由于本实施例的系统基本相似于方法实施例，所以本实施例描述过程比较简单，相关之处可以参见实施例一的部分说明即可，本发明提供的基于lstm的致密低渗油藏压裂井网气驱效果预测的系统的实施例仅仅是示意性的。
[0068]
本实施例提供的基于lstm的致密低渗油藏压裂井网气驱效果预测系统，包括：
[0069]
基础模型建立单元，被配置为建立致密低渗油藏压裂井网注气开发双重介质数值模拟基础模型；
[0070]
数据库建立单元，被配置为基于数值模拟基础模型，建立多套致密低渗油藏注气开发双重介质数值模拟模型，形成致密低渗油藏注气开发数据库；
[0071]
特征参数筛选单元，被配置为筛选用于构建日产油、地层压力、采收率代理模型相关的特征参数；
[0072]
预测模型建立单元，被配置为将筛选出的特征参数使用基于均方根反向传播的lstm算法，分别建立并训练日产油、地层压力及采收率的预测代理模型。
[0073]
实施例三
[0074]
本实施例提供一种与本实施例一所提供的基于lstm的致密低渗油藏压裂井网气驱效果预测的方法对应的电子设备，电子设备可以是用于客户端的电子设备，例如手机、笔记本电脑、平板电脑、台式机电脑等，以执行实施例一的方法。
[0075]
如图4所示，电子设备包括处理器、存储器、通信接口和总线，处理器、存储器和通信接口通过总线连接，以完成相互间的通信。总线可以是工业标准体系结构(isa，industry standard architecture)总线，外部设备互连(pci，peripheral component)总线或扩展工业标准体系结构(eisa，extended industry standard component)总线等等。存储器中存
储有可在所述处理器上运行的计算机程序，所述处理器运行所述计算机程序时执行本实施例一所提供的基于lstm的致密低渗油藏压裂井网气驱效果预测的方法。本领域技术人员可以理解，图4中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算设备的限定，具体的计算设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
[0076]
在一些实现中，上述的存储器中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本技术的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本技术各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read
‑
only memory)、随机存取存储器(ram，random access memory)、光盘等各种可以存储程序代码的介质。
[0077]
在另一些实现中，处理器可以为中央处理器(cpu)、数字信号处理器(dsp)等各种类型通用处理器，在此不做限定。
[0078]
实施例四
[0079]
本实施例一的基于lstm的致密低渗油藏压裂井网气驱效果预测的方法可被具体实现为一种计算机程序产品，计算机程序产品可以包括计算机可读存储介质，其上载有用于执行本实施例一所述的基于lstm的致密低渗油藏压裂井网气驱效果预测的方法的计算机可读程序指令。
[0080]
计算机可读存储介质可以是保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是但不限于电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意组合。
[0081]
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实现”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0082]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0083]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特
定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0084]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0085]
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。