摘要：本文针对现行给水管网电算程序的一些缺陷,作了若干改进，使其输入原始数据简单方便，并能自动圈划回路，适用于不同结构或多水源管网设计与运行管理,能输出管段流量、结点水压、水塔高度或水泵工况及其轴功率等结果，比现行程序更方便、通用、有效。简例介绍了新程序的功能。
关键词：给水系统 水力计算 解算程序
modification of calculating program for water supply pipe network
abstract: not being satisfactory at function, the present calculating program for water supply pipe network is modified in this paper.the new program,which selects independent meshes automatically,is more convenient to input original data,thus it is suitable to designing different structures of pipe networks and to running management of more than one water source,and it is capable of outputing pipe quantity,node pressure,water tower height or water pump working state and so on.in one word,the new program is more convenient and more effective than many other ones nowaday.simple examples are used to explain its function.
key words: water supply; pipe network; calculation; program
1 引 言
在现行给水管 网计算程序中，解环方程法管网平差程序存在如下缺点：（1）原始数据输入工作复杂，不仅要输入管段参数，而且要输入环数及管段所属环号，工作量大，自动化水平低，且极易出错。（2）通用性及适应性差，环状管网的程序无法应用于树状或混合状管网，对二级水泵直接供水或通过水塔供水，以及多泵站或泵与水塔、水塔与水塔联合供水等各种情况，现行很多程序不能通用，对给水管 网设计计算和运行管理计算，程序也只能分别处理。针对这些缺点，本文运用网络图论知识，对解环方程管网平差程序进行了改进，较好地满足了实际要求。
2 原始数据输入与水泵特性曲线拟合
将给水系统进行拓扑变换，使之成为强连通有向赋权图。对管网中各节点从1开始编号，每个节点流量视为一虚拟分支流量，虚拟分支阻力为r＝0，其流量终点即水泵吸水点。输入的原始数据（由数据库系统进行管理）有：各分支起、止节点编号，管段的管径、管长，分支中水泵型号，分支中恒压头（水塔和管段节点高差），以及分支流量。另外由一有格式数据文件输入分支、含泵站分支数、水泵特性曲线上五对流量与压力数据。
3 管道水头损失及节点水压计算
管道水头损失h按式（1）计算。
h=rq2 (1)
式中，q——管段流量，m3 /s；d——管径，m。r——摩阻系数，可按式（2）计算［1，2］ 。
r＝0．001629 l／d5．29 (2)
式中，l——管长，m；d——管径，m 。［3］ 。计算式为。
hbi =hai -hi +pfi +pvi (3)
式中hai 和hbi 分别为分支i中起、止节点的相对压力值；hi 为分支中的阻力；pfi 和pvi 分别为分支中水泵压头及恒压头。
4 阻力排序及独立网孔圈划
按图论知识，闭合链称为网孔或回路。一个包含网络所有节点但不形成回路相互连通的图称为树，树中分支称树支，其余分支称为弦。若以nb 、nj 、nm 分别表示网络的分支数、节点数、独立网孔数，则有nm =nb -nj +1。理论上，网孔圈划可由点边关联矩阵推演出回路矩阵来实现。但实际运算中此法较繁，故应用上多采用生成树原理进行圈划。先挑选出一棵最小生成树，再选出独立网孔［4］ 。q ＋nf ≤nm ，其中，nq ，nf 分别为管网中固定流量分支数及水泵分支数。输入的原始数据应符合该式。
5 迭代计算
采用hardy-cross方法，将回路压力平衡方程式按taylor级数展开，略去二阶以上的高阶微分项。利用方程式中一个已知根（即初拟流量）的近似值，逐次迭代计算，并校正其流量值。当迭代校正误差达到规定允许精度时，就求出了近似的真实值。具体步骤如下：
式中qj ——第j条分支流量； aij ——基本关联矩阵元素； nj ——网络的节点数；
nb ——网络的分支数。
式中，bij ——基本回路矩阵元素，aij 、bij 由程序根据输入的管网分支数自动生成。(0) i ，本文均取零。
(4) 对管网回路中各分支的初拟流量进行校正，如式
q(k+1) j =q(k) j +bij δq(k+1) j (7)
(5) 对所有独立回路，判断预期收敛精度，如式
(6) 当第（5）步不满足时，转第（3）步，直至满足式（8），得到管网计算结果。
6 程序框图及应用简例
本文编写的fortran程序包括主程序和4个子程序。其中子程序cvft求水泵特性能曲线，sort对分支阻力排序，mesh圈划独立网孔，itr进行迭代计算。见图1所示。
图1 改进后程序的框图
对所编程序采用图2所示简例进行调试和计算。已知管径d1 ～d8 ，各管长l，系数c，流量q9 ，q10 ，q11 ，q12 ，q13 ，q14 。
图2 计算简例采用的管网
考虑以下各种实际情况进行分别计算。v1 已知，可解算出分支流量及节点压力。v 。设q1 或q14 已知，计算得分支（1）中应设水塔的高度，分支1～8中的流量，节点1～7中压力。1 ，q14 ，pf1 ，pf14 ，管网总有效功率。1 或q14 已知，可算得1～7分支流量，1～7节点压力。还可计算出分支（1）中应设水塔的高度或设直接供水水泵的有效功率。1 ，pf1 ，功率，迭代次数。
7 结论
改进的给水管 网水力计算程序，输入原始数据简便，能自动圈划回路，适用于不同结构或多水源的管网设计或运行管理，能输出管段流量、节点水压、水塔高度或水泵工况及功率等结果。可见比现行很多程序更方便、更通用、更有效。
参考文献
［1］ 姚雨霖，周康伦，任周宇等. 城市给水排水. 北京：中国建筑工业出版社，1982. 109-119.
