吕秦，杨惠德

（四川染料厂，重庆 401256）

关键字：标准；天然气；压缩因子；密度；Octave软件；IEC 61131-3

1 为什么要对液化天然气压缩因子计算进行仿真

根据ISO 12213-2国际标准，我国对天然气压缩因子的计算颁布了相应的国家标准，即GB/T 17747.2[1]。该标准等效采用ISO国际标准，但删除了国际标准中附录F的有关程序实现部分。经与国际标准的程序对照，可以发现，国家标准中的（B.6）公式中累加项前的系数为1，而对照国际标准ISO 12213-2：2009[2]附录F的程序，其系数为2，因此，有必要确定其系数实际值，即需要对天然气压缩因子的计算进行仿真，确定其系数值。

国际标准ISO 12213-2附录F提供了Fortran编程语言的程序，但目前大多数科学计算都采用MATLAB实现。由于美国对我国有关单位使用MATLAB的限制，因此，需要有对应的仿真软件用于对天然气压缩因子计算进行仿真。为此，本文采用Octave软件进行天然气压缩因子的计算仿真。Octave软件是开源软件，其源于MATLAB。为检验其可替代性，本文将MATLAB程序和Octave程序进行了比较，实际结果显示，大多数MATLAB的语言可直接用于Octave程序，其主要区别是MATLAB程序中，一些关键字可用缩写，而Octave程序中，虽然可以使用，但会提示警告。例如，fonts关键字，在Octave程序中应列出全称fontsize。此外，MATLAB程序中的for…end循环语句，在Octave程序中应表示为for…endfor等。而且程序执行速度也比MATLAB的程序执行要慢些。

2 天然气压缩因子计算的仿真

根据国际标准附录F提供的Fortran程序编写了Octave程序，该程序采用友好的用户界面，可根据天然气的压力和温度以及各组分的体积百分比，计算该工况下天然气的压缩因子和密度。……