【文章內容簡介】 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, )‘校正后的偏差系數(shù) Z’ End If End Sub 模型 +WA 校正法 偏差系數(shù)計算模型 該模型與 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 為給定系數(shù)。 該法應用范圍是 :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 模型計算偏差系數(shù) 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 校正法 偏差系數(shù)計算模型 Dranchuk、 Purvis 和 Robinsion 根據(jù) BenedictWebbRubin 狀態(tài)方程 ,將偏差系數(shù)轉換為對比壓力和對比溫度的函數(shù) ,于 1974 年推導出了帶 8 個常數(shù)的經驗公式 ,其形式為 : 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 為給定系數(shù) 。 ppr 為擬對比壓力 ,無因次 。 Tpr 為擬對比溫度 ,無因次。 DPR 法使用 NewtonRaphson 迭代法解非線性問題可得到偏差系數(shù)的值。這種方法的使用范圍是 :1. 05 ≤ Tpr ≤ 3 。0. 2 ≤ ppr ≤ 30。 部分代碼如下：