1 相关概念
地质储量:地层原始条件下储层的蕴藏量,即静态储量。
动态储量:以现有的工艺技术和井网开采方式不变的条件下,所有井投入生产直至天然气产量和波及范围内的地层压力降为零时,可以从气藏中采出的天然气总量。**理论上,它等于现有井网控制条件下的地质储量。动储量,是储量的可流动部分,是波及范围内的可流动气量。**动态储量大小与生产井数、井网控制程度和波及系数有关。(它等于地质储量*波及系数)
可采储量:现有井网、工艺技术和经济极限条件下,能从气藏中采出的那一部分天然气量。(它等于地质储量*采收率)
因此,地质储量>动态储量>可采储量。
地质储量与动储量相同点在于与采气速度、开采时间和废弃条件无关,两者的前提都是废弃条件为零。
只有动储量和可采储量才能作为气井合理稳定产能和井网密度的重要依据。
2 动储量计算方法
2.1 压降法
物质平衡又称压降法,目前应用较为广泛且相对最精确的动态储量计算方法。目前物质平衡主要应用的气藏类型有:
- 定容封闭性气藏
- 水驱气藏
- 凝析气藏
- 异常高压气藏
特点:
- 准确性高,可靠性强,需要关井测压;
- 若全气藏关井,采出量应达到3%-5%;若分片关井,采出程度应达到10-15%;
- 地层压力恢复程度影响计算结果,低渗气藏压力恢复时间较长,开发早期误差大。
推导:
设气藏原始条件下,即原始地层压力 p i p_i pi和地层温度条件下,气藏内天然气原始地质储量(地面标准条件0.1Mpa和20 ° C °C °C)为G,它占有地下体积为 G B g i GB_{gi} GBgi。
在压力从 p i p_i pi下降到 p p p的过程中,累积采出气体和水的地面体积为 G p G_p Gp和 W p W_p Wp。
根据地下体积平衡原理可知:地层压力下降
Δ
p
\Delta{p}
Δp的过程中,累积产出天然气和水在压力
p
p
p下的体积为
G
p
B
g
+
W
p
B
w
G_pB_g+W_pB_w
GpBg+WpBw,应等于压力下降
Δ
p
\Delta{p}
Δp的过程中而引起的地下天然气的膨胀量
A
A
A、束缚水的膨胀和气藏孔隙体积的减少引起的含气孔隙体积的减少量
B
B
B、天然水侵量
W
e
W_e
We之和。
地下产出量=A+B+C
-
地下天然气膨胀量
天然气在 p i p_i pi下总体积为 G B g i GB_{gi} GBgi,地面体积为 G G G,在压力 p p p下地下体积为 G B g GB_g GBg。因此压力下降 Δ p \Delta{p} Δp所引起的地下天然气的膨胀量为:
A = G B g − G B g i (1) A=GB_g-GB_{gi}\tag1 A=GBg−GBgi(1) -
含气孔隙体积减小量
含气孔隙体积的减小量等于压力从原始地层压力 p i p_i pi降至压力 p p p时束缚水的膨胀量 d V w dV_w dVw和气藏孔隙体积的减小量 d V p dV_p dVp之和。
B = − d V w + d V p (2) B=-dV_w+dV_p\tag2 B=−dVw+dVp(2)
根据水和岩石有效压缩系数的定义,可分别写成如下形式:
C w = − 1 V w d V w d p C_w=-\frac{1}{V_w}\frac{dVw}{dp} Cw=−Vw1dpdVwC p = 1 V p d V p d p C_p=\frac{1}{V_p}\frac{dVp}{dp} Cp=Vp1dpdVp
式中, C w C_w Cw,水的压缩系数, 1 / M p a 1/Mpa 1/Mpa; C p C_p Cp,岩石的压缩系数, 1 / M p a 1/Mpa 1/Mpa; V w V_w Vw,束缚水的体积, m 3 m^3 m3;
V p V_p Vp,孔隙体积, m 3 m^3 m3;
根据原始条件下天然气的地下体积可分别计算出总孔隙体积 V p V_p Vp和束缚水体积 V w V_w Vw,即:
V p = G B g i 1 − S w i V_p=\frac{GB_{gi}}{1-S_{wi}} Vp=1−SwiGBgi
V w = V p ⋅ S w i = G B g i 1 − S w i ⋅ S w i V_w=V_p\cdot{S_{wi}}=\frac{GB_{gi}}{1-S_{wi}}\cdot{S_{wi}} Vw=Vp⋅Swi=1−SwiGBgi⋅Swi
因此:
B
=
(
C
w
V
w
+
C
p
V
p
)
Δ
p
=
G
B
g
i
1
−
S
w
i
(
C
W
S
w
i
+
C
p
)
Δ
p
(3)
B=(C_wV_w+C_pV_p)\Delta{p}=\frac{GB_{gi}}{1-S_{wi}}(C_WS_{wi}+C_p)\Delta{p}\tag3
B=(CwVw+CpVp)Δp=1−SwiGBgi(CWSwi+Cp)Δp(3)
最后将
A
,
B
,
C
A,B,C
A,B,C代入得到气藏的物质平衡通式:
(
G
B
g
−
G
B
g
i
)
+
G
B
g
i
1
−
S
w
i
(
C
W
S
w
i
+
C
p
)
Δ
p
+
W
e
=
G
p
B
g
+
W
p
B
w
(4)
\color{red} {(GB_g-GB_{gi})+\frac{GB_{gi}}{1-S_{wi}}(C_WS_{wi}+C_p)\Delta{p}+W_e=G_pB_g+W_pB_w}\tag4
(GBg−GBgi)+1−SwiGBgi(CWSwi+Cp)Δp+We=GpBg+WpBw(4)
整理得:
G
=
G
p
B
g
−
(
W
e
−
W
p
B
w
)
B
g
i
[
(
B
g
B
g
i
−
1
)
+
(
C
w
S
w
i
+
C
P
1
−
S
w
i
)
Δ
p
]
(5)
\color{red} { G=\frac{G_pB_g-{(W_e-W_pB_w)}}{B_{gi}[(\frac{B_g}{B_{gi}}-1)+(\frac{C_wS_{wi}+C_P}{1-S_{wi}})\Delta{p}]} }\tag 5
G=Bgi[(BgiBg−1)+(1−SwiCwSwi+CP)Δp]GpBg−(We−WpBw)(5)
无水驱气藏的物质平衡:
油藏没有水驱作用,即
W
e
=
0
,
W
p
=
0
W_e=0,W_p=0
We=0,Wp=0,可得无水驱气藏的物质平衡方程式:
G
p
B
g
=
G
(
B
g
−
B
g
i
)
+
G
B
g
i
(
C
w
S
w
i
+
C
p
1
−
S
w
i
)
Δ
p
(6)
\color{red}{G_pB_g=G(B_g-B_{gi})+GB_{gi}(\frac{C_wS_{wi}+C_p}{1-S_{wi}})\Delta{p}}\tag 6
GpBg=G(Bg−Bgi)+GBgi(1−SwiCwSwi+Cp)Δp(6)
上式右端第二项与右端第一项相比数值很小,可以忽略不计时,即认为开采过程中孔隙体积不变,即可转化为定容封闭气藏的物质平衡方程式:
G
p
B
g
=
G
(
B
g
−
B
g
i
)
(6)
\color{red}{G_pB_g=G(B_g-B_{gi})}\tag 6
GpBg=G(Bg−Bgi)(6)
天然气的目前和原始体积系数分别为:
B
g
=
p
s
c
Z
T
p
T
s
c
(7)
B_g=\frac{p_{sc}ZT}{pT_{sc}}\tag7
Bg=pTscpscZT(7)
B g i = p s c Z i T p i T s c (8) B_{gi}=\frac{p_{sc}Z_iT}{p_iT_{sc}}\tag8 Bgi=piTscpscZiT(8)
式中, T T T,地层温度, K K K;
T s c T_{sc} Tsc,地面标准状态下的温度, 293.15 K 293.15K 293.15K;
p p p,目前地层压力, M p a Mpa Mpa;
p s c p_{sc} psc,地面标准状态下压力, 0.101 M p a 0.101Mpa 0.101Mpa;
Z Z Z,压力 p p p下天然气的偏差因子, f f f;
Z i Z_i Zi, p i p_i pi下天然气的偏差因子, f f f。
将公式(7)和公式(8)代入公式(6)式整理得:
p
Z
=
p
i
Z
i
(
1
−
G
p
G
)
(9)
\color{red} \frac{p}{Z}=\frac{p_i}{Z_i}(1-\frac{G_p}{G})\tag9
Zp=Zipi(1−GGp)(9)
上式即为定容封闭气藏的压降方程式。该压降是在忽略束缚水膨胀和孔隙体积减小的情况下导出的,实际工作中,即使对于正常压力系统的气藏,(6)式右端第二项也可能达到右端第一项的十分之一,因此建议最好不要忽略(6)式右端第二项。
另外一种定容气藏推导:
G
=
G
p
+
G
B
g
i
B
g
B
=
V
g
V
g
s
=
Z
n
R
T
/
p
Z
s
c
n
R
T
s
c
/
p
s
c
=
p
s
c
Z
s
c
T
s
c
Z
T
p
B
g
i
=
V
g
i
V
g
s
=
p
s
c
Z
s
c
T
s
c
Z
i
T
i
p
i
B
g
i
B
g
=
p
/
T
Z
p
i
/
T
i
Z
i
=
p
/
Z
p
i
/
Z
i
G=G_p+\frac{GB_{gi}}{B_g} \\ B=\frac{V_g}{V_{gs}}=\frac{ZnRT/p}{Z_{sc}nRT_{sc}/p_{sc}}=\frac{p_{sc}}{Z_{sc}T_{sc}}\frac{ZT}{p} \\ B_{gi}=\frac{V_{gi}}{V_{gs}}=\frac{p_{sc}}{Z_{sc}T_{sc}}\frac{Z_iT_i}{p_i} \\ \frac{B_{gi}}{B_g}=\frac{p/TZ}{p_i/T_iZ_i}=\frac{p/Z}{p_i/Z_i}\\
G=Gp+BgGBgiB=VgsVg=ZscnRTsc/pscZnRT/p=ZscTscpscpZTBgi=VgsVgi=ZscTscpscpiZiTiBgBgi=pi/TiZip/TZ=pi/Zip/Z
因此:
G
=
G
p
+
G
p
/
Z
p
i
/
Z
i
p
Z
=
p
i
Z
i
(
1
−
G
p
G
)
\color{red}{G=G_p+G\frac{p/Z}{p_i/Z_i} \\ \frac{p}{Z}=\frac{p_i}{Z_i}(1-\frac{G_p}{G})}
G=Gp+Gpi/Zip/ZZp=Zipi(1−GGp)
拟压力:
p
p
=
p
Z
p_p=\frac{p}{Z}
pp=Zp
定容气藏的生产曲线:
p
p
=
p
p
i
(
1
−
G
p
G
)
\color{red}p_p=p_{pi}(1-\frac{G_p}{G})
pp=ppi(1−GGp),
p
p
=
a
−
b
G
p
,
a
=
p
p
i
,
b
=
p
p
i
G
p_p=a-bG_p,a=p_{pi}, b=\frac{p_{pi}}{G}
pp=a−bGp,a=ppi,b=Gppi
适应性:
物质平衡是利用气藏视地层压力( p Z \frac{p}{Z} Zp)与累积产气量( G p G_p Gp)所构成的“压降图”来确定气藏的储量的。压降法是物质平衡在封闭气藏应用中的一个特例。它可以用于:利用单井稳定生产后同期关井测压,获得稳定地层压力,计算单井储量;选用气井产量大,关井易稳定,井位处于构造中部的井作为气藏关井地层压力的代表,取等时关井时测得的地层压力与气藏累积产量作图,可计算出气藏储量;或气藏统一关井测压,获得各井稳定关井压力,选用算术平均及加权平均计算历次关井气藏平均地层压力及累积采气量也可计算气藏储量。
压降法不需要任何地质参数,对于地质结构复杂,无法求准储量空间的气藏,如碳酸盐岩裂缝性气藏,采用该法最好。
对于活跃的水压驱动气藏,由于开采中压力不下降(或下降不明显),对于异常高压气藏,由于不能忽略水的压缩性和岩石的变形,“压降图”在此时不再是一条直线。另外,在计算累积产气量时,由于实际放空量的估计往往误差较大,在一定程度上影响了储量的精度。
对于低渗透气藏,适应性也受到一定的限制。低渗透气井生产一定时间后关井,关井初期压力恢复较快,此后压力恢复越来越缓慢,需要相当长时间才能恢复至稳定(大量开发实践证明,一般需要0.5-1年)。为了不过多影响正常生产,通常关井时间有限(长庆气田一般1个月左右),测试的恢复压力因稳定而偏低。低渗透气藏渗流遵循“先易后难,先好后差”的自然规律,开发初期主要沿裂缝和相对高渗透带向生产井流动,此阶段测试得到的产气量和地层压力主要反映裂缝和相对高渗透带的天然气渗流特征,随着开发时间的延长,动态资料才逐渐反映低渗透带、更低渗透带的渗流特征。这导致得到的储量曲线呈“多段型”,每段对应不同物性气层参与渗流的开发阶段。
2.2 弹性第二相法
弹性第二相法是根据压降试井的压力变化而得出的一种方法,是目前计算动态储量的主要方法之一。对于一个有限封闭的气藏,当气井以稳定产量开井生产,有界封闭地层开井生产井底压力降落曲线一般可分为三个阶段:
-
第一阶段,不稳定早期,指压降漏斗没有传到边界之前的弹性第一阶段;
-
第二阶段,不稳定晚期,即压降漏斗传到边界之后;
-
第三阶段,拟稳定期,任一点压降速度相同,此阶段又称为弹性第二相过程,井底压力随时间变化关系为:
P w f 2 = P e 2 − 8.48 × 1 0 − 3 Q u K h p s c Z T T s c [ l g ( R e r w ) − 0.326 + 0.435 S ] − 2 Q p e t G C t P_{wf}^2=P_e^2-\frac{8.48\times10^{-3}Qu}{Kh}\frac{p_{sc}ZT}{T_{sc}}[lg(\frac{R_e}{r_w})-0.326+0.435S]-\frac{2Qp_et}{GC_t} Pwf2=Pe2−Kh8.48×10−3QuTscpscZT[lg(rwRe)−0.326+0.435S]−GCt2Qpet
令:
A = 2 Q p e G C t B = P e 2 − 8.48 × 1 0 − 3 Q u K h p s c Z T T s c [ l g ( R e r w ) − 0.326 + 0.435 S ] A=\frac{2Qp_e}{GC_t}\\ B=P_e^2-\frac{8.48\times10^{-3}Qu}{Kh}\frac{p_{sc}ZT}{T_{sc}}[lg(\frac{R_e}{r_w})-0.326+0.435S] A=GCt2QpeB=Pe2−Kh8.48×10−3QuTscpscZT[lg(rwRe)−0.326+0.435S]
则:
P
w
f
2
=
B
−
A
t
G
=
2
Q
p
e
A
C
t
P_{wf}^2=B-At\\ G=\frac{2Qp_e}{AC_t}
Pwf2=B−AtG=ACt2Qpe
特点:
- 定产生产达到拟稳态;
- 高精度仪表测流压;
- 测试过程中产量波动不超过10%;
- 修正等时试井的延续生产段的资料可以应用该方法。
2.3 压力恢复法
气井关井前需要长时间的稳产,且恢复探测到边界。
不稳定早期的气井压力恢复曲线为:
p
w
s
2
=
p
T
2
+
4.24
×
1
0
−
3
Q
μ
K
h
p
s
c
Z
T
T
s
c
[
l
g
8.085
K
t
ϕ
μ
C
t
r
w
2
+
0.875
S
]
p_{ws}^2=p_T^2+\frac{4.24\times10^{-3}Q\mu}{Kh}\frac{p_{sc}ZT}{T_{sc}}[lg\frac{8.085Kt}{\phi{\mu}C_tr_w^2}+0.875S]
pws2=pT2+Kh4.24×10−3QμTscpscZT[lgϕμCtrw28.085Kt+0.875S]
令:
m
=
4.24
×
1
0
−
3
Q
μ
K
h
p
s
c
Z
T
T
s
c
A
=
p
T
2
+
m
[
l
g
8.085
K
t
ϕ
μ
C
t
r
w
2
+
0.875
S
]
m=\frac{4.24\times10^{-3}Q\mu}{Kh}\frac{p_{sc}ZT}{T_{sc}}\\ A=p_T^2+m[lg\frac{8.085Kt}{\phi{\mu}C_tr_w^2}+0.875S]
m=Kh4.24×10−3QμTscpscZTA=pT2+m[lgϕμCtrw28.085Kt+0.875S]
即:
p
w
s
2
=
m
l
g
t
+
A
G
=
0.1077
Q
t
e
m
C
t
p
e
p_{ws}^2=mlgt+A\\ G=0.1077\frac{Qt_e}{mC_t}p_e
pws2=mlgt+AG=0.1077mCtQtepe
特点:
- 必须先求得地层平均压力和压力恢复达到边界时间以及斜率;
- 需关井前长时稳产,关井恢复需探测到边界。
2.4 压差曲线法
2.5 产量累计法
2.6 产量递减法
2.7 产量不稳定分析法
**动态物质平衡法:利用单井生产历史数据(产量和压力),不关井进行物质平衡分析,进而对气藏生产动态进行分析、评价和预测的新方法。**
特点:
- 利用单井的生产历史数据(产量和压力);
- 不必进行关井测压,不必定产或定压生产;
- 产量和井底压力都可以是变化的。
用途:
- 计算单井控制动储量;
- 判断是否达到泄流边界;
- 计算单井泄流面积和泄流半径;
- 计算目前地层压力;
- 计算渗透率和表皮系数,确定改造井;
- 判断是否存在复杂地质边界;
- 动态预测(配产、稳定期、采收率)
原理:
压降法需要关井至压力稳定,而生产过程就是物质平衡过程,如何在生产动态中运用压降法?
物质平衡方程:
P
‾
Z
‾
=
P
i
Z
i
(
1
−
G
p
G
)
d
d
t
(
P
‾
Z
‾
)
=
−
p
i
z
i
G
⋅
d
G
p
d
t
=
−
p
i
q
z
i
G
\frac{\overline{P}}{\overline{Z}}=\frac{P_i}{Z_i}(1-\frac{G_p}{G})\\ \frac{d}{dt}(\frac{\overline{P}}{\overline{Z}})=-\frac{p_i}{z_iG}\cdot{\frac{dG_p}{dt}}=-\frac{p_iq}{z_iG}
ZP=ZiPi(1−GGp)dtd(ZP)=−ziGpi⋅dtdGp=−ziGpiq
物质平衡方程:
P
‾
Z
‾
=
P
i
Z
i
(
1
−
G
p
G
)
d
d
t
(
P
‾
Z
‾
)
=
−
p
i
z
i
G
⋅
d
G
p
d
t
=
−
p
i
q
z
i
G
\frac{\overline{P}}{\overline{Z}}=\frac{P_i}{Z_i}(1-\frac{G_p}{G})\\ \frac{d}{dt}(\frac{\overline{P}}{\overline{Z}})=-\frac{p_i}{z_iG}\cdot{\frac{dG_p}{dt}}=-\frac{p_iq}{z_iG}
ZP=ZiPi(1−GGp)dtd(ZP)=−ziGpi⋅dtdGp=−ziGpiq
参考
- https://wenku.baidu.com/view/bb1ef9f4a22d7375a417866fb84ae45c3a35c250.html (田冷)
- http://www.doc88.com/p-118697856124.html