复杂燃气管网水力计算..docVIP

复杂燃气管网水力计算算法 Complex Gas Pipe Network Hydraulics Computation Method 苏州燃气集团有限公司 胡小龙 [摘要：] 本文论述了在管网网络图的基础上，用深度优先搜索遍历算法及线性化流量平衡方法进行多气源(同一种燃气,不同位置)、多环复杂管网的水力计算，此算法可用于中压、低压管网设计型问题及操作型问题计算(实际管网最大负荷计算)。 关健词：管网 水力 算法 流量 压力 一、算法概要： 本算法程序是以微软VC++6.0作为开发工具,以图作为系统数据结构,存贮系统相关数据。本算法程序运用的前题是首先录入管网基础数据构造管网网络图，即录入气源结点（开始结点）、普通结点（内部结点）、调压站或用户结点（终点结点）的结点座标（可用以计算边长度）、初始管径、壁厚等。管网图的形成，也即完成了管网水力算法所需数据的录入工作，由于网络图的可视性，不但可以用来检验数据录入的正确性，而且可用来初步检查管网结点计算压力和各边流量的准确性。形成管网图及水力计算所需参数后，用深度优先搜索遍历算法从气源点开始遍历管网所有结点，进行流量及压力的初始化，然后进行各结点流量平衡计算，直至各调压站站前压力及输出流量达计算要求。 二、水力算法模型及分析： 1、水力计算模型： 本算法的关健是管网各结点流量平衡计算,为了加快运算收敛速度,流量平衡算法采用线性化法,线性化法如下所述:(以中压管网,铸铁管为例) 中压管网水力计算公式为: Pa2-Pb2=127*λ*(i+1)*L*ρ*Q2/D5………….[1] (Pa、Pb:Kpa , L:cm,ρ:Kg/m3,Q: m3/h,D:cm,i为局部阻力当量系数,燃气处于标准状况,密度取平均值) λ=[2.089*10-7/D+2.3633/Re]0.284 i取10%-20% (依管网支点疏密等条件而定)。 Pa-Pb=Q*(127*λ*(i+1)*L*ρ*Q/D5)/ (Pa+Pb) 令: Cab= 127*λ*(i+1)*L*ρ*Q/D5)/ (Pa+Pb) 则: Pa-Pb= Cab*Q，即 Cab=（Pa-Pb）/Q--------------------------------------[2] 在进行循环计算时,Cab中的参数用上一`次计算的数值(主要指流量Q及结点压力Pa、Pb) 各结点（非终点结点）流量平衡时: ΣQji=0 结点压力为： Pi=（ΣPj/Cij）/（Σ1/Cij） 对于同一结点，二次相邻计算的压力值偏差相对误码差不得大于规定值一般为10-6左右。只有这样，才能保证公式[2]中等式左边参数中的Q及Pa、Pb同等式右边中Q及Pa、Pb 在误差范围内基本为同一值，流量平衡线性化方法才能保证同中压水力计算式相一致，也只有这样，才能保证各结点流量在误差范围内真正平衡。 计算流量时，阻力系数λ=[2.089*10-7/D+2.3633/Re]0.284 中中的Q用上一`次计算的Q值,由于本算法各结点压力及各边流量不断向目标值逼近,因此,此法计算所得阻力系数误差可忽略。 2、算法分析： 管网水力计算时，物性参数、粗糙度、管长、用气负荷一般为已知（可认高峰流量为已知），未知参数为各边管径D、流量Q、阻力系数λ及各结点压力P（气源始点压力可认为已知，如输送系统与管网一起进行优化计算，则气源始点压力为未知）设管网边数为M，结点数为N（终点结点数为N0），那么，系统变量为：3M+N（个）,系统方程数为：各边压差方程 M（个）、阻力系数方程M（个）结点流量平衡方程N-N0（个），因此，系统自由度为：3M+N-(2M+N- N0)=M+ N0 而这M+ N0个自由度一般为M个管径及N0-1个调压站高峰输出流量(即系统最大输送流量)及气源始点压力。以上M+ N0个自由度确定后，系统状态也就确定了。进行各结点流量平衡计算时，系统需初始化结点压力及各边流量，本算法流量初始化方法如下：各边最初流量初定为其下游各用户结点流量，然后，从气源点出发，将各结点上游各边流量按各结点下游各边初始流量所占比例重新分摊到各结点下游各边中。然后采用深度优先搜索遍历算法,从气源结点开始遍历各结点进行压力初始化及流量平衡计算。 三、算法流程： 流程图如图一所示。 2、流程简述： 管网基础数据录入完毕形成管网图后,初始化水力计算所需要的基础数据，主要为物性数据（密度、粘度），管网流体传输特性数据（粗糙度、局部阻力当量系数），管网负荷数据（气源厂高峰出口压力、各调压站用户数、户