图论法用于供水管网水力计算的研究
- 格式:doc
- 大小:36.00 KB
- 文档页数:17
图论法用于供水管网水力计算的研究摘要:图论理论是网络分析的主要工具,现用于管网的水力平衡计算,既充分发挥了图论理论的优势,使计算变得简便、迅捷,又可将管网附件加入计算,使结果更准确、更符合实际。
文中采用峰阵输入管网结构,使输入数据的工作量大大减少,易于编制程序,计算大型的复杂管网。
关键词:供水管网水力计算图论法前言供水管网的水力平衡计算是供水系统规划设计、经济评价和运行管理的基础。
水力平衡计算的目的就是在确定管径的情况下求出满足连续方程和能量方程的各节点压力水头和各管段流量。
目前常用的水力平衡计算方法有哈代-克罗斯法(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节点的用水量。
管段视为网络图中的对应边,管段的直径、管长、管道流量、摩损系数等作为管段对应边的权。
至此,与管网同构的网络图生成了。
图中箭头表示各条边的方向,即管段中水流方向。
网络图中节点与边的关联函数可以用完全关联矩阵I4×5表示如式(1)所示。
第3章 给水排水管网水力学基础 (2h)3.1 给水管网水流特征流态分析:<2000 层流雷诺数νVD=Re =2000~4000 过渡流水力光滑区eD80~4000 h f ∝V 1.75 >4000 紊流 过渡区85.0)2(4160~80eDe D hf ∝V 1.75~2阻力平方区 85.0)2(4160eD> h f ∝V 2紊流过渡区=过渡粗糙区 阻力平方区=紊流粗糙区恒定流与非恒定流:水力因素(水流参数)随时间变化 均匀流与非均匀流: 水力因素(水流参数)随空间变化 压力流与重力流:水流的水头:单位重量流体具有的机械能h / H (位置水头 位能Z)(压力水头 压能P/γ) (流速水头 动能V 2/2g)水头损失:流体克服流动阻力所消耗的机械能 (沿程阻力)(局部阻力)3.2 管渠水头损失计算沿程水头损失(frictional head loss):谢才(Chezy)公式 l RC v h f 22= (通用,R 水力半径=断面/湿周,C 谢才系数)达西-韦伯(Darcy-Weisbach)公式 gv D l h f 22λ= (适用于圆管满流,λ沿程阻力系数, )28Cg=λC 和λ的计算 ①科尔勃洛克-怀特公式:)Re53.38.14lg(7.17CR e C +-= )Re 51.27.3lg(21λλ+-=D e 简化 )Re 462.48.14lg(7.17875.0+-=R e C )Re462.47.3lg(21875.0+-=D e λ②海曾-威廉(Hazen-Williams)公式:148.0852.113.016.13qC gD W=λlDC q h Wf 87.4852.1852.167.10=(v=0.9m/s 时)注:81.000)(vvC C W W = (v 0=0.9m/s ) ③曼宁(Manning)公式:6/11R nC =(n 曼宁粗糙系数) lR v n h f 3/422=l D q n 333.52229.10=3/22/13/22/12/123/41)()(R i nn R lh ln R h v f f === ④巴普洛夫斯基公式:yR nC 1=(n 曼宁粗糙系数) 式中)10.0(75.013.05.2---=n R n y局部水头损失(local head loss ):gv h m 22ζ= (ζ局部阻力系数)水头损失公式指数形式:n f n m nf q s l aq l Dkq h === (a 比阻,s f 磨阻系数)n m m q s D g q g v h ===422282πζζ (s m 局部磨阻系数) 总:n m f m f g q s s h h h )(+=+= (s g 管道磨阻系数)3.3 非满流管渠水力计算满流:曼宁公式6/11R n C =,谢才公式l RC v h f 22=lR v n h f 3/422=,满流时l Dq n 333.52229.10= 2/13/23/22/12/123/41)()(I R nn R lh ln R h v f f === 2/13/2I R nA Av q == 非满流:充满度 y/D ,管中心到水面线夹角θ2/)2cos 1(/θ-=D y)21(cos 21Dy-=-θ)sin (82θθ-=D A)sin 1(4θθ-=D R则θθsin 10-=R R ,R 为非满流时水力半径,R 0为漫流时水力半径; πθθ2sin 0-=A A ,A 为非满流时过水断面,A 0为满流时过水断面; 323200)sin 1()(θθ-==R R v v ,v 为非满流时流速,v 0为满流时流速; 3235320002)sin ()(πθθθ-==R R A A q q ,q 为非满流时流量,q 0为满流时流量; (y/D=0.94时,q/q 0=1.08最大;y/D=0.81时,v/v 0=1.14最大)l D q n h f 333.520229.10=31620229.10D q n I l h f == nD I q 29.1038210= 2/32/13/83/516.20)sin (⎥⎦⎤⎢⎣⎡∙-=nq I D θθθ,23/83/53/2)sin (16.20⎥⎦⎤⎢⎣⎡-∙=D nq I θθθ例题:某污水管道设计流量q=100L/s ,采用水力坡度I=0.007,拟采用D=400mm 钢筋混凝土管,粗糙系数n=0.014,求充满度y/D 和流速v 。
第21卷 第3期2023年5月中国水利水电科学研究院学报(中英文)JournalofChinaInstituteofWaterResourcesandHydropowerResearchVol.21 No.3May,2023收稿日期:2022-10-25;网络首发时间:2023-05-09网络首发地址:https:??kns.cnki.net?kcms?detail?10.1788.TV.20230509.0947.001.html基金项目:国家自然科学基金项目(52079143);国家重点研发计划项目(2018YFC0407705,2022YFC3204601)作者简介:郭聪(1998-),硕士生,主要从事水资源配置研究。
E-mail:ln7xsz@163.com通讯作者:游进军(1977-),博士,正高级工程师,主要从事水资源配置研究。
E-mail:youjj@iwhr.com文章编号:2097-096X(2023)-03-0295-12图论在水资源系统建模与配置计算中的应用郭 聪,游进军,林鹏飞(中国水利水电科学研究院流域水循环模拟与调控国家重点实验室,北京 100038)摘要:为提高复杂水资源系统的供排水关系识别与模型构建计算能力,本文基于系统网络图探讨了图论在模拟水资源系统复杂水量传输过程、确定网络遍历顺序与系统配水计算中的应用。
利用图论中的拓扑排序算法与深度优先搜索算法对水资源系统网络关系进行不同目标的分析识别。
综合两种算法,进行一次供水、供水还原与二次供水的多次网络分水过程设计和迭代计算,结合最大最小公平算法等水量分配原则,设计基于网络系统的复杂水量分配计算过程,将用户需水、工程可供水量和供水网络关系结合进行供需自适应配置,解决多水源多用户对应关系下的水源供给不均衡问题。
以三亚市地表水资源配置开展实例研究,结果表明通过图论算法在系统拓扑关系遍历和配水计算过程中的应用,系统枯水年城镇缺水率从3.42%降至1.63%,供需均衡度从0.38提高至0.71,提高了配置的计算效率与合理性。
收稿日期:2000Ο07Ο14作者简介:郑苏娟(1965—),女,江苏南京人,讲师,硕士,主要从事图论研究.基于图论的城市供水管网抢修决策信息系统郑苏娟1,徐筱麟2,丁莲珍1(1.河海大学数理系,江苏南京 210098;2.解放军理工大学2系,江苏南京 210016)摘要:利用图论分析方法,建立了城市供水管网抢修决策数学模型,给出了该模型的系统软件.应用表明,该系统能迅速提高维修人员查找所要关闭阀门的速度.关键词:供水管网;抢修决策;数学模型中图分类号:T U991.61 文献标识码:A 文章编号:1000Ο1980(2001)05Ο0092Ο04城市自来水管网的管理和维护,是一项既复杂又重要的工作,它和城市人民的生活息息相关.在日常维修特别是在应急抢修中,为了使对用户的影响降低到最低程度,往往需要由具有多年工作经验的专门技术人员提出决策抢修方案.随着城市规模的不断扩大,供水线路的不断延伸,对决策人员的要求就会越来越高,抢修决策将会越来越复杂.因此,应用计算机实现工作中急需的管网抢修智能决策和对管网设施实行规范化管理,是城市自来水管网设施管理的必然趋势.1 自来水管网的基本组成和抢修决策的基本要求如图1所示,城市自来水管网的主要组成要素有:主管———自来水管网中主干管线;次管———自来水管网中通往用户的小口径管线;水源Q ———水厂及增压站(储水柜);用户水表T ———主管到次管压力转换的装置;阀门f ———自来水管网中,只需关闭少量阀门即可切断供水的最小区域.由于抢修决策主要针对主管和次管,所以图上没有区分主管和次管.自来水管网计算机抢修决策的基本要求是:(a )能在多水源情况下自动辨认最小封闭区的进水阀和出水阀;(b )在阀门失灵情况下具有延伸封闭区的功能;(c )能确定封闭区阀门关闭之后受影响的用户水表,输出用户信息;(d )可输出封闭区内的阀门单位图和所在范围内的管网地形图.由于城市中供水管网相互交错,液体具有双向可流动性,要采用计算机进行智能决策,必须建立合理的数学模型.2 数学模型由于供水管网的构造具有明显的点线特征,用图论的分析方法是最合理的.然而,诸如图1(a )所示的供水管网结构要用图论进行分析还比较困难,必须对其进行变形简化,以建立便于图论分析的数学模型.其基本思想是将封闭区收缩为一个点,将阀门拉长成一条线,并将调压站信息隐藏到封闭区,如图1(b )所示.3 供水管网抢修决策方案产生的原理在自来水管网中,只需关闭少量阀门即可切断供水源的最小区域.该区域称为封闭区.其中水源既可以是水厂,又可以是储水柜.同一封闭区中任意一点出现故障,故障区切断水源的抢修方案是相同的,每一封闭区对应一个抢修方案.因此,管网中有多少封闭区,就对应有多少个抢修方案;封闭区的确定取决于自来水管网的结构,管网结构一旦形成,封闭区的划分就不会改变.由于自来水管网有其相对稳定性,封闭区划分生成后可供长期使用.要使某一封闭区停水,所需关闭的阀门必定在围着这一封闭区内的几个阀门之中.这样,要关闭的只是第29卷第5期2001年9月河海大学学报JOURNA L OF H OH AI UNI VERSITY V ol.29N o.5Sep.2001其中与水源相连的阀门.一个封闭区被切断水源后,有时还会影响其他封闭区.所有影响的封闭区至少有一个与事故点相邻,它们是以故障点所在封闭区为“水源”的那些封闭区,可以根据这一特点向下检索直到找出所有受影响的封闭区.以图1(a )为例,该图G 的有序三元组为[B (G ),F (G ),ΨG ].其中:B (G )={b 1,b 2,b 3,b 4,b 5,b 6,b 7};F (G )={f 1,f 2,f 3,f 4,f 5,f 6,f 7};ΨG (f 1)=b 1b 2,ΨG (f 2)=b 2b 5,ΨG (f 3)=b 2b 3,ΨG (f 4)=b 3b 4,ΨG (f 5)=b 5b 6,ΨG (f 6)=b 3b 5,ΨG (f 7)=b 3b 7.图1 基于图论的供水管网数学模型Fig.1 Mathem atical model of w ater piping netw ork从中首先得到管网封闭区的关联矩阵M (G )和相邻矩阵A (G ). 得到关联矩阵M (G )和相邻矩阵A (G )后,将从每个水源所在的封闭区出发,根据相邻矩阵A (G )给出的关系生成相互供水关系矩阵,此举意味着有N 个水源就要搜索N 次.以图1为例,先从Q 1所在的封闭区b 1出发进行搜索,然后从Q 2所在的封闭区b 6出发搜索,则获得单个管网封闭区的供水关系矩阵(行为供水封闭区,列为被供水封闭区):将各个管网封闭区的供水关系矩阵相加,就可以得到整个管网的封闭区相互供水关系矩阵39第29卷第5期郑苏娟,等 基于图论的城市供水管网抢修决策信息系统例1 若封闭区b 5发生故障,b 5给b 2,b 3,b 6供水,但检索b 2,b 3,b 6的供水封闭区时,都存在其他封闭区供水,所以该故障只影响b 5封闭区的用户,即T 3,T 5.从M (G )中可知,与b 5相连的f 2,f 5,f 6均为进水阀.例2 若封闭区b 3发生故障,b 3给b 2,b 4,b 5,b 7供水,但检索b 2,b 4,b 5,b 7的供水封闭区时,b 4,b 7仅由封闭区b 3供水,所以该故障不仅影响b 3封闭区的用户,而且还影响封闭区b 4,b 7,即T 2,T 4,T 6,T 7.从M (G )中可知,与b 3相连的f 4,f 7为出水阀,f 3,f 6为进水阀.如果在要求关闭的阀门中,由于走丝、损坏及其他原因,该阀门无法关闭,就需进行二次决策,即通过将相邻的两个封闭区合并实现封闭区的延伸.在决策过程中,用户在第一次决策的基础上对无法正常关闭的阀门进行点取,就可完成二次决策.4 抢修决策信息系统总体结构根据实际工作对自来水管网抢修决策、管网维护及图档信息管理的基本需要,建立了城市供水管网抢修决策信息系统,如图2所示.图2 系统软件模块结构Fig.2 Module structure of system softw are5 结束语该系统已在南京市自来水管网的爆管分析和抢修决策中得到应用.由于采用本系统,维修管养部门计算受影响的用户范围,列出要关闭的阀门的全部信息,从接到电话到赶往现场止水过程仅1.5h.如果按以往的人工决策抢修过程,像这类大口径爆管的维修工作至少需要4h.显然,这大大提高了抢修速度.该图论分析方法也在城市煤气抢修决策中得到应用.如煤气中压管断裂抢修、中压管线改道等实际施工工程,由该系统产生的停气降压方案与人工产生的最佳方案完全一致,并准确给出了由于停气所影响的煤气用户.参考文献:[1]邦迪J A 主编.图论及其应用[M].吴望名等译.北京:科学出版社,1987.1~60.49河 海 大 学 学 报2001年9月Decision Making I nformation System for U rgent R epair of U rban W aterSupply Piping N etw ork B ased on G raph TheoryZHENG Su 2juan 1,XU Xiao 2ling 2,DING Lian 2zhen 1(1.Dept.o f Math.and Phys.,Hohai Univ.,Nanjing 210098,China ;2.Dept.o f Two ,the P LA Engineering Univ ,Nanjing 210016,China )Abstract :Based on the graph theory ,a mathematical m odel for urgent repair decision making of urban water supply pipng system is established ,and the system s oftware is discribed.Application shows that the system may help find quickly the valves to be closed.K ey w ords :water supply piping netw ork ;urgent repair decision making ;mathematical m odel《河海大学学报(自然科学版)》征订启事《河海大学学报(自然科学版)》是以水资源开发、利用与保护为重点的综合性学术期刊,主要刊登本校在水资源、水文、地质、测量、水利工程、水电工程、水运工程、海洋及海岸工程、水工结构、工程力学、水力学及河流动力学、岩土工程、计算机科学、电力工程、电子技术及自动化工程、工业与民用建筑、管理工程、水利经济、环境工程、机械工程等学科方面的科研成果学术论文、学术讨论、研究动态等学术性文章,可供上述有关专业的科技工作者及大专院校师生阅读和参考.本刊创办于1957年,是我国中文核心期刊,在国内工程技术界和学术界有较大影响.刊载的文章中,有不少国家科技攻关(重点)项目和各种科学基金资助项目的研究成果,部分达到了国内领先和国际先进水平,为我国水利、水电、水运工程及其他有关工程建设的规划、设计、施工和管理提供了科学理论、方法和具体建议,发挥了较大的社会效益和经济效益,深受工程界和科技界赞许,并获得全国优秀高校自然科学学报二等奖,以及全国水利系统优秀期刊、江苏省优秀期刊称号.本刊每逢单月出版,国内外公开发行,邮发代号:28Ο63,每期定价:8.00元,全年订费48.00元.欢迎广大读者在全国各地邮局订阅或直接与编辑部联系.联系地址:南京市西康路1号《河海大学学报》编辑部,邮政编码:210098.59第29卷第5期郑苏娟,等 基于图论的城市供水管网抢修决策信息系统。
图论法研究多水源管网供水分区阎立华;刘丽英;马婷婷【期刊名称】《沈阳建筑大学学报(自然科学版)》【年(卷),期】2007(023)001【摘要】目的为了从整体上把握供水管网各水源供水区域的范围,对多水源管网进行供水分区,使得各水源的供水区域可视化.方法从供水管网是一种图的观点出发,基于管网水力计算的平差结果,根据图论原理,生成供水管网流向图的邻接矩阵;在邻接矩阵的基础上搜索水源点的供水情况,形成判断出水源的供水分界点矩阵,进而确定了各水源的供水区域;利用Visual Basic和MapInfo联合编程技术将该方法进行程序编制.结果通过调用复杂多水源管网进行算例分析,实现了多水源管网供水分区的可视化显示.可以从管网总体上获得各个供水区域的范围、大小和分区的具体位置.结论对运行管理人员进行管网水量调配和设计人员进行管网改造具有指导作用.有效解决了多水源管网供水分区的问题,而且该方法具有算法简单可靠、程序编制简单、适用范围广、集成化程度高等特点.【总页数】4页(P113-116)【作者】阎立华;刘丽英;马婷婷【作者单位】沈阳建筑大学市政与环境工程学院,辽宁,沈阳,110168;沈阳建筑大学市政与环境工程学院,辽宁,沈阳,110168;沈阳建筑大学市政与环境工程学院,辽宁,沈阳,110168【正文语种】中文【中图分类】TU821.3【相关文献】1.城市供水管网分区供水技术研究 [J], 林英姿;赵学敏;赵泽明2.爆管事故对多水源管网供水分区的影响 [J], 戴成林;吕谋;董深;李璞3.城市供水管网分区供水技术研究 [J], 郭维剑4.城市供水管网分区供水技术研究 [J], 马彪5.多水源管网供水分区的可视化研究 [J], 刘丽英;阎立华;宋占龙因版权原因,仅展示原文概要,查看原文内容请购买。
图论法用于供水管网水力计算地研究河南省黄河河务局张丰周武汉水利电力大学魏永曜天津市大港油田集团公司设计院石继摘要图论理论是网络分析地主要工具,现用于管网地水力平衡计算,既充分发挥了图论理论地优势,使计算变得简便、迅捷,又可将管网附件加入计算,使结果更准确、更符合实际.文中采用峰阵输入管网结构,使输入数据地工作量大大减少,易于编制程序,计算大型地复杂管网.关键词供水管网,水力计算,图论法.前言供水管网地水力平衡计算是供水系统规划设计、经济评价和运行管理地基础.水力平衡计算地目地就是在确定管径地情况下求出满足连续方程和能量方程地各节点压力水头和各管段流量.目前常用地水力平衡计算方法有哈代-克罗斯法(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节点地用水量.管段视为网络图中地对应边,管段地直径、管长、管道流量、摩损系数等作为管段对应边地权.至此,与管网同构地网络图生成了.图中箭头表示各条边地方向,即管段中水流方向.网络图中节点与边地关联函数可以用完全关联矩阵I4×。
市政给水管网水力计算问题研究摘要:目前市场上出现的排水给水管材的规格和类别非常多,这给水力计算带来了很大的麻烦。
以往管理给水管网时基本属于经验式管理,存在科学性差。
随着测流点、测压点在市政给水管网中的设置,管网建模逐渐进入了实用化阶段。
通过介绍给水管网模型,介绍管网水力计算方程的研究问题。
关键词:水力计算;市政给排水;建模在市政给水管网的设计中,水力计算是管网设计的计算基础。
根据管网形状和管材不同,采用的参数或公式就不同。
随着管材市场的不断发展,目前市场上出现的给水排水水管的规格和类别越来越多,这给水力计算带来了很大的麻烦。
虽然有设计给水排水管道的相关设计手册中规定了针对各种管材的水力计算公式,但是还是不能够满足日益增多的管材规格,另外在查算时也非常不方便。
在目前的管网设计中,通常通过建立微观管网模型来获取动态水力信息,进而进行水力计算,但是由于技术限制,这种方法在使用过程中受到限制。
因此探究市政给水管网水力计算研究问题具有非常重要的意义。
供水管网模型就目前研究的供水管网模型类型来看,管网模型的类型包括了宏观和微观两种管网模型。
建立管网宏观模型时运用回归计算的方法,运用此方法的前提是基于大量的运行数据以及模型服从管网流量“比例负荷”。
通过这种计算方法,能够建立控制点压力分布以及在管网中各个水厂的供水压力的函数关系。
由于建立宏观模型是建立在统计的回归模型上,它的计算速度非常快,所以这种建模方法通常用在给水系统模块调度中,而在扩建、改建或者新建给水管网模块中并不适合。
根据实际的管网情况,管网中的管段、水泵以及阀门等全部的元素,不通过简化处理而建立的模型即为微观模型。
通过解环方程、解节点方程以及解管段方程能够将管网中节点以及管段的信息。
通过建立微观模型能够将给水管网中水力的全部运行状态准确表达出来,其重点表达的是水力实时状态和信息。
由于受到技术限制,一些管道的基础参数和拓扑关系的完整性很难获取,尤其是受到设备的限制,不能准确地将管网节点流量的动态数据准确获取。
图论法用于供水管网水力计算的研究图论法用于供水管网水力计算的研究摘要:图论理论是网络分析的主要工具,现用于管网的水力平衡计算,既充分发挥了图论理论的优势,使计算变得简便、迅捷,又可将管网附件加入计算,使结果更准确、更符合实际。
文中采用峰阵输入管网结构,使输入数据的工作量大大减少,易于编制程序,计算大型的复杂管网。
关键词:供水管网水力计算图论法前言供水管网的水力平衡计算是供水系统规划设计、经济评价和运行管理的基础。
水力平衡计算的目的就是在确定管径的情况下求出满足连续方程和能量方程的各节点压力水头和各管段流量。
目前常用的水力平衡计算方法有哈代-克罗斯法(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节点的用水量。
管段视为网络图中的对应边,管段的直径、管长、管道流量、摩损系数等作为管段对应边的权。
至此,与管网同构的网络图生成了。
图中箭头表示各条边的方向,即管段中水流方向。
网络图中节点与边的关联函数可以用完全关联矩阵I4×5表示如式(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。
网络图的生成树(全涉及树)可以有很多种,在计算时可以任选一种。
在本例中,选1、2、4这3条边为图的生成树,则补树(余树)的各边(弦)为3、5.各弦将与枝构成基本回路,一个基本回路中有且仅有1条弦。
用基本回路矩阵Bf表示则如式(2)所示。
枝1245 Bf= [ -1 1 0 1 0 ] (2)1-1-11式中每一行表示一个基本回路(环)。
环的方向以该环对应弦的方向为准。
“-1”表示管段中的流向与环中弦的方向相反,“1”表示相同,“0”表示该管段不在此环内。
Bf可用矩阵B和单位阵U表示为式(3)。
Bf=[B|U],其中B=[-11](3)1-1-1环阵与管段摩损列向量hf构成环方程如式(4)所示。
摩损向量的元素顺序与Bf中每行元素所对应的管段顺序相同。
Bf×hf=0。
其中hf=(h1,h2,h4,h3,h5)T(4)图论理论中,连续方程用割方程代替。
每个割方程只含一根枝,并和相关的弦构成割集,将图2分割成互不连通的脱离体。
这样,图中就有3个割集。
割集和割集阵Af如式(5)所示:割集K:割阵: Af= [ 枝1 2 4 弦3 5 ] (5)K1=(e1,e3,e5) 1 0 0 1 -1K2=(e2,e3,e5) 0 1 0 -1 1K3=(e4,e5) 0 0 1 0 1割阵Af中,每一行表示一个割集。
图中有3根枝,所以就有3个割集。
割阵中,“+1”表示该管段在此割集内,且管段流向与此割集内的枝中的流向相同,“-1”表示流向相反,“0”表示该管段不在此割集内。
式(5)的割阵Af和割集K一一对应。
割阵Af可用一个矩阵A 和一个单位阵U表示为:Af=[U|A],其中A=[1-1]-111割阵与流量列向量可构成割方程。
根据图论理论,割阵的行向量与环阵的行向量正交,这种关系可用式(6)表示。
[B|U]·[U|A]T=0或者[U|A]·[B|U]T=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个未知量[1]。
因此,必须借助P个管段摩损方程式。
管段摩损方程式线性化后的通式如(7)和(8)所示。
系数R称为阻尼系数,Y称为传导系数。
R和Y的具体形式与所选用的摩损公式有关,是D、C、L的函数。
摩损公式线性化后,R还是q的函数,Y还是h的函数。
不过,在求解过程中,总是把R和Y当作已知量来对待。
阻尼式:h=R×q(7)传导式:q=Y×h(8)式中R和Y是阻尼系数和传导系数矩阵。
如果摩损公式采用Hazen-William公式,则有:h=R×q=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为脱离体所含节点流量之和。
方环程: Bf×h=0,即 [B U]×[ h ] =0 (12)h′割方程: Af×q=Q,即 [UA]×[ q ] =Q (13)q′传导式:[ q ] = [ Y 0 ]×[ h ] (14)q′ 0 Y′ h′求枝摩损式(以管段摩损为未知量):首先将传导式(14)代入割方程(13)得:[U A]×[Y]×[h]=Q (15)Y′h′由环方程(12)可得Bh+h′=0,即h′=-Bh,代入式(15)得:[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)求得。
式(11)中C1.852×D4.87/10.68×L对某一管段来说是个常数,可用W 表示。
则传导系数Y可以表示为:Y=W×|q|-0.852(19)在迭代计算时,第一次可以直接用W代替Y进行计算,求出h\,q 后计算Y,再求新的q值,如此反复计算,直至前后两次的q值符合给定的误差标准为止。
为了避免可能出现的数值摆动现象,在第三次迭代时,用前两次迭代结果的流量平均值作为初始流量值[2],即:q=q(1)+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=2ξq/π2gD2;传导系数为:Y=π2gD2/2ξ×q-1计算时只需将闸板式阀门的R或Y值加入,即可计算。
蝶阀的计算方法与闸板式阀门类似3.2 逆止阀逆止阀是管网中最常见的设备之一,是水流方向控制设备,只允许水流单向通过。