这是一个使用VBA编写的宏，用于根据GB/T17747.3-2011标准计算天然气的压缩因子。代码中定义了两个函数`B`和`C`，分别计算系数，并在`ysyz`子程序中应用这些函数，根据输入的压力、温度和成分参数进行计算。宏通过迭代逐步逼近法求解，适用于Excel模板中进行天然气属性计算。
Function B(T, HCHu, x1, x2, x3, x4, x5) As Double
Dim Bij() As Double
ReDim Preserve Bij(1 To 5, 1 To 5)
On Error Resume Next
Bij(1, 1) = -0.425468 + 0.002865 * T - 0.00000462073 * T * T
Bij(1, 1) = Bij(1, 1) + (0.000877118 - 0.00000556281 * T + 0.0000000088151 * T * T) * HCHu
Bij(1, 1) = Bij(1, 1) + (-0.000000824747 + 0.00000000431436 * T - 6.08319E-12 * T * T) * HCHu * HCHu
Bij(2, 2) = -0.1446 + 0.00074091 * T - 0.00000091195 * T * T
Bij(3, 3) = -0.86834 + 0.0040376 * T - 0.0000051657 * T * T
Bij(4, 4) = -0.00110596 + 0.0000813385 * T - 0.000000098722 * T * T
Bij(5, 5) = -0.13082 + 0.00060254 * T - 0.0000006443 * T * T
Bij(1, 2) = (0.72 + 0.00001875 * (320 - T) * (320 - T)) * (Bij(1, 1) + Bij(2, 2)) / 2
Bij(1, 3) = -0.865 * (Bij(1, 1) * Bij(3, 3)) ^ (1 / 2)
Bij(1, 4) = -0.052128 + 0.00027157 * T - 0.00000025 * T * T
Bij(1, 5) = -0.068729 - 0.00000239381 * T + 0.000000518195 * T * T
Bij(2, 3) = -0.339693 + 0.00161176 * T - 0.00000204429 * T * T
Bij(2, 4) = 0.012
B = x1 * x1 * Bij(1, 1) + 2 * x1 * x2 * Bij(1, 2) + 2 * x1 * x3 * Bij(1, 3)
B = B + 2 * x1 * x4 * Bij(1, 4) + 2 * x1 * x5 * Bij(1, 5) + x2 * x2 * Bij(2, 2)
B = B + 2 * x2 * x3 * Bij(2, 3) + 2 * x2 * x4 * Bij(2, 4)
B = B + x3 * x3 * Bij(3, 3) + x4 * x4 * Bij(4, 4) + x5 * x5 * Bij(5, 5)
End Function
Function C(T, HCH, x1, x2, x3, x4, x5) As Double
Dim Cijk() As Double
ReDim Preserve Cijk(1 To 5, 1 To 5, 1 To 5)
On Error Resume Next
Cijk(1, 1, 1) = -0.302488 + 0.00195861 * T - 0.00000316302 * T * T
Cijk(1, 1, 1) = Cijk(1, 1, 1) + (0.000646422 - 0.00000422876 * T + 0.00000000688157 * T * T) * HCH
Cijk(1, 1, 1) = Cijk(1, 1, 1) + (-0.000000332805 + 0.0000000022316 * T - 3.67713E-12 * T * T) * HCH * HCH
Cijk(2, 2, 2) = 0.0078498 - 0.000039895 * T + 0.000000061187 * T * T
Cijk(3, 3, 3) = 0.0020513 + 0.000034888 * T - 0.000000083703 * T * T
最低0.47元/天 解锁文章
扫一扫
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
