天然气工程课程设计报告
姓名: 王岩
学号: 2008*******
班级: 054**-**
专业: 勘查技术与工程(油气井)
学院: 工程学院
指导老师: 宁伏龙老师
目录
程序运行图示
程序代码分析
该井节点分析
未来天然气工业的发展趋势和我国面临的一些瓶颈问题
一、程序运行图示
(HY模型/DPK模型/DPR模型与WA校正法)
(露点模型/泡点模型/等温闪蒸模型)
(Mukherjee-Brill模型)
二、程序代码分析
(HY模型,DPK模型,DPR模型)
+WA校正法
该法以 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 - 1/ Tpr )^2]
该法应用范围是:1. 2 ≤Tpr ≤ 3 ;0. 1 ≤ppr ≤24. 0。
部分代码如下:
Private mand1_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) =
……
Pc(11) =
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 >
F = (-) * (ppr / Tpr) * Exp((-) * ((1 - (1 / Tpr)) ^ 2)) + (x + x ^ 2 + x ^ 3 - x ^ 4) / ((1 - x) ^ 3) - (( / Tpr) - ( / (Tpr ^ 2)) + ( / (Tpr ^ 3))) * (x ^ 2) + (( / Tpr) - ( / (Tpr
乳液及膏霜类护发品的生产工艺 来自淘豆网www.taodocs.com转载请标明出处.