天然气工程课程设计报告.docVIP

天然气工程课程设计报告 姓 名： 王 岩 学 号： 2008******* 班 级： 054**-** 专 业： 勘查技术与工程（油气井） 学 院： 工程学院 指导老师： 宁伏龙老师 目录 程序运行图示 程序代码分析 该井节点分析 未来天然气工业的发展趋势和我国面临的一些瓶颈问题 一、程序运行图示 1.欢迎界面以及基本参数的输入 2.酸性气体偏差系数计算(HY模型/DPK模型/DPR模型与WA校正法) 3.天然气相平衡计算(露点模型/泡点模型/等温闪蒸模型) 4．天然气两相管流压力计算(Mukherjee－Brill模型) 5.天然气水合物生成预测及抑制剂用量计算 二、程序代码分析 1.气体偏差系数计算（HY模型，DPK模型，DPR模型） 1.1HY模型+WA校正法 1.1.1HY偏差系数计算模型 　 该法以 Starling-Carnahan 状态方程为基础 ,通过对 Standing-Katz 图版进行拟合 ,得到以下关系式: Z = 0. 06125*[ ppr / (ρ*Tpr ) ]*exp [ - 1. 2 *(1 - 1/ Tpr )^2] ρr 为拟对比密度 ,可用牛顿迭代法从如下方程求得: 0 = (ρr +ρr^2 +ρr^3 – ρr^4)/(1 - ρ r )^3 -(14. 76 / Tpr -9. 76 / Tpr^2+ 4. 58/Tpr^3)*ρr^2+(90. 7/Tpr-242. 2/Tpr^2+ 42. 4/Tpr^3)* ρr ^(2. 18+2. 82/ Tpr)-0. 01652( ppr / Tpr )* exp[ - 1.2*(1 - 1/ Tpr )^2] 该法应用范围是:1. 2 ≤Tpr ≤ 3 ;0. 1 ≤ppr ≤24. 0。 部分代码如下： Private Sub Command1_Click() Dim n(11) As Double -----------------------------------------------------------------------------------定义各组分 Dim i As Integer For i = 0 To 11 n(i) = Val(欢迎界面.摩尔分数(i)) / 100 Next i Dim p As Double, T As Double------------------------------------------------------------------定义压力、温度 p = Val(绝对工作压力) T = Val(绝对工作温度) Dim Tc(11) As Double, Pc(11) As Double----------------------------------------定义临界温度、临界压力 Tc(0) = 373.5833 …… Pc(11) = 295.5315 Dim Tpc As Double, Ppc As Double----------------------------------------计算拟临界压力、拟临界温度 Tpc = 0 Ppc = 0 For i = 0 To 11 Tpc = Tpc + n(i) * Tc(i) Ppc = Ppc + n(i) * Pc(i) Next i Dim Tpr As Double, ppr As Double-----------------------------------------计算拟对比压力、拟对比温度 Tpr = T / Tpc ppr = p / Ppc Dim ρr As Double------------------------------------------------------------------计算对比密度(牛顿迭代法) Dim x As Double, x1 As Double, F As Double, F1 As Double absolution = 1 x = 0 Do While absolution 0.00001 F = (-0.06152) * (ppr / Tpr) * Exp((-1.2) * ((1 - (1 / Tpr)) ^ 2)) + (x + x ^ 2 + x ^ 3 - x ^ 4) / ((1 - x) ^ 3) - ((14.76 / Tpr) - (9.76 / (Tpr ^ 2)) + (4.58 / (Tpr ^