• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于負(fù)梯度法的供水管網(wǎng)污染源識(shí)別

    2013-10-30 08:15:16信昆侖盛希夫項(xiàng)寧銀
    關(guān)鍵詞:差分法污染源步長(zhǎng)

    信昆侖,盛希夫,陶 濤,項(xiàng)寧銀

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

    近幾年,城市供水管網(wǎng)的突發(fā)污染事故屢有發(fā)生,除水源地污染引起管網(wǎng)飲用水問(wèn)題之外,如管道錯(cuò)接等供水管網(wǎng)內(nèi)部的原因也會(huì)造成水質(zhì)的局部乃至全局惡化.因此,如何快速準(zhǔn)確地定位污染源是飲用水安全保障急需解決的課題.關(guān)于污染源的定位研究,Laird等[1]提出了用非線性規(guī)劃方法來(lái)求解這一問(wèn)題.首先用源頭追蹤算法替代EPANET水質(zhì)模擬計(jì)算引擎從而產(chǎn)生必要的數(shù)據(jù),然后通過(guò)帶正則項(xiàng)的非線性規(guī)劃方法求出污染源侵入的時(shí)間和地點(diǎn).針對(duì)這一方法解的非唯一性問(wèn)題,Laird等[2]又提出了先用混合整數(shù)二次規(guī)劃的方法確定一個(gè)候選污染源節(jié)點(diǎn)集,然后在此污染源節(jié)點(diǎn)集的基礎(chǔ)上利用前一方法篩選出最有可能是污染源的節(jié)點(diǎn).但這個(gè)方法對(duì)于大中型管網(wǎng)的適用性和此模型與實(shí)際管網(wǎng)相符程度仍然有待研究.Guan等[3]利用EPANET作為內(nèi)嵌的水質(zhì)模擬器再結(jié)合負(fù)梯度法提出了模擬優(yōu)化算法.這一方法在已知管網(wǎng)模型的運(yùn)行模式和水質(zhì)監(jiān)測(cè)數(shù)據(jù)充足可靠的前提下能夠較好地解出污染源侵入的時(shí)間和地點(diǎn).但只考慮了水質(zhì)反應(yīng)為線性的情況,也沒(méi)有給出確定候選污染源節(jié)點(diǎn)的方法.Pries等[4]針對(duì)此問(wèn)題提出了關(guān)系樹(shù)-線性規(guī)劃算法.通過(guò)構(gòu)造上下游節(jié)點(diǎn)的污染物濃度關(guān)系樹(shù)并形成線性規(guī)則,再利用關(guān)系樹(shù)的線性規(guī)則反向求解線性規(guī)劃問(wèn)題,得到污染源節(jié)點(diǎn)位置及污染注入屬性.Cristo等[5]提出了離散最優(yōu)化的方法來(lái)求解污染源的識(shí)別問(wèn)題.主要通過(guò)構(gòu)造水質(zhì)污染矩陣來(lái)確定可能的污染源注入點(diǎn),再通過(guò)最小化模擬值和監(jiān)測(cè)值的差來(lái)從可能的污染源點(diǎn)中確定污染源的位置.另外,Liu等[6]同樣利用EPANET作為內(nèi)嵌的水質(zhì)模擬計(jì)算引擎用遺傳算法隨機(jī)搜索最優(yōu)解.由于以上方法均需要管網(wǎng)的水質(zhì)監(jiān)測(cè)點(diǎn)可以報(bào)告實(shí)時(shí)的污染物濃度,對(duì)于水質(zhì)監(jiān)測(cè)點(diǎn)僅能判斷是否出現(xiàn)污染物的情況,Yang等[7]提出用貝葉斯方法計(jì)算已檢出污染物的水質(zhì)監(jiān)測(cè)節(jié)點(diǎn)其上游節(jié)點(diǎn)為污染源的概率,然后根據(jù)概率的大小判斷出可能的污染源.

    本質(zhì)上來(lái)說(shuō),供水管網(wǎng)污染源的識(shí)別問(wèn)題是管網(wǎng)水質(zhì)模擬的反演問(wèn)題,從數(shù)學(xué)的角度來(lái)看也是一個(gè)優(yōu)化問(wèn)題.本文基于Guan等[3]利用負(fù)梯度法求解供水管網(wǎng)污染源識(shí)別問(wèn)題,結(jié)合Cristo等[5]提出的候選污染源節(jié)點(diǎn)集合方法,進(jìn)一步考慮供水管網(wǎng)水質(zhì)反應(yīng)的復(fù)雜性,用差分步長(zhǎng)法代替權(quán)重系數(shù)法求解梯度方向,采用EPANET作為內(nèi)嵌的水質(zhì)模擬計(jì)算引擎結(jié)合負(fù)梯度法來(lái)求解污染源侵入的時(shí)間、地點(diǎn)和侵入過(guò)程,并結(jié)合實(shí)例管網(wǎng)進(jìn)行算法求解性能的分析.

    1 模型構(gòu)建

    在城市供水管網(wǎng)擁有相當(dāng)精度的水力水質(zhì)模型且管網(wǎng)上安裝了若各個(gè)水質(zhì)監(jiān)測(cè)點(diǎn)的情況下,污染源的識(shí)別問(wèn)題可以轉(zhuǎn)化為一個(gè)數(shù)學(xué)上的最優(yōu)化問(wèn)題.本文以模擬時(shí)段內(nèi)各個(gè)水質(zhì)監(jiān)測(cè)點(diǎn)的實(shí)際污染物濃度值與模擬值差的平方和為目標(biāo)函數(shù),通過(guò)分析水質(zhì)監(jiān)測(cè)點(diǎn)不同時(shí)刻模擬濃度值與實(shí)測(cè)濃度值差來(lái)修正計(jì)算候選污染源在特定時(shí)刻的節(jié)點(diǎn)濃度.待模擬時(shí)段內(nèi)所有的候選污染源節(jié)點(diǎn)濃度更新后,再次進(jìn)行水質(zhì)模擬.當(dāng)模擬結(jié)果達(dá)到預(yù)設(shè)的條件后,輸出候選污染源節(jié)點(diǎn)濃度矩陣.由于求解所得的候選污染源污染物累計(jì)注入量與其作為真實(shí)污染源的可能性成正比,因此可以此作為確定污染源節(jié)點(diǎn)的依據(jù).

    首先定義

    式中:xi(t)表示候選污染源節(jié)點(diǎn)i在t時(shí)刻的污染物濃度值,對(duì)于k個(gè)時(shí)間段,則xi是一個(gè)行向量,即xi=[xi(t1),xi(t2),…,xi(tk)];yj(t)表示水質(zhì)監(jiān)測(cè)點(diǎn)j在t時(shí)刻的模擬污染物濃度值;n表示候選污染源節(jié)點(diǎn)總數(shù);m表示布置的水質(zhì)監(jiān)測(cè)點(diǎn)總數(shù);t表示模擬時(shí)刻.

    監(jiān)測(cè)點(diǎn)模擬濃度可以表示為

    結(jié)合監(jiān)測(cè)點(diǎn)的實(shí)際監(jiān)測(cè)結(jié)果,污染源識(shí)別問(wèn)題的目標(biāo)函數(shù)可以表示為[3]

    式中:X*為污染源識(shí)別問(wèn)題的最優(yōu)解,y′i(t)為水質(zhì)監(jiān)測(cè)點(diǎn)i在t時(shí)刻的實(shí)際濃度值.為了便于求解,這里需作幾點(diǎn)假設(shè):①污染物的流量相對(duì)于管網(wǎng)水量而言可忽略不計(jì),不影響管網(wǎng)原來(lái)的水力工況;②污染物在管網(wǎng)節(jié)點(diǎn)處投加,投加點(diǎn)可以是多個(gè),投加量可以隨時(shí)間變化,但在每一水質(zhì)時(shí)間步長(zhǎng)之內(nèi)是恒定不變的;③污染物侵入管網(wǎng)中每個(gè)節(jié)點(diǎn)的概率是相同的,也就是說(shuō)污染物是隨機(jī)侵入的;④管網(wǎng)中的水質(zhì)監(jiān)測(cè)點(diǎn)能夠監(jiān)測(cè)到任意濃度的污染物.

    2 模型求解

    由于式(3)表示的最優(yōu)化問(wèn)題目標(biāo)函數(shù)和約束條件均為非線性函數(shù),且約束條件是由供水管網(wǎng)恒定流方程、連續(xù)性方程、管道水質(zhì)反應(yīng)動(dòng)力學(xué)方程共同決定的復(fù)雜非線性等式約束,因此該問(wèn)題為高度復(fù)雜的非線性規(guī)劃問(wèn)題.考慮到約束條件可采用成熟的EPANET水力水質(zhì)分析模塊進(jìn)行隱式求解,故基于負(fù)梯度法,通過(guò)差分法求目標(biāo)函數(shù)梯度向量對(duì)此非線性規(guī)劃問(wèn)題進(jìn)行求解.求解過(guò)程如下:

    步驟1 給定初始近似點(diǎn)X(0)及精度ε>0,若‖▽f(X(0))‖2≤ε,則X(0)即為近似極小點(diǎn).需要說(shuō)明的是梯度向量的每一個(gè)元素均為目標(biāo)函數(shù)對(duì)X=[x1,x2,…,xn]相應(yīng)元素的偏導(dǎo)數(shù).該偏導(dǎo)數(shù)通過(guò)差分法近似求解

    式中:ΔX是第i個(gè)元素為δ且其余元素為0的列向量;δ為差分步長(zhǎng),其取值大小應(yīng)該合理確定;f即目標(biāo)函數(shù)(見(jiàn)式(3)).當(dāng)水質(zhì)反應(yīng)為線性關(guān)系時(shí),水質(zhì)監(jiān)測(cè)點(diǎn)所測(cè)得的污染物濃度值可以表達(dá)成污染物侵入節(jié)點(diǎn)的污染物濃度加權(quán)值[3],即

    偏導(dǎo)數(shù)的求法可以簡(jiǎn)化為

    式中:aj,i為水質(zhì)監(jiān)測(cè)點(diǎn)j對(duì)候選污染源節(jié)點(diǎn)i的權(quán)重系數(shù),即當(dāng)候選污染源節(jié)點(diǎn)i有一個(gè)單位的增量時(shí),在水質(zhì)監(jiān)測(cè)點(diǎn)j有相應(yīng)的污染物濃度增量.權(quán)重系數(shù)可以通過(guò)EPANET水質(zhì)分析中的源頭追蹤功能得到.

    步驟2 若‖▽f(X(0))‖2>ε,求步長(zhǎng)λ0,并計(jì)算

    求最佳步長(zhǎng)本文采用黃金分割法,最佳步長(zhǎng)即為使得目標(biāo)函數(shù)最小的步長(zhǎng).若‖▽f(X(1))‖2≤ε,則輸出X(1),否則轉(zhuǎn)入步驟3.

    步驟3 一般地說(shuō),若‖▽f(X(k))‖2≤ε,則X(k)即為所求的近似解;若‖▽f(X(k))‖2>ε,則求步長(zhǎng)λk,并確定下一個(gè)近似點(diǎn)

    如此迭代下去,直至達(dá)到所設(shè)定精度要求為止.

    3 算例研究

    以供水管網(wǎng)實(shí)例Net3為例進(jìn)行污染源的識(shí)別問(wèn)題求解.Net3管網(wǎng)模型含節(jié)點(diǎn)數(shù)92個(gè),水庫(kù)2座,網(wǎng)中水塔3個(gè),管段數(shù)117個(gè),水泵2臺(tái),如圖1所示.圖中虛線范圍內(nèi)的節(jié)點(diǎn)為候選污染源節(jié)點(diǎn).

    圖1 管網(wǎng)拓?fù)浣Y(jié)構(gòu)圖Fig.1 Topology layout of the case study network

    假定在節(jié)點(diǎn)101和119同時(shí)注入污染物,注入方式為設(shè)置點(diǎn)注入[8],注入時(shí)間段為3∶00am~10∶00am,污染物注入質(zhì)量濃度變化模式見(jiàn)圖2.管網(wǎng)水質(zhì)模擬時(shí)間總長(zhǎng)24h,水力時(shí)間步長(zhǎng)1h,水質(zhì)時(shí)間步長(zhǎng)5min.在管網(wǎng)中隨機(jī)選取6個(gè)節(jié)點(diǎn)作為水質(zhì)監(jiān)測(cè)點(diǎn)(分別為連接節(jié)點(diǎn)121,151,111,267,205,213),如圖1所示.通過(guò)水質(zhì)模擬,以這些節(jié)點(diǎn)模擬污染物質(zhì)量濃度數(shù)據(jù)作為監(jiān)測(cè)結(jié)果數(shù)據(jù)(見(jiàn)表1).

    表1 監(jiān)測(cè)點(diǎn)模擬污染物質(zhì)量濃度Tab.1 Simulated concentrations at monitoring sites mg·L-1

    3.1 確定候選節(jié)點(diǎn)集

    候選節(jié)點(diǎn)集是指根據(jù)水質(zhì)監(jiān)測(cè)點(diǎn)對(duì)污染物的監(jiān)測(cè)情況所確定的可能是管網(wǎng)污染事故污染源發(fā)生位置的節(jié)點(diǎn)集合.目前關(guān)于管網(wǎng)突發(fā)污染事故污染源追蹤的研究,仍然以通過(guò)經(jīng)驗(yàn)法確定候選節(jié)點(diǎn)集合為主[3],候選節(jié)點(diǎn)集合確定的范圍和準(zhǔn)確程度都會(huì)影響到污染源識(shí)別問(wèn)題最優(yōu)化算法的性能.為此,本研究根據(jù)文獻(xiàn)[5,9]提出的最小覆蓋集方法進(jìn)行候選節(jié)點(diǎn)集合的確定.首先求出節(jié)點(diǎn)的監(jiān)測(cè)時(shí)間矩陣,以24h為它們的服務(wù)水平,可以確定出每個(gè)監(jiān)測(cè)點(diǎn)所覆蓋的范圍.通過(guò)表1可知,本次污染事故假定模式下,實(shí)際上僅水質(zhì)監(jiān)測(cè)點(diǎn)151,111,267,205,213能監(jiān)測(cè)到污染物.因此可以把5個(gè)監(jiān)測(cè)點(diǎn)共同覆蓋的節(jié)點(diǎn)作為候選節(jié)點(diǎn)集(如圖1所示),即

    3.2 結(jié)果分析

    考慮到按照負(fù)梯度方向搜索得到的是局部最優(yōu)解,且絕大多數(shù)候選污染源節(jié)點(diǎn)均不是污染源注入點(diǎn),可將初始解向量的每個(gè)元素均設(shè)置為0.求解過(guò)程中,差分步長(zhǎng)δ取0.0100,最高迭代次數(shù)取360.圖2中顯示采用差分法計(jì)算梯度向量求解得到的結(jié)果,由于求解的多數(shù)候選點(diǎn)的質(zhì)量濃度值均很小,圖中僅顯示30個(gè)候選污染源節(jié)點(diǎn)中所求解的注入質(zhì)量濃度較大的4個(gè)節(jié)點(diǎn)的注入質(zhì)量濃度隨時(shí)間變化情況.候選污染源節(jié)點(diǎn)101是各時(shí)刻累計(jì)投加質(zhì)量濃度值最大的節(jié)點(diǎn),初步選定為污染源節(jié)點(diǎn).同時(shí),候選污染源節(jié)點(diǎn)10從3∶00am到10∶00am基本上以100mg·L-1的質(zhì)量濃度持續(xù)注入.圖1表明,節(jié)點(diǎn)10和節(jié)點(diǎn)101是同一個(gè)管段的起止節(jié)點(diǎn).該現(xiàn)象表明有一個(gè)污染源在以節(jié)點(diǎn)10和節(jié)點(diǎn)101為起止節(jié)點(diǎn)的管道中.事實(shí)上,候選污染源節(jié)點(diǎn)101的注入質(zhì)量濃度變化曲線基本與實(shí)際污染源節(jié)點(diǎn)101的質(zhì)量濃度變化曲線相吻合.此外,候選污染源節(jié)點(diǎn)119的投加質(zhì)量濃度呈拋物線型變化也與實(shí)際的污染源119投加質(zhì)量濃度變化相吻合,可以確定這里也是一個(gè)污染源注入點(diǎn).另外,候選污染源節(jié)點(diǎn)151作為污染源節(jié)點(diǎn)119的鄰近節(jié)點(diǎn),也有一定質(zhì)量濃度的污染物注入,這就更表明了另一個(gè)污染源就是節(jié)點(diǎn)119或者在其附近.

    圖2 污染物注入質(zhì)量濃度變化曲線(差分法)Fig.2 Injection concentrations of contamination source(finite difference method)

    3.3 影響求解結(jié)果的因素分析

    (1)差分步長(zhǎng)

    對(duì)于圖1所顯示的小型算例管網(wǎng),將迭代次數(shù)均設(shè)定為360次,差分步長(zhǎng)δ依次取值為0.1000,0.0100,0.0010和0.0001時(shí),計(jì)算結(jié)果的梯度向量模長(zhǎng)值依次為16.52,0.24,0.15和0.33.為了定量地衡量計(jì)算結(jié)果與實(shí)際結(jié)果的擬合程度,計(jì)算污染源節(jié)點(diǎn)101和119各個(gè)時(shí)段假定注入濃度和解出注入濃度差的平方和總和,依次為:1564.09,1914.74,7703.17,64201.58,如表2所示.

    由表2可知,差分步長(zhǎng)取0.0100時(shí)對(duì)應(yīng)的計(jì)算結(jié)果(梯度向量模長(zhǎng)為0.24,目標(biāo)函數(shù)值為1914.74,而梯度向量模長(zhǎng)和目標(biāo)函數(shù)值是反映梯度法求解質(zhì)量的兩個(gè)主要指標(biāo))最好,此時(shí)差分步長(zhǎng)取值正好與EPANET水質(zhì)選項(xiàng)中設(shè)定的水質(zhì)公差值相同.其原因在于水質(zhì)公差表示管網(wǎng)水質(zhì)的最小變化值,也即一部分水體與另一部分水體水質(zhì)的最小差值[9],這個(gè)值也表征了管網(wǎng)中水質(zhì)監(jiān)測(cè)點(diǎn)的污染物最低監(jiān)測(cè)限度,因此當(dāng)差分步長(zhǎng)取水質(zhì)公差值時(shí),計(jì)算結(jié)果最好.若差分步長(zhǎng)繼續(xù)減小,由于梯度計(jì)算將逐步超出管網(wǎng)模型的計(jì)算精度范圍,計(jì)算結(jié)果反而變差.

    表2 差分步長(zhǎng)求解結(jié)果對(duì)比Tab.2 Comparison of results with differentδs

    (2)梯度計(jì)算方法

    針對(duì)本文研究的算例管網(wǎng),分別分析了通過(guò)差分法(見(jiàn)式4)以及線性比例法(見(jiàn)式6)進(jìn)行梯度計(jì)算的求解效果.線性比例法計(jì)算后得到的結(jié)果如圖3所示,其中迭代次數(shù)為400次,收斂時(shí)間為1min,梯度向量的模長(zhǎng)為3230.83;通過(guò)差分法編碼計(jì)算后得到的結(jié)果如圖2所示,其中迭代次數(shù)為400次,收斂時(shí)間為5min,梯度向量的模長(zhǎng)為0.25.

    圖3 污染物注入質(zhì)量濃度變化曲線(線性比例法)Fig.3 Injection concentrations of contamination sources with linearly proportional method

    圖3表明在候選污染源節(jié)點(diǎn)115,117,119和120處有污染物注入,所定位出的污染源位置在假定污染源注入點(diǎn)附近區(qū)域,但對(duì)于污染源節(jié)點(diǎn)污染物濃度的變化識(shí)別精度較差;而通過(guò)差分法計(jì)算(見(jiàn)圖2)則能夠直接確定污染源注入的時(shí)間、地點(diǎn)和侵入過(guò)程.此外,迭代次數(shù)相同時(shí),采用線性比例法計(jì)算結(jié)果的梯度向量模長(zhǎng)遠(yuǎn)大于采用差分法計(jì)算的梯度向量模長(zhǎng).因而采用線性比例法計(jì)算的結(jié)果要比用差分步長(zhǎng)編碼計(jì)算的結(jié)果差,但用線性方式編碼計(jì)算的收斂速度快于用差分法計(jì)算的速度.在采用線性比例法進(jìn)行計(jì)算時(shí),還需要為迭代變量設(shè)置一個(gè)可取值區(qū)間,以保證解的非負(fù)性和所能取的最大值,否則將會(huì)出現(xiàn)無(wú)法收斂的情況.而用差分法計(jì)算時(shí),僅需要保證解的非負(fù)性,無(wú)需設(shè)置解的最大限值.分析其原因在于,線性比例法計(jì)算時(shí),由于水質(zhì)監(jiān)測(cè)點(diǎn)對(duì)于某一候選污染源節(jié)點(diǎn)的權(quán)重系數(shù),是取相應(yīng)時(shí)間步長(zhǎng)時(shí)段內(nèi)的權(quán)重系數(shù)平均值,因此所求得的梯度向量并不能精確地給出目標(biāo)函數(shù)的收斂方向,導(dǎo)致計(jì)算結(jié)果較差.因此,建議在需要快速確定污染源點(diǎn)位時(shí)可采用線性比例法,當(dāng)需要更精確地分析污染物注入時(shí)間、濃度及變化模式時(shí),采用基于管網(wǎng)水質(zhì)模擬的差分法進(jìn)行計(jì)算.

    4 結(jié)論

    負(fù)梯度法可有效解決供水管網(wǎng)單點(diǎn)污染物注入模式和多點(diǎn)污染物注入模式下的污染源定位問(wèn)題.相對(duì)于遺傳算法等隨機(jī)搜索算法,沿負(fù)梯度方向搜索具有更高的搜索效率,通過(guò)引入差分算法可有效降低目標(biāo)函數(shù)對(duì)于決策變量求導(dǎo)的難度.而通過(guò)將利用最小覆蓋集確定候選污染源節(jié)點(diǎn)集合的方法與梯度法相結(jié)合,可進(jìn)一步提高對(duì)污染源識(shí)別問(wèn)題的求解能力.此外,通過(guò)不同差分步長(zhǎng)下的計(jì)算結(jié)果表明,當(dāng)其值取水質(zhì)模擬所采用的水質(zhì)公差值時(shí),算法求解結(jié)果最好.而相對(duì)于用線性比例法求解梯度向量,差分法可獲得更高的求解質(zhì)量,不但可以定位污染源,并能以較高精度識(shí)別污染物注入的時(shí)間和濃度變化模式.

    由于污染源識(shí)別問(wèn)題是利用有限個(gè)水質(zhì)監(jiān)測(cè)點(diǎn)獲得的數(shù)據(jù)反演污染源的侵入過(guò)程,因此,水質(zhì)監(jiān)測(cè)點(diǎn)的數(shù)量及空間布置對(duì)于問(wèn)題的求解有直接的影響.因此,在將來(lái)的研究中,本方法應(yīng)進(jìn)一步與水質(zhì)監(jiān)測(cè)點(diǎn)的優(yōu)化布置結(jié)合,以更進(jìn)一步提高算法的求解效率.此外,對(duì)于無(wú)法有效獲取管網(wǎng)中水質(zhì)的精確濃度數(shù)據(jù)的情況,如何對(duì)突發(fā)污染事件進(jìn)行污染源的辨識(shí),是未來(lái)需要進(jìn)一步研究的問(wèn)題.

    [1]Laird C D,Biegler L T,van Bloemen Waanders B G,et al.Contamination source determination for water networks[J].Water Resources Planning and Management,2005,131(2):125.

    [2]Laird C D,Biegler L T,van Bloemen Waanders B G.Mixedinteger approach for obtaining unique solutions in source inversion of water networks[J].Water Resources Planning and Management,2006,132(4):242.

    [3]GUAN Jiabao,Aral M M,Maslia M L,et al.Identification of contaminant sources in water distribution systems using simulation optimization method:case study [J]. Water Resources Planning and Management,2006,132(4):252.

    [4]Pries A,Ostfeld A.Contamination source identification in water systems:a hybrid model trees-linear programming scheme[J].Water Resources Planning and Management,2006,132(4):263.

    [5]Cristo D C,Leopardi A.Pollution source identification of accidental contamination in water distribution networks[J].Water Resources Planning and Management,2008,134(2):197.

    [6]LIU Li,Ranji Ranjithan S,Mahinthakumar G.Contamination source identification in water distribution systems using an adaptive dynamic optimization procedure[J].Water Resources Planning and Management,2011,137(2):183.

    [7]YANG Xueyao,Boccelli D L,De Sanctis A E.A Bayesian approach for probabilistic contaminantion source identification[C]//World Environmental and Water Resource Congress.[S.l.]:ASCE,2011.

    [8]Rossman L A.EPANET2 users’manual[M].Cincinnati:U.S.Environmental Protection Agency,2000.

    [9]王鳳仙.供水管網(wǎng)水質(zhì)管理軟件應(yīng)用和預(yù)警技術(shù)研究[D].上海:同濟(jì)大學(xué)環(huán)境科學(xué)與工程學(xué)院,2010.WANG Fengxian.Study on application of water quality management software and early-warning technology for water distribution systems[D].Shanghai:College of Enviromental Science and Engineering of Tongji University,2010.

    猜你喜歡
    差分法污染源步長(zhǎng)
    二維粘彈性棒和板問(wèn)題ADI有限差分法
    基于Armijo搜索步長(zhǎng)的BFGS與DFP擬牛頓法的比較研究
    持續(xù)推進(jìn)固定污染源排污許可管理全覆蓋
    基于污染源解析的空氣污染治理對(duì)策研究
    十二五”期間佳木斯市污染源排放狀況分析
    看不見(jiàn)的污染源——臭氧
    基于逐維改進(jìn)的自適應(yīng)步長(zhǎng)布谷鳥(niǎo)搜索算法
    基于SQMR方法的三維CSAMT有限差分法數(shù)值模擬
    一種新型光伏系統(tǒng)MPPT變步長(zhǎng)滯環(huán)比較P&O法
    有限差分法模擬電梯懸掛系統(tǒng)橫向受迫振動(dòng)
    久久久久久久久久久久大奶| 久久精品国产亚洲av高清一级| 国产91精品成人一区二区三区 | 欧美亚洲 丝袜 人妻 在线| 又大又黄又爽视频免费| 脱女人内裤的视频| 欧美日韩黄片免| 久久毛片免费看一区二区三区| 亚洲欧美精品自产自拍| 亚洲精品国产一区二区精华液| 中文字幕av电影在线播放| 亚洲国产最新在线播放| 日本欧美国产在线视频| 手机成人av网站| 老司机亚洲免费影院| 一级毛片我不卡| 婷婷色综合www| 下体分泌物呈黄色| 精品少妇一区二区三区视频日本电影| 国产精品99久久99久久久不卡| 国语对白做爰xxxⅹ性视频网站| 好男人视频免费观看在线| av在线播放精品| 99久久综合免费| 亚洲精品日本国产第一区| 日本五十路高清| videosex国产| 最黄视频免费看| 大片电影免费在线观看免费| svipshipincom国产片| 91麻豆av在线| 日韩中文字幕视频在线看片| 国产精品一区二区在线不卡| 男人爽女人下面视频在线观看| 日韩伦理黄色片| 国产不卡av网站在线观看| 亚洲av日韩精品久久久久久密 | av视频免费观看在线观看| 日韩视频在线欧美| 一区福利在线观看| 人妻一区二区av| 国产成人一区二区在线| 亚洲精品美女久久av网站| 日韩熟女老妇一区二区性免费视频| 天天添夜夜摸| 老鸭窝网址在线观看| 麻豆av在线久日| 又大又黄又爽视频免费| 日本wwww免费看| 大片电影免费在线观看免费| 成人影院久久| 脱女人内裤的视频| 欧美日韩国产mv在线观看视频| 91九色精品人成在线观看| 精品福利观看| 狠狠精品人妻久久久久久综合| 国产野战对白在线观看| 亚洲国产看品久久| 一区二区日韩欧美中文字幕| 99re6热这里在线精品视频| 大香蕉久久网| 51午夜福利影视在线观看| 国产成人91sexporn| 丝袜美足系列| 天天影视国产精品| 一级片免费观看大全| 人人妻人人澡人人看| 在线观看人妻少妇| 一级毛片我不卡| 美女主播在线视频| 日韩伦理黄色片| 精品久久蜜臀av无| 久久久久精品国产欧美久久久 | 国产精品一国产av| 久久人人爽av亚洲精品天堂| 成年av动漫网址| 一本一本久久a久久精品综合妖精| 午夜福利在线免费观看网站| 国产免费福利视频在线观看| 电影成人av| 男人操女人黄网站| xxx大片免费视频| 一级毛片我不卡| 欧美av亚洲av综合av国产av| 日韩大码丰满熟妇| 男男h啪啪无遮挡| 99久久人妻综合| 免费高清在线观看日韩| 一边亲一边摸免费视频| 一二三四在线观看免费中文在| 国产日韩欧美在线精品| 另类亚洲欧美激情| 美女脱内裤让男人舔精品视频| 亚洲av欧美aⅴ国产| svipshipincom国产片| 99热全是精品| 亚洲av日韩精品久久久久久密 | 亚洲熟女毛片儿| 中文字幕人妻丝袜一区二区| 亚洲国产av新网站| 黄色视频在线播放观看不卡| 国产欧美日韩精品亚洲av| 久9热在线精品视频| 久久 成人 亚洲| 精品人妻一区二区三区麻豆| 国产女主播在线喷水免费视频网站| 精品国产一区二区久久| 大码成人一级视频| 男女边吃奶边做爰视频| 亚洲欧洲日产国产| 亚洲精品久久午夜乱码| 天天躁夜夜躁狠狠躁躁| 51午夜福利影视在线观看| av网站在线播放免费| 国产精品一区二区精品视频观看| 1024香蕉在线观看| 国产成人91sexporn| 热99国产精品久久久久久7| 国产又爽黄色视频| 国产黄频视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 天天影视国产精品| 大香蕉久久成人网| 国产在视频线精品| 18禁国产床啪视频网站| 后天国语完整版免费观看| xxxhd国产人妻xxx| 国产成人系列免费观看| 如日韩欧美国产精品一区二区三区| av国产精品久久久久影院| 亚洲av在线观看美女高潮| 成年美女黄网站色视频大全免费| 中国国产av一级| 老汉色∧v一级毛片| 亚洲成色77777| 成人国语在线视频| 日韩 亚洲 欧美在线| a级片在线免费高清观看视频| 欧美成狂野欧美在线观看| 久久青草综合色| 自线自在国产av| 狠狠婷婷综合久久久久久88av| 精品久久蜜臀av无| 妹子高潮喷水视频| 久久99一区二区三区| 久久午夜综合久久蜜桃| 韩国精品一区二区三区| 久久免费观看电影| 久久99热这里只频精品6学生| 老汉色∧v一级毛片| 国产三级黄色录像| 丝袜人妻中文字幕| 十分钟在线观看高清视频www| 精品少妇久久久久久888优播| 久久青草综合色| 午夜激情久久久久久久| 女人久久www免费人成看片| 蜜桃在线观看..| 性色av乱码一区二区三区2| 观看av在线不卡| 成人三级做爰电影| 亚洲精品av麻豆狂野| 老司机亚洲免费影院| 国产成人欧美| 精品视频人人做人人爽| 久久人妻熟女aⅴ| 中文字幕色久视频| 丝袜在线中文字幕| 精品免费久久久久久久清纯 | 美女午夜性视频免费| 最新的欧美精品一区二区| 黄色毛片三级朝国网站| 国产又色又爽无遮挡免| 少妇人妻 视频| 国产片特级美女逼逼视频| 亚洲,欧美,日韩| 亚洲av男天堂| 欧美亚洲日本最大视频资源| 亚洲国产欧美一区二区综合| 王馨瑶露胸无遮挡在线观看| 午夜免费观看性视频| 日本午夜av视频| 无限看片的www在线观看| 亚洲视频免费观看视频| 国产成人欧美在线观看 | 国产一区二区在线观看av| 亚洲欧美中文字幕日韩二区| 天天添夜夜摸| 50天的宝宝边吃奶边哭怎么回事| 精品一品国产午夜福利视频| 99香蕉大伊视频| 99国产精品一区二区三区| 一本一本久久a久久精品综合妖精| 日本vs欧美在线观看视频| 国产精品一区二区精品视频观看| 九草在线视频观看| 国产高清国产精品国产三级| 超碰97精品在线观看| 欧美人与善性xxx| 一区二区三区精品91| 在现免费观看毛片| 久久久久精品国产欧美久久久 | 精品一区二区三区四区五区乱码 | 中文字幕制服av| 亚洲av日韩精品久久久久久密 | 午夜福利视频在线观看免费| 精品熟女少妇八av免费久了| 国产精品99久久99久久久不卡| 国产精品.久久久| 国产精品三级大全| 免费黄频网站在线观看国产| 在线观看免费午夜福利视频| 日韩一区二区三区影片| 性色av乱码一区二区三区2| 久久久久久久久免费视频了| tube8黄色片| 亚洲美女黄色视频免费看| 最近手机中文字幕大全| 久久av网站| 性少妇av在线| 美女视频免费永久观看网站| 亚洲国产欧美一区二区综合| 超碰成人久久| 女人久久www免费人成看片| 国产亚洲午夜精品一区二区久久| 亚洲国产精品国产精品| 嫁个100分男人电影在线观看 | 欧美亚洲 丝袜 人妻 在线| 一本久久精品| 欧美黑人欧美精品刺激| 丰满人妻熟妇乱又伦精品不卡| 老熟女久久久| 中文字幕人妻熟女乱码| 国产精品一二三区在线看| 肉色欧美久久久久久久蜜桃| 午夜免费鲁丝| 午夜免费男女啪啪视频观看| 午夜激情久久久久久久| 十八禁高潮呻吟视频| 久久精品亚洲av国产电影网| 女性被躁到高潮视频| 国产三级黄色录像| 高潮久久久久久久久久久不卡| 亚洲欧美色中文字幕在线| 99精品久久久久人妻精品| 97人妻天天添夜夜摸| 精品久久久久久久毛片微露脸 | 一级毛片我不卡| 国产精品二区激情视频| av一本久久久久| 91成人精品电影| 香蕉国产在线看| 丝袜脚勾引网站| 黑丝袜美女国产一区| 中国国产av一级| 亚洲第一av免费看| 国产成人免费观看mmmm| 悠悠久久av| 日韩人妻精品一区2区三区| 欧美日韩国产mv在线观看视频| 亚洲成av片中文字幕在线观看| 熟女少妇亚洲综合色aaa.| 日韩制服丝袜自拍偷拍| 汤姆久久久久久久影院中文字幕| 欧美日韩av久久| 国产男女超爽视频在线观看| 少妇猛男粗大的猛烈进出视频| 热99久久久久精品小说推荐| 多毛熟女@视频| 97在线人人人人妻| 婷婷色综合大香蕉| 欧美激情极品国产一区二区三区| 黑人欧美特级aaaaaa片| 另类亚洲欧美激情| 国产精品欧美亚洲77777| 极品少妇高潮喷水抽搐| 一区在线观看完整版| 九色亚洲精品在线播放| 精品免费久久久久久久清纯 | 色精品久久人妻99蜜桃| 中文乱码字字幕精品一区二区三区| 国产免费现黄频在线看| 日日摸夜夜添夜夜爱| 欧美老熟妇乱子伦牲交| svipshipincom国产片| 天天躁狠狠躁夜夜躁狠狠躁| 少妇人妻久久综合中文| 亚洲av在线观看美女高潮| 成人黄色视频免费在线看| 久久久久国产精品人妻一区二区| 天天操日日干夜夜撸| 精品国产国语对白av| 欧美人与善性xxx| 亚洲精品美女久久久久99蜜臀 | 国产欧美亚洲国产| 精品少妇久久久久久888优播| 国产精品国产av在线观看| 久久久久视频综合| 国产精品麻豆人妻色哟哟久久| 亚洲色图 男人天堂 中文字幕| 一级a爱视频在线免费观看| 精品国产乱码久久久久久男人| 中文字幕av电影在线播放| 欧美黑人欧美精品刺激| 欧美日韩视频高清一区二区三区二| 久久亚洲精品不卡| 成年女人毛片免费观看观看9 | 国产国语露脸激情在线看| 又紧又爽又黄一区二区| 久久久久久久久久久久大奶| 欧美性长视频在线观看| 男女免费视频国产| www日本在线高清视频| 黑丝袜美女国产一区| 国产1区2区3区精品| 1024香蕉在线观看| 青春草亚洲视频在线观看| 精品一品国产午夜福利视频| av网站在线播放免费| 天天躁狠狠躁夜夜躁狠狠躁| 后天国语完整版免费观看| 国产在视频线精品| 美女午夜性视频免费| 日韩一本色道免费dvd| 日日爽夜夜爽网站| 精品亚洲乱码少妇综合久久| 最近中文字幕2019免费版| 国产精品免费大片| 久久精品久久久久久久性| 久久国产精品大桥未久av| 高清不卡的av网站| 美女高潮到喷水免费观看| 黑人欧美特级aaaaaa片| 中文字幕人妻丝袜一区二区| 久久青草综合色| 黄片小视频在线播放| 久久久久精品国产欧美久久久 | 老汉色av国产亚洲站长工具| av在线播放精品| 精品人妻熟女毛片av久久网站| 亚洲欧洲精品一区二区精品久久久| 亚洲国产看品久久| 亚洲欧洲日产国产| 亚洲视频免费观看视频| 国产一区二区 视频在线| 国产一卡二卡三卡精品| 黑人欧美特级aaaaaa片| 亚洲精品美女久久久久99蜜臀 | 一级毛片黄色毛片免费观看视频| 99久久综合免费| 女性生殖器流出的白浆| 午夜福利视频精品| 午夜福利在线免费观看网站| 国语对白做爰xxxⅹ性视频网站| 成在线人永久免费视频| 欧美乱码精品一区二区三区| 少妇裸体淫交视频免费看高清 | 精品国产超薄肉色丝袜足j| 国产精品国产三级专区第一集| 97精品久久久久久久久久精品| 美女大奶头黄色视频| 男人爽女人下面视频在线观看| 日本a在线网址| 高清欧美精品videossex| 丝袜美腿诱惑在线| 中文字幕精品免费在线观看视频| 男女免费视频国产| 国产黄色视频一区二区在线观看| 超色免费av| 国产精品国产三级国产专区5o| 亚洲人成电影观看| 亚洲精品中文字幕在线视频| 日韩伦理黄色片| 男人舔女人的私密视频| e午夜精品久久久久久久| av天堂在线播放| 亚洲av男天堂| 国产xxxxx性猛交| 一边亲一边摸免费视频| 免费看不卡的av| 国产精品一区二区在线不卡| 亚洲精品美女久久久久99蜜臀 | 天堂中文最新版在线下载| 不卡av一区二区三区| 久久99一区二区三区| 曰老女人黄片| 麻豆av在线久日| 国产日韩欧美视频二区| 如日韩欧美国产精品一区二区三区| 蜜桃国产av成人99| 黄色片一级片一级黄色片| www.自偷自拍.com| 亚洲成国产人片在线观看| 国产日韩欧美在线精品| svipshipincom国产片| 国产不卡av网站在线观看| 久久国产精品人妻蜜桃| 岛国毛片在线播放| 欧美日韩国产mv在线观看视频| 久久中文字幕一级| 手机成人av网站| 日韩av在线免费看完整版不卡| 日韩一区二区三区影片| 99香蕉大伊视频| 丝瓜视频免费看黄片| 亚洲国产av影院在线观看| 国产日韩欧美亚洲二区| 高清av免费在线| 国产主播在线观看一区二区 | 性高湖久久久久久久久免费观看| 一本综合久久免费| 亚洲av欧美aⅴ国产| 久久人人爽人人片av| 久久性视频一级片| 18禁观看日本| 爱豆传媒免费全集在线观看| 欧美变态另类bdsm刘玥| 国产成人av激情在线播放| 九草在线视频观看| av电影中文网址| 高清黄色对白视频在线免费看| 日韩 亚洲 欧美在线| 嫁个100分男人电影在线观看 | av又黄又爽大尺度在线免费看| 在线观看免费视频网站a站| 丰满人妻熟妇乱又伦精品不卡| 尾随美女入室| 午夜福利视频精品| 51午夜福利影视在线观看| 午夜激情av网站| 嫁个100分男人电影在线观看 | 国产精品人妻久久久影院| 久久精品熟女亚洲av麻豆精品| 免费观看a级毛片全部| 悠悠久久av| 精品久久久久久电影网| 久久狼人影院| 中文精品一卡2卡3卡4更新| 晚上一个人看的免费电影| 少妇粗大呻吟视频| 国产在线视频一区二区| 国产片内射在线| 精品一区二区三卡| bbb黄色大片| 999精品在线视频| 视频在线观看一区二区三区| 精品久久久精品久久久| 亚洲精品美女久久久久99蜜臀 | av天堂在线播放| 69精品国产乱码久久久| 咕卡用的链子| 亚洲精品美女久久av网站| 中文字幕最新亚洲高清| 国产xxxxx性猛交| 99热网站在线观看| 丝袜在线中文字幕| 男的添女的下面高潮视频| 无限看片的www在线观看| 色网站视频免费| 国语对白做爰xxxⅹ性视频网站| 人成视频在线观看免费观看| 婷婷色av中文字幕| 两人在一起打扑克的视频| 日韩精品免费视频一区二区三区| 色视频在线一区二区三区| 丝袜美腿诱惑在线| 亚洲欧美精品综合一区二区三区| 亚洲天堂av无毛| 黑人巨大精品欧美一区二区蜜桃| 婷婷丁香在线五月| 国产91精品成人一区二区三区 | 国产高清国产精品国产三级| 欧美激情极品国产一区二区三区| 亚洲欧美清纯卡通| 丝袜人妻中文字幕| 成年av动漫网址| 久久久久视频综合| 亚洲精品一二三| 亚洲色图综合在线观看| 国产亚洲精品久久久久5区| 精品福利观看| 麻豆国产av国片精品| 亚洲av在线观看美女高潮| 国产男人的电影天堂91| 亚洲男人天堂网一区| 热re99久久国产66热| 亚洲av日韩在线播放| 亚洲熟女毛片儿| 国产精品秋霞免费鲁丝片| 日韩电影二区| xxxhd国产人妻xxx| 日本a在线网址| 天堂俺去俺来也www色官网| 国产成人一区二区三区免费视频网站 | 国产精品九九99| 99国产精品一区二区三区| 只有这里有精品99| 精品亚洲成a人片在线观看| 国产精品欧美亚洲77777| 19禁男女啪啪无遮挡网站| 欧美人与性动交α欧美软件| 超色免费av| av天堂在线播放| 一级毛片电影观看| 国产免费现黄频在线看| 国产精品99久久99久久久不卡| 观看av在线不卡| 在线观看一区二区三区激情| 亚洲情色 制服丝袜| 亚洲成色77777| 黄色视频在线播放观看不卡| 日韩伦理黄色片| 男人添女人高潮全过程视频| 国产99久久九九免费精品| 男女之事视频高清在线观看 | 亚洲国产欧美一区二区综合| 天天躁夜夜躁狠狠躁躁| √禁漫天堂资源中文www| 国产一区二区激情短视频 | 人妻 亚洲 视频| 久久国产亚洲av麻豆专区| 欧美在线一区亚洲| 在线亚洲精品国产二区图片欧美| 亚洲精品国产一区二区精华液| 久久久久精品人妻al黑| 五月天丁香电影| 波多野结衣一区麻豆| 免费少妇av软件| 久久99热这里只频精品6学生| 欧美xxⅹ黑人| 亚洲国产欧美日韩在线播放| 大码成人一级视频| 日韩,欧美,国产一区二区三区| 夜夜骑夜夜射夜夜干| 成人手机av| 丁香六月天网| 91老司机精品| 成年美女黄网站色视频大全免费| 蜜桃国产av成人99| 777米奇影视久久| 欧美日本中文国产一区发布| 色婷婷久久久亚洲欧美| 亚洲天堂av无毛| 日韩 欧美 亚洲 中文字幕| av在线播放精品| 交换朋友夫妻互换小说| 成人国产一区最新在线观看 | 老汉色av国产亚洲站长工具| 中文乱码字字幕精品一区二区三区| 欧美人与性动交α欧美精品济南到| 日韩 欧美 亚洲 中文字幕| 成人影院久久| 亚洲久久久国产精品| 脱女人内裤的视频| 国产精品一区二区在线观看99| 亚洲自偷自拍图片 自拍| 精品福利永久在线观看| 97精品久久久久久久久久精品| 嫁个100分男人电影在线观看 | 免费av中文字幕在线| 国产免费现黄频在线看| 麻豆乱淫一区二区| 久久久久久亚洲精品国产蜜桃av| 午夜免费男女啪啪视频观看| 国语对白做爰xxxⅹ性视频网站| 中文欧美无线码| 亚洲精品成人av观看孕妇| 久久精品国产亚洲av高清一级| 亚洲欧美清纯卡通| 国产精品一区二区在线观看99| 丝袜美足系列| 在线看a的网站| 亚洲欧美精品综合一区二区三区| 久久性视频一级片| 日韩制服丝袜自拍偷拍| 高清av免费在线| 一本一本久久a久久精品综合妖精| 中文字幕高清在线视频| 色网站视频免费| 又黄又粗又硬又大视频| 麻豆av在线久日| 飞空精品影院首页| 成人国语在线视频| 欧美黑人精品巨大| 久久人人爽av亚洲精品天堂| 亚洲成色77777| 亚洲精品国产av成人精品| 国产成人精品久久久久久| 麻豆乱淫一区二区| 欧美日韩一级在线毛片| 久久久久久久国产电影| 国产成人欧美在线观看 | 亚洲精品自拍成人| 日韩中文字幕欧美一区二区 | 每晚都被弄得嗷嗷叫到高潮| 国产麻豆69| 一个人免费看片子| 免费av中文字幕在线| 国产成人一区二区在线| 午夜91福利影院| 男女之事视频高清在线观看 | 成年人免费黄色播放视频| 国产免费视频播放在线视频| 丝袜美腿诱惑在线| 日本欧美国产在线视频| 免费人妻精品一区二区三区视频| 亚洲图色成人| 一边亲一边摸免费视频| 十八禁网站网址无遮挡| 热re99久久国产66热| 亚洲黑人精品在线| 蜜桃国产av成人99| 亚洲精品乱久久久久久| 两个人免费观看高清视频| 色视频在线一区二区三区|