下载文档原格式
摘要天然气的压缩因子、体积系数、压缩系数、粘度等高压物性参数随气藏压力和温度的变化而变化,定量描述和预测这些参数的变化规律具有十分重要的实际意义。
通过电脑程序将天然气高压物性的相关经验公式转变为程序计算，能够很快的得到计算结果以及对计算结果的图形分析，通过最后的图形分析我们能很直观的看出高压物性参数之间的关系，有利于确定合理的开发速度和规模，节省投资，将资金投向回报率高的方案上。
本文中首先介绍了天然气高压物性参数的相关经验及半经验公式，再利用excelVBA实现公式的程序计算，只要输入原始数据，点击相应的高压物性计算按钮就能得到计算结果，数据分析窗体能够自动输出高压物性与相关参数的图形。
关键字：高压物性偏差系数粘度压力温度ABSTRACTGas compressibility factor, volume factor, compressibility, viscosity and other physical parameters with high pressure gas reservoir pressure and temperature changes, quantitatively describe and predict the variation of these parameters has a very important practical significance. Through a computer program related to the physical properties of high pressure natural gas into a program to calculate the empirical formula, can quickly get the results and the calculation results of the graphical analysis, graphical analysis through the last we can see the pressure very intuitive relationship between the physical parameters will help determine a reasonable pace and scale of development, reduce investment, high return on the capital investment program on.This paper first introduces the gas pressure in the physical parameters relevant experience and semi.empirical formulas to achieve reuse excelVBA program calculates the formula, as long as the input raw data, click on the appropriate button to be able to calculate the physical properties of high pressure to get results, analysis of data form can be automatically output pressure properties and related parameters graphics.Keywords: PVT variation ；coefficient of viscosity；pressure ；temperature ；coefficient of volume.重庆科技学院本科生毕业设计目录目录摘要 (I)ABSTRACT (II)1 绪论 (1)1.1国内外天然气高压物性参数计算发展历史 (1)1.2国内天然气分布 (2)1.3天然气高压物性参数计算的意义 (2)2 计算方法介绍 (4)2.1天然气临界参数计算 (4)2.1.1 天然气平均分子量 (4)2.1.2 天然气的相对密度 (4)2.1.3 拟临界压力P PC和拟临界温度T PC (4)2.1.4 拟对比压力P Pr和拟对比温度T Pr的计算 (5)2.2天然气的偏差因子Z计算 (6)2.2.1 Pong.Robinson方程法 (6)2.2.2 Cranmer方法 (7)2.2.3 DPR法 (7)2.2.4 DAK法 (8)2.2.5 平均值法 (9)2.3天然气压缩因子计算 (9)2.4天然气体积系数计算 (10)2.5天然气膨胀系数的计算 (10)2.6天然气的粘度计算 (10)2.6.1 Lee关系式 (11)2.6.2 Dempsey关系式 (11)3 EXCELVBA程序计算 (13)3.1VBA简介 (13)3.2界面设计 (13)3.3操作步骤 (14)3.3.1 原始数据 (15)3.3.2 拟临界压力、温度，拟对比压力、温度的计算 (16)3.3.3 天然气高压物性的计算 (17)3.3.4 图形分析 (22)3.3.5 数据查询 (25)4 结论 (26)参考文献 (27)致谢 (28)1重庆科技学院本科生毕业设计 1 绪论1 绪论1.1 国内外天然气高压物性参数计算发展历史天然气高压物性参数计算问题早在20世纪40年代就有人提出了。
原油高压物性（PVT）实验描述在不同压力下，油藏流体的相平衡状态会发生变化。
一个油田在开发早期，最好就抓紧取样，开展原油高压物性（PVT）实验，使样品能尽量接近原始油藏流体。
通过PVT实验，掌握油藏流体及其在不同压力下的体积特征，为我们对油田动态预测奠定一个坚实的基础。
原油高压物性（PVT）实验有两类：
一类是等组分膨胀实验，它是把烃类流体样品在油藏温度及超过油藏原始压力下放入PVT容器中，在等温条件下逐步减少容器的压力，测量烃类体积在每个压力下的变化。
这项实验的目的在于确定:
（1）饱和压力（泡点压力，原油内的溶解气开始分离出去时的压力）；
（2）高于饱和压力时在油藏温度条件下的单相流体的压缩系数；
（3）总烃类体积与压力的函数关系。
另一类是差异分离实验，它是在油藏开发过程中，随着压力降低，从原油中分离出来的溶解气不断地被采出来，在油藏中气与液相也不断重新建立新的平衡，这项实验的目的在于确定:
（1）溶解气与压力的函数关系；
（2）原油体积的收缩率与压力的关系；
（3）分离气体的组成、压缩系数和相对密度；
（4）剩余油的密度、黏度与压力的函数关系。
高压物质的性质与相变高压物质指的是在极端压力下表现出非常特殊性质的物质。
它们在普通大气压下可能表现得和其他物质无异，但当经受到极高压力的作用后，它们的行为和性质就开始呈现出与众不同的特点。
本文将讨论高压物质的性质和相变，揭示这些物质背后的奇妙之处。
一、高压引发的结构和性质变化当物质受到高压的作用时，其原子或分子之间的距离会变得更近，导致原子核之间的排斥力增加。
在这种情况下，物质中的原子或分子将产生更强的相互作用力，从而引发结构和性质的变化。
以晶体为例，普通情况下晶体的原子或离子排列有一定的规律，呈现出比较稳定的结构。
然而，当晶体所受的压力达到一定程度后，其结构会发生重塑。
一些晶体可能会变得更加紧密，形成更密实的晶格结构，而另一些晶体则可能出现亚稳态结构。
高压还可能影响物质的电学、磁学以及导电性等性质。
许多物质在高压下的电导率会发生显著的变化，部分物质甚至会表现出与室温下截然不同的电学性质。
此外，高压也会改变物质的磁性质，使得原本不带磁矩的物质产生磁矩，或改变已存在的磁矩的方向。
二、高压下的相变现象高压条件下物质的相变行为也非常引人注目。
通常情况下，固体在受到足够高的压力时，会转变为液体或气体。
这种相变被称为熔化相变和汽化相变。
然而，在高压下，相变现象会变得更为复杂。
对于某些物质，高压可能导致它们由液体直接转变为固体，跳过熔化的过程，这被称为“超熔相变”。
这种相变行为在很多有机化合物以及金属中都被观察到。
在某些情况下，高压下物质的相变行为变得极其不稳定。
例如，一些晶体在高压下会转变为金属液态，这就是所谓的“金属-非金属转变”。
这种相变行为常常伴随着电学、磁学性质的剧变，同时也引发了材料学和地球科学领域的广泛研究。
三、高压物质的应用前景高压物质的研究不仅仅具有理论意义，还有着广泛的实际应用前景。
首先，高压能够改变物质的化学性质，促使一些反应更快、更有效地进行。
通过对高压下化学反应的研究，我们能够探索新的反应途径和合成方法，为新材料的开发和生产提供更广阔的空间。
中国石油大学渗流物理实验报告
实验日期: 成绩: 班级: 学号: 姓名: 教师:
同组者:
不可压缩流体平面径向稳定渗流实验
一、实验目的
1、平面径向渗流实验是达西定律在径向渗流方式下的体现，通过本实验加深对达西定律的理解；
2、要求熟悉平面径向渗流方式下的压力降落规律，并深刻理解该渗流规律与单向渗流规律的不同，进而对渗透率突变地层、非均质地层等复杂情况下的渗流问题及其规律深入分析和理解。
二、实验原理
平面径向渗流实验以稳定渗流理论为基础，采用圆形填砂模型，以流体在模型中的流动模拟水平均质地层中不可压缩流体平面径向稳定渗流过程。
保持填砂模型内、外边缘压力恒定，改变出口端流量，在稳定条件下测量填砂模型不同位置处的水头高度，可绘制水头高度或压力随位置的变化曲线(压降漏斗曲线)；根据平面径向稳定渗流方程的解计算填砂模型的流动系数及渗透率。
三、实验流程
图2-1 平面径向流实验流程图
1－测压管（模拟井）；2～16－测压管（共16根）；18―圆形边界（填砂模型）；19－排液管（生产井筒）；20—量筒； 21—进水管线；22—供液筒；23－溢流管；24—排水阀；25—进水阀；26—供水阀。
四、实验步骤
1、记录填砂模型半径、填砂模型厚度，模拟井半径、测压管间距等数据。
2、打开供水阀“26”，打开管道泵电源，向供液筒注水，通过溢流管使供液筒内液面保持恒定。
3、关闭排水阀“24”，打开进水阀“25”向填砂模型注水。
4、当液面平稳后，打开排水阀“24”，控制一较小流量。
5、待液面稳定后，测试一段时间内流入量筒的水量，重复三次。
；
6、记录液面稳定时各测压管内水柱高度。
7、调节排水阀，适当放大流量，重复步骤5、6；在不同流量下测量流量及各测压管高度，共测三组流量。
8、关闭排水阀24、进水阀25，结束实验。
注：待学生全部完成实验后，先关闭管道泵电源，再关闭供水阀26。
五、实验要求及数据处理
1、实验要求
(1)将原始数据记录于测试数据表中，根据记录数据将每组的3个流量求平均值，并计算测压管高度、压力；绘制三个流量下压力随位置的变化曲线（压降漏斗曲线），说明曲线形状及其原因。
(2)根据平面径向稳定渗流方程，计算填砂模型平均渗透率、不同半径范围的渗透率，评价砂体的均匀性。
(3)写出填砂模型流量与总压差的关系表达式，并绘出流量与总压差的关系曲线。
2、数据处理
流量与总压差的关系表达式：
()
w
e
w e R R P P Kh ln 2Q μπ-=
（2-1）
任意半径范围的渗透率计算公式：
()212
12ln
P P h r r Q K -=
πμ （2-2）
式中：e P —模型外边缘压力，10-1MPa ；
w P —模型出口端面压力，10-1MPa ；
e R —供给边缘半径，cm ； w R —井筒半径，cm ；
h —地层厚度，cm ； μ—流体粘度，s mPa ⋅；
1P 、2P —任意半径1r 、2r 处的压力，10-1MPa 。
3、平面径向流实验数据与处理结果记录表
实验设备编号： 径3#
表2 测压管液面读数与高度、压力计算结果记录表
其中：H 为实际读数；h 为计算高度，P 为计算压力。
填砂模型（内）半径＝ 18.0 cm ， 填砂厚度＝ 2.5 cm ，
中心孔（内）半径＝ 0.3 cm ， 相邻两测压管中心间距＝ 4.44 cm ，
4、单向流实验数据处理和曲线绘制
(1)举例说明水头高度、压力、和流量的计算过程，并将计算结果填在表1中；绘制三个流量下压力随位置的变化曲线（压降漏斗曲线），说明曲线形状及其原因。
以H=66.5cm 为例进行计算举例：
h=H-2.5/2-0.4=66.5-1.25-0.4=64.85cm p=ρgh=1000*9.8*0.6458=6355.3Pa
s)(cm3/ 78.23
84
.267.283.23321=++=++=
Q Q Q Q
图2 压降漏斗曲线
对于每个流量下的压降曲线，其基本关于纵轴对称，表明距离井筒相同的点压差基本一样；越靠近井筒压差下降越快，距离增大时压差增加缓慢，基本不怎么变化，原因在于井筒的压力降对远处的影响减弱；曲线基本呈对数形态，流量越大，曲线越靠下。
(2)根据平面径向稳定渗流方程，计算填砂模型平均渗透率、不同半径范围的渗透率，评价砂体的均匀性。
①以流量1为例：cm h 65.652
85
.6545.651=+=
ε )1.0(0643.07.64336565.0*8.9*1000
11MPa Pa h g P ====εερ (0.1MPa)
0564.09.56395755.0*8.9*100012====Pa gh P w ρ )(72.91)
0564.00643.0(*5.2*23.018
ln
*1*78.2)(2ln
2212111m P P h r r Q K μππμ=-=-=
同理：cm h 75.642=ε，)1.0(0635.02MPa P =ε，)1.0(0494.02MPa P w =
)(37.83)
0494.00635.0(*5.2*23.018
ln
*1*51.4)(2ln
2212122m P P h r r Q K μππμ=-=-=
cm h 05.603=ε，)1.0(0588
.03MPa P =ε，)1.0(0345.03MPa P w = )(94.72)
0345.00588.0(*5.2*23.018
ln
*1*80.6)(2ln
2212133m P P h r r Q K μππμ=-=-=
所以模型的平均渗透率为：
)(68.823
94
.7237.8372.9132321m K K K K μ=++=++=
以半径范围为0至4.44cm 的地层平均渗透率的计算为例： 在流量1的情况下：
)1.0(0635.04.63502
6485
.06475.0*
8.9*100011MPa Pa h g P ==+==ρ
(0.1MPa)
0564.09.56395755.0*8.9*100022====Pa h g P ρ )(17.67)
0564.00635.0(*5.2*23.044
.4ln
*1*78.2)(2ln
22121111m P P h r r Q K μππμ=-=-=
同理在流量2的情况下：
)1.0(0622
.0'1MPa P =，)1.0(0494.0'2MPa P = )(44.60)0494
.00622.0(*5.2*23.044
.4ln
*1*51.4)(2ln
22'
1'21212m P P h r r Q K μππμ=-=-= )1.0(0570
.0'
'1MPa P =，)1.0(0345.0''2MPa P = )(84.51)0345
.00570.0(*5.2*23.044
.4ln
*1*80.6)(2ln
22'
'1''21313m P P h r r Q K μππμ=-=-=、
所以，此半径范围内的地层平均渗透率为：
)(82.583
84
.5144.6017.67321312111m K K K K μ=++=++=
以此类推，在计算半径范围为4.44cm 至8.88cm 的地层平均渗透率时，1r 取8.88cm ，2r 取4.44cm ；在计算半径范围为8.88cm 至13.32cm 的地层平均渗透率时，1r 取13.32cm ，2r 取8.88cm ；在计算半径范围为13.32cm 至17.76cm 的地层平均渗透率时，1r 取17.76cm ，
2r 取13.32cm ；其中所有对应的1P 、2P 应当取半径相对应的测压管内所测压力。
(3)写出填砂模型流量与总压差的关系表达式，并绘出流量与总压差的关系曲线。
图3 流量与总压差的关系曲线
流量与总压差的关系表达式：
()w
e
w e R R P P Kh ln 2Q μπ-=
由上式可得出流量与总压差之间呈线性关系，由流量与总压差的关系曲线可知实验数据基本满足线性关系。
六．实验总结
由于实验设备的长期使用填砂管的均质性变差，并非严格的均质。
其次，由于测得的流量数据组数不多，即使图像近似线性，由曲线作出的图像并不能充分证明流量与压差的线性关系。
再次，实验过程读数必须等到稳定时进行，且读数必须按照读数规则进行。
