• <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
    小说图片视频综合网站| 老司机深夜福利视频在线观看| 夜夜夜夜夜久久久久| 俄罗斯特黄特色一大片| 日本五十路高清| 在线播放无遮挡| 美女高潮的动态| 99国产精品一区二区蜜桃av| 色噜噜av男人的天堂激情| 成人国产综合亚洲| 国产成人啪精品午夜网站| 亚洲aⅴ乱码一区二区在线播放| 日本五十路高清| 黄色一级大片看看| 日本在线视频免费播放| 欧美激情在线99| 精品一区二区免费观看| 亚洲人成网站在线播放欧美日韩| 国产精品一区二区免费欧美| av女优亚洲男人天堂| 成人特级av手机在线观看| 2021天堂中文幕一二区在线观| 激情在线观看视频在线高清| 又爽又黄a免费视频| 亚洲国产色片| 深爱激情五月婷婷| 日本 av在线| 永久网站在线| 欧美高清成人免费视频www| 中出人妻视频一区二区| 国产成人aa在线观看| 国产麻豆成人av免费视频| 51午夜福利影视在线观看| 国产免费av片在线观看野外av| 欧美日韩黄片免| 9191精品国产免费久久| 午夜福利成人在线免费观看| 一边摸一边抽搐一进一小说| 久久人人爽人人爽人人片va | 91麻豆精品激情在线观看国产| 国产精品一区二区性色av| 在线国产一区二区在线| 日本黄色片子视频| 欧美黑人巨大hd| 亚洲熟妇熟女久久| 自拍偷自拍亚洲精品老妇| 国产精品国产高清国产av| 久久国产精品人妻蜜桃| 非洲黑人性xxxx精品又粗又长| av视频在线观看入口| 国产精品一及| 国产高清视频在线播放一区| 欧美高清成人免费视频www| 色尼玛亚洲综合影院| 在线观看舔阴道视频| 中文字幕av在线有码专区| 婷婷亚洲欧美| 91午夜精品亚洲一区二区三区 | 香蕉av资源在线| a级毛片免费高清观看在线播放| 九色成人免费人妻av| 99国产综合亚洲精品| 亚洲最大成人手机在线| 精品99又大又爽又粗少妇毛片 | 亚洲午夜理论影院| 久久精品国产亚洲av天美| 亚洲最大成人手机在线| 国内久久婷婷六月综合欲色啪| 欧美乱妇无乱码| 在线观看美女被高潮喷水网站 | 一区二区三区四区激情视频 | 亚洲国产欧洲综合997久久,| 国产午夜精品久久久久久一区二区三区 | 99热6这里只有精品| 午夜福利在线观看吧| 嫩草影院精品99| 亚洲激情在线av| 久久欧美精品欧美久久欧美| 日本免费一区二区三区高清不卡| 99精品在免费线老司机午夜| 九色国产91popny在线| 国产 一区 欧美 日韩| 欧美精品啪啪一区二区三区| 亚洲第一电影网av| 亚洲一区高清亚洲精品| 美女 人体艺术 gogo| 黄色丝袜av网址大全| 99久久精品一区二区三区| 亚洲成av人片在线播放无| 久久国产精品人妻蜜桃| 蜜桃亚洲精品一区二区三区| 九色国产91popny在线| 日本三级黄在线观看| 老女人水多毛片| 欧美最黄视频在线播放免费| 脱女人内裤的视频| 午夜福利成人在线免费观看| 极品教师在线视频| 简卡轻食公司| 国产麻豆成人av免费视频| 十八禁人妻一区二区| 国内毛片毛片毛片毛片毛片| 国产精品,欧美在线| 国产三级在线视频| 女同久久另类99精品国产91| 色5月婷婷丁香| 十八禁国产超污无遮挡网站| 美女高潮的动态| 久久久色成人| 日韩av在线大香蕉| 色综合站精品国产| 国产私拍福利视频在线观看| 3wmmmm亚洲av在线观看| 国内少妇人妻偷人精品xxx网站| 在线天堂最新版资源| 国内精品一区二区在线观看| 真人做人爱边吃奶动态| 高潮久久久久久久久久久不卡| 色综合欧美亚洲国产小说| 老女人水多毛片| 99久久成人亚洲精品观看| 成人亚洲精品av一区二区| 精品熟女少妇八av免费久了| 少妇的逼好多水| 日韩精品中文字幕看吧| 在线播放无遮挡| 啦啦啦韩国在线观看视频| 亚洲人成网站在线播| 色视频www国产| 亚洲人成电影免费在线| 国产精华一区二区三区| 国产精品久久视频播放| 观看美女的网站| 国产精品精品国产色婷婷| 长腿黑丝高跟| 国产国拍精品亚洲av在线观看| 高清日韩中文字幕在线| www.色视频.com| 最新中文字幕久久久久| 欧美bdsm另类| 久久国产乱子免费精品| 人人妻人人澡欧美一区二区| 97碰自拍视频| 97碰自拍视频| 又粗又爽又猛毛片免费看| 国产免费一级a男人的天堂| 在线天堂最新版资源| 在线a可以看的网站| 欧美不卡视频在线免费观看| 少妇熟女aⅴ在线视频| 一本综合久久免费| 人人妻人人看人人澡| 日日夜夜操网爽| 精品一区二区三区av网在线观看| 国产在视频线在精品| 91麻豆av在线| 桃色一区二区三区在线观看| av女优亚洲男人天堂| 亚洲欧美清纯卡通| 天天躁日日操中文字幕| 国产精品女同一区二区软件 | 亚洲精品成人久久久久久| 日本与韩国留学比较| av中文乱码字幕在线| 免费av毛片视频| 亚洲国产精品成人综合色| 在线免费观看不下载黄p国产 | 中亚洲国语对白在线视频| 亚洲国产精品成人综合色| 69av精品久久久久久| 91久久精品国产一区二区成人| 天天躁日日操中文字幕| 国模一区二区三区四区视频| 亚洲五月婷婷丁香| 精品欧美国产一区二区三| 亚洲成人久久爱视频| 狂野欧美白嫩少妇大欣赏| 国产麻豆成人av免费视频| 18禁黄网站禁片午夜丰满| 国产成人aa在线观看| 日韩有码中文字幕| 麻豆成人午夜福利视频| 欧美日韩瑟瑟在线播放| 成人永久免费在线观看视频| 亚洲国产日韩欧美精品在线观看| 色综合婷婷激情| 国产欧美日韩精品一区二区| 怎么达到女性高潮| 精品一区二区三区av网在线观看| 精品人妻一区二区三区麻豆 | 亚洲无线观看免费| 中文资源天堂在线| 18+在线观看网站| 少妇的逼水好多| 免费一级毛片在线播放高清视频| 可以在线观看毛片的网站| 国产激情偷乱视频一区二区| 国产激情偷乱视频一区二区| 99久久99久久久精品蜜桃| 亚洲一区二区三区不卡视频| 亚洲欧美日韩东京热| 国产毛片a区久久久久| 亚洲av二区三区四区| 在线天堂最新版资源| 亚洲专区中文字幕在线| 亚洲乱码一区二区免费版| 亚洲精品日韩av片在线观看| 国内揄拍国产精品人妻在线| 美女黄网站色视频| 久久人人爽人人爽人人片va | 亚洲精华国产精华精| 欧美成人一区二区免费高清观看| 亚洲精品久久国产高清桃花| 国产亚洲av嫩草精品影院| 国产淫片久久久久久久久 | 直男gayav资源| 又爽又黄a免费视频| 精品久久久久久久久久久久久| 尤物成人国产欧美一区二区三区| 国产在视频线在精品| 日本一本二区三区精品| 国产色爽女视频免费观看| 久久久精品大字幕| 真人做人爱边吃奶动态| 搡老熟女国产l中国老女人| 成人永久免费在线观看视频| 老熟妇乱子伦视频在线观看| 欧美区成人在线视频| 老熟妇仑乱视频hdxx| 色在线成人网| 国产欧美日韩一区二区三| 国产国拍精品亚洲av在线观看| 欧美性感艳星| 国产精品久久久久久人妻精品电影| 在线观看午夜福利视频| 亚洲人成电影免费在线| 欧美中文日本在线观看视频| 久久精品国产亚洲av天美| 成人一区二区视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 欧美最新免费一区二区三区 | 久久国产乱子伦精品免费另类| 午夜a级毛片| 五月玫瑰六月丁香| 久久久精品欧美日韩精品| 日日干狠狠操夜夜爽| 午夜福利在线观看吧| 午夜精品在线福利| 亚洲成人久久爱视频| 女人被狂操c到高潮| 精品久久久久久成人av| 日韩欧美国产在线观看| 日日夜夜操网爽| 色5月婷婷丁香| 又紧又爽又黄一区二区| 悠悠久久av| 日韩成人在线观看一区二区三区| avwww免费| 国产欧美日韩精品亚洲av| 欧美性感艳星| 亚洲五月婷婷丁香| 成年免费大片在线观看| www.熟女人妻精品国产| 色综合婷婷激情| 精品不卡国产一区二区三区| 国产三级中文精品| 国产蜜桃级精品一区二区三区| 久久久久久大精品| 久久午夜福利片| 无人区码免费观看不卡| 免费黄网站久久成人精品 | 一区福利在线观看| 欧美精品国产亚洲| 天堂网av新在线| 亚洲av第一区精品v没综合| 女生性感内裤真人,穿戴方法视频| 国产av不卡久久| 我要看日韩黄色一级片| 亚洲性夜色夜夜综合| 欧美性感艳星| 亚洲av熟女| 欧美成狂野欧美在线观看| 国产69精品久久久久777片| 国产精品亚洲av一区麻豆| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产黄色小视频在线观看| 一区二区三区免费毛片| 久久久国产成人免费| 久99久视频精品免费| 麻豆av噜噜一区二区三区| 蜜桃久久精品国产亚洲av| 亚洲三级黄色毛片| 亚洲国产日韩欧美精品在线观看| 又黄又爽又免费观看的视频| 成人特级黄色片久久久久久久| 午夜福利在线观看免费完整高清在 | 午夜激情欧美在线| 欧美成人免费av一区二区三区| 国产色爽女视频免费观看| netflix在线观看网站| 一个人看的www免费观看视频| 欧美一区二区精品小视频在线| 亚洲片人在线观看| 国产一区二区在线av高清观看| 欧美日韩乱码在线| 亚洲国产精品成人综合色| 亚洲最大成人手机在线| 国产精品伦人一区二区| 真实男女啪啪啪动态图| 18+在线观看网站| 亚洲经典国产精华液单 | 亚洲五月天丁香| 真人一进一出gif抽搐免费| 亚洲欧美精品综合久久99| 啦啦啦观看免费观看视频高清| 757午夜福利合集在线观看| 中文亚洲av片在线观看爽| 国产男靠女视频免费网站| 欧美精品国产亚洲| 有码 亚洲区| 久久久久九九精品影院| 久久久久精品国产欧美久久久| 成人欧美大片| 亚洲av二区三区四区| 国产男靠女视频免费网站| 美女高潮喷水抽搐中文字幕| 一区二区三区四区激情视频 | 国产真实乱freesex| 色尼玛亚洲综合影院| 国产熟女xx| 每晚都被弄得嗷嗷叫到高潮| 亚洲综合色惰| 51午夜福利影视在线观看| 午夜激情欧美在线| 久久久久久久亚洲中文字幕 | 淫秽高清视频在线观看| 国产欧美日韩一区二区精品| 欧美国产日韩亚洲一区| 搞女人的毛片| 成人国产一区最新在线观看| 亚洲av成人av| 日本一二三区视频观看| 色尼玛亚洲综合影院| 亚洲欧美日韩东京热| 最好的美女福利视频网| 国产精品日韩av在线免费观看| 国产亚洲精品久久久com| 丰满人妻熟妇乱又伦精品不卡| 国产探花极品一区二区| 男人舔奶头视频| 欧美精品国产亚洲| 男插女下体视频免费在线播放| 校园春色视频在线观看| 一个人看视频在线观看www免费| 日韩欧美国产一区二区入口| 亚洲av成人av| 免费观看的影片在线观看| 国产精品日韩av在线免费观看| 国产亚洲欧美在线一区二区| 国产亚洲av嫩草精品影院| 色尼玛亚洲综合影院| 亚洲va日本ⅴa欧美va伊人久久| 成人国产一区最新在线观看| 深夜精品福利| 国产精品久久视频播放| 欧美3d第一页| 午夜日韩欧美国产| 日本免费一区二区三区高清不卡| 国产蜜桃级精品一区二区三区| 精品免费久久久久久久清纯| 老司机午夜十八禁免费视频| 成人欧美大片| 国产一区二区三区视频了| 最新中文字幕久久久久| 观看美女的网站| 97超级碰碰碰精品色视频在线观看| 欧美最黄视频在线播放免费| 日本黄色视频三级网站网址| 黄色女人牲交| 欧美黑人巨大hd| 午夜福利成人在线免费观看| 搞女人的毛片| 2021天堂中文幕一二区在线观| 亚洲无线在线观看| 简卡轻食公司| 99热精品在线国产| 麻豆久久精品国产亚洲av| 97超视频在线观看视频| 中国美女看黄片| 直男gayav资源| .国产精品久久| 国产蜜桃级精品一区二区三区| 免费人成视频x8x8入口观看| 极品教师在线视频| 91在线观看av| 午夜福利高清视频| 久久精品影院6| 久9热在线精品视频| 日本一本二区三区精品| 亚洲成人免费电影在线观看| 精品乱码久久久久久99久播| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美精品综合久久99| 熟妇人妻久久中文字幕3abv| 无遮挡黄片免费观看| 欧美日韩乱码在线| 国产视频内射| 久久久久久久久中文| 99久久久亚洲精品蜜臀av| 国产精品一区二区三区四区久久| 免费观看的影片在线观看| 真实男女啪啪啪动态图| 亚洲第一区二区三区不卡| 亚洲精品色激情综合| 日韩亚洲欧美综合| 精品国内亚洲2022精品成人| 国产免费一级a男人的天堂| 色哟哟·www| x7x7x7水蜜桃| 国产午夜精品论理片| 丁香欧美五月| 狠狠狠狠99中文字幕| 老司机福利观看| 又紧又爽又黄一区二区| 国产高清视频在线观看网站| 日韩有码中文字幕| 中文资源天堂在线| 久久精品久久久久久噜噜老黄 | 国产欧美日韩精品一区二区| 神马国产精品三级电影在线观看| 午夜亚洲福利在线播放| 免费黄网站久久成人精品 | 超碰av人人做人人爽久久| 欧美丝袜亚洲另类 | 淫妇啪啪啪对白视频| 日本免费一区二区三区高清不卡| 综合色av麻豆| 久久久国产成人精品二区| 黄色一级大片看看| 精品一区二区免费观看| 男人的好看免费观看在线视频| 亚洲欧美日韩高清在线视频| 人妻夜夜爽99麻豆av| 成人永久免费在线观看视频| 亚洲乱码一区二区免费版| 69av精品久久久久久| 中文字幕人妻熟人妻熟丝袜美| 亚洲中文字幕日韩| 亚洲内射少妇av| 夜夜夜夜夜久久久久| 国产精品国产高清国产av| 他把我摸到了高潮在线观看| 国产精品嫩草影院av在线观看 | 成人午夜高清在线视频| 国产在视频线在精品| 一区二区三区高清视频在线| 99国产综合亚洲精品| 亚洲国产精品久久男人天堂| 亚洲电影在线观看av| a级毛片免费高清观看在线播放| 成年女人毛片免费观看观看9| 麻豆成人av在线观看| 国产亚洲欧美在线一区二区| 亚洲第一欧美日韩一区二区三区| 一级黄色大片毛片| 国产一区二区亚洲精品在线观看| 最新在线观看一区二区三区| 国产精品1区2区在线观看.| 久久精品影院6| 成年女人毛片免费观看观看9| 熟女电影av网| 国产精品爽爽va在线观看网站| 亚洲人成网站在线播| 亚洲成人久久性| 欧美日韩乱码在线| av在线老鸭窝| 精品一区二区三区av网在线观看| 男人舔奶头视频| 亚洲欧美日韩卡通动漫| 午夜精品一区二区三区免费看| 午夜福利18| 美女cb高潮喷水在线观看| 91狼人影院| 国产成年人精品一区二区| 欧美丝袜亚洲另类 | 亚洲中文日韩欧美视频| 亚洲成人中文字幕在线播放| 男人和女人高潮做爰伦理| 久久久久免费精品人妻一区二区| 他把我摸到了高潮在线观看| 国产中年淑女户外野战色| 国产精品1区2区在线观看.| 国产亚洲欧美在线一区二区| 精品久久久久久成人av| 日韩免费av在线播放| 久久亚洲真实| 国产精品野战在线观看| 午夜a级毛片| 亚洲人与动物交配视频| 国产黄片美女视频| 亚洲国产欧美人成| 日本三级黄在线观看| 亚洲精品在线观看二区| 十八禁人妻一区二区| 欧美+日韩+精品| 波多野结衣高清作品| 成熟少妇高潮喷水视频| 国产高清视频在线播放一区| 男女视频在线观看网站免费| 午夜免费男女啪啪视频观看 | 18禁黄网站禁片午夜丰满| 国产白丝娇喘喷水9色精品| 成年女人永久免费观看视频| 午夜福利成人在线免费观看| 中文字幕av成人在线电影| 嫩草影院入口| 欧美日本视频| 亚洲激情在线av| 性欧美人与动物交配| 人妻夜夜爽99麻豆av| 亚洲七黄色美女视频| 日日夜夜操网爽| 亚洲美女视频黄频| 免费在线观看亚洲国产| 亚洲第一区二区三区不卡| 久久6这里有精品| 久久久久久久久中文| 美女高潮喷水抽搐中文字幕| 成熟少妇高潮喷水视频| 好男人电影高清在线观看| 桃红色精品国产亚洲av| 黄色配什么色好看| 国产一区二区亚洲精品在线观看| www.熟女人妻精品国产| 亚洲成人精品中文字幕电影| 国产真实乱freesex| 精品福利观看| av在线天堂中文字幕| 国产探花在线观看一区二区| 国产免费av片在线观看野外av| xxxwww97欧美| 精品一区二区三区视频在线观看免费| 极品教师在线免费播放| av在线老鸭窝| 亚洲片人在线观看| 国产精品亚洲av一区麻豆| 两个人的视频大全免费| 欧美成人一区二区免费高清观看| 色5月婷婷丁香| 国产成人福利小说| 高清在线国产一区| 美女大奶头视频| 18禁黄网站禁片午夜丰满| 夜夜看夜夜爽夜夜摸| 亚洲aⅴ乱码一区二区在线播放| 午夜视频国产福利| 色吧在线观看| 老司机福利观看| 亚洲乱码一区二区免费版| 简卡轻食公司| 亚洲中文字幕一区二区三区有码在线看| 欧美成人a在线观看| 天天一区二区日本电影三级| 乱人视频在线观看| 真人做人爱边吃奶动态| 男女床上黄色一级片免费看| 一本精品99久久精品77| 亚洲美女视频黄频| 一级av片app| 亚洲国产精品sss在线观看| 一进一出抽搐动态| 国产高清视频在线播放一区| 99riav亚洲国产免费| 又黄又爽又免费观看的视频| 免费观看人在逋| 成人精品一区二区免费| 九色成人免费人妻av| 又黄又爽又刺激的免费视频.| 久久婷婷人人爽人人干人人爱| 特大巨黑吊av在线直播| 亚洲美女黄片视频| 两人在一起打扑克的视频| 国产精品久久久久久亚洲av鲁大| 熟女人妻精品中文字幕| 麻豆久久精品国产亚洲av| 免费观看的影片在线观看| 亚洲人成网站在线播放欧美日韩| 久久午夜福利片| 婷婷精品国产亚洲av| 亚洲aⅴ乱码一区二区在线播放| 天天一区二区日本电影三级| 日韩欧美免费精品| 99久久99久久久精品蜜桃| 麻豆国产97在线/欧美| 免费人成在线观看视频色| 国产欧美日韩精品亚洲av| 午夜精品一区二区三区免费看| 亚洲av美国av| 国产高清三级在线| 久久久久久久久久黄片| 国产大屁股一区二区在线视频| 久久久久久久久大av| www.熟女人妻精品国产| 他把我摸到了高潮在线观看| www.熟女人妻精品国产| 黄色女人牲交| 黄色丝袜av网址大全| 香蕉av资源在线| 天堂影院成人在线观看| 香蕉av资源在线| 观看免费一级毛片| 国产精品免费一区二区三区在线| 国产aⅴ精品一区二区三区波| 美女黄网站色视频| 真实男女啪啪啪动态图|