优化深度学习降维重构参数的油藏自动历史拟合方法技术

本发明专利技术涉及一种优化深度学习降维重构参数的自动历史拟合方法及系统，包括以下步骤：步骤一、构建初始油藏模型，通过自动编码器方法训练深度自编码模型，使用解码器进行降维参数化表征，步骤二、利用粒子群优化算法优化自编码模型的神经网络结构，步骤三、ES‑MDA方法迭代同化数据更新潜变量，自编码模型解码器重构油藏模型并进行数值模拟，步骤四、多次迭代进行观测数据拟合，最终输出拟合模型，完成自动历史拟合流程。相对于现有技术，本发明专利技术通过深度自编码模型、粒子群优化算法的结合可将高维地质特征如渗透率映射到低维连续高斯空间进行重参数化表征，结合集合光滑多次数据同化(ES‑MDA)方法吸收生产历史实现复杂大规模油藏自动历史拟合高效求解。

【技术实现步骤摘要】
优化深度学习降维重构参数的油藏自动历史拟合方法
本专利技术属于地球物理领域,具体地，涉及一种优化深度学习降维重构参数的自动历史拟合方法及系统。
技术介绍
油藏自动历史拟合是油藏开发管理的重要组成部分，基于油藏数值模拟技术通过拟合油气井的生产观测数据对储层数值模型的不确定性参数如渗透率、孔隙度等进行反演调整，为井网规划、生产优化及方案设计等提供可靠的储层数值模型。由于储层地质特征的复杂性，往往需要足够精细的网格构建储层数值模型，导致自动历史拟合需要调整的参数通常达数十万乃至数百万。油藏数值模型降维重构方法通常被用于解决高维网格变量难以直接反演调整的难题。近年来，随着高性能计算机、传感器等广泛应用，计算能力、实时数据等显著增加。因此机器学习方法在图像处理、计算机视觉等领域获得了巨大的成功，在工程领域机器学习方法也显示出巨大的应用潜力。数据驱动求解和机器学习回归预测模型等研究证明了机器学习方法能够有效的应用于复杂系统数值建模、参数反演和优化决策。目前机器学习方法在自动历史拟合领域的应用尚处于起步阶段，但机器学习领域近年来的变革性突破为自动历史拟合技术的更新换代带来了新的契机。
技术实现思路
为克服现有技术存在的缺陷，本专利技术提供一种优化深度学习降维重构参数的油藏自动历史拟合方法，采用自动编码器(Autoencoder,AE)对大规模网格油藏静态参数进行降维与重构处理，同时利用进化算法优化自编码模型的结构，对复杂地质特征进行特征提取及降维分解，特征压缩的同时保留了足够的信息实现保持与复杂地质特征一致的油藏模型重构，结合降维后的低维空间连续特征与数据光滑多次数据同化方法(ES-MDA)进行自动历史拟合，解决了数据同化方法对复杂离散特征难以有效应用的问题。为实现上述目的，本专利技术采用下述方案：优化深度学习降维重构参数的油藏自动历史拟合方法，包括以下步骤：步骤一、构建初始油藏模型，通过自动编码器方法训练深度自编码模型，使用解码器进行降维参数化表征步骤二、利用粒子群优化算法优化自编码模型的神经网络结构步骤三、ES-MDA方法迭代同化数据更新潜变量，自编码模型解码器重构油藏模型并进行数值模拟步骤四、多次迭代进行观测数据拟合，最终输出拟合模型，完成自动历史拟合流程。相对于现有技术，本专利技术的有益效果如下：通过深度自编码模型、粒子群优化算法的结合可将高维地质特征如渗透率映射到低维连续高斯空间进行重参数化表征，结合集合光滑多次数据同化(ES-MDA)方法吸收生产历史实现复杂大规模油藏自动历史拟合高效求解。附图说明图1是优化深度学习降维重构参数的自动历史拟合方法流程示意图；图2A是连续相油藏算例的真实油藏模型；图2B是连续相油藏算例的随机初始油藏模型集合；图3A是优化前后深度自编码(DAE)模型生成渗透率场效果对比的初始油藏模型；图3B为是优化前深度自编码模型重构结果(隐层节点数为519-391-60)；图3C是优化后深度自编码模型重构结果(隐层节点数为3905-256-41)；图4A是ES-MDA结合深度自动编码器-粒子群优化方法生产观测数据拟合结果的生产井1的日产油量，其中圆点表示观测数据点，细实线是初始油藏模型集合数值模拟预测结果，虚线是历史拟合更新模型集合的数值模拟预测结果；图4B是生产井2的日产油量的历史拟合结果；图4C是生产井3的日产油量的历史拟合结果；图4D是生产井4的日产油量的历史拟合结果；具体实施方式参照图1，优化深度学习降维重构参数的油藏自动历史拟合方法，包括以下步骤：步骤一、构建初始油藏模型，通过自动编码器方法训练深度自编码模型，使用解码器进行降维参数化表征，具体过程如下：1.1：构建初始油藏模型将油藏开发生产动态响应dobs和油藏模型参数m看作是油藏系统的输入与输出：dobs＝g(m)+ε(1)其中g(·)表示油藏系统的数值模型或油藏数值模拟器，ε表示生产动态响应的观测误差，使用贝叶斯概率推理方法建立油藏模型参数估计数学模型：其中，p(m)称为先验概率，表示模型参数的已知信息；p(dobs|m)称为似然函数，用于量化观测数据与估计模型参数进行数值模拟预测之间的拟合误差；p(m|dobs)即后验概率，表示在已知观测数据dobs的情形下模型参数m的概率分布；根据油藏的静态观测数据，使用随机地质统计建模方法可以建立大量的随机初始油藏模型mi(i×1，2，...，Ne)，结合高斯概率模型建立油藏模型参数的先验概率表达：假设观测误差ε服从高斯分布，即ε～N(0,CD)，根据式(1)可以推得似然函数：结合式(2)、(3)、(4)可得到油藏模型参数的后验概率表达：1.2：训练自编码模型，使用解码器进行降维参数化表征读取大规模网格油藏静态参数即高维油藏静态参数,并采用自动编码器对所述大规模网格油藏静态参数进行降维,得到降维后的油藏静态参数。用初始油藏模型数据训练深度自编码模型，应根据不同模型设置不同的网络层数，隐层数通常设为3-5层，根据训练效果设置训练迭代次数，通常设为100-200次；构造自动编码器目标函数,并根据所述自动编码器目标函数将自动编码器输入层中的所述高维油藏静态参数压缩至隐藏层并去掉数据中的冗余信息,然后在输出层中对压缩至隐藏层中的数据进行降维,得到降维后的油藏静态参数,其中,所述高维油藏静态参数具体包括各划分网格的渗透率以及孔隙度。步骤二、利用粒子群优化算法优化自编码模型的神经网络结构，过程如下：将生成的深度AE模型网络各层节点个数xi作为粒子群优化算法的优化参数，模型训练后的训练误差作为优化目标函数：粒子群优化算法初始化为一群随机粒子(随机解)，各层节点的初始化范围如表1中所示，然后通过迭代找到最优解。当粒子的训练误差更小时，则该粒子位置作为局部极值pbest,所有粒子的最小训练误差所在位置作为全局极值gbest。在每一次的迭代中，粒子通过跟踪两个“极值”(pbest，gbest)来更新自己。在找到这两个最优值后，粒子通过下面的公式来更新自己的速度和位置。设置适当的终止条件，在到达预设的迭代次数后,优化结束并得到最优组合。vi＝vi+c1×rand()×(phesti-xi)+c2×rand()×(gbesti-xi)(6)xi＝xi+vi(7)其中，i＝1,2,…,N,N是此群中粒子的总数。vi：粒子的速度；rand()：介于(0,1)之间的随机数；xi：粒子的当前位置；c1和c2：学习因子，通常c1＝c2＝2；利用优化后得到的结果作为隐层节点数组合训练深度自编码模型，将模型参数输入进行降维处理得到潜变量。步骤三、ES-MDA方法迭代同化数据更新潜变量，自编码模型解码器重构油藏模型并进行数值模拟，过程如下：3.1ES-MDA方法迭代同化数据更新潜变量：使用E本文档来自技高网...

【技术保护点】
1.优化深度学习降维重构参数的油藏自动历史拟合方法，其特征在于，包括以下步骤：/n步骤一、构建初始油藏模型，通过自动编码器方法训练深度自编码模型，使用解码器进行降维参数化表征/n步骤二、利用粒子群优化算法优化自编码模型的神经网络结构/n步骤三、ES-MDA方法迭代同化数据更新潜变量，自编码模型解码器重构油藏模型并进行数值模拟/n步骤四、多次迭代进行观测数据拟合，最终输出拟合模型，完成自动历史拟合流程。/n

【技术特征摘要】
1.优化深度学习降维重构参数的油藏自动历史拟合方法，其特征在于，包括以下步骤：
步骤一、构建初始油藏模型，通过自动编码器方法训练深度自编码模型，使用解码器进行降维参数化表征
步骤二、利用粒子群优化算法优化自编码模型的神经网络结构
步骤三、ES-MDA方法迭代同化数据更新潜变量，自编码模型解码器重构油藏模型并进行数值模拟
步骤四、多次迭代进行观测数据拟合，最终输出拟合模型，完成自动历史拟合流程。


2.根据权利要求1所述的优化深度学习降维重构参数的油藏自动历史拟合方法，其特征在于，步骤一具体过程如下：
1.1：构建初始油藏模型
将油藏开发生产动态响应dobs和油藏模型参数m看作是油藏系统的输入与输出：
dobs＝g(m)+ε(1)
其中g(·)表示油藏系统的数值模型或油藏数值模拟器，ε表示生产动态响应的观测误差，使用贝叶斯概率推理方法建立油藏模型参数估计数学模型：



其中，p(m)称为先验概率，表示模型参数的已知信息；p(dobs|m)称为似然函数，用于量化观测数据与估计模型参数进行数值模拟预测之间的拟合误差；p(m|dobs)即后验概率，表示在已知观测数据dobs的情形下模型参数m的概率分布；
根据油藏的静态观测数据，使用随机地质统计建模方法可以建立大量的随机初始油藏模型ni(i＝1，2，...，Ne)，结合高斯概率模型建立油藏模型参数的先验概率表达：



假设观测误差ε服从高斯分布，即ε～N(0,CD)，根据式(1)可以推得似然函数：



结合式(2)、(3)、(4)可得到油藏模型参数的后验概率表达：



1.2：训练自编码模型，使用解码器进行降维参数化表征
读取大规模网格油藏静态参数即高维油藏静态参数,并采用自动编码器对所述大规模网格油藏静态参数进行降维,得到降维后的油藏静态参数；用初始油藏模型数据训练深度自编码模型，应根据不同模型设置不同的网络层数，隐层数通常设为3-5层，根据训练效果设置训练迭代次数，通常设为100-200次；
构造自动编码器目标函数,并根据所述自动编码器目标函数将自动编码器输入层中的所述高维油藏静态参数压缩至隐藏层并去掉数据中的冗余信息,然后在输出层中对压缩至隐藏层中的数据进行降维,得到降维后的油藏静态参数,其中,所述...

【专利技术属性】
技术研发人员：张凯，陈昕晟，张世明，刘晨，姚军，吴义志，张黎明，刘均荣，姚传进，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东;37