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

    基于自動生成策略的供水管網(wǎng)抗震優(yōu)化算法

    2016-07-26 09:46:16元宏偉
    關(guān)鍵詞:模擬退火算法供水管網(wǎng)遺傳算法

    劉 威, 元宏偉, 徐 良

    (1.同濟大學(xué) 土木工程防災(zāi)國家重點實驗室,上海 200092;2. 同濟大學(xué) 土木工程學(xué)院,上海 200092;3. 北京市建筑設(shè)計研究院,北京 100045;4. 上海核工程研究設(shè)計院,上海 200233)

    ?

    基于自動生成策略的供水管網(wǎng)抗震優(yōu)化算法

    劉威1,2, 元宏偉3, 徐良4

    (1.同濟大學(xué) 土木工程防災(zāi)國家重點實驗室,上海 200092;2. 同濟大學(xué) 土木工程學(xué)院,上海 200092;3. 北京市建筑設(shè)計研究院,北京 100045;4. 上海核工程研究設(shè)計院,上海 200233)

    摘要:以管網(wǎng)年費用折算值為優(yōu)化目標(biāo)、管網(wǎng)拓撲結(jié)構(gòu)與管徑為優(yōu)化參數(shù)、管網(wǎng)節(jié)點抗震可靠度為約束條件, 建立了供水管網(wǎng)抗震優(yōu)化設(shè)計模型.基于自動生成策略,并結(jié)合環(huán)形管網(wǎng)判斷方法,分別利用遺傳算法、遺傳-模擬退火算法和微粒群算法進行了供水管網(wǎng)的抗震拓撲優(yōu)化分析.利用3種優(yōu)化方法對2個典型供水管網(wǎng)進行了對比分析.對比分析表明,遺傳-模擬退火算法具有最好的優(yōu)化能力.

    關(guān)鍵詞:供水管網(wǎng); 自動生成策略; 遺傳算法; 遺傳-模擬退火算法; 微粒群算法

    生命線工程系統(tǒng)是指維系現(xiàn)代城市功能與區(qū)域經(jīng)濟功能的基礎(chǔ)性工程設(shè)施系統(tǒng)[1].作為城市生命線系統(tǒng)的重要組成部分,城市供水管網(wǎng)對城市的日常運行和災(zāi)后搶險救災(zāi)具有重要意義.進行城市供水管網(wǎng)抗震優(yōu)化設(shè)計,確保城市供水系統(tǒng)在地震中和地震后安全可靠運行,減少地震災(zāi)害的損失,加快災(zāi)后城市恢復(fù)的時間,是亟待解決的問題.在這方面的研究中,Chen等[2]首次以管網(wǎng)建設(shè)造價為優(yōu)化目標(biāo),利用遺傳算法并結(jié)合啟發(fā)式管徑優(yōu)化算法進行了基于功能可靠度的供水管網(wǎng)抗震拓撲優(yōu)化和管徑優(yōu)化研究.進而,李杰等[3]、邢燕等[4]分別利用模擬退火算法和遺傳-模擬退火混合算法對供水管網(wǎng)進行了抗震拓撲優(yōu)化設(shè)計.徐良等[5]以管網(wǎng)年費用折算值為優(yōu)化目標(biāo)、以管網(wǎng)拓撲結(jié)構(gòu)與管段管徑為優(yōu)化參數(shù)、以管網(wǎng)節(jié)點最低可靠度為約束條件,建立供水管網(wǎng)抗震優(yōu)化設(shè)計模型,并將微粒群算法應(yīng)用到管網(wǎng)系統(tǒng)抗震優(yōu)化問題中.上述研究都從一個初始經(jīng)驗管網(wǎng)出發(fā)來進行簡化,優(yōu)化結(jié)果對初始經(jīng)驗管網(wǎng)依賴很大.2010年,李杰等[6]首先提出了基于自動生成策略的供水管網(wǎng)抗震優(yōu)化設(shè)計方法.這一方法只需要管網(wǎng)用戶節(jié)點信息就可以生成滿足抗震要求的供水管網(wǎng),但這一方法在生成管線時不考慮實際工程背景,而是直接認為管線是連接管網(wǎng)節(jié)點的直線,與實際工程尚有差距.

    在上述工作基礎(chǔ)上,本文將供水管網(wǎng)所在區(qū)域的道路信息引入管網(wǎng)優(yōu)化設(shè)計中,并結(jié)合環(huán)形網(wǎng)絡(luò)判斷方法,提出了新的自動生成網(wǎng)絡(luò)拓撲結(jié)構(gòu)的策略.基于此策略,本文利用遺傳算法、遺傳-模擬退火算法和微粒群算法進行了城市供水管網(wǎng)系統(tǒng)抗震拓撲優(yōu)化分析,并結(jié)合算例進行了對比分析.

    1供水管網(wǎng)的抗震優(yōu)化模型

    以城市供水管網(wǎng)拓撲結(jié)構(gòu)和管徑為優(yōu)化參數(shù)、抗震可靠度作為約束條件、管網(wǎng)建造運營費用最小為設(shè)計目標(biāo)來進行供水管網(wǎng)優(yōu)化設(shè)計,可建立如下供水管網(wǎng)抗震優(yōu)化模型[5]:

    minW

    (1)

    式中:βmin為供水管網(wǎng)所有節(jié)點的抗震可靠度指標(biāo)的最小值,可以采用文獻[7]介紹的均值一次二階矩方法來求解;β0為供水管網(wǎng)設(shè)計時允許的抗震可靠指標(biāo)限值;W為供水管網(wǎng)的年費用折算值,可用下式表示[8]:

    (2)

    式中:p為管網(wǎng)的每年折舊和大修的百分率;T為管網(wǎng)建設(shè)投資償還期,年;n為泵站數(shù)目;WPi為泵站i的單位運行電費指標(biāo),元·(m-3·s·m·年);qi為泵站i的設(shè)計揚水流量,m3·s-1;hpi為泵站i的最大揚程,m;WN為管網(wǎng)造價,可用下式計算得到:

    (3)

    式中:lj為管線j的長度,m;dj為管線j的管徑,m;a1,a2,a3為造價經(jīng)驗系數(shù),可以分別取62.105 1,1 979.7,1.486[9];γj為管線j的連通系數(shù),鋪設(shè)管線時取為1,不鋪設(shè)時取為0;m為管網(wǎng)中管線的數(shù)目.

    2供水管網(wǎng)自動生成策略

    自動生成技術(shù)的含義是從某種已有拓撲結(jié)構(gòu)(包括空白結(jié)構(gòu))出發(fā),按照一定的規(guī)則,將單元增加到結(jié)構(gòu)中去,并不斷地進行調(diào)整,直至最后形成滿足預(yù)先給定要求的拓撲構(gòu)形為止.在供水管網(wǎng)抗震優(yōu)化設(shè)計中,采用自動生成技術(shù)可以生成優(yōu)化算法的初始種群.考慮到在實際工程中供水管網(wǎng)應(yīng)盡可能沿道路延伸,本文將供水管網(wǎng)所在區(qū)域的道路信息引入管網(wǎng)優(yōu)化設(shè)計中,建立了基于道路信息的供水管網(wǎng)自動生成策略.

    相比于樹狀供水管網(wǎng),環(huán)狀供水管網(wǎng)有著更好的可靠性.因此,本文將環(huán)形管網(wǎng)生成方法引入自動生成策略中,實現(xiàn)了供水管網(wǎng)的優(yōu)化設(shè)計.

    2.1建立管網(wǎng)的初始信息

    將管網(wǎng)所在區(qū)域的道路信息用道路節(jié)點間的路徑來表示,同時將管網(wǎng)節(jié)點放置于路徑之中.這樣在初始狀態(tài)下,管網(wǎng)信息由道路節(jié)點、道路節(jié)點間的路徑以及管網(wǎng)節(jié)點信息組成.同時,根據(jù)實際情況可以給出道路節(jié)點以及管網(wǎng)節(jié)點間的路徑長度.

    2.2判斷能否生成環(huán)形網(wǎng)

    2.3管線生成規(guī)則

    管線的生成規(guī)則為:①當(dāng)管線兩節(jié)點之間沒有通過其他管網(wǎng)節(jié)點或只通過道路節(jié)點,則生成1條管線;②當(dāng)兩節(jié)點之間的管線通過其他管網(wǎng)節(jié)點,則以這些節(jié)點為端點生成多條管線.

    上述規(guī)則保證任意一條管線以節(jié)點為端點,中間不通過其他節(jié)點.以圖1d中的網(wǎng)絡(luò)為例(節(jié)點8,9和10為道路節(jié)點),若節(jié)點7,5間生成管線,即為管線7-5,管線的連接方式為7-5,長度為500 m;若節(jié)點7,6間生成管線,管線的連接方式為7-9-6,由于9是道路節(jié)點,不是管網(wǎng)節(jié)點,故生成管線7-6,長度為900 m;若節(jié)點7,1間生成管線,管線的連接方式為7-3-1,由于節(jié)點3為管網(wǎng)節(jié)點,則生成管線7-3和3-1,長度分別為500 m和400 m.

    2.4生成初始管網(wǎng)拓撲結(jié)構(gòu)

    在優(yōu)化算法中,首先要生成一系列初始管網(wǎng).在生成初始管網(wǎng)時,按照管網(wǎng)節(jié)點的編號順序逐個遍歷所有節(jié)點.對每個管網(wǎng)節(jié)點,生成它與其他節(jié)點之間的管線.生成方法為:對處于管網(wǎng)邊緣的節(jié)點,利用Dijkstra算法[10]找出與其路徑最短的3個節(jié)點;對處于管網(wǎng)中心的節(jié)點,找出與其路徑最短的4個節(jié)點,同時以一定概率生成它們之間的管線.

    a管網(wǎng)1b管網(wǎng)2

    c 管網(wǎng)3

    d 管網(wǎng)4

    圖2 待建供水管網(wǎng)

    如圖2所示待建供水管網(wǎng),4為水源點,1-3和5-12為需水節(jié)點, 13-19為道路節(jié)點.以邊緣節(jié)點1為例,如圖3a所示,與其路徑最短的3個節(jié)點為節(jié)點3,4和2.若節(jié)點1和3間生成管線,即為管線1-3,管線的連接方式為1-19-3.節(jié)點4為中間節(jié)點,如圖3b所示,與其路徑最短的4個節(jié)點分別為節(jié)點9,2,1和6.若節(jié)點4和2間生成管線,即為管線4-2,管線的連接方式為4-13-2;若節(jié)點4,9間生成管線,即為管線4-9,管線的連接方式為4-15-9.對每個節(jié)點進行一次操作,得到的一個可能的管網(wǎng)拓撲結(jié)構(gòu)如圖4所示.

    2.5環(huán)狀修補過程

    需要指出的是,在環(huán)狀管網(wǎng)設(shè)計過程中,可能會產(chǎn)生一些無意義的解,即:出現(xiàn)節(jié)點度為零的節(jié)點,從而造成1個管網(wǎng)斷成2個或多個子網(wǎng);或者生成的管網(wǎng)為非環(huán)形網(wǎng).因此要對每一次生成的管網(wǎng)進行環(huán)狀修補,修補過程如下:

    (1)將初始管網(wǎng)個體中生成的一條管線暫時刪除,同時暫時刪除這條管線所經(jīng)過道路的道路信息.

    (2)對此個體在此道路信息下進行連通性修補.

    a生成與節(jié)點1連接的管線b生成與節(jié)點4連接的管線圖3 管網(wǎng)生成過程Fig.3 Processofnetworkgeneration圖4 管網(wǎng)拓撲結(jié)構(gòu)Fig.4 Networktopology

    首先,采用廣度優(yōu)先搜索,判斷生成的管網(wǎng)拓撲結(jié)構(gòu)的連通性,若存在與水源點不連通的節(jié)點,則對其進行連通性修補.進行修補時,對每個與水源點不連通的節(jié)點按前述自動生成管線的方式來增加管線,即:對處于管網(wǎng)邊緣的不連通節(jié)點,找出與其連接路徑最短的3個與水源點連通的節(jié)點,若此時與水源點連通的節(jié)點數(shù)小于3,則找出所有與水源點連通的節(jié)點;對處于管網(wǎng)中心的不連通節(jié)點,找出與其連接路徑最短的4個與水源點連通的節(jié)點,若此時與水源點連通的節(jié)點數(shù)小于4,則找出所有與水源點連通的節(jié)點.同時,以一定概率生成它們間的管線.

    (3)將暫時刪除的那條管線和道路信息還原.

    (4)對個體中每條已生成的管線依次重復(fù)步驟(1)~(3).

    (5)對修補后的個體判斷是否為環(huán)狀網(wǎng).若非環(huán)狀網(wǎng),則刪除此個體,重新生成.

    圖5是對圖4進行環(huán)狀修補后得到的環(huán)狀管網(wǎng).

    圖5 修補生成的環(huán)狀網(wǎng)

    3基于自動生成策略的優(yōu)化算法在供水管網(wǎng)抗震優(yōu)化中的應(yīng)用

    在管網(wǎng)優(yōu)化過程中,可以采用生物進化的原則構(gòu)造管網(wǎng)拓撲優(yōu)化問題的優(yōu)化算法.

    3.1遺傳算法

    遺傳算法是借鑒生物進化原則在20世紀70年代初期由美國密西根大學(xué)的Holland[11]教授提出的一種自適應(yīng)并行全局優(yōu)化概率搜索算法.其計算步驟包括編碼、生成初始種群、個體評價、選擇操作、交叉操作、變異操作和收斂性判斷等.

    (1)編碼.把一個問題的解從其解空間轉(zhuǎn)換到遺傳算法所能處理的搜索空間的過程稱為編碼.假定一根管線管徑可取100,150,200,250,300,350,400,450,500 mm,則對應(yīng)的基因可取1,2,3,4,5,6,7,8,9,而基因取0值則表示這根管線不鋪設(shè).所有優(yōu)化參數(shù)對應(yīng)的基因按照指定的次序排列起來,就構(gòu)成一條染色體.在管網(wǎng)拓撲優(yōu)化模型中一個染色體對應(yīng)管網(wǎng)的一種拓撲結(jié)構(gòu)方案.多條染色體構(gòu)成遺傳算法的一個種群.

    (2)生成初始種群.初始種群的生成由自動生成策略來實現(xiàn).

    (3)個體評價.通過計算個體適應(yīng)度來評價個體,這里定義個體s(個體也即優(yōu)化的一個解)的適應(yīng)度函數(shù)為

    (4)

    式中:M為預(yù)先指定的一個較大的數(shù)值;W(s)為個體s對應(yīng)管網(wǎng)的年費用折算值;P1(s),P2(s)為懲罰函數(shù),可采用下式來計算懲罰函數(shù):

    徐曉春等(2004)對民樂銅礦區(qū)的輝銅礦化英安斑巖等樣品進行了Sm-Nd同位素年齡測定,等時線年齡為228±56 Ma,說明與該套中三疊世火山巖地層的時代屬于同一時代。

    (5)

    P2(s)=

    (6)

    (4)選擇操作.選擇操作是根據(jù)計算得到的個體適應(yīng)度選擇優(yōu)勝的個體進入下一代,并淘汰劣質(zhì)個體.本文采用最佳個體保存方法和適應(yīng)度比例方法.最佳個體保存方法是將上一代最優(yōu)的個體直接選擇進入下一代,并將最優(yōu)個體從群體中刪除; 適應(yīng)度比例方法即按群體中各個體的適應(yīng)度來確定個體被選擇的概率.當(dāng)選擇出來的個體達到規(guī)定的數(shù)目時,選擇操作結(jié)束,本文取為群體個體總數(shù)的1/3.

    (5)交叉操作.交叉操作是遺傳算法中的重要操作之一,其以一定的交叉概率把2個父代個體的部分結(jié)構(gòu)加以替換重組從而生成新個體.這里采用單點交叉方法,其操作過程是在個體串中隨機設(shè)定一個交叉點,對隨機選擇出來的2個父代個體在該點前的部分結(jié)構(gòu)進行互換,從而生成2個新個體.這里對每一次交叉操作生成的新個體進行環(huán)狀網(wǎng)判斷及相應(yīng)的環(huán)狀修補.

    (6)變異操作.變異操作是指以一定的變異概率將某個染色體編碼串中的基因用其他等位基因進行替換.這里對交叉獲得的新管網(wǎng)通過對其中的管線以一定的概率用其他直徑的管線代替來形成一個新的個體.在這過程中,為了引導(dǎo)管網(wǎng)個體向較好的方向發(fā)展,對滿足最小可靠度指標(biāo)要求的節(jié)點,以一定概率刪除以此節(jié)點為端點的管線;對不滿足最小可靠度指標(biāo)要求的節(jié)點,按生成管線的方式以一定概率生成與其他節(jié)點連接的管線.對新個體進行環(huán)狀網(wǎng)判斷及相應(yīng)的環(huán)狀修補.

    (7)收斂判斷.本文采用固定迭代次數(shù)作為算法終止條件.

    3.2遺傳-模擬退火算法

    將模擬退火算法[12]嵌入到遺傳算法中發(fā)展出來的遺傳-模擬退火混合算法[13]可有效地增強遺傳算法的搜索能力.混合算法的計算過程與遺傳算法基本一樣,只是在變異操作時采用模擬退火操作.下面結(jié)合管網(wǎng)優(yōu)化問題來介紹模擬退火操作的步驟.

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

    模擬退火操作針對交叉操作得到的解來進行,每次選擇一個解作為當(dāng)前的初始解,然后對當(dāng)前解進行隨機擾動產(chǎn)生一個新解,根據(jù)式(7)確定接受新解的概率P,也即用解s2代替解s1的概率.

    (7)

    式中:t為當(dāng)前的溫度;f(s)為解s的目標(biāo)函數(shù).

    (8)

    式中:W(s)為解s對應(yīng)管網(wǎng)的年費用折算值;P3(s)為懲罰函數(shù),可以采用下式來計算:

    (9)

    式中:b1,b2為常數(shù),其值的大小需要根據(jù)具體的管網(wǎng)來確定.

    3.2.2隨機擾動模型

    3.2.3溫度參數(shù)的控制

    溫度參數(shù)是模擬退火過程中最關(guān)鍵的參數(shù),主要包括起始溫度的選取、溫度的下降方法、每一溫度迭代次數(shù)的確定等.

    一般采用t0=Cδ來估計初始溫度.其中,C為通過試算得到的一個充分大的數(shù),δ=Wmax-Wmin,Wmax,Wmin為父代群體中個體年費用折算值的最大值和最小值.

    模擬退火操作要求溫度下降到零,整個系統(tǒng)以概率1收斂到全局最優(yōu)解.本文所用溫度下降的方法為tr′+1=αtr′,其中:r′為模擬退火操作的計算代數(shù);α為降溫速率,0<α<1,α越接近于1溫度下降越慢,這種方法簡單易行,它的每一步以相同的比率降溫.

    本文采用固定長度法來給定每一溫度的迭代長度.

    3.3微粒群算法

    微粒群算法[14]是通過模擬鳥群的捕食行為而建立的優(yōu)化算法.在微粒群算法中,設(shè)待優(yōu)化問題的解在一個d維的搜索空間中,群體中的第s個微粒位置表示待優(yōu)化問題的一個解,可表示為一個d維矢量,Xs=(xs1,xs2,…,xsd)T;微粒s的速度表示其位置的改變,可用Vs=(vs1,vs2,…,vsd)T來表示.當(dāng)前群體的最好位置的微粒l用Kl=(kl1,kl2,…,kld)T表示,整個演化過程中歷史最好位置的微粒g用Kg=(kg1,kg2,…,kgd)T表示.則微粒群中,對每一代微粒s的第u維的進化方程為

    (10)

    (11)

    式中:γ為慣性權(quán)重,取值越大全局搜索能力越強,取值越小局部搜索能力越好;r為進化代數(shù);R1()和R2()為2個隨機函數(shù),產(chǎn)生0~1之間的隨機數(shù);c1,c2為加速度常數(shù),取較大值可以拓展微粒的搜索空間,取較小值則可提高微粒的搜索精度.

    將微粒群算法應(yīng)用到管網(wǎng)抗震優(yōu)化設(shè)計中時,每個微粒的位置對應(yīng)一個管網(wǎng)結(jié)構(gòu)方案,即問題的一個解.其編碼方式與遺傳算法相同.在進化過程中,當(dāng)位置矢量中某一維的值小于零時,則取零;當(dāng)其值大于最大管徑編碼值時,則取最大管徑編碼值;處于兩者間時,則選擇與其最接近的管徑編碼值.

    管網(wǎng)抗震優(yōu)化問題的優(yōu)化目標(biāo)是為了獲得在滿足節(jié)點最低抗震可靠度條件下管網(wǎng)最低的年費用折算值,于是對不滿足約束條件的解,通過增加年費用折算值對其懲罰.解s的適應(yīng)值函數(shù)取為

    (12)

    式中:P4為懲罰函數(shù),可取足夠大的常數(shù).

    管網(wǎng)抗震優(yōu)化的微粒群算法的步驟如下:①微粒群個體的初始位置由自動生成技術(shù)來實現(xiàn),并隨機產(chǎn)生初始速度;②計算每個微粒的適應(yīng)值,對不滿足約束條件的微粒采用罰函數(shù)進行處理;③對每個微粒,將其適應(yīng)值與其當(dāng)前群體最好位置的適應(yīng)值進行比較,若較好,則將其作為當(dāng)前群體最好位置;④將當(dāng)前群體最好位置的適應(yīng)值與整個演化過程中歷史最好位置的適應(yīng)值進行比較,若較好,則將其作為歷史最好位置;⑤分別對每個微粒的各維的速度和位置進行更新,并對無效的管網(wǎng)結(jié)構(gòu)進行修補,如未達到一個預(yù)設(shè)的最大迭代數(shù),則返回步驟②.

    3.4供水管網(wǎng)抗震優(yōu)化過程

    城市供水管網(wǎng)系統(tǒng)的抗震優(yōu)化設(shè)計過程如圖6所示.整個設(shè)計過程依托于優(yōu)化算法的大框架,將自動生成策略嵌入此大框架之中,進而完成整個設(shè)計過程,最終給出優(yōu)化的設(shè)計管網(wǎng).

    圖6 供水管網(wǎng)抗震設(shè)計過程

    4實例分析

    4.1實例1

    圖7為一個20個節(jié)點的待建供水管網(wǎng),其中節(jié)點1~19為用戶節(jié)點,20為水源點,21~33為道路節(jié)點,虛線給出了待建供水管網(wǎng)所在區(qū)域的道路分布狀況.考慮地震烈度為8度進行管網(wǎng)抗震優(yōu)化設(shè)計.管材為鑄鐵.約束條件取最低可靠度指標(biāo)為0.84(對應(yīng)節(jié)點最小可靠度0.80).根據(jù)工程實際,選定管網(wǎng)中管線的管徑種類和管徑值,分別為160,200,250,300,350,400,450,500,550 mm 這9 種,對應(yīng)的管徑編碼取為1,2,3,4,5,6,7,8,9.取種群規(guī)模為60,迭代次數(shù)為50,管線生成概率取0.9,各優(yōu)化算法基本參數(shù)如下.

    圖7 實例1待建供水管網(wǎng)

    遺傳算法:交叉概率采用動態(tài)自適應(yīng)算法來確定,即:若個體適應(yīng)度大于群體的平均適應(yīng)度,則交叉率取0.90;反之,則交叉率取0.99.變異概率取為0.2.

    遺傳-模擬退火算法:交叉概率同遺傳算法,初始溫度為100度,當(dāng)前循環(huán)溫度為上一代的0.4倍,內(nèi)循環(huán)次數(shù)為10次.

    微粒群算法:慣性權(quán)重γ從0.9至0.4線性減小,加速度常數(shù)c1和c2均從1.5至0.4線性減小.

    利用前述3種方法分別對待建管網(wǎng)進行抗震優(yōu)化設(shè)計.圖8給出了3種方法的供水管網(wǎng)優(yōu)化方案.表1給出了3種方法所得優(yōu)化方案的管網(wǎng)年費用折算值,表2給出了3種優(yōu)化算法給出的優(yōu)化結(jié)果中各個節(jié)點的可靠度指標(biāo).從表1可以看出,遺傳-模擬退火算法的費用最低,微粒群算法次之,遺傳算法最高,但遺傳算法與微粒群算法的費用差別不大.從表2可以看出,遺傳-模擬退火算法得到的結(jié)果中節(jié)點可靠度指標(biāo)最高,平均值達到2.012 3;遺傳算法次之,節(jié)點可靠度指標(biāo)平均值為1.504 2;微粒群算法最低,為1.343 9.因此,遺傳-模擬退火算法得到的優(yōu)化結(jié)果年費用折算值低,節(jié)點可靠度指標(biāo)均值高,表現(xiàn)最好.遺傳算法雖然年費用折算值比微粒群算法略高,但是得到的節(jié)點可靠度指標(biāo)均值比微粒群算法好.因此本例中,遺傳-模擬退火算法表現(xiàn)最好,遺傳算法和微粒群算法表現(xiàn)相當(dāng).

    a 基于微粒群算法的管網(wǎng)優(yōu)化方案

    b 基于遺傳算法的管網(wǎng)優(yōu)化方案

    c 基于遺傳-模擬退火算法的管網(wǎng)優(yōu)化方案

    圖8 管網(wǎng)優(yōu)化結(jié)果

    表1 實例1三種優(yōu)化算法的管網(wǎng)年費用折算值

    表2 實例1不同優(yōu)化算法的節(jié)點可靠度指標(biāo)

    4.2實例2

    圖9為綿竹市城區(qū)供水管網(wǎng),這里僅保留了各需水點、水源點和道路信息,其中節(jié)點1~58為需水點,59~62為水源點,63~64為道路節(jié)點,虛線給出了待建供水管網(wǎng)所在區(qū)域的道路分布狀況.利用前述3種方法分別對管網(wǎng)重新進行抗震優(yōu)化設(shè)計.考慮到本例中存在與管網(wǎng)由單條管線相連的水源點,因而在管網(wǎng)優(yōu)化設(shè)計中加入判定規(guī)則,以使管網(wǎng)中此類管線在優(yōu)化進程中的環(huán)狀網(wǎng)判斷中不被刪除,從而保證所有水源點都與管網(wǎng)相連.考慮地震烈度為8度進行管網(wǎng)抗震優(yōu)化設(shè)計.管材為鑄鐵.約束條件取最低可靠度指標(biāo)0.84(對應(yīng)節(jié)點最小可靠度0.8).各算法基本參數(shù)選取與實例1同.表3給出3種方法所得優(yōu)化方案的管網(wǎng)年費用折算值.從表中同樣可見,基于自動生成策略的遺傳-模擬退火算法具有最好的優(yōu)化性能,圖10為其優(yōu)化結(jié)果.

    圖9 綿竹市待建供水管網(wǎng)

    優(yōu)化算法年費用折算值/萬元節(jié)點可靠度指標(biāo)最小值微粒群算法276.30.8617遺傳算法272.80.9584遺傳模擬退火算法271.71.0634

    5結(jié)論

    結(jié)合環(huán)形供水管網(wǎng)的判斷方法并將道路信息引入管網(wǎng)抗震優(yōu)化設(shè)計中,建立了供水管網(wǎng)的自動生成策略.同時,將自動生成策略嵌入到3種優(yōu)化算法中,實現(xiàn)遺傳算法、遺傳-模擬退火混合算法和微粒群算法在供水管網(wǎng)的抗震優(yōu)化設(shè)計中的應(yīng)用.利用上述3種算法對2個供水管網(wǎng)進行了實例分析,結(jié)果表明,遺傳-模擬退火算法具有最好的優(yōu)化能力.

    參考文獻:

    [1]李杰. 生命線工程抗震—基礎(chǔ)理論與應(yīng)用[M]. 北京: 科學(xué)出版社, 2005.

    LI Jie. Lifeline earthquake engineering—Basic methods and applications [M]. Beijing: Science Press, 2005.

    [2]Chen L L, Li J, Xu C C. Seismic reliability-based optimization for water supply network[C]//Proc. the 9th International Conference on Computing in Civil and Building Engineering. Taipei: Taiwan University, 2002:793-798.

    [3]李杰,衛(wèi)書麟,劉威. 基于模擬退火算法的供水管網(wǎng)抗震優(yōu)化設(shè)計[J].地震工程與工程振動,2009,29(3):108.

    LI Jie, WEI Shulin, LIU Wei. Seismic reliability optimization of water distribution networks with simulated annealing algorithms[J]. Earthquake Engineering and Engineering Vibration. 2009, 29(3):108.

    [4]邢燕,李杰. 混合優(yōu)化策略在生命線管網(wǎng)拓撲設(shè)計中的應(yīng)用[J].同濟大學(xué)學(xué)報:自然科學(xué)版,2008,36(5):569.

    圖10 基于遺傳-模擬退火算法的綿竹管網(wǎng)優(yōu)化方案

    XING Yan, LI Jie. Hybrid optimization strategy for topology design of lifeline networks[J]. Journal of Tongji University: Natural Science, 2008,36(5):569.

    [5]徐良,劉威,李杰. 基于微粒群算法的供水管網(wǎng)抗震優(yōu)化設(shè)計[J].防災(zāi)減災(zāi)工程學(xué)報,2010,30(3):269.

    XU Liang, LIU Wei, LI Jie. Seismic reliability optimization of water distribution network with particle swarm algorithm[J]. Journal of Disaster Prevention and Mitigation Engineering, 2010,30(3):269.

    [6]李杰,邢燕. 基于可靠度的生命線工程網(wǎng)絡(luò)抗震設(shè)計[J]. 同濟大學(xué)學(xué)報:自然科學(xué)版,2010,38(6):783.

    LI Jie, XING Yan. Reliability-based seismic design of lifeline networks[J]. Journal of Tongji University: Natural Science, 2010,38(6):783.

    [7]陳玲俐,李杰. 城市供水管網(wǎng)系統(tǒng)抗震功能可靠度分析[J]. 工程力學(xué),2004,21(4):45.

    CHEN Lingli, LI Jie. Aseismatic serviceability analysis of water supply network[J]. Engineering Mechanics,2004,21(4):45.

    [8]嚴煦世, 劉遂慶. 給水排水管網(wǎng)系統(tǒng)[M]. 北京: 中國建筑工業(yè)出版社, 2003.

    YAN Xushi, LIU Suiqing. Water and wastewater pipeline network systems[M]. Beijing: China Architecture and Building Press, 2003.

    [9]邵知宇. 給水管網(wǎng)分段線性優(yōu)化模型[D]. 上海: 同濟大學(xué),2001.

    SHAO Zhiyu. Piecewise linear optimization model of water pipeline network[D]. Shanghai: Tongji University, 2001.

    [10]Dijkstra E W. A note on two problems in connexion with graphs[J]. Numerische Mathematlk, 1959,1(1):269.

    [11]Holland J H. Adaptation in natural and artificial systems[M]. ANN Arbor: University of Michigan Press, 1975.

    [12]Kirkpatrick S, Gelatt C D, Vecchi M P. Optimization by simulated annealing[J]. Science, 1983, 220(4958):671.

    [13]Jeong I, Lee J. Adaptive simulated annealing genetic algorithm for system identifictation [J]. Engineering Applications of Artificial Intelligence, 1996, 9(5):523.

    [14]Kennedy J. A new optimizer using particle swarm theory [C]//The Sixth International Symposium on Micro Machine and Human Science. Nagoya: Nagoya Municipal Industrial Research Institute, 1995: 39-43.

    收稿日期:2015-07-07

    基金項目:國家自然科學(xué)基金面上項目(51278380)

    中圖分類號:P315.9

    文獻標(biāo)志碼:A

    Seismic Topology Optimization Algorithms for Water Distribution Networks based on Automatic Generation Strategy

    LIU Wei1,2, YUAN Hongwei3, XU Liang4

    (1. State Key Laboratory for Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China; 2. College of Civil Engineering, Tongji University, Shanghai 200092, China; 3. Beijing Institute of Architectural Design, Beijing 100045, China; 4.Shanghai Nuclear Engineering Research & Design Institute, Shanghai 200233, China)

    Abstract:Taking water distribution network’s annual reduced cost as optimization object and seismic reliability as optimization restriction, a network topology optimization model is established. Then, an automatic generation strategy based on the existed road information is proposed. Combining with ring network generation method and automatic generation strategy, three approaches, a genetic algorithm, a simulated annealing genetic algorithm and a particle swarm algorithm, are employed to solve the seismic optimization model. Two networks are optimized using three algorithms and a comparative study is performed. The comparison indicates that the simulated annealing genetic algorithm performs the best.

    Key words:water distribution network; automatic generation strategy; genetic algorithm; simulated annealing genetic algorithm; particle swarm algorithm

    第一作者: 劉威(1976—),男,副教授,工學(xué)博士,主要研究方向為生命線工程防災(zāi).E-mail:liuw@#edu.cn

    猜你喜歡
    模擬退火算法供水管網(wǎng)遺傳算法
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應(yīng)用
    分區(qū)計量(DMA)在實際工程的應(yīng)用
    綠色科技(2017年8期)2017-05-22 14:41:31
    基于遺傳算法和LS-SVM的財務(wù)危機預(yù)測
    智能傳感器中的算法應(yīng)用
    城鎮(zhèn)供水管網(wǎng)多措并舉漏損控制措施探析
    科技傳播(2016年19期)2016-12-27 16:09:53
    供水管網(wǎng)管道爆管事故發(fā)生原因及控制措施
    城市地下供水管網(wǎng)事故應(yīng)急處置
    科技視界(2016年23期)2016-11-04 13:26:29
    改進的模擬退火算法及其在裝填問題中的應(yīng)用
    基于BP人工神經(jīng)網(wǎng)絡(luò)的離散型車間生產(chǎn)調(diào)度指標(biāo)預(yù)測模型的研究
    科技視界(2016年3期)2016-02-26 09:45:54
    乱系列少妇在线播放| 少妇的逼水好多| 国产一区二区激情短视频| 色噜噜av男人的天堂激情| 97人妻精品一区二区三区麻豆| 久久国产乱子免费精品| 99久久中文字幕三级久久日本| 美女xxoo啪啪120秒动态图| 偷拍熟女少妇极品色| 最近中文字幕高清免费大全6 | 能在线免费观看的黄片| 亚洲av.av天堂| 免费观看精品视频网站| 亚洲久久久久久中文字幕| 欧美日韩亚洲国产一区二区在线观看| 亚洲不卡免费看| av在线老鸭窝| 又爽又黄a免费视频| 欧美国产日韩亚洲一区| 国产精品久久电影中文字幕| 亚洲人成伊人成综合网2020| 久久精品国产99精品国产亚洲性色| 国产成人aa在线观看| 成年免费大片在线观看| 中出人妻视频一区二区| 内射极品少妇av片p| 两个人的视频大全免费| 国产亚洲91精品色在线| 最后的刺客免费高清国语| 精品日产1卡2卡| 国产亚洲精品av在线| 一进一出抽搐gif免费好疼| 99热网站在线观看| 一边摸一边抽搐一进一小说| 亚洲久久久久久中文字幕| 人人妻人人看人人澡| av专区在线播放| 一区福利在线观看| 国产精品亚洲一级av第二区| 免费无遮挡裸体视频| 成人精品一区二区免费| 18+在线观看网站| 亚洲av中文字字幕乱码综合| 亚洲国产日韩欧美精品在线观看| 国产探花极品一区二区| 动漫黄色视频在线观看| 亚洲精品一区av在线观看| 亚洲精品粉嫩美女一区| 三级毛片av免费| 欧美一级a爱片免费观看看| 久久香蕉精品热| 五月伊人婷婷丁香| 免费不卡的大黄色大毛片视频在线观看 | 男女视频在线观看网站免费| 欧美性猛交黑人性爽| 床上黄色一级片| or卡值多少钱| 99国产精品一区二区蜜桃av| 国产精华一区二区三区| 非洲黑人性xxxx精品又粗又长| 校园春色视频在线观看| 国内精品久久久久精免费| 中国美白少妇内射xxxbb| 黄色配什么色好看| 亚洲人成网站在线播放欧美日韩| 精品一区二区三区视频在线观看免费| 此物有八面人人有两片| 五月伊人婷婷丁香| 免费高清视频大片| 观看美女的网站| 麻豆成人午夜福利视频| 国产在线男女| 国产爱豆传媒在线观看| 色综合站精品国产| 国模一区二区三区四区视频| 日本精品一区二区三区蜜桃| 亚洲专区国产一区二区| 免费搜索国产男女视频| 亚洲国产欧美人成| 岛国在线免费视频观看| 男女边吃奶边做爰视频| 亚洲av成人av| a级毛片a级免费在线| 国产一区二区亚洲精品在线观看| 男女之事视频高清在线观看| 一本久久中文字幕| 蜜桃亚洲精品一区二区三区| 成人精品一区二区免费| 国产 一区精品| 在线观看美女被高潮喷水网站| 精品久久久久久久久久久久久| av福利片在线观看| 精品乱码久久久久久99久播| 成人性生交大片免费视频hd| 欧美国产日韩亚洲一区| 99热这里只有是精品50| 成人综合一区亚洲| 天堂av国产一区二区熟女人妻| 欧美日本视频| 不卡视频在线观看欧美| 久久久久免费精品人妻一区二区| 在线观看舔阴道视频| 日韩国内少妇激情av| 最近最新中文字幕大全电影3| 在线免费十八禁| 日本a在线网址| 国产国拍精品亚洲av在线观看| 精品人妻1区二区| 欧美不卡视频在线免费观看| 蜜桃亚洲精品一区二区三区| 91久久精品国产一区二区成人| 国产一区二区在线观看日韩| av.在线天堂| 久久人妻av系列| 999久久久精品免费观看国产| 99精品在免费线老司机午夜| 老司机午夜福利在线观看视频| 精品久久久久久久末码| 亚洲七黄色美女视频| 国产淫片久久久久久久久| 久久精品国产鲁丝片午夜精品 | 淫秽高清视频在线观看| 欧美xxxx性猛交bbbb| 国产在线精品亚洲第一网站| 国产乱人伦免费视频| 亚洲av.av天堂| 狠狠狠狠99中文字幕| 波多野结衣高清无吗| 久久久久久九九精品二区国产| 国产乱人视频| 男人舔女人下体高潮全视频| 国产一区二区亚洲精品在线观看| 欧美xxxx黑人xx丫x性爽| 干丝袜人妻中文字幕| 成人一区二区视频在线观看| av天堂中文字幕网| 精品久久久久久,| 国产一区二区在线av高清观看| 亚洲精品久久国产高清桃花| 午夜精品久久久久久毛片777| 国产一级毛片七仙女欲春2| 亚洲综合色惰| 亚洲黑人精品在线| 国产亚洲精品久久久com| 国产69精品久久久久777片| 成人精品一区二区免费| 久久精品人妻少妇| 亚洲18禁久久av| 精品午夜福利在线看| 成人av在线播放网站| 成人精品一区二区免费| 亚洲电影在线观看av| 最新中文字幕久久久久| 国产一区二区三区在线臀色熟女| 亚洲va在线va天堂va国产| 女人被狂操c到高潮| 久久久久久久久久久丰满 | 哪里可以看免费的av片| 国产精品自产拍在线观看55亚洲| 亚洲最大成人手机在线| 无人区码免费观看不卡| 精品久久久久久成人av| 韩国av在线不卡| 精品一区二区三区视频在线| 丝袜美腿在线中文| 美女xxoo啪啪120秒动态图| 舔av片在线| 国产精品自产拍在线观看55亚洲| 99九九线精品视频在线观看视频| 男人的好看免费观看在线视频| 国产亚洲欧美98| 啦啦啦韩国在线观看视频| 成人综合一区亚洲| 国产美女午夜福利| 成年版毛片免费区| 成人精品一区二区免费| 成年人黄色毛片网站| 免费看美女性在线毛片视频| 春色校园在线视频观看| 亚洲欧美精品综合久久99| 欧美+日韩+精品| 国产精品久久久久久亚洲av鲁大| 91在线精品国自产拍蜜月| 级片在线观看| 日韩欧美精品v在线| 亚洲精品影视一区二区三区av| 黄色欧美视频在线观看| 久久精品国产自在天天线| 春色校园在线视频观看| 日本熟妇午夜| 精品久久国产蜜桃| 亚洲av电影不卡..在线观看| 老熟妇仑乱视频hdxx| 欧美性感艳星| 久久精品国产亚洲av天美| 人妻夜夜爽99麻豆av| 一级毛片久久久久久久久女| 欧美最黄视频在线播放免费| 国产 一区精品| 99riav亚洲国产免费| 麻豆一二三区av精品| 搡老岳熟女国产| 亚洲欧美清纯卡通| 久久精品国产亚洲av涩爱 | 黄色日韩在线| 听说在线观看完整版免费高清| 国产精品久久电影中文字幕| 又粗又爽又猛毛片免费看| 国产高清有码在线观看视频| 大又大粗又爽又黄少妇毛片口| 国产精品福利在线免费观看| 国产真实乱freesex| 可以在线观看毛片的网站| 亚洲性久久影院| 国产日本99.免费观看| 亚洲不卡免费看| 亚洲性久久影院| 老熟妇仑乱视频hdxx| 国产精品人妻久久久影院| 国模一区二区三区四区视频| 亚洲国产精品久久男人天堂| 欧美不卡视频在线免费观看| 国产精品一及| 一本一本综合久久| 乱码一卡2卡4卡精品| 亚洲电影在线观看av| 99精品在免费线老司机午夜| 国产伦在线观看视频一区| 国产免费一级a男人的天堂| 日本五十路高清| 亚洲av日韩精品久久久久久密| 国产人妻一区二区三区在| 日本三级黄在线观看| 搡老熟女国产l中国老女人| 波多野结衣高清无吗| 色综合站精品国产| bbb黄色大片| 日日摸夜夜添夜夜添小说| 国产爱豆传媒在线观看| 99在线人妻在线中文字幕| 99久国产av精品| 悠悠久久av| 欧美在线一区亚洲| 在线播放无遮挡| 一区福利在线观看| 日本一二三区视频观看| 亚洲成a人片在线一区二区| 欧美日韩中文字幕国产精品一区二区三区| 国产精品自产拍在线观看55亚洲| 国产成人aa在线观看| 直男gayav资源| 日韩,欧美,国产一区二区三区 | 好男人在线观看高清免费视频| 最后的刺客免费高清国语| bbb黄色大片| 亚洲av不卡在线观看| 婷婷丁香在线五月| 日本-黄色视频高清免费观看| 九九在线视频观看精品| 国产亚洲精品久久久久久毛片| 免费在线观看影片大全网站| 亚洲国产精品sss在线观看| 自拍偷自拍亚洲精品老妇| 亚洲精品一区av在线观看| 婷婷色综合大香蕉| 波多野结衣高清作品| 超碰av人人做人人爽久久| 在线观看午夜福利视频| 欧美色欧美亚洲另类二区| 日日夜夜操网爽| 国产高清激情床上av| aaaaa片日本免费| 欧美三级亚洲精品| 69人妻影院| 黄色丝袜av网址大全| 午夜福利视频1000在线观看| 99热这里只有是精品在线观看| 99精品在免费线老司机午夜| 神马国产精品三级电影在线观看| 我的女老师完整版在线观看| 性插视频无遮挡在线免费观看| 国产成人aa在线观看| 成人精品一区二区免费| 又粗又爽又猛毛片免费看| 欧美日韩乱码在线| 亚洲美女搞黄在线观看 | 中国美白少妇内射xxxbb| 久久精品国产亚洲av香蕉五月| 亚洲成a人片在线一区二区| 我的老师免费观看完整版| 一个人观看的视频www高清免费观看| 久久久久久久久中文| 色噜噜av男人的天堂激情| 国产精品一区www在线观看 | 日本撒尿小便嘘嘘汇集6| 欧美三级亚洲精品| 亚洲欧美清纯卡通| 欧美区成人在线视频| 亚洲av熟女| 韩国av一区二区三区四区| 亚洲精品粉嫩美女一区| 99精品在免费线老司机午夜| 久久精品影院6| 中文在线观看免费www的网站| 18禁黄网站禁片免费观看直播| 非洲黑人性xxxx精品又粗又长| 国产精品一区二区性色av| 性色avwww在线观看| 成人国产综合亚洲| 两个人视频免费观看高清| 国产亚洲精品av在线| 色哟哟·www| 99热这里只有是精品50| 99久久中文字幕三级久久日本| 亚洲人成网站在线播| 亚洲真实伦在线观看| 99在线视频只有这里精品首页| 国产成人一区二区在线| 成人三级黄色视频| 小说图片视频综合网站| 国产精品人妻久久久影院| 国产成人av教育| 欧美不卡视频在线免费观看| 国产精品,欧美在线| 亚洲人成伊人成综合网2020| 看免费成人av毛片| 丰满人妻一区二区三区视频av| 午夜激情福利司机影院| 成人国产综合亚洲| 免费高清视频大片| 老司机午夜福利在线观看视频| 日本黄色视频三级网站网址| 我要看日韩黄色一级片| 国产v大片淫在线免费观看| 亚洲中文字幕一区二区三区有码在线看| 高清日韩中文字幕在线| 99国产精品一区二区蜜桃av| 久久久久久久久久久丰满 | 亚洲精品一区av在线观看| 中亚洲国语对白在线视频| 精品国内亚洲2022精品成人| 国产爱豆传媒在线观看| 乱人视频在线观看| 非洲黑人性xxxx精品又粗又长| 搡老熟女国产l中国老女人| 日韩欧美一区二区三区在线观看| 能在线免费观看的黄片| 亚洲国产精品合色在线| 波多野结衣高清作品| 国产单亲对白刺激| 亚洲中文字幕日韩| 日韩av在线大香蕉| 国内精品一区二区在线观看| 亚洲av熟女| 成人高潮视频无遮挡免费网站| 男女之事视频高清在线观看| 婷婷亚洲欧美| 中亚洲国语对白在线视频| 一个人看的www免费观看视频| 久久精品综合一区二区三区| 特大巨黑吊av在线直播| 小说图片视频综合网站| 久久久久久国产a免费观看| 男人狂女人下面高潮的视频| 免费观看在线日韩| 欧美zozozo另类| 欧美日本视频| 久久久久久国产a免费观看| 国产在视频线在精品| 搡老岳熟女国产| aaaaa片日本免费| 久久热精品热| 成年人黄色毛片网站| 精品福利观看| 九九热线精品视视频播放| 色视频www国产| 亚洲精品一卡2卡三卡4卡5卡| 国产黄a三级三级三级人| 亚洲精品456在线播放app | 国产一区二区三区在线臀色熟女| 九九爱精品视频在线观看| 亚洲 国产 在线| 国模一区二区三区四区视频| 中国美白少妇内射xxxbb| 88av欧美| 亚洲av不卡在线观看| 精品人妻偷拍中文字幕| 男女边吃奶边做爰视频| 男人舔女人下体高潮全视频| 舔av片在线| 国产在线精品亚洲第一网站| 亚洲国产精品sss在线观看| 最近最新免费中文字幕在线| 久久久精品大字幕| av国产免费在线观看| 国产精品美女特级片免费视频播放器| 在线a可以看的网站| 色精品久久人妻99蜜桃| 亚洲电影在线观看av| 久久久色成人| 色哟哟哟哟哟哟| 国产精品一区www在线观看 | 亚洲成人久久爱视频| 日韩精品青青久久久久久| xxxwww97欧美| 亚洲真实伦在线观看| 成人二区视频| 国产乱人视频| 久久午夜亚洲精品久久| 欧美黑人巨大hd| 久久人妻av系列| 精品人妻视频免费看| 午夜福利视频1000在线观看| 亚洲精华国产精华液的使用体验 | 好男人在线观看高清免费视频| 日韩欧美在线乱码| www.色视频.com| 国产成人aa在线观看| 欧美区成人在线视频| 日本免费a在线| 内射极品少妇av片p| 欧美激情在线99| 在线观看66精品国产| 俺也久久电影网| 麻豆久久精品国产亚洲av| 18禁黄网站禁片午夜丰满| 神马国产精品三级电影在线观看| 91精品国产九色| 久久这里只有精品中国| 十八禁国产超污无遮挡网站| 床上黄色一级片| netflix在线观看网站| 欧美日韩乱码在线| av天堂中文字幕网| 国产淫片久久久久久久久| 亚洲av中文av极速乱 | 国产蜜桃级精品一区二区三区| 久久久久精品国产欧美久久久| 亚洲自偷自拍三级| 国产亚洲精品综合一区在线观看| 搡老妇女老女人老熟妇| 欧美区成人在线视频| 欧美色视频一区免费| 午夜福利高清视频| 最新在线观看一区二区三区| 国产伦精品一区二区三区视频9| 十八禁国产超污无遮挡网站| 久久久久久大精品| 黄色配什么色好看| 久久久久九九精品影院| 国产爱豆传媒在线观看| 一级av片app| 一区二区三区激情视频| 精品久久久久久久久久免费视频| 一级a爱片免费观看的视频| 最近最新免费中文字幕在线| 日韩一区二区视频免费看| 色哟哟·www| 91狼人影院| 免费搜索国产男女视频| 午夜福利在线观看免费完整高清在 | 午夜影院日韩av| 精品日产1卡2卡| 我的女老师完整版在线观看| 赤兔流量卡办理| 日本爱情动作片www.在线观看 | 日日撸夜夜添| 国产伦精品一区二区三区四那| 国内毛片毛片毛片毛片毛片| 国产高清不卡午夜福利| 熟女人妻精品中文字幕| 久久国产乱子免费精品| 九色国产91popny在线| 天堂√8在线中文| 搡老熟女国产l中国老女人| 少妇猛男粗大的猛烈进出视频 | 日韩,欧美,国产一区二区三区 | 国产男靠女视频免费网站| 精品午夜福利视频在线观看一区| 日日撸夜夜添| 人人妻人人看人人澡| 成人午夜高清在线视频| 国产一区二区三区视频了| x7x7x7水蜜桃| 日韩精品有码人妻一区| 欧美一级a爱片免费观看看| 久久久精品欧美日韩精品| 婷婷亚洲欧美| av女优亚洲男人天堂| 国产一区二区三区av在线 | 成熟少妇高潮喷水视频| 五月玫瑰六月丁香| 国产伦一二天堂av在线观看| 久久午夜福利片| 老熟妇仑乱视频hdxx| 91精品国产九色| 99热只有精品国产| 中文资源天堂在线| 在线观看午夜福利视频| 亚洲精品亚洲一区二区| 日韩欧美在线二视频| 亚洲内射少妇av| 日韩国内少妇激情av| 午夜福利视频1000在线观看| 日本与韩国留学比较| 别揉我奶头 嗯啊视频| 精品久久国产蜜桃| 少妇高潮的动态图| 99精品久久久久人妻精品| 少妇人妻一区二区三区视频| 国产在线精品亚洲第一网站| 偷拍熟女少妇极品色| 在线观看免费视频日本深夜| 国产成人av教育| 日韩精品有码人妻一区| av专区在线播放| 乱人视频在线观看| 色5月婷婷丁香| АⅤ资源中文在线天堂| 听说在线观看完整版免费高清| 国产av麻豆久久久久久久| 久久精品影院6| 天堂av国产一区二区熟女人妻| 亚洲av日韩精品久久久久久密| 窝窝影院91人妻| 九色成人免费人妻av| 欧美一级a爱片免费观看看| 在线免费观看的www视频| av在线老鸭窝| 欧美黑人欧美精品刺激| 精品久久久噜噜| 欧美精品国产亚洲| 亚洲欧美日韩高清在线视频| 高清日韩中文字幕在线| 真实男女啪啪啪动态图| 欧美黑人欧美精品刺激| 麻豆国产av国片精品| 可以在线观看毛片的网站| 国产av不卡久久| 国产成人a区在线观看| 日韩欧美三级三区| 最近最新中文字幕大全电影3| www.色视频.com| 国产 一区 欧美 日韩| 国产精品免费一区二区三区在线| 久久久久久久久久久丰满 | a级毛片a级免费在线| 人人妻人人看人人澡| 制服丝袜大香蕉在线| 少妇高潮的动态图| 日本一二三区视频观看| 国产一区二区三区av在线 | 日本黄色片子视频| 搡老岳熟女国产| 精品久久国产蜜桃| 精品99又大又爽又粗少妇毛片 | 亚洲av免费高清在线观看| 国产免费一级a男人的天堂| 精品人妻熟女av久视频| 午夜福利在线在线| 啦啦啦啦在线视频资源| av国产免费在线观看| 人妻久久中文字幕网| 少妇丰满av| 国产在视频线在精品| 五月玫瑰六月丁香| 免费一级毛片在线播放高清视频| 丰满乱子伦码专区| 国产精品野战在线观看| 免费高清视频大片| 亚洲天堂国产精品一区在线| 日日啪夜夜撸| 三级毛片av免费| 欧美又色又爽又黄视频| 成熟少妇高潮喷水视频| 久久人人爽人人爽人人片va| 看黄色毛片网站| 欧美日韩乱码在线| 免费观看在线日韩| 变态另类丝袜制服| 成人一区二区视频在线观看| 亚洲最大成人av| 俄罗斯特黄特色一大片| 赤兔流量卡办理| 91久久精品国产一区二区三区| 18+在线观看网站| 国产精品女同一区二区软件 | 久久婷婷人人爽人人干人人爱| 国产高清三级在线| 亚洲av.av天堂| 精品久久久噜噜| 有码 亚洲区| 麻豆一二三区av精品| 噜噜噜噜噜久久久久久91| 午夜日韩欧美国产| 我要搜黄色片| 网址你懂的国产日韩在线| 亚洲欧美清纯卡通| 在线观看av片永久免费下载| 国产69精品久久久久777片| 国产成年人精品一区二区| 婷婷亚洲欧美| 在线看三级毛片| 美女高潮喷水抽搐中文字幕| 不卡视频在线观看欧美| bbb黄色大片| 免费电影在线观看免费观看| 人妻制服诱惑在线中文字幕| 美女 人体艺术 gogo| 美女高潮喷水抽搐中文字幕| 伦理电影大哥的女人| xxxwww97欧美| 亚洲熟妇熟女久久| 国产老妇女一区| 精品日产1卡2卡|