CN106127604A - 一种动态储量计算方法及装置 - Google Patents
一种动态储量计算方法及装置 Download PDFInfo
- Publication number
- CN106127604A CN106127604A CN201610436211.9A CN201610436211A CN106127604A CN 106127604 A CN106127604 A CN 106127604A CN 201610436211 A CN201610436211 A CN 201610436211A CN 106127604 A CN106127604 A CN 106127604A
- Authority
- CN
- China
- Prior art keywords
- pressure drop
- gas production
- target area
- unit
- drop gas
- 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
- 238000000205 computational method Methods 0.000 title claims abstract description 12
- 238000004519 manufacturing process Methods 0.000 claims abstract description 273
- 238000000034 method Methods 0.000 claims abstract description 22
- 230000001186 cumulative effect Effects 0.000 claims description 30
- 230000035699 permeability Effects 0.000 claims description 23
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 239000011148 porous material Substances 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 3
- 238000012360 testing method Methods 0.000 abstract description 3
- 230000007246 mechanism Effects 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 241001269238 Data Species 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000000630 rising effect Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000007788 liquid Substances 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
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000004836 empirical method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000003475 lamination Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- 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
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
- Filling Or Discharging Of Gas Storage Vessels (AREA)
Abstract
本发明涉及一种动态储量计算方法及装置,其中,方法包括:确定目标区域的单位压降采气量公式;根据单位压降采气量公式,利用单位压降采气量与单位套压降采气量之间的相互关系,绘制单位套压降采气量曲线;利用单位套压降采气量曲线判断目标区域的储量是否已达到整体动用阶段;分已达到整体动用阶段的目标区域的动态储量、未达到整体动用阶段的目标区域的动态储量这两种情况进行计算。本技术方案根据单位套压降采气量曲线计算的动态储量实践检验效果较好,误差在可接受范围内,尤其对于生产时间较长的气田,结果更加可靠。
Description
技术领域
本发明涉及石油勘探技术领域,特别涉及一种动态储量计算方法及装置。
背景技术
低渗气田资源潜力巨大,已成为我国一种主要的目标区域类型,天然气产能增长的40%以上来自于低渗气田。由于低渗气田储层具有微观孔隙结构复杂、非均质性较强、有效砂体连通性差等特点,其渗流特征与常规中高渗目标区域有明显不同,体现为一般需要较长时间才能达到拟稳定渗流状态,因而气田、气井的动态储量计算比较困难,尤其在气田开发早期,由于时间较短,尚未达到泄流边界,气田整体储量尚未完全动用,气田早期体现的动态储量往往偏小。
目前动态储量的计算方法主要有3种,一种是物质平衡法,主要依据累计产量与地层压力数据,该方法局限在于需要准确的平均地层压力数据和较长的生产时间,对于大面积低渗透目标区域适用性较差;二是弹性二相法等利用试井数据的相关方法,该类方法主要适用于小型的有界封闭弹性气驱目标区域;三是产量累积法、递减曲线法等根据实践所总结的经验法及各种预测方法,一般利用生产数据等动态资料,这类方法存在经验公式的取值准确性,典型曲线匹配多解性等问题。
总体来说,目前对低渗气田、气井的动态储量均要求有一定生产年限,使气田、气井体现一定动态特征的情况下,对其进行动态拟合和储量求取,另外,对于地层压力获取困难,整体处于储量动用早期的大规模低渗透气田,进行动态储量计算存在较大困难。
发明内容
为解决现有技术的问题,本发明提出一种动态储量计算方法及装置,根据单位套压降采气量曲线计算的动态储量实践检验效果较好,误差在可接受范围内,尤其对于生产时间较长的气田,结果更加可靠。
为实现上述目的,本发明提供了一种动态储量计算方法,包括:
确定目标区域的单位压降采气量公式;
根据所述单位压降采气量公式,利用单位压降采气量与单位套压降采气量之间的相互关系,绘制单位套压降采气量曲线;
利用所述单位套压降采气量曲线判断目标区域的储量是否已达到整体动用阶段;
当目标区域的储量达到整体动用阶段时,利用在地面标准条件下单位套压降采气量曲线达到稳定时目标区域的累积产气量、单位套压降采气量曲线达到稳定时目标区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时目标区域的视井口套压确定目标区域的动态储量。
优选地,还包括:
当目标区域的储量未达到整体动用阶段时,利用单位套压降采气量曲线预测目标区域的最大单位套压降采气量;
选择一参照区域,利用在地面标准条件下单位套压降采气量曲线达到稳定时参照区域的累积产气量、单位套压降采气量曲线达到稳定时参照区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时参照区域的视井口套压获得所述参照区域的动态储量;其中,所述参照区域为目标区域内储量达到整体动用阶段的区域或与所述目标区域相邻的储量达到整体动用阶段的区域;
利用参照区域的动态储量、参照区域的最大单位套压降采气量以及目标区域的最大单位套压降采气量确定目标区域的动态储量。
优选地,所述目标区域为低渗透气田区域或低渗透气井区域。
优选地,所述目标区域的单位压降采气量的公式为:
其中,Ep为单位压降采气量;Gp为目标区域在地面标准条件下的累计产气量;ψ(pi)为原始地层条件下的目标区域视地层压力;pi为目标区域原始地层压力;ψ(p)为目前地层压力下的目标区域视地层压力;p为目标区域目前地层压力;G为目标区域在地面标准条件下的原始地质储量;Ω为目标区域的地下连通孔隙体积;Bgi为原始地层条件下的气体体积系数。
优选地,所述单位压降采气量与单位套压降采气量之间的相互关系为:目标区域的单位套压降采气量与目标区域的视井口套压降的乘积等于目标区域的单位压降采气量与目标区域的视地层压降的乘积。
优选地,所述单位压降采气量与单位套压降采气量之间的相互关系为:在定产条件下,目标区域的视地层压降等于视井底流压降等于视井口套压降。
优选地,所述参照区域的动态储量利用在地面标准条件下单位套压降采气量曲线达到稳定时参照区域的累积产气量、单位套压降采气量曲线达到稳定时参照区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时参照区域的气田视井口套压获得。
对应地,为实现上述目的,本发明还提供了一种动态储量计算装置,包括:
单位压降采气量公式确定单元,用于确定目标区域的单位压降采气量公式;
绘制单元,用于根据所述单位压降采气量公式,利用单位压降采气量与单位套压降采气量之间的相互关系,绘制单位套压降采气量曲线;
判断单元,用于利用所述单位套压降采气量曲线判断目标区域的储量是否已达到整体动用阶段;
计算单元,用于当目标区域的储量达到整体动用阶段时,利用在地面标准条件下单位套压降采气量曲线达到稳定时目标区域的累积产气量、单位套压降采气量曲线达到稳定时目标区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时目标区域的视井口套压确定目标区域的动态储量。
优选地,所述计算单元还用于当目标区域的储量未达到整体动用阶段时,利用单位套压降采气量曲线预测目标区域的最大单位套压降采气量;选择一参照区域,利用在地面标准条件下单位套压降采气量曲线达到稳定时参照区域的累积产气量、单位套压降采气量曲线达到稳定时参照区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时参照区域的视井口套压获得所述参照区域的动态储量;利用参照区域的动态储量、参照区域的最大单位套压降采气量以及目标区域的最大单位套压降采气量确定目标区域的动态储量;其中,所述参照区域为目标区域内储量达到整体动用阶段的区域或与所述目标区域相邻的储量达到整体动用阶段的区域。
优选地,所述目标区域为低渗透气田区域或低渗透气井区域。
优选地,所述单位压降采气量公式确定单元确定的目标区域的单位压降采气量公式为:
其中,Ep为单位压降采气量;Gp为目标区域在地面标准条件下的累计产气量;ψ(pi)为原始地层条件下的目标区域视地层压力;pi为目标区域原始地层压力;ψ(p)为目前地层压力下的目标区域视地层压力;p为目标区域目前地层压力;G为目标区域在地面标准条件下的原始地质储量;Ω为目标区域的地下连通孔隙体积;Bgi为原始地层条件下的气体体积系数。
优选地,所述绘制单元涉及的单位压降采气量与单位套压降采气量之间的相互关系为:目标区域的单位套压降采气量与目标区域的视井口套压降的乘积等于目标区域的单位压降采气量与目标区域的视地层压降的乘积。
优选地,所述单位压降采气量与单位套压降采气量之间的相互关系为:在定产条件下,目标区域的视地层压降等于目标区域的视井底流压降等于目标区域的视井口套压降。
优选地,所述计算单元利用在地面标准条件下单位套压降采气量曲线达到稳定时参照区域的累积产气量、单位套压降采气量曲线达到稳定时参照区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时参照区域的视井口套压获得参照区域的动态储量。
上述技术方案具有如下有益效果:
1、降低了对数据的要求,避免了对精确的全气田地层平均压力的需求,尤其是对于已达到整体动用的气田、气井,从理论上支撑了使用套压数据代替地层压力数据的合理性,可直接使用套压等生产数据进行动态储量的精确计算;
2、对于尚未达到整体动用阶段、地层压力测取困难的气田、气井,由于其处于相对较早的生产阶段,则可以近似认为属于定产生产,即地层压降近似等于套压降,使用单位套压降采气量近似代替单位压降采气量,通过与气田内部或邻区已达到整体动用阶段的气井进行类比计算,计算该气田或气井的动态储量;
3、对于处于储量动用早期的大规模低渗气田,常规手段拟合和预测动态储量具有较大困难,本发明基于单位压降采气量图版,提供了一种动态储量计算方法,使研究人员可以在较早阶段开始套合曲线,并预测最大单位压降采气量或单位套压降采气量,避免了对较长生产时间的要求;
4、计算误差较小,提出了精确的动态储量计算公式,对于达到整体动用的气田,可由单位套压降采气量直接计算;对于尚未整体动用的气田,采用套合单位压降采气量曲线的方法,预测时间推前到整体动用时的节点,而非其他方法的气田压力为零的时间节点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提出一种动态储量计算方法流程图;
图2为本发明提出一种动态储量计算装置框图;
图3为本实施例的机理模型示意图;
图4为本实施例的机理模型单位套压降采气量变化曲线图;
图5为本实施例的低渗气田单位套压降采气量曲线图;
图6为拟稳态压力分布图;
图7为本实施例的苏里格气田单位套压降采气量变化曲线图;
图8为本实施例的苏38-16-5单位套压降采气量变化曲线图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本技术方案的工作原理为:本技术方案通过大量实例分析和机理模拟,系统地分析了低渗气田全生命周期单位压降采气量的多阶段变化规律,编制了单位套压降采气量曲线,并对处于不同储量动用阶段的目标区域进行动态储量计算。对于储量已全部处于整体动用阶段的气田、气井来说,可通过生产数据精确计算动态储量;对于储量未全部处于整体动用阶段的气田、气井来说,可根据单位套压降采气量曲线进行单位压降采气量预测,然后通过与其内部气井或周围气井的类比,求取其动态储量。本技术方案可通过单位套压降采气量曲线预测,较好地估算大规模低渗透气田或气井的动态储量,而这类气田或气井由于整体动用需要很长时间,在开发早期通过常规计算方法往往不能通过拟合计算其动态储量。
基于上述工作原理,本发明提出一种动态储量计算方法,如图1所示。包括:
步骤101):确定目标区域的单位压降采气量公式;
步骤102):根据所述单位压降采气量公式,利用单位压降采气量与单位套压降采气量之间的相互关系,绘制单位套压降采气量曲线;
步骤103):利用所述单位套压降采气量曲线判断目标区域的储量是否已达到整体动用阶段;
步骤104):当目标区域的储量达到整体动用阶段时,利用在地面标准条件下单位套压降采气量曲线达到稳定时目标区域的累积产气量、单位套压降采气量曲线达到稳定时目标区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时目标区域的视井口套压确定目标区域的动态储量。
优选地,还包括:
当目标区域的储量未达到整体动用阶段时,利用单位套压降采气量曲线预测目标区域的最大单位套压降采气量;
选择一参照区域,利用在地面标准条件下单位套压降采气量曲线达到稳定时参照区域的累积产气量、单位套压降采气量曲线达到稳定时参照区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时参照区域的视井口套压获得所述参照区域的动态储量;其中,所述参照区域为目标区域内储量达到整体动用阶段的区域或与所述目标区域相邻的储量达到整体动用阶段的区域;
利用参照区域的动态储量、参照区域的最大单位套压降采气量以及目标区域的最大单位套压降采气量确定目标区域的动态储量。
基于上述工作原理,本发明提出一种动态储量计算装置,如图2所示。包括:
单位压降采气量公式确定单元201,用于确定目标区域的单位压降采气量公式;
绘制单元202,用于根据所述单位压降采气量公式,利用单位压降采气量与单位套压降采气量之间的相互关系,绘制单位套压降采气量曲线;
判断单元203,用于利用所述单位套压降采气量曲线判断目标区域的储量是否已达到整体动用阶段;
计算单元204,用于当目标区域的储量达到整体动用阶段时,利用在地面标准条件下单位套压降采气量曲线达到稳定时目标区域的累积产气量、单位套压降采气量曲线达到稳定时目标区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时目标区域的视井口套压确定目标区域的动态储量。
优选地,所述计算单元204还用于当目标区域的储量未达到整体动用阶段时,利用单位套压降采气量曲线预测目标区域的最大单位套压降采气量;选择一参照区域,利用在地面标准条件下单位套压降采气量曲线达到稳定时参照区域的累积产气量、单位套压降采气量曲线达到稳定时参照区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时参照区域的视井口套压获得所述参照区域的动态储量;利用参照区域的动态储量、参照区域的最大单位套压降采气量以及目标区域的最大单位套压降采气量确定目标区域的动态储量;其中,所述参照区域为目标区域内储量达到整体动用阶段的区域或与所述目标区域相邻的储量达到整体动用阶段的区域。
实施例:
步骤A):单位压降采气量的定义及公式推导;
单位压降采气量是指目标区域单位压降的采出气量。
根据定义,可得
其中,ψ(pi)=pi/Zi,ψ(p)=p/Z。式中:Ep为单位压降采气量,108m3/MPa;Gp为目标区域在地面标准条件(0.101MPa,20℃)下的累计产气量,108m3;ψ(Δp)为目标区域的视地层压降,MPa;ψ(pi)为原始地层条件下的目标区域视地层压力,MPa;pi为目标区域原始地层压力,MPa;Zi为原始地层条件下的气体偏差因子;ψ(p)为目前地层压力下的目标区域视地层压力,MPa;p为目标区域目前地层压力,MPa;Z为目前地层条件下的气体偏差因子。
定容目标区域是一种理想化的目标区域,其容积在天然气开采过程中不发生变化。对于定容目标区域,有:
式中:G为目标区域在地面标准条件下的原始地质储量,108m3。
将式(2)代入式(1),得
式中:Ω为目标区域的地下连通孔隙体积,108m3;Bgi为原始地层条件下的气体体积系数,m3/m3。
由式(3)可见,单位压降采气量与采气速度等相关开发技术措施关系不大,主要与气田/气井动用储量的规模有关,即单位压降采气量与气田/气井的动用储量正相关,对于特定的目标区域,原始地层压力一定,Bgi、ψ(pi)均为定值,单位压降采气量主要与目标区域的地下连通孔隙的体积有关。同时应该注意到,在目标区域开采过程中,低渗目标区域单位压降采气量不断增大,体现了目标区域内部连通孔隙的不断增加。
步骤B):绘制单位套压降采气量曲线;
针对地层压力求取有困难的目标区域,进一步研究了目标区域的单位压降采气量与单位套压降采气量之间的相互关系,有:
即Epψ(Δp)=Ecapψ(Δpca) (5)
式中,Ecap为目标区域的单位套压降采气量,108m3/MPa;ψ(Δpca)为目标区域的视井口套压降,MPa。
另外,根据渗流力学原理,在一定范围的排气面积内,气井定产生产一段时间,流动达到拟稳态后,层内各点视压力随时间变化基本相同,不同时间的视压力分布曲线互成一组平行的曲线簇,如图6所示。此时,只要产气量和渗流阻力不变,则视地层压力下降值(简称视地层压降)等于视井底流动压力下降值(简称视井底流压降)。同时,由于视井底流动压力与视井口套压具有稳定的差值关系,因此,视地层压降等于视井口套压下降值(简称视井口套压降),故而,对于定产井,有
ψ(Δp)=ψ(Δpwf)=ψ(Δpca) (6)
式中:pwf为井底流压,pca为井口套压,ψ(Δpwf)为视井底流压降,ψ(Δpca)为目标区域的视井口套压降。
将(6)式代入(1)式,可得
即对于定产气田、气井,可用单位套压降采气量代替单位压降采气量。
目前我国苏里格等大型低渗透气田多采用配产的生产制度,在初期相当一段时间内可以看做定产,利用气井生产过程中累积产量、套压数据,计算单位套压降采气量并绘制成图。通过下文苏里格气田大量实例计算发现,单位套压降采气量求取的动态储量,与实际符合较好,可用于低渗气田、气井的动态储量计算。
步骤C):根据单位套压降采气量曲线图,确定目标区域的储量是否处于整体动用阶段;
根据国外Hugoton、国内靖边、卧龙河嘉五一、相国寺石炭系目标区域等大量实例和机理模型数值模拟计算结果,在低渗区域的全生命周期中,单位套压降采气量曲线具有特定曲线形态,如图4、图5所示。可分为四个典型阶段:缓慢上升阶段、快速上升阶段、平稳上升阶段和稳定阶段,对应了气田储量动用的四个阶段。四个阶段特征描述如下:
1)缓慢上升阶段。在该阶段内,气田动态连通的孔隙缓慢上升,对应生产特征为气井地层压力快速下降,虽然累积产气量快速增加,但是总体上,单位套压降采气量缓慢增加,体现了随着产能建设的进行,不断增加供气源以增加动用储量的过程。
2)快速上升阶段。在该阶段内,气田连通的孔隙体积快速上升,生产上随着钻井继续不断增加,累积产气量不断增加,同时压降速度有所放缓。因此,单位套压降采气量快速上升,该阶段体现了气田储量的快速动用。
3)平稳上升阶段。该阶段气田连通孔隙平稳上升,体现了虽然地层压力下降进一步放缓,但产量增长逐步减小或不增长,累积采气量增长逐步减小。因此,单位套压降采气量呈现平稳上升的特征,体现了气田储量在经历快速上升之后的平稳推进。该阶段的出现,代表两种情况:一是产能建设的放缓,二是对于已结束产能建设的气田,随着地层压力下降,外围更加低渗气源的逐步动用,它们的贡献使整体动用储量平稳上升。
4)整体动用阶段。该阶段气田地下连通孔隙已达到上限,生产特征为气田没有办法通过钻井等手段维持产量,气田产量下降,累积产气量增长缓慢,地层压力下降非常缓慢,气田的单位压降采气量保持不变,气田动态储量已经全部动用,经过一段时间的生产,随着产量趋小和压力趋小,最终气田结束其生命周期。
由于单位压降采气量体现的是气田储量动用程度,显然地,气田单位套压降采气量曲线的变化趋势同单位压降采气量曲线是一致的。根据气田单位套压降采气量曲线,可以确定气田所处的储量动用阶段,判断气田储量是否已经整体动用。随后,可以分以下两种情况,进行气田、气井动态储量的计算。
步骤D)已达到整体动用阶段的气田、气井动态储量的计算
对于单位套压降采气量曲线已趋平稳,储量已达到整体动用的气田、气井,根据动态储量定义,可使用下式计算动态储量:
Gr=Gps+Esψ(Δps)=Gps+Es(ψ(ps)-0)=Gps+Esψ(ps) (8)
式中,Gr为气田动态储量,108m3;Gps为在地面标准条件(0.101MPa,20℃)下,单位套压降采气量曲线达到稳定时的累积产气量,108m3;Es为最大单位压降采气量,即单位压降采气量曲线达到稳定时的数值,108m3/MPa;ψ(Δps)为单位压降采气量曲线达到稳定时的视地层压降,MPa;ψ(ps)为单位压降采气量曲线达到稳定时的视地层压力,MPa。
根据式(5),单位压降采气量与单位套压降采气量满足以下关系式:
Esψ(Δps)=Ecasψ(Δpcas) (9)
式中,Ecas为最大单位套压降采气量,即单位套压降采气量曲线达到稳定时的数值,108m3/MPa,ψ(Δpcas)为单位套压降采气量曲线达到稳定时的视井口套压,MPa。
将式(9)代入式(8),可得:
Gr=Gps+Esψ(Δps)=Gps+Ecasψ(Δpcas) (10)
此式可以精确计算已达到整体动用的低渗气田的动态储量,并且仅需要产量和套压两种生产数据,对于生产时间较长,地层压力测取有困难的气井、气田,是十分方便的。
步骤E)未达到整体动用阶段的气田、气井动态储量的计算
对于尚未到达整体动用阶段的,通过以下步骤计算气田动态储量:
①如图5所示,根据单位套压降采气量曲线形态,预测该目标气田的单位套压降采气量曲线达到稳定时的数值Ecas;
②选取目标气田中生产时间较长、已到达整体动用阶段的气井,或气井周围生产时间较长、已到达整体动用的邻井,读取其趋稳的单位套压降采气量Ewcas;
③根据式(10)计算选择的气井动态储量Gwr;
④通过以下步骤,计算目标气田的动态储量Gr。
对于定容目标区域,根据式(3),有
同时,对于尚未达到整体动用阶段的气田、气井,可以近似认为定产生产阶段,基本符合式(6)、(7)的前提条件,可以用单位套压降采气量代替单位压降采气量,因此,式(11)可变换为:
变换可得
根据式(13),可以计算得到尚未到达整体动用阶段的目标区域的动态储量。
1、机理模型动态储量计算
设计机理模型一个,储量3×108m3,渗透率0.3mD,不考虑压裂,600米井距,共钻9口井,分3年钻井,每年钻3口井,单井日产气8000m3,生产20年,见附表1。
如图3所示,机理模型经过20年的生产,单位压降采气量曲线由平稳上升到全部动用阶段的拐点已经出现,气田已到达全部动用状态,此时,达到稳定时的累积产气量Gps为1.92×108m3,气田最大单位压降采气量Es为0.1079×108m3/MPa,气田单位压降采气量曲线达到稳定时的视地层压力ψ(ps)为10.24MPa,按照式(8),可得机理模型的动态储量
Gr=Gps+Esψ(ps))=1.92+0.1079×10.24=3.0(108m3)
由于机理模型为均质模型,其动态模型理论数值即地质储量3.0×108m3,因此,本发明计算方法准确,误差为0。
表1机理模型生产数据表
2、苏里格气田动态储量计算
苏里格气田为典型的大规模低渗透气田,气田整体探明储量目前为4.2×108m3,且尚在进一步增加,缺乏整体动态储量规模认识;同时该气田地层压力数据较少,平均地层压力求取困难,运用套压数据计算苏里格气田单位压降采气量曲线图,数据信息见附表2。如图7所示,苏里格气田目前单位套压降采气量曲线处于快速上升阶段,尚有很大的开发潜力。
首先,将图7覆于附图5上,拉伸使曲线吻合,则苏里格气田单位套压降采气量稳定时可达756亿方/MPa。
然后,选取2003年10月投产的苏38-16-5井,如图8所示。该井生产时间超过10年,数据信息见表3,单位套压降采气量曲线已到达整体动用阶段,达到稳定时的累积产气量Gps为0.5003×108m3,气井最大单位套压降采气量Ewcas为0.0255×108m3/MPa,气井单位套压降采气量曲线达到稳定时的视井口套压降ψ(Δpcas)为3.4MPa,按照式(10),可得该井的动态储量:
Gwr=Gps+Esψ(Δps)=Gps+Ecasψ(Δpcas)=0.5003+0.0255×3.4=0.587(108m3)
最后,根据式(13),计算苏里格气田动态储量:
对比采用相同数据用累积产量法计算气田动态储量为19500亿方,该发明与压降法误差10.7%,对于一个储量动用尚处于早期快速动用阶段的大型低渗气田,误差在可接受范围之内,该发明计算符合较好,结果较为可靠。
附表2苏里格气田生产运行简表
附表3苏38-16-5气井生产运行简表
本发明通过大量实例分析和机理模拟,系统地分析了低渗气田全生命周期单位压降采气量的多阶段变化规律,编制了低渗目标区域单位压降采气量曲线变化图版,并对处于不同储量动用阶段的目标区域进行动态储量计算。该发明对于储量已全部动用的气田、气井,可通过生产数据精确计算动态储量;对于储量未全部动用的气田、气井,可根据图版进行单位压降采气量预测,然后通过与其内部气井或周围气井的类比,求取其动态储量。该发明可通过图版预测,较好地估算大规模低渗透气田的动态储量,而这类气田由于整体动用需要很长时间,在开发早期通过常规计算方法往往不能通过拟合计算其动态储量。
以上具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (14)
1.一种动态储量计算方法,其特征在于,包括:
确定目标区域的单位压降采气量公式;
根据所述单位压降采气量公式,利用单位压降采气量与单位套压降采气量之间的相互关系,绘制单位套压降采气量曲线;
利用所述单位套压降采气量曲线判断目标区域的储量是否已达到整体动用阶段;
当目标区域的储量达到整体动用阶段时,利用在地面标准条件下单位套压降采气量曲线达到稳定时目标区域的累积产气量、单位套压降采气量曲线达到稳定时目标区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时目标区域的视井口套压确定目标区域的动态储量。
2.如权利要求1所述的方法,其特征在于,还包括:
当目标区域的储量未达到整体动用阶段时,利用单位套压降采气量曲线预测目标区域的最大单位套压降采气量;
选择一参照区域,利用在地面标准条件下单位套压降采气量曲线达到稳定时参照区域的累积产气量、单位套压降采气量曲线达到稳定时参照区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时参照区域的视井口套压获得所述参照区域的动态储量;其中,所述参照区域为目标区域内储量达到整体动用阶段的区域或与所述目标区域相邻的储量达到整体动用阶段的区域;
利用参照区域的动态储量、参照区域的最大单位套压降采气量以及目标区域的最大单位套压降采气量确定目标区域的动态储量。
3.如权利要求1或2所述的方法,其特征在于,所述目标区域为低渗透气田区域或低渗透气井区域。
4.如权利要求1~3任一权利要求所述的方法,其特征在于,所述目标区域的单位压降采气量公式为:
其中,Ep为单位压降采气量;Gp为目标区域在地面标准条件下的累计产气量;ψ(pi)为原始地层条件下的目标区域视地层压力;pi为目标区域原始地层压力;ψ(p)为目前地层压力下的目标区域视地层压力;p为目标区域目前地层压力;G为目标区域在地面标准条件下的原始地质储量;Ω为目标区域的地下连通孔隙体积;Bgi为原始地层条件下的气体体积系数。
5.如权利要求1~3任一权利要求所述的方法,其特征在于,所述单位压降采气量与单位套压降采气量之间的相互关系为:目标区域的单位套压降采气量与目标区域的视井口套压降的乘积等于目标区域的单位压降采气量与目标区域的视地层压降的乘积。
6.如权利要求1~3任一权利要求所述的方法,其特征在于,所述单位压降采气量与单位套压降采气量之间的相互关系为:在定产条件下,目标区域的视地层压降等于视井底流压降等于视井口套压降。
7.如权利要求1~3任一权利要求所述的方法,其特征在于,所述参照区域的动态储量利用在地面标准条件下单位套压降采气量曲线达到稳定时参照区域的累积产气量、单位套压降采气量曲线达到稳定时参照区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时参照区域的视井口套压获得。
8.一种动态储量计算装置,其特征在于,包括:
单位压降采气量公式确定单元,用于确定目标区域的单位压降采气量公式;
绘制单元,用于根据所述单位压降采气量公式,利用单位压降采气量与单位套压降采气量之间的相互关系,绘制单位套压降采气量曲线;
判断单元,用于利用所述单位套压降采气量曲线判断目标区域的储量是否已达到整体动用阶段;
计算单元,用于当目标区域的储量达到整体动用阶段时,利用在地面标准条件下单位套压降采气量曲线达到稳定时目标区域的累积产气量、单位套压降采气量曲线达到稳定时目标区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时目标区域的视井口套压确定目标区域的动态储量。
9.如权利要求8所述的装置,其特征在于,所述计算单元还用于当目标区域的储量未达到整体动用阶段时,利用单位套压降采气量曲线预测目标区域的最大单位套压降采气量;选择一参照区域,利用在地面标准条件下单位套压降采气量曲线达到稳定时参照区域的累积产气量、单位套压降采气量曲线达到稳定时参照区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时参照区域的视井口套压获得所述参照区域的动态储量;利用参照区域的动态储量、参照区域的最大单位套压降采气量以及目标区域的最大单位套压降采气量确定目标区域的动态储量;其中,所述参照区域为目标区域内储量达到整体动用阶段的区域或与所述目标区域相邻的储量达到整体动用阶段的区域。
10.如权利要求8或9所述的装置,其特征在于,所述目标区域为低渗透气田区域或低渗透气井区域。
11.如权利要求8~10任一权利要求所述的装置,其特征在于,所述单位压降采气量公式确定单元确定的目标区域的单位压降采气量公式为:
其中,Ep为单位压降采气量;Gp为目标区域在地面标准条件下的累计产气量;ψ(pi)为原始地层条件下的目标区域视地层压力;pi为目标区域原始地层压力;ψ(p)为目前地层压力下的目标区域视地层压力;p为目标区域目前地层压力;G为目标区域在地面标准条件下的原始地质储量;Ω为目标区域的地下连通孔隙体积;Bgi为原始地层条件下的气体体积系数。
12.如权利要求8~10任一权利要求所述的装置,其特征在于,所述绘制单元涉及的单位压降采气量与单位套压降采气量之间的相互关系为:目标区域的单位套压降采气量与目标区域的视井口套压降的乘积等于目标区域的单位压降采气量与目标区域的视地层压降的乘积。
13.如权利要求8~10任一权利要求所述的装置,其特征在于,所述单位压降采气量与单位套压降采气量之间的相互关系为:在定产条件下,目标区域的视地层压降等于视井底流压降等于视井口套压降。
14.如权利要求8~10任一权利要求所述的装置,其特征在于,所述计算单元利用在地面标准条件下单位套压降采气量曲线达到稳定时参照区域的累积产气量、单位套压降采气量曲线达到稳定时参照区域的最大单位套压降采气量以及单位套压降采气量曲线达到稳定时参照区域的视井口套压获得参照区域的动态储量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610436211.9A CN106127604B (zh) | 2016-06-17 | 2016-06-17 | 一种动态储量计算方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610436211.9A CN106127604B (zh) | 2016-06-17 | 2016-06-17 | 一种动态储量计算方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106127604A true CN106127604A (zh) | 2016-11-16 |
CN106127604B CN106127604B (zh) | 2019-05-07 |
Family
ID=57470954
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610436211.9A Active CN106127604B (zh) | 2016-06-17 | 2016-06-17 | 一种动态储量计算方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106127604B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108612525A (zh) * | 2018-04-19 | 2018-10-02 | 重庆科技学院 | 一种气藏动态储量计算方法 |
CN110067555A (zh) * | 2018-01-24 | 2019-07-30 | 中国石油天然气股份有限公司 | 碳酸盐岩油井的最小动态储量的确定方法和装置 |
CN112049624A (zh) * | 2019-06-06 | 2020-12-08 | 中国石油天然气股份有限公司 | 油井动态储量的预测方法、装置、设备及存储介质 |
CN113515833A (zh) * | 2020-04-09 | 2021-10-19 | 中国石油化工股份有限公司 | 一种控制气井压降速度的配产设计方法 |
CN114382465A (zh) * | 2020-10-21 | 2022-04-22 | 中国石油天然气股份有限公司 | 确定气井的稳产年限的方法、装置、终端和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104832156A (zh) * | 2015-05-05 | 2015-08-12 | 中国石油天然气股份有限公司 | 一种预估气井产量的方法 |
CN105243182A (zh) * | 2015-09-08 | 2016-01-13 | 中国石油天然气股份有限公司 | 致密油压裂水平井的动态储量计算方法 |
CN105545275A (zh) * | 2015-12-31 | 2016-05-04 | 中国石油天然气股份有限公司 | 气田气井部署方法和装置 |
CN105649616A (zh) * | 2015-12-29 | 2016-06-08 | 中国石油天然气股份有限公司 | 一种低渗气藏井下节流条件下气井动储量评价方法 |
-
2016
- 2016-06-17 CN CN201610436211.9A patent/CN106127604B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104832156A (zh) * | 2015-05-05 | 2015-08-12 | 中国石油天然气股份有限公司 | 一种预估气井产量的方法 |
CN105243182A (zh) * | 2015-09-08 | 2016-01-13 | 中国石油天然气股份有限公司 | 致密油压裂水平井的动态储量计算方法 |
CN105649616A (zh) * | 2015-12-29 | 2016-06-08 | 中国石油天然气股份有限公司 | 一种低渗气藏井下节流条件下气井动储量评价方法 |
CN105545275A (zh) * | 2015-12-31 | 2016-05-04 | 中国石油天然气股份有限公司 | 气田气井部署方法和装置 |
Non-Patent Citations (1)
Title |
---|
季丽丹 等: ""气藏开发阶段划分新方法"", 《断块油气田》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110067555A (zh) * | 2018-01-24 | 2019-07-30 | 中国石油天然气股份有限公司 | 碳酸盐岩油井的最小动态储量的确定方法和装置 |
CN108612525A (zh) * | 2018-04-19 | 2018-10-02 | 重庆科技学院 | 一种气藏动态储量计算方法 |
CN108612525B (zh) * | 2018-04-19 | 2021-05-28 | 重庆科技学院 | 一种气藏动态储量计算方法 |
CN112049624A (zh) * | 2019-06-06 | 2020-12-08 | 中国石油天然气股份有限公司 | 油井动态储量的预测方法、装置、设备及存储介质 |
CN112049624B (zh) * | 2019-06-06 | 2024-04-30 | 中国石油天然气股份有限公司 | 油井动态储量的预测方法、装置、设备及存储介质 |
CN113515833A (zh) * | 2020-04-09 | 2021-10-19 | 中国石油化工股份有限公司 | 一种控制气井压降速度的配产设计方法 |
CN114382465A (zh) * | 2020-10-21 | 2022-04-22 | 中国石油天然气股份有限公司 | 确定气井的稳产年限的方法、装置、终端和存储介质 |
CN114382465B (zh) * | 2020-10-21 | 2024-04-30 | 中国石油天然气股份有限公司 | 确定气井的稳产年限的方法、装置、终端和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106127604B (zh) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106127604A (zh) | 一种动态储量计算方法及装置 | |
CN110334431B (zh) | 一种低渗透致密气藏单井控制储量计算及剩余气分析方法 | |
CN105243182B (zh) | 致密油压裂水平井的动态储量计算方法 | |
CN105587298B (zh) | 多流管模式的水驱油井含水率反演方法 | |
CN104632187B (zh) | 一种确定产水煤层气井动态储量的方法 | |
CN106481332B (zh) | 用于确定页岩气多段压裂水平井内外区动态储量的方法 | |
CN104018829B (zh) | 一种利用煤层气井生产数据测量气水相渗曲线的方法 | |
CN106484933A (zh) | 一种用于确定页岩气井井控动态储量的方法及系统 | |
CN105631184B (zh) | 一种油气井产量递减分析方法和系统 | |
CN108229051B (zh) | 一种预测油藏空气泡沫驱采收率的方法 | |
CN107038268A (zh) | 一种确定非均质储层五点井网水驱波及系数的方法 | |
CN104453834A (zh) | 一种井组注采关系优化调整方法 | |
CN104695950A (zh) | 火山岩油藏产能预测方法 | |
CN106545336A (zh) | 考虑致密气藏渗流机理的产能计算方法 | |
CN104612635A (zh) | 煤层群瓦斯联合抽采达标预评判方法 | |
CN109242364B (zh) | 一种高温高压气井模拟井筒体积置换产能评价法 | |
CN110929462B (zh) | 一种使用生产数据计算低渗气藏真实压力及气藏储量的方法 | |
CN114427432B (zh) | 一种气藏剩余气开发潜力确定方法 | |
CN106640021B (zh) | 压后放喷参数的计算方法及装置 | |
CN106194164B (zh) | 边底水油藏开发岩心实验模拟方法 | |
CN104975827B (zh) | 预测二氧化碳驱油藏指标的物质平衡方法 | |
CN111794740A (zh) | 一种适用在缝洞型碳酸盐岩油藏动态储量计算的方法 | |
CN107066718A (zh) | 一种四维动态地应力模拟方法 | |
CN110750930A (zh) | 一种基于裂缝连续体模型预测裂缝性储层应力演化的方法 | |
CN107130955B (zh) | 井底流压的确定方法及储集体天然能量的确定方法 |
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 |