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

    基于遺傳算法的供水管網(wǎng)需求乘數(shù)因子校核

    2019-10-11 11:24:36毛潤(rùn)康杜坤周明陳攀雷雨晴丁榕藝楊佳莉
    軟件導(dǎo)刊 2019年7期
    關(guān)鍵詞:校核遺傳算法

    毛潤(rùn)康 杜坤 周明 陳攀 雷雨晴 丁榕藝 楊佳莉

    基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(51608242);云南省應(yīng)用基礎(chǔ)研究青年項(xiàng)目(2017FD094)

    作者簡(jiǎn)介:毛潤(rùn)康(1995-),男,昆明理工大學(xué)建筑工程學(xué)院碩士研究生,研究方向?yàn)槭姓こ?杜坤(1986-),男,博士,昆明理工大學(xué)建筑工程學(xué)院講師,研究方向?yàn)槭姓こ獭1疚耐ㄓ嵶髡撸憾爬ぁ?/p>

    摘 要:水力模型已廣泛應(yīng)用于給水管網(wǎng)設(shè)計(jì)、分析與運(yùn)行中。在所有水力模型中,需水量是導(dǎo)致模型輸出最不確定的參數(shù)之一。因用水情況不確定,使得管網(wǎng)中的節(jié)點(diǎn)需水量變得異常復(fù)雜。在大多數(shù)實(shí)際管網(wǎng)中,用于校核節(jié)點(diǎn)需水量的監(jiān)測(cè)設(shè)備數(shù)量有限,且小于未知量個(gè)數(shù),使得節(jié)點(diǎn)需水量校核作為欠定問題,令節(jié)點(diǎn)需水量校準(zhǔn)產(chǎn)生較大誤差,并且傳統(tǒng)遺傳算法校核節(jié)點(diǎn)需水量的方法是假定所有節(jié)點(diǎn)的需求乘數(shù)因子一致,這也導(dǎo)致校核后的模型無法接近真實(shí)運(yùn)行情況,因此提出在欠定條件下用遺傳算法解決需求乘數(shù)因子的校核問題。通過對(duì)一個(gè)實(shí)際案例多次運(yùn)行并取平均值作為結(jié)果進(jìn)行驗(yàn)證,結(jié)果表明,遺傳算法的校核結(jié)果不僅能夠與被測(cè)位置的實(shí)際值相擬合,而且可以得到非測(cè)量位置的管道流量和節(jié)點(diǎn)水頭,其中校核后的節(jié)點(diǎn)水頭和管道流量誤差較小,平均誤差分別為1.78%、4.05%。該方法相比于傳統(tǒng)校核方法具有更高精度,且更能反映出管網(wǎng)真實(shí)運(yùn)行情況,同時(shí)還避免了傳統(tǒng)校核方法中因遺傳算法產(chǎn)生局部最優(yōu)解而導(dǎo)致誤差偏大的問題,對(duì)于大型管網(wǎng)模型校核也具有一定參考價(jià)值。

    關(guān)鍵詞:遺傳算法;校核;需求乘數(shù)因子校核;欠定問題

    DOI:10. 11907/rjdk. 191608 開放科學(xué)(資源服務(wù))標(biāo)識(shí)碼(OSID):

    中圖分類號(hào):TP301文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1672-7800(2019)007-0060-05

    Calibration of Water Demand Multipliers in Water

    Distribution Systems Using Genetic Algorithms

    MAO Run-kang, DU Kun, ZHOU Ming, CHEN Pan, LEI Yu-qing, DING Rong-yi, YANG Jia-li

    (Faculty of Civil Engineering and Mechanics, Kunming University of Science and Technology, Kunming 650500, China)

    Abstract:Hydraulic models have been widely used for design, analysis, and operation of water distribution systems. As with all hydraulic models, water demands are one of the main parameters that cause the most uncertainty to the model outputs. However, estimation of the demand parameters is usually complicated due to the stochastic behavior of the water consumptions. This is an underdetermined system where the number of measurements is limited and less than the number of unknows,so the calibration of water demand will produce a large error,and the traditional genetic algorithm to calibrate the water demand of nodes is to assume that the demand multiplier factors of all nodes are the same, which leads to the problem that the calibration of model can not be close to the real operation. This paper presents an approach to calibration of the demand multiplier factors under an ill-posed condition where the number of measurements is less than the number of parameter variables. The problem is solved by using a genetic algorithm (GA). A practical case was run multiple times was taken and the average to validate that not only is the GA able to match the calibrated values at measured locations, but by using multiple runs of the GA model, the flow rates and nodal heads at nonmeasured locations can be estimated. The error of node head and pipeline flow is small after check, and their average error is 1.78% and 4.05%, respectively. Compared with the traditional optimized method, this method has higher accuracy and can reflect the real operation of the pipe network. At the same time, it also avoids the problem of large error caused by the local optimal solution produced by genetic algorithm in the traditional optimized method. It also has a certain reference for the calibration of large pipe network model.

    Key Words:genetic algorithms; calibration; calibration of water demand multipliers; underdetermined problem

    0 引言

    供水系統(tǒng)(WDS)作為城市基礎(chǔ)設(shè)施不可缺少的組成部分,每天都要容納大量的水量,以確保提供穩(wěn)定用水。供水管網(wǎng)隨著人口增長(zhǎng)而變得更加復(fù)雜,因此需要加強(qiáng)管理。在現(xiàn)有管網(wǎng)中確定管段流量和節(jié)點(diǎn)水頭是優(yōu)化WDS的必要步驟,該任務(wù)可以通過使用傳感器等測(cè)量設(shè)備完成。然而,傳感器只能捕獲系統(tǒng)中某些節(jié)點(diǎn)位置的狀態(tài)。因此,如何利用有限個(gè)數(shù)的設(shè)備校核整個(gè)WDS模型需要作進(jìn)一步研究。

    WDS模型校核是指對(duì)節(jié)點(diǎn)需水量和管段粗糙系數(shù)取值進(jìn)行調(diào)整,使模型預(yù)測(cè)節(jié)點(diǎn)水壓和管段流量的模擬數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)在一定誤差精度范圍內(nèi)相匹配的過程。Ormsbee&Lingireddy很好地解決了管網(wǎng)校核問題[1],其提出一個(gè)七步校準(zhǔn)過程,包括:①模型預(yù)定用途確定;②模型參數(shù)確定;③校核數(shù)據(jù)收集;④模型結(jié)果評(píng)價(jià);⑤模型宏觀校正;⑥敏感性分析;⑦模型微觀校正。

    Savic等[2]對(duì)管網(wǎng)校核進(jìn)行了完整的文獻(xiàn)綜述,管網(wǎng)校核可以作如下分類:狀態(tài)(瞬態(tài)/靜態(tài))、計(jì)算方法(迭代/顯式/隱式)或優(yōu)化方法(傳統(tǒng)/改進(jìn))。

    針對(duì)管網(wǎng)節(jié)點(diǎn)流量校核,吳學(xué)偉等[3]將節(jié)點(diǎn)水壓作為已知參數(shù)對(duì)節(jié)點(diǎn)流量進(jìn)行計(jì)算,但由于驗(yàn)證的管網(wǎng)是實(shí)驗(yàn)室的小型管網(wǎng),因此對(duì)大型管網(wǎng)還需作進(jìn)一步研究;王榮和[4]提出可以用遺傳算法校核模型的節(jié)點(diǎn)流量與管段粗糙系數(shù);陳宇輝[5]研究了遺傳算法、模擬退火算法、自適應(yīng)遺傳算法在求解目標(biāo)函數(shù)時(shí)的應(yīng)用效果;鄭飛飛等[6]提出一種基于非交叉抖動(dòng)爬行突變的遺傳算法(CMBGA)校核管網(wǎng),該算法沒有交叉操作,只有變異操作,并通過兩個(gè)案例得到驗(yàn)證,因而相比傳統(tǒng)算法在結(jié)果上得到了很大改善,且需要調(diào)整的參數(shù)較少;劉佳明等[7]針對(duì)標(biāo)準(zhǔn)遺傳算法收斂速度慢,且容易陷入局部最優(yōu)解的問題,對(duì)遺傳算法進(jìn)行改進(jìn),提高了算法速度與精度;Yu等[8]提出一種基于貪婪算法的抽樣設(shè)計(jì)法校核水力模型,該方法以現(xiàn)有傳感器為基礎(chǔ),在每個(gè)優(yōu)化仿真步驟中依次添加一個(gè)新傳感器,從而盡可能使校核中的不確定性最小化,并取得了較好的效率和精度。

    管網(wǎng)模型的不確定性是由多種因素引起的,如管道粗糙度、節(jié)點(diǎn)需水量或閥門狀態(tài)等。由于管道粗糙系數(shù)和節(jié)點(diǎn)需水量對(duì)管網(wǎng)不確定性的影響較大,因此其通常用于靜態(tài)校核。然而,粗糙度系數(shù)通常只在長(zhǎng)期內(nèi)變化。例如,Haddad等[9]考慮了該參數(shù)的年度變化。因此,節(jié)點(diǎn)需水量是短期內(nèi)影響管網(wǎng)模型不確定性的主要參數(shù)(例如每小時(shí)、每天)。

    在許多文獻(xiàn)中,研究人員采用不同方法對(duì)需水量估算進(jìn)行研究,例如擴(kuò)展卡爾曼濾波、跟蹤狀態(tài)校核和卡爾曼濾波、遺傳算法和粒子濾波。但這些模型是建立在給定框架基礎(chǔ)上的,在該框架中,測(cè)量位置是預(yù)先確定的,校準(zhǔn)數(shù)量小于測(cè)量點(diǎn)數(shù)量。只有少數(shù)文獻(xiàn)直接討論了過欠定系統(tǒng),如比例需求法與奇異值分解(SVD)法。

    在以往研究中,當(dāng)測(cè)量數(shù)量小于校核參數(shù)變量數(shù)量時(shí),對(duì)供水系統(tǒng)節(jié)點(diǎn)需水量的校核是一個(gè)非線性欠定問題。通過QR分解、SVD、牛頓—拉夫森方法中的Moore-Penrose偽逆矩陣等局部線性化方法,可以得到問題局部解。但是由于解的不唯一性,計(jì)算結(jié)果不僅離實(shí)際解差距很遠(yuǎn),甚至在某些節(jié)點(diǎn)上會(huì)產(chǎn)生負(fù)值。顯然,如果測(cè)量設(shè)備數(shù)據(jù)被認(rèn)為是唯一的已知輸入,那么尋找可靠的需水量校準(zhǔn)方法仍然是水力研究人員面臨的挑戰(zhàn)。

    鑒于此,本文提出一種校核欠定系統(tǒng)需水量乘數(shù)因子的方法,即測(cè)量數(shù)量小于需求參數(shù)變量數(shù)量。EPANET工具包用于求解管網(wǎng)方程組,遺傳算法(GAS)用于尋找已知測(cè)量輸入與其模擬值之間的最優(yōu)解。將遺傳算法多次運(yùn)行后的平均值作為校核后的最優(yōu)解,通過對(duì)管網(wǎng)不同情景進(jìn)行測(cè)試以對(duì)管網(wǎng)可靠性進(jìn)行評(píng)估。此外,本研究還探討了使用Piller抽樣設(shè)計(jì)技術(shù)選擇最佳測(cè)量位置,以提高校準(zhǔn)模型的精確性[10],然后將該方法應(yīng)用于實(shí)際案例研究,最后得出結(jié)論并給出相關(guān)建議。

    1 GA校準(zhǔn)模型

    該模型將隱式技術(shù)應(yīng)用于穩(wěn)態(tài)水力模擬,將校核過程描述為一個(gè)優(yōu)化問題。目標(biāo)函數(shù)是使模型中模擬值與其相應(yīng)實(shí)測(cè)值之間的差值最小化。決策變量是“決策變量”章節(jié)描述的節(jié)點(diǎn)需求乘數(shù)因子(DMFs)。

    2 目標(biāo)函數(shù)

    本文采用最小二乘法求解目標(biāo)函數(shù),該方法是將測(cè)量位置上管道流量和節(jié)點(diǎn)水頭的測(cè)量值與模擬值之間的平方差之和最小化[11]。

    [MinF=i=1nwHH0i-Hi2+j=1mwqQ0j-Qj2]? ? (1)

    式中,n為水壓監(jiān)測(cè)點(diǎn)數(shù),i為第i個(gè)水壓監(jiān)測(cè)點(diǎn),wH、wq為水壓與流量值的權(quán)重系數(shù),[H0i、Hi]為第i點(diǎn)的壓力監(jiān)測(cè)值與模型值,[Q0j、Qj]為第j管道的流量監(jiān)測(cè)值與模型值。

    管道流量、節(jié)點(diǎn)水頭和節(jié)點(diǎn)需水量均與時(shí)間有關(guān)。然而,在短時(shí)間內(nèi),如30min或1h,可以被認(rèn)為是恒定的,因此在公式(1)中沒有明確給出時(shí)間依賴符號(hào)。測(cè)量值是從現(xiàn)場(chǎng)測(cè)量設(shè)備中獲得的,或者可以通過運(yùn)行諸如EPANET等水力模擬工具包生成這些值。加權(quán)因子wH、wq可以用監(jiān)測(cè)值的平方倒數(shù)進(jìn)行計(jì)算[12]。

    3 決策變量

    在WDS模型中,每個(gè)節(jié)點(diǎn)的需水量是通過將每個(gè)時(shí)間t上的基本需水量與其相應(yīng)的DMF相乘進(jìn)行計(jì)算的。

    [Dk,t=D0,k×fk,t]? ? ? ? ? ? ?(2)

    D0,k為第K節(jié)點(diǎn)基本需水量,使用季度/年度用水帳單信息,fk,t為第K節(jié)點(diǎn)在T時(shí)刻的需水量乘數(shù)因子。

    優(yōu)化問題的決策變量是在每個(gè)時(shí)刻節(jié)點(diǎn)需水量的需求乘數(shù)因子fk,t,其范圍如下:

    [fminkfk,tfmaxk]? ? ? ? ? ? (3)

    其中[fmink]必須大于零,[fmaxk]是根據(jù)峰值需求因素的典型值進(jìn)行選擇的[13],本文中[fmaxk]取1.5。

    表1給出了一次求解6個(gè)需求乘數(shù)因子的GA(即染色體)案例。它們是從0~1.50的范圍內(nèi)選取的,增量步長(zhǎng)為0.02(值隨機(jī)產(chǎn)生)。每個(gè)需求乘數(shù)因子有一個(gè)浮點(diǎn)數(shù)編碼,從0~75不等。通過使用這些編碼信息,將遺傳算法產(chǎn)生的染色體解碼成一組需求乘數(shù)因子,并乘以基本需求得到的實(shí)際需水量,即可用于水力模擬過程。

    表1 需求校核問題的GA染色體及解碼實(shí)例

    4 校核精度要求

    校核后的模型必須達(dá)到一定精度,模型精度是由實(shí)際測(cè)量值與模擬值差值進(jìn)行衡量的。評(píng)價(jià)管網(wǎng)模型是否符合實(shí)際,國(guó)內(nèi)和國(guó)際上都沒有出臺(tái)相應(yīng)技術(shù)標(biāo)準(zhǔn)或行業(yè)標(biāo)準(zhǔn)。但是一些研究機(jī)構(gòu)根據(jù)多年的科研和工程經(jīng)驗(yàn),給出了相應(yīng)建議值,僅供參考。英國(guó)水研究中心(WRC)校核標(biāo)準(zhǔn)包括以下幾部分[14]:

    (1)流量監(jiān)測(cè)點(diǎn):當(dāng)主干管流量大于總用水量的10%時(shí),誤差取測(cè)量值的±5%;否則,誤差取測(cè)量值的±10%。

    (2)壓力監(jiān)測(cè)點(diǎn):85%的監(jiān)測(cè)點(diǎn)壓力偏差在±0.5m 或? ? ±5%的整個(gè)系統(tǒng)最大水頭損失;95%的監(jiān)測(cè)點(diǎn)壓力偏差在 ±0.75m 或±7.5%的整個(gè)系統(tǒng)最大水頭損失;100%的監(jiān)測(cè)點(diǎn)壓力偏差在±2m或±15%的整個(gè)系統(tǒng)最大水頭損失。

    (3)分界線:模擬計(jì)算得到的管網(wǎng)壓力分界線應(yīng)與實(shí)際情況相吻合。

    (4)供水趨勢(shì):模擬計(jì)算得到的供水趨勢(shì)應(yīng)與實(shí)際情況相吻合。

    (5)壓力分布:模擬計(jì)算得到的各節(jié)點(diǎn)水壓分布情況應(yīng)與實(shí)際情況相吻合,計(jì)算得到的高壓區(qū)和低壓區(qū)等應(yīng)與實(shí)際情況相吻合。

    根據(jù)國(guó)內(nèi)供水管網(wǎng)專家趙洪賓教授[15]提出的符合中國(guó)實(shí)際情況的管網(wǎng)水力模型校核精度建議值:①計(jì)算出的各水源出廠供水量、供水壓力與實(shí)測(cè)記錄相吻合;②計(jì)算出的各壓力監(jiān)測(cè)點(diǎn)水壓與實(shí)測(cè)記錄吻合程度:100%監(jiān)測(cè)點(diǎn)的壓力實(shí)測(cè)值與計(jì)算值之差≤±40 kPa,80%監(jiān)測(cè)點(diǎn)的壓力實(shí)測(cè)值與計(jì)算值之差≤±20 kPa,50%監(jiān)測(cè)點(diǎn)的壓力實(shí)測(cè)值與計(jì)算值之差≤±10 kPa;③對(duì)于管段流量占管網(wǎng)總供水量1%以上的管段,誤差<±5%,對(duì)于管段流量占管網(wǎng)總供水量0.5%以上的管段,誤差<±10%;④計(jì)算得到的各節(jié)點(diǎn)水壓分布情況應(yīng)與實(shí)際情況相吻合,計(jì)算得到的低壓區(qū)應(yīng)與實(shí)際情況相吻合。

    5 GA過程與算子

    本研究實(shí)現(xiàn)的遺傳算法校核是用MATLAB編寫的,該算法流程如圖2所示。隨機(jī)生成染色體的初始群體,并將其解碼為對(duì)應(yīng)每個(gè)節(jié)點(diǎn)的DMF值,給管網(wǎng)每個(gè)節(jié)點(diǎn)分配DMF中的一個(gè)值,然后調(diào)用EPANET模擬管網(wǎng)的穩(wěn)態(tài)水力學(xué)。通過計(jì)算目標(biāo)函數(shù),得到校核后測(cè)量點(diǎn)的管道流量和水頭(Qj,Hi),并與實(shí)際測(cè)量值進(jìn)行比較。應(yīng)用目標(biāo)函數(shù)的逆(1/(F+1)表達(dá)式,以避免F為零時(shí)的不確定形式)定義GA每個(gè)成員的適應(yīng)度函數(shù)。這是衡量每個(gè)成員質(zhì)量的標(biāo)準(zhǔn),用來決定每個(gè)成員的生存機(jī)會(huì)。

    通過應(yīng)用遺傳算法的選擇、交叉與變異,產(chǎn)生繼承前幾代特征的新一代,然后反復(fù)校核過程,直到達(dá)到標(biāo)準(zhǔn)后停止。

    對(duì)于選擇算子,Goldberg &Deb [16]的一項(xiàng)研究建議使用錦標(biāo)賽選擇算子,因?yàn)槠渚哂懈玫氖諗啃?。此外,Goldberg還提出概率相對(duì)較高的兩點(diǎn)交叉算子(Pc=0.6~1.0)與突變概率為Pm≈1/編碼長(zhǎng)度[17]的位變異算子,以提高GA模型性能。因此,在本文研究中將競(jìng)賽選擇、兩點(diǎn)交叉、位變異等方法應(yīng)用于遺傳算法模型校核。

    圖1 需求乘數(shù)因子遺傳算法校核流程

    6 測(cè)量地點(diǎn)選擇

    校準(zhǔn)模型的準(zhǔn)確性不僅取決于現(xiàn)場(chǎng)測(cè)量數(shù)量,還取決于測(cè)量位置。Nhu&Angu等[18]使用遺傳算法對(duì)管網(wǎng)進(jìn)行校核,發(fā)現(xiàn)在相同管網(wǎng)上如果測(cè)量點(diǎn)位置不同,得到的結(jié)果也有所不同。許多研究人員利用各種數(shù)學(xué)和統(tǒng)計(jì)方法研究了最優(yōu)測(cè)量位置問題,例如Yu、Powell、Vítkovsk、Berry、Propato以及Krause等在校準(zhǔn)模型中采用抽樣設(shè)計(jì)(SD)方法選擇測(cè)量地點(diǎn)。Piller使用SD方法將狀態(tài)向量估計(jì)中測(cè)量誤差的影響降到最低;Bush&Uber提出3種基于D-最優(yōu)性準(zhǔn)則的SD方法[19]:最大和法、加權(quán)和法和最大最小法,這3種方法都是以雅可比矩陣為基礎(chǔ)的;Kapelan等[20]發(fā)現(xiàn)在SD模型中使用遺傳算法可以得到最優(yōu)的壓力測(cè)量位置,其中第一個(gè)模型被描述為一個(gè)單目標(biāo),第二個(gè)模型被描述為一個(gè)多目標(biāo)優(yōu)化問題。

    在本文中,Pliier提出的基于貪婪算法的SD方法被運(yùn)用于校核需水量乘數(shù)因子模型中,用于尋找最佳測(cè)量位置。測(cè)量位置對(duì)校核結(jié)果的影響通過評(píng)價(jià)遺傳算法模型收斂性確定。

    7 SD方法

    管網(wǎng)中的水力穩(wěn)態(tài)能夠解決節(jié)點(diǎn)連續(xù)性方程與管道能量方程的非線性問題,節(jié)點(diǎn)水頭和管道流量敏感度對(duì)節(jié)點(diǎn)需水量的影響可表示為:

    [AT?q?θ+GD=0ns,ndD?q?θ-A?h?θ=0np,nd]? ? ?(4)

    其中,A、AT為未知節(jié)點(diǎn)水頭的影響矩陣及其轉(zhuǎn)置矩陣,其能反映出節(jié)點(diǎn)及管道情況;q、h為管道流量與節(jié)點(diǎn)水頭向量;GD為分配后的節(jié)點(diǎn)需求矩陣。在本文中,節(jié)點(diǎn)需水量的聚合不用考慮。因此,GD為節(jié)點(diǎn)基本需水量的對(duì)角矩陣,D為管道信息的對(duì)角矩陣,其中對(duì)角元素是各管道水頭損失方程的導(dǎo)數(shù),公式(4)的結(jié)果是管道流量和節(jié)點(diǎn)水壓與節(jié)點(diǎn)需水量相關(guān)的雅克比矩陣。

    [?q?θ=-D-1AATD-1A-1GD?h?θ=-ATD-1A-1GD]? (5)

    由于未知需求參數(shù)f0需要估計(jì),Piller提出的基于貪婪算法的SD方法能使測(cè)量誤差對(duì)狀態(tài)向量估計(jì)的影響最小化,并迭代地選擇測(cè)量位置(S矩陣)。測(cè)量位置的選擇矩陣S已被確定,因此ST0矩陣是滿秩矩陣,而且其偽逆矩陣[(ST0)*∞]的值中具有最小值。矩陣T0=E1J0E2表示雅克比矩陣的平衡矩陣,在公式中E1表示左乘對(duì)角矩陣以表示節(jié)點(diǎn)與管道測(cè)量精度,E2表示右乘與參數(shù)變量變化對(duì)應(yīng)的矩陣,J0= J(yf0)表示公式(5)中在f0處計(jì)算的雅克比矩陣。本文中由于只有需水量參數(shù)需要校核,因此E1 表示單位矩陣,即所有測(cè)量都具有相同精度,E2如公式(6)所示。

    [E2=diag(f0)]? ? ? ? ? (6)

    當(dāng)節(jié)點(diǎn)或管道被選擇作為測(cè)量位置時(shí),其值為1,否則為0,文獻(xiàn)[7]中對(duì)該方法進(jìn)行了詳細(xì)解釋。

    8 案例

    用該案例測(cè)試上述方法在更大管網(wǎng)中的性能,該管網(wǎng)水力模型由EPANET提供,包括26個(gè)管道、27個(gè)節(jié)點(diǎn)、2個(gè)儲(chǔ)水池及3個(gè)泵站。假設(shè)管網(wǎng)中最多可安裝3種測(cè)量設(shè)備,此外假定泵的流量已知,管網(wǎng)如圖4所示。

    圖2 案例管網(wǎng)

    基于貪婪算法的測(cè)量位置選擇模型發(fā)現(xiàn),對(duì)于這種管網(wǎng),管道中的管道流量比節(jié)點(diǎn)水頭更加敏感,則在管道5、12和22處設(shè)置3個(gè)流量測(cè)量點(diǎn),如圖1所示。該遺傳算法的校核模型具有以下特點(diǎn):①?zèng)Q策變量為22個(gè),對(duì)應(yīng)于管網(wǎng)的22個(gè)需求乘數(shù)因子;②決策變量選擇范圍為0~1.5;③選取群體大小N = 500、交叉概率Pc = 0.8、變異概率Pm = 0.0 4和代數(shù)大小N = 1 000作為遺傳算法參數(shù);④為評(píng)價(jià)決策變量增量的影響,分別考察了[Δθ]=0.000 5、0.005、0.02和0.1的不同增量。

    在Intel Corei5(2.9GHz)的計(jì)算機(jī)上運(yùn)行100次GA模型,總運(yùn)行時(shí)間約為16.3h,每一次運(yùn)行約10min。表2總結(jié)了不同增長(zhǎng)步長(zhǎng)對(duì)模型的影響。對(duì)于決策變量4個(gè)不同的增長(zhǎng)步長(zhǎng),誤差大致相等。由此也可得出一個(gè)結(jié)論,即遺傳算法校核模型對(duì)決策變量增長(zhǎng)步長(zhǎng)不敏感。

    表2 不同增長(zhǎng)步長(zhǎng)引起的誤差

    圖3-圖5給出了在決策變量增長(zhǎng)步長(zhǎng)為0.02時(shí)的結(jié)果,分別為節(jié)點(diǎn)水頭、節(jié)點(diǎn)需水量、管道流量。圖的下部表示水壓、需水量、流量真實(shí)值和模擬值,圖的上部表示真實(shí)值與模擬值誤差,另外圖中部分節(jié)點(diǎn)編號(hào)被省略。

    圖3 平均校準(zhǔn)需水量(運(yùn)行GA? ? ?圖4 平均校準(zhǔn)節(jié)點(diǎn)水壓(運(yùn)行GA100次)與管網(wǎng)實(shí)際需水量比較? ? ? ? ?100次)與管網(wǎng)實(shí)際水壓比較

    將節(jié)點(diǎn)需水量分成3組,分別為:0~30L/s、30L/s~60L/s、60L/s以上。圖3顯示出:在需水量較小的第一組中,校核后的誤差很大,為0.67%~61.33%;在第二組中,誤差為14.33%~38.89%;第三組中,誤差為0.53%~9.69%。

    如圖4所示,實(shí)際水頭大的誤差反而小,校核前后各節(jié)點(diǎn)水頭誤差均小于2.5%。其中在節(jié)點(diǎn)15處的實(shí)際水頭為44.66,校核后的水頭為45.69,誤差為2.31%。在節(jié)點(diǎn)2處的實(shí)際水頭為115.84,校核后的水頭為116.83,誤差為0.85%,平均誤差為1.78。

    相應(yīng)地,由圖5可以看出,校核后高流量的管段誤差比低流量的誤差小,最大誤差出現(xiàn)在編號(hào)為24的管道,實(shí)際流量為8L/s,校核流量為10.4%,對(duì)應(yīng)誤差為30%。管網(wǎng)中管道流量校核后的平均誤差為4.05%。其中管段5、12、22為測(cè)量點(diǎn)(已知)。

    圖5 平均校準(zhǔn)管道流量(運(yùn)行GA 100次)與管網(wǎng)實(shí)際流量比較

    一般通過比較監(jiān)測(cè)點(diǎn)的實(shí)測(cè)值與模型計(jì)算值是否吻合衡量模型校核準(zhǔn)確度。由圖4、圖5可以看出,管道流量和節(jié)點(diǎn)水頭誤差均滿足誤差要求,與實(shí)測(cè)值較吻合,說明校核后節(jié)點(diǎn)的需求乘數(shù)因子更符合實(shí)際情況。

    9 結(jié)語

    傳統(tǒng)遺傳算法校核節(jié)點(diǎn)需求量會(huì)出現(xiàn)局部最優(yōu)解,從而導(dǎo)致模擬值與實(shí)際值出現(xiàn)較大偏差,而且由于該傳統(tǒng)方法是假定節(jié)點(diǎn)需求乘數(shù)因子一致,從而導(dǎo)致校核結(jié)果不能反映實(shí)際情況。為解決上述問題,本文提出用遺傳算法校核需求乘數(shù)因子,將遺傳算法運(yùn)行多次并取平均值作為校核結(jié)果,避免了單次運(yùn)行出現(xiàn)局部最優(yōu)解而影響結(jié)果,并通過校核需求乘數(shù)因子使管網(wǎng)中的節(jié)點(diǎn)情況更接近真實(shí)情況,使得校核結(jié)果更加真實(shí)可靠。

    本方法雖在案例中取得了較好效果,為了更好地加以運(yùn)用,未來研究工作還需要尋找更先進(jìn)的方法校核節(jié)點(diǎn)需水量,并減少計(jì)算時(shí)間。此外,由于測(cè)量數(shù)據(jù)中存在誤差而導(dǎo)致校核模型不準(zhǔn)確的問題也需要重視。最后,解決管網(wǎng)漏損問題對(duì)于校核結(jié)果的可靠性也十分重要。

    參考文獻(xiàn):

    [1] BOULOS P F, ORMSBEE L E. Explicit network calibration for multiple loading conditions[J]. Civil Engineering Systems, 1991, 8(3):8.

    [2] KAPELAN Z S, SAVIC D A, WALTERS G A. Calibration of water distribution hydraulic models using a Bayesian-type procedure[J]. J. Hydraul. Eng.,2007,133(8): 927-936,

    [3] 吳學(xué)偉,趙洪賓. 給水管網(wǎng)狀態(tài)估計(jì)方法的研究[J]. 哈爾濱建筑大學(xué)學(xué)報(bào),1995,28(6):60-64.

    [4] 范江,杜坤. 基于加權(quán)最小二乘法的供水管網(wǎng)節(jié)點(diǎn)流量校核[J]. 土木建筑與環(huán)境工程,2016,38(3):73-79.

    [5] 王榮和,姚仁忠,潘建華. 遺傳算法在給水管網(wǎng)現(xiàn)狀分析中的應(yīng)用[J]. 給水排水,2000,26(9):31-36.

    [6] ZHENG F,ZECCHIN A C,SIMPSON A R,et al. Noncrossover dither creeping mutation-based genetic algorithm for pipe network optimization[J]. Journal of Water Resources Planning and Management, 2014,140(4):553-557.

    [7] 劉佳明,王寧,徐浩廣. 面向供水管網(wǎng)水力模型自動(dòng)校核問題的改進(jìn)遺傳算法[J]. 計(jì)算機(jī)系統(tǒng)應(yīng)用,2017,26(12):104-109.

    [8] YU S, SHIPENG C,TUQIAO Z,et al. A greedy sampling design algorithm for the modal calibration of nodal demand in water distribution systems[J]. Mathematical Problems in Engineering,2019, 3917571:1-11.

    [9] SEIFOLLAHI-AGHMIUNI S,HADDAD O B. Effects of pipe roughness uncertainty on water distribution network performance during its operational period[J]. Water Resources Management,2013,27(5):1581-1599.

    [10] PILLER O. Modeling the behavior of a network—Hydraulic analysis and sampling procedures for parameter estimation[D]. Bordeaux;Univ.of Bordeaux,1995.

    [11] 杜坤. 供水管網(wǎng)水力模型校核與漏損定位研究[D]. 重慶:重慶大學(xué),2014.

    [12] DI N A,DI N M,GISONNI C,et al. A genetic algorithm for demand pattern and leakage estimation in a water distribution network[J]. Journal of Water Supply Research and Technology-Aqua,2014,64(1):35-46.

    [13] BEAL C,STEWART R. Identifying residential water end uses underpinning peak day and peak hour demand[J]. J.Water Resour.Plann. Manage,2014,140(7):1-10.

    [14] 李成林. 基于遺傳算法城市供水管網(wǎng)模型校核與優(yōu)化調(diào)度研究[D]. 合肥:安徽建筑大學(xué),2017.

    [15] 趙洪賓. 供水管網(wǎng)水力、水質(zhì)模型校核標(biāo)準(zhǔn)探討[J]. 中國(guó)給水排水,2008,24(18):104-106.

    [16] GOLDBERG D E. A comparative analysis of selection schemes used in genetic algorithms[J]. Foundations of genetic algorithms,1991,1:69-93.

    [17] GOLDBERG D E. Genetic algorithms in search, optimization and machine learning[D]. Tuscaloosa:The University of Alabama,1989.

    [18] NHU D,ANGUS S,JOCHEN D,et al. Demand estimation in water distribution systems:solving underdetermined problems using genetic algorithms[C]. In XVIII International Conference on Water Distribution Systems,WDSA2016,Procedia Engineering,2017,186:193-201.

    [19] BUSH C A, UBER J G. Sampling Design Methods for Water Distribution Model Calibration[J]. Journal of Water Resources Planning and Management, 1998, 124(6):334-344.

    [20] HUTTON C,KAPELAN Z,VAMVAKERIDOU-LYROUDIA L,et al. Real-time demand estimation in water distribution systems under uncertainty[C]. WDSA 2012,14th Water Distribution Systems Analysis Conference, 2012: 1374-1385.

    猜你喜歡
    校核遺傳算法
    發(fā)電機(jī)進(jìn)相試驗(yàn)失磁保護(hù)與低勵(lì)限制配合關(guān)系的校核
    通行費(fèi)校核分析方法探究
    FPSO火炬塔海運(yùn)綁扎方案及強(qiáng)度校核
    遺傳算法對(duì)CMAC與PID并行勵(lì)磁控制的優(yōu)化
    LKJ徑路數(shù)據(jù)校核系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應(yīng)用
    基于遺傳算法和LS-SVM的財(cái)務(wù)危機(jī)預(yù)測(cè)
    協(xié)同進(jìn)化在遺傳算法中的應(yīng)用研究
    基于改進(jìn)的遺傳算法的模糊聚類算法
    午夜精品在线福利| 亚洲成人久久爱视频| 日本a在线网址| 日韩视频一区二区在线观看| 哪里可以看免费的av片| 精品国产一区二区三区四区第35| 国产高清视频在线播放一区| 亚洲五月色婷婷综合| 日本熟妇午夜| 看黄色毛片网站| 国产精品亚洲av一区麻豆| 在线视频色国产色| 一区二区三区国产精品乱码| 中文在线观看免费www的网站 | 99re在线观看精品视频| www.自偷自拍.com| 日本在线视频免费播放| 90打野战视频偷拍视频| 国产精品日韩av在线免费观看| 国产私拍福利视频在线观看| 亚洲成人精品中文字幕电影| 精品人妻1区二区| 国产精品99久久99久久久不卡| 欧美日韩瑟瑟在线播放| netflix在线观看网站| 亚洲午夜理论影院| 久久久久久久久免费视频了| 精品久久蜜臀av无| cao死你这个sao货| 亚洲熟妇熟女久久| 国产精品九九99| 99在线人妻在线中文字幕| 亚洲精品美女久久久久99蜜臀| 亚洲欧美精品综合一区二区三区| 国产片内射在线| 欧洲精品卡2卡3卡4卡5卡区| 精品久久蜜臀av无| 国产成人一区二区三区免费视频网站| 91国产中文字幕| a在线观看视频网站| av视频在线观看入口| 日韩成人在线观看一区二区三区| 黄色毛片三级朝国网站| 丁香六月欧美| 人人妻人人澡人人看| 中文字幕人妻熟女乱码| 国产高清有码在线观看视频 | 精品国产国语对白av| 国产精品电影一区二区三区| 国产精品 欧美亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 国产不卡一卡二| 日韩中文字幕欧美一区二区| 欧美性长视频在线观看| 亚洲天堂国产精品一区在线| 亚洲国产欧美一区二区综合| 色播亚洲综合网| a级毛片a级免费在线| 日本精品一区二区三区蜜桃| 巨乳人妻的诱惑在线观看| 12—13女人毛片做爰片一| 少妇熟女aⅴ在线视频| 国产精品国产高清国产av| 天天一区二区日本电影三级| 脱女人内裤的视频| 中文字幕最新亚洲高清| 久久久国产成人精品二区| 久热这里只有精品99| 国产精品久久久久久亚洲av鲁大| 国产一区二区三区视频了| 久久久久久免费高清国产稀缺| 婷婷精品国产亚洲av| 欧洲精品卡2卡3卡4卡5卡区| 亚洲五月天丁香| 久久精品国产亚洲av香蕉五月| 免费一级毛片在线播放高清视频| 国产私拍福利视频在线观看| 黑人巨大精品欧美一区二区mp4| 少妇的丰满在线观看| 天堂√8在线中文| 欧美中文综合在线视频| 欧美黑人欧美精品刺激| 少妇 在线观看| 十分钟在线观看高清视频www| www.999成人在线观看| 亚洲一区二区三区不卡视频| 最新在线观看一区二区三区| 首页视频小说图片口味搜索| 给我免费播放毛片高清在线观看| 久久精品成人免费网站| 欧美激情高清一区二区三区| 婷婷精品国产亚洲av在线| 一级毛片精品| 欧美最黄视频在线播放免费| 精品久久久久久久久久久久久 | 不卡一级毛片| 在线观看66精品国产| 国产亚洲精品av在线| 国产精品乱码一区二三区的特点| www日本在线高清视频| 制服人妻中文乱码| 亚洲国产中文字幕在线视频| 老司机午夜福利在线观看视频| 在线观看免费日韩欧美大片| 欧美日韩亚洲综合一区二区三区_| 久久精品夜夜夜夜夜久久蜜豆 | 久久精品亚洲精品国产色婷小说| 中文字幕最新亚洲高清| 成年免费大片在线观看| 青草久久国产| 母亲3免费完整高清在线观看| 亚洲五月天丁香| 91麻豆精品激情在线观看国产| 国产不卡一卡二| 一a级毛片在线观看| 白带黄色成豆腐渣| 欧美性长视频在线观看| 欧美最黄视频在线播放免费| 搡老妇女老女人老熟妇| 国产又黄又爽又无遮挡在线| 一个人观看的视频www高清免费观看 | 男女下面进入的视频免费午夜 | 国产精品久久久久久亚洲av鲁大| 亚洲色图av天堂| 欧美日韩瑟瑟在线播放| 国产精品98久久久久久宅男小说| 成人特级黄色片久久久久久久| av免费在线观看网站| 国产亚洲精品久久久久久毛片| 亚洲国产高清在线一区二区三 | 巨乳人妻的诱惑在线观看| 欧美久久黑人一区二区| 午夜久久久在线观看| 老司机福利观看| √禁漫天堂资源中文www| 91老司机精品| 午夜久久久久精精品| 精品国产乱子伦一区二区三区| 国产成人精品久久二区二区免费| 国产亚洲精品一区二区www| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩成人在线观看一区二区三区| 国产精品香港三级国产av潘金莲| 99热只有精品国产| 亚洲av日韩精品久久久久久密| 一个人免费在线观看的高清视频| 欧美性长视频在线观看| 中文字幕最新亚洲高清| 久久香蕉激情| 国产黄色小视频在线观看| 久久精品国产清高在天天线| 亚洲 国产 在线| 欧美不卡视频在线免费观看 | 黄频高清免费视频| 亚洲精品美女久久久久99蜜臀| 日本 av在线| 女性生殖器流出的白浆| 麻豆av在线久日| 国产国语露脸激情在线看| 亚洲男人的天堂狠狠| 窝窝影院91人妻| 婷婷亚洲欧美| 欧美+亚洲+日韩+国产| 久久亚洲精品不卡| 午夜福利免费观看在线| 国产高清有码在线观看视频 | 女生性感内裤真人,穿戴方法视频| 黄频高清免费视频| 麻豆成人午夜福利视频| av片东京热男人的天堂| 国产又黄又爽又无遮挡在线| 亚洲一区中文字幕在线| 男人舔女人的私密视频| 午夜免费观看网址| 欧美激情高清一区二区三区| 国产亚洲av嫩草精品影院| 99久久综合精品五月天人人| 免费在线观看日本一区| 精品少妇一区二区三区视频日本电影| 国产精品久久久久久亚洲av鲁大| 国产一卡二卡三卡精品| 91av网站免费观看| 99久久国产精品久久久| 久久久久久久久免费视频了| 国产精品99久久99久久久不卡| 午夜福利在线在线| 久久久久久久久久黄片| 最近最新中文字幕大全免费视频| 99久久国产精品久久久| 人人澡人人妻人| 激情在线观看视频在线高清| 人人妻人人澡欧美一区二区| 日本精品一区二区三区蜜桃| 国产成人精品久久二区二区91| 免费在线观看亚洲国产| 亚洲人成网站高清观看| 成人国产综合亚洲| www.自偷自拍.com| 美女高潮喷水抽搐中文字幕| а√天堂www在线а√下载| 97超级碰碰碰精品色视频在线观看| 性欧美人与动物交配| 欧美色欧美亚洲另类二区| 在线观看免费午夜福利视频| 成年人黄色毛片网站| 国产精品久久视频播放| 亚洲av熟女| 日本免费a在线| 听说在线观看完整版免费高清| 淫秽高清视频在线观看| 亚洲av日韩精品久久久久久密| 亚洲第一青青草原| 香蕉av资源在线| 国产高清激情床上av| 亚洲第一av免费看| 女人高潮潮喷娇喘18禁视频| 精品电影一区二区在线| 身体一侧抽搐| 成人特级黄色片久久久久久久| 性色av乱码一区二区三区2| 精品国产一区二区三区四区第35| 久久精品aⅴ一区二区三区四区| 俄罗斯特黄特色一大片| 在线观看日韩欧美| 日本三级黄在线观看| 搡老熟女国产l中国老女人| 日本免费一区二区三区高清不卡| 久久天堂一区二区三区四区| 成人午夜高清在线视频 | 精品一区二区三区视频在线观看免费| 国产精品久久久av美女十八| www.999成人在线观看| 国产成人影院久久av| 老司机福利观看| 国产亚洲精品综合一区在线观看 | 黄网站色视频无遮挡免费观看| 欧美av亚洲av综合av国产av| 久久久水蜜桃国产精品网| 日韩大码丰满熟妇| 熟女少妇亚洲综合色aaa.| 亚洲中文日韩欧美视频| 亚洲三区欧美一区| 一a级毛片在线观看| 99国产综合亚洲精品| 日本一区二区免费在线视频| xxx96com| 国产乱人伦免费视频| 午夜福利18| 欧美性猛交黑人性爽| 欧美乱妇无乱码| 色综合欧美亚洲国产小说| 又黄又粗又硬又大视频| 日韩精品青青久久久久久| 亚洲av成人av| 两性夫妻黄色片| 久久午夜综合久久蜜桃| 亚洲专区字幕在线| 18禁观看日本| 日韩三级视频一区二区三区| 国产精品久久视频播放| 精品久久蜜臀av无| 成人永久免费在线观看视频| 国产精品久久久久久亚洲av鲁大| www.999成人在线观看| 久久精品91无色码中文字幕| 国产激情久久老熟女| 久久 成人 亚洲| 麻豆久久精品国产亚洲av| 亚洲第一青青草原| 欧美成人性av电影在线观看| 桃色一区二区三区在线观看| 一级毛片高清免费大全| 亚洲一区二区三区色噜噜| 亚洲,欧美精品.| 18禁裸乳无遮挡免费网站照片 | 老熟妇乱子伦视频在线观看| 在线视频色国产色| 听说在线观看完整版免费高清| 无人区码免费观看不卡| 免费在线观看日本一区| 在线观看www视频免费| 亚洲成人国产一区在线观看| 男女视频在线观看网站免费 | av片东京热男人的天堂| 中文字幕高清在线视频| 天天躁夜夜躁狠狠躁躁| 色综合站精品国产| 18禁黄网站禁片午夜丰满| 亚洲精品av麻豆狂野| 少妇粗大呻吟视频| 波多野结衣高清作品| 一本久久中文字幕| 日韩欧美国产在线观看| 日韩中文字幕欧美一区二区| 国产91精品成人一区二区三区| 一本久久中文字幕| 精品第一国产精品| 亚洲国产高清在线一区二区三 | av超薄肉色丝袜交足视频| 亚洲男人的天堂狠狠| 色老头精品视频在线观看| avwww免费| 欧美zozozo另类| 欧美中文日本在线观看视频| 搞女人的毛片| 嫩草影视91久久| 午夜福利欧美成人| 日日干狠狠操夜夜爽| 校园春色视频在线观看| 成人国产综合亚洲| 99re在线观看精品视频| 免费人成视频x8x8入口观看| 婷婷丁香在线五月| 久久久久国产精品人妻aⅴ院| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲成av片中文字幕在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 一级毛片高清免费大全| 国产精品影院久久| 日本免费a在线| 亚洲av电影不卡..在线观看| 久久香蕉精品热| 国产午夜精品久久久久久| 国产区一区二久久| 亚洲色图av天堂| 99热只有精品国产| 18禁黄网站禁片午夜丰满| 在线观看舔阴道视频| 欧美av亚洲av综合av国产av| 色播亚洲综合网| 在线播放国产精品三级| 欧美三级亚洲精品| 视频区欧美日本亚洲| 久久草成人影院| 777久久人妻少妇嫩草av网站| 91麻豆精品激情在线观看国产| 午夜免费激情av| 人妻丰满熟妇av一区二区三区| 91字幕亚洲| 免费观看人在逋| 最新在线观看一区二区三区| 黄色毛片三级朝国网站| 久久国产乱子伦精品免费另类| 成人午夜高清在线视频 | 波多野结衣av一区二区av| 一级片免费观看大全| 九色国产91popny在线| 自线自在国产av| 精品国产乱码久久久久久男人| 欧美成人免费av一区二区三区| 成人特级黄色片久久久久久久| 99在线人妻在线中文字幕| 国产亚洲欧美98| 国产精品99久久99久久久不卡| 欧美黄色片欧美黄色片| 日本撒尿小便嘘嘘汇集6| 精品福利观看| 久久精品影院6| 国产区一区二久久| av免费在线观看网站| 一二三四社区在线视频社区8| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品合色在线| 香蕉av资源在线| 午夜精品久久久久久毛片777| 国产熟女xx| 色综合亚洲欧美另类图片| 中文字幕人妻熟女乱码| 无限看片的www在线观看| 久久精品国产亚洲av香蕉五月| 亚洲色图 男人天堂 中文字幕| 男男h啪啪无遮挡| 国产又色又爽无遮挡免费看| 18美女黄网站色大片免费观看| 亚洲 欧美一区二区三区| 无遮挡黄片免费观看| 男人的好看免费观看在线视频 | 精品高清国产在线一区| 亚洲国产毛片av蜜桃av| 色精品久久人妻99蜜桃| 一级毛片高清免费大全| 日韩 欧美 亚洲 中文字幕| 成年人黄色毛片网站| netflix在线观看网站| 国产高清激情床上av| 亚洲成人久久爱视频| 搡老熟女国产l中国老女人| 国产av一区二区精品久久| 亚洲成人久久爱视频| 欧美性长视频在线观看| 日韩欧美一区视频在线观看| 99精品久久久久人妻精品| 99国产精品一区二区三区| 宅男免费午夜| 国产视频内射| 久久香蕉精品热| 久久中文字幕人妻熟女| 午夜福利成人在线免费观看| 精品不卡国产一区二区三区| 久久久水蜜桃国产精品网| 老司机午夜福利在线观看视频| 女同久久另类99精品国产91| 久久久久久大精品| 国产精品久久电影中文字幕| 12—13女人毛片做爰片一| 日日干狠狠操夜夜爽| 久久青草综合色| 国产高清激情床上av| 99热6这里只有精品| 女人爽到高潮嗷嗷叫在线视频| 日韩精品免费视频一区二区三区| 亚洲国产欧洲综合997久久, | 国产成人av激情在线播放| 久久欧美精品欧美久久欧美| 日韩大码丰满熟妇| 一本综合久久免费| 欧美黑人欧美精品刺激| 可以免费在线观看a视频的电影网站| 亚洲三区欧美一区| 欧美日韩乱码在线| 日本一本二区三区精品| 精品国产一区二区三区四区第35| 亚洲色图 男人天堂 中文字幕| 久久中文字幕一级| √禁漫天堂资源中文www| 老汉色∧v一级毛片| 男女之事视频高清在线观看| 久久久久精品国产欧美久久久| 一区二区三区高清视频在线| 人成视频在线观看免费观看| 日日摸夜夜添夜夜添小说| 亚洲第一av免费看| 在线观看免费午夜福利视频| 亚洲欧美精品综合一区二区三区| 亚洲av日韩精品久久久久久密| 中文在线观看免费www的网站 | 欧美日本亚洲视频在线播放| 精品欧美国产一区二区三| 激情在线观看视频在线高清| 成人亚洲精品一区在线观看| 级片在线观看| 最好的美女福利视频网| cao死你这个sao货| 色播在线永久视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产极品粉嫩免费观看在线| 久热爱精品视频在线9| 黄色视频,在线免费观看| 精品一区二区三区av网在线观看| 国产精品久久久久久亚洲av鲁大| 久久九九热精品免费| 欧美日韩瑟瑟在线播放| 啦啦啦 在线观看视频| 色尼玛亚洲综合影院| 亚洲国产精品sss在线观看| 欧美人与性动交α欧美精品济南到| 国产黄片美女视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲国产中文字幕在线视频| 男人舔奶头视频| 欧美成狂野欧美在线观看| 色老头精品视频在线观看| 俄罗斯特黄特色一大片| 亚洲欧美日韩高清在线视频| 一区福利在线观看| 国产精品亚洲av一区麻豆| 久久久久久九九精品二区国产 | 国产亚洲av嫩草精品影院| 成人永久免费在线观看视频| 欧美激情极品国产一区二区三区| 精品久久久久久久久久久久久 | 国产黄色小视频在线观看| 精品久久久久久久毛片微露脸| 人妻久久中文字幕网| 亚洲成人久久性| 十八禁人妻一区二区| 免费无遮挡裸体视频| 一个人免费在线观看的高清视频| 精品久久久久久久久久免费视频| 中文字幕人成人乱码亚洲影| 欧美日韩一级在线毛片| 欧美三级亚洲精品| 亚洲国产欧美一区二区综合| 亚洲va日本ⅴa欧美va伊人久久| 淫秽高清视频在线观看| 高潮久久久久久久久久久不卡| 欧美黑人精品巨大| 久久精品国产99精品国产亚洲性色| 每晚都被弄得嗷嗷叫到高潮| 十八禁人妻一区二区| 美女国产高潮福利片在线看| 欧美黑人欧美精品刺激| 欧美人与性动交α欧美精品济南到| 一级毛片女人18水好多| 久久久久久免费高清国产稀缺| 亚洲专区国产一区二区| 老司机靠b影院| 亚洲国产精品sss在线观看| 国产蜜桃级精品一区二区三区| 婷婷六月久久综合丁香| 国产国语露脸激情在线看| 在线观看免费视频日本深夜| 亚洲av第一区精品v没综合| 精品国产国语对白av| 久久香蕉国产精品| 欧美久久黑人一区二区| 香蕉av资源在线| 国内久久婷婷六月综合欲色啪| 美女免费视频网站| 99久久精品国产亚洲精品| 最近最新中文字幕大全免费视频| 亚洲精品国产一区二区精华液| 美女大奶头视频| 国产成人一区二区三区免费视频网站| 亚洲精品国产区一区二| 91国产中文字幕| 侵犯人妻中文字幕一二三四区| 51午夜福利影视在线观看| 在线观看舔阴道视频| 亚洲狠狠婷婷综合久久图片| 丁香六月欧美| 变态另类丝袜制服| 人人澡人人妻人| 亚洲国产毛片av蜜桃av| 视频在线观看一区二区三区| 老司机在亚洲福利影院| 欧美激情极品国产一区二区三区| 久久久国产成人精品二区| 欧美一区二区精品小视频在线| 麻豆成人午夜福利视频| 老汉色∧v一级毛片| 成年人黄色毛片网站| 动漫黄色视频在线观看| 香蕉国产在线看| 亚洲 欧美 日韩 在线 免费| a在线观看视频网站| 婷婷精品国产亚洲av| 亚洲国产欧美一区二区综合| 白带黄色成豆腐渣| 久久亚洲真实| 亚洲午夜精品一区,二区,三区| 色播亚洲综合网| 亚洲午夜精品一区,二区,三区| av视频在线观看入口| 老司机深夜福利视频在线观看| 国产精品精品国产色婷婷| 看黄色毛片网站| 露出奶头的视频| 91在线观看av| 制服丝袜大香蕉在线| 99国产精品一区二区三区| 国产精品98久久久久久宅男小说| 日本免费a在线| 久久精品亚洲精品国产色婷小说| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品一区二区三区| 最近最新中文字幕大全免费视频| 女性被躁到高潮视频| 夜夜爽天天搞| 亚洲五月色婷婷综合| 成人手机av| 国产午夜福利久久久久久| 后天国语完整版免费观看| 日韩av在线大香蕉| 欧美激情久久久久久爽电影| 最近最新免费中文字幕在线| 亚洲欧美精品综合一区二区三区| 99久久无色码亚洲精品果冻| 欧美不卡视频在线免费观看 | 久久亚洲真实| 熟女电影av网| 国产色视频综合| 每晚都被弄得嗷嗷叫到高潮| 可以在线观看的亚洲视频| 国产成人欧美| 热99re8久久精品国产| 99国产极品粉嫩在线观看| 首页视频小说图片口味搜索| 亚洲欧美日韩无卡精品| 精品久久蜜臀av无| 男人舔女人下体高潮全视频| 天堂影院成人在线观看| 中文字幕人妻丝袜一区二区| 我的亚洲天堂| 亚洲av日韩精品久久久久久密| 欧美日本视频| 国产精品av久久久久免费| 中文字幕精品免费在线观看视频| 久久久久久久精品吃奶| 久久香蕉精品热| 久久久国产成人免费| av天堂在线播放| 91麻豆av在线| 国产精品久久久人人做人人爽| 中文字幕久久专区| 日韩视频一区二区在线观看| 亚洲最大成人中文| 亚洲精华国产精华精| 黑人巨大精品欧美一区二区mp4| 亚洲色图 男人天堂 中文字幕| 妹子高潮喷水视频| 欧美成人一区二区免费高清观看 | 女人高潮潮喷娇喘18禁视频| 欧美日韩一级在线毛片| 亚洲精品在线美女| 亚洲精品久久国产高清桃花| 亚洲色图av天堂| 久久久精品欧美日韩精品| 成人一区二区视频在线观看| 久久久久久久久久黄片| 国内久久婷婷六月综合欲色啪| 国产精品久久久av美女十八| 国产视频内射| 免费看a级黄色片| 欧美亚洲日本最大视频资源|