CN110363325B - 复杂断块油藏多目标注采调控优化方法 - Google Patents
复杂断块油藏多目标注采调控优化方法 Download PDFInfo
- Publication number
- CN110363325B CN110363325B CN201910370374.5A CN201910370374A CN110363325B CN 110363325 B CN110363325 B CN 110363325B CN 201910370374 A CN201910370374 A CN 201910370374A CN 110363325 B CN110363325 B CN 110363325B
- Authority
- CN
- China
- Prior art keywords
- production
- injection
- scheme
- oil reservoir
- optimization
- 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.)
- Active
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 97
- 238000005457 optimization Methods 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 35
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 50
- 238000011161 development Methods 0.000 claims abstract description 48
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 7
- 230000008901 benefit Effects 0.000 claims abstract description 6
- 238000002347 injection Methods 0.000 claims description 25
- 239000007924 injection Substances 0.000 claims description 25
- 230000000694 effects Effects 0.000 claims description 17
- 239000003129 oil well Substances 0.000 claims description 12
- 238000004088 simulation Methods 0.000 claims description 12
- 238000006073 displacement reaction Methods 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 9
- 239000000243 solution Substances 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 238000011160 research Methods 0.000 claims description 6
- 230000035772 mutation Effects 0.000 claims description 5
- 230000003068 static effect Effects 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 230000000903 blocking effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims 1
- 238000009412 basement excavation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 11
- 238000000605 extraction Methods 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000004065 wastewater treatment Methods 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Animal Husbandry (AREA)
- General Health & Medical Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Educational Administration (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Primary Health Care (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种复杂断块油藏多目标注采调控优化方法,包括:步骤1,收集目标区块信息,对目标区块进行建模;步骤2,将目标区块按照特征及矛盾划分区域;步骤3,根据各个区块/区域的主要矛盾,对区块设计相应的目标函数;步骤4,利用多目标优化算法对目标函数进行多目标优化;步骤5,在获得的多目标优化结果集合中根据目标函数以及经济效益这些指标选取部分结果作为最终解决方案。该复杂断块油藏多目标注采调控优化方法充分利用现有井网,变“局部高效挖潜”为“分区高效调控”,分区域进行井网、注采调控优化,实现均衡水驱开发,从而解决目前技术中存在的无法改善平面矛盾的问题。
Description
技术领域
本发明涉及油气田开发技术领域,特别是涉及到一种复杂断块油藏多目标注采调控优化方法。
背景技术
胜利油区复杂断块油藏主力断块具有一定含油面积,能形成不规则注采井网,进入特高含水期后,由于断层遮挡作用、油藏非均质性强,层间、平面动用不均衡,剩余油分布差异大,特别是在低油价形势下,目前胜利油区复杂断块油藏已经经历多轮次的层系井网调整,缺乏进一步层系细分、井网加密调整的物质基础和经济可行性,油藏调整优化从层系细分重组、井网加密完善提高储量动用程度转变为充分利用目前层系、井网,注采液量调控优化努力实现均衡驱替开发。另外,长期以来,复杂断块油藏开发调整优化过程中主要是根据油藏工程评价认识设置3-4套优化方案,利用油藏数值模拟技术模拟计算确定最优方案,以层系、井网整体最优单一目标为主要的优化目标,没有考虑断块油藏不同构造部位地质、开发特点差异导致的开发不均衡性,由此会加剧驱替不均衡问题。同时人工优化方案设计优化方案数量有限,方案成败与否依赖油藏工程师的专业水平,很难做到科学决策。因此,针对复杂断块油藏面临的平面开发矛盾、目前低油价的形势,我们发明了一种新的复杂断块油藏多目标注采调控优化方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种实现了油藏不同区域均达到最优的目标效果的复杂断块油藏多目标注采调控优化方法。
本发明的目的可通过如下技术措施来实现:复杂断块油藏多目标注采调控优化方法,该复杂断块油藏多目标注采调控优化方法包括:步骤1,收集目标区块信息,对目标区块进行建模;步骤2,将目标区块按照特征及矛盾划分区域;步骤3,根据各个区块/区域的主要矛盾,对各个区块设计相应的目标函数;步骤4,利用多目标优化算法对目标函数进行多目标优化;步骤5,在所获得的多目标优化结果集合中,根据目标函数以及经济效益指标,选取部分结果作为最终解决方案。
本发明的目的还可通过如下技术措施来实现:
步骤1包括:
(1)根据断块油藏的特点,确定实施的断块油藏区块;
(2)分析油藏注采情况,确定区块适合进行注采调控区域;
(3)收集油藏静态数据;
(4)根据地质参数构建油藏静态模型;
(5)收集油藏动态生产数据;
(6)建立数值模拟模型,进行油藏数值模拟开发计算。
在步骤1中,利用精细地质研究成果对区块建模,之后基于油水井动态数据进行精细数值模拟研究,得到目标区块剩余油分布规律,对目标区块平面注采井生产情况、压力状况、水淹规律、剩余油分布情况进行分析。
在步骤2中,根据油藏的水驱动用情况以及剩余油分布状况,将油藏划分成不同的区域;复杂断块油藏划分为三个区域,分别是高部位断层遮挡难动用带、腰部井间动用不均衡带、低部位高水淹区。
在步骤3中,对于高部位断层遮挡难动用带,该区域油井含水率较低,水驱动用程度较低,在生产开发过程中先不关注产水情况,目标函数确保产油量最大化,见公式3:
在步骤3中,对于腰部井间动用不均衡带,该区域油井含水率中等,剩余储量较大,需均匀驱替以降低该区域内部矛盾,尽可能采出更多的剩余油;定义流场强度的均匀程度为流场变异系数,见式(5),流场变异系数越大代表油藏驱替越不均衡,系数越小代表油藏驱替程度越均匀;目标函数为流场变异系数最小,见式(4),约束条件为单井极限注采量与注采平衡:
在步骤3中,对于低部位高水淹带,该区域油井含水率较高,需要提高注采比,保持油藏整体压力以支撑其他区域的开采情况,保证水线均匀分布有助于抑制异常高低压区的产生,实现正常的水驱开发效率,约束条件为单井极限注采量与累计注采比2:1以上
步骤4包括:
(4.1)确定最大迭代次数maxInter、种群数Pop、交叉概率Pc和变异概率Pm;
(4.2)根据油藏对各个注采井的约束条件,随机生成Pop套生产方案;
(4.3)构建外部存档集,确保能够实时的保留优化过程中的最优结果,外部存档集的大小是Pop或者其他数值;
(4.4)将初始化的Pop套方案存入外部存档集;
(4.5)将初始方案中的各个方案随机选择进行交叉变异操作得到新的方案,将新的方案再次存入外部存档集,对存档集进行拥挤距离排序操作后再进行非支配排序操作;
(4.6)删除外部存档集中较差的部分,以使外部存档集内数据不超过外部存档集的大小;
(4.7)多次重复(4.5)(4.6)步,直到达到最大迭代次数maxInter 后停止并输出外部存档集。
在步骤4中,步骤4.5中的非支配排序和拥挤距离排序采用以下方法:
非支配等级:
设xa、xb均为油藏开发可行方案,当满足下列条件时,即称xa与xb相比是Pareto占优的;
其中f1~f3在本发明中是指三个区块的目标函数;
①将不被除了它之外任何一个开发方案支配的方案找出,认为它的支配等级为1;
②将所有的支配等级为1的方案从原方案集合中剔除,剩余的开发方案中不被除了它之外任何一个开发方案支配的方案设为2;
③持续进行以上步骤直到将所有的方案全部分配支配等级;
拥挤距离:
①将现存档集中的开发方案按照第一目标f1从小到大排序;
②将排序后的结果中第一个和最后一个拥挤距离设为无穷大,其余的按照以下公式计算:
di=abs(fi-1-fi+1)/abs(f1-fmax) (2)
其中i=2...max;abs代表绝对值;fi代表第i个目标;
③不断重复以上方法,将三个目标全部计算完,然后将三个目标的d 全部相加作为该开发方案的拥挤距离。
在步骤5中,根据油田实际生产需求从外部存档集中选择合适的解决方案;选取经济效益最大化结果,或选择某一个目标函数优化效果更好的优化效果,或者各个目标函数均衡优化的生产方案。
本发明中的复杂断块油藏多目标注采调控优化方法,涉及复杂断块油藏高含水期多目标注采调控优化,在给定初始条件下,结合断块油藏地质及开发特点,根据动用差异化程度及影响差异动用的主控因素不同,通过设置多目标优化算法优化目标函数,结合油藏数值模拟技术,在不增加额外投入的前提下,让油藏整体达到最优的注采优化目标效果。
本发明中的复杂断块油藏多目标注采调控优化方法,引入多目标优化机制,首先根据油藏地质、开发特点及开发面临的主要矛盾差异,按照优化目标差异将油藏分成不同差异动用区域,然后利用多目标优化算法对这些互相影响的区域进行优化。本方法充分利用现有井网,变“局部高效挖潜”为“分区高效调控”,分区域进行井网、注采调控优化,实现均衡水驱开发,从而解决目前技术中存在的无法改善平面矛盾的问题。本发明将油藏数值模拟方法与最优化理论方法结合,针对不同主控因素主导下的动用不均衡区设置目标函数及约束条件,自动搜索求解最优解,有效动用非主流线剩余油,实现了油藏不同区域均达到最优的目标效果。
附图说明
图1为本发明的复杂断块油藏多目标注采调控优化方法的一具体实施例的流程图;
图2为本发明的一具体实施例中目标断块油藏数值模拟含油饱和度图;
图3为本发明的一具体实施例中多目标优化流程图;
图4为本发明的一具体实施例中理论区块渗透率场图;
图5为本发明的一具体实施例中特征区块剖面图;
图6为本发明的一具体实施例中特征区块开始优化时含有饱和度的示意图;
图7为本发明的一具体实施例中特征区块分区示意图;
图8为本发明的一具体实施例中特征区块多目标优化前沿示意图;
图9为本发明的一具体实施例中特征区块优化前后累计产油量对比图;
图10为本发明的一具体实施例中特征区块优化前后含水率对比图;
图11为本发明的一具体实施例中特征区块优化前后含水率-累计产油量对比图;
图12为本发明的一具体实施例中特征区块优化前后注入量-产油量对比图;
图13为本发明的一具体实施例中特征区块优化前后压力对比图;
图14为本发明的一具体实施例中特征区块优化前后第一目标对比图;
图15为本发明的一具体实施例中特征区块优化前后第三目标对比图;
图16为本发明的一具体实施例中特征区块第二目标优化前效果图;
图17为本发明的一具体实施例中特征区块第二目标人工调节效果图;
图18为本发明的一具体实施例中特征区块第二目标优化后效果图;
图19为本发明的一具体实施例中特征区块优化前剩余油饱和度图;
图20为本发明的一具体实施例中特征区块人工调节剩余油饱和度图;
图21为本发明的一具体实施例中特征区块优化后剩余油饱和度图;
图22为本发明的一具体实施例中特征区块优化前流线场图;
图23为本发明的一具体实施例中特征区块人工调节流线场图;
图24为本发明的一具体实施例中特征区块优化后流线场图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的复杂断块油藏多目标注采调控优化方法的流程图。
在步骤101中,利用eclipse油藏模拟软件对目标区块进行数值模拟,其中图2为典型断块模型区块示意图。利用精细地质研究成果对区块建模,其中图4为所构建的理论区块示意图,根据断块油藏的特点,地层有一定的倾斜角,所以如图5所示,区块剖面图有一定的倾斜角。为了能够体现油藏开发后期的特点,让油藏生产一段时间,实现部分井含水率到达85%以上,其中剩余油饱和度如图6所示。之后基于油水井动态数据进行精细数值模拟研究,得到目标区块剩余油分布规律,对目标区块平面注采井生产情况、压力状况、水淹规律、剩余油分布情况进行分析,流程进入步骤 102。
在步骤102中,根据目标区块水淹及水驱动用情况,将油藏分为不同区域,本目标区块划分为高中低三个主要的部位,其中高部遮挡区平均含油饱和度为0.47,剩余可采储量为22.9万吨,占全部可采储量的28.2%;中部滞留区平均含油饱和度为0.4,剩余可采储量为36万吨,占全部可采储量的44.2%;低部水淹区平均含油饱和度为0.28,剩余可采储量为22.5 万吨,占全部可采储量的27.6%。不同区域含油饱和度情况如图7示,流程进入步骤103。
在步骤103中,根据步骤102所确定的划分区域部位,确定三个区域各自的目标函数。
(1)高部位断层遮挡难动用带,该区域油井含水率较低,水驱动用程度较低,在生产开发过程中先不关注产水情况,目标函数确保产油量最大化,见下式,约束条件为单井极限注采量与最小井底流压。
(2)腰部井间动用不均衡带,该区域油井含水率中等,剩余储量较大,需均匀驱替以降低该区域内部矛盾,尽可能采出更多的剩余油。定义流场强度的均匀程度为流场变异系数,见式(5),流场变异系数越大代表油藏驱替越不均衡,系数越小代表油藏驱替程度越均匀。目标函数为流场变异系数最小,见式(4),约束条件为单井极限注采量与注采平衡。
(3)低部位高水淹带,该区域油井含水率较高,需要提高注采比,保持油藏整体压力以支撑其他区域的开采情况,保证水线均匀分布有助于抑制异常高低压区的产生,实现正常的水驱开发效率,约束条件为单井极限注采量与累计注采比2:1以上,流程进入104。
以上公式符号含义分别为:
Ti为总的时间步,t为时间步长,Np为生产井数,NI为注入井数,为单位时间步第i口生产井产油量,为单位时间步第j口生产井注水量,rate为年利率,Wij为流场变异系数,Swcj和Swkj分别为生产前后的含水饱和度,为孔隙体积,ɑ与b分别为权重系数,本目标区块分别取值为1500与10。
在步骤104,利用多目标优化算法对断块油藏分区优化模型进行优化。具体步骤为:
(1)确定最大迭代次数为50次,交叉概率Pc为0.6,变异概率Pm为 0.1。随机生成50套开发方案,并标号分别为1~50号开发方案,利用油藏数值模拟软件对目标区块进行模拟,模拟完成后计算不同目标的数值,并存入存档集。
(2)将所有数据均转换为二进制数据,对于所有的第m个开发方案,随机选取第n个开发方案(n与m不相等),将两个方案进行交叉变异生成新的解决方案并存入外部存档集。
(3)将外部存档集按拥挤距离排序后按非支配等级排序,然后删除超出存档集数量的较差的解。
(4)往复迭代第(4.2)与(4.3)步到达最大迭代次数50代,输出外部存档集。将外部存档集中,NPV值最大的开发方案,三个目标分别最优的开发方案,以及三个目标效果均衡的开发方案取出作为最终方案的备选方案,流程进入步骤105。
在步骤105中,根据油田实际生产需求从外部存档集中选择合适的解决方案。可以选取经济效益最大化结果,也可以选择某一个目标函数优化效果更好的优化效果,或者各个目标函数均衡优化的生产方案。
具体实例:
1.首先根据断块油藏的特点建立特征模型。
为了有更强烈的平面矛盾,建立的渗透率场如图4所示,其中深色区域为高渗区域,浅色为低渗区域。根据断块油藏的特点,地层有一定的倾斜角,所以如图5所示,区块剖面图有一定的倾斜角。
为了能够体现油藏开发后期的特点,让油藏生产一段时间,实现部分井含水率到达85%以上,其中剩余油饱和度如图6所示。
2.根据油藏的平面矛盾,对油藏进行分区,按照上中下分为三个部分,三个部分如图7所示。
其中高部位为断层遮挡部位仅有一口水井含水率较低,需要提高采油量;腰部井间动用不均衡带井间弱驱部位有大量聚集于不同区域的剩余油,需要注采流线调整实现均衡水驱;低部位高水淹带为强驱部位含水率较高,需要人工边水驱提高驱油效率。
3.使用多目标算法对目标区块进行多目标优化,具体步骤如下:
(1)确定最大迭代次数为50次,交叉概率Pc为0.6,变异概率Pm为0.1。随机生成50套开发方案,并标号分别为1~50号开发方案,利用油藏数值模拟软件对目标区块进行模拟,模拟完成后计算不同目标的数值,并存入存档集。
(2)将所有数据均转换为二进制数据,rand是一个0到1之间的随机数,对于所有的第m个开发方案,随机选取第n个开发方案(n与m不相等),将两个方案进行交叉变异生成新的解决方案并存入外部存档集。
(3)将外部存档集按拥挤距离排序后按非支配等级排序,然后删除超出存档集数量的较差的解。
(4)往复迭代第(2)与(3)步到达最大迭代次数50代,输出外部存档集。得到的存档集在笛卡尔坐标系中如图8。将外部存档集中,NPV 值最大的开发方案,三个目标分别最优的开发方案,以及三个目标效果均衡的开发方案取出作为最终方案的备选方案。
其中Δt为时间步长,d为通货膨胀率,本次实验分别设为60天,通货膨胀设为0。中括号内为所有油井产收益减去污水处理花费再减去水井注入花费。变量含义分别为,a为油价,Qo为油井产油量,j为第j口油井,b 为污水处理单价,Qw为产水量,c为水价,Qwi为注入量。其中a为1500,b 为20,c为10。
(5)选取三个目标优化效果均衡的开发方案,该方案与人工调节方案以及初始方案对比结果如图9-13所示;优化前后不同目标的优化情况如图14-18所示,其中图16-18是优化前后中部区域流场强度对比图,优化后流场强度变得更加均匀;图19-24为优化前后剩余油与流线场对比图,优化过后生产效果有着显著改善。其中人工调节方案为提高INJ-02的注入量,降低PRO-05的采出量,提高PRO-02,PRO-04以及PRO-07的采出量。表1给出了三个方案的目标函数值,由图表可以看出多目标算法优化过后,不同目标均有着显著的优化,高部产油量明显提升,中部流场强度变得更加均匀,低部目标函数也有显著改善。
表1优化前后不同区域目标函数优化效果表
本发明中的特高含水期断块油藏多目标注采调控优化方法,首次提出了针对特高含水期断块油藏注采程度不均衡的差异化的多目标注采优化方法,根据油藏特征及剩余油富集模式差异及导致注采不均衡主控因素不同,合理设置分块分区、针对性建立注采调控目标函数,注采调控适配剩余油分布,实现均衡开发,进一步提高特高含水期断块油藏采收率工作,特别是针对已开发断块油田,具有一定的技术引领作用,为断块油藏高含水期的高效挖潜与提高采收率技术应用提供决策依据,其推广应用前景广阔,经济社会效益显著。
Claims (4)
1.复杂断块油藏多目标注采调控优化方法,其特征在于,该复杂断块油藏多目标注采调控优化方法包括:
步骤1,收集目标区域信息,对目标区域进行建模;
步骤2,将目标区域按照特征及矛盾划分区域;
步骤3,根据各个区域的主要矛盾,对区域设计相应的目标函数;
步骤4,利用多目标优化算法对目标函数进行多目标优化;
步骤5,在获得的多目标优化结果集合中选取经济效益最大化结果,或选择某一个目标函数优化效果更好的优化效果,或者各个目标函数均衡优化的生产方案作为最终解决方案;
在步骤2中,根据油藏的水驱动用情况以及剩余油分布状况,将油藏划分成不同的区域;复杂断块油藏划分为三个区域,分别是高部位断层遮挡难动用带、腰部井间动用不均衡带、低部位高水淹区;
在步骤3中,对于高部位断层遮挡难动用带,该区域油井含水率低,水驱动用程度低,在生产开发过程中先不关注产水情况,目标函数确保产油量最大化,见公式(3):
对于腰部井间动用不均衡带,该区域油井含水率中等,剩余储量大,需均匀驱替以降低该区域内部矛盾;定义流场强度的均匀程度为流场变异系数,见公式(5),流场变异系数越大代表油藏驱替越不均衡,系数越小代表油藏驱替程度越均匀;目标函数为流场变异系数最小,见公式(4),约束条件为单井极限注采量与注采平衡:
其中,f2为腰部位目标函数,m为网格数,k为网格编号,Wk为第k个网格的流场变异系数,ΔRwk是流过网格k过水倍数,Swck和Swk分别为网格k生产前后的含水饱和度,QINk是流过网格k的流体流量,Vk为网格k的体积,φk为网格k的孔隙度;
对于低部位高水淹区,该区域油井含水率高,需要提高注采比,保持油藏整体压力以支撑其他区域的开采情况,保证水线均匀分布有助于抑制异常高低压区的产生,实现正常的水驱开发效率,约束条件为单井极限注采量与累计注采比2:1以上
步骤4包括:
(4.1)确定最大迭代次数maxInter、种群数Pop、交叉概率Pc和变异概率Pm;
(4.2)根据油藏对各个注采井的约束条件,随机生成Pop套生产方案;
(4.3)构建外部存档集,确保能够实时的保留优化过程中的最优结果,外部存档集的大小是Pop或者其他数值;
(4.4)将初始化的Pop套方案存入外部存档集;
(4.5)将初始方案中的各个方案随机选择进行交叉变异操作得到新的方案,将新的方案再次存入外部存档集,对存档集进行拥挤距离排序操作后再进行非支配排序操作,删除外部存档集中排序最末的部分,以使外部存档集内数据不超过外部存档集的大小;
(4.6)多次重复(4.5)步,直到达到最大迭代次数maxInter后停止并输出外部存档集。
2.根据权利要求1所述的复杂断块油藏多目标注采调控优化方法,其特征在于,步骤1包括:
(1)根据断块油藏的特点,确定实施的断块油藏区域;
(2)分析油藏注采情况,确定区域适合进行注采调控区域;
(3)收集油藏静态数据;
(4)根据地质参数构建油藏静态模型;
(5)收集油藏动态生产数据;
(6)建立数值模拟模型,进行油藏数值模拟开发计算。
3.根据权利要求1所述的复杂断块油藏多目标注采调控优化方法,其特征在于,在步骤1中,利用精细地质研究成果对区域建模,之后基于油水井动态数据进行精细数值模拟研究,得到目标区域剩余油分布规律,对目标区域平面注采井生产情况、压力状况、水淹规律、剩余油分布情况进行分析。
4.根据权利要求1所述的复杂断块油藏多目标注采调控优化方法,其特征在于,在步骤4中,步骤4.5中的非支配排序和拥挤距离排序采用以下方法:
非支配等级:
设xa、xb为两种开发方案,fu(xa)是xa的第u个目标函数,fu(xb)是xb的第u个目标函数,fv(xa)是xa的第v个目标函数,fv(xb)是xb的第v个目标函数,当满足下列条件时,即称xa与xb相比是Pareto占优的;
其中f1~f3是指三个区域的目标函数;
①将不被除了它之外任何一个开发方案支配的方案找出,认为它的支配等级为1;
②将所有的支配等级为1的方案从原方案集合中剔除,剩余的开发方案中不被除了它之外任何一个开发方案支配的方案设为2;
③持续进行以上步骤直到将所有的方案全部分配支配等级;
拥挤距离:
①将现存档集中的开发方案按照第一目标f1从小到大排序;
②将排序后的结果中第一个和最后一个拥挤距离设为无穷大,其余的按照以下公式计算:
dr=abs(f1,r-1-f1,r+1)/abs(f1,r-f1,max) (2)
其中r为外部存档集中的排序后的第r个方案,最后一个设为max,r=2…max;abs代表绝对值;f1,r代表第一目标f1第r个方案结果;
③不断重复以上方法,将三个目标全部计算完,然后将三个目标的拥挤距离全部相加作为该开发方案的拥挤距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910370374.5A CN110363325B (zh) | 2019-05-06 | 2019-05-06 | 复杂断块油藏多目标注采调控优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910370374.5A CN110363325B (zh) | 2019-05-06 | 2019-05-06 | 复杂断块油藏多目标注采调控优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110363325A CN110363325A (zh) | 2019-10-22 |
CN110363325B true CN110363325B (zh) | 2022-04-19 |
Family
ID=68215261
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910370374.5A Active CN110363325B (zh) | 2019-05-06 | 2019-05-06 | 复杂断块油藏多目标注采调控优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110363325B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110795893A (zh) * | 2019-11-07 | 2020-02-14 | 中国石油化工股份有限公司 | 注水开发油田注采系统能耗整体优化方法 |
CN111364955A (zh) * | 2020-03-17 | 2020-07-03 | 中国石油化工股份有限公司 | 注采井间流场演化模拟方法 |
CN111173507B (zh) * | 2020-03-17 | 2023-04-07 | 中国石油化工股份有限公司 | 高含水油田剩余油预测方法 |
CN111861129B (zh) * | 2020-06-22 | 2022-07-01 | 中国石油大学(华东) | 一种基于多尺度油藏模型的多保真注采优化方法 |
CN113969768B (zh) * | 2020-07-23 | 2024-05-31 | 中国石油化工股份有限公司 | 一注多采井组定向赋能-差异释放式体积水驱方法 |
CN114165227B (zh) * | 2020-09-11 | 2024-03-01 | 中国石油化工股份有限公司 | 基于油藏均衡驱替的流场调整可行性定量评价方法 |
CN113240161B (zh) * | 2021-04-27 | 2024-06-14 | 中国石油天然气股份有限公司 | 净现值预测模型建立方法、装置、存储介质及电子设备 |
CN115478821B (zh) * | 2022-07-27 | 2024-07-12 | 中国石油化工股份有限公司 | 一种边水油藏立体抑水的方法 |
CN116882323B (zh) * | 2023-09-07 | 2023-11-28 | 中国石油大学(华东) | 一种考虑时序性及细分任务的自适应代理策略优化方法 |
CN117076956B (zh) * | 2023-10-16 | 2024-01-26 | 西安石油大学 | 缝洞油藏物理模型相似准数优化方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101657890B1 (ko) * | 2015-04-06 | 2016-09-20 | 서울대학교산학협력단 | 다목적 유전 알고리즘과 실물옵션에 기반한 저류층의 생산성 평가 방법 |
CN109165783A (zh) * | 2018-08-15 | 2019-01-08 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 油田开发方案优化方法及装置 |
CN109281648A (zh) * | 2018-09-27 | 2019-01-29 | 中国石油天然气股份有限公司 | 确定油藏的合理井网密度的方法和装置 |
CN109577968A (zh) * | 2018-12-03 | 2019-04-05 | 中国石油大学(华东) | 一种水平井流入控制装置限流强度多目标优化方法 |
-
2019
- 2019-05-06 CN CN201910370374.5A patent/CN110363325B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101657890B1 (ko) * | 2015-04-06 | 2016-09-20 | 서울대학교산학협력단 | 다목적 유전 알고리즘과 실물옵션에 기반한 저류층의 생산성 평가 방법 |
CN109165783A (zh) * | 2018-08-15 | 2019-01-08 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 油田开发方案优化方法及装置 |
CN109281648A (zh) * | 2018-09-27 | 2019-01-29 | 中国石油天然气股份有限公司 | 确定油藏的合理井网密度的方法和装置 |
CN109577968A (zh) * | 2018-12-03 | 2019-04-05 | 中国石油大学(华东) | 一种水平井流入控制装置限流强度多目标优化方法 |
Non-Patent Citations (1)
Title |
---|
"断块油藏高部位剩余油径向钻孔参数优化";杨志军等;《特种油气藏》;20180228;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110363325A (zh) | 2019-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110363325B (zh) | 复杂断块油藏多目标注采调控优化方法 | |
CN111441767B (zh) | 油藏生产动态预测方法及装置 | |
US8005658B2 (en) | Automated field development planning of well and drainage locations | |
CN108843286B (zh) | 一种分层采油选井的技术方法 | |
CN103617147A (zh) | 一种矿井突水水源层识别方法 | |
CN108416475A (zh) | 一种页岩气产能非确定性预测方法 | |
CN109209312B (zh) | 一种适用聚合物驱的资源开发方法 | |
CN109063403B (zh) | 一种滑溜水压裂优化设计方法 | |
CN106503792A (zh) | 一种基于自适应模块化神经网络的瓦斯浓度预测方法 | |
CN108280289B (zh) | 基于局部加权c4.5算法的冲击地压危险等级预测方法 | |
CN115438823A (zh) | 一种井壁失稳机制分析与预测方法及系统 | |
CN111350485B (zh) | 井网调整方法及装置 | |
Ogbeiwi et al. | Robust optimisation of water flooding using an experimental design-based surrogate model: A case study of a Niger-Delta oil reservoir | |
CN117291125B (zh) | 一种水驱油藏井网与射孔层段综合调整优化方法及系统 | |
CN109555518B (zh) | 一种基于聚类和判别算法的冲积扇单井构型识别方法 | |
CN117574755B (zh) | 页岩储层井工厂水平井压裂参数分层多级优化方法 | |
CN113431534A (zh) | 一种低渗致密油藏co2吞吐选井方法 | |
CN110924935B (zh) | 致密油藏井底流压调控方案确定方法、装置和设备 | |
Saafan et al. | A multi-objective optimization approach to groundwater management using genetic algorithm | |
CN113610446B (zh) | 一种复杂分散断块油田群投产顺序的决策方法 | |
CN112348350B (zh) | 一种油田开发后期层系重组方法、计算机设备及存储介质 | |
Al-Shamma et al. | History matching using hybrid parameterisation and optimisation methods | |
CN115705452A (zh) | 整装砂岩油藏开发中后期的新型采收率预测方法 | |
CN113837890B (zh) | 一种基于数据统计的页岩气水平井井型优选方法 | |
Al-Aghbari et al. | Hybrid approach of using bi-objective genetic programming in well control optimization of waterflood management |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |