基于IAOA-PNN模型的天然气压缩因子计算方法研究

采输技术
DOI ：10.3969/j.issn.1001-2206.2023.06.012
基于IAOA-PNN 模型的天然气压缩因子计算方法研究
孙玮
中国石油华北油田分公司检验检测中心，河北任丘062552摘
要：针对部分压气站未设置气相色谱分析仪，无法获取天然气压缩因子的现状，通过拉丁超立方抽样获
取虚拟天然气组分样本，随后以准确度较高的GERG-2008方程为基础，计算天然气密度、热值和压缩因子，形成具有热力学性质的天然气数据库，最后搭建概率神经网络（PNN）模型用于数据的训练、验证和预测，并对预测模型结果进行现场验证。

结果表明，IAOA 算法在收敛速度、训练精度和稳定性上优于AOA 算法、PSO 算法和GWO 算法，证明了算法从种群初始化和密度因子方面进行优化的有效性和科学性；现场校验时本文模型的相对误差维持在-1%~2%之间，且对于组分含量和工况的变化不敏感，可适用于大部分管输气的工况条件；密度是影响算法精度的重要参数，在参数缺失的情况下，应优先保证现场具有监测温度、压力和密度的仪器设备。

关键词：压缩因子；AOA；PNN；GERG-2008方程；相对误差
Research on the calculation method of natural gas compression factor based on IAOA-PNN model
SUN Wei
Inspection and Testing Center of Petrochina Huabei Oilfield Company,Renqiu 062552,China
Abstract：In response to the situation where some pressurized gas stations do not have gas chromatography analyzers and cannot obtain natural gas compression factors,this study used Latin hypercube sampling to obtain virtual natural gas component samples.Subsequently,based on the highly accurate GERG-2008equation,natural gas density,calorific value and compression factor of natural gas were calculated to create a thermodynamic natural gas database.Finally,a probabilistic neural network (PNN)model was built for data training,verification and prediction,and the results of the prediction model are verified on site.The results show that IAOA algorithm is superior to AOA algorithm,PSO algorithm and GWO algorithm in convergence speed,training accuracy and stability,which proves that the algorithm is effective and scientific in population initialization and density factor optimization.The relative error of this model is between -1%and 2%,and it is not sensitive to the change of component content and working conditions,making it suitable for most gas transmission conditions.Density is an important parameter that affects the accuracy of the algorithm,and in cases of missing parameters,it is important to prioritize the presence of equipment and instruments for monitoring temperature,pressure,and density on-site.Keywords：compression factor;AOA;PNN;GERG-2008equation;relative error
天然气作为清洁能源，在我国实现“碳达峰”
“碳中和”的目标中占有重要地位[1]。

在管输气贸易结算的过程中，常使用的计量方式有：体积计量、能量计量和质量计量[2-3]。

虽然能量计量是国际贸易中天然气现货和期权交易的通行手段，但
受我国计量标准、计量技术、配套检测技术和溯源性的影响，我国仍以体积计量为主，即计量单位时间内流过标准孔板的天然气流量。

其中，天然气压缩因子计算结果的准确性是影响贸易交接完备性的重要因素之一。

目前，关于压缩因子的计算有实验法、图版
法和状态方程法等3种。

实验法的结果最为准确[4]，但存在测试周期长、实验费用高的缺点，现场天然气组分多变，实验不可能覆盖所有的工况条件；图版法大多基于Standing-Katz 的标准图版进行回归拟合[5]，但当混合气体中含有重烃组分或较多非烃组分时，误差明显增大，无法满足工程需求；状态方程法以GB/T 17747.2—2011中推荐的AGA8-92DC 方程[6]和ISO 20765-2:2015中推荐的GERG-2008方程为主[7]
，但在计算时需要已知天然
采输技术
气的组分、温度和压力，天然气组分需要通过在线或离线的气相色谱分析仪获取，这对于大多数的中间压气站而言，存在设备费用昂贵的问题，通常压气站并不具备条件。

随着机器学习和计算机技术的发展，人工神经网络已被广泛应用于持液率预测、两相流流型预测和腐蚀速率预测等领域，但在压缩因子计算上还鲜有报道。

基于此，通过拉丁超立方抽样获取虚拟天然气组分样本，随后以准确度较高的GERG-2008方程为基础，计算天然气密度、热值和压缩因子，形成具有热力学性质的天然气数据库，最后搭建概率神经网络（PNN ）模型用于数据的训练、验证和预测，并采用改进的阿基米德优化算法（IAOA ）调整预测效果。

研究结果可为未设置气相色谱分析站点的压缩因子计算提供新的思路和方法。

1PNN 模型
PNN 模型属于前馈性神经网络的一种，其泛
化和学习能力较强，可用于非线性函数的回归和固定模式的分类识别等领域。

不同于其余神经网络，PNN 的求和层采用群体竞争的方式，神经元通过估计不同类别的概率，从而获得竞争输出的机会，最后一次迭代中仅有一个神经元可以输出。

对于模式层，可通过式（1）计算特征向量与训练集中天然气压缩因子的相互关系。

f ()
K ,W i =exp éëêêù
û
úú-()K -W i T ()
K -W i
2δ2
（1）式中：K 为输入特征向量；W i 为输入层与模式层的
权重；δ为光滑因子，是影响PNN 模型预测结果的重要参数。

优化PNN 模型的重点在于寻找最优的光滑因子，确保模型预测结果符合预期。

2IAOA 算法
AOA 算法是由Hashim 等在2020年提出[8]，
与其他元启式优化算法相似，该算法将流体中的物体视为种群，通过不断调整个体的密度、
体积和加速度，实现个体的状态平衡，从而完成寻优过程。

虽然AOA 算法的设置参数少、模
型简单，但仍存在局部搜索能力较弱、寻优精度较低的缺陷。

AOA 算法的原理见文献[9]。

2.1
基于佳点集的种群初始化
在初始化种群阶段，AOA 算法容易受到随机分配的影响，发生群体聚集现象，这种形式不利于提高算法在解空间上的遍历性。

文献[10]采用佳点集的方式获得点分布序列，并将解序列映射到欧式空间。

该方法的构造与AOA 算法的维度无关。

采用随机法和佳点集法在二维种群上的初始化结果见图1。

在同样生成100个种群的条件下，采用佳点集方法产生的初始种群分布更为均匀，可提高搜索效率。

1.00.80.60.40.20.00.00.2
0.4
0.6
0.8
1.0
（a ）随机法
（b ）佳点集法
0.0
0.2
0.4
0.6
0.8
1.0
1.00.80.60.40.20.0图1
不同方法在二维种群上的初始化结果
2.2优化自适应密度因子
在全局搜索和局部搜索阶段，密度因子d 对位
置信息的更新起到关键作用，原始算法见式（2）。

d t +1=exp æèççöø÷÷t max -t t max -æèççö
ø
÷÷t t max （2）
式中：d t +1
为t +1次迭代时的密度因子，t 为当前迭代次数，t max 为最大迭代次数。

采输技术
d 随迭代次数的增加近似呈线性递减趋势，这不利于算法的全局搜索和局部开发。

因此，引入tan 函数进行修正，见式（3）。

对比两种密度因子的变化趋势，见图2。

改进后密度因子先缓慢减小
后快速降低，形成了动态自适应调节极值，有利于保留更多的自身信息，平衡寻优能力。

d t +1=d start -()
d start -d end ×tan éë
êêù
ûúúπω×æèççöø÷÷t t max 2
（3）式中：d start 、d end 分别为迭代开始、迭代结束时的初
值，通过实验验证，得到d start =2.7、d end =0.01的效果较好；ω为控制曲线平滑程度的权值，取ω=4。

图2
密度因子优化对比
迭代次数/次
500
100
150
200
密度因子d 3.02.52.01.51.00.50.0改进前改进后
3
模型构建及分析
3.1
数据来源
首先，调研目前国内管输气的组成（见表1），
利用拉丁超立方抽样的采样方式获取10000个虚拟天然气组分样本，部分数据见表2。

以甲烷为例，具体步骤如下：其一，将（60,100）的区间分成10000段，在每一段中随机抽取一个数值；其二，将抽取的数值通过标准正态分布反函数映射到标准正态分布；其三，利用Matlab 软件中的sort 函数，打乱抽取样本，形成数据集，同时保证所有组分的摩尔分数之和为100%。

此外，根据压气站的实际运行情况，通过抽样方式生成工况温度和压力，温度范围250～400K ，压力小于等于20MPa 。

王国云等[10]、张镨等[11]、王辉等[12]均证明了
GERG-2008状态方程在计算天然气物性上的优越性，故后续采用该状态方程计算压缩因子、密度和热值等热力学参数。

参照ISO 20765.2:2015的公式、方法和流程编制计算程序，求解热力学参数
的迭代过程即为求解一阶非线性方程，常用的方法有二分法、简单迭代法、迭代加速法和Newton 迭代法。

以ISO 20765.2:2015附录G 中的Gas1、Gas2、Gas3作为标准气样，通过对比附录中的实验数据（以压缩因子为例）验证不同求解方法的
准确性，见图3。

其中，二分法的计算结果与实验数据的吻合性较好，其余方法在某些温度和压力条件下出现了较大偏离，故后续采用该方法计算相关参数。

3.2
数据训练及预测
将3.1节得到的数据集按照8∶1∶1的比例分为训练集、验证集和测试集。

训练集用于训练模型超参数，验证集用于预估模型性能和微调模型超参数，测试集用于评估模型的泛化能力和最终性能。

将温度、压力、密度和热值等现场可测量参数作为PNN 模型的输入，将压缩因子作为
PNN
孙玮：基于IAOA-PNN 模型的天然气压缩因子计算方法研究
采输技术
模型的输出，并通过IAOA 算法对光滑因子进行寻优。

在数据训练前，为消除量纲和数据取值范围对预测结果的影响，先对数据进行离散标准化处理，将其映射到[0,1]区间。

为验证IAOA 算法的优越性，将其与AOA 算法、PSO （粒子群）算法和GWO （灰狼）算法进行对比，以均方根误差（RMSE ）为适应度函数，对比不同算法在训练过程中的适应度曲线变化情况，见图4。

其中，各算法的种群规模为30，空间维度为1，最大迭代次数为100，各算法的基本参数取值参照文献[9]。

IAOA 算法在迭代至25次时即达
到收敛状态，适应度值最小为0.00258；PSO 算法的适应度值呈先降低后上下波动的趋势，说明该算法已陷入局部最优解的困扰，增加迭代次数才有收敛的可能；AOA 算法和GWO 算法分别在59次和61次时达到收敛状态，但适应度值较IAOA 算法相比明显偏大。

综上，IAOA 算法在收敛速度、训练精度和稳定性上优于其余算法，证明了算法改进的必要性和合理性，此时的光滑因子取值为0.65218。

图4
不同算法的适应度曲线
迭代次数/次20406080100
0适应度函数值
0.250.200.150.100.050.00
IAOA AOA PSO GWO
模型的训练、验证及预测结果见图5。

压缩因子在训练集、验证集和测试集上的决定系数分别为0.9997、0.9989、0.9984，RMSE 分别为0.00258、0.00756、0.00851，所有数据均在回归线附近分布，只有少量压缩因子较低的数据发生了偏移，但偏移量较小，这种现象可能与虚拟组分中重烃含量较高或温度接近临界凝析温度等因素有关。

实际值模型值0.6
0.7
0.8
0.9
1.0
1.00.90.80.70.6（b ）验证集
实际值模型值
0.6
0.7
0.8
0.9
1.0
1.00.90.80.70.6（c ）测试集
实际值模型值0.6
0.7
0.8
0.9
1.0
（a ）训练集
1.00.90.80.70.6图5模型的训练、验证和预测结果
温度/K
240260280300320340360380400420
压缩因子实验值二分法
简单迭代法迭代加速法Newton
1.000.950.900.850.800.750.700.650.60
图3
不同求解方法对比
采输技术
3.3现场校验结果
选取某输气管道上的3个压气站为例进行模型
校验，温度测量采用Rosemount 214C 型传感器，不确定度±0.1℃；压力和流量测量采用Rosemount
2051CFP 型一体化孔板流量计，不确定度最高0.9%；密度测量采用横河GD402G 型密度分析仪，当温度波动不超过10℃时，不确定度最高
0.001kg/m 3；天然气热值通过SIEMENS SITRANS-CV 型热值分析仪获取。

将以上可测量参数输入训练好的IAOA-PNN 模型获得工况下的压缩因子，同时结合工况流量Q L ，计算标况流量Q b ，以获取实时在线的天然气流量，现场布置的天然气测量系
统见图6。

图6天然气测量系统示意
热值
密度
压缩因子和标况流量
天然气出口节流阀
流量计
空冷器
天然气进口
P
T
IAOA-PNN 模型修正因子
将气相色谱分析仪得到的天然气组分和温度、压力作为已知参数，代入GERG-2008状态方程求解压缩因子，将该值作为实际值，与模型值进行对比，结果见图7（只标出了前三个最大的组分含量）。

实际值与模型值的曲线一致性较好，两者的吻合程度较高；随着压力的增加，压缩因子先减小后增大，在增大区间，压缩因子与压力呈线性关系，压缩因子随温度的增加逐渐增大，但增大趋势变缓，预测结果与Standing-Katz 标准图版描绘的变化趋势相符，再次证明了本文算法的有效性。

此外，1#、2#、3#压气站的分子量分别为16.356、17.793和19.091，虽然分子量逐渐增大，组分中除C 1以外的重烃含量逐渐增多，但从箱线图（见图8）中可知，相对误差并没有明显增大，不同组分和温度下的相对误差始终维持在-1%~2%之间，说明IAOA-PNN 模型对于组分含量和工况的变化不敏感，可适用于大部分管输气的工况条件。

3.4
不同可测量参数对预测结果的影响
考虑到现场实际工况中的温度、压力数据已并入SCADA 系统中，而密度和热值并非必要测试
项目，存在一定的数据缺失风险。

将本文模型定为IAOA-PNN-1模型，将包含温度、压力和热值数据的模型定为IAOA-PNN-2模型，将包含温度、压力和密度数据的模型定为IAOA-PNN-3模型。

以1#压气站为例，观察不同温度下的相对误差变化，见图9。

IAOA-PNN-2模型的相对误差范围与IAOA-
PNN-1模型相比略微增大，IAOA-PNN-3模型的
相对误差范围较其余模型有大幅增加，说明密度
压力/MPa 101520
（a ）1#压气站
5
压缩因子
GERG-2008（250K ）IAOA-PNN （250K ）GERG-2008（300K ）IAOA-PNN （300K ）GERG-2008（350K ）IAOA-PNN （350K ）
98.54%CH 4+0.647%C 2H 6+0.5%N 2
1.051.000.950.900.850.800.750.700.650.600.550.50
（b ）2#压气站
90.04%CH 4+3.69%C 2H 6+4.48%N 2
压力/MPa
101520
5
压缩因子
GERG-2008（250K ）IAOA-PNN （250K ）GERG-2008（300K ）IAOA-PNN （300K ）GERG-2008（350K ）IAOA-PNN （350K ）
1.051.000.950.900.850.800.750.700.650.600.550.50
85.01%CH 4+9.38%C 2H 6+3.49%C 3H 8
压力/MPa
10
1520
05
压缩因子
GERG-2008（250K ）
IAOA-PNN （250K ）GERG-2008（300K ）IAOA-PNN （300K ）GERG-2008（350K ）IAOA-PNN （350K ）
1.00.90.80.70.60.50.40.3
（c ）3#压气站
图7
不同压气站上的现场校验结果
孙玮：基于IAOA-PNN 模型的天然气压缩因子计算方法研究
采输技术
250K 300K 350K 250K 300K350K 250K 300K 350K IAOA-PNN-1相对误差/%
543210
-1-2-3
IAOA-PNN-2IAOA-PNN-3
图9
不同参数模型下的压缩因子相对误差箱线图
是影响算法精度和稳定性的重要参数。

在参数缺失的情况下，应优先保证现场具有监测温度、压力和密度的仪器设备。

4结论
1）针对部分压气站未设置气相色谱分析仪，
无法获取天然气压缩因子的现状，以温度、压
力、密度和热值为可测量参数，以压缩因子作为未知参数，通过拉丁超立方抽样方式构建虚拟天然气组分样本，将样本代入IAOA-PNN 模型进行训练和预测，得到一种简便易行的压缩因子计算方法。

2）IAOA 算法在收敛速度、训练精度和稳定性
上优于AOA 算法、PSO 算法和GWO 算法，证明了算法改进的有效性和科学性。

3）通过在压气站进行现场校验，本文模型的
相对误差维持在-1%~2%之间，且对于组分含量和工况的变化不敏感，可适用于大部分管输气的工况条件。

当参数缺失时，应优先保证现场具有监测温度、压力和密度的仪器设备。

参考文献
[1]熊治富，邓金花，陈丹.川东北地区天然气偏差系数计算方法
评价[J].西南石油大学学报（自然科学版），2015，37（4）：60-66.
[2]王辉，周柳玲，陈晓科.高含硫原料天然气贸易交接计量方法
研究[J].工业计量，2023，33（3）：10-14.
[3]王艺晨，叶继根，吴淑红.基于对勾函数的新型天然气压缩因
子计算关联式[J].石油与天然气化工，2021，50（5）：38-43.
[4]SOAVE G.Equilibrium constants from a modified Redlich-Kwong equation of state[J].Chemical Engineering Science，1972，27（6）：1197-1203.
[5]KATZ D L，CORNELL D，ROBAYASHI R，et al.Handbook
of natural gas engineering[M].New York：McGraw-Hill Book Company，1959.
[6]GB/T 17747.2—2011，天然气压缩因子的计算第２部分：用
摩尔组成进行计算[S].
[7]UNE-EN ISO 20765-2:2019，Natural gas -Calculation of
thermodynamic properties -Part 2:Single-phase properties
（gas，liquid，and dense fluid）for extended ranges of ap-plication[S].
[8]HASHIM F A，HUSSAIN K，HOUSSEIN E H，et al.Archi-medes optimization algorithm：a new metaheuristic algo-rithm for solving optimization problems[J].Applied Intelli-gence，2021，51（3）：1531-1551.
[9]罗仕杭，何庆.融合Sin 混沌和分段权值的阿基米德优化算法
[J].计算机工程与应用，2022，58（14）：63-72.
[10]王国云，李长俊，贾文龙.基于GERG-2008方程对CH 4-H 2
混合物压缩因子的计算与因素分析[J].天然气化工-C1化学与化工，2022，47（2）：89-94.
[11]张镨，周理，陶建，等.物性值法计算天然气压缩因子适应性
分析[J].石油与天然气化工，2023，52（2）：48-54，64.[12]王辉，李木盛，陈荟宇.准确度更高的高含硫天然气压缩因子
计算方法[J].石油与天然气化工，2022，51（6）：50-60.
作者简介：
孙
玮（1987—），女，河北任丘人，工程师，2010年毕业于长江大学油气储运专业，现主要从事检验检测技术工作。

Email：***********************.cn
收稿日期：2023-09-16；修回日期：2023-10-23
250K 300K 350K 250K 300K 350K 250K 300K 350K
11#压气站2#压气站3#压气站相对误差/%
2.01.51.00.50.0-0.5-1.0
图8
不同压气站的压缩因子相对误差箱线图。