高金良，趙洪賓，陳兵
（哈爾濱工業(yè)大學(xué)給排水系統(tǒng)研究室）

　　摘要： 本文提出了一個管網(wǎng)簡化的算法，具有可以適應(yīng)各種管網(wǎng)拓?fù)浣Y(jié)構(gòu)、計算結(jié)果準(zhǔn)確、可任意確定管網(wǎng)簡化程度的特點。經(jīng)比較管網(wǎng)簡化前后的模型，簡化精度完全滿足管網(wǎng)計算的要求。該算法的實現(xiàn)，為供水系統(tǒng)實時優(yōu)化調(diào)度奠定了堅實基礎(chǔ)。
　　關(guān)鍵詞：復(fù)雜管網(wǎng)；簡化；線性模型；非線性模型

Calculation and Analysis of Complex Water Supply Network Simplification
Gao Jinliang, Zhao Hongbin, Chen Bing
(Workgroup of Water & Wastewater System, Harbin Institute of Technology)

　　Abstract: This paper describes a algorithm of network simplification. It has some advantages that it can adapt many kinds of topology network, get accurate result and determine the extent of network simplification. The precision satisfy the requirement of network calculation by comparing the model of network simplification before and after. We realize the algorithm, so build the stabilized foundation of real-time optimal dispatch in water supply system.
　　Keywords: complex network; simplification; linear model; non-linear model

一、管網(wǎng)簡化的目的、意義

　　我國城市配水管網(wǎng)是保證城市居民生活、生產(chǎn)用水的重要市政設(shè)施。管網(wǎng)輸送水到千家萬戶，隨著城市的發(fā)展，逐漸形成密如蛛網(wǎng)的配水管網(wǎng)。這造成了管網(wǎng)計算分析的極大困難。為此我們必須對極其復(fù)雜龐大的配水管網(wǎng)進(jìn)行簡化。管網(wǎng)簡化的結(jié)果主要有三個用途：加速管網(wǎng)平差、得出宏觀管網(wǎng)、實現(xiàn)管網(wǎng)優(yōu)化調(diào)度。
　　近年來，計算機(jī)技術(shù)迅猛發(fā)展，計算速度迅速提高，普通微機(jī)的時鐘頻率已達(dá)1GHz。同時管網(wǎng)平差計算在引入了序列分析、壓縮矩陣、有限元等技術(shù)后算法計算效率也大大提高。如對有4000個節(jié)點，5000條管段的特大城市復(fù)雜管網(wǎng)平差，20秒內(nèi)即可完成。目前可不需要事先進(jìn)行管網(wǎng)簡化即可實現(xiàn)復(fù)雜管網(wǎng)平差模擬計算。同時供水綜合信息系統(tǒng)、管網(wǎng)GIS技術(shù)日臻成熟^[1]，現(xiàn)已不難得出管網(wǎng)的宏觀干線拓?fù)溥B接圖形。
　　但對于管網(wǎng)的優(yōu)化調(diào)度這一重要課題，由于至今對大規(guī)?；旌想x散非線性模型求解十分困難，耗時驚人。優(yōu)化問題的數(shù)學(xué)復(fù)雜性要遠(yuǎn)比管網(wǎng)簡化更大，因此需要簡化后的管網(wǎng)模型，以使得優(yōu)化調(diào)度計算的時間能令人接受。目前，供水系統(tǒng)優(yōu)化技術(shù)應(yīng)用進(jìn)展緩慢的一個主要原因是缺少有效的管網(wǎng)簡化方法。故管網(wǎng)簡化對于供水系統(tǒng)計算分析仍是十分必要的。而且由于目前管網(wǎng)簡化的目標(biāo)單一，不須用于前二個用途，也使得管網(wǎng)簡化算法的考慮可擺脫很多限制，更易于成功實現(xiàn)。

二、管網(wǎng)簡化的現(xiàn)有方法

　　管網(wǎng)簡化的現(xiàn)有方法主要有：管線省略、合并、分解。
　　⒈管線省略認(rèn)為，通常情況下，管徑的大小與管線在供配水中的作用和水力條件的影響成正比。在管網(wǎng)計算允許的范圍內(nèi)，逐步去掉小管徑的管線，同時對省略管線附近的管線的水力參數(shù)進(jìn)行處理。但需要人為地判斷管線的對管網(wǎng)工作的影響程度（如干線間的小連通管），且不可避免產(chǎn)生計算誤差，尤其對于成環(huán)的小管徑的管線很難省略。
　?、苍诠芫€合并處理中，并聯(lián)、串聯(lián)管線的合并在理論上可以得到精確的計算結(jié)果，在大中城市的給水管網(wǎng)中，有30%的管線可以做此簡化處理。此外平行管線合并和管線省略一樣需要考慮對其四周管線的處理，有一定誤差。
　?、彻芫W(wǎng)分解是當(dāng)一個供水區(qū)域僅由少數(shù)幾條管線供水和其他管網(wǎng)相連時，可以將連接的管線斷開，將一個管網(wǎng)分解為兩個獨立的管網(wǎng)，對兩個管網(wǎng)分別單獨計算。這在連接管線的流向和流量易于確定的時候可以分解，但在大多數(shù)管網(wǎng)成環(huán)布置、多處連接的情況下難以實現(xiàn)。

三、管網(wǎng)簡化的理論

　　對于給水管網(wǎng)，節(jié)點和管段之間關(guān)聯(lián)可以由關(guān)聯(lián)矩陣Λ來表示。矩陣的每行對應(yīng)著每一個節(jié)點n=1,2,Λ,N；每列對應(yīng)著每一條管段j=1,2,_Λ,J。。在矩陣Λ中，每一列僅有二個非零元1、-1，分別對應(yīng)于該管段的起始和終止節(jié)點，其余元素為0；每一行有1、-1，分別對應(yīng)于流入和流出該節(jié)點的管段。則給水管網(wǎng)的水頭損失方程；連續(xù)性方程和能量方程可表示為：

　　q=q(△h)=g│△h│^0.54sign(△h)......................(1)

　　q^nod=Λq............................................(2)

　　△h=Λ^Th............................................(3)

　　q=Q(△h)............................................(4)

　　其中，q-----管段流量向量；
　　　　　△h-----管段水頭損失；
　　　　　g-----管段關(guān)聯(lián)系數(shù)，由管長、管徑、阻力系數(shù)確定；
　　　　　q^nod------節(jié)點流量向量；
　　　　　h-----節(jié)點水壓向量；
　　　　　Q(△h)-----向量函數(shù)，Q(△h)=(q₁(△h₁),Λ,q_j(△h_J))^T。

　　由⑵⑶⑷式可得，

　　ΛQ(Λ^Th)=q^nod.....................................(5)

⒈管網(wǎng)模型線形化
　　將⑸式線性化對△h求導(dǎo)可得，

　　Λ(dQ/dΔh)Λ^T×δh=δq^nod...........................(6)

　　式中，設(shè)A=Λ(dQ/dΔh)Λ^T，此矩陣為雅可比矩陣。

　　且有。則A可得，

　　 ........⑺

⒉化簡線形化模型

　　對⑺式可以多次重復(fù)應(yīng)用高斯約當(dāng)消去過程。如將矩陣中消去前r行，即從管網(wǎng)模型中簡化掉r個節(jié)點。得新的簡化后的模型可寫為，

　　 ........⑻

　　對節(jié)點流量q也作相應(yīng)處理。
⒊化簡后模型恢復(fù)為非線形模型
　　經(jīng)過高斯約當(dāng)消去過程，得到了化簡后的管網(wǎng)模型A^(r)。新的矩陣A^(r)包含了原有線形模型的全部特性，即形成了一個新的管段聯(lián)系的新管網(wǎng)模型。其管段關(guān)聯(lián)系數(shù)g為，

　　 ........⑼

　　對于新的管網(wǎng)模型，可假定新生成的管段的阻力系數(shù)C值為100，管長等于管徑，來生成實際具有物理意義的管段。這就得到了新的簡化后的非線性的管網(wǎng)模型。
　　該管網(wǎng)簡化過程的物理解釋是，通過高斯約當(dāng)消去刪除一部分原有的管段，合成新的管段中保存了原有管段間的關(guān)聯(lián)。簡化的非線形模型很好的保留了原始模型的特性。該方法雖然是通過在一個確定的工況點求導(dǎo)來線性化模型，但在一個很寬的管網(wǎng)工作運(yùn)行條件范圍內(nèi)，模型的精度完全可以滿足要求。

四、管網(wǎng)簡化的實現(xiàn)

　　對上述化簡理論過程的實現(xiàn)如圖1所示，化簡程序如圖2所示。

五、應(yīng)用實例及結(jié)論

　　應(yīng)用上述管網(wǎng)簡化理論，以北方某特大城市實際復(fù)雜管網(wǎng)為例進(jìn)行管網(wǎng)模型簡化計算。該原始管網(wǎng)模型共有3439個節(jié)點，4514條管段，3個供水水源。按第12個時段的運(yùn)行工況點進(jìn)行簡化計算，去除了3407個節(jié)點，形成了僅有32個節(jié)點（其中包括10個水源節(jié)點）， 422條管段的新管網(wǎng)模型。限于篇幅，簡化前后的管網(wǎng)圖省略。通過對簡化前后的管網(wǎng)模型進(jìn)行平差模擬計算，發(fā)現(xiàn)簡化過程對管網(wǎng)計算的誤差影響很小。部分節(jié)點、管段的誤差分析見圖3、4。

　　

　　

　　本文提出的管網(wǎng)簡化算法，由于是程序自動處理，可以適應(yīng)各種拓?fù)浣Y(jié)構(gòu)、復(fù)雜程度的供水管網(wǎng)，可任意確定管網(wǎng)簡化程度，在簡化的過程中不需用戶任何干預(yù)處理，簡化精度完全滿足管網(wǎng)計算的要求。管網(wǎng)經(jīng)過簡化后，可大大減少供水管網(wǎng)優(yōu)化調(diào)度計算的耗時，克服了實現(xiàn)管網(wǎng)優(yōu)化調(diào)度實時性差的難題。

參考文獻(xiàn)：

　　給水管網(wǎng)理論和計算 中國建筑工業(yè)出版社 嚴(yán)煦世、趙洪賓

---

作者簡介：
高金良. 男 講師 哈爾濱工業(yè)大學(xué)市政環(huán)境工程學(xué)院
通 訊 處：0451-6282281 150090 哈爾濱工業(yè)大學(xué)（二區(qū)）624^#信箱
Email: gaojinliang@mail.hrbucea.edu.cn