CN104504472A - 储气库气水交互区井产能的预测方法及装置 - Google Patents
储气库气水交互区井产能的预测方法及装置 Download PDFInfo
- Publication number
- CN104504472A CN104504472A CN201410809184.6A CN201410809184A CN104504472A CN 104504472 A CN104504472 A CN 104504472A CN 201410809184 A CN201410809184 A CN 201410809184A CN 104504472 A CN104504472 A CN 104504472A
- Authority
- CN
- China
- Prior art keywords
- gas
- well
- air water
- gas storage
- reservoir
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 201
- 238000003860 storage Methods 0.000 title claims abstract description 146
- 238000000034 method Methods 0.000 title claims abstract description 103
- 230000003993 interaction Effects 0.000 title abstract 6
- 238000004519 manufacturing process Methods 0.000 claims abstract description 92
- 238000006073 displacement reaction Methods 0.000 claims abstract description 83
- 230000035699 permeability Effects 0.000 claims abstract description 76
- 230000008569 process Effects 0.000 claims abstract description 68
- 238000011160 research Methods 0.000 claims abstract description 47
- 238000002474 experimental method Methods 0.000 claims abstract description 38
- 238000011161 development Methods 0.000 claims abstract description 20
- 238000012360 testing method Methods 0.000 claims abstract description 20
- 238000013461 design Methods 0.000 claims abstract description 18
- 230000008859 change Effects 0.000 claims abstract description 10
- 230000001603 reducing effect Effects 0.000 claims abstract description 6
- 230000002452 interceptive effect Effects 0.000 claims description 47
- 238000005325 percolation Methods 0.000 claims description 35
- 230000004941 influx Effects 0.000 claims description 33
- 239000004576 sand Substances 0.000 claims description 24
- 230000006870 function Effects 0.000 claims description 20
- 238000009933 burial Methods 0.000 claims description 12
- 238000004088 simulation Methods 0.000 claims description 7
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 238000012407 engineering method Methods 0.000 claims description 6
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 claims description 6
- 238000013508 migration Methods 0.000 claims description 6
- 230000005012 migration Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 abstract description 6
- 239000007789 gas Substances 0.000 description 220
- 238000010586 diagram Methods 0.000 description 11
- 238000004590 computer program Methods 0.000 description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- 230000008901 benefit Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 239000003345 natural gas Substances 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005303 weighing 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
- 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
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种储气库气水交互区井产能的预测方法及装置,其中方法包括:确定气藏净水侵量与压力的函数关系,根据所述函数关系确定在改建储气库运行过程中处于气水交互区的研究对象井;处理所述研究对象井的气藏开发数据或试井资料,得到气井二项式产能方程,确定气水交互驱替降低产能的影响因素;根据所述气水交互驱替降低产能的影响因素,设计储气库建库过程的渗流实验,测试储气库设计寿命每一个循环周期的气相渗透率变化情况;用渗流实验结果修正气井二项式产能方程,得到扣除储气库运行过程中由于气水交互驱替而损失渗流能力的所述研究对象井的产能。本发明可以提高水淹气藏改建储气库的井产能预测的准确性。
Description
技术领域
本发明涉及天然气地下储存、气藏工程、室内实验技术领域,尤其涉及储气库气水交互区井产能的预测方法及装置。
背景技术
储气库是调节天然气供需平衡的有力手段,衡量其调节能力的重要指标是采气能力,采气能力通过储气库井产能实现,准确预测井产能是困扰相关科研人员的技术难题。
目前,预测水淹气藏改建储气库井产能时,大都没有考虑气水交驱替作用对产能的影响。如文献2012年2月第32卷第2期《天然气工业》的文章“相国寺石炭系气藏改建地下储气库运行参数设计”,直接将气藏开发阶段的井产能等同于改建储气库后的井产能;文献2011年2月第32卷第1期《新疆石油地质》的文章“永安油田永21块地下储气库气井产能的确定”,所描述的方法仅适用于储层有一次水侵经历的情况,无法准确描述储气库运行过程中气水多次交互驱替的物理过程。
因此,当前水淹气藏改建储气库的井产能预测方法存在局限性,特别是位于气水交互区储层的井产能预测方法报道较少,需要相对准确的计算方法出现。
发明内容
本发明实施例提供一种储气库气水交互区井产能的预测方法,用以提高水淹气藏改建储气库的井产能预测的准确性,该方法包括:
确定气藏净水侵量与压力的函数关系,根据所述函数关系确定在改建储气库运行过程中处于气水交互区的研究对象井;
处理所述研究对象井的气藏开发数据或试井资料,得到气井二项式产能方程,确定气水交互驱替降低产能的影响因素;
根据所述气水交互驱替降低产能的影响因素,设计储气库建库过程的渗流实验,测试储气库设计寿命每一个循环周期的气相渗透率变化情况;
用渗流实验结果修正气井二项式产能方程,得到扣除储气库运行过程中由于气水交互驱替而损失渗流能力的所述研究对象井的产能。
一个实施例中,确定气藏净水侵量与压力的函数关系,根据所述函数关系确定在改建储气库运行过程中处于气水交互区的研究对象井,包括:
通过气藏工程方法或数值模拟方法得到净水侵量W,与储层平均压力P回归成函数关系:
W=f(P);
代入储气库上限压力Pmax,得到储气库运行至上限压力时的水侵量:
Wmax=f(Pmax);
代入储气库下限压力Pmin,得到储气库运行至下限压力时的水侵量:
Wmin=f(Pmin);
根据研究对象井所处的埋藏深度,计算储气库气水界面运移至该埋藏深度时的压力P井,进而得到储气库运行至压力P井时的水侵量:
W井=f(P井);
若Wmax<W井<Wmin,则确定该井为所述研究对象井。
一个实施例中,处理所述研究对象井的气藏开发数据或试井资料,得到气井二项式产能方程,确定气水交互驱替降低产能的影响因素,包括:
选择气藏开发阶段两个生产稳定数据点:
生产数据点1(Pe1,Pwf1,Q1);
生产数据点2(Pe2,Pwf2,Q2);
其中,Pe1是生产数据点1的地层压力,Pwf1是生产数据点1的井底流压,Q1是生产数据点1的气体流速;Pe2是生产数据点2的地层压力,Pwf2是生产数据点2的井底流压,Q2是生产数据点2的气体流速;
依据所述两个生产稳定数据点,处理得到拟压力气井产能二项式方程:
Pe'-Pwf′=AQ+BQ2;
其中,P'为真实气体的拟压力:
Pe'为地层压力的拟压力,Pwf′为井底流压的拟压力,Q为气体流速,μ为气体粘度,Z为气体压缩因子,A和B是常数,通过所述两个生产稳定数据点联立求解,是渗透率K的函数;
确定气水交互驱替降低产能的影响因素为储层渗透率K。
一个实施例中,根据所述气水交互驱替降低产能的影响因素,设计储气库建库过程的渗流实验,测试储气库设计寿命每一个循环周期的气相渗透率变化情况,包括:
根据所述气水交互驱替降低产能的影响因素,设计气驱水-水驱气实验,描述气水交互区气水互驱替的物理过程,得到气水交互区的端点的气相相对渗透率:
K1/K:第一次气水交互驱替后气相相对渗透率,描述储气库第一周期储层经过一次气水互区造成渗流能力降低后剩余的相对渗透率;
K2/K:第二次气水交互驱替后气相相对渗透率,描述储气库第二周期储层经二次气水互区造成渗流能力降低后剩余的相对渗透率;
K3/K:第三次气水交互驱替后气相相对渗透率,描述储气库第三周期储层经过三次气水互区造成渗流能力降低后剩余的相对渗透率;
直至Kn/K:第n次气水交互驱替后气相相对渗透率,描述储气库第n周期储层经过n次气水互区造成渗流能力降低后剩余的相对渗透率。
一个实施例中,用渗流实验结果修正气井二项式产能方程,得到扣除储气库运行过程中由于气水交互驱替而损失渗流能力的所述研究对象井的产能,包括:
将渗流实验结果代入真实气体的拟压力函数,得到修正的拟压力函数表达式:
将修正后的拟压力表达式代入气井二项式产能方程,得到修正后二项式产能方程:
Pen″-Pewf″=AQ+BQ2;
其中,Pen″为地层压力的修正拟压力,Pewf″为井底流压的修正拟压力。。
本发明实施例还提供一种储气库气水交互区井产能的预测装置,用以提高水淹气藏改建储气库的井产能预测的准确性,该装置包括:
研究对象井确定模块,用于确定气藏净水侵量与压力的函数关系,根据所述函数关系确定在改建储气库运行过程中处于气水交互区的研究对象井;
影响因素确定模块,用于处理所述研究对象井的气藏开发数据或试井资料,得到气井二项式产能方程,确定气水交互驱替降低产能的影响因素;
渗流实验模块,用于根据所述气水交互驱替降低产能的影响因素,设计储气库建库过程的渗流实验,测试储气库设计寿命每一个循环周期的气相渗透率变化情况;
产能预测模块,用于用渗流实验结果修正气井二项式产能方程,得到扣除储气库运行过程中由于气水交互驱替而损失渗流能力的所述研究对象井的产能。
一个实施例中,研究对象井确定模块具体用于:
通过气藏工程方法或数值模拟方法得到净水侵量W,与储层平均压力P回归成函数关系:
W=f(P);
代入储气库上限压力Pmax,得到储气库运行至上限压力时的水侵量:
Wmax=f(Pmax);
代入储气库下限压力Pmin,得到储气库运行至下限压力时的水侵量:
Wmin=f(Pmin);
根据研究对象井所处的埋藏深度,计算储气库气水界面运移至该埋藏深度时的压力P井,进而得到储气库运行至压力P井时的水侵量:
W井=f(P井);
若Wmax<W井<Wmin,则确定该井为所述研究对象井。
一个实施例中,影响因素确定模块具体用于:
选择气藏开发阶段两个生产稳定数据点:
生产数据点1(Pe1,Pwf1,Q1);
生产数据点2(Pe2,Pwf2,Q2);
其中,Pe1是生产数据点1的地层压力,Pwf1是生产数据点1的井底流压,Q1是生产数据点1的气体流速;Pe2是生产数据点2的地层压力,Pwf2是生产数据点2的井底流压,Q2是生产数据点2的气体流速;
依据所述两个生产稳定数据点,处理得到拟压力气井产能二项式方程:
Pe'-Pwf′=AQ+BQ2;
其中,P'为真实气体的拟压力:
Pe'为地层压力的拟压力,Pwf′为井底流压的拟压力,Q为气体流速,μ为气体粘度,Z为气体压缩因子,A和B是常数,通过所述两个生产稳定数据点联立求解,是渗透率K的函数;
确定气水交互驱替降低产能的影响因素为储层渗透率K。
一个实施例中,渗流实验模块具体用于:
根据所述气水交互驱替降低产能的影响因素,设计气驱水-水驱气实验,描述气水交互区气水互驱替的物理过程,得到气水交互区的端点的气相相对渗透率:
K1/K:第一次气水交互驱替后气相相对渗透率,描述储气库第一周期储层经过一次气水互区造成渗流能力降低后剩余的相对渗透率;
K2/K:第二次气水交互驱替后气相相对渗透率,描述储气库第二周期储层经二次气水互区造成渗流能力降低后剩余的相对渗透率;
K3/K:第三次气水交互驱替后气相相对渗透率,描述储气库第三周期储层经过三次气水互区造成渗流能力降低后剩余的相对渗透率;
直至Kn/K:第n次气水交互驱替后气相相对渗透率,描述储气库第n周期储层经过n次气水互区造成渗流能力降低后剩余的相对渗透率。
一个实施例中,产能预测模块具体用于:
将渗流实验结果代入真实气体的拟压力函数,得到修正的拟压力函数表达式:
将修正后的拟压力表达式代入气井二项式产能方程,得到修正后二项式产能方程:
Pen″-Pewf″=AQ+BQ2;
Pen″为地层压力的修正拟压力,Pewf″为井底流压的修正拟压力。
本发明实施例在预测气水交互区井产能时,没有像现有方法直接取气藏开发阶段处理得到的恒定数值,而是扣除由于气水交互驱替而损失的产能,进而给出了储气库寿命内每一个循环周期的产能变化,与现有方法相比,本发明实施例的计算结果准确度大幅度提高,而且随着储气库运行周期的延长,技术优势进越发凸显。本发明实施例的预测方案,将为储气库前期评价提供相对准确的井产能技术参数。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中储气库气水交互区井产能的预测方法的示意图;
图2为本发明实施例中气水交互区位置的示意图;
图3为本发明实施例中储气库气水交互区井产能的预测装置的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
发明人发现,现有针对水淹气藏改建储气库时气水交互区的井产能预测方法存在局限性的原因在于,没有扣除储气库运行过程中气水交互驱替对产能的降低作用,更不能预测储气库逐个周期扣除气水交互驱替影响后的井产能,基于此,针对现有气水交互区井产能预测方法的不足,在本发明实施例中扣除储气库运行过程中由于气水交互驱替而损失的渗流能力,进而得到一种相对准确预测储气库气水交互区井产能的方法。
具体的,本发明实施例依据建库储层中的气水界面,确定储气库井所在的气水交互储层区域,针对气水交互区的特点设计物理模拟实验,利用实验结果描述由于气水交互驱替而损失的渗流能力,最终得到气水交互区井产能的连续变化预测方法。
图1为本发明实施例中储气库气水交互区井产能的预测方法的示意图。如图1所示,本发明实施例中储气库气水交互区井产能的预测方法可以包括:
步骤101、确定气藏净水侵量与压力的函数关系,根据所述函数关系确定在改建储气库运行过程中处于气水交互区的研究对象井;
步骤102、处理所述研究对象井的气藏开发数据或试井资料,得到气井二项式产能方程,确定气水交互驱替降低产能的影响因素;
步骤103、根据所述气水交互驱替降低产能的影响因素,设计储气库建库过程的渗流实验,测试储气库设计寿命每一个循环周期的气相渗透率变化情况;
步骤104、用渗流实验结果修正气井二项式产能方程,得到扣除储气库运行过程中由于气水交互驱替而损失渗流能力的所述研究对象井的产能。
具体实施时,先研究气藏净水侵量-压力的函数关系,确定某井在改建储气库运行过程中处于气水交互区。
实施例中,净水侵量W可以通过气藏工程方法或数值模拟方法得到,与储层平均压力P回归成函数关系:
W=f(P);
代入储气库上限压力Pmax,得到储气库运行至上限压力时的水侵量:
Wmax=f(Pmax);
代入储气库下限压力Pmin,得到储气库运行至下限压力时的水侵量:
Wmin=f(Pmin);
根据研究对象井所处的埋藏深度,计算储气库气水界面运移至该埋藏深度时的压力P井,进而得到储气库运行至压力P井时的水侵量:
W井=f(P井);
如果Wmax<W井<Wmin,则该井适于本发明实施例中来预测产能,确定该井为所述研究对象井。
在确定研究对象之后,处理气藏开发数据或试井资料,得到气井二项式产能方程,明确气水交互驱替降低产能的影响因素。实施例中,可以优选气藏开发阶段两个生产稳定数据点:
生产数据点1(Pe1,Pwf1,Q1);
生产数据点2(Pe2,Pwf2,Q2);
其中,Pe1是生产数据点1的地层压力,Pwf1是生产数据点1的井底流压,Q1是生产数据点1的气体流速;Pe2是生产数据点2的地层压力,Pwf2是生产数据点2的井底流压,Q2是生产数据点2的气体流速;
依据这两个生产稳定数据点,处理得到拟压力气井产能二项式方程:
Pe'-Pwf′=AQ+BQ2;
其中,P'为真实气体的拟压力:
Pe'为地层压力的拟压力,Pwf′为井底流压的拟压力,Q为气体流速,μ为气体粘度,Z为气体压缩因子,A和B是常数,通过所述两个生产稳定数据点联立求解,是渗透率K的函数;
因此,影响产能变化的参数主要为储层渗透率K,需要着重研究储层在改建储气库后的变化规律。
在确定气水交互驱替降低产能的影响因素后,根据所述气水交互驱替降低产能的影响因素,即气水交互驱替的特点,设计储气库建库过程的渗流实验,测试储气库设计寿命n周期中每一个循环周期的气相渗透率变化情况。
图2为本发明实施例中气水交互区位置的示意图,如图2所示,气水交互区,位于储气库上限压力与下限压力对应的气水界面之间的储层,即储气库注满气时为气区、采气末期为水区,在储气库运行过程中处于气水交互驱替状态。
根据气水交互驱替的特点,设计气驱水-水驱气实验,描述气水交互区气水互驱替的物理过程,得到气水交互区的端点的气相相对渗透率:
K1/K:第一次气水交互驱替后气相相对渗透率,描述储气库第一周期储层经过一次气水互区造成渗流能力降低后剩余的相对渗透率;
K2/K:第二次气水交互驱替后气相相对渗透率,描述储气库第二周期储层经二次气水互区造成渗流能力降低后剩余的相对渗透率;
K3/K:第三次气水交互驱替后气相相对渗透率,描述储气库第三周期储层经过三次气水互区造成渗流能力降低后剩余的相对渗透率;
……
Kn/K:第n次气水交互驱替后气相相对渗透率,描述储气库第n周期储层经过n次气水互区造成渗流能力降低后剩余的相对渗透率。
最后,用渗流实验结果修正气井二项式产能方程,进而得到扣除储气库运行过程中由于气水交互驱替而损失渗流能力的井产能。
实施例中,可以将渗流实验结果代入真实气体的拟压力函数,得到修正的拟压力函数表达式:
将修正后的拟压力表达式代入气井二项式产能方程,得到修正后二项式产能方程:
Pen″-Pewf″=AQ+BQ2;
Pen″为地层压力的修正拟压力,Pewf″为井底流压的修正拟压力。
可以得知,产能方程中的A和B修正前后保持不变,进而保证了气藏开发阶段计算得到的A、B系数的继续有效。
可见,通过实验结果,修正产能方程左端的方便计算的拟压力,保持等式右端不易计算的A、B系数不变,即可实现产能方程的修正,获得任意注采周期气水交互驱替后井产能的表达式。
基于同一发明构思,本发明实施例中还提供了一种储气库气水交互区井产能的预测装置,如下面的实施例所述。由于该装置解决问题的原理与储气库气水交互区井产能的预测方法相似,因此该装置的实施可以参见储气库气水交互区井产能的预测方法的实施,重复之处不再赘述。
图3为本发明实施例中储气库气水交互区井产能的预测装置的示意图。如图3所示,本发明实施例中储气库气水交互区井产能的预测装置可以包括:
研究对象井确定模块301,用于确定气藏净水侵量与压力的函数关系,根据所述函数关系确定在改建储气库运行过程中处于气水交互区的研究对象井;
影响因素确定模块302,用于处理所述研究对象井的气藏开发数据或试井资料,得到气井二项式产能方程,确定气水交互驱替降低产能的影响因素;
渗流实验模块303,用于根据所述气水交互驱替降低产能的影响因素,设计储气库建库过程的渗流实验,测试储气库设计寿命每一个循环周期的气相渗透率变化情况;
产能预测模块304,用于用渗流实验结果修正气井二项式产能方程,得到扣除储气库运行过程中由于气水交互驱替而损失渗流能力的所述研究对象井的产能。
实施例中,研究对象井确定模块301具体可以用于:
通过气藏工程方法或数值模拟方法得到净水侵量W,与储层平均压力P回归成函数关系:
W=f(P);
代入储气库上限压力Pmax,得到储气库运行至上限压力时的水侵量:
Wmax=f(Pmax);
代入储气库下限压力Pmin,得到储气库运行至下限压力时的水侵量:
Wmin=f(Pmin);
根据研究对象井所处的埋藏深度,计算储气库气水界面运移至该埋藏深度时的压力P井,进而得到储气库运行至压力P井时的水侵量:
W井=f(P井);
若Wmax<W井<Wmin,则确定该井为所述研究对象井。
实施例中,影响因素确定模块302具体可以用于:
选择气藏开发阶段两个生产稳定数据点:
生产数据点1(Pe1,Pwf1,Q1);
生产数据点2(Pe2,Pwf2,Q2);
其中,Pe1是生产数据点1的地层压力,Pwf1是生产数据点1的井底流压,Q1是生产数据点1的气体流速;Pe2是生产数据点2的地层压力,Pwf2是生产数据点2的井底流压,Q2是生产数据点2的气体流速;
依据所述两个生产稳定数据点,处理得到拟压力气井产能二项式方程:
Pe'-Pwf′=AQ+BQ2;
其中,P'为真实气体的拟压力:
Pe'为地层压力的拟压力,Pwf′为井底流压的拟压力,Q为气体流速,μ为气体粘度,Z为气体压缩因子,A和B是常数,通过所述两个生产稳定数据点联立求解,是渗透率K的函数;
确定气水交互驱替降低产能的影响因素为储层渗透率K。
实施例中,渗流实验模块303具体可以用于:
根据所述气水交互驱替降低产能的影响因素,设计气驱水-水驱气实验,描述气水交互区气水互驱替的物理过程,得到气水交互区的端点的气相相对渗透率:
K1/K:第一次气水交互驱替后气相相对渗透率,描述储气库第一周期储层经过一次气水互区造成渗流能力降低后剩余的相对渗透率;
K2/K:第二次气水交互驱替后气相相对渗透率,描述储气库第二周期储层经二次气水互区造成渗流能力降低后剩余的相对渗透率;
K3/K:第三次气水交互驱替后气相相对渗透率,描述储气库第三周期储层经过三次气水互区造成渗流能力降低后剩余的相对渗透率;
直至Kn/K:第n次气水交互驱替后气相相对渗透率,描述储气库第n周期储层经过n次气水互区造成渗流能力降低后剩余的相对渗透率。
实施例中,产能预测模块304具体可以用于:
将渗流实验结果代入真实气体的拟压力函数,得到修正的拟压力函数表达式:
将修正后的拟压力表达式代入气井二项式产能方程,得到修正后二项式产能方程:
Pen″-Pewf″=AQ+BQ2;
其中,Pen″为地层压力的修正拟压力,Pewf″为井底流压的修正拟压力。
综上所述,本发明实施例在预测气水交互区井产能时,没有像现有方法直接取气藏开发阶段处理得到的恒定数值,而是扣除由于气水交互驱替而损失的产能,进而给出了储气库寿命内每一个循环周期的产能变化,与现有方法相比,本发明实施例的计算结果准确度大幅度提高,而且随着储气库运行周期的延长,技术优势进越发凸显。本发明实施例的预测方案,将为储气库前期评价提供相对准确的井产能技术参数。
本领域内的技术人员应明白,本发明的实施例可提供为方法、***、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种储气库气水交互区井产能的预测方法,其特征在于,包括:
确定气藏净水侵量与压力的函数关系,根据所述函数关系确定在改建储气库运行过程中处于气水交互区的研究对象井;
处理所述研究对象井的气藏开发数据或试井资料,得到气井二项式产能方程,确定气水交互驱替降低产能的影响因素;
根据所述气水交互驱替降低产能的影响因素,设计储气库建库过程的渗流实验,测试储气库设计寿命每一个循环周期的气相渗透率变化情况;
用渗流实验结果修正气井二项式产能方程,得到扣除储气库运行过程中由于气水交互驱替而损失渗流能力的所述研究对象井的产能。
2.如权利要求1所述的方法,其特征在于,确定气藏净水侵量与压力的函数关系,根据所述函数关系确定在改建储气库运行过程中处于气水交互区的研究对象井,包括:
通过气藏工程方法或数值模拟方法得到净水侵量W,与储层平均压力P回归成函数关系:
W=f(P);
代入储气库上限压力Pmax,得到储气库运行至上限压力时的水侵量:
Wmax=f(Pmax);
代入储气库下限压力Pmin,得到储气库运行至下限压力时的水侵量:
Wmin=f(Pmin);
根据研究对象井所处的埋藏深度,计算储气库气水界面运移至该埋藏深度时的压力P井,进而得到储气库运行至压力P井时的水侵量:
W井=f(P井);
若Wmax<W井<Wmin,则确定该井为所述研究对象井。
3.如权利要求2所述的方法,其特征在于,处理所述研究对象井的气藏开发数据或试井资料,得到气井二项式产能方程,确定气水交互驱替降低产能的影响因素,包括:
选择气藏开发阶段两个生产稳定数据点:
生产数据点1(Pe1,Pwf1,Q1);
生产数据点2(Pe2,Pwf2,Q2);
其中,Pe1是生产数据点1的地层压力,Pwf1是生产数据点1的井底流压,Q1是生产数据点1的气体流速;Pe2是生产数据点2的地层压力,Pwf2是生产数据点2的井底流压,Q2是生产数据点2的气体流速;
依据所述两个生产稳定数据点,处理得到拟压力气井产能二项式方程:
P′e-P′wf=AQ+BQ2;
其中,P'为真实气体的拟压力:
P′e为地层压力的拟压力,P′wf为井底流压的拟压力,Q为气体流速,μ为气体粘度,Z为气体压缩因子,A和B是常数,通过所述两个生产稳定数据点联立求解,是渗透率K的函数;
确定气水交互驱替降低产能的影响因素为储层渗透率K。
4.如权利要求3所述的方法,其特征在于,根据所述气水交互驱替降低产能的影响因素,设计储气库建库过程的渗流实验,测试储气库设计寿命每一个循环周期的气相渗透率变化情况,包括:
根据所述气水交互驱替降低产能的影响因素,设计气驱水-水驱气实验,描述气水交互区气水互驱替的物理过程,得到气水交互区的端点的气相相对渗透率:
K1/K:第一次气水交互驱替后气相相对渗透率,描述储气库第一周期储层经过一次气水互区造成渗流能力降低后剩余的相对渗透率;
K2/K:第二次气水交互驱替后气相相对渗透率,描述储气库第二周期储层经二次气水互区造成渗流能力降低后剩余的相对渗透率;
K3/K:第三次气水交互驱替后气相相对渗透率,描述储气库第三周期储层经过三次气水互区造成渗流能力降低后剩余的相对渗透率;
直至Kn/K:第n次气水交互驱替后气相相对渗透率,描述储气库第n周期储层经过n次气水互区造成渗流能力降低后剩余的相对渗透率。
5.如权利要求4所述的方法,其特征在于,用渗流实验结果修正气井二项式产能方程,得到扣除储气库运行过程中由于气水交互驱替而损失渗流能力的所述研究对象井的产能,包括:
将渗流实验结果代入真实气体的拟压力函数,得到修正的拟压力表达式:
将修正后的拟压力表达式代入气井二项式产能方程,得到修正后二项式产能方程:
P″en-P″ewf=AQ+BQ2;
其中,P″en为地层压力的修正拟压力,P″ewf为井底流压的修正拟压力。
6.一种储气库气水交互区井产能的预测装置,其特征在于,包括:
研究对象井确定模块,用于确定气藏净水侵量与压力的函数关系,根据所述函数关系确定在改建储气库运行过程中处于气水交互区的研究对象井;
影响因素确定模块,用于处理所述研究对象井的气藏开发数据或试井资料,得到气井二项式产能方程,确定气水交互驱替降低产能的影响因素;
渗流实验模块,用于根据所述气水交互驱替降低产能的影响因素,设计储气库建库过程的渗流实验,测试储气库设计寿命每一个循环周期的气相渗透率变化情况;
产能预测模块,用于用渗流实验结果修正气井二项式产能方程,得到扣除储气库运行过程中由于气水交互驱替而损失渗流能力的所述研究对象井的产能。
7.如权利要求6所述的装置,其特征在于,研究对象井确定模块具体用于:
通过气藏工程方法或数值模拟方法得到净水侵量W,与储层平均压力P回归成函数关系:
W=f(P);
代入储气库上限压力Pmax,得到储气库运行至上限压力时的水侵量:
Wmax=f(Pmax);
代入储气库下限压力Pmin,得到储气库运行至下限压力时的水侵量:
Wmin=f(Pmin);
根据研究对象井所处的埋藏深度,计算储气库气水界面运移至该埋藏深度时的压力P井,进而得到储气库运行至压力P井时的水侵量:
W井=f(P井);
若Wmax<W井<Wmin,则确定该井为所述研究对象井。
8.如权利要求7所述的装置,其特征在于,影响因素确定模块具体用于:
选择气藏开发阶段两个生产稳定数据点:
生产数据点1(Pe1,Pwf1,Q1);
生产数据点2(Pe2,Pwf2,Q2);
其中,Pe1是生产数据点1的地层压力,Pwf1是生产数据点1的井底流压,Q1是生产数据点1的气体流速;Pe2是生产数据点2的地层压力,Pwf2是生产数据点2的井底流压,Q2是生产数据点2的气体流速;
依据所述两个生产稳定数据点,处理得到拟压力气井产能二项式方程:
P′e-P′wf=AQ+BQ2;
其中,P'为真实气体的拟压力:
P′e为地层压力的拟压力,P′wf为井底流压的拟压力,Q为气体流速,μ为气体粘度,Z为气体压缩因子,A和B是常数,通过所述两个生产稳定数据点联立求解,是渗透率K的函数;
确定气水交互驱替降低产能的影响因素为储层渗透率K。
9.如权利要求8所述的装置,其特征在于,渗流实验模块具体用于:
根据所述气水交互驱替降低产能的影响因素,设计气驱水-水驱气实验,描述气水交互区气水互驱替的物理过程,得到气水交互区的端点的气相相对渗透率:
K1/K:第一次气水交互驱替后气相相对渗透率,描述储气库第一周期储层经过一次气水互区造成渗流能力降低后剩余的相对渗透率;
K2/K:第二次气水交互驱替后气相相对渗透率,描述储气库第二周期储层经二次气水互区造成渗流能力降低后剩余的相对渗透率;
K3/K:第三次气水交互驱替后气相相对渗透率,描述储气库第三周期储层经过三次气水互区造成渗流能力降低后剩余的相对渗透率;
直至Kn/K:第n次气水交互驱替后气相相对渗透率,描述储气库第n周期储层经过n次气水互区造成渗流能力降低后剩余的相对渗透率。
10.如权利要求9所述的装置,其特征在于,产能预测模块具体用于:
将渗流实验结果代入真实气体的拟压力函数,得到修正的拟压力函数表达式:
将修正后的拟压力表达式代入气井二项式产能方程,得到修正后二项式产能方程:P″en-P″ewf=AQ+BQ2;
P″en为地层压力的修正拟压力,P″ewf为井底流压的修正拟压力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410809184.6A CN104504472B (zh) | 2014-12-23 | 2014-12-23 | 储气库气水交互区井产能的预测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410809184.6A CN104504472B (zh) | 2014-12-23 | 2014-12-23 | 储气库气水交互区井产能的预测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104504472A true CN104504472A (zh) | 2015-04-08 |
CN104504472B CN104504472B (zh) | 2017-09-26 |
Family
ID=52945867
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410809184.6A Active CN104504472B (zh) | 2014-12-23 | 2014-12-23 | 储气库气水交互区井产能的预测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104504472B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104964904A (zh) * | 2015-06-16 | 2015-10-07 | 中国石油天然气股份有限公司 | 储气库水淹气区储层渗透率的测试方法及装置 |
CN107435528A (zh) * | 2016-05-20 | 2017-12-05 | 中国石油天然气股份有限公司 | 火山岩气藏气井配产的方法 |
CN108629471A (zh) * | 2017-03-20 | 2018-10-09 | 中国石油天然气股份有限公司 | 储气库气井交互注采的方法 |
CN113250749A (zh) * | 2021-06-22 | 2021-08-13 | 中国石油化工股份有限公司西北油田分公司 | 一种凝析气藏型储气库的模拟方法及其*** |
CN114756826A (zh) * | 2022-06-16 | 2022-07-15 | 中海油田服务股份有限公司 | 气井产能预测方法、装置、电子设备和存储介质 |
CN115828625A (zh) * | 2022-12-21 | 2023-03-21 | 西南石油大学 | 考虑储层结盐影响的储气库气井产能计算方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101187961A (zh) * | 2007-11-20 | 2008-05-28 | 中国海洋石油总公司 | 射孔完井产能预测*** |
CN101806224A (zh) * | 2010-03-09 | 2010-08-18 | 煤炭科学研究总院重庆研究院 | 井下煤层气抽采产能预测方法 |
CN101906966A (zh) * | 2010-07-16 | 2010-12-08 | 中国石油天然气股份有限公司 | 储层产能的预测方法及装置 |
FR2984562A1 (fr) * | 2011-12-15 | 2013-06-21 | Terra 3E | Procede et systeme de modelisation dynamique d'un ecoulement de fluide polyphasique |
CN103352693A (zh) * | 2013-06-20 | 2013-10-16 | 中国石油天然气集团公司 | 一种基于含油气含量的产能预测***及其方法 |
CN103577886A (zh) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | 一种低渗气藏水平井分段压裂产能预测方法 |
CN104134101A (zh) * | 2014-07-23 | 2014-11-05 | 中国石油集团川庆钻探工程有限公司 | 低渗透储层天然气产能预测方法 |
-
2014
- 2014-12-23 CN CN201410809184.6A patent/CN104504472B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101187961A (zh) * | 2007-11-20 | 2008-05-28 | 中国海洋石油总公司 | 射孔完井产能预测*** |
CN101806224A (zh) * | 2010-03-09 | 2010-08-18 | 煤炭科学研究总院重庆研究院 | 井下煤层气抽采产能预测方法 |
CN101906966A (zh) * | 2010-07-16 | 2010-12-08 | 中国石油天然气股份有限公司 | 储层产能的预测方法及装置 |
FR2984562A1 (fr) * | 2011-12-15 | 2013-06-21 | Terra 3E | Procede et systeme de modelisation dynamique d'un ecoulement de fluide polyphasique |
CN103577886A (zh) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | 一种低渗气藏水平井分段压裂产能预测方法 |
CN103352693A (zh) * | 2013-06-20 | 2013-10-16 | 中国石油天然气集团公司 | 一种基于含油气含量的产能预测***及其方法 |
CN104134101A (zh) * | 2014-07-23 | 2014-11-05 | 中国石油集团川庆钻探工程有限公司 | 低渗透储层天然气产能预测方法 |
Non-Patent Citations (2)
Title |
---|
范新文等: ""涩北气田气井合理配产的综合技术对策"", 《天然气工业》 * |
郭平等: ""储气库多次气水互驱实验"", 《大庆石油地质与开发》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104964904A (zh) * | 2015-06-16 | 2015-10-07 | 中国石油天然气股份有限公司 | 储气库水淹气区储层渗透率的测试方法及装置 |
CN107435528A (zh) * | 2016-05-20 | 2017-12-05 | 中国石油天然气股份有限公司 | 火山岩气藏气井配产的方法 |
CN108629471A (zh) * | 2017-03-20 | 2018-10-09 | 中国石油天然气股份有限公司 | 储气库气井交互注采的方法 |
CN113250749A (zh) * | 2021-06-22 | 2021-08-13 | 中国石油化工股份有限公司西北油田分公司 | 一种凝析气藏型储气库的模拟方法及其*** |
CN113250749B (zh) * | 2021-06-22 | 2024-02-27 | 中国石油化工股份有限公司西北油田分公司 | 一种凝析气藏型储气库的模拟方法及其*** |
CN114756826A (zh) * | 2022-06-16 | 2022-07-15 | 中海油田服务股份有限公司 | 气井产能预测方法、装置、电子设备和存储介质 |
CN115828625A (zh) * | 2022-12-21 | 2023-03-21 | 西南石油大学 | 考虑储层结盐影响的储气库气井产能计算方法 |
CN115828625B (zh) * | 2022-12-21 | 2023-11-21 | 西南石油大学 | 考虑储层结盐影响的储气库气井产能计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104504472B (zh) | 2017-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104504472A (zh) | 储气库气水交互区井产能的预测方法及装置 | |
Simons et al. | A model for overland flow and associated processes within the Hydroinformatics Modelling System | |
Sepehri et al. | Prediction of mining-induced surface subsidence and ground movements at a Canadian diamond mine using an elastoplastic finite element model | |
CN105401927B (zh) | 油藏井网井距确定方法及装置 | |
CN104636818A (zh) | 储气库气油交互区井产能的预测方法及装置 | |
Wu et al. | Flow characteristics and escape-route optimization after water inrush in a backward-excavated karst tunnel | |
Maihemuti et al. | Numerical simulation analysis of reservoir bank fractured rock-slope deformation and failure processes | |
CN104392109A (zh) | 一种油藏的动态离散裂缝模拟方法和装置 | |
Tal et al. | An improved numerical manifold method for simulation of sequential excavation in fractured rocks | |
US10650175B2 (en) | Hydrocarbon migration and accumulation methods and systems | |
CN103469780A (zh) | 一种重力坝坝基深层抗滑稳定临界滑动面的计算方法 | |
CN104346511A (zh) | 一种全新的油藏动态监测方法及装置 | |
Huang et al. | Simulating groundwater inflow in the underground tunnel with a coupled fracture-matrix model | |
Wang et al. | Simulated crack and slip plane propagation in soil slopes with embedded discontinuities using XFEM | |
Beer et al. | Isogeometric boundary element method for the simulation of underground excavations | |
CN104680244A (zh) | 一种油气田产能预测方法及装置 | |
Deng et al. | 3D finite element modeling of directional hydraulic fracturing based on deformation reinforcement theory | |
Fomychov et al. | Determination of parameters of non-uniform fractured rock massif in computing experiment | |
CN103942430A (zh) | 一种基于组合模型的建筑物沉降预测方法 | |
Han et al. | Analytic solutions of the forces and displacements for multicentre circular arc tunnels | |
Jiang et al. | Investigation of coordinated development of coal mining and pipeline protection under boundary effect of thick unconsolidated layer | |
Liu et al. | Analysis of external water pressure for a tunnel in fractured rocks | |
Zhang et al. | Distribution for hydraulic head on tunnel structures in water-rich mountainous region considering influences of fault geology using virtual image technique | |
Zhang et al. | Research on the limit operating pressure of a bedded salt gas storage cavern: A case study | |
Rastorguev et al. | Software DHI feflow and processing modflow in solving groundwater flow problems for HPP projects design and reconstruction in a complex geological and hydrogeological environment |
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 |