中崛为您服务：戴先生 13739056275
图论法用于供水管网水力计算
摘要：图论理论是网络分析的主要工具，现用于管网的水力平衡计算 ，既充 分发挥了图论理论的优势，使计算变得简便、迅捷，又可将管网附件加入计算， 使结果更准确、更符合实际。文中采用峰阵输入管网结构，使输入数据的工作量 大大减少，易于编制程序，计算大型的复杂管网。
其中，a 表示管段中过水断面的高度，d 表示管段直径，a/d 表示阀门开关。 当开度为0时，阀门完全关闭，没有流量通过；当开度为1时，阀门完全打开，对 水流不产生影响。
将阀门水头损失公式用流量表示为：hf=ξ×v2/2g=ξ×2q2/π2gD2 则阻尼系数 R 为：R=2ξq/π2gD2;传导系数为：Y=π2gD2/2ξ×q-1 计算时只需将闸板式阀门的 R 或 Y 值加入，即可计算。 蝶阀的计算方法与闸板式阀门类似 3.2逆止阀 逆止阀是管网中最常见的设备之一，是水流方向控制设备，只允许水 流单向通过。
中崛为您服务：戴先生 13739056275
图2
图1中有一水库 A，三个给水点 B、C、D，Q1表示水库节点供水量，Q2\,Q3\,Q4 分别表示 B、C、D 节点的用水量。管段视为网络图中的对应边，管段的直径、管 长、管道流量、摩损系数等作为管段对应边的权。至此，与管网同构的网络图生 成了。图中箭头表示各条边的方向，即管段中水流方向。
中崛为您服务：戴先生 13739056275
割集 K：
割阵： Af= [ 枝1 2 4 弦3 5 ]
(5)
K1=(e1,e3,e5)
1 0 0 1 -1
K2=(e2,e3,e5)
0 1 0 -1 1
K3=(e4,e5)
0 010 1
割阵 Af 中，每一行表示一个割集。图中有3根枝，所以就有3个割集。割阵中，“+1” 表示该管段在此割集内，且管段流向与此割集内的枝中的流向相同，“-1”表示 流向相反，“0”表示该管段不在此割集内。式(5)的割阵 Af 和割集 K 一一对应。 割阵 Af 可用一个矩阵 A 和一个单位阵 U 表示为：
为了避免可能出现的数值摆动现象，在第三次迭代时，用前两次迭代结果的 流量平均值作为初始流量值［２］，即：
q=q（１）+q(2)2
(20)
中崛为您服务：戴先生 13739056275
中崛为您服务：戴先生 13739056275
求得 q(3)，……，这样收敛速度加快。 3管网附件
实际管网中，有许多控制、安全、量测设施，如加压泵、控制阀、逆止阀、 减压阀等附件，对管网运行产生重要影响。传统计算方法都未涉及到管网附件问 题，不仅使计算准确性受损，而且其计算程序无法用于日常管理工作。
图论法处理管网附件时，将附件所在管段视为特殊管段，这些管段的摩损式 要根据其附件的水力学特征计算摩损值，再加入到管网中进行水力平衡计算。本 文给出几种较常见管网附件的处理方法。对于其它附件，具体问题具体处理，在 此就不一一详述了。 3.1普通阀门 闸板式阀门是用得最多的一种阀门，在一 般的水力计算过程中， 闸板式阀门的水头损失计算一般引用公式 hf=ξ×v2/2g,ξ值见文献[3]。
将供水管网中的管段概化成一条线段(即图中的边)，将有附件的管段看成图 中的特殊管段，边与边由节点相连。这样，一个供水系统的管网图就转化为图论 中的网络图。而且管道中的水流是有方向的，所以管网图是有向图。
根据以上所述原则，可将图1所示管网系统，转化为图2所示的网络图。
图1
中崛为您服务：戴先生 13739056275
1 -1 Af=阵与流量列向量可构成割方程。
根据图论理论，割阵的行向量与环阵的行向量正交，这种关系可用式(6)表 示。
[B|U]·[U|A]T=0或者[U|A]·[B|U]T=0
(6)
所以有 B=-AT 或者 A=-BT。这样，环阵可以由割阵求出，反之亦然。
［U A］×［ Y 0 ］ × ［
h］ =Q (16)
0 Y′
-Bh
即
h×［Y-AY′B］=Q
(17)
根据正交定理得：
h×［Y+AY′AT］=Q
(18)
这就是图论法的求枝摩损式计算公式。h 即为枝管段的摩损向量。解得枝摩 损值 h 后，其余变量可由相应的公式求出。由环方程可得 h′=-B×h，即可求出 弦摩损向量 h′,q、q′向量可以由式(14)求得。
1 0
01]
(2)
中崛为您服务：戴先生 13739056275
中崛为您服务：戴先生 13739056275
式中每一行表示一个基本回路(环)。环的方向以该环对应弦的方向为准。“-1” 表示管段中的流向与环中弦的方向相反，“1”表示相同，“0”表示该管段不在此 环内。Bf 可用矩阵 B 和单位阵 U 表示为式(3)。
任何管道的水力计算都可以用管段流量 q\,水头损失 h\,管径 D\,管长 L 和 管壁条件 C 等5个因素来描述。一般 D、L 和 C 为已知条件，只有 q 和 h 未知。因 此，求解一个管网的水力平衡问题，可从两方面考虑：一是利用 q 和 h 的关系， 消去 h，以 q 为未知量计算，求出 q 后，反求 h；二是首先消去 q，以 h 为未知 量计算；解出 h 之后，再反求。图论法也可从这两方面入手，即求弦流量式和求 枝摩损式。前者只适用于环状网，而后者则适用于所有类型的管网，所以本文着 重介绍后者。
网络图的生成树(全涉及树)可以有很多种，在计算时可以任选一种。在本例 中，选1、2、4这3条边为图的生成树，则补树(余树)的各边(弦)为3、5.各弦将 与枝构成基本回路，一个基本回路中有且仅有1条弦。用基本回路矩阵 Bf 表示则 如式(2)所示。
枝1 Bf= [
2 4 弦3 5
-1 1 0 1 -1 -1
式(11)中 C1.852×D4.87/10.68×L 对某一管段来说是个常数，可用 W 表示。则传导 系数 Y 可以表示为：
Y=W×｜q｜-0.852
(19)
在迭代计算时，第一次可以直接用 W 代替 Y 进行计算，求出 h\,q 后计算 Y， 再求新的 q 值，如此反复计算，直至前后两次的 q 值符合给定的误差标准为止。
-1 1 0 Bf=［B｜U］,其中 B= [ 1 -1 -1 ] (3)
环阵与管段摩损列向量 hf 构成环方程如式(4)所示。摩损向量的元素顺序与 Bf 中每行元素所对应的管段顺序相同。
Bf×hf=0。其中 (4)
hf=(h1,h2,h4,h3,h5)T
图论理论中，连续方程用割方程代替。每个割方程只含一根枝，并和相关的 弦构成割集，将图2分割成互不连通的脱离体。这样，图中就有3个割集。割集和 割集阵 Af 如式(5)所示：
(9) (10) (11)
用 h 向量表示管段摩损：h 表示枝摩损，h′表示弦摩损；
用 q 向量表示管段流量：q 枝管段流量，q′表示弦管段流量。
割方程的右端项 Q 为脱离体所含节点流量之和。
方环程： Bf×h=0,即 [B U]×
［ h］
=0 (12)
中崛为您服务：戴先生 13739056275
中崛为您服务：戴先生 13739056275
设一管网有 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 当作已知量来对待。
关联矩阵通过选主元初等行变换即可得到割阵：先选关联阵第一行中一非零
中崛为您服务：戴先生 13739056275
中崛为您服务：戴先生 13739056275
枝元素为主元，并使其为+1，消去其它各行中此主元；再选第二行、第三行、… 的主元，最后即得割阵 Af。因此，可以由关联矩阵导出割阵和环阵。2图论法模 型
网络图中节点与边的关联函数可以用完全关联矩阵 I4×5表示如式(1)所示。
顶点
边的编号
(1)
式中：Iij= { 1，表示 j 管段与 i 节点相连，且管内水流流离该节点； 0，表示此管段不与该节点关联；
-1，表示 j 管段与 i 节点相连，且管内水流流入该节点。
完全关联矩阵与管段流量列向量 q 以及节点流量列向量 Q 可组成管网节点方 程(即连续方程)Iij×q+Q=0,q=(q1,q2,q3,q4,q5)T,Q=(Q1,Q2,Q3,Q4)T。
关键词：供水管网 水力计算 图论法
前言 供水管网的水力平衡计算是供水系统规划设计、经济评价和运行管理的基
础。水力平衡计算的目的就是在确定管径的情况下求出满足连续方程和能量方程 的各节点压力水头和各管段流量。目前常用的水力平衡计算方法有哈代-克罗斯 法(Hardy-Cross)，牛顿-莱福逊法(New ton-Raphson)，线性理论法 (Linear-Theory)，有限元法(Finite Element)等等。所有这些方法各有所长， 适用范围各不相同，有的还需人工假设管段流量，使输入数据工作量增大，且未 考虑管网附件的影响。本文介绍的图论法将复杂的管网处理为相应的“网络图”， 并建立相应的数学模型，用峰阵输入原始数据来描述管网结构，输入的数据量最 少，不易出错，易于计算大型的复杂管网。其计算过程可同时考虑管网附件，如 控制阀、加压泵、逆止阀、减压阀等，使计算结果更符合实际。 1图论原理
阻尼 h=R×
式：
q
(7)
传导 q=Y
式：
×h (8)