图论法用于供水管网水力计算的研究.doc

图论法用于供水管网水力计算的研究河南省黄河河务局 张丰周武汉水利电力大学 魏永曜天津市大港油田集团公司设计院 石继摘 要 图论理论是网络分析的主要工具，现用于管网的水力平衡计算 ，既充分发挥了图论理论的优势，使计算变得简便、迅捷，又可将管网附件加入计算，使结果更准确、更符合实际。文中采用峰阵输入管网结构，使输入数据的工作量大大减少，易于编制程序，计算大型的复杂管网。关键词 供水管网，水力计算，图论法。前 言供水管网的水力平衡计算是供水系统规划设计、经济评价和运行管理的基础。水力平衡计算的目的就是在确定管径的情况下求出满足连续方程和能量方程的各节点压力水头和各管段流量。目前常用的水力平衡计算方法有哈代-克罗斯法(Hardy-Cross)，牛顿-莱福逊法(New ton-Raphson)，线性理论法(Linear-Theory)，有限元法(Finite Element)等等。所有这些方法各有所长，适用范围各不相同，有的还需人工假设管段流量，使输入数据工作量增大，且未考虑管网附件的影响。本文介绍的图论法将复杂的管网处理为相应的“网络图”，并建立相应的数学模型，用峰阵输入原始数据来描述管网结构，输入的数据量最少，不易出错，易于计算大型的复杂管网。其计算过程可同时考虑管网附件，如控制阀、加压泵、逆止阀、减压阀等，使计算结果更符合实际。1 图论原理将供水管网中的管段概化成一条线段(即图中的边)，将有附件的管段看成图中的特殊管段，边与边由节点相连。这样，一个供水系统的管网图就转化为图论中的网络图。而且管道中的水流是有方向的，所以管网图是有向图。根据以上所述原则，可将图1所示管网系统，转化为图2所示的网络图。 图1图2图1中有一水库A，三个给水点B、C、D，Q1表示水库节点供水量，Q2,Q3,Q4分别表示B、C、D节点的用水量。管段视为网络图中的对应边，管段的直径、管长、管道流量、摩损系数等作为管段对应边的权。至此，与管网同构的网络图生成了。图中箭头表示各条边的方向，即管段中水流方向。网络图中节点与边的关联函数可以用完全关联矩阵I45表示如式(1)所示。顶点边的编号(1)式中：Iij=1，表示j管段与i节点相连，且管内水流流离该节点；0，表示此管段不与该节点关联；-1，表示j管段与i节点相连，且管内水流流入该节点。完全关联矩阵与管段流量列向量q以及节点流量列向量Q可组成管网节点方程(即连续方程)Iijq+Q=0,q=(q1,q2,q3,q4,q5)T,Q=(Q1,Q2,Q3,Q4)T。网络图的生成树(全涉及树)可以有很多种，在计算时可以任选一种。在本例中，选1、2、4这3条边为图的生成树，则补树(余树)的各边(弦)为3、5.各弦将与枝构成基本回路，一个基本回路中有且仅有1条弦。用基本回路矩阵Bf表示则如式(2)所示。 枝124弦35Bf=-11010(2)1-1-101式中每一行表示一个基本回路(环)。环的方向以该环对应弦的方向为准。“-1”表示管段中的流向与环中弦的方向相反，“1”表示相同，“0”表示该管段不在此环内。Bf可用矩阵B和单位阵U表示为式(3)。Bf=BU,其中B=-110(3)1-1-1环阵与管段摩损列向量hf构成环方程如式(4)所示。摩损向量的元素顺序与Bf中每行元素所对应的管段顺序相同。Bfhf=0。其中hf=(h1,h2,h4,h3,h5)T(4)图论理论中，连续方程用割方程代替。每个割方程只含一根枝，并和相关的弦构成割集，将图2分割成互不连通的脱离体。这样，图中就有3个割集。割集和割集阵Af如式(5)所示：割集K：割阵：Af=枝124弦35(5)K1=(e1,e3,e5)1001-1K2=(e2,e3,e5)010-11K3=(e4,e5)00101割阵Af中，每一行表示一个割集。图中有3根枝，所以就有3个割集。割阵中，“+1”表示该管段在此割集内，且管段流向与此割集内的枝中的流向相同，“-1”表示流向相反，“0”表示该管段不在此割集内。式(5)的割阵Af和割集K一一对应。割阵Af可用一个矩阵A和一个单位阵U表示为：Af=UA,其中A=1-1-1101割阵与流量列向量可构成割方程。根据图论理论，割阵的行向量与环阵的行向量正交，这种关系可用式(6)表示。B|UU|AT=0或者U|AB|UT=0(6)所以有B=-AT或者A=-BT。这样，环阵可以由割阵求出，反之亦然。关联矩阵通过选主元初等行变换即可得到割阵：先选关联阵第一行中一非零枝元素为主元，并使其为+1，消去其它各行中此主元；再选第二行、第三行、的主元，最后即得割阵Af。因此，可以由关联矩阵导出割阵和环阵。2 图论法模型任何管道的水力计算都可以用管段流量q,水头损失h,管径D,管长L和管壁条件C等5个因素来描述。一般D、L和C为已知条件，只有q和h未知。因此，求解一个管网的水力平衡问题，可从两方面考虑：一是利用q和h的关系，消去h，以q为未知量计算，求出q后，反求h；二是首先消去q，以h为未知量计算；解出h之后，再反求。图论法也可从这两方面入手，即求弦流量式和求枝摩损式。前者只适用于环状网，而后者则适用于所有类型的管网，所以本文着重介绍后者。设一管网有J个节点，P条管段，L个环，则三者满足L=P-J+1的关系。管网的每一管段都有q和h两个未知量，因而未知量的个数为2P。但管网环方程有L个，线性无关的连续性方程有J-1个，总数为L+J-1=P个，不能求解2P个未知量。因此，必须借助P个管段摩损方程式。管段摩损方程式线性化后的通式如(7)和(8)所示。系数R称为阻尼系数，Y称为传导系数。R和Y的具体形式与所选用的摩损公式有关，是D、C、L的函数。摩损公式线性化后，R还是q的函数，Y还是h的函数。不过，在求解过程中，总是把R和Y当作已知量来对待。阻尼式：h=Rq(7)传导式：q=Yh(8)式中R和Y是阻尼系数和传导系数矩阵。如果摩损公式采用Hazen-William公式，则有：h=Rq=10.68q1.852L/(C1.852D4.87)=10.68L|q|0.852/(C1.852D4.87)q(9)R=10.68L|q|0.852/(C1.852D4.87)(10)Y=1/R=C1.852D4.87/(10.68L|q|0.852)=C1.852D4.87/(10.68L)|q|-0.852(11)用h向量表示管段摩损：h表示枝摩损，h表示弦摩损；用q向量表示管段流量：q枝管段流量，q表示弦管段流量。割方程的右端项Q为脱离体所含节点流量之和。方环程：Bfh=0,即 B Uh=0(12)h割方程：Afq=Q,即 UAq=Q(13)q传导式：q=Y0h(14)q0Yh求枝摩损式(以管段摩损为未知量)：首先将传导式(14)代入割方程(13)得：U AY0h=Q(15)0Yh由环方程(12)可得Bh+h=0，即h=-Bh，代入式(15)得：U AY0h=Q(16)0Y-Bh即hY-AYB=Q(17)根据正交定理得：hY+AYAT=Q(18)这就是图论法的求枝摩损式计算公式。h即为枝管段的摩损向量。解得枝摩损值h后，其余变量可由相应的公式求出。由环方程可得h=-Bh，即可求出弦摩损向量h,q、q向量可以由式(14)求得。式(11)中C1.852D4.87/10.68L对某一管段来说是个常数，可用W表示。则传导系数Y可以表示为：Y=Wq-0.852(19)在迭代计算时，第一次可以直接用W代替Y进行计算，求出h,q后计算Y，再求新的q值，如此反复计算，直至前后两次的q值符合给定的误差标准为止。为了避免可能出现的数值摆动现象，在第三次迭代时，用前两次迭代结果的流量平均值作为初始流量值，即：q=q（）+q(2)2(20)求得q(3)，这样收敛速度加快。3 管网附件实际管网中，有许多控制、安全、量测设施，如加压泵、控制阀、逆止阀、减压阀等附件，对管网运行产生重要影响。传统计算方法都未涉及到管网附件问题，不仅使计算准确性受损，而且其计算程序无法用于日常管理工作。图论法处理管网附件时，将附件所在管段视为特殊管段，这些管段的摩损式要根据其附件的水力学特征计算摩损值，再加入到管网中进行水力平衡计算。本文给出几种较常见管网附件的处理方法。对于其它附件，具体问题具体处理，在此就不一一详述了。3.1 普通阀门 闸板式阀门是用得最多的一种阀门，在一 般的水力计算过程中，闸板式阀门的水头损失计算一般引用公式hf=v2/2g,值见文献3。其中，a表示管段中过水断面的高度，d表示管段直径，a/d表示阀门开关。当开度为0时，阀门完全关闭，没有流量通过；当开度为1时，阀门完全打开，对水流不产生影响。将阀门水头损失公式用流量表示为：hf=v2/2g=2q2/2gD2则阻尼系数R为：R=2q/2gD2;传导系数为：Y=2gD2/2q-1计算时只需将闸板式阀门的R或Y值加入，即可计算。蝶阀的计算方法与闸板式阀门类似3.2 逆止阀 逆止阀是管网中最常见的设备之一，是水流方向控制设备，只允许水流单向通过。设第K根管段上装有一个逆止阀，管段的传导系数为yk。并且，逆止阀的工作方向为从节点i到节点j，如图3(a)所示。则逆止阀的工作状态可以分为以下两种：1当节点水压HiHj，则水流可从K管段中流过。此时就相当于一根普通管段，逆止阀 就象不存在一样。管段的传导系数就是普通管段的传导系数，如图3(b)所示。图32当节点水压HiHj，逆止阀将处于工作状态，阀门自动关闭，水流无法通过，如图3(c)所示。相当于这条管段不存在，yk=0，流量qk=yk(hk)n=0。3.3 减压阀 减压阀(Pressure-Reducing ValvesPRV)可控制管网中该阀门下游端压力水头值保持在某一范围内，而不致压力太高，且兼有逆止阀的作用。减压阀有一额定工作压力，当其上游压力超过工作压力时，其下游端会维持一恒定的压力。若上游压力低于工作压力，则减压阀不起作用。若下游水头高于上游，PRV相当于一处于工作状态的逆止阀，管段流量为零。图4如图4所示，上游节点i,下游节点j，管长L，下游端管长，整个管段的传导系数yk，PRV下游端管段的传导系数，上游端水头值He，额定工作水头Hset,管段上游节点水头值Hi，管段下游节点水头值Hj。减压阀的工作状态如下： 1当HeHset且HiHj，减压阀对水流没有影响，可视为普通管段。传导系数就是yk,如图5(a)所示。图42当HeHset且HiHj，减压阀处于工作状态，下游端水压力等于Hset管段流量qk是水头Hset和Hj以及传导系数的函数。这种情况就好象在减压阀处有一固定水头为Hset的水库在向供水，如图5(b)。管段流量可表示为：qk=-(Hset-Hj)n。3当HiHj，减压阀就相当于一个逆止阀，管段中没有流量通过，传导系数为零，yk=0,如图5(c)。图53.4 加压泵 由于水泵的能量输入，水流流过水泵之后，水头不仅未象普通管段那样减少，反而增加了，这就是水泵扬程。但也可把扬程视为“管段”的“水头损失”，只是这“水头损失”是负值。加压泵一般选用离心泵，其流量和扬程关系可由水泵特性曲线H=Hx-SxQ2表示，Hx和Sx由水泵特性曲线求出。对水泵来说，有实际意义的是落在水泵高效区的一小段曲线。而管道计算常用传导公式一般形式为q=h,、为常数。因此可用适合管网计算的曲线(21)拟合水泵特性曲线，使Q可用一常数项加一个关于H的非常数项来表示。q=q0-ah2(21)在水泵特性曲线高效区靠近两端的位置各选一具有代表性的点坐标值，代入式(21)，即可求 出系数a和常数项q0。这样，拟合曲线就确定了。水泵特性拟合曲线式(21)由于有一常数项q0，仍无法直接给出此“管段”的传导系数。因此可把式(21)分成2部分：第1部分为q=q0，是一个恒定流量；第2部分为q=-ah2，是一与扬程H有关的流量.所以水泵这一特殊“管段”可视为两根管段的复合，即一管段中的流量始终是q0，另一管段的传导系数为“-a”。可在水泵上游节点加上一个值为“-q0”的节点输入流量，下游节点加上一个值为“q0”的节点输出流量。水泵管段的传导公式就变为式(22)，即可按普通管段计算。计算结果加上流量q0即得水泵流量。 q=-ah2(22)3.5 水塔 水塔节点的压力水头可认为是恒定值，即H(i)=constant。当对置水塔向管网供水时，可以把水塔的供水量作为节点入流量，纳入管网进行计算。4、算例 作者用此法计算武汉市汉阳区自来水管网，得到了较满意的结果。该系统有330个节(其中3个水源节点)，395根管段，66个环，见图6。当不考虑附件时，部分计算结果如表所示，经验证，结果完全满足节点方程和环方程。若考虑节点56、57之间有一逆止阀，只允许从节点57流向56。考虑附件后的部分计算结果如表所示。图6 武汉市汉阳区自来水水管网简图 无附件有附件管段上游下游流量流速上游节下游节流量流速上游节下游节编号节点节点点水头点水头点水头点水头m3/Sm/smmm3/sm/smm11220.5840.7466.2966.090.5840.7466.2966.092121.1801.0466.2966.191.1801.0466.2966.191905857-0.339-1.6850.1250.830.0050.0243.0543.001945756-0.405-1.0550.8352.050043.0056.801955655-0.440-5.5652.0552.980.0120.0456.8057.503941173290.0070.3952.5452.360.0070.3950.32500290.4152.0551.710.0290.4157.5056.84从对照表可以看出，虽然管网中只出现了一个逆止阀，但整个管网的水力状态受到很大影响。节点1是供水节点，离逆止阀较远，所以没有受到大的影响；节点55和56的节点水头普通升高，增幅达9%；节点57和58的水头则大幅度下降，下降幅度达15%；节点119和329的水头也有所下降，但下降幅度相对较小，是因为这两个节点离