CN117371361B - 一种致密油气藏压裂水平井变制度生产动态模拟监测方法 - Google Patents
一种致密油气藏压裂水平井变制度生产动态模拟监测方法 Download PDFInfo
- Publication number
- CN117371361B CN117371361B CN202311675971.1A CN202311675971A CN117371361B CN 117371361 B CN117371361 B CN 117371361B CN 202311675971 A CN202311675971 A CN 202311675971A CN 117371361 B CN117371361 B CN 117371361B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- production
- tight
- dimensionless
- horizontal 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.)
- Active
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 229
- 230000009466 transformation Effects 0.000 title claims abstract description 29
- 238000012544 monitoring process Methods 0.000 title claims abstract description 22
- 238000005094 computer simulation Methods 0.000 title claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 175
- 239000007789 gas Substances 0.000 claims abstract description 136
- 239000004215 Carbon black (E152) Substances 0.000 claims abstract description 47
- 229930195733 hydrocarbon Natural products 0.000 claims abstract description 47
- 150000002430 hydrocarbons Chemical class 0.000 claims abstract description 47
- 238000004088 simulation Methods 0.000 claims abstract description 35
- 238000006243 chemical reaction Methods 0.000 claims abstract description 34
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 24
- 239000003345 natural gas Substances 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 150000001875 compounds Chemical class 0.000 claims abstract description 4
- 230000010354 integration Effects 0.000 claims abstract description 3
- 230000008569 process Effects 0.000 claims description 66
- 230000035699 permeability Effects 0.000 claims description 40
- 238000009792 diffusion process Methods 0.000 claims description 16
- 230000001186 cumulative effect Effects 0.000 claims description 14
- 238000012986 modification Methods 0.000 claims description 14
- 230000004048 modification Effects 0.000 claims description 14
- 239000002131 composite material Substances 0.000 claims description 10
- 230000006835 compression Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 8
- 230000000704 physical effect Effects 0.000 claims description 8
- 238000002407 reforming Methods 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 7
- 238000012795 verification Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 abstract description 9
- 239000003209 petroleum derivative Substances 0.000 abstract description 2
- 238000005206 flow analysis Methods 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 8
- 238000011161 development Methods 0.000 description 7
- 238000005192 partition Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
- E21B2200/20—Computer models or simulations, e.g. for reservoirs under production, drill bits
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种致密油气藏压裂水平井变制度生产动态模拟监测方法,具体涉及石油天然气开采技术领域。本发明基于致密油气藏的现场数据,划分生产阶段后,利用致密气藏五线性流模型进行反演模拟获取未知参数组中各参数的最优值后,再建立生产制度转换模型,结合致密油气藏的现场数据和未知参数组反演结果,利用致密气藏五线性流模型模拟日产量并进行拟合,得到压裂水平井日产量预测模型并验证其准确性后,通过时间积分得到压裂水平井累计产量预测模型,将压裂水平井日产量预测模型和压裂水平井累计产量预测模型用于产能预测。本发明方法有效克服了传统复合线性流解析方法模拟监测精度不足的问题,计算效率高,为后续生产制度优化调整提供了依据。
Description
技术领域
本发明涉及石油天然气开采技术领域,具体涉及一种致密油气藏压裂水平井变制度生产动态模拟监测方法。
背景技术
致密油气藏具有储层物性差且自然产能低的特点,必须通过水平井多级压裂才能实现经济开采。压后储层分为压裂裂缝、改造储层和未改造储层等多个连续流动区域,非均质特征显著。目前,主要通过简化物理模型与求解过程建立复合分区模型,利用解析方法或半解析方法描述压裂水平井的性能,研究非常规油气藏的开发动态。
虽然,解析方法或半解析方法的计算效率高,可用于储层或裂缝等未知参数的快速反演,科学认识油气藏生产动态。但是,现有基于解析方法或半解析方法所构建的模型主要针对单一生产制度条件(例如,定产量或定井底流压),但是,在实际生产过程中,特别是油藏开发初期,工作制度经常会发生调整,利用现有模型模拟预测变生产制度条件下压裂水平井的生产动态精度不足,严重制约了后续对生产制度优化调整的理论指导。因此,亟需提供一种致密油气藏压裂水平井变制度生产动态模拟监测方法,为致密油气藏压裂水平井的生产制度优化调整提供理论指导。
发明内容
本发明旨在解决现有技术难以准确预测致密油气藏压裂水平井变生产制度实际数据的问题,提出了一种致密油气藏压裂水平井变制度生产动态模拟监测方法,为指导致密油气藏压裂水平井生产制度的优化调整奠定了基础。
为实现上述目的,本发明采用如下技术方案:
一种致密油气藏压裂水平井变制度生产动态模拟监测方法,具体包括以下步骤:
步骤1,选取致密油气藏,获取致密油气藏的储层物性参数、压裂作业参数和生产历史,进行数据整理准备;
步骤2,根据生产历史中的井底流压和累计产量,将致密油气藏的生产历史划分为多个生产阶段;
步骤3,根据压裂水平井的特征,选取致密气藏五线性流模型作为复合线性流模型,结合压裂水平井的邻井资料和现场经验,利用致密气藏五线性流模型进行反演模拟获取由裂缝渗透率、裂缝半长、储层改造区域宽度、储层改造区域渗透率所组成的未知参数组,得到未知参数组中各参数的最优值;
步骤4,将致密油气藏的各生产阶段分别划分为两个部分,建立生产制度转换模型,结合步骤1中所获取的致密油气藏储层物性参数、压裂作业参数、生产历史和步骤3中所确定的未知参数组反演结果,利用致密气藏五线性流模型模拟得到各生产阶段中不同生产部分的日产量模拟值并进行数据拟合,得到致密油气藏的压裂水平井日产量预测模型,验证压裂水平井日产量预测模型的预测效果后,对验证后的压裂水平井日产量预测模型在时间上进行积分,得到致密油气藏的压裂水平井累计产量预测模型;
步骤5,利用验证后致密油气藏的压裂水平井日产量预测模型和压裂水平井累计产量预测模型预测变制度条件下致密油气藏内压裂水平井的生产动态,确定致密油气藏中压裂水平井的日产量预测值和累计产量预测值。
优选地,所述步骤1中,所述致密油气藏的储层物性参数包括储层的厚度、初始压力、温度、储层渗透率、孔隙度、饱和度以及储层内的油气密度;
所述致密油气藏的压裂作业参数包括压裂水平井的长度、裂缝簇数和井间距;
所述致密油气藏的生产历史包括生产时间、井底流压、产出速度、日产量和累计产量。
优选地,所述步骤2中,同一生产阶段内压裂水平井的生产制度不变且压裂水平井的日产量和井底流压恒定。
优选地,所述步骤3中,具体包括以下步骤:
步骤3.1,选取致密气藏五线性流模型作为复合线性流模型;
所述致密气藏五线性流模型包括第一区域、第二区域、第三区域、第四区域和第五区域,其中,第一区域为储层改造区域,第二区域、第三区域和第四区域均为储层未改造区域,第五区域为水力压裂区域,内部存在水力裂缝;
步骤3.2,建立天然气拟压力函数,将无因次变量引入致密气藏五线性流模型中,得到井底流压和产出量的计算公式;
步骤3.3,设置未知参数组中各参数的取值范围和最大反演次数,利用致密气藏五线性流模型进行反演,得到未知参数组中各参数的最优值。
优选地,所述天然气拟压力函数为:
(1)
式中,为储层压力,单位为Pa;为天然气压缩因子;为粘度,单位为cp;
所述无因次变量包括无因次拟压力、无因次时间、无因次产量、无因次扩散系数和无因次距离;
所述无因次拟压力为:
(2)
式中,为第一区域的渗透率,单位为m2;为储层厚度,单位为m;为标况温度,
取值为273.15K;为储层的初始压力,单位为Pa;为标况压力,取值为1.01325×105Pa;为地面产量,单位为m3/s;为储层温度,单位为K;
所述无因次时间为:
(3)
式中,为时间,单位为s;为第一区域的孔隙度;为第一区域的综合压缩系
数,单位为Pa-1;为裂缝长度,单位为m;
所述无因次产量为:
(4)
式中,为井底流压,单位为Pa;
所述无因次扩散系数为:
(5)
其中,
(6)
式中,为区域编号,当时,表示第一区域,当时,表示第二区域,当时,表示第三区域,当时,表示第四区域,当时,表示第五区域;为第
区域的扩散系数;为第区域的渗透率,单位为m2;为第区域的孔隙度;为第区域
的综合压缩系数,单位为Pa-1。
优选地,所述无因次距离包括无因次宽度、储层改造区域的无因次宽度、裂缝半间距的无因次宽度、无因次长度、储层改造区域的无因次长度、无因次井间半间距长度和无因次裂缝导流能力,其中,所述无因次宽度为:
(7)
式中,为无因次宽度;为距离压裂裂缝的距离,单位为m;
所述储层改造区域的无因次宽度为:
(8)
式中,为储层改造区域的无因次宽度;为储层改造区域的宽度,单位为m;
所述裂缝半间距的无因次宽度为:
(9)
式中,为裂缝半间距无因次宽度;为裂缝半间距,单位为m;
所述无因次长度为:
(10)
式中,为无因次长度;为距离水平井筒的距离,单位为m;
所述储层改造区域的无因次长度为:
(11)
式中,为储层改造区域的无因次长度;为储层改造区域的长度,单位为m;
所述无因次井间半间距长度为:
(12)
式中,为无因次井间半间距长度;为井间半间距的长度,单位为m;
所述无因次裂缝导流能力为:
(13)
式中,为无因次裂缝导流能力;为裂缝渗透率,单位为m2;为裂缝宽度,单
位为m。
优选地,将无因次变量引入致密气藏五线性流模型中,所述第三区域和第四区域的线性流模型均为:
(14)
其中,
(15)
(16)
式中,为拉普拉斯空间第区域的无因次拟压力;为拉普拉斯变换参数;为拉普拉斯空间第区域的无因次拟压力;为双曲余弦函数;为第区域复合变量;
所述第二区域的线性流模型为:
(17)
其中,
(18)
(19)
式中,为拉普拉斯空间第一区域的无因次拟压力,为拉普拉斯空间第二
区域的无因次拟压力,为拉普拉斯空间第四区域的无因次拟压力;为第二区域渗透
率;为第四区域渗透率;为第二区域扩散系数;为拉普拉斯参数的相关函数;为双曲正切函数;为第四区域复合变量;
所述第一区域的线性流模型为:
(20)
其中,
(21)
(22)
(23)
(24)
式中,为拉普拉斯空间第三区域的无因次拟压力,为拉普拉斯空间水力
压裂区域的无因次拟压力;为第一区域扩散系数;为无因次裂缝宽度;为第三区
域渗透率;、、均为拉普拉斯参数的相关函数;为双曲正弦函数;
为第三区域复合变量。
优选地,所述水力压裂区域的线性流模型为:
(25)
其中,
(26)
(27)
式中,、均为拉普拉斯参数的相关函数;
根据水力压裂区域的线性流模型,得到定产量条件下的井底流压,如公式(28)所示:
(28)
式中,为定产量条件下的井底流压,无因次;
基于Duhamel原理,根据定产量条件下的井底流压,得到定井底流压条件下的
产出量,如公式(29)所示:
(29)
式中,为定压力条件下的产出量,无因次。
优选地,所述步骤3.3中,具体包括以下步骤:
步骤3.3.1,根据压裂水平井的邻井资料和现场经验,分别设置未知参数组中各参数的取值范围,并预设最大反演次数和未知参数组中各参数的初始值;
步骤3.3.2,将步骤1中所获取致密油气藏的储层物性参数、压裂作业参数和生产历史和预设的未知参数组的初始值输入至致密气藏五线性流模型中,利用致密气藏五线性流模型模拟致密油气藏的生产过程得到日产量的模拟值以及井底流压的模拟值,并根据生产制度决定计算致密油气藏日产量或致密油气藏井底流压的目标函数值;
步骤3.3.3,判断当前反演次数是否超过预设最大反演次数,若当前迭代次数未超过预设最大反演次数,则更新反演次数和未知参数组的取值,返回步骤3.3.2中继续利用致密气藏五线性流模型进行模拟,否则,则进入步骤3.3.4中;
步骤3.3.4,获取致密气藏五线性流模型各次反演模拟得到的致密油气藏日产量目标函数值,确定致密油气藏日产量目标函数值的最小值所对应未知参数组中各参数的取值,得到未知参数组的反演结果,包括裂缝渗透率最优值、裂缝半长最优值、储层改造区域宽度最优值、储层改造区域渗透率最优值。
优选地,所述致密油气藏日产量的目标函数,如公式(30)所示:
(30)
式中,为致密油气藏日产量的目标函数;为最大模拟时间,单位为s;为致密油气藏日产量的模拟值,单位为m3/s;为致密油气藏日产量的真实值,单
位为m3/s,根据致密油气藏的历史生产历史确定;为致密油气藏日产量目标函数的比率系
数;
所述致密油气藏井底流压的目标函数,如公式(31)所示:
(31)
其中,
(32)
(33)
(34)
(35)
式中,为致密油气藏井底流压的目标函数;为致密油气藏产量规整
化压力的模拟值,为致密油气藏模拟得到的产量规整化压力的导数;为
致密油气藏产量规整化压力的真实值,为致密油气藏产量规整化压力导数的真实
值;为致密油气藏井底流压目标函数的比率系数;为产量规整化压力的权重系数,
为产量规整化压力导数的权重系数。
优选地,所述步骤4中,具体包括以下步骤:
步骤4.1,将致密油气藏的各生产阶段均划分为第一生产部分和第二生产部分,建立生产制度转换模型;
步骤4.2,根据步骤1中所获取致密油气藏的储层物性参数、压裂作业参数、生产历史以及步骤3中所确定未知参数组的最优值重新设置致密气藏五线性流模型,利用致密气藏五线性流模型结合生产制度转换模型模拟,得到各生产阶段中第一生产部分和第二生产部分的日产量模拟值并进行对数曲线拟合,得到致密油气藏的压裂水平井日产量预测模型;
步骤4.3,基于致密油气藏的生产历史和生产制度,验证日产量预测模型的预测效果;
通过计算压裂水平井日产量模拟值与压裂水平井日产量实测值之间的精度值,若精度值小于预设精度值,则调整生产制度转换模型的制度转换参数,更新压裂水平井日产量预测模型后,返回步骤4.2中,否则,结束验证并进入步骤4.4;
步骤4.4,获取验证后的压裂水平井日产量预测模型,并进行时间上积分得到压裂水平井累计产量预测模型。
优选地,所述生产制度转换模型,如公式(36)所示:
(36)
式中,为第个生产阶段中第一生产部分生产制度的初始化压力,单位为Pa;为第个生产阶段中第一生产部分生产制度的初始化压力,单位为Pa;为第个生产阶段的平均压力,单位为Pa;为第个生产阶段中第二生产部分生产制度的初
始化压力,单位为Pa;为第一生产部分和第二生产部分之间生产制度的制度转换时间,单
位为s;、、均为制度转换参数。
优选地,所述步骤5中,通过改变生产制度,利用生产制度转换模型确定生产制度中各生产阶段之间的制度转换时间,结合致密气藏五线性流模型模拟得到致密油气藏内压裂水平井在变制度条件下的生产动态,利用验证后的压裂水平井日产量预测模型和压裂水平井累计产量预测模型进行预测,得到不同生产制度条件下致密油气藏中压裂水平井的日产量预测值和累计产量预测值。
本发明所带来的有益技术效果:
本发明提出了一种致密油气藏压裂水平井变制度生产动态模拟监测方法,通过将解析模型与智能优化方法相结合,快速精确地预测致密油气藏中压裂水平井的生产动态,既克服了现有解析模型应用于致密油气藏压裂水平井变生产制度问题的局限性,又避免了智能优化方法应用于致密油气藏参数优化时缺少物理背景的问题。
本发明通过构建生产制度转换模型,引入制度转换参数,结合不同初始化压力条件下的日产量,利用复合线性流模型模拟得到变制度条件下致密油气藏中压裂水平井的日产量、累计产量、井底流压等生产动态,建立致密油气藏的压裂水平井日产量预测模型和压裂水平井累计产量预测模型并进行验证,实现了对连续变化生产制度条件下致密油气藏压裂水平井油气产量的准确预测。
同时,本发明方法还具有强扩展性,适用于致密油气藏各类复合线性流分区模型中,有效克服了传统解析模型的局限性,能够快速模拟不同生产制度条件下压裂水平井的开发动态,准确预测压裂水平井的产量,为指导致密油气藏水平井的生产制度优化奠定了基础。
附图说明
图1为本发明致密油气藏压裂水平井变制度生产动态模拟监测方法的流程图。
图2为致密气藏五线性流模型的示意图,图中,为储层改造区域的宽度,为裂
缝半间距,为储层改造区域的长度,为井间半间距的长度,为裂缝长度。
图3为同一生产阶段中第一生产部分和第二生产部分的压力剖面图。
图4为对数曲线拟合所确定的致密油气藏压裂水平井日产量预测模型。
图5为变制度条件下采用本发明预测结果与实际产量数据的对比图。
具体实施方式
下面结合附图和某致密气藏中的压裂水平井为例,对本发明作进一步详细描述。
由于现有的压裂水平井复合分区模型应用于致密油气藏压裂水平井变制度生产动态预测时,预测过程中假设生产制度不变(即定井底流压或定产量),未考虑实际生产过程中变制度产生对致密油气藏压裂水平井日产量和累计产量的影响,存在压裂水平井产能预测准确性较差的问题,所以本发明公开了一种致密油气藏压裂水平井变制度生产动态模拟监测方法。
本实施例采用本发明所提出的致密油气藏压裂水平井变制度生产动态模拟监测方法,如图1所示,具体包括以下步骤:
步骤1,选取致密油气藏,获取致密油气藏的储层物性参数、压裂作业参数和生产历史,进行数据整理准备,其中,所述致密油气藏的储层物性参数包括储层的厚度、初始压力、温度、储层渗透率、孔隙度、饱和度以及储层内的油气密度;所述致密油气藏的压裂作业参数包括压裂水平井的长度、裂缝簇数和井间距;所述致密油气藏的生产历史包括生产时间、井底流压、产出速度、日产量和累计产量。
本实施例中,致密油气藏具体为鄂尔多斯盆地内的某致密气藏,获取致密气藏的厚度为15m,初始压力为18 Mpa,温度为67℃,孔隙度为3.5%,含气饱和度为60%,气体比重为0.7,渗透率为0.0005 mD,裂缝半间距为33.33m,井间半间距为150m,致密气藏的中部深度为2067m。
步骤2,根据生产历史中的井底流压和累计产量,将致密油气藏的生产历史划分为11个生产阶段,同一生产阶段内压裂水平井的生产制度不变且压裂水平井的日产量和井底流压恒定。
步骤3,根据压裂水平井的特征,选取致密气藏五线性流模型作为复合线性流模型,结合压裂水平井的邻井资料和现场经验,利用复合线性流模型进行反演模拟获取由裂缝渗透率、裂缝半长、储层改造区域宽度、储层改造区域渗透率所组成的未知参数组,反演模拟得到未知参数组中各参数的最优值,具体包括以下步骤:
步骤3.1,选取致密气藏五线性流模型作为复合线性流模型,致密气藏五线性流模型包括第一区域、第二区域、第三区域、第四区域和第五区域,如图2所示,其中,第一区域为储层改造区域,第二区域、第三区域和第四区域均为储层未改造区域,第五区域为水力压裂区域,内部存在水力裂缝。
步骤3.2,由于气体物性参数具有强非线性,为了方便求解裂缝渗透率、裂缝半长、储层改造区域宽度和储层改造区域渗透率,建立天然气拟压力函数,将无因次变量引入致密气藏五线性流模型中,得到井底流压和产出量的计算公式。
本实施例中,所述天然气拟压力函数如公式(1)所示:
(1)
式中,为天然气拟压力函数,单位为Pa2/cp;为储层压力,单位为Pa;为天然
气压缩因子;为粘度,单位为cp。
所述无因次变量包括无因次拟压力、无因次时间、无因次产量、无因次扩散系数和无因次距离。
所述无因次拟压力的计算公式,如公式(2)所示:
(2)
式中,为无因次压力;为第一区域的渗透率,单位为m2;为储层厚度,单位
为m;为标况温度,取值为273.15K;为储层的初始压力,单位为Pa;为标况压力,取
值为1.01325×105Pa;为地面产量,单位为m3/s;为储层温度,单位为K。
所述无因次时间的计算公式,如公式(3)所示:
(3)
式中,为无因次时间;为时间,单位为s;为第一区域的孔隙度;为第一区
域的综合压缩系数,单位为Pa-1;为裂缝长度,单位为m。
所述无因次产量的计算公式,如公式(4)所示:
(4)
式中,为无因次产量;为井底流压,单位为Pa。
所述无因次扩散系数的计算公式,如公式(5)所示:
(5)
其中,
(6)
式中,为区域编号,当时,表示第一区域,当时,表示第二区域,当时,表示第三区域,当时,表示第四区域,当时,表示第五区域;为第
区域的无因次扩散系数;为第区域的扩散系数;为第区域的渗透率,单位为m2;为
第区域的孔隙度;为第区域的综合压缩系数,单位为Pa-1。
所述无因次距离包括无因次宽度、储层改造区域的无因次宽度、裂缝半间距的无因次宽度、无因次长度、储层改造区域的无因次长度、无因次井间半间距长度和无因次裂缝导流能力,其中,所述无因次宽度如公式(7)所示:
(7)
式中,为无因次宽度;为距离压裂裂缝的距离,单位为m。
所述储层改造区域的无因次宽度如公式(8)所示:
(8)
式中,为储层改造区域的无因次宽度;为储层改造区域的宽度,单位为m。
所述裂缝半间距的无因次宽度如公式(9)所示:
(9)
式中,为裂缝半间距无因次宽度;为裂缝半间距,单位为m。
所述无因次长度如公式(10)所示:
(10)
式中,为无因次长度;为距离水平井筒的距离,单位为m。
所述储层改造区域的无因次长度如公式(11)所示:
(11)
式中,为储层改造区域的无因次长度;为储层改造区域的长度,单位为m。
所述无因次井间半间距长度如公式(12)所示:
(12)
式中,为无因次井间半间距长度;为井间半间距的长度,单位为m。
所述无因次裂缝导流能力如公式(13)所示:
(13)
式中,为无因次裂缝导流能力;为裂缝渗透率,单位为m2;为裂缝宽度,单
位为m。
将无因次变量引入致密气藏五线性流模型中,确定致密气藏五线性流模型中各区域的线性流模型,得到致密气藏五线性流模型中各区域的无因次拟压力解。
所述致密气藏五线性流模型中第三区域和第四区域的线性流模型均设置为:
(14)
其中,
(15)
(16)
式中,为拉普拉斯空间第区域的无因次拟压力;为拉普拉斯变换参数;为拉普拉斯空间第区域的无因次拟压力;为双曲余弦函数;为第区域复合变量。
所述致密气藏五线性流模型中第二区域的线性流模型设置为:
(17)
其中,
(18)
(19)
式中,为拉普拉斯空间第一区域的无因次拟压力,为拉普拉斯空间第二
区域的无因次拟压力,为拉普拉斯空间第四区域的无因次拟压力;为第二区域渗透
率,单位为m2;为第四区域渗透率,单位为m2;为第二区域扩散系数;为拉普拉
斯参数的相关函数;为双曲正切函数;为第四区域复合变量。
所述致密气藏五线性流模型中第一区域的线性流模型设置为:
(20)
其中,
(21)
(22)
(23)
(24)
式中,为拉普拉斯空间第三区域的无因次拟压力,为拉普拉斯空间水力
压裂区域的无因次拟压力;为第一区域扩散系数;为无因次裂缝宽度;为第三区
域渗透率,单位为m2;、、均为拉普拉斯参数的相关函数;为双曲正
弦函数;为第三区域复合变量。
所述水力压裂区域的线性流模型设置为:
(25)
其中,
(26)
(27)
式中,、均为拉普拉斯参数的相关函数。
根据水力压裂区域的线性流模型,得到定产量条件下的井底流压,如公式(28)所示:
(28)
式中,为定产量条件下的井底流压,无因次。
基于Duhamel原理,根据定产量条件下的井底流压,得到定井底流压条件下的
产出量,如公式(29)所示:
(29)
式中,为定压力条件下的产出量,无因次。
由此确定拉普拉斯空间内井底流压和产出量的函数关系,利用拉普拉斯逆变换数值模拟方法能够确定真实井筒空间内定产量生产时压裂水平井的井底流压以及定井底流压生产时压裂水平井的产量变化信息。
步骤3.3,设置未知参数组中各参数的取值范围和最大反演次数,利用致密气藏五线性流模型进行反演,得到未知参数组中各参数的最优值,具体包括以下步骤:
步骤3.3.1,根据压裂水平井的邻井资料和现场经验,分别设置未知参数组中各参数的取值范围,并预设最大反演次数和未知参数组中各参数的初始值,本实施例中最大反演次数设置为100。
步骤3.3.2,利用致密气藏五线性流模型进行反演,将步骤1中所获取致密油气藏的储层物性参数、压裂作业参数和生产历史输入至致密气藏五线性流模型中,再将预设的未知参数组中各参数的初始值输入至致密气藏五线性流模型中,利用致密气藏五线性流模型模拟致密油气藏的生产过程得到日产量的模拟值以及井底流压的模拟值,并根据生产制度决定计算致密油气藏日产量或致密油气藏井底流压的目标函数值。
本实施例中,所述致密油气藏日产量的目标函数,如公式(30)所示:
(30)
式中,为致密油气藏日产量的目标函数;为最大模拟时间,单位为s;为致密油气藏日产量的模拟值,单位为m3/s;为致密油气藏日产量的真实值,单
位为m3/s,根据致密油气藏的历史生产历史确定;为致密油气藏日产量目标函数的比率系
数。
所述致密油气藏井底流压的目标函数,如公式(31)所示:
(31)
其中,
(32)
(33)
(34)
(35)
式中,为致密油气藏井底流压的目标函数;为致密油气藏产量规整
化压力的模拟值,为致密油气藏模拟得到的产量规整化压力的导数;为
致密油气藏产量规整化压力的真实值,为致密油气藏产量规整化压力导数的真实
值;为致密油气藏井底流压目标函数的比率系数;为产量规整化压力的权重系数,
为产量规整化压力导数的权重系数。
步骤3.3.3,判断当前反演次数是否超过预设的最大反演次数,若当前迭代次数未超过预设的最大反演次数,则更新反演次数和未知参数组中各参数的取值,返回步骤3.3.2中继续利用致密气藏五线性流模型进行模拟,否则,则进入步骤3.3.4中。
步骤3.3.4,获取致密气藏五线性流模型各次反演模拟得到的致密油气藏日产量目标函数值,确定致密油气藏日产量目标函数值的最小值所对应未知参数组中各参数的取值,得到未知参数组的反演结果。
本实施例中,反演发现当裂缝渗透率取值为0.1348 mD、裂缝半长取值为45 m、储层改造区域宽度取值为3.63m、储层改造区域渗透率取值为0.0037 mD时致密油气藏日产量的目标函数值最小,即裂缝渗透率最优值为0.1348 mD、裂缝半长最优值为45 m、储层改造区域宽度最优值3.63m、储层改造区域渗透率最优值为0.0037 mD。
步骤4,将致密油气藏的各生产阶段分别划分为两个部分,建立生产制度转换模型,结合步骤1中所获取的致密油气藏储层物性参数、压裂作业参数、生产历史和步骤3中所确定的未知参数组反演结果,利用致密气藏五线性流模型模拟模拟得到各生产阶段中不同生产部分的日产量模拟值并进行数据拟合,得到致密油气藏的压裂水平井日产量预测模型,验证压裂水平井日产量预测模型的预测效果后,对验证后的压裂水平井日产量预测模型在时间上进行积分,得到致密油气藏的压裂水平井累计产量预测模型,具体包括以下步骤:
步骤4.1,由于常规油气藏变制度生产模拟时基于物质平衡方程确定平均压力以及初始化泄流面积区域内的压力,但是压裂后的储层具有强非均质特征,储层改造区域和储层未改造区域的之间的物性差异大,致密油气藏开发过程中两者的压力传播速度不同,使得常规反演方法进行模拟时仅采用平均压力初始化油气藏条件将产生严重误差。
所以,本实施例中针对致密油气藏的各生产阶段,将各生产阶段均划分为第一生产部分和第二生产部分,各生产阶段中的第一生产部分和第二生产部分采用不同的初始化压力进行开采,如图3所示。
本实施例中,建立生产制度转换模型如公式(36)所示:
(36)
式中,为第个生产阶段中第一生产部分生产制度的初始化压力,单位为Pa;为第个生产阶段中第一生产部分生产制度的初始化压力,单位为Pa;为第个生产阶段的平均压力,单位为Pa;为第个生产阶段中第二生产部分生产制度的初
始化压力,单位为Pa;为第一生产部分和第二生产部分之间生产制度的制度转换时间,单
位为s;、、均为制度转换参数。
步骤4.2,根据步骤1中所获取致密油气藏的储层物性参数、压裂作业参数、生产历史以及步骤3中所确定的裂缝渗透率最优值、裂缝半长最优值、储层改造区域宽度最优值、储层改造区域渗透率最优值重新设置致密气藏五线性流模型,利用致密气藏五线性流模型结合生产制度转换模型进行模拟,得到各生产阶段中第一生产部分和第二生产部分的日产量模拟值,对各生产阶段中第一生产部分的日产量模拟值和第二生产部分的日产量模拟值进行对数曲线拟合,如图4所示,得到致密油气藏的压裂水平井日产量预测模型。
步骤4.3,验证日产量预测模型的预测效果。
基于致密油气藏的生产历史和生产制度,利用生产制度转换模型计算得到致密油气藏各生产阶段的初始化压力和制度转换时间,再根据各生产阶段的初始化压力和制度转换时间,利用致密气藏五线性流模型模拟配合采用压裂水平井日产量预测模型计算得到变制度条件下致密油气藏中压裂水平井日产量的模拟值,并与生产历史中压裂水平井日产量的实测值进行对比,计算压裂水平井日产量的模拟值与生产历史中压裂水平井日产量实测值之间的精度值,若精度值小于10%,则调整生产制度转换模型中的各制度转换参数,更新致密油气藏的压裂水平井日产量预测模型后,返回步骤4.2中,否则,结束对压裂水平井日产量预测模型预测效果的验证,并进入步骤4.4。
步骤4.4,获取验证后致密油气藏的压裂水平井日产量预测模型,通过对致密油气藏的压裂水平井日产量预测模型在时间上进行积分,得到致密油气藏的压裂水平井累计产量预测模型。
本实施例中配合利用商业数值模拟软件CMG检验本发明中压裂水平井日产量预测
模型的可行性,利用致密气藏五线性流模型分别将井底流压设置为10MPa、8MPa、6MPa和
5MPa模拟压裂水平井的定压生产,选用井底流压为10MPa时压裂水平井的生产数据进行历
史拟合,获取致密气藏五线性流模型的裂缝渗透率、裂缝半长、储层改造区域宽度、储层改
造区域渗透率,利用井底流压为8MPa时压裂水平井的生产数据产量优选制度转换参数,优
选后确定制度转换参数取值为0.51124、制度转换参数取值为0.615586、制度转换参数取值为0.3647,此时致密油气藏的压裂水平井日产量预测模型精度值最高。
本实施例中,致密油气藏的压裂水平井共计生产576天,共划分为11个生产阶段,利用气藏的压裂水平井日产量预测模型结合生产制度转换模型进行模拟,预测得到压裂水平井日产量的模拟日产量和模拟累计产量,并与生产历史中压裂水平井的实际日产量和实际累计产量进行对比,如图5所示,对比后发现,本发明方法所建立的压裂水平井日产量预测模型预测结果准确,适用于致密油气藏的产能预测。
步骤5,利用验证后致密油气藏的压裂水平井日产量预测模型和压裂水平井累计产量预测模型进行产能预测,模拟不同生产制度条件下压裂水平井的开发动态,确定致密油气藏中压裂水平井的日产量预测值和累计产量预测值,从而得到致密油气藏的可采储量。
假设本实施例致密气藏中的压裂水平井生产20年后不再具有经济开采价值,利用验证后致密油气藏的压裂水平井累计产量预测模型进行产能预测,确定压裂水平井的累计产量为569万方,从而确定致密气藏中压裂水平井的可采储量为569万方,为后续致密气藏的生产制度调整优化提供了依据。
综上可得,采用本发明方法有效克服了传统解析模型的局限性,能够快速预测不同生产制度条件下压裂水平井的开发动态,准确致密油气藏压裂水平井的产量,有利于指导致密油气藏水平井开发。
在本发明描述中,需要说明的是,除非另有明确的规定和限定,术语“设置”、“安装”、“相连”、“连接”、“固定”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以通过具体情况理解上述术语在本发明中的具体含义。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (10)
1.一种致密油气藏压裂水平井变制度生产动态模拟监测方法,其特征在于,具体包括以下步骤:
步骤1,选取致密油气藏,获取致密油气藏的储层物性参数、压裂作业参数和生产历史,进行数据整理准备;
步骤2,根据生产历史中的井底流压和累计产量,将致密油气藏的生产历史划分为多个生产阶段;
步骤3,根据压裂水平井的特征,选取致密气藏五线性流模型作为复合线性流模型,结合压裂水平井的邻井资料和现场经验,利用致密气藏五线性流模型进行反演模拟获取由裂缝渗透率、裂缝半长、储层改造区域宽度、储层改造区域渗透率所组成的未知参数组,得到未知参数组中各参数的最优值;
步骤4,将致密油气藏的各生产阶段分别划分为两个部分,建立生产制度转换模型,结合步骤1中所获取的致密油气藏储层物性参数、压裂作业参数、生产历史和步骤3中所确定的未知参数组反演结果,利用致密气藏五线性流模型模拟得到各生产阶段中不同生产部分的日产量模拟值并进行数据拟合,得到致密油气藏的压裂水平井日产量预测模型,验证压裂水平井日产量预测模型的预测效果后,对验证后的压裂水平井日产量预测模型在时间上进行积分,得到致密油气藏的压裂水平井累计产量预测模型;
步骤5,利用验证后致密油气藏的压裂水平井日产量预测模型和压裂水平井累计产量预测模型预测变制度条件下致密油气藏内压裂水平井的生产动态,确定致密油气藏中压裂水平井的日产量预测值和累计产量预测值;
所述步骤3中,具体包括以下步骤:
步骤3.1,选取致密气藏五线性流模型作为复合线性流模型;
所述致密气藏五线性流模型包括第一区域、第二区域、第三区域、第四区域和第五区域,其中,第一区域为储层改造区域,第二区域、第三区域和第四区域均为储层未改造区域,第五区域为水力压裂区域,内部存在水力裂缝;
步骤3.2,建立天然气拟压力函数,将无因次变量引入致密气藏五线性流模型中,得到井底流压和产出量的计算公式;
步骤3.3,设置未知参数组中各参数的取值范围和最大反演次数,利用致密气藏五线性流模型进行反演,得到未知参数组中各参数的最优值;
所述步骤4中,具体包括以下步骤:
步骤4.1,将致密油气藏的各生产阶段均划分为第一生产部分和第二生产部分,建立生产制度转换模型;
步骤4.2,根据步骤1中所获取致密油气藏的储层物性参数、压裂作业参数、生产历史以及步骤3中所确定未知参数组的最优值重新设置致密气藏五线性流模型,利用致密气藏五线性流模型结合生产制度转换模型模拟,得到各生产阶段中第一生产部分和第二生产部分的日产量模拟值并进行对数曲线拟合,得到致密油气藏的压裂水平井日产量预测模型;
步骤4.3,基于致密油气藏的生产历史和生产制度,验证日产量预测模型的预测效果;
通过计算压裂水平井日产量模拟值与压裂水平井日产量实测值之间的精度值,若精度值小于预设精度值,则调整生产制度转换模型的制度转换参数,更新压裂水平井日产量预测模型后,返回步骤4.2中,否则,结束验证并进入步骤4.4;
步骤4.4,获取验证后的压裂水平井日产量预测模型,并进行时间上积分得到压裂水平井累计产量预测模型;
所述生产制度转换模型,如公式(36)所示:
(36)
式中,为第/>个生产阶段中第一生产部分生产制度的初始化压力,单位为Pa;/>为第/>个生产阶段中第一生产部分生产制度的初始化压力,单位为Pa;/>为第/>个生产阶段的平均压力,单位为Pa;/>为第/>个生产阶段中第二生产部分生产制度的初始化压力,单位为Pa;/>为第一生产部分和第二生产部分之间生产制度的制度转换时间,单位为s;、/>、/>均为制度转换参数。
2.根据权利要求1所述的致密油气藏压裂水平井变制度生产动态模拟监测方法,其特征在于,所述步骤1中,所述致密油气藏的储层物性参数包括储层的厚度、初始压力、温度、储层渗透率、孔隙度、饱和度以及储层内的油气密度;
所述致密油气藏的压裂作业参数包括压裂水平井的长度、裂缝簇数和井间距;
所述致密油气藏的生产历史包括生产时间、井底流压、产出速度、日产量和累计产量。
3.根据权利要求1所述的致密油气藏压裂水平井变制度生产动态模拟监测方法,其特征在于,所述步骤2中,同一生产阶段内压裂水平井的生产制度不变且压裂水平井的日产量和井底流压恒定。
4.根据权利要求1所述的致密油气藏压裂水平井变制度生产动态模拟监测方法,其特征在于,所述天然气拟压力函数为:
(1)
式中,为储层压力,单位为Pa;/>为天然气压缩因子;/>为粘度,单位为cp;
所述无因次变量包括无因次拟压力、无因次时间、无因次产量、无因次扩散系数和无因次距离;
所述无因次拟压力为:
(2)
式中,为第一区域的渗透率,单位为m2;/>为储层厚度,单位为m;/>为标况温度,取值为273.15K;/>为储层的初始压力,单位为Pa;/>为标况压力,取值为1.01325×105Pa;/>为地面产量,单位为m3/s;/>为储层温度,单位为K;
所述无因次时间为:
(3)
式中,为时间,单位为s;/>为第一区域的孔隙度;/>为第一区域的综合压缩系数,单位为Pa-1;/>为裂缝长度,单位为m;
所述无因次产量为:
(4)
式中,为井底流压,单位为Pa;
所述无因次扩散系数为:
(5)
其中,
(6)
式中,为区域编号,当/>时,表示第一区域,当/>时,表示第二区域,当/>时,表示第三区域,当/>时,表示第四区域,当/>时,表示第五区域;/>为第/>区域的扩散系数;/>为第/>区域的渗透率,单位为m2;/>为第/>区域的孔隙度;/>为第/>区域的综合压缩系数,单位为Pa-1。
5.根据权利要求4所述的致密油气藏压裂水平井变制度生产动态模拟监测方法,其特征在于,所述无因次距离包括无因次宽度、储层改造区域的无因次宽度、裂缝半间距的无因次宽度、无因次长度、储层改造区域的无因次长度、无因次井间半间距长度和无因次裂缝导流能力,其中,所述无因次宽度为:
(7)
式中,为无因次宽度;/>为距离压裂裂缝的距离,单位为m;
所述储层改造区域的无因次宽度为:
(8)
式中,为储层改造区域的无因次宽度;/>为储层改造区域的宽度,单位为m;
所述裂缝半间距的无因次宽度为:
(9)
式中,为裂缝半间距无因次宽度;/>为裂缝半间距,单位为m;
所述无因次长度为:
(10)
式中,为无因次长度;/>为距离水平井筒的距离,单位为m;
所述储层改造区域的无因次长度为:
(11)
式中,为储层改造区域的无因次长度;/>为储层改造区域的长度,单位为m;
所述无因次井间半间距长度为:
(12)
式中,为无因次井间半间距长度;/>为井间半间距的长度,单位为m;
所述无因次裂缝导流能力为:
(13)
式中,为无因次裂缝导流能力;/>为裂缝渗透率,单位为m2;/>为裂缝宽度,单位为m。
6.根据权利要求5所述的致密油气藏压裂水平井变制度生产动态模拟监测方法,其特征在于,将无因次变量引入致密气藏五线性流模型中,所述第三区域和第四区域的线性流模型均为:
(14)
其中,
(15)
(16)
式中,为拉普拉斯空间第/>区域的无因次拟压力;/>为拉普拉斯变换参数;/>为拉普拉斯空间第/>区域的无因次拟压力;/>为双曲余弦函数;/>为第/>区域复合变量;
所述第二区域的线性流模型为:
(17)
其中,
(18)
(19)
式中,为拉普拉斯空间第一区域的无因次拟压力,/>为拉普拉斯空间第二区域的无因次拟压力,/>为拉普拉斯空间第四区域的无因次拟压力;/>为第二区域渗透率;为第四区域渗透率;/>为第二区域扩散系数;/>为拉普拉斯参数的相关函数;为双曲正切函数;/>为第四区域复合变量;
所述第一区域的线性流模型为:
(20)
其中,
(21)
(22)
(23)
(24)
式中,为拉普拉斯空间第三区域的无因次拟压力,/>为拉普拉斯空间水力压裂区域的无因次拟压力;/>为第一区域扩散系数;/>为无因次裂缝宽度;/>为第三区域渗透率;/>、/>、/>均为拉普拉斯参数的相关函数;/>为双曲正弦函数;/>为第三区域复合变量。
7.根据权利要求6所述的致密油气藏压裂水平井变制度生产动态模拟监测方法,其特征在于,所述水力压裂区域的线性流模型为:
(25)
其中,
(26)
(27)
式中,、/>均为拉普拉斯参数/>的相关函数;
根据水力压裂区域的线性流模型,得到定产量条件下的井底流压,如公式(28)所示:
(28)
式中,为定产量条件下的井底流压,无因次;
基于Duhamel原理,根据定产量条件下的井底流压,得到定井底流压条件下的产出量,如公式(29)所示:
(29)
式中,为定压力条件下的产出量,无因次。
8.根据权利要求7所述的致密油气藏压裂水平井变制度生产动态模拟监测方法,其特征在于,所述步骤3.3中,具体包括以下步骤:
步骤3.3.1,根据压裂水平井的邻井资料和现场经验,分别设置未知参数组中各参数的取值范围,并预设最大反演次数和未知参数组中各参数的初始值;
步骤3.3.2,将步骤1中所获取致密油气藏的储层物性参数、压裂作业参数和生产历史和预设的未知参数组的初始值输入至致密气藏五线性流模型中,利用致密气藏五线性流模型模拟致密油气藏的生产过程得到日产量的模拟值以及井底流压的模拟值,并根据生产制度决定计算致密油气藏日产量或致密油气藏井底流压的目标函数值;
步骤3.3.3,判断当前反演次数是否超过预设最大反演次数,若当前迭代次数未超过预设最大反演次数,则更新反演次数和未知参数组的取值,返回步骤3.3.2中继续利用致密气藏五线性流模型进行模拟,否则,则进入步骤3.3.4中;
步骤3.3.4,获取致密气藏五线性流模型各次反演模拟得到的致密油气藏日产量目标函数值,确定致密油气藏日产量目标函数值的最小值所对应未知参数组中各参数的取值,得到未知参数组的反演结果,包括裂缝渗透率最优值、裂缝半长最优值、储层改造区域宽度最优值、储层改造区域渗透率最优值。
9.根据权利要求8所述的致密油气藏压裂水平井变制度生产动态模拟监测方法,其特征在于,所述致密油气藏日产量的目标函数,如公式(30)所示:
(30)
式中,为致密油气藏日产量的目标函数;/>为最大模拟时间,单位为s;/>为致密油气藏日产量的模拟值,单位为m3/s;/>为致密油气藏日产量的真实值,单位为m3/s,根据致密油气藏的历史生产历史确定;/>为致密油气藏日产量目标函数的比率系数;
所述致密油气藏井底流压的目标函数,如公式(31)所示:
(31)
其中,
(32)
(33)
(34)
(35)
式中,为致密油气藏井底流压的目标函数;/>为致密油气藏产量规整化压力的模拟值,/>为致密油气藏模拟得到的产量规整化压力的导数;/>为致密油气藏产量规整化压力的真实值,/>为致密油气藏产量规整化压力导数的真实值;为致密油气藏井底流压目标函数的比率系数;/>为产量规整化压力的权重系数,/>为产量规整化压力导数的权重系数。
10.根据权利要求1所述的致密油气藏压裂水平井变制度生产动态模拟监测方法,其特征在于,所述步骤5中,通过改变生产制度,利用生产制度转换模型确定生产制度中各生产阶段之间的制度转换时间,结合致密气藏五线性流模型模拟得到致密油气藏内压裂水平井在变制度条件下的生产动态,利用验证后的压裂水平井日产量预测模型和压裂水平井累计产量预测模型进行预测,得到不同生产制度条件下致密油气藏中压裂水平井的日产量预测值和累计产量预测值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311675971.1A CN117371361B (zh) | 2023-12-08 | 2023-12-08 | 一种致密油气藏压裂水平井变制度生产动态模拟监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311675971.1A CN117371361B (zh) | 2023-12-08 | 2023-12-08 | 一种致密油气藏压裂水平井变制度生产动态模拟监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117371361A CN117371361A (zh) | 2024-01-09 |
CN117371361B true CN117371361B (zh) | 2024-03-19 |
Family
ID=89389592
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311675971.1A Active CN117371361B (zh) | 2023-12-08 | 2023-12-08 | 一种致密油气藏压裂水平井变制度生产动态模拟监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117371361B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105003239A (zh) * | 2014-04-15 | 2015-10-28 | 中国海洋石油总公司 | 海上压裂水平井压后效果评价方法 |
WO2016192077A1 (zh) * | 2015-06-04 | 2016-12-08 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种致密气压裂水平井数值试井模型建立求解方法 |
CN107506948A (zh) * | 2017-09-27 | 2017-12-22 | 科为联合(青岛)能源技术研究院有限公司 | 一种基于动态泄流体积的页岩油气综合产量分析方法 |
CN111428425A (zh) * | 2020-03-19 | 2020-07-17 | 西南石油大学 | 一种页岩油藏变裂缝渗透率分段压裂水平井产能计算方法 |
CN111980654A (zh) * | 2020-10-12 | 2020-11-24 | 西南石油大学 | 一种非均匀页岩油藏分段压裂水平井产能计算方法 |
AU2020103953A4 (en) * | 2020-12-08 | 2021-02-11 | Southwest Petroleum University | Method and system for predicting production of fractured horizontal well in shale gas reservoir |
CN112360422A (zh) * | 2020-12-08 | 2021-02-12 | 西南石油大学 | 一种页岩气藏压裂水平井产量预测方法及系统 |
CN116291331A (zh) * | 2023-05-05 | 2023-06-23 | 中国石油大学(北京) | 油气藏产能预测方法、装置、设备和存储介质 |
CN116384267A (zh) * | 2023-02-21 | 2023-07-04 | 中国石油大学(北京) | 致密储层压裂水平井的最终可采储量的确定方法及设备 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110374572B (zh) * | 2019-08-12 | 2021-01-05 | 西南石油大学 | 非常规油气藏多段压裂水平井各改造段特性参数确定方法 |
-
2023
- 2023-12-08 CN CN202311675971.1A patent/CN117371361B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105003239A (zh) * | 2014-04-15 | 2015-10-28 | 中国海洋石油总公司 | 海上压裂水平井压后效果评价方法 |
WO2016192077A1 (zh) * | 2015-06-04 | 2016-12-08 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种致密气压裂水平井数值试井模型建立求解方法 |
CN107506948A (zh) * | 2017-09-27 | 2017-12-22 | 科为联合(青岛)能源技术研究院有限公司 | 一种基于动态泄流体积的页岩油气综合产量分析方法 |
CN111428425A (zh) * | 2020-03-19 | 2020-07-17 | 西南石油大学 | 一种页岩油藏变裂缝渗透率分段压裂水平井产能计算方法 |
CN111980654A (zh) * | 2020-10-12 | 2020-11-24 | 西南石油大学 | 一种非均匀页岩油藏分段压裂水平井产能计算方法 |
AU2020103953A4 (en) * | 2020-12-08 | 2021-02-11 | Southwest Petroleum University | Method and system for predicting production of fractured horizontal well in shale gas reservoir |
CN112360422A (zh) * | 2020-12-08 | 2021-02-12 | 西南石油大学 | 一种页岩气藏压裂水平井产量预测方法及系统 |
CN116384267A (zh) * | 2023-02-21 | 2023-07-04 | 中国石油大学(北京) | 致密储层压裂水平井的最终可采储量的确定方法及设备 |
CN116291331A (zh) * | 2023-05-05 | 2023-06-23 | 中国石油大学(北京) | 油气藏产能预测方法、装置、设备和存储介质 |
Non-Patent Citations (4)
Title |
---|
致密气藏压裂水平井产能分析及压裂优化设计;孙元伟;程远方;时凤霞;王艳丽;印树明;;新疆石油地质(第06期);全文 * |
致密气藏压裂水平井反卷积试井模型;王飞;潘子晴;;石油学报(第07期);全文 * |
致密油藏体积压裂水平井分区复合产能模型;任岚;黄鑫;;大庆石油地质与开发(第02期);全文 * |
页岩气藏多段压裂水平井产能预测模型;谢亚雄;刘启国;王卫红;李若莹;胡小虎;;大庆石油地质与开发(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN117371361A (zh) | 2024-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2007265695B2 (en) | Method for comparing and back allocating production | |
CN107506948B (zh) | 一种基于动态泄流体积的页岩油气综合产量分析方法 | |
CN108694254B (zh) | 一种变产变压生产气井产量经验递减曲线分析方法 | |
CN108804819A (zh) | 一种低渗气藏动态储量评价方法 | |
CN111456709B (zh) | 一种基于测井曲线的水平井多级压裂分段分簇方法 | |
CN114427432B (zh) | 一种气藏剩余气开发潜力确定方法 | |
CN109426672B (zh) | 基于不确定地质模型的油藏注采参数优化方法 | |
CN109488282A (zh) | 可动储量物性下限的确定方法 | |
CN109057786A (zh) | 一种用于非常规油气藏可采储量估算的自适应模型 | |
CN111927417A (zh) | 一种页岩气分段压裂水平井组储量动用状况评价方法 | |
CN111287740A (zh) | 基于真实应变下的异常高压气藏动态储量计算方法 | |
CN106640042A (zh) | 一种评价气井单井产能装置及方法 | |
CN111502652A (zh) | 一种三孔介质气藏水平井产量递减及生产动态预测方法 | |
CN117371361B (zh) | 一种致密油气藏压裂水平井变制度生产动态模拟监测方法 | |
CN112610205B (zh) | 一种强非均质储层气井早期产能评价方法 | |
CN111950112A (zh) | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 | |
CN113445988A (zh) | 一种低渗碳酸盐岩气藏气井产能评价方法 | |
CN111950111B (zh) | 一种适用于底部开放的碳酸盐岩储层动态分析方法 | |
CN116050626B (zh) | 一种多级压裂水平井生产分级优化与效果预测方法 | |
CN116291331A (zh) | 油气藏产能预测方法、装置、设备和存储介质 | |
CN112035993A (zh) | 一种底部定压的碳酸盐岩储层测试评价方法 | |
CN111222261B (zh) | 动静参数结合的产量劈分方法 | |
CN110984976B (zh) | 多层合采气井产量劈分方法 | |
CN113153220A (zh) | 一种气井井筒疏通效果评价方法 | |
CN113901378A (zh) | 压力梯度影响因素权重系数确定方法 |
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 |