首页期刊简介编委会编辑部过刊浏览征稿简则期刊信息下载中心欢迎订阅English
引 en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。

基于主成分分析和代理模型的油藏生产注采优化方法

  • 张凯
  • ， 陈国栋
  • ， 薛小明
  • ， 张黎明
  • ， 孙 海
  • ， 姚传进

1. 中国石油大学(华东)石油工程学院,山东青岛 266580；

A reservoir production optimization method based on principal component analysis and surrogate model

  • ZHANG Kai
  • ， CHEN Guodong
  • ， XUE Xiaoming
  • ， ZHANG Liming
  • ， SUN Hai
  • ， YAO Chuanjin

1. School of Petroleum Engineering in China University of Petroleum(East China), Qingdao 266580 , China；

作者简介:

张凯(1980-),男,教授,博士,博士生导师,国家自然科学基金优秀青年基金获得者,山东省"泰山学者"特聘专家,研究方向为油气田开发、油藏工程及优化。E-mail:zhangkai@upc.edu.cn。

中图分类号:TE34

文献标识码:A

文章编号:1673-5005(2020)03-0090-08

DOI:10.3969/j.issn.1673-5005.2020.03.010

参考文献 1
ZHANG Kai,ZHANG Liming,YAO Jun,et al.Water flooding optimization with adjoint model under control constraints [J].Journal of Hydrodynamics,2014,26(1):75-85. 
查找原文
参考文献 2
张凯,李阳,姚军,等.油藏生产优化理论研究[J].石油学报,2010,31(1):78-83.ZHANG Kai,LI Yang,YAO Jun,et al.Theoretical re-search on production optimization of oil reservoirs [J].Acta Petrolei Sinica,2010,31(1):78-83. 
查找原文
参考文献 3
张凯,路然然,张黎明,等.基于序列二次规划算法的油藏动态配产配注优化 [J].油气地质与采收率,2014,21(1):45-50.ZHANG Kai,LU Ranran,ZHANG Liming,et al.Opti-mization of dynamic distribution and distribution of reser-voirs based on sequential quadratic programming algo-rithm[J].Petroleum Geology and Recovery Efficiency,2014,21(1):45-50. 
查找原文
参考文献 4
DEHDARI V,OLIVER D S.Sequential quadratic pro-gramming for solving constrained production optimization:case study from Brugge field[J].SPE Journal,2012,17(3):874-884. 
查找原文
参考文献 5
LI L,JAFARPOUR B,MOHAMMAD-KHANINEZHAD M R.A simultaneous perturbation stochastic approxima-tion algorithm for coupled well placement and control op-timization under geologic uncertainty[J].Computational Geoences,2013,17(1):167-188. 
查找原文
参考文献 6
CHEN Y,OLIVER D S,ZHANG D.Efficient ensemble-based closed-loop production optimization[J].SPE Jour-nal,2009,14(4):634-645. 
查找原文
参考文献 7
TRELEA I C.The particle swarm optimization algorithm:convergence analysis and parameter selection[J].Infor-mation Processing Letters,2003,85(6):317-325. 
查找原文
参考文献 8
ZHANG L,WANG S,ZHANG K,et al.Cooperative ar-tificial bee colony algorithm with multiple populations for interval multiobjective optimization problems [J].IEEE Transactions on Fuzzy Systems,2019,27(5):1052-1065. 
查找原文
参考文献 9
赵辉.油藏开发闭合生产优化理论研究[D].青岛:中国石油大学(华东),2011.ZHAO Hui.Theoretical research on reservoir closed-loop production optimization [D].Qingdao:China University of Petroleum(East China),2011. 
查找原文
参考文献 10
张凯,张秀清,张黎明,等.基于改进的同步扰动随机逼近算法的聚表二元驱优化[J].中国石油大学学报(自然科学版),2017,41(5):102-109.ZHANG Kai,ZHANG Xiuqing,ZHANG Liming,et al.A novel approach for optimization of polymer-surfactant flooding based on simultaneous perturbation stochastic approximation algorithm[J].Journal of China University of Petroleum(Edition of Natural Science),2017,41(5):102-109. 
查找原文
参考文献 11
ZHANG K,ZHANG X,NI W,et al.Nonlinear con-strained production optimization based on augmented La-grangian function and stochastic gradient[J].Journal of Petroleum Science & Engineering,2016,146:418-431. 
查找原文
参考文献 12
ZHOU Z,ONG Y S,NAIR P B,et al.Combining glob-al and local surrogate models to accelerate evolutionary optimization[J].IEEE Transactions on Systems Man & Cybernetics Part C,2006,37(1):66-76. 
查找原文
参考文献 13
HE J,XIE J,WEN X H,et al.An /html/zgsydxxb/20200310/alternative proxy for history matching using proxy-for-data approach and re-duced order modeling[J].Journal of Petroleum Science & Engineering,2016,146:392-399. 
查找原文
参考文献 14
HE J,XIE J,SARMA P,et al.Proxy-based work flow for a priori evaluation of data-acquisition programs[J].SPE Journal,2016,21(4):1400-1412. 
查找原文
参考文献 15
LIAN Y,LIOU M.Multiobjective optimization using coupled response surface model and evolutionary algo-rithm[J].European Radiology,2005,43(6):1316-1325. 
查找原文
参考文献 16
HOROWITZ B,AFONSO S M B,de MENDONA C V P.Surrogate based optimal waterflooding management [J].Journal of Petroleum Science & Engineering,2013,112(sup):206-219. 
查找原文
参考文献 17
GOLZARI A,HAGHIGHAT S M,JAMSHIDI S.Devel-opment of an adaptive surrogate model for production op-timization[J].Journal of Petroleum Science and Engi-neering,2015,133:677-688. 
查找原文
参考文献 18
LUO C,ZHANG S L,WANG C,et al.A metamodel-assisted evolutionary algorithm for expensive optimization [J].Journal of Computational & Applied Mathematics,2011,236(5):759-764. 
查找原文
参考文献 19
LIU B,ZHANG Q,GIELEN G G E.A Gaussian process surrogate model assisted evolutionary algorithm for medium scale expensive optimization problems[J].IEEE Transac-tions on Evolutionary Computation,2014,18(2):180-192. 
查找原文
参考文献 20
REGIS R G.An initialization strategy for high-dimen-sional surrogate-based expensive black-box optimization [J].Springer Proceedings in Mathematics & Statistics,2013,62(5):837-853. 
查找原文
参考文献 21
张凯,姚军,刘顺,等.埕岛油田 6A+B 区块油藏动态调控[J].中国石油大学学报(自然科学版),2009,33(6):71-76.ZHANG Kai,YAO Jun,LIU Shun,et al.Reservoir dy-namic control in block 6A+B of Chengdao Oilfield [J].Journal of China University of Petroleum(Edition of Natural Science),2009,33(6):71-76. 
查找原文
参考文献 22
张凯,路然然,周文胜,等.无梯度多参数自动历史拟合方法[J].中国石油大学学报(自然科学版),2014,38(5):109-115.ZHANG Kai,LU Ranran,ZHOU Wensheng,et al.Multi-parameter gradient-free automatic history matching method [J].Journal of China University of Petroleum(Edition of Natural Science),2014,38(5):109-115. 
查找原文
参考文献 23
LIU Bo,WANG Ling,JIN Yihui.Advances in differen-tial evolution[J].Control & Decision,2007,22(7):721-729. 
查找原文
参考文献 24
QUEIPO N V,HAFTKA R T,SHYY W,et al.Surro-gate-based analysis and optimization [J].Progress in Aerospace Sciences,2005,41(1):1-28. 
查找原文
参考文献 25
BHOSEKAR A,IERAPETRITOU M.Advances in sur-rogate based modeling,feasibility analysis,and optimi-zation:a review[J].Computers & Chemical Engineer-ing,2018,108(4):250-267. 
查找原文
参考文献 26
LIU B,GROUT V,NIKOLAEVA A.Efficient global optimization of actuator based on a surrogate model as-sisted hybrid algorithm[J].IEEE Transactions on In-dustrial Electronics,2018,65:5712-5721. 
查找原文
参考文献 27
张天宇.数据降维算法研究及其应用[D].上海:上海交通大学,2008.ZHANG Tianyu.Research on dimensionality reduction algorithms and its applications [D].Shanghai:Shang-hai Jiao Tong University,2008. 
查找原文
参考文献 28
JANSEN J D,FONSECA R M,KAHROBAEI S,et al.The Egg model:a geological ensemble for reservoir sim-ulation[J].Geoscience Data Journal,2014,1(2):192-195. 
查找原文
参考文献 29
SOREK N,GILDIN E,BOUKOUVALA F,et al.Di-mensionality reduction for production optimization using polynomial approximations [J].Computational Geosci-ences,2017,21(2):247-266.
查找原文
目录contents

    摘要

    针对油藏注采优化变量多,直接使用代理模型优化结果变差的问题,提出基于主成分分析和代理模型的油藏生产注采优化方法,代理模型通过对油藏数值模拟器采样,构建速度快、逼近原数值模拟模型精度的数学模型。 通过主成分分析对数据进行降维,而后借助基于代理模型的差分进化算法进行优化,最终得到实际生产优化问题的调控方案,并应用 SADE-PCA 方法对 Egg 模型进行实例计算。 结果表明:提出的基于主成分分析的代理辅助进化算法,在保证精度的前提下,可有效降低变量的维数,解决油藏大规模变量注采优化的难题;优化的注采方案能够有效增油控水,提高油田的经济效益;插值控制方法(ICM)在实现平滑的调控方案和提高计算效率方面具有良好性能。

    Abstract

    When dealing with reservoir production optimization problems with many variables, directly using surrogate model may deteriorate the optimization results. In this study, a surrogate-assisted evolutionary algorithm based on principal component analysis was proposed for production optimization. Surrogate model is a mathematical model that is computationally fast with an accuracy close to the original numerical simulation model. The principal component analysis was used to reduce the dimension- ality of the data, and then it was optimized by using a differential evolutionary algorithm based on the surrogate model. Conse- quently, the optimal control scheme of an actual production optimization problem was obtained, and the Egg model was adopted to test the proposed method via the SADE-PCA method. The results show that the proposed surrogate-assisted evolutionary algo- rithm based on principal component analysis can effectively reduce the dimension of variables and solve large-scale production optimization of reservoirs efficiently. The optimized injection-production control scheme can effectively increase cumulative oil production and reduce water production to improve the economic benefits of oilfields. In addition, the interpolation control method (ICM) has great performance in achieving smooth control schemes and improving calculation efficiency.

  • 目前,新油田发现的数量逐年递减,世界各地许多大型油田产量也在逐年下降。 如何立足现有条件,通过对各单井流量进行调控,使注入水流入未波及的区域,以提高油藏采收率、最大化油藏开发潜力至关重要。 近年来生产优化领域算法可以分为梯度优化算法和无梯度算法。 基于梯度的优化方法[1-2](如序列二次规划[3-4])使用伴随方法来获得梯度信息。 然而,基于梯度的优化方法容易陷入局部最优, 梯度获取困难。 无梯度方法可以进一步划分为近似梯度算法(如SPSA算法[5] 和En-Opt算法[6])和启发式算法( 如粒子群算法[7] 和人工蜂群算法(ABC) [8])。 近似梯度算法[9-11] 可以计算近似的数值梯度,不需要使用伴随方法。 启发式算法全局性好,可获得更高的油藏开发经济效益,但需要大量数值模拟调用次数。 基于代理模型的进化算法在油藏生产优化问题中可以减少数值模拟调用次数降低时间成本[12]。 因此近年来研究者提出了很多基于代理模型的算法来解决生产优化问题[13-14],如多项式响应面模型[15]、克里金模型[16]、人工神经网络[17]和径向基函数[18]。 然而,代理模型在构建模型时计算量大,处理高维问题时性能恶化[19]。 Regis [20] 提出基于克里金模型构建代理的方法主要应用于少于15 个决策变量的问题。 为解决上述问题,笔者提出结合代理模型和降维算法的新框架解决生产优化问题,并将油藏模型实例结果与现有进化算法和直接代理优化进行对比分析。

  • 1 油藏生产优化问题

  • 油藏生产优化通过调整油水井的注采参数改变地下流体的分布状态,提高油藏波及系数,以最大化目标函数(如净现值(NPV)或累积产油量),进而实现生产效益最大化[21-22]。 由于决策变量和目标函数之间的非线性关系,找到油水井最优注采方案是一项具有挑战性的工作。 生产优化问题的目标是找到最优控制变量组xm×n,其中m为井的数量,n为控制时间步数。 因此变量的数量等于待优化井数乘以控制时间步数。

  • 油田开发的主要目的是经济效益最大化,因此结合最优控制理论,以净现值作为生产优化性能指标,得到简化的目标函数表达式为

  • f(x,y)=NPV(x,y)=∑t=1nΔtQo,tro-Qw,trw-Qi,tri(1+b)pt
    (1)

  • 约束条件为

  • c(x)⩽0
    (2)

  • x1⩽x⩽xu
    (3)

  • 式中,x为要优化的变量组,即井控参数向量(如井底压力或油水井流量);y为状态参数矩阵,如油藏模型所有网格的压力、饱和度和组分等;f(x,y)为目标函数;n为控制时间步数;Qo,t、Q和Qi,t分别为时间步t对应的产油速度、产水速度和注水速度,m 3/d;ro 为原油价格,元/m 3;r和ri 分别为水处理成本和注水成本,元/m 3;b为平均年利率,%;p为经过的时间,a;c(x)为决策变量约束(如现场操作约束、 物理约束等);xl 和xu 分别为下边界和上边界。

  • 2 算法机制

  • 代理模型方法通过少量样本点构建目标函数的近似模型,从而在计算目标函数值时使用代理模型而不是直接调用油藏数值模拟器。 由于代理模型计算目标函数值所需计算代价远小于直接调用油藏数值模拟器计算,因此可以有效提高优化效率。 基于代理模型的方法通常建立在具有强大全局搜索能力的进化算法的基础上。 本文中采用流行的差分进化算法作为优化求解器。

  • 2.1 差分进化算法

  • 差分进化算法(DE) [23]是一种流行的启发式算法,用于解决全局优化问题,并广泛应用于工程优化问题。 经典的DE算法如下:

  • 变异。 本文中选择“DE/best/1” 作为变异策略,

  • vi=xbest+Mu(xi1-xi2)
    (4)

  • 式中,xbest为目前找到的最优解,即最优的井控参数向量;i 1,i 2∈[1,NP ] 为相互不同的随机整数;vi 为第i个个体变异向量;Mu 为一个通常设置为0.5 的比例因子。

  • 交叉。

  • ui(j)={vi(j),Uj(0,1)⩽CR或j=jrandxi(j),Uj(0,1)>CR且j≠jrand
    (5)

  • 式中,Uj(0,1)为从0 到1 均匀分布的任意数;j rand为确保试验向量ui 不重复而随机选择的索引;CR 为交叉率。

  • 选择。

  • xi'={ui,f(ui)<f(xi)xi,f(ui)⩾f(xi)
    (6)

  • 式中,x＇i为新的井控参数目标向量。

  • 关于差分进化算法的更多介绍可以参阅文献[23]。

  • 2.2 代理模型方法

  • 近年来研究者已经提出了许多近似函数构造代理模型[13-14],如多项式响应面模型[15]、克里金模型[16]、人工神经网络[17] 和径向基函数[18]。 克里金模型可以从数学角度提供模型的不确定性[12,16,18], 因此本文中选择克里金模型作为代理模型。

  • 克里金模型(普通克里金法[24])假定确定性函数为

  • fp(x)=μ+ε(x)
    (7)

  • 式中,μ 为样本点的平均常数;ε( x) 为具有零均值的自相关误差项,是样本点之间的广义距离的函数。

  • ε(x)用高斯相关函数作为变差函数的相关结构给出:

  • cov⁡(ε(x(i)),ε(x(j)))=σ2exp⁡(-∑k=1Nφk(x(i)-1)x(j))2)
    (8)

  • 式中,N为x的维数;σ 2为设计点的方差;φk 为第k维的相关系数。

  • 具体而言,在给定一组Ns 数据(x,f)之后,通过最大似然函数估计参数 μ、σ 和 φ [12]。 最优无偏估计为

  • E(fp(x))=μ¯+rTR-1(f-Iμ¯)
    (9)

  • 其中

  • r(x)=[R(x,x(1)),⋯,R(x,x(Ns))]

  • 式中,R为随机过程的相关矩阵; μ¯为最大似然函数估计得到的参数;I为Ns 个元素为1 的向量。

  • 均方误差为

  • s2(x)=σ2[1-rTR-1r+(1-ITR-1r)ITR-1I]
    (10)

  • 选择校正个体时个体的目标函数和代理模型的预测精度都很重要。 为了平衡这两个因素的影响, 已经提出了很多标准,例如预期改进(ExI) [25]、改进概率( PoI) [12] 和最低置信区间( LCB) [26]。 本文中选择LCB策略。 给定预测分布, 预筛选的值定义[26]为

  • flcb(x)=f^(x)-ωs(x)
    (11)

  • 式中, ^ f(x) 为克里金模型的预测值;f lcb( x) 为基于LCB准则的筛选值;ω为平衡个体适应度和不确定性的常数。

  • 基于LCB动态加点策略的代理模型最小化目标函数过程如图1 所示。 首先通过试验设计方法取几个样本点作为初始样本点,然后构建代理模型来近似目标函数,根据LCB策略使用差分进化算法筛选出下一步加点的所在位置,逐步动态取样完善代理模型精度,逼近全局最优点。

  • 图1 基于LCB动态加点策略的代理模型优化过程示意图

  • Fig.1 Schematic diagram of surrogate-assisted optimization process based on LCB strategy

  • 2.3 降维方法主成分分析

  • 由于油藏注采优化问题变量多,直接用克里金插值构建的近似模型误差大,难以实现对数值模拟器进行有效的近似。 因此采用主成分分析方法对优化问题进行降维处理。 主成分分析( PCA) [27] 是一种常用的线性降维方法,可以将高维数据转换为低维数据。 在克里金代理模型中,输入变量过多不仅会增加代理计算的复杂性,代理模型的精度也会下降,影响优化的效率和精度。 因此本文中采用主成分分析方法对数据进行降维。

  • PCA假设数据之间为线性关系,在保存原始数据的协方差结构的基础上计算低维表达,即最大化总体方差。 PCA的目标函数为

  • UPCA=argmaxOPCA∑i=1n∥Yi-Y¯∥2=arg⁡maxDPCA∑i=1N∥UPCAT(Xi-X¯)∥2
    (12)

  • 其中

  • X¯=1N∑i=1NXi,Y¯=1N∑i=1NYi，

  • Xi=(xi1,xi2,⋯,xir),Yi=(yi1,yi2,⋯,yir).

  • 转换矩阵存在尺度约束,U T PCAUPCA=Id,Id 为d ×d单位矩阵。 将原来的数据集X=(X1,X2,…,XN) T 转换为新的数据集Y=(Y1,Y2,…,YN) T,根据在总方差中的贡献率选择Yi,从而可以选择前m个主成分代替原数据集。 主成分分析通过对高维优化空间进行降维,将高维油藏优化数据在低维表示,进而提高代理模型精度及优化效率。

  • 2.4 SADE-PCA算法优化计算步骤

  • 代理模型的优化过程比传统的进化算法更有效,但当处理油藏注采优化问题时,由于优化变量过多导致参数优化搜索空间巨大,代理模型对油藏数值模拟器的近似能力大大降低。 当直接使用基于代理模型的差分进化算法( SADE) 进行优化时,由于变量过多,直接采用代理模型进行优化会导致参数优化搜索空间巨大难以构建高精度的代理模型,即少量样本点无法建立足够精度的代理模型。 为了利用低维空间中的基于代理模型的优化方法,在建立代理模型之前,利用主成分分析将数据从高维空间转换为低维空间。 随后,利用低维数据构建更高精度的克里金代理模型。 因此搜索可以集中在由代理模型引导的适应度高的小区域。 SADE-PCA算法的流程如图2 所示。

  • 图2 SADE-PCA流程

  • Fig.2 Flow diagram of SADE-PCA

  • 3 计算实例

  • 基于油藏数值模拟软件Eclipse,以油藏实例Egg模型[28] 为例,应用SADE-PCA算法进行测试, 并与传统进化算法和直接代理优化算法进行比较。 Egg模型是一个在水驱生产下的非均质油藏模型, 有18553 个有效网格,共有8 口注水井和4 口生产井。 模型渗透率及井位分布如图3 所示。 模型详细说明见文献[28]。 生产井以39.5 MPa的井底压力定压生产。 需要优化的控制变量为8 口注水井的注水速率。 每口井的最大允许注水速度为79.5 m 3/d,最小速度为0 m 3/d。 原油价格为80.0 元/m 3,注水成本为5.0 元/m 3,水处理成本为19.0 元/m 3,平均年利率为0%。

  • 图3 Egg模型渗透率分布

  • Fig.3 Permeability distribution of Egg model

  • 要优化的变量数为注水井数乘以控制步数。 一般控制步数越多,优化越有效,即获得更高的净现值。 控制步数越多,则具有更大的自由度,优化问题能够实现更高的净现值。 然而,控制变量数增多意味着迭代次数增加才能收敛,优化算法需要更多的计算资源才能收敛。 Sorek等[29]引入了一种函数控制方法( FCM) 和插值控制方法( ICM) 以减少计算量,实现适合油田实际应用的平滑注采调控。 ICM方法通过选择几个等距插值点作为优化变量以表示沿时间范围的井控变量,其他井控变量通过三次样条插值确定,ICM方法示意图如图4 所示。 本文中利用插值控制方法缓解由控制时间步数量引起的优化效率和精度的矛盾,减少优化变量数目,并提供平滑的井控方案。 Egg模型的总生产周期为3600 d, 模拟生产时间步长设为90 d。 因此需要优化的变量总数为8×40=320。 选择对每口井用10 个插值点进行控制优化,控制变量的数量减少到8×10=80。

  • 图4 插值控制方法示意图

  • Fig.4 Schematic diagram of interpolation control method

  • 为了比较不同方法的结果,分别用SPSA、EnOpt、DE、DE+ICM、SADE+ICM和SADE-PCA+ICM方法对Egg模型进行优化,其中SPSA和En-Opt算法属于近似梯度类算法。 由于近似梯度类算法的优化结果极大取决于初始点的选取,对程序进行5 次运行,取平均结果。 净现值随数值模拟调用次数变化结果如图5 所示。 在有限的数值模拟调用次数下,近似梯度类方法(SPSA、En-Opt算法)相比差分进化算法表现出明显的优势。 DE方法经过800 次迭代优化后,从初始的1.867×10 7 元增长到2.267× 10 7 元。 使用ICM方法后,净现值增长到2.503×10 7 元,相比直接用DE方法优化,提高了10.41%。 由于代理模型在变量多的情况下无法实现对油藏数值模拟器进行有效近似,SADE+ICM方法甚至比DE+ ICM方法的最终优化结果差。 这说明代理模型在优化过程中起到了反作用,代理模型在筛选过程中误导了差分进化算法产生的子解。 这是因为代理模型在高维空间中失去了准确性,误差大。 使用降维方法后,SADE-PCA+ICM方法与DE+ICM方法相比优化效果改善,说明代理模型在优化过程中对子解的筛选过程促进了优化进程,改善了优化效率。 与DE优化结果相比,SADE-PCA+ICM在相同的数值模拟调用次数下优化出的净现值增长到2.720×10 7 元, 提高了19.98%。 因此从优化结果和优化效率来看,SADE-PCA+ICM方法收敛性最好,能够满足油藏生产优化要求,同时为三次采油油藏开发方案的优化提供了新思路。

  • 图5 各算法的净现值优化结果

  • Fig.5 Optimization result of net present value with different algorithms

  • 各优化方法最终的油藏累积产油量、累积注水量、累积产水量和含水率随时间变化如图6 所示。 对比几种不同方法优化结果,SADE-PCA+ICM提供的方案有效地降低了油井水的产出量,提高了波及系数,有效改善了油田经济开发效益。 各优化方法计算所得的对8 口注水井的最终生产调控如图7 所示。 图7 中1 控制步为90 d,不使用ICM方法的DE的注水井调控方案如图7(a)所示,其同一口井相邻两控制步的数值相差较大,调控方案不稳定且振荡, 不利于现场实施,且注水井注入量变化过大会导致井底压力变化过大,损害储层,导致出砂等现象。 利用ICM的优化算法可以提供平滑的生产调控方案。

  • 图6 各算法优化后油藏累积产油量、累积注水量、累积产水量和油田含水率随时间变化Fig.6 Variation of accumulative oil production, water injection, water production and water cut of reservoir with time after different algorithms optimization

  • 图7 各算法注水井调控

  • Fig.7 Water injection well controls for different algorithms

  • 4 结束语

  • 针对油藏注采优化变量多,直接使用代理模型优化结果变差的问题,提出了SADE-PCA方法,通过主成分分析对数据进行降维,而后借助基于代理模型的差分进化算法进行优化,最终得到实际生产优化问题的调控方案。 应用SADE-PCA方法对Egg模型进行了实例计算。 优化所得的最优调控方案通过对生产制度进行调控,有效地降低了油井水的产出量,增大了原油产量,提高了油田经济开发效益, 进一步验证了方法的有效性。 插值控制方法(ICM) 在实现平滑的调控方案和提高计算效率方面具有良好性能。 使用ICM方法后,代理模型结合主成分分析的方法可以更有效地解决实际油藏注采优化问题。

  • 参考文献

    • [1] ZHANG Kai,ZHANG Liming,YAO Jun,et al.Water flooding optimization with adjoint model under control constraints [J].Journal of Hydrodynamics,2014,26(1):75-85. 

    • [2] 张凯,李阳,姚军,等.油藏生产优化理论研究[J].石油学报,2010,31(1):78-83.ZHANG Kai,LI Yang,YAO Jun,et al.Theoretical re-search on production optimization of oil reservoirs [J].Acta Petrolei Sinica,2010,31(1):78-83. 

    • [3] 张凯,路然然,张黎明,等.基于序列二次规划算法的油藏动态配产配注优化 [J].油气地质与采收率,2014,21(1):45-50.ZHANG Kai,LU Ranran,ZHANG Liming,et al.Opti-mization of dynamic distribution and distribution of reser-voirs based on sequential quadratic programming algo-rithm[J].Petroleum Geology and Recovery Efficiency,2014,21(1):45-50. 

    • [4] DEHDARI V,OLIVER D S.Sequential quadratic pro-gramming for solving constrained production optimization:case study from Brugge field[J].SPE Journal,2012,17(3):874-884. 

    • [5] LI L,JAFARPOUR B,MOHAMMAD-KHANINEZHAD M R.A simultaneous perturbation stochastic approxima-tion algorithm for coupled well placement and control op-timization under geologic uncertainty[J].Computational Geoences,2013,17(1):167-188. 

    • [6] CHEN Y,OLIVER D S,ZHANG D.Efficient ensemble-based closed-loop production optimization[J].SPE Jour-nal,2009,14(4):634-645. 

    • [7] TRELEA I C.The particle swarm optimization algorithm:convergence analysis and parameter selection[J].Infor-mation Processing Letters,2003,85(6):317-325. 

    • [8] ZHANG L,WANG S,ZHANG K,et al.Cooperative ar-tificial bee colony algorithm with multiple populations for interval multiobjective optimization problems [J].IEEE Transactions on Fuzzy Systems,2019,27(5):1052-1065. 

    • [9] 赵辉.油藏开发闭合生产优化理论研究[D].青岛:中国石油大学(华东),2011.ZHAO Hui.Theoretical research on reservoir closed-loop production optimization [D].Qingdao:China University of Petroleum(East China),2011. 

    • [10] 张凯,张秀清,张黎明,等.基于改进的同步扰动随机逼近算法的聚表二元驱优化[J].中国石油大学学报(自然科学版),2017,41(5):102-109.ZHANG Kai,ZHANG Xiuqing,ZHANG Liming,et al.A novel approach for optimization of polymer-surfactant flooding based on simultaneous perturbation stochastic approximation algorithm[J].Journal of China University of Petroleum(Edition of Natural Science),2017,41(5):102-109. 

    • [11] ZHANG K,ZHANG X,NI W,et al.Nonlinear con-strained production optimization based on augmented La-grangian function and stochastic gradient[J].Journal of Petroleum Science & Engineering,2016,146:418-431. 

    • [12] ZHOU Z,ONG Y S,NAIR P B,et al.Combining glob-al and local surrogate models to accelerate evolutionary optimization[J].IEEE Transactions on Systems Man & Cybernetics Part C,2006,37(1):66-76. 

    • [13] HE J,XIE J,WEN X H,et al.An /html/zgsydxxb/20200310/alternative proxy for history matching using proxy-for-data approach and re-duced order modeling[J].Journal of Petroleum Science & Engineering,2016,146:392-399. 

    • [14] HE J,XIE J,SARMA P,et al.Proxy-based work flow for a priori evaluation of data-acquisition programs[J].SPE Journal,2016,21(4):1400-1412. 

    • [15] LIAN Y,LIOU M.Multiobjective optimization using coupled response surface model and evolutionary algo-rithm[J].European Radiology,2005,43(6):1316-1325. 

    • [16] HOROWITZ B,AFONSO S M B,de MENDONA C V P.Surrogate based optimal waterflooding management [J].Journal of Petroleum Science & Engineering,2013,112(sup):206-219. 

    • [17] GOLZARI A,HAGHIGHAT S M,JAMSHIDI S.Devel-opment of an adaptive surrogate model for production op-timization[J].Journal of Petroleum Science and Engi-neering,2015,133:677-688. 

    • [18] LUO C,ZHANG S L,WANG C,et al.A metamodel-assisted evolutionary algorithm for expensive optimization [J].Journal of Computational & Applied Mathematics,2011,236(5):759-764. 

    • [19] LIU B,ZHANG Q,GIELEN G G E.A Gaussian process surrogate model assisted evolutionary algorithm for medium scale expensive optimization problems[J].IEEE Transac-tions on Evolutionary Computation,2014,18(2):180-192. 

    • [20] REGIS R G.An initialization strategy for high-dimen-sional surrogate-based expensive black-box optimization [J].Springer Proceedings in Mathematics & Statistics,2013,62(5):837-853. 

    • [21] 张凯,姚军,刘顺,等.埕岛油田 6A+B 区块油藏动态调控[J].中国石油大学学报(自然科学版),2009,33(6):71-76.ZHANG Kai,YAO Jun,LIU Shun,et al.Reservoir dy-namic control in block 6A+B of Chengdao Oilfield [J].Journal of China University of Petroleum(Edition of Natural Science),2009,33(6):71-76. 

    • [22] 张凯,路然然,周文胜,等.无梯度多参数自动历史拟合方法[J].中国石油大学学报(自然科学版),2014,38(5):109-115.ZHANG Kai,LU Ranran,ZHOU Wensheng,et al.Multi-parameter gradient-free automatic history matching method [J].Journal of China University of Petroleum(Edition of Natural Science),2014,38(5):109-115. 

    • [23] LIU Bo,WANG Ling,JIN Yihui.Advances in differen-tial evolution[J].Control & Decision,2007,22(7):721-729. 

    • [24] QUEIPO N V,HAFTKA R T,SHYY W,et al.Surro-gate-based analysis and optimization [J].Progress in Aerospace Sciences,2005,41(1):1-28. 

    • [25] BHOSEKAR A,IERAPETRITOU M.Advances in sur-rogate based modeling,feasibility analysis,and optimi-zation:a review[J].Computers & Chemical Engineer-ing,2018,108(4):250-267. 

    • [26] LIU B,GROUT V,NIKOLAEVA A.Efficient global optimization of actuator based on a surrogate model as-sisted hybrid algorithm[J].IEEE Transactions on In-dustrial Electronics,2018,65:5712-5721. 

    • [27] 张天宇.数据降维算法研究及其应用[D].上海:上海交通大学,2008.ZHANG Tianyu.Research on dimensionality reduction algorithms and its applications [D].Shanghai:Shang-hai Jiao Tong University,2008. 

    • [28] JANSEN J D,FONSECA R M,KAHROBAEI S,et al.The Egg model:a geological ensemble for reservoir sim-ulation[J].Geoscience Data Journal,2014,1(2):192-195. 

    • [29] SOREK N,GILDIN E,BOUKOUVALA F,et al.Di-mensionality reduction for production optimization using polynomial approximations [J].Computational Geosci-ences,2017,21(2):247-266.

  • 基本信息

    中图分类号: TE34
    文献标识码: A
    DOI: 10.3969/j.issn.1673-5005.2020.03.010
    文章编号: 1673-5005(2020)03-0090-08

    基金信息

    国家自然科学基金项目(51722406,51874335,51674280)；
    山东省自然科学基金项目(JQ201808,ZR2019JQ21)；
    中石油重大科技 项目(ZD2019-183-008)；
    山东省高等学校青创科技支持计划(2019KJH002)；
    国家油气重大专项(2016ZX05025001-006)；

    引用信息

    张凯,陈国栋,薛小明,等. 基于主成分分析和代理模型的油藏生产注采优化方法[J]. 中国石油大学学报 (自然科学版),2020,44(3):90-97.

    ZHANG Kai, CHEN Guodong, XUE Xiaoming, et al. A reservoir production optimization method based on principal compo- nent analysis and surrogate model[J]. Journal of China University of Petroleum(Edition of Natural Science),2020,44(3): 90-97.

    稿件历史

    收稿日期: 2019-07-25

  • 参考文献

    • [1] ZHANG Kai,ZHANG Liming,YAO Jun,et al.Water flooding optimization with adjoint model under control constraints [J].Journal of Hydrodynamics,2014,26(1):75-85. 

    • [2] 张凯,李阳,姚军,等.油藏生产优化理论研究[J].石油学报,2010,31(1):78-83.ZHANG Kai,LI Yang,YAO Jun,et al.Theoretical re-search on production optimization of oil reservoirs [J].Acta Petrolei Sinica,2010,31(1):78-83. 

    • [3] 张凯,路然然,张黎明,等.基于序列二次规划算法的油藏动态配产配注优化 [J].油气地质与采收率,2014,21(1):45-50.ZHANG Kai,LU Ranran,ZHANG Liming,et al.Opti-mization of dynamic distribution and distribution of reser-voirs based on sequential quadratic programming algo-rithm[J].Petroleum Geology and Recovery Efficiency,2014,21(1):45-50. 

    • [4] DEHDARI V,OLIVER D S.Sequential quadratic pro-gramming for solving constrained production optimization:case study from Brugge field[J].SPE Journal,2012,17(3):874-884. 

    • [5] LI L,JAFARPOUR B,MOHAMMAD-KHANINEZHAD M R.A simultaneous perturbation stochastic approxima-tion algorithm for coupled well placement and control op-timization under geologic uncertainty[J].Computational Geoences,2013,17(1):167-188. 

    • [6] CHEN Y,OLIVER D S,ZHANG D.Efficient ensemble-based closed-loop production optimization[J].SPE Jour-nal,2009,14(4):634-645. 

    • [7] TRELEA I C.The particle swarm optimization algorithm:convergence analysis and parameter selection[J].Infor-mation Processing Letters,2003,85(6):317-325. 

    • [8] ZHANG L,WANG S,ZHANG K,et al.Cooperative ar-tificial bee colony algorithm with multiple populations for interval multiobjective optimization problems [J].IEEE Transactions on Fuzzy Systems,2019,27(5):1052-1065. 

    • [9] 赵辉.油藏开发闭合生产优化理论研究[D].青岛:中国石油大学(华东),2011.ZHAO Hui.Theoretical research on reservoir closed-loop production optimization [D].Qingdao:China University of Petroleum(East China),2011. 

    • [10] 张凯,张秀清,张黎明,等.基于改进的同步扰动随机逼近算法的聚表二元驱优化[J].中国石油大学学报(自然科学版),2017,41(5):102-109.ZHANG Kai,ZHANG Xiuqing,ZHANG Liming,et al.A novel approach for optimization of polymer-surfactant flooding based on simultaneous perturbation stochastic approximation algorithm[J].Journal of China University of Petroleum(Edition of Natural Science),2017,41(5):102-109. 

    • [11] ZHANG K,ZHANG X,NI W,et al.Nonlinear con-strained production optimization based on augmented La-grangian function and stochastic gradient[J].Journal of Petroleum Science & Engineering,2016,146:418-431. 

    • [12] ZHOU Z,ONG Y S,NAIR P B,et al.Combining glob-al and local surrogate models to accelerate evolutionary optimization[J].IEEE Transactions on Systems Man & Cybernetics Part C,2006,37(1):66-76. 

    • [13] HE J,XIE J,WEN X H,et al.An /html/zgsydxxb/20200310/alternative proxy for history matching using proxy-for-data approach and re-duced order modeling[J].Journal of Petroleum Science & Engineering,2016,146:392-399. 

    • [14] HE J,XIE J,SARMA P,et al.Proxy-based work flow for a priori evaluation of data-acquisition programs[J].SPE Journal,2016,21(4):1400-1412. 

    • [15] LIAN Y,LIOU M.Multiobjective optimization using coupled response surface model and evolutionary algo-rithm[J].European Radiology,2005,43(6):1316-1325. 

    • [16] HOROWITZ B,AFONSO S M B,de MENDONA C V P.Surrogate based optimal waterflooding management [J].Journal of Petroleum Science & Engineering,2013,112(sup):206-219. 

    • [17] GOLZARI A,HAGHIGHAT S M,JAMSHIDI S.Devel-opment of an adaptive surrogate model for production op-timization[J].Journal of Petroleum Science and Engi-neering,2015,133:677-688. 

    • [18] LUO C,ZHANG S L,WANG C,et al.A metamodel-assisted evolutionary algorithm for expensive optimization [J].Journal of Computational & Applied Mathematics,2011,236(5):759-764. 

    • [19] LIU B,ZHANG Q,GIELEN G G E.A Gaussian process surrogate model assisted evolutionary algorithm for medium scale expensive optimization problems[J].IEEE Transac-tions on Evolutionary Computation,2014,18(2):180-192. 

    • [20] REGIS R G.An initialization strategy for high-dimen-sional surrogate-based expensive black-box optimization [J].Springer Proceedings in Mathematics & Statistics,2013,62(5):837-853. 

    • [21] 张凯,姚军,刘顺,等.埕岛油田 6A+B 区块油藏动态调控[J].中国石油大学学报(自然科学版),2009,33(6):71-76.ZHANG Kai,YAO Jun,LIU Shun,et al.Reservoir dy-namic control in block 6A+B of Chengdao Oilfield [J].Journal of China University of Petroleum(Edition of Natural Science),2009,33(6):71-76. 

    • [22] 张凯,路然然,周文胜,等.无梯度多参数自动历史拟合方法[J].中国石油大学学报(自然科学版),2014,38(5):109-115.ZHANG Kai,LU Ranran,ZHOU Wensheng,et al.Multi-parameter gradient-free automatic history matching method [J].Journal of China University of Petroleum(Edition of Natural Science),2014,38(5):109-115. 

    • [23] LIU Bo,WANG Ling,JIN Yihui.Advances in differen-tial evolution[J].Control & Decision,2007,22(7):721-729. 

    • [24] QUEIPO N V,HAFTKA R T,SHYY W,et al.Surro-gate-based analysis and optimization [J].Progress in Aerospace Sciences,2005,41(1):1-28. 

    • [25] BHOSEKAR A,IERAPETRITOU M.Advances in sur-rogate based modeling,feasibility analysis,and optimi-zation:a review[J].Computers & Chemical Engineer-ing,2018,108(4):250-267. 

    • [26] LIU B,GROUT V,NIKOLAEVA A.Efficient global optimization of actuator based on a surrogate model as-sisted hybrid algorithm[J].IEEE Transactions on In-dustrial Electronics,2018,65:5712-5721. 

    • [27] 张天宇.数据降维算法研究及其应用[D].上海:上海交通大学,2008.ZHANG Tianyu.Research on dimensionality reduction algorithms and its applications [D].Shanghai:Shang-hai Jiao Tong University,2008. 

    • [28] JANSEN J D,FONSECA R M,KAHROBAEI S,et al.The Egg model:a geological ensemble for reservoir sim-ulation[J].Geoscience Data Journal,2014,1(2):192-195. 

    • [29] SOREK N,GILDIN E,BOUKOUVALA F,et al.Di-mensionality reduction for production optimization using polynomial approximations [J].Computational Geosci-ences,2017,21(2):247-266.

    • 版权所有 中国石油大学学报（自然科学版）编辑部 Copyright©2008 All Rights Reserved
    主管单位：中华人民共和国教育部 主办单位：中国石油大学（华东）
    地址： 青岛市黄岛区长江西路66号中国石油大学期刊社 邮编：266580 电话：0532-86983553 E-mail： journal@upc.edu.cn
    本系统由：北京勤云科技发展有限公司设计