CN105913152A - 基于混合遗传算法的油藏生产过程设定点优化方法 - Google Patents
基于混合遗传算法的油藏生产过程设定点优化方法 Download PDFInfo
- Publication number
- CN105913152A CN105913152A CN201610226697.3A CN201610226697A CN105913152A CN 105913152 A CN105913152 A CN 105913152A CN 201610226697 A CN201610226697 A CN 201610226697A CN 105913152 A CN105913152 A CN 105913152A
- Authority
- CN
- China
- Prior art keywords
- oil reservoir
- production process
- oil
- process set
- producing well
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Genetics & Genomics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Development Economics (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于混合遗传算法的油藏生产过程设定点优化方法,属于油藏生产过程优化领域,由两层循环形式组成,内层循环中由改进遗传算法对界约束子问题进行全局寻优以得到下一迭代点,外层循环则修正乘子向量和罚参数向量、检查收敛准则是否满足、重新构造界约束子问题或在收敛准则满足时终止。本发明能有效地处理复杂的非线性约束优化问题,能够对油藏生产过程设定点更加有效的进行优化。
Description
技术领域
本发明属于油藏生产过程优化领域,特别是涉及一种基于混合遗传算法的油藏生产过程设定点优化方法。
背景技术
包括油藏生产过程设定点优化在内的复杂工业生产过程优化与控制方法研究一直是一个热门领域,蕴藏着极大的经济效益,具有广阔的应用前景,其发展对我国国民经济意义重大。
注水是油藏采油过程中最经济地提高采收率的开采方式之一,同时也是我国大部分油藏主要采油方式之一。在油藏注水采油生产过程中,一般在油藏地面的适当位置挖掘有若干个注水井和生产井,每个注水井上设有一套注水控制系统,每个生产井上也设有一套生产控制系统,以保证向各个注水井以适当速率注水并从各个生产井以适当速率抽取油水混合物,从而将油藏中的石油开采出来。
过程优化控制理论与方法一直是一个热门领域,如何寻求工艺变量的最优设定值并实时地对系统进行设定以使生产过程最优地运行,是当今国际复杂工业生产过程优化与控制领域的研究热点之一。油藏生产过程长期以来积累了大量的注水井的注入率和生产井的产出率等历史数据,通过这些数据的研究与分析,可以得到有关各注水井和生产井之间的连通信息,实时优化各注水井的注入率和各生产井的产出率,并及时调整各注水井和各生产井对应控制系统的设定点,可使石油开采过程保持在最优工况下进行。
然而,很少有文献报道对油藏生产过程中各注入率和产出率的设定点优化,传统的基于梯度信息的优化方法难以对油藏生产过程设定点优化问题进行有效的求解。
发明内容
有鉴于现有技术的上述缺陷,本发明所要解决的技术问题是提供一种更加有效的油藏生产过程设定点优化方法。
为实现上述目的,本发明提供了一种基于混合遗传算法的油藏生产过程设定点优化方法,按以下步骤进行:
步骤一、设定在油藏地面挖掘m个注水井和n个生产井,设定生产井j的总产出率为qj,m、n均为正整数;
建立油藏容量模型:其中ik是注水井k的注水率,τj为生产井j的时间常数,λkj为注水井k与生产井j之间的权重,
设定生产井j在t时刻的石油产出率为qcj,建立油藏生产率模型求取aj和bj;
步骤二、建立油藏生产过程设定点优化模型:p0当前石油销售价格,pwj为从油水混合液中分离水的单位费用,qwj为生产井j在t时刻的水产出率;
步骤三、以油藏容量模型、油藏生产率模型为约束条件,以所述R为目标函数,得到油藏生产过程设定点优化问题:
步骤四、对于油藏生产过程设定点优化问题,将其离散化为一般的非线性约束优化问题:
步骤五、对z′(x),设定Lagrange乘子向量λk、和罚参数向量σk、增广Lagrange乘子法第k步所求得子问题为:
其中是修改增广Lagrange函数;x为决策变量,λ为拉格朗日乘子参数;σ为罚函数参数;
其中,
步骤六、由遗传算法对步骤五所述子问题进行全局寻优以得到下一迭代点,然后修正乘子向量和罚参数向量,检查收敛准则是否满足,当收敛准则不满足时重新构造界约束子问题,收敛准则满足时终止算法。
较佳的,步骤六所述遗传算法在进化过程中随机组合一个交叉算子和一个变异算子产生若干个新的个体。
本发明的有益效果是:本发明能有效地处理复杂的非线性约束优化问题,能够对油藏生产过程设定点更加有效的进行优化。
附图说明
图1是本发明一具体实施方式的流程示意图。
图2是组合交叉算子和变异算子示意图。
图3是第一种生产井的实际石油产出率与石油生产率模型计算所得到的石油产出率的比较示意图。
图4是第二种生产井的实际石油产出率与石油生产率模型计算所得到的石油产出率的比较示意图。
图5是第三种生产井的实际石油产出率与石油生产率模型计算所得到的石油产出率的比较示意图。
图6是第四种生产井的实际石油产出率与石油生产率模型计算所得到的石油产出率的比较示意图。
图7是第一种生产井的实际总产出率和容量模型计算值的比较示意图。
图8是第二种生产井的实际总产出率和容量模型计算值的比较示意图。
图9是第三种生产井的实际总产出率和容量模型计算值的比较示意图。
图10是第四种生产井的实际总产出率和容量模型计算值的比较示意图。
图11是p0=60,pw=1,r=0.0时的注入率最优设定值示意图。
图12是p0=60,pw=1,r=0.2时的注入率最优设定值示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明:
如图1至图12所示,一种基于混合遗传算法的油藏生产过程设定点优化方法,由两层循环形式组成,内层循环中由改进遗传算法对子问题进行全局寻优以得到下一迭代点,外层循环则修正乘子向量和罚参数向量、检查收敛准则是否满足、重新构造界约束子问题或在收敛准则满足时终止算法,按以下步骤进行:
步骤一、设定在油藏地面挖掘m个注水井和n个生产井,设定生产井j的总产出率为qj,m、n均为正整数;
建立油藏容量模型:其中ik是注水井k的注水率,τj为生产井j的时间常数,λkj为注水井k与生产井j之间的权重,
设定生产井j在t时刻的石油产出率为qcj,建立油藏生产率模型求取aj和bj;
步骤二、建立油藏生产过程设定点优化模型:p0当前石油销售价格,pwj为从油水混合液中分离水的单位费用,qwj为生产井j在t时刻的水产出率;
步骤三、以油藏容量模型、油藏生产率模型为约束条件,以所述R为目标函数,得到油藏生产过程设定点优化问题:
步骤四、对于油藏生产过程设定点优化问题,将其离散化为一般的非线性约束优化问题:
步骤五、对z′(x),设定Lagrange乘子向量λk、和罚参数向量σk、增广Lagrange乘子法第k步所求得子问题为:
其中是修改增广Lagrange函数;
其中,
步骤六、由遗传算法对步骤五所述子问题进行全局寻优以得到下一迭代点,然后修正乘子向量和罚参数向量,检查收敛准则是否满足,当收敛准则不满足时重新构造界约束子问题,收敛准则满足时终止算法。
本实施例中,步骤六所述遗传算法在进化过程中随机组合一个交叉算子和一个变异算子产生若干个新的个体。
本发明由软件Eclipse对非均质油藏Synfield进行模拟所得到的数据,用来验证本文所提出的油藏生产过程设定点优化模型以及所提出的混合优化方法的有效性。
Synfield油藏是一个有5个注水井(分别用I01,I02,I03,I04和I05表示)和4个生产井(分别用P01,P02,P03和P04表示)的多层各向异性油藏,油藏中水平方向的渗透率是竖直方向的3倍。用Eclipse对Synfield模拟时,时间长度为100个月,即用5个注水井的100个月注入率模拟产生4个生产井相应的100个月的产出率数据.将1-50个月的数据获得设定点优化模型中的参数,然后用混合优化方法求解51-100个月时段的油藏生产过程设定点优化模型,得到5个注水井在该时段内每个月的最优注入率设定值,再以这些最优设定值用Eclipse对Synfield进行模拟产生4个生产井对应的新50组产出率数据,最后将新数据与对比数据进行比较,分析油藏生产过程设定值优化前后的经济效益.
由油藏容量模型和1-50个月的数据得到了油藏中各个井之间的连通性权重及每个生产井的时间常数,如表1所示.
表1由模型和数据所得权重及时间常数
从表1可知,对于非均质油藏,对称注水井和生产井之间的权重是不对称的,权重的大小恰当地量化了各注水井与生产井之间连通性大小。
由图7至图10可知,容量模型能较好地描述注水油藏中水的衰减和滞后特性,以及各注水井注入率和各生产井产出率之间的相互关系.生产井P01,P02,P03和P04的总产出率和容量模型计算量之间拟合的确定性系数R2分别为0.9959、0.9970、0.9973和0.9991。基于前面所得到的权重λ和模型数据,对油藏生产率模型(3)中参数aj和bj的估计结果如表2.
表2由模型数据所得参数aj和bj
由图3至图6可知,4个生产井的石油产出率都拟合的很好,由石油生产率模型计算所得到的石油产出率与其实际石油产出率拟合的确定性系数R2分别为0.9629、0.9590、0.9608和0.9698,说明所建立的石油生产率模型较准确地描述了它们的石油产出率,可对它们的石油产出率较准确预测。利用油藏生产过程设定点优化模型求各注水井的注入率和各生产井的产出率的最优设定值,对每个注水井,本实施例设其注入率下限为零,上限分别为其100个月数据中的最大值(RB/天),即lk=0,k=1,2,3,4,5,u1=3603.739,u2=2000.666,u3=1464.216,u4=1351.601,u5=2220.153.
本实施例优化时间段从第51个月开始到第100个月,以月为步长离散化问题z(x)得到下面非线性约束优化问题:
其中,优化变量为:ik(l),qj(l),k=1,2,3,4,5,j=1,2,3,4,l=51,52,…,100,共450个,其参数设置为:在外层迭代中,容许误差ε=10-6,乘子向量的初始值λ0,罚参数向量的初始值σ0,罚参数的最大容许值σu=1010,罚参数的递增因子γ=10,不可行性的递减因子ζ=0.25,外层最大迭代次数Km=50.在内层迭代中,改进遗传算法的种群规模N=300,交叉概率Pc=0.7,变异概率Pm=0.1,内层最大迭代次数为5000。
将油水混合液分离水的价格pw固定为1$/RB,石油价格po以每次增加1美元的幅度由1$/RB提高到450$/RB,对每组(po,pw)及折旧率r分别取0、0.1、0.15和0.2所对应的优化问题进行求解,得到了不同价格和折旧率下各注水井注入率和各生产井产出率从第51~100个月的最优设定值.图11和图12给出了po=60$/RB,pw=1$/RB,折旧率r分别取0和0.2时,对应的5个注水井从第51到100个月的注入率的设定值最优方案.在相同价格和不同折旧率下,每个注水井从第51~100个月的注入率最优设定值是不同的.在两种不同折旧率的情况下,注入井I01基本上都是以其最大注入率注水;在折旧率r=0.0的情况下,注水井I02一直以最大注入率注水,注水井I03、I04和I05的注入率设定值开始有所变化,但稳定后也基本上以其最大注入率进行注水;在折旧率r=0.2的情况下,注水井I02和I03的注入率设定值有变化,注水井I04和I05的注入率最优设定值一直以其最大注水率进行注水.对油藏从第51~100个月的石油开采过程,各注水井注入率和各生产井产出率的不同设定值对应有不同的经济效益,利用问题(10)的目标函数进行评价,石油价格po=60$/RB和油水分离价格pw=1$/RB时得到的寻优结果如表3所示.
表3优化前后的经济收益比较
由表3的结果可以清晰地看出,对油藏生产过程中注水井注入率和生产井生产率的设定值优化带来了显著的经济效益,这充分说明了本发明提出的设定点优化模型和优化方法的有效性。
以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。
Claims (2)
1.一种基于混合遗传算法的油藏生产过程设定点优化方法,其特征在于按以下步骤进行:
步骤一、设定在油藏地面挖掘m个注水井和n个生产井,设定生产井j的总产出率为qj,m、n均为正整数;
建立油藏容量模型:其中ik是注水井k的注水率,τj为生产井j的时间常数,为注水井k与生产井j之间的权重,
设定生产井j在t时刻的石油产出率为qoj,建立油藏生产率模型求取参数aj和bj;
步骤二、建立油藏生产过程设定点优化模型:p0当前石油销售价格,pwj为从油水混合液中分离水的单位费用,qwj为生产井j在t时刻的水产出率;
步骤三、以油藏容量模型、油藏生产率模型为约束条件,以所述R为目标函数,得到油藏生产过程设定点优化问题:
步骤四、对于油藏生产过程设定点优化问题,将其离散化为一般的非线性约束优化问题:
步骤五、对z'(x),设定Lagrange乘子向量λk、和罚参数向量σk、增广Lagrange乘子法第k步所求得子问题为:
其中是修改增广Lagrange函数;
其中,
步骤六、由遗传算法对步骤五所述子问题进行全局寻优以得到下一迭代点,然后修正乘子向量和罚参数向量,检查收敛准则是否满足,当收敛准则不满足时重新构造界约束子问题,收敛准则满足时终止算法。
2.如权利要求1所述的基于混合遗传算法的油藏生产过程设定点优化方法,其特征是:步骤六所述遗传算法在进化过程中随机组合一个交叉算子和一个变异算子产生若干个新的个体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610226697.3A CN105913152B (zh) | 2016-04-13 | 2016-04-13 | 基于混合遗传算法的油藏生产过程设定点优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610226697.3A CN105913152B (zh) | 2016-04-13 | 2016-04-13 | 基于混合遗传算法的油藏生产过程设定点优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105913152A true CN105913152A (zh) | 2016-08-31 |
CN105913152B CN105913152B (zh) | 2019-09-24 |
Family
ID=56746885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610226697.3A Expired - Fee Related CN105913152B (zh) | 2016-04-13 | 2016-04-13 | 基于混合遗传算法的油藏生产过程设定点优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105913152B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101725346A (zh) * | 2009-12-15 | 2010-06-09 | 中国石油大学(华东) | 油藏井间动态连通性反演方法 |
CN102419786A (zh) * | 2011-10-13 | 2012-04-18 | 中国石油大学(华东) | 聚合物驱提高原油采收率的动态规划方法 |
CN104216341A (zh) * | 2013-05-31 | 2014-12-17 | 中国石油化工股份有限公司 | 一种基于改进随机扰动近似算法的油藏生产实时优化方法 |
CN104615862A (zh) * | 2015-01-14 | 2015-05-13 | 中国石油天然气股份有限公司 | 基于进化算法的高含水油田确定井位的方法 |
US20150160369A1 (en) * | 2013-12-09 | 2015-06-11 | King Fahd University Of Petroleum And Minerals | Method for well placement |
-
2016
- 2016-04-13 CN CN201610226697.3A patent/CN105913152B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101725346A (zh) * | 2009-12-15 | 2010-06-09 | 中国石油大学(华东) | 油藏井间动态连通性反演方法 |
CN102419786A (zh) * | 2011-10-13 | 2012-04-18 | 中国石油大学(华东) | 聚合物驱提高原油采收率的动态规划方法 |
CN104216341A (zh) * | 2013-05-31 | 2014-12-17 | 中国石油化工股份有限公司 | 一种基于改进随机扰动近似算法的油藏生产实时优化方法 |
US20150160369A1 (en) * | 2013-12-09 | 2015-06-11 | King Fahd University Of Petroleum And Minerals | Method for well placement |
CN104615862A (zh) * | 2015-01-14 | 2015-05-13 | 中国石油天然气股份有限公司 | 基于进化算法的高含水油田确定井位的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105913152B (zh) | 2019-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bao et al. | Fully implicit simulation of polymer flooding with MRST | |
Shirangi et al. | Closed-loop field development under uncertainty by use of optimization with sample validation | |
CN103422849B (zh) | 注水开发老油田井网重建方法及装置 | |
CN104750896B (zh) | 一种缝洞型碳酸盐岩油藏数值模拟方法 | |
CN105842753B (zh) | 对烃源岩内残余有机碳含量进行恢复的方法及装置 | |
CN104615862A (zh) | 基于进化算法的高含水油田确定井位的方法 | |
CN110032782A (zh) | 一种城市级智能交通信号控制系统及方法 | |
CN104868465B (zh) | 基于模糊机会约束的电力系统网架重构优化方法 | |
CN109543828A (zh) | 一种基于小样本条件下的吸水剖面预测方法 | |
CN106022463A (zh) | 基于改进粒子群算法的个性化学习路径优化方法 | |
CN103790565A (zh) | 复杂断块油藏开发井网优化设计方法 | |
CN111625922A (zh) | 一种基于机器学习代理模型的大规模油藏注采优化方法 | |
CN107239845A (zh) | 油藏开发效果预测模型的构建方法 | |
CN103984788A (zh) | 一种煤巷锚杆支护自动智能设计与优化系统 | |
Alenezi et al. | A data-driven smart proxy model for a comprehensive reservoir simulation | |
CN104216341A (zh) | 一种基于改进随机扰动近似算法的油藏生产实时优化方法 | |
CN103198354A (zh) | 一种利用bp神经网络与nsga2实现油田抽油机采油节能增产优化方法 | |
CN103886209A (zh) | 基于马尔科夫的跳变电力系统时滞稳定性分析系统及方法 | |
CN103177155A (zh) | 一种基于bp神经网络与spea2算法的油田抽油机采油节能增产优化方法 | |
CN106545319A (zh) | 断块油藏水驱后转人工气顶-边水双向驱提高采收率方法 | |
CN106505624A (zh) | 确定配电网分布式电源最优吸纳能力的调控系统及方法 | |
CN109800454A (zh) | 一种基于卡尔曼滤波的煤层气采集方法、系统及存储介质 | |
CN104466979B (zh) | 一种输电网无功优化调控系统及方法 | |
CN102352743A (zh) | 基于最优控制理论的聚合物驱方案优化设计方法 | |
CN109870909A (zh) | 一种基于rbf神经网络和自适应搜索的人工免疫算法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190924 Termination date: 20200413 |