天然气工程课程报告

天然气工程课程报告内容摘要：

Double If T = Then Dim T1 As Double T1 = T + * (p / 2760 ( * (10 ^ (8))) * (p ^ 2)) Tpr1 = T1 / Tpc1 Ppr1 = p / Ppc1 z1 = (1 + ρr + ρr ^ 2 ρr ^ 3) / ((1 ρr) ^ 3) ( / Tpr1 / (Tpr1 ^ 2) + / (Tpr1 ^ 3)) * ρr + ( / Tpr1 / (Tpr1 ^ 2) + / (Tpr1 ^ 3)) * (ρr ^ ( + ( / Tpr1))) HY 模型 WAZ = Format(z1, ) Else Tpr1 = T / Tpc1 Ppr1 = p / Ppc1 z1 = (1 + ρr + ρr ^ 2 ρr ^ 3) / ((1 ρr) ^ 3) ( / Tpr1 / (Tpr1 ^ 2) + / (Tpr1 ^ 3)) * ρr + ( / Tpr1 / (Tpr1 ^ 2) + / (Tpr1 ^ 3)) * (ρr ^ ( + ( / Tpr1))) HY 模型 WAZ = Format(z1, )‘校正后的偏差系数 Z’ End If End Sub 模型 +WA 校正法 偏差系数计算模型 该模型与 DranchukPurvisRobinsion 计算法相同 ,但相对密度应采用牛顿迭代法从下式求出 : 0=1 + (A1 + A2/Tpr+ A3/Tpr^3+ A4/Tpr^4+A5/Tpr^5)* ρr +(A6+A7/Tpr+ A8/Tpr^2)*ρr^2 A9*(A7/Tpr+ A8/Tpr^3) ρr^5 +A10/Tpr^3*ρr ^2*(1 + A11*ρr2 )* exp ( A11*ρr^2 ) 0. 27*ppr/(ρr *Tpr) 式中 : Ai 为给定系数。 该法应用范围是 :1. 0 ≤ Tpr ≤ 3 ,0. 2 ≤ ppr ≤ 30。 或者 0. 7 ≤ Tpr ≤ 1. 0 , ppr 1. 0。 部分代码如下： 39。 计算对比密度（牛顿迭代法） 39。 Dim A1 As Double, A2 As Double, A3 As Double, A4 As Double, A5 As Double, A6 As Double, A7 As Double, A8 As Double, A9 As Double, A10 As Double, A11 As Double A1 = A2 = A3 = A4 = A5 = A6 = A7 = A8 = A9 = A10 = A11 = Dim ρr As Double Dim f11 As Double Dim f12 As Double Dim x As Double, x1 As Double, F As Double, F1 As Double absolution = 1 x = 0 Do While absolution F = () * ppr / Tpr + x + (A1 + A2 / Tpr + A3 / (Tpr ^ 3) + A4 / (Tpr ^ 4) + A5 / (Tpr ^ 5)) * (x ^ 2) + (A6 + A7 / Tpr + A8 / (Tpr ^ 2)) * (x ^ 3) A9 * (A7 / Tpr + A8 / (Tpr ^ 2)) * (x ^ 6) + A10 * (1 + A11 * (x ^ 2)) * ((x ^ 3) / (Tpr ^ 3)) * Exp(A11 * (x ^ 2)) f11 = 1 + 2 * (A1 + A2 / Tpr + A3 / (Tpr ^ 3) + A4 / (Tpr ^ 4) + A5 / (Tpr ^ 5)) * x + 3 * (A6 + A7 / Tpr + A8 / (Tpr ^ 2)) * (x ^ 2) 6 * A9 * (A7 / Tpr + A8 / (Tpr ^ 2)) * (x ^ 5) f12 = (A10 / (Tpr ^ 3)) * (3 * (x ^ 2) + A11 * (3 * (x ^ 4) 2 * A11 * (x ^ 6))) * Exp(A11 * (x ^ 2)) F1 = f11 + f12 x1 = x F / F1 absolution = Abs(x1 x) x = x1 Loop ρ r = x Dim z As Double39。 用 DAK 模型计算偏差系数 39。 z = 1 + (A1 + A2 / Tpr + A3 / (Tpr ^ 3) + A4 / (Tpr ^ 4) + A5 / (Tpr ^ 5)) * ρr + (A6 + A7 / Tpr + A8 / (Tpr ^ 2)) * (ρr ^ 2) A9 * (A7 / Tpr + A8 / (Tpr ^ 2)) * (ρr ^ 5) + A10 * (1 + A11 * (x ^ 2)) * ((x ^ 2) / (Tpr ^ 3)) * Exp((A11) * (x ^ 2)) DAK 模型 Z = Format(z, )‘DPK 模型求得未校正的 Z’ 校正法同上，略。 模型 +WA 校正法 偏差系数计算模型 Dranchuk、 Purvis 和 Robinsion 根据 BenedictWebbRubin 状态方程 ,将偏差系数转换为对比压力和对比温度的函数 ,于 1974 年推导出了带 8 个常数的经验公式 ,其形式为 : Z = 1 +( A1 + A2/ Tpr+ A3/TPr^3)* ρr+ (A4 + A5/Tpr)* ρr^2+(A5 * A6/Tpr)* ρr ^5 A7/TPr^3* ρr^2*(1 + A8* ρr^2 )* exp ( A8* ρr^2) 式中 : Ai 为给定系数。 ppr 为拟对比压力 ,无因次。 Tpr 为拟对比温度 ,无因次。 DPR 法使用 NewtonRaphson 迭代法解非线性问题可得到偏差系数的值。 这种方法的使用范围是 :1. 05 ≤ Tpr ≤ 3。 0. 2 ≤ ppr ≤ 30。 部分代码如下：。