第1页共4页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第1页共4页燃气管网的邻接矩阵描述及相应的水力计算算法李其1，祝达2(1.天津市煤气工程设计院，天津市300381；2.天津大地燃气技术开发有限公司，天津市300061）摘要：介绍了用邻接矩阵描述燃气管网的方法和对应的水力计算算法，并总结了用邻接矩阵描述燃气管网方法的特点。关键词：燃气管网；水力计算；邻接矩阵；节点法燃气管网的水力计算一般采用节点法，节点法首先要有对燃气管网的准确描述，然后构建相应得数学模型进而求解。本文就对用邻接矩阵描述燃气管网的数学模型及其求解算法进行介绍和分析。1.燃气管网的描述用节点法描述燃气管网需要确定如下参数：pi---各节点压力qi---各节点负荷Li,j---各管段管长Di,j---各管段管径图1表示一个由6个节点和6个管段构成的燃气管网，其节点已编号，其中节点1，4为已知压力节点（气源点），节点2，3，5，6为已知负荷节点。各节点压力可表示为向量P=(p1,p2,,…,p6)T各节点负荷可表示为向量Q=(q1,q2,…,q6)T从图1可以看出燃气管网可以描述为一个连通的无向图，则各管段管长和管径可分别表示为图的邻接矩阵第2页共4页第1页共4页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第2页共4页L=(L1,1L1,2¿⋅¿L1,6L2,1L2,2¿⋅¿L2,6¿⋅¿¿⋅¿¿⋅¿¿⋅¿L6,1L6,2¿⋅¿L6,6)和D=(D1,1D1,2¿⋅¿D1,6D2,1D2,2¿⋅¿D2,6¿⋅¿¿⋅¿¿⋅¿¿⋅¿D6,1D6,2¿⋅¿D6,6)因为管段无方向所以矩阵L，D均为以主对角线为中心的对称矩阵，因此可用L，D的半三角矩阵来描述管网，可利用矩阵的对称特性将两矩阵合并为矩阵LD=(LD1,1L1,2L2,3¿⋅¿L1,6D2,1LD2,2L2,3¿⋅¿L2,6D3,1D3,2LD3,3¿⋅¿L3,6¿⋅¿¿⋅¿¿⋅¿¿⋅¿¿⋅¿D6,1D6,2D6,3¿⋅¿LD6,6)LD主对角线上的元素LDii无实际意义，可规定其值为0。两节点不连通则规定相应得元素值为0，例如图1中节点1，4不连通则有：L1,4=0D4,1=0通过上面的分析可知图1所示的燃气管网可用向量p，q和矩阵LD来描述，很显然，这种描述方法可推广至有n个节点的任意燃气管网。由燃气输配的连续性和能量规律可知向量p，q和矩阵LD是相关的，即已知其中任何两个量就可以计算确定第三个量，因此只需向量p，q和矩阵LD的任何两个量即可确定地描述一个燃气管网。2.邻接矩阵节点法的算法燃气管网中，对于任意节点流量是连续的，流出节点的流量=流入节点的流量，即有连续性方程：ΣQi=0（1）对于有n个节点的燃气管网，可列n个连续性方程，可解得n个未知数。应用节点法时将力量和压力降的非线性关系：ΔP=aQn（2）转化成线性关系：Q=CΔP（3）式中C为系数，即C=QΔP（4）上式中Q与ΔP均为绝对值，节点流量的符号规定如下：流向节点的流量为负，因为该管段试点的压力大于该节点的压力；流离节点的管段流量为正，因为该管段在节点处的第3页共4页第2页共4页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第3页共4页压力大于其终点的压力。压力降值的符号与流量相同。对于任意节点i应用公式（3）和连续性方程（1）：Ci,1(pi−p1)+Ci,2(pi−p2)+¿⋅¿+Ci,n(pi−pn)+qi=0展开为：(−Ci,1)p1+(−Ci,2)p2+¿⋅¿+(Ci,1+Ci,2+¿⋅¿+Ci,n)pi+¿⋅¿+(−Ci,n)pn=qi构造系数矩阵C=(C1,1+C1,2+¿⋅¿+C1,n−C1,2¿⋅¿−C1,n−C2,1C2,1+C2,2+¿⋅¿+C2,n¿⋅¿−C2,n¿⋅¿¿⋅¿¿⋅¿¿⋅¿−Cn,1−Cn，2¿⋅¿Cn,1+Cn,2+¿⋅¿+Cn,n)则节点法求解即转化为解线性方程组Cp=q（5）对于有n个节点的燃气管网，设气源节点数量为m个，由于气源点的节点压力已知，而节点负荷未知则向量p中有m个分量为已知，n-m个分量未知，而向量q中则有m个分量为未知，因此线性方程组（5）中恰有(n-m)+m=n个未知数，又因为该线性方程组线性无关，所以对于任意燃气管网方程组（5）有惟一解。开始计算时，系数矩阵可按燃气管网各管段的雷诺数确定，于是初拟的各管段流量为Q=WF=ReνdF=Reπdν4雷诺数Re可任取一个数值，但为使迭代算法快速收敛，应取尽量接近实际的数值。对于低压燃气管网可取Re=2×104，对于高中压管网可取Re=3×105。再由式（2）求Δ...
