龍?zhí)煊澹?王軍慧

(重慶大學(xué) 城市建設(shè)與環(huán)境工程學(xué)院,重慶400045)

供水管網(wǎng)的水力模型在工程應(yīng)用中比較多，水質(zhì)模型由于其復(fù)雜性，工程應(yīng)用還比較少。導(dǎo)致水質(zhì)模型的不確定性的因素眾多，不僅包括水力輸入?yún)?shù)：節(jié)點流量Q和管段水損系數(shù)C；還包括水質(zhì)輸入?yún)?shù)：主體水反應(yīng)系數(shù)kb和管壁反應(yīng)系數(shù)kw。對供水管網(wǎng)模型的不確定性研究，以水力的不確定性的研究見多，而水質(zhì)不確定研究的研究較少。Kapelan等[1]在僅考慮節(jié)點流量不確定性的條件下，提出了一個隨機(jī)優(yōu)化設(shè)計模型，并用遺傳算法進(jìn)行了求解。劉孟軍[2]用一階泰勒展開式導(dǎo)出的節(jié)點壓力和管段流量對節(jié)點流量的近似偏導(dǎo)矩

金溪[3]用線性化供水管網(wǎng)概率水力模型，研究了節(jié)點流量的隨機(jī)性對管網(wǎng)水力特性的影響，研究發(fā)現(xiàn)當(dāng)節(jié)點流量的變異系數(shù)較大時，進(jìn)行線性近似的隨機(jī)性分析會導(dǎo)致較大的誤差。Barkdoll[4]在一個較小的供水管網(wǎng)之上，分析穩(wěn)態(tài)流和非穩(wěn)態(tài)流條件下了節(jié)點流量對于水壓和水質(zhì)的影響，這里的節(jié)點流量是基本用水量，給定節(jié)點流量的分布特征是高斯分布和均勻分布。

分析以上研究以及相關(guān)研究可以發(fā)現(xiàn)，目前有關(guān)供水管網(wǎng)模型的不確定性的研究，主要針對節(jié)點流量的隨機(jī)性或輸入數(shù)據(jù)的誤差引起的管網(wǎng)系統(tǒng)模型參數(shù)的不確定性，而針對供水管網(wǎng)模型主要校核參數(shù)的不確定性對于水質(zhì)的影響的研究極少。為此，在假定Q(基本用水量)、C、kb和kw的隨機(jī)性都服從正態(tài)分布的前提下，使用蒙特卡羅隨機(jī)抽樣法，對所獲的輸入?yún)?shù)的抽樣值，應(yīng)用動態(tài)水質(zhì)模型計算供水管網(wǎng)的節(jié)點水質(zhì)，研究了基本用水量、水損系數(shù)、主體水反應(yīng)系數(shù)和管壁反應(yīng)系數(shù)單獨(dú)作用及共同作用下對供水管網(wǎng)節(jié)點水質(zhì)不確定性的影響，并給出了應(yīng)用實例。

1 研究方法

因為不確定性本身是一個概率的現(xiàn)象，輸入?yún)?shù)是隨機(jī)取值的，在供水管網(wǎng)的不確定性的分析中應(yīng)采用概率的方法。蒙特卡羅法(MCS)作為一種隨機(jī)抽樣算法，在供水管網(wǎng)模型的不確定性分析中得到了廣泛的應(yīng)用，并且認(rèn)為以MCS計算得到的模型輸出量不確定性是正確的標(biāo)準(zhǔn)，其他方法在不確定性研究中應(yīng)用時應(yīng)先與MCS計算的結(jié)果進(jìn)行比較以確定方法是否可行。MCS是一種用數(shù)值模擬來解決與隨機(jī)變量有關(guān)的實際工程問題的方法，其基本思路是，對一個具體的問題，首先建立描述該問題的概率分布模型，然后根據(jù)該概率模型進(jìn)行隨機(jī)抽樣生成一組隨機(jī)數(shù)，以隨機(jī)數(shù)作為輸入值計算問題的目標(biāo)解，重復(fù)抽樣計算過程得到目標(biāo)解的統(tǒng)計規(guī)律，并以此作為問題的解[5]。

2 案例分析

這里采用EPANET軟件安裝目錄下提供的案例管網(wǎng)Net2，Net2管網(wǎng)的拓?fù)浣Y(jié)構(gòu)如圖2所示。Net2管網(wǎng)中包含一個水源節(jié)點1和一個水塔，水源節(jié)點1的供水是斷續(xù)的，消毒劑都是從水源節(jié)點1進(jìn)入管網(wǎng)中。模型輸入?yún)?shù)：水力輸入?yún)?shù)節(jié)點流量(包括節(jié)點基本用水量和用水模式)和水損系數(shù)不做改變。水質(zhì)輸入?yún)?shù)認(rèn)為所有管道的主體水反應(yīng)系數(shù)相同，kb= -0.8 d-1；供水管網(wǎng)的管壁余氯衰減系數(shù)分成兩類，管徑DN300的管壁系數(shù)為kw1= -0.648 d-1，管徑DN200的管壁系數(shù)為kw2= -0.792 d-1(供水管網(wǎng)只有這兩種管徑的管段)，余氯反應(yīng)為一級反應(yīng)。對Net2管網(wǎng)執(zhí)行延時水力水質(zhì)分析的時間為55 h，水質(zhì)時間步長5 min，報告水質(zhì)時間步長1 h，鐘表計時開始時間12 am。

圖1 算法流程圖

圖2 Net2管網(wǎng)拓?fù)浣Y(jié)構(gòu)

為了研究輸入?yún)?shù)對于動態(tài)水質(zhì)的不確定性的影響，這里水源節(jié)點1的余氯濃度隨模擬時間是不斷變化的，水源節(jié)點1的源頭水質(zhì)為1.0 mg/L，水質(zhì)變化的時間模式如圖3。Net2管網(wǎng)用水節(jié)點的初始余氯濃度都設(shè)為0.3 mg/L。水源節(jié)點的余氯濃度隨時間不斷變化，且水源節(jié)點的供水是斷續(xù)的，顯然模擬時段內(nèi)各節(jié)點余氯濃度也是隨時間不斷變化的。

認(rèn)為Net2所有的輸入?yún)?shù)都是服從正太分布的，量化輸入?yún)?shù)的不確定性，輸入?yún)?shù)的均值為模型中的值，取水損系數(shù)的變異系數(shù)為0.05，節(jié)點基本用水量的變異系數(shù)為0.1，參考Pasha[7]所有反應(yīng)系數(shù)的變異系數(shù)取0.1。對Net2進(jìn)行抽樣模擬中，每次利用MCS生成的隨機(jī)數(shù)個數(shù)為節(jié)點數(shù)、管段數(shù)、主體水反應(yīng)系數(shù)(等于管段數(shù))和管壁反應(yīng)系數(shù)(等于管段數(shù))之和。為了保證余氯濃度值統(tǒng)計量的充分收斂，在每次不確定性分析中都進(jìn)行了3 000次的抽樣計算。對Net2的余氯模型進(jìn)行水質(zhì)不確定性分析，通過模擬計算可以得到整個時段內(nèi)節(jié)點余氯濃度的分布特征。

圖4給出了在所有輸入?yún)?shù)的共同影響下節(jié)點7在40時余氯濃度的標(biāo)準(zhǔn)差和MCS抽樣計算次數(shù)的關(guān)系。節(jié)點7余氯濃度的標(biāo)準(zhǔn)差在0.003左右，節(jié)點余氯濃度的不確定度很小，這也說明了動態(tài)水質(zhì)模型節(jié)點水質(zhì)的魯棒性。從節(jié)點余氯濃度的統(tǒng)計量和抽樣次數(shù)的關(guān)系可知，即便在對這么的變量的抽樣模擬時，在約1 000次的模擬計算時余氯濃度的統(tǒng)計值已經(jīng)開始收斂。同時由于MCS在量化不確定性因素的方法簡單，因此在評價諸多因素對供水管網(wǎng)模型的不確定性影響時，MCS不失為一種很好的方法。

圖3 水源節(jié)點1余氯時間模式

圖4 40時節(jié)點7余氯濃度標(biāo)準(zhǔn)差

圖5和6分別給出了節(jié)點7和33在40 h余氯濃度的箱形圖，其中7是靠近供水點節(jié)點，節(jié)點33是末端節(jié)點。從節(jié)點余氯濃度的箱形圖中可以看出，所有不確定輸入?yún)?shù)共同作用下，節(jié)點余氯的不確定性程度最大。單個不確定參數(shù)對節(jié)點7余氯濃度的影響程度從大到小依次是：kw、kb、Q和C，單個不確定參數(shù)對節(jié)點33余氯濃度的影響程度從大到小依次是：kb、kw、Q和C?？傮w來說余氯反應(yīng)系數(shù)造成的不確定性程度最大，因此在水質(zhì)模型校核中對主體水反應(yīng)系數(shù)和管壁反應(yīng)系數(shù)進(jìn)行校核，是建立準(zhǔn)確的水質(zhì)模型的關(guān)鍵所在。比較水力模型輸入?yún)?shù)基本用水量和水損系數(shù)對余氯濃度的不確定性影響程度，都有基本用水量的影響程度要大于水損系數(shù)，這一點和節(jié)點流量和水損系數(shù)對于節(jié)點壓力的不確定影響程度的大小關(guān)系是一樣的[8]。

圖5 節(jié)點7余氯濃度箱形圖

圖6 節(jié)點33余氯濃度箱形圖

圖7和圖8分別給出了在所有不確定輸入?yún)?shù)的作用下，節(jié)點7和節(jié)點33在整個模擬時段上余氯濃度的誤差圖。節(jié)點7的余氯濃度曲線，有4個波谷。這4個出現(xiàn)的時間對應(yīng)著節(jié)點1不供水的4個時間段，這時候節(jié)點7的用水是由水塔提供，節(jié)點7使用由水塔流過來的水，其水質(zhì)要比由水源直接流過來要差，這是因為水在水塔中停留余氯濃度在不斷衰減。同時也可以看出，在由水源節(jié)點供水向水塔供水轉(zhuǎn)變時節(jié)點7余氯濃度的標(biāo)準(zhǔn)差相對較大，節(jié)點7余氯濃度標(biāo)準(zhǔn)差的最大值出現(xiàn)在47時的0.013 9 mg/L。

節(jié)點33的余氯濃度曲線和節(jié)點7的余氯曲線形式完全不同，從模擬時間開始節(jié)點33的余氯濃度一直在下降，這是因為節(jié)點33靠近管網(wǎng)末端，其水齡比較大，從起始開始由水源節(jié)點1所供給的水流到節(jié)點33需要21.89 h(節(jié)點水齡值)。那么在這個時間之間，節(jié)點33的余氯濃度只與管網(wǎng)中初始濃度0.3 mg/L有關(guān)，而和水源節(jié)點1的余氯濃度無關(guān)，所以曲線也一直在下降。在水源節(jié)點1流至節(jié)點33之后的3個不供水的時間段，節(jié)點33的余氯濃度曲線有上升的趨勢，其中在最后一個上升趨勢很明顯。這個原因和節(jié)點7余氯濃度處于波谷的原因是一樣的，都是因為給節(jié)點的水質(zhì)供給是由一個水源切換到另外一個水源。同樣的對于末端節(jié)點33而言有，由水源節(jié)點供水向水塔供水轉(zhuǎn)變時節(jié)點33余氯濃度的標(biāo)準(zhǔn)差相對較大，節(jié)點33余氯濃度標(biāo)準(zhǔn)差的最大值出現(xiàn)在47時的0.016 2 mg/L。

圖7 節(jié)點7余氯濃度誤差

圖8 節(jié)點33余氯濃度誤差

比較節(jié)點7和節(jié)點33對應(yīng)時刻節(jié)點余氯濃度的標(biāo)準(zhǔn)差，大部分時刻節(jié)點33余氯濃度標(biāo)準(zhǔn)差都是要大于節(jié)點7，這可能與節(jié)點在管網(wǎng)中的相對位置有關(guān)，節(jié)點33位于管網(wǎng)末端，其節(jié)點水齡比較大，同一時刻流至末端節(jié)點的水受到的不確定因素的影響比前端水齡較小的節(jié)點要多，所以其受不確定性參數(shù)的影響較大。

和水力模型的不確定分析類似，水質(zhì)模型的不確定分析可以用水質(zhì)的不確定區(qū)間作為水質(zhì)模型計算的值，這樣更能準(zhǔn)確的把握供水管網(wǎng)中的水質(zhì)情況。通過節(jié)點水質(zhì)的分布情況，對供水管網(wǎng)的水質(zhì)可靠度進(jìn)行評估。在實際的供水管網(wǎng)模型的應(yīng)用中，如果發(fā)現(xiàn)在出廠水水質(zhì)滿足標(biāo)準(zhǔn)的情況下，監(jiān)測節(jié)點的余氯濃度值長時間小于水力模型計算的節(jié)點水質(zhì)不確定區(qū)間下限值，則供水管網(wǎng)中可能發(fā)生了如下情況：(1) 監(jiān)測設(shè)備故障，導(dǎo)致所監(jiān)測的余氯值與實際情況不符；(2) 供水管網(wǎng)發(fā)生污染事故，監(jiān)測節(jié)點的余氯濃度值無法保障供水安全。建議建模工作者可以將監(jiān)測節(jié)點水質(zhì)不確定區(qū)間的下限值作為管網(wǎng)事故預(yù)警的標(biāo)準(zhǔn)，當(dāng)監(jiān)測節(jié)點的監(jiān)測余氯濃度值頻繁小于標(biāo)準(zhǔn)值時發(fā)出預(yù)警信號，提醒相關(guān)工作人員對供水管網(wǎng)中存在的問題進(jìn)行核實和處理。同時水質(zhì)不確定區(qū)間可以為水廠消毒劑的投加劑量提供參考等。

3 結(jié)論

(1) 由輸入?yún)?shù)所導(dǎo)致動態(tài)水質(zhì)模型節(jié)點余氯濃度的不確定性是很小的，動態(tài)水質(zhì)模型具有魯棒性。

(2) 在水力和水質(zhì)不確定輸入?yún)?shù)的共同作用下，節(jié)點余氯濃度的不確定性最大；余氯反應(yīng)系數(shù)是對余氯濃度影響程度最大的輸入?yún)?shù)，在水質(zhì)模型校核時應(yīng)針對主體水反應(yīng)系數(shù)和管壁反應(yīng)系數(shù)進(jìn)行校核；水力輸入?yún)?shù)節(jié)點基本用水對水質(zhì)的不確定性影響的程度要高于水損系數(shù)。

(3) 同一時刻一般有末端節(jié)點余氯濃度的不確定性要大于前端節(jié)點，且當(dāng)節(jié)點用水由一水源切換到另外一水源時，節(jié)點余氯濃度的不確定性程度最大。

參考文獻(xiàn)：

[1]KAPELAN Z, SAVIC D, WALTERS G A. Robust least cost design of water distribution systems using Gas[J]. In Proceedings of the computer control for water industry CCWI, 2003(2):101-105

[2]劉孟軍,鄒平華,何鐘怡. 供水管網(wǎng)中節(jié)點流量隨機(jī)性的影響分析[J].給水排水，2007,33:350-354

[3]金溪. 供水管網(wǎng)概率水力模型線性化求解方法[J].河海大學(xué)學(xué)報，2011,39(4)：458-463

[4]BARKDOLL B D, DIDGAM H. Effect of user demand on water quality and hydraulics of distribution systems[J].Proc World Water & Environmental Resources Congress， 2003(6)：107-109

[5]王中宇,劉智敏，夏新濤. 測量誤差與不確定度評定[M]. 北京：科學(xué)出版社,2008

[6]謝善斌, 金波, 劉遂慶. 基于不確定度量的給水管網(wǎng)用戶的用水模式分析[J].中國給水排水,2011,27(7): 59-62

[7]PASHA M F K. Uncertainty Analysis and Calibration of Water Distribution Quality Models[D]. Tucson Arizona:Universityof Arizona,2006

[8]PASHA M F, LANSEY K. Analysis of uncertainty on waterdistribution hydraulics and water quality Proc[J]. World Water Congress, 2005(6)：207-209