CN108694254A - 一种变产变压生产气井产量经验递减曲线分析方法 - Google Patents
一种变产变压生产气井产量经验递减曲线分析方法 Download PDFInfo
- Publication number
- CN108694254A CN108694254A CN201710218995.2A CN201710218995A CN108694254A CN 108694254 A CN108694254 A CN 108694254A CN 201710218995 A CN201710218995 A CN 201710218995A CN 108694254 A CN108694254 A CN 108694254A
- Authority
- CN
- China
- Prior art keywords
- gas well
- fbhp
- production
- flowing bottomhole
- bottomhole pressure
- 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
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 79
- 230000007423 decrease Effects 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000009466 transformation Effects 0.000 title claims abstract description 23
- 238000013211 curve analysis Methods 0.000 title claims abstract description 20
- 239000000463 material Substances 0.000 claims abstract description 37
- 230000003247 decreasing effect Effects 0.000 claims abstract description 28
- 238000011161 development Methods 0.000 claims abstract description 15
- 230000008859 change Effects 0.000 claims abstract description 5
- 230000001186 cumulative effect Effects 0.000 claims description 8
- 238000004088 simulation Methods 0.000 claims description 4
- 238000011958 production data acquisition Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 28
- 238000011156 evaluation Methods 0.000 abstract description 6
- 238000007405 data analysis Methods 0.000 abstract description 2
- 238000011160 research Methods 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 126
- 239000007788 liquid Substances 0.000 description 14
- 238000004364 calculation method Methods 0.000 description 13
- 230000018109 developmental process Effects 0.000 description 13
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 238000007796 conventional method Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 239000003345 natural gas Substances 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000005325 percolation Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000003795 desorption Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000009776 industrial production Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000005514 two-phase flow 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
- 230000001133 acceleration Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011067 equilibration Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 235000000334 grey box Nutrition 0.000 description 1
- 244000085685 grey box Species 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002699 waste 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
-
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (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
本发明提供一种变产变压生产气井产量经验递减曲线分析方法,该方法引入拟压力规整化日产气量和物质平衡时间,建立起一种生产数据分析方法,解决了目前国内大多数页岩气井在投产初期要求配产有一定稳产期因而无法通过常规递减分析方法进行产量预测研究的技术问题,考虑到页岩气井独特的产量变化规律,引入幂律指数递减模型,给出了一种页岩气井初期产能分析及后期产量预测方法。本发明在对生产数据进行分析评价时不必要求气井处于定压递减生产阶段,在根据分析模型进行未来产能预测时可以同时预测气井产量、井底流压随时间变化。本发明的评价结果可以用于页岩气井产能评价和开发指标预测。
Description
技术领域
本发明属于致密气、页岩气勘探开发领域,主要涉及一种变产变压生产气井产量经验递减曲线分析方法。
背景技术
页岩气是我国天然气勘探开发的重要新领域之一,目前已经有许多页岩气井投产,并积累了大量的试采动态数据,这些数据的有效分析是指导气井下一步生产极为重要的依据。
对于一口气井来说,在生产进入递减阶段以后,可以在分析生产历史数据的基础上评价确定生产气井未来最可能的寿命,评估未来产量。外推法是最简单和最直接的方式,具体做法是绘制出产量相对于时间或累计产量的曲线,并将曲线进行外推,递减曲线分析和指标预测是其中非常重要的一项内容。气井产量递减分析是通过分析气井生产数据,确定产量递减模式、初始递减率以及递减指数,并进行气井动态预测,准确分析气井产量变化规律是评价气井开发效果及预测最终可采储量的重要依据。
页岩气藏属于低孔、特低渗气藏,有着特殊的储集运移机理,需要大型水力压裂之后才有工业产能,这决定了其动态分析要比常规气藏复杂得多,表现为气井生产初期阶段产出的主要是裂缝系统中的游离气,气井初始产量较高但是递减较快;之后的产量主要来自基质孔隙中的自由气和吸附气,由于基质中的天然气解吸扩散-渗流速度要慢得多,气井逐渐进入低产稳产阶段。这也导致页岩气井生产动态区别于常规气井,一般具有很长的产量递减期,即使生产数年后也很难达到拟稳定流,因而适用于常规气藏的传统产量递减分析方法在应用时受到限制。
目前国外学者考虑到页岩气藏独特的渗流规律,在常规分析方法的基础上引入了幂律指数、扩展指数等其他递减模式,并在美国德克萨斯州中北部沃思堡(Fort Worth)盆地巴奈特(Barnett)页岩气藏以及路易斯安那州海恩斯维尔(Haynesville)页岩气藏数口井的生产动态分析中进行了应用。国内方面页岩气井产量递减分析研究一直以常规的Arps(Arps提出了三种递减规律,分别为指数递减、双曲线递减和调和递减)递减分析为主,近几年也开始关注并尝试使用幂律、扩展指数等新型递减模型进行气井产量变化规律分析和预测。
但不管是常规的还是新型的递减模型在应用时均要求气井达到边界控制流动,进入递减阶段半年以上、且井底压力变化不大,否则无法评价气井产能,并且当利用通过历史生产阶段产量分析而建立起的递减模型进行预测时,井底流压也要基本保持不变,否则会引起严重误差。而目前国内的页岩气井多采用限产方式投产,且工作制度不稳定。除此之外,页岩气藏的特殊性也决定了其难以达到拟稳态流动阶段。这一方面造成传统方法在页岩气井产量递减分析时不再适用,另一方面,使得目前现场已经积累的大量的生产动态数据无法用于动态预测,造成资源的浪费。
分析常规的幂律指数等页岩气藏产量递减分析技术可以发现,针对大量变产变压生产的页岩气井,目前缺乏有效准确的产量递减分析方法,尤其是在页岩气藏开发的早期阶段。
发明内容
为解决上述现有技术中,特别是在页岩气藏开发的早期阶段,缺乏有效准确的页岩气井产量递减分析方法的技术问题,本发明提供一种变产变压生产气井产量经验递减曲线分析方法,具体方案如下:
一种变产变压生产气井产量经验递减曲线分析方法,包括以下步骤:
生产数据采集、计算,并对日产气量、井底流压和生产时间数据点进行规整化处理,以获得拟压力规整化日产气量和物质平衡时间;
使用幂律指数递减模型拟合拟压力规整化日产气量和物质平衡时间散点数据,得到产量经验递减模型的参数,以确定所述产量经验递减模型;
根据预设工作制度预测气井未来开发指标变化。
进一步的,在拟合过程中,所述拟压力规整化日产气量散点数据包括以井口油套压计算的所述井底流压进行规整化处理而得到的数据。
进一步的,在拟合过程中,当存在实测井底流压,则所述拟压力规整化日产气量散点数据还包括以实测井底流压进行规整化处理而得到的数据。
进一步的,所述预设工作制度包括设置预测的最大时间范围,设置气井稳产期日产气量,设置最低井底流压,设置模拟时步步长。
进一步的,以现场统计的最后一个数据点对应的物质平衡时间为起点,通过时步步长推算其后的物质平衡时间,根据经验产量递减模型计算所述推算的物质平衡时间对应的拟压力规整化日产气量,并得到相应的拟井底流压,通过查表得出所述拟井底流压对应的井底流压。
进一步的,当经查表获得的所述井底流压高于所述最低井底流压,则该气井可以按照定产方式生产。
进一步的,当经查表获得的所述井底流压低于或等于所述最低井底流压,则该气井以最低井底流压生产,并在井底流压等于所述最低井底流压条件下重新计算气井日产气量,进而得到当前时步的气井累积产气量。
进一步的,根据所述推算的物质平衡时间计算相应的实际生产时间。
进一步的,当所述实际生产时间小于或等于所述预测的最大时间范围,则进行下一时步的指标预测。
进一步的,当所述实际生产时间大于所述预测的最大时间范围,则计算结束,并输出页岩气井指标预测结果。
与现有技术相比,本发明提供的一种变产变压生产气井产量经验递减曲线分析方法通过引入拟压力规整化日产气量的概念,消除了常规方法在进行递减分析时要求气井井底流压波动不大的限制,因而可以进行变产变压生产气井的递减曲线分析。另外,该方法使用幂律指数递减模型拟合分析数据,由于该模型的分段特性,可以考虑页岩气井日产气量早期递减快、后期递减慢的分阶段递减特征。再者,常规方法不考虑井底流压的变化对产量递减分析的影响,要求在分析阶段与未来产量预测阶段气井的井底流压变化不大,而本方法通过引入的拟压力规整化日产气量可以考虑井底流压变化对气井日产气量的影响。因此,本发明通过上述技术方案消除了常规气井产量递减分析方法应用的局限性,很好的解决早期变产变压生产气井产量数据分析及生产指标预测的问题。
附图说明
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中:
图1为本发明的变产变压生产气井产量经验递减曲线分析、预测流程图;
图2为本发明实施例中A-1井的生产动态数据;
图3为本发明实施例中A-1井的井底流压、套压变化曲线;
图4为本发明实施例中散点拟合得到的幂律指数递减公式;
图5为本发明实施例中预测得到的拟压力规整化日产气量曲线;
图6为本发明实施例中气井未来开发指标变化的预测结果。
在附图中,相同的部件采用相同的附图标记,附图并未按实际比例绘制。
具体实施方式
下面将结合附图对本发明作进一步说明。
页岩气藏属于低孔、特低渗气藏,有着特殊的储集运移机理,需要大型水力压裂之后才有工业产能,这决定了其动态分析要比常规气藏复杂得多,表现为气井生产初期阶段产出的主要是裂缝系统中的游离气,气井初始产量较高但是递减较快;之后的产量主要来自基质孔隙中的自由气和吸附气,由于基质中的天然气解吸扩散-渗流速度要慢得多,气井逐渐进入低产稳产阶段。这也导致页岩气井生产动态区别于常规气井,一般具有很长的产量递减期,即使生产数年后也很难达到拟稳定流,因而适用于常规气藏的传统产量递减分析方法在应用时受到限制。本实施例提供一种变产变压生产气井产量经验递减曲线分析方法,通过引入拟压力规整化日产气量,利用幂律指数递减模型进行散点拟合,并预测气井未来开发指标,很好地解决了上述问题。下文将以A-1井为例进行详细分析。
本实施例各参数的单位中的d为天数,注意区分管子内径的表示符号d。如无特别说明,本实施例中所提到的时间,均以气井开始投产之日为时间起算的参考日,例如,评价时间30年,是指自气井开始投产之日起30年。
收集A-1井的生产参数;收集并设置地质参数;测试并设置页岩气PVT参数(PVT参数,即天然气高压物性参数);收集气井管柱及井眼轨迹数据;得到本实施例A-1井的相关参数如表1所示,本实施例A-1井的气井生产动态数据如图2所示。
表1:A-1井相关参数表
根据井筒管柱结构、气井生产方式由井口油套压、产气量、产水量计算储层中深位置的井底流压。此处用Gray模型(即灰箱模型)计算井底流压,当气井不产水时,持液率参数Hl等于0,Gray模型以干气方式计算井底流压;当气井产水时,持液率参数Hl大于0,Gray模型以两相流模型来计算井底流压。压力梯度计算公式如下:
式(1)中,ρm=ρlHl+ρg(1-Hl) (2)
式(2)中,
式(3)中,
式(4)中,
式(4)中,
式(4)中,
式(2)、(4)中,
式(1)--(8)中,z:井筒上某一点的深度,单位m;p:深度z处对应的井筒压力,单位Pa;f:两相摩阻系数,无因次;vm:气液混合物表观速度,单位m/s;d:管子内径,单位m;ρNS:无滑脱气液混合物密度,单位kg/m3;ρm:混合物密度,单位kg/m3;g:重力加速度,单位m/s2;ρg:气体密度,kg/m3;ρl:液体密度,kg/m3;Hl:持液率,无因次;Vsl气井产液速度,单位m/s、Vsg:气井产气速度,单位m/s;σl:气液间的界面张力,单位mN/m。
本实施例A-1井,井筒管柱结构为2500m下深的油管,内径62mm,外径73mm,其生产方式为油管生产,且产水,所述Gray模型以气液两相流模型,并根据井口油套压采用上述压力梯度计算公式计算井底流压,计算结果如图3所示。
对日产气量、井底流压和生产时间数据点进行规整化处理,计算拟压力规整化日产气量和物质平衡时间,相应计算公式如下:
拟压力规整化日产气量计算公式:
式(9)中,
式(9)--(10)中,PNR为拟压力规整化日产气量,单位104m3/(d·MPa2/mPa·s);t为实际生产时间,单位d;qg(t)为t时刻的气井日产气量,单位104m3/d;ψi原始地层压力对应的拟压力,单位MPa2/mPa·s;ψ(t)是t时刻的井底流压对应的拟压力,单位MPa2/mPa·s;p0是参考压力,可取为一个标准大气压,单位MPa;p是井底流压,单位MPa;μ是气体粘度,单位mPa.s;z是气体压缩系数,无因次。
通过引入拟压力规整化日产气量的概念,可以考虑井底流压变化对气井日产气量的影响,消除了常规方法在进行递减分析时要求气井井底流压波动不大的限制,可以进行变产变压生产气井的递减曲线分析。
物质平衡时间计算公式:
式(11)中tmb为物质平衡时间,单位d;qg(t)为t时刻的气井日产气量,单位104m3/d;t为实际生产时间,单位d。
将A-1井收集到的数据,及式(1)--(8)计算所得的井底流压数据,代入式(9)--(11)中,得到一系列的拟压力规整化日产气量与物质平衡时间的散点数据,该散点数据如图4中的散点所示。
使用幂律指数递减模型拟合拟压力规整化日产气量和物质平衡时间散点数据,确定递减模型参数。
该幂律指数递减模型的拟合目标函数如下:
式(12)中,qa,D∞,Da,n为待求的拟合参数,obj(qa,D∞,Da,n)表示包含qa,D∞,Da,n这四个参数的目标函数,参数qa表示初始时刻(即气井开始投产的时刻)的拟压力规整化日产气量,单位104m3/(d·MPa2/mPa·s);参数D∞表示生产后期的递减率,单位d-1,参数Da表示生产初期的递减率,单位d-1,n为递减指数,无因次;n0表示实测井底流压点个数;PNR(tmb,i)表示第i个实测流压点对应的拟压力规整化日产气量值,单位104m3/(d·MPa2/mPa·s);tmb,i表示第i个实测流压点对应的物质平衡时间;单位d;m0表示井口油套压折算井底流压点的个数;PNRz(tzmb,j)表示第j个折算井底流压点对应的拟压力规整化日产气量值,单位104m3/(d·MPa2/mPa·s);tzmb,j表示第j个折算井底流压点对应的物质平衡时间。单位d;公式(12)表示寻找最优的参数组合,使平均单点拟合误差平方值最小,也就是通常所说的最小二乘法。
由于幂律指数递减模型的分段特性,因此,该模型可以考虑页岩气井日产气量早期递减快、后期递减慢的分阶段递减特征,符合气井预测的现实需求。
优选的,在拟合过程中,主要是以井口油套压计算的井底流压规整化而得到的拟压力规整化日产气量为主,如果有实测井底流压,则该实测井底流压规整化而得到的拟压力规整化日产气量也可以放入拟合目标函数来降低井底流压计算偏差对拟合结果影响。
将图4中的散点数据利用拟合目标函数进行散点拟合,通过规划求解,拟合得到A-1井的产量经验递减模型(也即是幂律指数递减模型)为:
PNR(tmb)=51.89exp(-1.76tmb 0.126) (13)
式(13)中,tmb为物质平衡时间,单位d;PNR(tmb)为物质平衡时间tmb所对应的拟压力规整化日产气量,单位104m3/(d·MPa2/mPa·s)。
在拟合得到的A-1井的产量经验递减模型的基础上,根据预设的工作制度就可以预测A-1井未来的开发指标。下文将详细阐明A-1井的预测步骤。
步骤(1):
A-1井预设工作制度如下:
设置气井的稳产产量为收集到的气井生产末期产量,即q0=6.1×104m3/d,最低井底流压限定值为pmin=7MPa,设置预测的最大时间范围,tend=30年,设置模拟时步步长Δtc=30天;
工作制度的预设为后续的预测步骤提供了基本的前提条件,在不同的工作制度下可以获得不同的预测结果。
步骤(2):
利用模拟时步步长,以现场统计的A-1井最后一个生产数据点对应的物质平衡时间为指标预测起点,即tmbend=818天,往后推算下一个物质平衡时间,计算公式如下:
tmbn=tmbend+n×Δtc,n=1,2,3,… (14)
式(14)中,n为自tmbend=818天为起点,往后推算的物质平衡时间的个数,无因次;tmbn为自tmbend=818天为起点,往后推算的第n个物质平衡时间,单位d。
步骤(3):
利用A-1井的产量经验递减模型,计算tmbn时刻对应的拟压力规整化日产气量,计算公式如下:
PNR(tmbn)=51.89exp(-1.76tmbn 0.126) (15)
式(15)中,tmbn为自tmbend=818天为起点,往后推算的第n个物质平衡时间,单位d;PNR(tmbn)为物质平衡时间tmbn所对应的拟压力规整化日产气量,单位104m3/(d·MPa2/mPa·s)。
由式(15)计算物质平衡时间tmbn时刻对应的拟压力规整化日产气量,所得数据如图5所示;
步骤(4):
按A-1井稳产产量为q0,计算任意时刻对应的拟井底流压ψ(tmbn),计算公式如下:
式(16)中,A-1井稳产产量q0=6.1×104m3/d;PNR(tmbn)为物质平衡时间tmbn所对应的拟压力规整化日产气量,单位104m3/(d·MPa2/mPa·s)。ψ(tmbn)为物质平衡时间tmbn时刻对应的拟井底流压,单位MPa2/mPa·s。
步骤(5):
根据压力-拟压力表(该表为常规参数表,可通过式(10)计算获得)插值求出拟井底流压ψ(tmbn)对应的井底流压pwf(tmbn);若井底流压值pwf(tmbn)高于规定的最低井底流压值pmin,则该井可以按照定产方式生产,q(tmbn)=q0;若井底流压值pwf(tmbn)低于或等于规定的最低井底流压值pmin,表明此时气井以最低井底流压生产,此时对应的时间即是气井的稳产期ts,再以pwf(tmbn)=pmin重新计算气井日产气量,计算公式如下:
q(tmbn)=PNR(tmbn)×(ψi-ψmin) (17)
式(17)中,q(tmbn)为物质平衡时间tmbn时刻对应的气井日产气量,单位104m3/d;PNR(tmbn)为物质平衡时间tmbn时刻所对应的拟压力规整化日产气量,单位104m3/(d·MPa2/mPa·s)。ψi为原始地层压力对应的拟压力,单位MPa2/mPa·s;ψmin为最低井底流压所对应的拟井底流压,单位MPa2/mPa·s;
本实施例中根据式(16)计算得到的A-1井的拟井底流压ψ(tmbn),通过压力-拟压力表插值求出拟井底流压ψ(tmbn)对应的井底流压pwf(tmbn),计算结果表明当tmbn=707d时,计算出的井底压力pwf(tmbn)<7MPa,故当tmbn>707d时,气井应该按照定井底流压的方式生产,并根据pwf(tmbn)=7MPa,应用式(17)重新计算该时刻气井日产气量q(tmbn),计算所得数据如图6所示;
步骤(6):
计算当前时步累积产气量。累积产气量可以通过气井日产气量和物质平衡时间进行计算,公式如下:
Qc(tn)=tmbn×q(tmbn) (18)
式(18)中,Qc(tn)表示第n(此处n与式(14)中的n含义一致)个数据点对应的气井累积产量,单位108m3;tmbn与式(14)中的tmbn含义一致;q(tmbn)与式(17)中的q(tmbn)含义一致。
根据A-1井计算所得当前时步气井日产气量q(tmbn)和物质平衡时间tmbn,应用式(18)计算A-1井当前时步的累积产气量,计算所得数据如图6所示。
步骤(7):
根据物质平衡时间tmbn计算对应的实际生产时间tn,计算公式如下:
式(19)中,tn为当前时步的实际生产时间,单位d;tn-1为前一时步的实际生产时间,单位d;q(tmbn)与式(17)中的q(tmbn)含义一致;tmbn与式(14)中的tmbn含义一致;q(tmbn-1)为前一时步物质平衡时间对应的气井日产气量,单位104m3/d;tmbn-1为前一时步物质平衡时间,单位d。
步骤(8):
若根据式(19)算得的A-1井的实际生产时间tn大于tend=30年,则计算结束,否则按照步骤(2)~(7)进行下一时步指标预测;
步骤(9):
输出页岩气井指标预测结果,包括气井日产气量、井底流压、累积产气量。模拟气井生产30年,预测得到的A-1井各时间点气井日产气量、井底流压、累积产气量如图6所示。A-1井稳产期为1.94年,30年末气井日产气量为0.44×104m3/d,井底流压为7Mpa,累产气量为1.93×108m3/d。
本实施例分析所得预测结果准确可靠,为A-1井开发效果及最终可采储量的评价提供了重要的依据。
虽然已经参考优选实施例对本发明进行了描述,但在不脱离本发明的范围的情况下,可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是,只要不存在结构冲突,各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文中公开的特定实施例,而是包括落入权利要求的范围内的所有技术方案。
Claims (10)
1.一种变产变压生产气井产量经验递减曲线分析方法,其特征在于,包括以下步骤:
生产数据采集、计算,并对日产气量、井底流压和生产时间数据点进行规整化处理,以获得拟压力规整化日产气量和物质平衡时间;
使用幂律指数递减模型拟合拟压力规整化日产气量和物质平衡时间散点数据,得到产量经验递减模型的参数,以确定所述产量经验递减模型;
根据预设工作制度预测气井未来开发指标变化。
2.根据权利要求1所述的变产变压气井产量经验递减曲线分析方法,其特征在于,在拟合过程中,所述拟压力规整化日产气量散点数据包括以井口油套压计算的所述井底流压进行规整化处理而得到的数据。
3.根据权利要求2所述的变产变压气井产量经验递减曲线分析方法,其特征在于,在拟合过程中,当存在实测井底流压,则所述拟压力规整化日产气量散点数据还包括以实测井底流压进行规整化处理而得到的数据。
4.根据权利要求1所述的变产变压气井产量经验递减曲线分析方法,其特征在于,所述预设工作制度包括设置预测的最大时间范围,设置气井稳产期日产气量,设置最低井底流压,设置模拟时步步长。
5.根据权利要求4所述的变产变压气井产量经验递减曲线分析方法,其特征在于,以现场统计的最后一个生产数据点对应的物质平衡时间为起点,通过设置的时步步长推算其后的物质平衡时间,根据经验产量递减模型计算所述推算的物质平衡时间对应的拟压力规整化日产气量,并得到相应的拟井底流压,通过查表得出所述拟井底流压对应的井底流压。
6.根据权利要求5所述的变产变压气井产量经验递减曲线分析方法,其特征在于,当经查表获得的所述井底流压高于所述最低井底流压,则该气井可以按照定产方式生产。
7.根据权利要求5所述的变产变压气井产量经验递减曲线分析方法,其特征在于,当经查表获得的所述井底流压低于或等于所述最低井底流压,则该气井以最低井底流压生产,并在井底流压等于所述最低井底流压条件下重新计算气井日产气量,进而得到当前时步的气井累积产气量。
8.根据权利要求6或7所述的变产变压气井产量经验递减曲线分析方法,其特征在于,根据所述推算的物质平衡时间计算相应的实际生产时间。
9.根据权利要求8所述的变产变压气井产量经验递减曲线分析方法,其特征在于,当所述实际生产时间小于或等于所述预测的最大时间范围,则进行下一时步的指标预测。
10.根据权利要求9所述的变产变压气井产量经验递减曲线分析方法,其特征在于,当所述实际生产时间大于所述预测的最大时间范围,则计算结束,并输出页岩气井指标预测结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710218995.2A CN108694254B (zh) | 2017-04-06 | 2017-04-06 | 一种变产变压生产气井产量经验递减曲线分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710218995.2A CN108694254B (zh) | 2017-04-06 | 2017-04-06 | 一种变产变压生产气井产量经验递减曲线分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108694254A true CN108694254A (zh) | 2018-10-23 |
CN108694254B CN108694254B (zh) | 2021-10-08 |
Family
ID=63842068
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710218995.2A Active CN108694254B (zh) | 2017-04-06 | 2017-04-06 | 一种变产变压生产气井产量经验递减曲线分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108694254B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109681184A (zh) * | 2018-12-05 | 2019-04-26 | 中石化重庆涪陵页岩气勘探开发有限公司 | 页岩气井生产动态预测方法及设备 |
CN110633848A (zh) * | 2019-09-03 | 2019-12-31 | 中国石油天然气股份有限公司 | 煤层气井稳产时间的预测方法及装置 |
CN111506865A (zh) * | 2020-07-01 | 2020-08-07 | 西南石油大学 | 一种页岩气井生产拟合与预测方法 |
CN111764892A (zh) * | 2020-07-03 | 2020-10-13 | 中国石油天然气集团有限公司 | 页岩气井变压变产阶段生产制度适应性判别方法 |
CN111852466A (zh) * | 2020-08-28 | 2020-10-30 | 四川长宁天然气开发有限责任公司 | 一种用于页岩气井规模配产与管网运行优化的方法 |
CN113236207A (zh) * | 2021-07-13 | 2021-08-10 | 西南石油大学 | 一种强非均质性储层中产水气井的定产量递减预测方法 |
CN113323656A (zh) * | 2021-06-16 | 2021-08-31 | 中海石油(中国)有限公司 | 封闭凝析气藏的开发指标预测方法及计算机可读存储介质 |
CN113550744A (zh) * | 2020-04-26 | 2021-10-26 | 中国石油天然气股份有限公司 | 一种页岩气排水采气效果评价方法 |
CN115199240A (zh) * | 2022-08-25 | 2022-10-18 | 西南石油大学 | 一种页岩气井产量预测方法、装置及存储介质 |
CN116085685A (zh) * | 2023-03-15 | 2023-05-09 | 上海叁零肆零科技有限公司 | 保证用气高峰时段天然气稳定供气的方法和系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130346040A1 (en) * | 2012-06-20 | 2013-12-26 | Schlumberger Technology Corporation | Shale gas production forecasting |
CN104500030A (zh) * | 2014-11-21 | 2015-04-08 | 中国石油天然气股份有限公司 | 超高压气藏生产动态异常数据诊断及修正方法 |
CN106056252A (zh) * | 2016-06-13 | 2016-10-26 | 西南石油大学 | 一种页岩气田井组划分三级优化方法 |
CN106481332A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 用于确定页岩气多段压裂水平井内外区动态储量的方法 |
CN106484933A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 一种用于确定页岩气井井控动态储量的方法及系统 |
-
2017
- 2017-04-06 CN CN201710218995.2A patent/CN108694254B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130346040A1 (en) * | 2012-06-20 | 2013-12-26 | Schlumberger Technology Corporation | Shale gas production forecasting |
CN104500030A (zh) * | 2014-11-21 | 2015-04-08 | 中国石油天然气股份有限公司 | 超高压气藏生产动态异常数据诊断及修正方法 |
CN106481332A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 用于确定页岩气多段压裂水平井内外区动态储量的方法 |
CN106484933A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 一种用于确定页岩气井井控动态储量的方法及系统 |
CN106056252A (zh) * | 2016-06-13 | 2016-10-26 | 西南石油大学 | 一种页岩气田井组划分三级优化方法 |
Non-Patent Citations (2)
Title |
---|
方圆: "JY页岩气藏单井产量递减规律研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑(月刊 )》 * |
王董东 等: "页岩气藏可采储量计算新方法", 《大庆石油地质与开发》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109681184A (zh) * | 2018-12-05 | 2019-04-26 | 中石化重庆涪陵页岩气勘探开发有限公司 | 页岩气井生产动态预测方法及设备 |
CN110633848A (zh) * | 2019-09-03 | 2019-12-31 | 中国石油天然气股份有限公司 | 煤层气井稳产时间的预测方法及装置 |
CN110633848B (zh) * | 2019-09-03 | 2022-03-29 | 中国石油天然气股份有限公司 | 煤层气井稳产时间的预测方法及装置 |
CN113550744A (zh) * | 2020-04-26 | 2021-10-26 | 中国石油天然气股份有限公司 | 一种页岩气排水采气效果评价方法 |
CN113550744B (zh) * | 2020-04-26 | 2023-08-22 | 中国石油天然气股份有限公司 | 一种页岩气排水采气效果评价方法 |
CN111506865A (zh) * | 2020-07-01 | 2020-08-07 | 西南石油大学 | 一种页岩气井生产拟合与预测方法 |
CN111764892B (zh) * | 2020-07-03 | 2023-05-12 | 中国石油天然气集团有限公司 | 页岩气井变压变产阶段生产制度适应性判别方法 |
CN111764892A (zh) * | 2020-07-03 | 2020-10-13 | 中国石油天然气集团有限公司 | 页岩气井变压变产阶段生产制度适应性判别方法 |
CN111852466A (zh) * | 2020-08-28 | 2020-10-30 | 四川长宁天然气开发有限责任公司 | 一种用于页岩气井规模配产与管网运行优化的方法 |
CN113323656A (zh) * | 2021-06-16 | 2021-08-31 | 中海石油(中国)有限公司 | 封闭凝析气藏的开发指标预测方法及计算机可读存储介质 |
CN113236207B (zh) * | 2021-07-13 | 2021-09-10 | 西南石油大学 | 一种强非均质性储层中产水气井的定产量递减预测方法 |
CN113236207A (zh) * | 2021-07-13 | 2021-08-10 | 西南石油大学 | 一种强非均质性储层中产水气井的定产量递减预测方法 |
CN115199240A (zh) * | 2022-08-25 | 2022-10-18 | 西南石油大学 | 一种页岩气井产量预测方法、装置及存储介质 |
CN115199240B (zh) * | 2022-08-25 | 2023-05-19 | 西南石油大学 | 一种页岩气井产量预测方法、装置及存储介质 |
CN116085685A (zh) * | 2023-03-15 | 2023-05-09 | 上海叁零肆零科技有限公司 | 保证用气高峰时段天然气稳定供气的方法和系统 |
CN116085685B (zh) * | 2023-03-15 | 2023-06-23 | 上海叁零肆零科技有限公司 | 保证用气高峰时段天然气稳定供气的方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108694254B (zh) | 2021-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108694254A (zh) | 一种变产变压生产气井产量经验递减曲线分析方法 | |
CN109594968B (zh) | 一种页岩气多段压裂水平井压后裂缝参数评价方法及系统 | |
CN105569646B (zh) | 一种油气井技术可采储量预测方法 | |
CN107563899B (zh) | 油气井产能预测方法及装置 | |
CN1966934B (zh) | 一种随钻预测钻头底下地层坍塌压力和破裂压力的方法 | |
CN105134196B (zh) | 一种缝洞型碳酸盐岩凝析气井开发指标分析方法及装置 | |
CN107608940B (zh) | 一种油井间抽周期确定方法 | |
CN104504604B (zh) | 一种定性气井井筒积液的方法 | |
CN111794740B (zh) | 一种适用在缝洞型碳酸盐岩油藏动态储量计算的方法 | |
CN101253402A (zh) | 评价欠饱和煤层气储层的方法 | |
CN106194154B (zh) | 一种非常规油气藏中长期产能预测方法 | |
CN105178943A (zh) | 一种实时校正井筒压力的方法 | |
CN106127599A (zh) | 用于致密气藏在集气站生产模式下的气井产量劈分方法 | |
Liu et al. | Dynamic behavior of coalbed methane flow along the annulus of single-phase production | |
CN102587887B (zh) | 一种气井井况预测方法 | |
CN102425386B (zh) | 一种符合幂律模式的钻井液流变参数控制方法 | |
CN115587674B (zh) | 油藏改建储气库扩容达产过程气井动态产能预测方法 | |
CN105672997A (zh) | 钻井液地层漏失量监测方法 | |
CN109815543A (zh) | 计算气田动态储量的方法 | |
Ibrahim et al. | A new look at reserves estimation of unconventional gas reservoirs | |
CN106150489A (zh) | 一种地层承压能力动态测试方法及固井方法 | |
CN106761680A (zh) | 一种化学降粘辅助螺杆泵举升稠油工艺的判断方法 | |
CN107605474B (zh) | 一种随钻预测气层产量的方法及装置 | |
CN116384267A (zh) | 致密储层压裂水平井的最终可采储量的确定方法及设备 | |
Toczek et al. | Total Specific Energy (TSE) as a new indicator of drilling efficiency |
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 |