天然气PVT参数计算

A1 2 .9 7 0 5 4 7 1 4
A 2 0 .2 8 6 2 6 4 0 5 4
A3 0 .0 0 8 0 5 4 2 0 5 2 2
A 7 0 .0 1 0 4 4 3 2 4 1 3
A1 1 0 .0 0 4 4 1 0 1 5 5 1 2
A1 5 0 .0 0 0 6 0 9 5 7 9 2 6 3
式中： Bg为天然气体积系数 VR为天然气地下体积量，m3 Vsc为天然气地面标准条件下体积量，m3 将天然气状态方程、地面条件代入上式，可得天然气体积 系数计算公式：
B g 3 .4 4 7 1 0
4
ZT p
5.实例计算 根据实验室数据，天然气的临界温度为193K，临界压力为4.6MPa，气层温 度为110 C ，气层压力为30MPa，气体相对密度为0.6，求天然气的pVT数 据。
A 2 1 .0 4 6 7
A7 T
r
3
r
2
(1 A8 ) e x p ( A8 ) 0 .2 7 r
2 2
r r
g (r ) f (r )
'
求得 r 后计算Z值：
A3 0 .5 7 8 3
A 4 0 .5 3 5 3
A5 0 .6 1 2 3
式中:
g
N
2
A 4 2 .8 0 8 6 0 9 4 9
A5 3 .4 9 8 0 3 3 0 5
A 6 0 .3 6 0 3 7 3 0 2
mPa s 为天然气粘度， l
mPa s C 为氮气粘度校正值， O 2 mPa s N 为硫化氢粘度校正值，
2
为在气层温度和大气压下的气体粘度， mPa s
天然气PVT参数计算
• • • • • 1.天然气偏差系数 2.天然气等温压缩系数 3.天然气的粘度 4.天然气的体积系数 5.实例计算
1.天然气偏差系数
天然气是一种真实气体，其状态方程式为
在已知 p r 和 T r 的情况下，需要经过一个迭代过程求解 r ：
r ( i 1) r ( i )
r
2
) Tr r
2
r
A5 A6 r5 A7 r2 A A A Z 1 ( A1 2 3 ) r ( A4 5 ) r2 3 (1 A8 r2 ) exp( A8 r2 ) Tr Tr3 Tr Tr Tr
式中：
A1 0 .3 1 5 1
lg ( g ) 6 .2 4 1 0
3
3
)
r
r
r
r
C O C O 2 (9 .0 8 1 0
2
3
)
式中:
H
2S
H 2 (8 .4 9 1 0
3
lg ( g ) 3 .7 3 1 0
3
)
A 0 2 .4 6 2 1 1 8 2
H S
2
，
，
分别为其含量，%
A1 2 0 .0 8 3 9 3 8 7 1 7 8
g
为天然气相对密度（空气=1.0）
C T 为地层温度，
公式中的 l 可以通过如下的迭代关系式求取
4.天然气的体积系数 天然气体积系数是天然气地下体积量与天然气地面标准条件下体积量 的比值，表示为：
Bg VR V sc
3
r
5
6
3
3
r
T ( A 8 A9 p r A1 0 p A1 1 p ) T ( A1 2 A1 3 p r A1 4 p A1 5 p )
2 2 3
r
3
2
N N 2 (8 .4 8 1 0
2
3
lg ( g ) 9 .5 9 1 0
6
pV ZnRT
式中的Z叫做气体偏差系数，表示在某一温度和压力条件下，同一质量气体的 真实体积与理想体积之比值。利用Dranchuk，Purvis和Robinson（1974）的 相关经验公式计算：
f (r )
3
g (r )
A3 T
r r
f ( r ) A 5 A 6 ( A 4 T r A 5 ) ( A1T r A 2
求得天然气pVT参数为：
气体偏差系数
Z 1 .0 0 3 1 5 1
mPa s
气体粘度
g 0 .0 2 0 1 6 7 3
气体体积系数
B g 0 .0 0 4 2 4 9 9 9
气体等温压缩系数
C g 0 .0 2 5 8 8
M Pa
1
A 6 0 .1 0 4 9
pr p
A 7 0 .6 8 1 6
A8 0 .6 8 4 5
T Tc
Z
0 .2 7 p r
rTr
r
0 .2 7 p r Z Tr
pc
Tr
这种算法收敛很快，只需迭代5次或更少次数，即可达到令人满意的精度。
2.天然气等温压缩系数
天然气等温压缩系数，是指在恒温条件下随压力变化的单位体积变化量，可写为:
mPa s 为二氧化碳粘度校正值，
C O 2 H 2S
A8 0 .7 9 3 3 8 5 6 8 4
A9 1 .3 9 6 4 3 3 0 6
A1 3 0 .1 8 6 4 0 8 8 4 6
A1 0 0 .1 4 9 1 4 4 9 2 5
A1 4 0 .0 2 0 3 3 6 7 8 8 1
Cg 1 V ( V p )T
根据对比压力的定义，则天然气的等温压缩系数可以简化为：
Cg Cr pc
将上式代入天然气状态方程可以得到：
Cg 1 p 1 Z Z p
式中:
A1 0 .3 1 5 1
A 2 1 .0 4 6 7
为了简化方程，这里引入对比压缩系数，其定义如下 C r
1 Cr
A3 0 .5 7 8 3
A 4 0 .5 3 5 3
r Z
Z r pr
A5 0 .6 1 2 3
A 6 0 .1 0 4 9
A7 0 .6 8 1 6
A8 0 .6 8 4 5
式中:
dZ d r 2 1 Tr r
2
r
[5 A5 A 6 2 ( A 4 T r A5 ) ( A1T r A 2
2 2
2S
g Tr l
) A 0 A1 p r A 2 p A3 p T r ( A4 A5 p r A6 p A7 p )
3
r
3
r
2
3
r
l (1.709 10 2.062 10 g )(1.8T 32) 8.188 10 6.15 10 lg( g )
5 2
r r
A3 T
r
2
)r
A7 T
r
2
(1 A8 A ) ex p ( A8 )]
2 2 4 2
r 8 r r
3.天然气的粘度
在地层条件下，天然气粘度是压力、温度和气体组分的函数，根据Carr的关系式有 :
ln (
l l N CO H