基于SHBWR状态方程求解城镇燃气压缩因子的计算方法

随着城镇燃气气源的天然气化和管网的大型化,燃气管网水热力计算时为满足计算精度要求必须考虑燃气的压缩性。本文采用基于SHBWR状态方程求解城镇燃气的压缩因子,数值计算方法求解以燃气密度为自变量的非线性方程。本文通过分析燃气密度函数在邻近有效解的范围内曲线的变化趋势,选用牛顿迭代法、弦截法及二分法求解该方程,并编写相应的计算程序。通过算例分析,牛顿迭代法和弦截法是较好的求解方法,选取适当的初值时,计算很少的次数(算例中约为2-6次)就可以达到同样的精度。