曾 翰，陶 濤

(同濟(jì)大學(xué)環(huán)境與工程學(xué)院，上海 200092)

近年來(lái)，隨著國(guó)內(nèi)用水需求的不斷增大和水資源的日益短缺，供水管網(wǎng)漏損問(wèn)題越發(fā)突出。管網(wǎng)的物理漏損隨壓力的增大而增大，對(duì)管網(wǎng)進(jìn)行壓力管理、降低管網(wǎng)剩余壓力是降低管網(wǎng)漏失的有效途徑。為了加強(qiáng)對(duì)管網(wǎng)壓力及漏損的控制，可以采用安裝閥門或流量計(jì)的方法將供水管網(wǎng)劃分成若干個(gè)分區(qū)，再安裝調(diào)壓閥實(shí)現(xiàn)壓力控制。

關(guān)于供水管網(wǎng)的分區(qū)方法，目前國(guó)際上尚未有統(tǒng)一的定論，實(shí)踐中仍以經(jīng)驗(yàn)方法為主，也有部分學(xué)者開展了自動(dòng)分區(qū)方法的研究，常用的分區(qū)方法有聚類分析[1,5]、社區(qū)發(fā)現(xiàn)[2-3,5]、廣度優(yōu)先搜素[4-5]、最短路徑分析[6]等。國(guó)內(nèi)對(duì)分區(qū)的研究多側(cè)重于計(jì)量分區(qū)(DMA)，分區(qū)之間主要通過(guò)流量計(jì)隔離，而本文提出的分區(qū)方法，將管網(wǎng)劃分為干管管段和若干連接在干管上的分區(qū)，各分區(qū)之間通過(guò)閥門截?cái)喽恢苯酉噙B。對(duì)于這種分區(qū)結(jié)構(gòu)，各分區(qū)僅與干管相連，分區(qū)間水力聯(lián)系較少，在具備區(qū)域計(jì)量的功能之余，便于分別對(duì)各個(gè)區(qū)域進(jìn)行壓力控制。為了獲得這樣的分區(qū)方案，本文的方法首先從供水管網(wǎng)中定義干管，挑選干管上的若干支管作為起點(diǎn)，計(jì)算管網(wǎng)各個(gè)節(jié)點(diǎn)到這些支管的最短路徑，并把每個(gè)節(jié)點(diǎn)分配給最短路徑權(quán)重和最小的支管，每根支管和分配給它的節(jié)點(diǎn)即組成了一個(gè)分區(qū)，隨后，在滿足分區(qū)規(guī)模要求的前提下合并部分小分區(qū)，形成最終的分區(qū)方案。供水管網(wǎng)分區(qū)需在滿足供水可靠性的前提下，兼顧壓力管理、管網(wǎng)水質(zhì)安全性以及分區(qū)改造工程的經(jīng)濟(jì)性。因此，管網(wǎng)分區(qū)問(wèn)題是一個(gè)多目標(biāo)問(wèn)題，針對(duì)這一點(diǎn)，本文的方法以分區(qū)內(nèi)壓力標(biāo)準(zhǔn)差總和比、節(jié)點(diǎn)平均水齡以及管網(wǎng)改造費(fèi)用為目標(biāo)函數(shù)，采用非支配排序遺傳算法(NSGA-Ⅱ)來(lái)優(yōu)化支管的挑選過(guò)程，從而獲得分區(qū)方案的帕累托前沿。

1 多目標(biāo)供水管網(wǎng)分區(qū)方法

本文提出的多目標(biāo)供水管網(wǎng)分區(qū)方法的流程如圖1所示。

圖1 多目標(biāo)供水管網(wǎng)分區(qū)方法流程圖Fig.1 Flow Chart of Multi-Objective Optimization of Sectorization for Water Supply Distribution Network

1.1 獲取管網(wǎng)拓?fù)鋱D并設(shè)置參數(shù)

從EPANET的input文件中讀取管網(wǎng)數(shù)據(jù)，構(gòu)建管網(wǎng)鄰接矩陣A。鄰接矩陣的權(quán)值設(shè)置將在章節(jié)1.3.2進(jìn)行討論。

設(shè)置的參數(shù)包括分區(qū)最小規(guī)模、分區(qū)最大規(guī)模和節(jié)點(diǎn)服務(wù)水頭。關(guān)于DMA分區(qū)的規(guī)模，主要依據(jù)經(jīng)驗(yàn)確定，500～5 000戶是DMA分區(qū)的常用規(guī)模[4]。確定好分區(qū)的最小和最大規(guī)模后，可估算出分區(qū)的最大和最小數(shù)目。節(jié)點(diǎn)服務(wù)水頭可以參考相關(guān)地方規(guī)范進(jìn)行設(shè)置。

1.2 定義干管和候選分區(qū)入口

在實(shí)際管網(wǎng)中，干管和分配管的區(qū)別是比較模糊的，因?yàn)楦晒芤部梢猿袚?dān)配水的功能，而分配管也有可能向干管輸送流量?？梢钥紤]依據(jù)管徑來(lái)篩選干管，因?yàn)楣芏蔚墓軓皆酱?，則管段能夠輸送的流量越大，其作為輸水干管的重要性也越大。因此挑選管徑和流量都較大的管段作為干管。對(duì)于多水源管網(wǎng)，為了提高供水可靠性，發(fā)揮多水源的聯(lián)合調(diào)度功能，還可以挑選一些連接各個(gè)水源點(diǎn)的管段作為干管。

挑選出干管后，沿著干管找出干管上的所有支管，每根支管即視為一個(gè)候選的分區(qū)入口。挑選出的干管數(shù)量應(yīng)足夠多，從而保證干管上的候選入口數(shù)目大于分區(qū)最小數(shù)量，如果候選入口數(shù)不足，則應(yīng)進(jìn)一步挑選更多的干管。此外，并非所有支管都適合作為候選分區(qū)入口，當(dāng)支管所供給的節(jié)點(diǎn)較少(達(dá)不到分區(qū)最小規(guī)模)，且這些節(jié)點(diǎn)和管網(wǎng)的其他區(qū)域不相連時(shí)，該支管不作為候選分區(qū)入口，其供給的節(jié)點(diǎn)視為次要節(jié)點(diǎn)，不作為分區(qū)的考慮對(duì)象。

1.3 優(yōu)化分區(qū)方案

1.3.1 初始化

隨機(jī)生出第一代種群P0，種群由N個(gè)個(gè)體組成，每個(gè)個(gè)體的染色體包含一種分區(qū)入口挑選方案以及分區(qū)的合并等級(jí)。

分區(qū)入口挑選方案可以用一串長(zhǎng)度等于候選分區(qū)入口數(shù)目的二進(jìn)制數(shù)來(lái)表達(dá)，“1”表示某分區(qū)入口被選中，需在該入口安裝流量計(jì)監(jiān)測(cè)其流量；“0”表示某分區(qū)入口未被選中，需在該入口安裝閥門將其截?cái)唷?/p>

合并等級(jí)是用于限制各個(gè)分區(qū)合并程度的變量。由合并等級(jí)計(jì)算出分區(qū)的合并閾值，合并閾值是一個(gè)介于分區(qū)最小和最大規(guī)模的數(shù)值，如果兩個(gè)分區(qū)的總規(guī)模超出合并閾值，則兩個(gè)分區(qū)不合并，反之則合并。假設(shè)合并等級(jí)是一個(gè)三位二進(jìn)制數(shù)，則能夠表達(dá)的合并等級(jí)為0～7，合并閾值按式(1)計(jì)算。

(1)

1.3.2 節(jié)點(diǎn)分組

對(duì)于每一個(gè)個(gè)體，根據(jù)分區(qū)入口挑選方案，進(jìn)行節(jié)點(diǎn)分組。假定干管上一共有K個(gè)候選分區(qū)入口，某一個(gè)個(gè)體共有K1(K1≥K)個(gè)入口被選中，則給水管網(wǎng)中除干管節(jié)點(diǎn)和次要節(jié)點(diǎn)外的每個(gè)節(jié)點(diǎn)i都有K1個(gè)可能的供水源。對(duì)于一個(gè)給定的候選分區(qū)入口k和節(jié)點(diǎn)i，存在有限條連接這兩者的路徑，在這些路徑中，定義路徑權(quán)重和最小的路徑為入口k到節(jié)點(diǎn)i的最短路徑，其權(quán)重和定義為入口k到節(jié)點(diǎn)i的距離，如式(2)。

(2)

其中：{pathki}—入口k到節(jié)點(diǎn)i的最短路徑的管段集合；

wp—管段p的權(quán)值。

本文的最短路徑分析采用迪杰斯特拉(Dijstra)[7]算法。此處的最短路徑是對(duì)從入口k到節(jié)點(diǎn)i的最優(yōu)供水路徑的近似，考慮到分區(qū)前后管網(wǎng)拓?fù)浒l(fā)生了改變，管段流向、管段流量、管段水頭損失可能隨之改變，因此，本文選取分區(qū)前后均不發(fā)生變化的參數(shù)作為權(quán)值。EPANET中常采用海曾-威廉公式計(jì)算管段水頭損失，采用該公式對(duì)應(yīng)的管段阻力系數(shù)作為管段權(quán)值，假設(shè)管網(wǎng)共有n根管段，則管段權(quán)值的計(jì)算如式(3)。

(3)

其中：Cl—管段l的海曾-威廉粗糙系數(shù)；

dl—管段l的管徑，m；

Ll—管段l的管長(zhǎng)，m。

對(duì)于一個(gè)給定的節(jié)點(diǎn)i，計(jì)算它到K1個(gè)候選分區(qū)入口的距離，得到如式(4)的距離集合。

distancesi={distance1i,…,distanceki,…,

distanceK1i}

(4)

根據(jù)上述距離集合，取距離最小的入口為節(jié)點(diǎn)i的供水源。計(jì)算每個(gè)節(jié)點(diǎn)i到K1個(gè)入口的最短路徑，并確定每個(gè)節(jié)點(diǎn)i的供水源，每個(gè)供水源和分配給它的節(jié)點(diǎn)即組成了一個(gè)分區(qū)，從而把整個(gè)管網(wǎng)劃分為干管部分和K1個(gè)單入口供水的分區(qū)。

1.3.3 分區(qū)合并

經(jīng)過(guò)步驟1.3.2，從供水管網(wǎng)中劃分出個(gè)分區(qū)，每個(gè)分區(qū)都僅由一個(gè)入口供水，在這種情況下，如果分區(qū)較多，則分區(qū)邊界也較多，需在管網(wǎng)中截?cái)噍^多的管段，改造成本較高，且對(duì)壓力的影響也較大；如果分區(qū)較少，雖然能降低改造成本，但是每個(gè)入口的供水范圍都較大，難以保證每個(gè)分區(qū)的供水壓力達(dá)標(biāo)。因此引入分區(qū)合并過(guò)程，對(duì)小分區(qū)進(jìn)行合并，在降低改造成本之余，適當(dāng)增加部分分區(qū)的供水入口，避免出現(xiàn)分區(qū)壓力偏低的現(xiàn)象。

分析K1個(gè)分區(qū)的鄰接關(guān)系，如果兩個(gè)分區(qū)之間有管段相連，且兩者的規(guī)模總和不超過(guò)合并閾值，則合并這兩個(gè)分區(qū)。不斷重復(fù)這個(gè)過(guò)程，直至所有的分區(qū)都不能再合并為止，得到該個(gè)體的最終分區(qū)方案。

完成分區(qū)合并后，K1個(gè)被選上的入口即為分區(qū)的入口，在水力模型中設(shè)置這些管段的狀態(tài)為開啟；K～K1個(gè)未被選上的入口是需要安裝閥門截?cái)嗟墓芏?，在水力模型中設(shè)置這些管段的狀態(tài)為關(guān)閉；分區(qū)之間相連的管段即分區(qū)的邊界，需要安裝閥門截?cái)噙@些管段，在水力模型中設(shè)置這些管段的狀態(tài)為關(guān)閉，從而實(shí)現(xiàn)各個(gè)分區(qū)的相互隔離。

1.3.4 計(jì)算目標(biāo)函數(shù)

用EPANET對(duì)分區(qū)后的管網(wǎng)水力模型執(zhí)行延時(shí)模擬，考察分區(qū)后的管網(wǎng)表現(xiàn)，考慮的目標(biāo)函數(shù)包括壓力標(biāo)準(zhǔn)差總和比、節(jié)點(diǎn)平均水齡以及管網(wǎng)改造費(fèi)用。

(1)壓力標(biāo)準(zhǔn)差總和比

供水管網(wǎng)分區(qū)后的壓力均勻性是評(píng)價(jià)分區(qū)方案優(yōu)劣的重要指標(biāo)之一，節(jié)點(diǎn)壓力相差較小的分區(qū)具有更大的降壓潛力。以分區(qū)內(nèi)24 h的節(jié)點(diǎn)壓力標(biāo)準(zhǔn)差的均值表征分區(qū)內(nèi)的壓力均勻性，所有分區(qū)的節(jié)點(diǎn)壓力標(biāo)準(zhǔn)差總和代表整個(gè)管網(wǎng)分區(qū)后的壓力均勻性，總和越小，則均勻性越好。供水管網(wǎng)分區(qū)的目的是為了使管網(wǎng)局部的壓力變得更均勻，因此以管網(wǎng)分區(qū)后與分區(qū)前的分區(qū)壓力標(biāo)準(zhǔn)差總和之比作為目標(biāo)函數(shù)，假設(shè)經(jīng)過(guò)合并后，共有個(gè)分區(qū)，則按式(5)計(jì)算。

(5)

其中：nk2—第k2個(gè)分區(qū)的節(jié)點(diǎn)數(shù)目；

式(5)中，分區(qū)前的壓力是指不對(duì)水力模型進(jìn)行任何修改，直接運(yùn)行而得到的壓力。

(2)節(jié)點(diǎn)平均水齡

以管網(wǎng)24 h的節(jié)點(diǎn)平均水齡WA來(lái)表征管網(wǎng)的水質(zhì)，平均水齡越小，則說(shuō)明分區(qū)對(duì)管網(wǎng)水質(zhì)的影響越小，假設(shè)供水管網(wǎng)共有m個(gè)節(jié)點(diǎn)，則WA按式(6)計(jì)算。

(6)

(3)管網(wǎng)改造費(fèi)用

供水管網(wǎng)分區(qū)方案應(yīng)使安裝的水表和閥門數(shù)量都較少，以降低管網(wǎng)改造費(fèi)用，改造費(fèi)用C按式(7)計(jì)算。

(7)

K3—需要安裝閥門的管段總數(shù)。

1.3.5 非支配優(yōu)選

根據(jù)各個(gè)個(gè)體的目標(biāo)函數(shù)值和是否破壞約束條件進(jìn)行非支配排序，其中破壞約束條件的個(gè)體將被任一不破壞約束條件的個(gè)體所支配。本文考慮的約束條件包括以下幾點(diǎn)。

(1)供水范圍必須覆蓋所有節(jié)點(diǎn)。

(2)分區(qū)規(guī)模約束，如式(8)。

sizemin

(8)

(3)節(jié)點(diǎn)壓力約束，如式(9)。

(9)

(4)水力約束，如式(10)。

G(h,q)=0

(10)

其中：size(dmak)—第k個(gè)分區(qū)的規(guī)模；

sizemin,sizemax—允許的分區(qū)最小、最大規(guī)模；

Hmin—允許的最小節(jié)點(diǎn)壓力，m；

h,q—節(jié)點(diǎn)水頭向量和節(jié)點(diǎn)需水量向量；

G—水力模型非線性方程組，包括流量連續(xù)性方程和環(huán)能量方程組。

按非支配等級(jí)從高到低挑選個(gè)體直到個(gè)體數(shù)超過(guò)N，在最后一個(gè)非支配等級(jí)中的個(gè)體挑選中，擁擠距離大的個(gè)體將被優(yōu)先選擇。N個(gè)被選擇的個(gè)體組成第二代父代種群P1。

1.3.6 循環(huán)

父代種群P1通過(guò)選擇、交叉和變異得到第二代子代種群Q1，重復(fù)1.3.2～1.3.5的過(guò)程得到P2，不斷循環(huán)直到滿足算法終止條件為止，得到第t代種群Pt。種群Pt即為Pareto最優(yōu)解或近似最優(yōu)解。

2 算例分析

2.1 算例概況

選取ZZ市供水管網(wǎng)的水力模型作為實(shí)例研究對(duì)象，該水力模型包含17 282個(gè)節(jié)點(diǎn)，18 518根管段，4個(gè)水源，總供水量為11 610 L/s，9個(gè)流量控制閥，管網(wǎng)拓?fù)淙鐖D2所示。

圖2 ZZ市供水管網(wǎng)水力模型Fig.2 Hydraulic Model of Water Supply Distribution Network in ZZ City

挑選部分管徑大于等于800的管段作為干管，為了保證供水的可靠性，挑選了一些連接各個(gè)水源點(diǎn)的管段作為干管，并根據(jù)步驟1.2所述識(shí)別出連接于干管上的次要節(jié)點(diǎn)，如圖3中的黑色管段所示。

預(yù)計(jì)劃分10～100個(gè)分區(qū)，由于缺乏ZZ市管網(wǎng)的節(jié)點(diǎn)用戶數(shù)量信息，本文以分區(qū)的節(jié)點(diǎn)流量總和表征分區(qū)的規(guī)模，設(shè)置分區(qū)最小規(guī)模和最大規(guī)模分別為116.1 L/s和1 161 L/s。節(jié)點(diǎn)服務(wù)水頭參考最新的規(guī)范，設(shè)置為16 m。由于ZZ市局部地勢(shì)較高，該市的管網(wǎng)在分區(qū)前就有部分節(jié)點(diǎn)的水壓達(dá)不到16 m，預(yù)計(jì)分區(qū)后這些節(jié)點(diǎn)的水壓也難以達(dá)標(biāo)，因此算法的節(jié)點(diǎn)壓力約束如下設(shè)置：對(duì)于分區(qū)前節(jié)點(diǎn)水壓不小于16 m的節(jié)點(diǎn)，約束其水壓在分區(qū)后不小于16 m；對(duì)于分區(qū)前節(jié)點(diǎn)水壓小于16 m的節(jié)點(diǎn)，約束其水壓在分區(qū)后大于0。在運(yùn)算前設(shè)置9個(gè)流量控制閥的狀態(tài)為開啟，即閥門全開而不限制管段流過(guò)的流量，同時(shí)也不產(chǎn)生水頭損失，流量控制閥在鄰接矩陣A中對(duì)應(yīng)的管段的權(quán)值設(shè)置為0。水力模型的延時(shí)模擬時(shí)間為48 h，讀取后24 h的節(jié)點(diǎn)壓力和節(jié)點(diǎn)水齡數(shù)據(jù)，每小時(shí)讀取一次。

圖3 ZZ市供水管網(wǎng)的干管及次要節(jié)點(diǎn)Fig.3 Water Mains and Minor Nodes of Water Supply Distribution Network in ZZ City

計(jì)算管網(wǎng)改造成本所采用的流量計(jì)與閥門的價(jià)格，如表1所示。

表1 閥門及流量計(jì)參考價(jià)格

采用matlab工具箱自帶的多目標(biāo)遺傳算法實(shí)現(xiàn)多目標(biāo)優(yōu)化，設(shè)置個(gè)體數(shù)為100，代數(shù)為500，其余參數(shù)采用默認(rèn)值。

2.2 計(jì)算結(jié)果

采用單臺(tái)i7核心、16 GB內(nèi)存的臺(tái)式計(jì)算機(jī)運(yùn)算該算例，36 h運(yùn)算結(jié)束，得到的Pareto前沿共包含32個(gè)解，壓力標(biāo)準(zhǔn)差總和比在0.92～1.22，節(jié)點(diǎn)平均水齡在20.4～23.2 h，改造費(fèi)用在364.6～504.7萬(wàn)元，分區(qū)方案的入口數(shù)為45～97個(gè)，截?cái)喙芏螖?shù)為112～168根，分區(qū)數(shù)目為12～16個(gè)。從分區(qū)方案的入口數(shù)和分區(qū)數(shù)目的范圍來(lái)看，平均每個(gè)分區(qū)都有3～6個(gè)入口，即大多數(shù)分區(qū)都有一根或一根以上的管段與干管相連。這可能是因?yàn)?，原始管網(wǎng)本身就存在壓力達(dá)不到16 m的節(jié)點(diǎn)，使得部分區(qū)域降壓空間較小，因此必須設(shè)法增加這些區(qū)域與干管的連接管段以獲得區(qū)域壓力符合約束條件的解。

根據(jù)所得結(jié)果可知Pareto前沿為一三維曲面，32個(gè)解大致均勻地分散在其中，三個(gè)目標(biāo)函數(shù)值兩兩間沒有呈現(xiàn)出顯著的趨勢(shì)和相關(guān)性，將32個(gè)解按管網(wǎng)改造費(fèi)用劃分為5組，如圖4所示。由圖4可知：當(dāng)改造費(fèi)用較高時(shí)，分區(qū)方案的節(jié)點(diǎn)平均水齡和壓力均勻性都較小(如圖中的紫色和綠色點(diǎn)，主要聚集在左下角，對(duì)應(yīng)的改造費(fèi)用為450～480萬(wàn)元和480～510萬(wàn)元)；把解分為5組后，每一組解中，壓力標(biāo)準(zhǔn)差總和比與節(jié)點(diǎn)平均水齡大致呈負(fù)相關(guān)。

圖4 多目標(biāo)管網(wǎng)分區(qū)方法的Pareto前沿Fig.4 Pareto Front of Multi-Objective Optimization of Sectorization for Water Supply Distribution Network

表2羅列了32個(gè)解中10個(gè)解的3個(gè)目標(biāo)函數(shù)值以及分區(qū)后的分區(qū)數(shù)、用水最高時(shí)平均壓力、用水最低時(shí)平均壓力、全管網(wǎng)最低壓力、全管網(wǎng)最高壓力以及低壓節(jié)點(diǎn)(壓力低于16 m的節(jié)點(diǎn))數(shù)目。由表2可知，大多數(shù)分區(qū)方案的節(jié)點(diǎn)平均水齡低于原始管網(wǎng)，這可能是因?yàn)樵诜謪^(qū)前，由于各區(qū)域間有較多的連接管相連，每個(gè)區(qū)域相當(dāng)于由多個(gè)入口供水，而分區(qū)后，每個(gè)區(qū)域由較少的幾個(gè)入口供水，導(dǎo)致整個(gè)區(qū)域的水量集中由少數(shù)的幾根管道供應(yīng)，使得這些管道流量流速增大，水齡減小，從而使得整個(gè)管網(wǎng)的節(jié)點(diǎn)平均水齡有所減少。此外，分區(qū)后管網(wǎng)用水最低時(shí)的平均壓力相比原始管網(wǎng)有所上升，而最高時(shí)的平均壓力有所下降，說(shuō)明分區(qū)對(duì)最高時(shí)工況的影響較大。盡管如此，分區(qū)后全管網(wǎng)最低壓力均高于分區(qū)前，且分區(qū)后低壓節(jié)點(diǎn)數(shù)比分區(qū)前少，說(shuō)明分區(qū)一定程度上能夠改善管網(wǎng)的低壓情況。壓力標(biāo)準(zhǔn)差總和比的數(shù)值均不小于0.92，部分解的數(shù)值大于1，說(shuō)明分區(qū)后的分區(qū)內(nèi)壓力均勻性的改善程度并不高，甚至沒有改善，這可能是因?yàn)樵诒舅惴ǖ姆謪^(qū)方案中，各分區(qū)僅和干管相連，分區(qū)之間不相連，因此分區(qū)的供水入口較少，加劇了分區(qū)壓力的不均勻。

表2 10個(gè)解的目標(biāo)函數(shù)值以及分區(qū)后的壓力表現(xiàn)

好的分區(qū)方案需兼顧壓力均勻性、水質(zhì)安全性以及經(jīng)濟(jì)性，解s-3的改造費(fèi)用為413.5萬(wàn)元，在32個(gè)解中處于較低的水平，節(jié)點(diǎn)平均水齡為23.1 h，相比分區(qū)前的23.0 h變化不大，壓力標(biāo)準(zhǔn)差總和比為0.95，相比分區(qū)前有小幅減小，說(shuō)明局部壓力均勻性有一定的改善，全管網(wǎng)最低壓力和低壓節(jié)點(diǎn)數(shù)分別為4.4 m和612，說(shuō)明低壓狀況相比分區(qū)前(3.0 m和724)也有所改善，可見解s-3能較好地兼顧改造成本、管網(wǎng)水質(zhì)和管網(wǎng)壓力，該解對(duì)應(yīng)的分區(qū)方案如圖5所示。

圖5 解s-3的分區(qū)方案Fig.5 Layout of Sectorization for Water Supply Distribution Network of Solution s-3

3 結(jié)論

通過(guò)上述的算例分析可知，本文的算法能夠得到分區(qū)方案的Pareto前沿，管理者可以根據(jù)管網(wǎng)實(shí)際運(yùn)行條件和當(dāng)?shù)亟?jīng)濟(jì)水平從Pareto前沿中選擇合適的分區(qū)方案。本算法的分區(qū)方案一定程度上能夠改善管網(wǎng)局部的壓力均勻性，且能夠優(yōu)化管網(wǎng)水齡和改善管網(wǎng)低壓情況。

然而，本文所采用的算例中不包含泵站，閥門為流量控制閥，因此在管段權(quán)值的設(shè)置方面未能對(duì)泵站和其他閥門情況進(jìn)行討論。本算法也未考慮節(jié)點(diǎn)用戶類型，分區(qū)方案中有可能把幾種類型的用戶劃分到一個(gè)區(qū)域中。因此算法中的權(quán)值設(shè)置、節(jié)點(diǎn)分組和分區(qū)合并策略需更深入的研究。

此外，算例計(jì)算結(jié)果表明，大多數(shù)的分區(qū)都有一根或一根以上的管段與干管相連，在這種情況下，有可能出現(xiàn)某些管段是分區(qū)入口，而某些管段是分區(qū)出口的現(xiàn)象，給壓力管理帶來(lái)不便。盡管如此，平均每個(gè)分區(qū)的入口數(shù)較少，只需監(jiān)測(cè)較少的管段流量就能夠監(jiān)測(cè)分區(qū)的流量，可以考慮在某些分區(qū)內(nèi)加設(shè)泵站，提高分區(qū)壓力，從而進(jìn)一步減少分區(qū)的入口數(shù)。

[1]葉健.結(jié)合圖論和評(píng)價(jià)體系的城市供水管網(wǎng)PMA分區(qū)優(yōu)化研究[D].哈爾濱: 哈爾濱工業(yè)大學(xué),2015.

[2]DIAO K,ZHOU Y,RAUCH W.Automated creation of district metered area boundaries in water distribution systems[J].Journal of Water Resources Planning & Management Asce,2013,139(2):184-190.

[3]CAMPBELL E,IZQUIERDO J,MONTALVO I,et al.A novel water supply network sectorization methodology based on a complete economic analysis,including uncertainties[J].Water,2016,8(5):179.

[4]HAJEBI S,ROSHANI E,CARDOZO N,et al.Water distribution network sectorisation using graph theory and many-objective optimisation[J].Journal of Hydroinformatics,2015,18(1): 77-95.

[5]PERELMAN L S,ALLEN M,PREIS A,et al.Automated sub-zoning of water distribution systems[J].Environmental Modelling & Software,2015,19(3):1-14.

[6]ZHENG F,SIMPSON A R,ZECCHIN A C.A decomposition and multistage optimization approach applied to the optimization of water distribution systems with multiple supply sources[J].Water Resources Research,2013,49(1):380-399.

[7]NARDO A D,NATALE M D,SANTONASTASO G F,et al.Water network sectorization based on genetic algorithm and minimum dissipated power paths[J].Water Science & Technology Water Supply,2013,13(4): 951-957.