[发明专利]一种压缩因子计算方法

权利要求书

1.一种压缩因子计算方法，其特征在于，包括训练阶段和测试阶段。

2.如权利要求1所述的压缩因子计算方法，其特征在于，所述训练阶段包括以下步骤：

S1，利用AGA8号报告获得1200组不同气体组分、不同工作温度和不同工作压力的压缩因子数据；

S2，计算训练样本的混合气体临界压力与临界温度；

S3，计算训练样本的混合气体拟对比压力与拟对比温度；

S4，将Ppr与Tpr以式(5)为分界线划分成两部分；

P_pr＝-18.2882T_pr³+32.3514T_pr²-14.9162T_pr+2.0352 (5)

S5，建立GMDH网络输出输入的关系式；

S6，将式作为目标函数，分别训练已被划分的两部分；

S7，训练得到两个GMDH网络，利用AICc优化网络模型；

AICc＝(2k-2L)/n+2k(k+1)/(n-k-1) (8)

其中k代表模型参数数量，L代表对数似然值，n代表网络节点数。

S8，训练完成。

3.如权利要求2所述的压缩因子计算方法，其特征在于，所述压缩因子数据，其中70％作为训练样本，30％作为测试样本。

4.如权利要求2所述的压缩因子计算方法，其特征在于，所述混合气体临界压力与临界温度的计算公式为；

其中xi为气体摩尔组分，Pci与Tci分别为对应的纯气体临界压力与临界温度。

5.如权利要求2所述的压缩因子计算方法，其特征在于，所述混合气体拟对比压力与拟对比温度的计算公式为；T_pr＝T/T_c，P_pr＝P/P_c；其中Pc与Tc分别为混合气体的临界压力与临界温度，T与P分别为管道内气体的工作温度与工作压力(绝压)，Ppr与Tpr分别为混合气体的拟对比压力与拟对比温度。

6.如权利要求2所述的压缩因子计算方法，其特征在于，所述GMDH网络输出输入的关系式为Z(x,y)＝a_x,yx²+b_x,yxy+c_x,yy²+d_x,yx+e_x,yy+f_x,y，其中x和y代表网络节点输入，Z代表网络节点输出，a_x,y,...,f_x,y代表多项式系数。

7.如权利要求1所述的压缩因子计算方法，其特征在于，所述测试阶段包括计算测试样本的P_pr与T_pr。

8.如权利要求7所述的压缩因子计算方法，其特征在于，所述测试阶段还包括将测试样本划分P_pr与T_pr成两部分。

9.如权利要求7所述的压缩因子计算方法，其特征在于，所述测试阶段还包括将划分后的测试样本P_pr传入对应部分进行测试，得到对应的P_pr压缩因子估计值。

10.如权利要求7所述的压缩因子计算方法，其特征在于，所述测试阶段还包括将划分后的测试样本T_pr传入对应部分进行测试，得到对应的T_pr压缩因子估计值。