天然氣工程課程報(bào)告-文庫吧資料
2024-10-16 14:09本頁面
【正文】 2 * (A1 + A2 / Tpr + A3 / (Tpr ^ 3)) * x + 3 * (A4 + A5 / Tpr) * (x ^ 2) + 6 * A5 * A6 / Tpr * (x ^ 5) + A7 / (Tpr ^ 3) * (3 * (x ^ 2) + A8 * (3 * (x ^ 4) 2 * A8 * (x ^ 6))) * Exp(A8 * (x ^ 2)) x1 = x F / F1 absolution = Abs(x1 x) x = x1 Loop ρ r = x Dim z As Double39。 部分代碼如下： 39。這種方法的使用范圍是 :1. 05 ≤ Tpr ≤ 3 。 Tpr 為擬對(duì)比溫度 ,無因次。 模型 +WA 校正法 偏差系數(shù)計(jì)算模型 Dranchuk、 Purvis 和 Robinsion 根據(jù) BenedictWebbRubin 狀態(tài)方程 ,將偏差系數(shù)轉(zhuǎn)換為對(duì)比壓力和對(duì)比溫度的函數(shù) ,于 1974 年推導(dǎo)出了帶 8 個(gè)常數(shù)的經(jīng)驗(yàn)公式 ,其形式為 : 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ù) 。用 DAK 模型計(jì)算偏差系數(shù) 39。計(jì)算對(duì)比密度（牛頓迭代法） 39?；蛘?0. 7 ≤ Tpr ≤ 1. 0 , ppr 1. 0。 nhc = n(0) + n(2) nh = n(0) ε = 15 * (nhc nhc ^ 2) + * (nh ^ nh ^ 2) For i = 0 To 11 Tc1(i) = Tc(i) ε Pc1(i) = Pc(i) * Tc1(i) / Tc(i) Next i Dim Tpc1 As Double, Ppc1 As Double, z1 As Double Tpc1 = 0 Ppc1 = 0 For i = 0 To 11 Tpc1 = Tpc1 + n(i) * Tc1(i) Ppc1 = Ppc1 + n(i) * Pc1(i) Next i Dim Tpr1 As Double, Ppr1 As 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ù)計(jì)算模型 該模型與 DranchukPurvisRobinsion 計(jì)算法相同 ,但相對(duì)密度應(yīng)采用牛頓迭代法從下式求出 : 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ù)。 部分 代碼如下： Dim nhc As Double, nh As Double, ε As Double, Tc1(11) As Double, Pc1(11) As Double39。 Tci′為 i 組分的校正臨界溫度 , K。 根據(jù) Wichert2Aziz 的觀點(diǎn) ,每個(gè)組分的臨界溫度和臨界壓力都應(yīng)與參數(shù)ε有關(guān) ,臨界參數(shù)的校正關(guān)系式如下所示 : T ci′ = Tci ε p ci′ = pci * T ci′ / Tci 式中 : Tci 為 i 組分的臨界溫度 ,K。參數(shù)ε的關(guān)系式如下 : ε = 15 (M M2) + 4. 167 ( N0. 5 N2) 式中 :M 為氣體混合物中 H2 S 與 CO2 的摩爾分?jǐn)?shù)之和 。用 HY 模型計(jì)算偏差系數(shù) 39。計(jì)算對(duì)比密度 (牛頓迭代法 )39。計(jì)算擬對(duì)比壓力、擬對(duì)比 溫度 39。計(jì)算擬臨界壓力、擬臨界溫度 39。定義臨界溫度、臨界壓力 39。定義壓力、溫度 39。定義各組分 39。0. 1 ≤ ppr ≤ 24. 0。 /cp 288 208000 根據(jù)： 得 q=1653Mscf/d Forchheimer 方程求解 5012 根據(jù)： 得： q=1645 Mscf/d q=1656Mscf/d ： 1656 Mscf/d 問題 Pwf（ psia） q（ Mscf/d） 壓力法 壓力平方法 15 1994 1067 245 1906 1064 475 1817 1057 704 1726 1044 934 1635 1026 1164 1543 1004 1394 1450 976 1624 1355 943 1854 1259 905 2084 1163 862
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1