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

    考慮可靠性的交通網(wǎng)絡(luò)應(yīng)急資源布局優(yōu)化

    2018-03-07 01:48:49秦進(jìn)葉勇申純燕張威洪健
    關(guān)鍵詞:儲(chǔ)備染色體可靠性

    秦進(jìn),葉勇,申純燕,張威,洪健

    (中南大學(xué) 交通運(yùn)輸工程學(xué)院,湖南 長(zhǎng)沙 410075)

    近年來(lái),世界范圍內(nèi)頻發(fā)的各類自然災(zāi)害事件和人為災(zāi)害事故,造成了大量的人員傷亡及經(jīng)濟(jì)損失,其中我國(guó)是世界上自然災(zāi)害最嚴(yán)重的國(guó)家之一。應(yīng)急資源優(yōu)化布局是防災(zāi)減災(zāi)的重要策略和措施,應(yīng)急資源在災(zāi)后能否提供快速有效的服務(wù),將對(duì)實(shí)施及時(shí)有效的救援具有重大意義。而交通網(wǎng)絡(luò)是實(shí)施災(zāi)后救援活動(dòng)的重要物質(zhì)基礎(chǔ),對(duì)災(zāi)后應(yīng)急救援能否達(dá)到高效化和災(zāi)區(qū)社會(huì)經(jīng)濟(jì)活動(dòng)能否盡快恢復(fù)正常化,交通網(wǎng)絡(luò)都是其主導(dǎo)性的影響因素之一。由于應(yīng)急資源布局對(duì)應(yīng)急救援的重要性以及求解的困難性,其長(zhǎng)期以來(lái)都受到國(guó)內(nèi)外相關(guān)學(xué)者的關(guān)注。Masood等[1]建立了消防站選址的多目標(biāo)數(shù)學(xué)模型。付德強(qiáng)等[2]建立了多目標(biāo)選址決策模型,并結(jié)合模型特點(diǎn)采用非支配排序多目標(biāo)遺傳算法。丁雪楓等[3]考慮應(yīng)急設(shè)施選址的公平性、效率性、成本等因素,建立多目標(biāo)規(guī)劃決策模型,提出了基于模擬植物生長(zhǎng)思路的求解方法。楊金順等[4]從建立了公路網(wǎng)應(yīng)急救援點(diǎn)的多目標(biāo)選址優(yōu)化模型。王旭坪等[5]建立了應(yīng)急響應(yīng)時(shí)間的感知滿意度函數(shù)以衡量災(zāi)民對(duì)救援響應(yīng)時(shí)間的滿意程度,構(gòu)建了一個(gè)突發(fā)事件發(fā)生后的初始階段應(yīng)急物資分配問(wèn)題。張玲等[6]考慮災(zāi)害發(fā)生時(shí)需求、運(yùn)輸代價(jià)等不確定因素,建立基于情景的應(yīng)急資源布局兩階段數(shù)學(xué)規(guī)劃模型來(lái)解決應(yīng)急資源布局問(wèn)題。Ceyhun等[7]在綜合考慮了應(yīng)急救援過(guò)程中的時(shí)間限制和資金限制的基礎(chǔ)上,建立了一個(gè)基于覆蓋集的多目標(biāo)應(yīng)急車輛選址模型。艾云飛等[8]基于集合覆蓋理論和引力模型,建立水上應(yīng)急物資儲(chǔ)備庫(kù)選址-分配優(yōu)化模型。俞武揚(yáng)[9]以最大化應(yīng)急服務(wù)質(zhì)量與數(shù)量的綜合指標(biāo)為目標(biāo),研究了應(yīng)急點(diǎn)具有最低服務(wù)質(zhì)量和數(shù)量要求下的應(yīng)急設(shè)施選址問(wèn)題。盡管目前已有較多研究人員對(duì)應(yīng)急資源布局問(wèn)題開(kāi)展了具體研究,但是多是從人道救援的角度出發(fā)研究應(yīng)急物資布局問(wèn)題。交通網(wǎng)絡(luò)作為災(zāi)后人道救援的最基本物質(zhì)保障,其應(yīng)急預(yù)防護(hù)方面的研究卻還較少,尤其是既有研究中很少考慮對(duì)應(yīng)急服務(wù)的可靠性要求。本文正是針對(duì)交通網(wǎng)絡(luò)基礎(chǔ)設(shè)施的應(yīng)急預(yù)防護(hù)工作,在考慮服務(wù)可靠性和資源多樣性等的前提下,建立優(yōu)化模型進(jìn)行深入的分析和研究。

    1 變量定義與假設(shè)

    本文所研究的問(wèn)題描述如下:針對(duì)災(zāi)害風(fēng)險(xiǎn)區(qū)域內(nèi)的交通網(wǎng)絡(luò),災(zāi)后其基礎(chǔ)設(shè)施有多個(gè)可能的應(yīng)急維護(hù)點(diǎn)(下稱應(yīng)急點(diǎn)),每個(gè)應(yīng)急點(diǎn)均對(duì)多種應(yīng)急設(shè)備資源提出需求?,F(xiàn)需在多個(gè)備選點(diǎn)上,選擇建設(shè)一定數(shù)量的應(yīng)急設(shè)備資源儲(chǔ)備點(diǎn)(下稱應(yīng)急儲(chǔ)備點(diǎn)),并合理配置各類設(shè)備資源,使得其在規(guī)劃期內(nèi)能以最小成本為所有應(yīng)急點(diǎn)提供具有一定可靠度的應(yīng)急救援服務(wù)。

    為簡(jiǎn)化問(wèn)題分析,首先做如下合理假設(shè):

    1) 應(yīng)急儲(chǔ)備點(diǎn)的備選位置及其相應(yīng)的建設(shè)成本等均已知;

    2) 應(yīng)急點(diǎn)對(duì)應(yīng)急設(shè)備資源的需求,可以由不同應(yīng)急儲(chǔ)備點(diǎn)進(jìn)行提供,即任一應(yīng)急點(diǎn)均可被多個(gè)應(yīng)急儲(chǔ)備點(diǎn)覆蓋并提供服務(wù);

    3) 災(zāi)后各應(yīng)急點(diǎn)對(duì)應(yīng)急設(shè)備需求的類型和數(shù)量,根據(jù)其基礎(chǔ)設(shè)施類型、重要程度和所在位置等情況的不同而不同;

    4) 災(zāi)后各類應(yīng)急設(shè)備都須在規(guī)定時(shí)間內(nèi)到達(dá)應(yīng)急點(diǎn)??紤]到災(zāi)后交通組織的特殊性,可將此時(shí)間約束轉(zhuǎn)化為距離約束,即應(yīng)急點(diǎn)與為其服務(wù)的應(yīng)急儲(chǔ)備點(diǎn)之間的距離,不應(yīng)大于給定的最大值。

    為便于描述問(wèn)題,定義變量如下:I表示災(zāi)后交通網(wǎng)絡(luò)中應(yīng)急點(diǎn)的集合,i∈I;J表示交通網(wǎng)絡(luò)中應(yīng)急儲(chǔ)備點(diǎn)的集合,j∈J;M 表示應(yīng)急設(shè)備資源類別的集合,主要為各類應(yīng)急救援車輛設(shè)備,包括清障車、托盤車、起吊車和消防車等,m∈M;為災(zāi)害發(fā)生后應(yīng)急點(diǎn) i對(duì)應(yīng)急設(shè)備 m的最少需求數(shù)量,與應(yīng)急點(diǎn)的類型、重要程度和所在位置等有關(guān);為應(yīng)急儲(chǔ)備點(diǎn)j上允許儲(chǔ)備的應(yīng)急設(shè)備m的最大數(shù)量;為應(yīng)急儲(chǔ)備點(diǎn) j上儲(chǔ)備的應(yīng)急設(shè)備 m的實(shí)際數(shù)量;Uj為應(yīng)急儲(chǔ)備點(diǎn)j上的儲(chǔ)備能力;μm為應(yīng)急設(shè)備m的儲(chǔ)備能力占用率;dij為應(yīng)急點(diǎn)i至應(yīng)急儲(chǔ)備點(diǎn)j之間的實(shí)際距離;Di為應(yīng)急點(diǎn)i至能為其提供服務(wù)的應(yīng)急儲(chǔ)備點(diǎn)之間的最大距離,該值與應(yīng)急點(diǎn)所處區(qū)域的交通道路條件有關(guān);Si為與應(yīng)急點(diǎn)i之間的距離不大于Di的所有應(yīng)急儲(chǔ)備點(diǎn)的集合,即能有效覆蓋應(yīng)急點(diǎn) 的所有應(yīng)急儲(chǔ)備點(diǎn)的集合,為受災(zāi)區(qū)域內(nèi)應(yīng)急儲(chǔ)備點(diǎn)上的設(shè)備m處于繁忙的概率,該概率值可以根據(jù)文獻(xiàn)[10]中所提出的公式進(jìn)行估算;fj為選擇在備選點(diǎn)j建設(shè)應(yīng)急儲(chǔ)備點(diǎn)時(shí)的建設(shè)成本,千元; rjm為在應(yīng)急存儲(chǔ)點(diǎn) j處的應(yīng)急設(shè)備 m的單位購(gòu)置成本(千元),同類設(shè)備考慮到市場(chǎng)差異和運(yùn)輸條件的不同,購(gòu)置成本也有所不同; vmj為在應(yīng)急存儲(chǔ)點(diǎn)j處的設(shè)備m的單位庫(kù)存成本(千元/年),包括人員工資和維護(hù)保養(yǎng)費(fèi)用等;W為規(guī)劃期內(nèi)允許建設(shè)的應(yīng)急儲(chǔ)備點(diǎn)的最大數(shù)量;α為應(yīng)急儲(chǔ)備點(diǎn)的服務(wù)水平,即其能為應(yīng)急點(diǎn)提供服務(wù)的可靠性水平,0<α<1;ε為銀行利率水平;H為系統(tǒng)規(guī)劃期時(shí)間長(zhǎng)度,a;

    同時(shí)定義決策變量如下:

    2 優(yōu)化模型的建立

    2.1 應(yīng)急服務(wù)的可靠性約束

    以應(yīng)急點(diǎn)i∈I對(duì)設(shè)備m∈M的需求為例,應(yīng)急服務(wù)可靠性約束可表述為:針對(duì)應(yīng)急點(diǎn) i,系統(tǒng)必須保證不低于事先給定的服務(wù)水平α為其提供不少于的應(yīng)急設(shè)備進(jìn)行服務(wù),即覆蓋應(yīng)急點(diǎn)i的所有應(yīng)急儲(chǔ)備點(diǎn),能以概率 α保證至少有數(shù)量為的設(shè)備m為該應(yīng)急點(diǎn)提供服務(wù)。

    對(duì)于應(yīng)急點(diǎn) i∈I,記覆蓋其的所有應(yīng)急儲(chǔ)備點(diǎn)上的設(shè)備m的總數(shù)為,則以概率α保證至少有數(shù)量為的設(shè)備m能即時(shí)為應(yīng)急點(diǎn)i提供所需服務(wù),可理解為當(dāng)災(zāi)后該應(yīng)急點(diǎn)向相鄰應(yīng)急儲(chǔ)備點(diǎn)提出服務(wù)請(qǐng)求時(shí),上述所有設(shè)備中處于空閑狀態(tài)的設(shè)備數(shù)量的概率應(yīng)不小于α。

    根據(jù)前面的變量定義,災(zāi)后區(qū)域內(nèi)應(yīng)急設(shè)備m的繁忙概率為 pm,即其處于空閑狀態(tài)的概率為(1-pm),因此對(duì)于應(yīng)急點(diǎn)i,其能及時(shí)得到數(shù)量不少于的應(yīng)急設(shè)備進(jìn)行支援的概率為:

    由于式(1)中,pm和都是事先給定的常數(shù),因此概率實(shí)際上可以視為是設(shè)備總數(shù)的函數(shù),記為=f(Nm)。i

    由此,應(yīng)急服務(wù)的可靠性約束可描述為:

    定理:若 0 <p<1,q為給定的正整數(shù),正整數(shù)n≥q,記:

    則f(n)為關(guān)于n的單調(diào)遞增函數(shù)。

    證明:

    令g(n)=1-f(n),可得:

    由此,可知:

    根據(jù)組合計(jì)算公式,可知:

    根據(jù)式(8),可很容易的發(fā)現(xiàn)有:

    此即證明g(n)是關(guān)于n的單調(diào)減函數(shù)。

    考慮到f(n)=1-g(n),因此f(n)為關(guān)于n的單調(diào)增函數(shù)。

    證畢。

    由此,根據(jù)前面的變量定義,應(yīng)急服務(wù)的可靠性約束可以描述為:

    2.2 成本分析

    系統(tǒng)總成本涉及3個(gè)方面,即儲(chǔ)備點(diǎn)的建設(shè)成本、設(shè)備購(gòu)置成本和設(shè)備庫(kù)存成本。對(duì)于應(yīng)急儲(chǔ)備點(diǎn) j∈J,其建設(shè)成本 fj,顯然與備選點(diǎn) j所處位置的土地價(jià)格和建設(shè)規(guī)模等有關(guān),這里假設(shè)為已知;其設(shè)備購(gòu)置成本rj和設(shè)備庫(kù)存成本vj,根據(jù)變量定義可知:

    在規(guī)劃期內(nèi)建設(shè)成本和設(shè)備購(gòu)置成本是一次性投入,而設(shè)備庫(kù)存成本則一般是按年支出。為在時(shí)間單位上保持一致性,可利用年金現(xiàn)值公式將相關(guān)一次性支出成本分?jǐn)偟揭?guī)劃期的每一年。年金現(xiàn)值轉(zhuǎn)化系數(shù)γ的計(jì)算如下:

    根據(jù)以上分析,應(yīng)急儲(chǔ)備點(diǎn)j上的總成本cj,可表示為

    則系統(tǒng)總成本C可以表示為:

    2.3 優(yōu)化模型

    根據(jù)以上分析,考慮服務(wù)可靠性的交通網(wǎng)絡(luò)應(yīng)急資源布局優(yōu)化模型,可建立如下:

    上述優(yōu)化模型中,目標(biāo)函數(shù)(16)是要將應(yīng)急系統(tǒng)的總成本最小化,包括建設(shè)成本、設(shè)備購(gòu)置成本和儲(chǔ)存成本;約束條件(17)是應(yīng)急服務(wù)的可靠性約束,其中是滿足可靠性約束(2)的參數(shù)的最小值;式(18)是應(yīng)急儲(chǔ)備點(diǎn)的儲(chǔ)備能力約束;式(19)是應(yīng)急儲(chǔ)備點(diǎn)上對(duì)其單一應(yīng)急資源的數(shù)量約束;式(20)是允許建設(shè)的應(yīng)急儲(chǔ)備點(diǎn)的最大量約束;式(21)~(22)是變量的取值約束。

    3 基于矩陣實(shí)數(shù)編碼的遺傳算法

    應(yīng)急資源布局優(yōu)化問(wèn)題早已經(jīng)被證明是一個(gè)典型的NP-Hard問(wèn)題[11-12],一般的數(shù)學(xué)方法很難在有效時(shí)間內(nèi)求解上述問(wèn)題。遺傳算法(Genetic Algorithm,GA)是一種現(xiàn)代啟發(fā)式算法,算法計(jì)算過(guò)程簡(jiǎn)單且求解效率高,對(duì)搜索空間有廣泛的適用性[13],因此本文選擇基于GA設(shè)計(jì)優(yōu)化模型的求解方法。

    GA中解的編碼形式對(duì)算法的效率和性能有本質(zhì)影響。矩陣形式編碼的GA對(duì)于求解高維、多峰和非線性的優(yōu)化問(wèn)題具有突出的優(yōu)勢(shì),尤其是針對(duì)本文所研究問(wèn)題,使用實(shí)數(shù)矩陣編碼形式,可有效避免編碼過(guò)長(zhǎng)或編碼、解碼復(fù)雜等現(xiàn)象,同時(shí)矩陣染色體可以直接應(yīng)用遺傳操作進(jìn)行尋優(yōu),無(wú)須將問(wèn)題分解成2層優(yōu)化問(wèn)題進(jìn)行求解,從而提高遺傳算法的速度和尋優(yōu)精度[14-16]。因此本文基于矩陣實(shí)數(shù)編碼的遺傳算法(Matrix Real-coded Genetic Algorithm,MRGA)設(shè)計(jì)模型的求解算法并進(jìn)行計(jì)算分析。

    在使用 MRGA進(jìn)行模型求解時(shí),首先需計(jì)算求解滿足式(2)約束的最小設(shè)備數(shù)量,由于此前已經(jīng)證明函數(shù) f ()為單調(diào)遞增函數(shù),因此可選擇二分法進(jìn)行求取。

    3.1 求解的二分法

    步驟 1:令n1=0,分別計(jì)算 f ′(n1)和 f ′(n2),轉(zhuǎn)步驟2;

    否 則 , 若 f′(n1) ? f ′( n ) > 0 , 則 令 n1=n; 若f′(n1) ? f ′( n) < 0 ,則令n2=n;

    步驟4:重復(fù)步驟2和3,直至算法終止。

    3.2 MRGA算法

    1) 編碼和初始化

    矩陣編碼方法是將矩陣整體作為遺傳子代個(gè)體,不需將矩陣展開(kāi)成一串元素,能確保子代個(gè)體基因的完整性[13],以矩陣作為個(gè)體,即矩陣染色體。每個(gè)個(gè)體矩陣中的元素即為矩陣染色體的基因。種群由多個(gè)個(gè)體矩陣組成。

    用Y1×J表示變量yj的矩陣,XM×J表示 xm的矩

    j陣。用s表示初始種群個(gè)體的個(gè)數(shù),隨機(jī)生成s個(gè)滿足解條件的矩陣Y和X,由于模型有2個(gè)矩陣變作為初始群體的矩陣染色體,將生成的s個(gè)矩陣染色體作為初始種群。矩陣為(1+M)×J階矩陣,遺傳子代第k代種群Pk個(gè)體的數(shù)目為 s,令 Pk={R1, R2,…,Rs}表示該種群,其中:

    XM×J的行表示設(shè)備資源m在各備選點(diǎn)的數(shù)量,列表示在備選點(diǎn)j處的各種設(shè)備的數(shù)量。Rl表示k種群中第l個(gè)個(gè)體,即矩陣染色體,為矩陣染色體的基因元素。

    在進(jìn)行選擇、交叉、變異等遺傳運(yùn)算時(shí)對(duì)合并矩陣R操作,在計(jì)算適應(yīng)度時(shí)將合并矩陣分解。

    2) 適應(yīng)度計(jì)算

    定義適應(yīng)度函數(shù)如下:

    計(jì)算時(shí)將種群中的矩陣染色體分解為矩陣X和Y代入目標(biāo)函數(shù),求解總成本C,MAX表示一個(gè)足夠大的正數(shù),C越小,適應(yīng)度值Fit越大,表明該個(gè)體的適應(yīng)性越強(qiáng)。

    3) 選擇復(fù)制

    對(duì)群體進(jìn)行輪盤賭選擇操作,且保留最優(yōu)個(gè)體。具體而言,尋找種群中最大適應(yīng)度對(duì)應(yīng)的矩陣染色體,并用之取代適應(yīng)度最小的矩陣染色體和當(dāng)前種群第一個(gè)矩陣染色體。在實(shí)際操作運(yùn)算中,交叉、變異等操作可能會(huì)改變所有子代矩陣染色體導(dǎo)致適應(yīng)度退化,為了防止此類現(xiàn)象的發(fā)生,確保最優(yōu)個(gè)體的遺傳,本算法的交叉、變異運(yùn)算不對(duì)種群第一個(gè)矩陣染色體實(shí)施,即采用最優(yōu)保存策略[13]。

    4) 交叉

    對(duì)群體除了第一個(gè)矩陣染色體以外的其他染色體兩兩交叉其對(duì)應(yīng)的列基因元素,交叉的概率設(shè)為pc,交叉的列是隨機(jī)選取的[14]。

    5) 變異

    采用在個(gè)體的行向量上進(jìn)行多窗口變異算子操作[15],使群體除第一個(gè)矩陣染色體外的其他矩陣染色體的行基因元素,按一定變異規(guī)則選取其他值,變異窗口寬度為w,隨機(jī)確定變異操作的行向量和窗口變異的起始位,變異的概率設(shè)為pm,每次變異中進(jìn)行2種不同窗口寬度(分別為1和2)的窗口變異操作,具體方法可參見(jiàn)文獻(xiàn)[15]。不同之處在于第一行向量變異為0-1變異,其余行向量為實(shí)整數(shù)變異。

    4 算例分析

    4.1 基礎(chǔ)數(shù)據(jù)準(zhǔn)備

    待選的應(yīng)急儲(chǔ)備點(diǎn)和應(yīng)急點(diǎn)分別為5和15個(gè),考慮2種應(yīng)急設(shè)備,即M={1,2};其應(yīng)急設(shè)備儲(chǔ)備能力占用率分別為 μ1=0.8,μ2=1.4;設(shè)備繁忙概率δ1=0.4,δ2=0.3;允許建設(shè)的最多儲(chǔ)備點(diǎn)數(shù)W=4,應(yīng)急服務(wù)的可靠性約束統(tǒng)一設(shè)置為 α=0.8,銀行利率水平ε=0.35%,規(guī)劃期為H=5年。備選的應(yīng)急儲(chǔ)備點(diǎn)信息、設(shè)備的單位購(gòu)置成本和儲(chǔ)存成本分別如下表1和2所示,所有應(yīng)急點(diǎn)的相關(guān)信息如表3所示。

    MRGA算法相關(guān)參數(shù)設(shè)置如下:種群個(gè)數(shù)為20,交叉概率pc=0.8,變異概率pm=0.1,種群進(jìn)化代 500。本文所提出的所有算法,均是在 Intel G3250,4G Memory的個(gè)人電腦上,使用MATLAB R2013a編程實(shí)現(xiàn)。

    表1 應(yīng)急儲(chǔ)備點(diǎn)基本信息Table 1 Information of emergency reserve points

    表2 應(yīng)急設(shè)備資源信息Table 2 Information of emergency resource

    4.2 計(jì)算分析

    首先分析可靠性約束對(duì)設(shè)備數(shù)量的影響。利用提出的二分法,很容易的可計(jì)算得到當(dāng)滿足可靠性約束α=0.8時(shí),覆蓋應(yīng)急點(diǎn)i的所有應(yīng)急儲(chǔ)備點(diǎn)上的設(shè)備m∈M的最小數(shù)量,具體如下表4所示。

    從表4中可發(fā)現(xiàn),由于對(duì)應(yīng)急點(diǎn)服務(wù)的可靠性約束,導(dǎo)致其周邊應(yīng)急儲(chǔ)備點(diǎn)為該應(yīng)急點(diǎn)準(zhǔn)備的服務(wù)設(shè)備數(shù)量,要明顯大于其所需的實(shí)際數(shù)量。

    表3 應(yīng)急點(diǎn)相關(guān)信息Table 3 Information of emergency demand points

    考慮到啟發(fā)式算法在求解時(shí),所得最優(yōu)解的質(zhì)量可能會(huì)受到初始解的影響,因此在算法的實(shí)現(xiàn)中,分別計(jì)算并記錄了 10組在不同初始解條件下所求取的最優(yōu)解,并以此對(duì)前文所提出的優(yōu)化模型和求解算法進(jìn)行驗(yàn)證。表5給出了不同初始解下,MRGA算法所求取的最優(yōu)解情況。

    表4 最小設(shè)備數(shù)量Table 4 Minimum numbers of machine resource

    表4 最小設(shè)備數(shù)量Table 4 Minimum numbers of machine resource

    應(yīng)急點(diǎn)i 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 N m=2 10 5 8 5 7 12 19 0 13 10 4 13 12 7 13 m i m=1 10 14 12 14 10 19 17 12 16 14 8 19 0 6 2

    圖1 N 1隨可靠性水平α變化情況iFig. 1 Minimum numberN 1 vs reliability αi

    圖2 N 2 隨可靠性水平α變化情況iFig. 2 Minimum numberN 2 vs reliability αi

    從表5可看到,在初始解不同的前提下,MRGA算法 10次計(jì)算所求得的最優(yōu)解是很接近的,最大值(213 270)與最小值(212 080)的相對(duì)偏差僅為0.56%,且相對(duì)與初始解,使用算法求解的最優(yōu)解目標(biāo)值均有12.74%~21.98%的節(jié)省,每次求解運(yùn)算時(shí)間絕大部分都在2 s以內(nèi),這也證明了本文提出的模型和MRGA算法的有效性。

    圖3描繪了當(dāng)給定的最大覆蓋距離Di發(fā)生變化時(shí),總費(fèi)用 C的變化情況(假設(shè)是所有距離都按同一幅度進(jìn)行同時(shí)變化)。從圖中可看到,隨著 Di的不斷增大,C隨之呈現(xiàn)持續(xù)降低的趨勢(shì)。但是當(dāng)Di過(guò)小時(shí),會(huì)導(dǎo)致問(wèn)題無(wú)最優(yōu)解(圖中虛線部分)。這主要是由于隨著Di的增加,同一應(yīng)急點(diǎn)可被更多的應(yīng)急儲(chǔ)備點(diǎn)覆蓋,因此各儲(chǔ)備點(diǎn)上可儲(chǔ)存相對(duì)更少的設(shè)備來(lái)保證達(dá)到所需的服務(wù)水平,從而有效降低設(shè)備購(gòu)置成本和庫(kù)存成本。而無(wú)可行解的情況,則主要是因?yàn)榇藭r(shí)給定的Di過(guò)小,導(dǎo)致某些應(yīng)急點(diǎn)始終不能被任何儲(chǔ)備點(diǎn)所覆蓋,從而使原問(wèn)題無(wú)可行解。

    表5 不同初始解對(duì)應(yīng)的最優(yōu)解情況Table 5 Initial solutions vs the optimal solutions

    圖3 系統(tǒng)總費(fèi)用C隨最大覆蓋距離Di變化情況Fig. 3 Total cost C vs the maximum distance Di

    圖4 總費(fèi)用隨可靠性水平α變化情況Fig. 4 Total cost C vs reliability α

    圖4描繪了C隨可靠性水平α變化的情形。隨著α的增加,C也隨之增大,且增長(zhǎng)速度不斷加快。這是由于要維持較高的服務(wù)可靠性,就必須儲(chǔ)存更多的應(yīng)急設(shè)備以備不時(shí)之需,從而導(dǎo)致購(gòu)置成本和庫(kù)存成本會(huì)顯著上升,尤其是當(dāng)α達(dá)到0.7以后,隨著α的進(jìn)一步增大,C的增長(zhǎng)速度明顯加快。

    5 結(jié)論

    1) 在不確定環(huán)境下,將應(yīng)急需求點(diǎn)對(duì)應(yīng)急資源服務(wù)的可靠性約束,轉(zhuǎn)化為應(yīng)急資源儲(chǔ)備點(diǎn)可提供救援服務(wù)應(yīng)急資源數(shù)量的約束,在此基礎(chǔ)上以建立了交通網(wǎng)絡(luò)應(yīng)急資源布局優(yōu)化模型,并結(jié)合二分法和矩陣編碼遺傳算法設(shè)計(jì)求解方法。算例分析表明,所提出的優(yōu)化方法計(jì)算快捷,所得結(jié)果合理。

    2) 可靠性水平和服務(wù)覆蓋距離都會(huì)對(duì)總費(fèi)用產(chǎn)生影響,且服務(wù)覆蓋距離的變化與總費(fèi)用的變化呈現(xiàn)出反比關(guān)系,可靠性水平與總費(fèi)用的變化呈現(xiàn)出正比關(guān)系,即系統(tǒng)總費(fèi)用會(huì)隨著覆蓋距離的增大而減小,隨著可靠性水平的增大而增大。

    3) 決策者若需要控制系統(tǒng)總成本,可以考慮設(shè)法增大應(yīng)急點(diǎn)的有效覆蓋距離,在降低成本的同時(shí)不會(huì)降低系統(tǒng)的服務(wù)水平,或設(shè)置合理恰當(dāng)?shù)目煽啃运揭蟆?/p>

    [1] Masood A B, Mortagy K, Alsayed A C. A multi-objective model for locating fire stations[J]. European Journal of Operational Research, 1998, 110(2): 243-260.

    [2] 付德強(qiáng), 張偉, 吳先鋒, 等. 基于NSGA-Ⅱ的災(zāi)害隨機(jī)環(huán)境下應(yīng)急儲(chǔ)備庫(kù)多目標(biāo)選址決策模型研究[J]. 軟科學(xué), 2014, 23(4): 64-69.FU Deqiang, ZHANG Wei, WU Xianfeng, et al. The decision model of the multi-objective supplies reserve storage location problem under the environment of random risks based NSGA-Ⅱ[J]. Soft Science, 2014,23(4): 64-69.

    [3] 丁雪楓, 尤建新, 王洪豐. 突發(fā)事件應(yīng)急設(shè)施選址問(wèn)題的模型及優(yōu)化算法[J]. 同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版),2012, 40(9): 1428-1433.DING Xuefeng, YOU Jianxin, WANG Hongfeng. A model of emergency rescue facility location problem for emergency incidents and its optimal solution algorithm[J].Journal of Tongji University (Natural Science), 2012,40(9): 1428-1433.

    [4] 楊金順, 孫洪運(yùn), 李林波, 等. 公路網(wǎng)災(zāi)害應(yīng)急救援點(diǎn)多目標(biāo)選址模型及算法[J]. 同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 41(12): 1843-1848.YANG Jinshun, SUN Hongyun, LI Linbo, et al.Multi-objectives location planning model and solution algorithm for disasters rescue site of highway network [J].Journal of Tongji University (Natural Science), 2013,41(12): 1843-1848.

    [5] 王旭坪, 董莉, 陳明天. 考慮感知滿意度的多受災(zāi)點(diǎn)應(yīng)急資源分配模型[J]. 系統(tǒng)管理學(xué)報(bào), 2013, 22(2): 251-256.WANG Xuping, DONG Li, CHEN Mingtian.Multiple-area post-disaster resource distribution model considering perception satisfaction[J]. Journal of Systems& Management, 2013, 22(2): 251-256.

    [6] 張玲, 董銀紅, 張敏. 基于情景分析的應(yīng)急資源布局決策[J]. 系統(tǒng)工程, 2014(3): 137-142.ZHANG Ling, DONG Yinhong, ZHANG Min.Decision-making of emergency resource distribution based on scenario analysis[J]. Systems Engineering,2014(3): 137-142.

    [7] Ceyhun A, Hasan S, Irem O. A fuzzymulti-objective covering-based vehicle location model for emergency services[J]. Services Original Research Article Computers & Operations Research, 2007, 34(3):705-726.

    [8] 艾云飛, 呂靖, 張麗麗. 水上應(yīng)急物資儲(chǔ)備庫(kù)選址-分配優(yōu)化模型[J]. 大連海事大學(xué)學(xué)報(bào), 2015(2): 62-66.AI Yunfei, Lü Jing, ZHANG Lili. The optimization location-allocation model of water emergency supplies repertories[J]. Journal of Dalian Maritime University,2015(2): 62-66.

    [9] 俞武揚(yáng). 服務(wù)能力受損情景下的應(yīng)急設(shè)施選址模型[J].控制與決策, 2016, 31(11): 1979-1984.YU Wuyang. Emergency facility location model under service ability damage scenarios[J]. Control and Decision,2016, 31(11): 1979-1984.

    [10] Revelle C, Hogan K. A reliability constrained siting model with localestimates of busy fractions[J].Environment and Planning B: Planning and Design, 2015,15(2): 143-152.

    [11] Wlodzimierz O. On the distribution approach to location problems[J]. Computers & Industrial Engineering, 1999,37(3): 595-612.

    [12] Luis G, Acosta E, Roberto D G. Dual-based heuristics for a hierarchical covering location problem[J]. Computer &Operations Research, 2003, 30(2): 165-180.

    [13] 邢文訓(xùn), 謝金星. 現(xiàn)代優(yōu)化計(jì)算方法[M]. 北京: 清華大學(xué)出版社, 2001.XING Wenxun, XIE Jinxing. Modern optimum computation method[M]. Beijing: The Tsinghua Press,2001.

    [14] 劉鯖潔, 陳桂明, 劉小方. 基于矩陣編碼的遺傳算法研究[J]. 計(jì)算機(jī)工程, 2011, 37(13): 160-162.LIU Qingjie, CHEN Guiming, LIU Xiaofang. Research on genetic algorithm based on matrix coding[J].Computer Engineering, 2011, 37(13): 160-162.

    [15] 杜永貴, 陳鑫. 矩陣編碼的遺傳算法[J]. 太原理工大學(xué)學(xué)報(bào), 2012, 43(2): 111-113, 118.DU Yonggui, CHEN Xin. Agenetic algorithm with matrix coding[J]. Journal of Taiyuan University of Technology,2012, 43(2): 111-113, 118.

    [16] 孫力勇, 張焰, 蔣傳文. 基于矩陣實(shí)數(shù)編碼遺傳算法求解大規(guī)模機(jī)組組合問(wèn)題[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2006,26(2): 82-87.SUN Liyong, ZHANG Yan, JIANG Chuanwen. A solution to the unit commitment problem based on matrix real-coded genetic algorithm[J]. Proceedings of the CSEE,2006, 26(2): 82-87.

    猜你喜歡
    儲(chǔ)備染色體可靠性
    釋放鉀肥儲(chǔ)備正當(dāng)時(shí)
    國(guó)家儲(chǔ)備林:為未來(lái)儲(chǔ)備綠色寶藏
    可靠性管理體系創(chuàng)建與實(shí)踐
    多一條X染色體,壽命會(huì)更長(zhǎng)
    為什么男性要有一條X染色體?
    電子制作(2017年2期)2017-05-17 03:55:06
    外匯儲(chǔ)備去哪兒了
    支點(diǎn)(2017年3期)2017-03-29 08:31:38
    能忍的人壽命長(zhǎng)
    再論高等植物染色體雜交
    基于可靠性跟蹤的薄弱環(huán)節(jié)辨識(shí)方法在省級(jí)電網(wǎng)可靠性改善中的應(yīng)用研究
    欧美极品一区二区三区四区| 肉色欧美久久久久久久蜜桃 | 少妇猛男粗大的猛烈进出视频 | 免费av毛片视频| 青春草亚洲视频在线观看| 大片免费播放器 马上看| 99久久精品国产国产毛片| 五月天丁香电影| 久久人人爽人人片av| 国产成人a区在线观看| 黄色日韩在线| 97人妻精品一区二区三区麻豆| 国产高清国产精品国产三级 | av专区在线播放| 亚洲成人av在线免费| 男女下面进入的视频免费午夜| 亚洲精品国产成人久久av| 老司机影院成人| av在线app专区| 伦理电影大哥的女人| 只有这里有精品99| 激情 狠狠 欧美| 熟女人妻精品中文字幕| 久久这里有精品视频免费| 午夜福利在线观看免费完整高清在| 深夜a级毛片| 亚洲av电影在线观看一区二区三区 | 内地一区二区视频在线| 国产精品一区二区在线观看99| 日本午夜av视频| 免费大片18禁| 欧美 日韩 精品 国产| 欧美高清成人免费视频www| 久久久久久九九精品二区国产| 日韩三级伦理在线观看| 婷婷色综合www| 极品少妇高潮喷水抽搐| 欧美xxⅹ黑人| 精品久久久久久久久av| 久久精品夜色国产| 免费av观看视频| 日韩人妻高清精品专区| 国产精品三级大全| 少妇 在线观看| 精品久久久久久久人妻蜜臀av| 欧美精品一区二区大全| 久久ye,这里只有精品| 日韩亚洲欧美综合| 日韩一区二区视频免费看| 国产女主播在线喷水免费视频网站| 国产精品人妻久久久影院| 高清欧美精品videossex| 国产综合懂色| 三级国产精品欧美在线观看| 亚洲内射少妇av| 亚洲精品日本国产第一区| 国产一区有黄有色的免费视频| 熟女电影av网| 亚州av有码| 日韩免费高清中文字幕av| 乱码一卡2卡4卡精品| 99视频精品全部免费 在线| 美女国产视频在线观看| 国产精品熟女久久久久浪| 最近最新中文字幕大全电影3| 你懂的网址亚洲精品在线观看| 高清日韩中文字幕在线| 国产老妇伦熟女老妇高清| 免费少妇av软件| 人妻夜夜爽99麻豆av| 久久久国产一区二区| 男人爽女人下面视频在线观看| 男人添女人高潮全过程视频| 午夜福利在线在线| 亚洲色图综合在线观看| 伊人久久国产一区二区| .国产精品久久| 久久久成人免费电影| 欧美成人一区二区免费高清观看| 新久久久久国产一级毛片| 九九在线视频观看精品| 国产黄a三级三级三级人| 久久久久精品性色| 亚洲在久久综合| av免费观看日本| av国产精品久久久久影院| 欧美成人精品欧美一级黄| 精品酒店卫生间| 视频区图区小说| 亚洲婷婷狠狠爱综合网| 久久精品熟女亚洲av麻豆精品| 麻豆久久精品国产亚洲av| 国产黄a三级三级三级人| 免费黄色在线免费观看| 日本熟妇午夜| 菩萨蛮人人尽说江南好唐韦庄| 男男h啪啪无遮挡| 97在线人人人人妻| 纵有疾风起免费观看全集完整版| 亚洲av一区综合| 国产精品久久久久久久久免| 插逼视频在线观看| 亚洲国产欧美在线一区| 国产精品99久久99久久久不卡 | 色播亚洲综合网| 久久久久久久久久久丰满| 插逼视频在线观看| 九色成人免费人妻av| 国产成人aa在线观看| 中国美白少妇内射xxxbb| 最新中文字幕久久久久| 春色校园在线视频观看| 成人亚洲精品av一区二区| 一个人看视频在线观看www免费| 黄片wwwwww| 欧美xxxx黑人xx丫x性爽| 老女人水多毛片| 两个人的视频大全免费| 校园人妻丝袜中文字幕| 久久精品国产亚洲av天美| 大陆偷拍与自拍| xxx大片免费视频| 亚洲欧美一区二区三区国产| 国产精品一区二区三区四区免费观看| av免费在线看不卡| 最近的中文字幕免费完整| 久久精品国产亚洲av涩爱| 26uuu在线亚洲综合色| 成人无遮挡网站| 国产高清三级在线| 2021天堂中文幕一二区在线观| 亚洲,一卡二卡三卡| 亚洲精品一区蜜桃| 国产精品一区二区性色av| 久久精品综合一区二区三区| 国产爽快片一区二区三区| 国产老妇女一区| 久久人人爽人人爽人人片va| 91精品国产九色| 欧美日韩综合久久久久久| 久久人人爽人人爽人人片va| 欧美 日韩 精品 国产| 最近中文字幕2019免费版| 又粗又硬又长又爽又黄的视频| 久久精品夜色国产| 欧美最新免费一区二区三区| 国产精品爽爽va在线观看网站| 极品教师在线视频| 国产精品一区二区性色av| 激情 狠狠 欧美| 在线亚洲精品国产二区图片欧美 | 国产黄色视频一区二区在线观看| 少妇高潮的动态图| 午夜激情久久久久久久| 欧美 日韩 精品 国产| 国产av码专区亚洲av| 亚洲内射少妇av| 亚洲精品中文字幕在线视频 | 涩涩av久久男人的天堂| 神马国产精品三级电影在线观看| 六月丁香七月| 亚洲国产精品999| 亚洲美女视频黄频| 51国产日韩欧美| 亚洲国产日韩一区二区| 日本猛色少妇xxxxx猛交久久| 天堂俺去俺来也www色官网| 伦理电影大哥的女人| 大码成人一级视频| 晚上一个人看的免费电影| 欧美区成人在线视频| 色哟哟·www| 婷婷色av中文字幕| 高清午夜精品一区二区三区| 一级二级三级毛片免费看| 丝瓜视频免费看黄片| 亚洲欧美一区二区三区黑人 | 久久精品综合一区二区三区| 亚洲av不卡在线观看| 精品熟女少妇av免费看| 综合色丁香网| 中文字幕久久专区| 国产一区二区三区av在线| 欧美日韩国产mv在线观看视频 | 亚洲精品久久久久久婷婷小说| 国产成人免费无遮挡视频| 亚洲在线观看片| 国产精品偷伦视频观看了| 欧美xxxx黑人xx丫x性爽| 日韩欧美一区视频在线观看 | 免费看a级黄色片| 国产免费一区二区三区四区乱码| 亚洲精品一二三| 成人国产麻豆网| 在线免费十八禁| 搞女人的毛片| 亚洲国产欧美人成| 熟妇人妻不卡中文字幕| 亚洲一区二区三区欧美精品 | 亚洲av欧美aⅴ国产| 久久久久久久精品精品| 成人欧美大片| 亚洲美女视频黄频| 国产视频内射| 三级国产精品欧美在线观看| 丝袜美腿在线中文| 日韩欧美一区视频在线观看 | 久久99蜜桃精品久久| 国语对白做爰xxxⅹ性视频网站| 国产在线男女| 国产精品国产三级国产av玫瑰| 国产中年淑女户外野战色| 人人妻人人澡人人爽人人夜夜| 大片电影免费在线观看免费| 精品99又大又爽又粗少妇毛片| 在线看a的网站| 国产精品女同一区二区软件| 亚洲在线观看片| 日本午夜av视频| 国产片特级美女逼逼视频| 日韩国内少妇激情av| 五月玫瑰六月丁香| 国产久久久一区二区三区| 欧美精品人与动牲交sv欧美| 精品久久久久久久末码| 91狼人影院| 亚洲av免费在线观看| 韩国高清视频一区二区三区| 麻豆久久精品国产亚洲av| 亚洲,一卡二卡三卡| 美女被艹到高潮喷水动态| 少妇丰满av| 日韩 亚洲 欧美在线| 五月开心婷婷网| 九色成人免费人妻av| 蜜桃久久精品国产亚洲av| 天天躁夜夜躁狠狠久久av| av国产久精品久网站免费入址| 小蜜桃在线观看免费完整版高清| 黑人高潮一二区| 一级av片app| 身体一侧抽搐| 成人无遮挡网站| 欧美人与善性xxx| 欧美精品国产亚洲| 老师上课跳d突然被开到最大视频| 国产精品福利在线免费观看| 久久久久国产精品人妻一区二区| 亚洲av日韩在线播放| 一二三四中文在线观看免费高清| 欧美人与善性xxx| 欧美精品国产亚洲| 亚洲不卡免费看| 神马国产精品三级电影在线观看| 国产精品久久久久久精品电影小说 | 久久久久久久久久人人人人人人| 91在线精品国自产拍蜜月| 国产探花在线观看一区二区| 九九久久精品国产亚洲av麻豆| 麻豆精品久久久久久蜜桃| av.在线天堂| 亚洲色图综合在线观看| 伦精品一区二区三区| 国产亚洲91精品色在线| 国产大屁股一区二区在线视频| av播播在线观看一区| av女优亚洲男人天堂| 久久久久久久久久人人人人人人| 欧美激情在线99| 人人妻人人看人人澡| 各种免费的搞黄视频| 91精品国产九色| 亚洲精品视频女| 免费观看的影片在线观看| 五月开心婷婷网| 我要看日韩黄色一级片| 欧美变态另类bdsm刘玥| 久久久国产一区二区| 久久精品国产鲁丝片午夜精品| 亚洲成人精品中文字幕电影| 麻豆成人午夜福利视频| 国产成人精品久久久久久| 少妇猛男粗大的猛烈进出视频 | 能在线免费看毛片的网站| 国产毛片在线视频| 国产亚洲91精品色在线| 国产日韩欧美在线精品| 最近最新中文字幕免费大全7| 欧美 日韩 精品 国产| 婷婷色麻豆天堂久久| 国产精品嫩草影院av在线观看| 三级经典国产精品| 亚洲国产精品成人综合色| 欧美日本视频| 狂野欧美激情性xxxx在线观看| 欧美97在线视频| 99久久精品热视频| 免费黄网站久久成人精品| 国产精品偷伦视频观看了| 乱系列少妇在线播放| 你懂的网址亚洲精品在线观看| 一二三四中文在线观看免费高清| 日韩,欧美,国产一区二区三区| 综合色av麻豆| 午夜免费观看性视频| 亚州av有码| 伊人久久国产一区二区| 日韩制服骚丝袜av| 国产在线一区二区三区精| 亚洲熟女精品中文字幕| 性色av一级| 午夜福利高清视频| 亚洲成人久久爱视频| 在线天堂最新版资源| av卡一久久| 免费看光身美女| 偷拍熟女少妇极品色| 亚洲av中文字字幕乱码综合| 国国产精品蜜臀av免费| 午夜视频国产福利| 久久精品综合一区二区三区| 80岁老熟妇乱子伦牲交| 欧美国产精品一级二级三级 | 日日啪夜夜爽| 亚洲怡红院男人天堂| 亚洲国产成人一精品久久久| 久久精品国产自在天天线| 国产精品一二三区在线看| 精品人妻偷拍中文字幕| 国产成人aa在线观看| 国产女主播在线喷水免费视频网站| 国产在视频线精品| 欧美区成人在线视频| 天天躁日日操中文字幕| 看免费成人av毛片| 国产成人免费无遮挡视频| 国产一区二区三区综合在线观看 | 欧美日韩亚洲高清精品| 3wmmmm亚洲av在线观看| 亚洲精品国产av蜜桃| 欧美日韩视频精品一区| 国产精品久久久久久久久免| 久久精品国产自在天天线| 国产一级毛片在线| 亚洲精品成人久久久久久| 韩国高清视频一区二区三区| 国产一区二区三区av在线| 建设人人有责人人尽责人人享有的 | 在线a可以看的网站| 亚洲精品国产成人久久av| 久久精品国产鲁丝片午夜精品| 精品国产一区二区三区久久久樱花 | 日本与韩国留学比较| 一级毛片电影观看| 女人被狂操c到高潮| 色网站视频免费| 欧美性感艳星| 精品一区二区三卡| 搡女人真爽免费视频火全软件| 精品人妻熟女av久视频| 亚洲精华国产精华液的使用体验| 80岁老熟妇乱子伦牲交| 99久久中文字幕三级久久日本| 成年av动漫网址| 熟女电影av网| 亚洲av一区综合| 丝袜喷水一区| 只有这里有精品99| 亚洲av免费高清在线观看| 一级毛片aaaaaa免费看小| 亚洲激情五月婷婷啪啪| 在线观看免费高清a一片| 黄色一级大片看看| 97热精品久久久久久| 一个人看的www免费观看视频| 在线观看免费高清a一片| 久久久久久久久久久免费av| 97热精品久久久久久| 99视频精品全部免费 在线| 日韩免费高清中文字幕av| 肉色欧美久久久久久久蜜桃 | 亚洲在久久综合| 日本爱情动作片www.在线观看| 在线精品无人区一区二区三 | 久久久久久久国产电影| 国产伦精品一区二区三区视频9| 亚洲图色成人| 免费不卡的大黄色大毛片视频在线观看| 久久久a久久爽久久v久久| 精品久久久久久电影网| 又黄又爽又刺激的免费视频.| 爱豆传媒免费全集在线观看| 可以在线观看毛片的网站| 亚洲国产欧美人成| 亚洲性久久影院| 国产黄片视频在线免费观看| 波野结衣二区三区在线| 99久国产av精品国产电影| 成人毛片60女人毛片免费| 国产 精品1| 男人舔奶头视频| 麻豆乱淫一区二区| av免费观看日本| 熟妇人妻不卡中文字幕| 亚洲一级一片aⅴ在线观看| 99热这里只有是精品50| 国产精品三级大全| 人人妻人人爽人人添夜夜欢视频 | 精品久久久久久电影网| 99热这里只有是精品50| 久久久精品94久久精品| 欧美人与善性xxx| 免费观看无遮挡的男女| 一级毛片黄色毛片免费观看视频| 久久久久网色| 亚洲精品久久久久久婷婷小说| 啦啦啦啦在线视频资源| 少妇人妻 视频| 中文在线观看免费www的网站| 777米奇影视久久| 99久久精品一区二区三区| 嫩草影院精品99| 欧美精品国产亚洲| 日韩一区二区视频免费看| 舔av片在线| 日韩国内少妇激情av| av卡一久久| 亚洲国产高清在线一区二区三| 国产淫片久久久久久久久| 久久精品人妻少妇| 免费看日本二区| 亚洲美女视频黄频| 国产色婷婷99| 成人特级av手机在线观看| 欧美人与善性xxx| 国产成人免费无遮挡视频| 大陆偷拍与自拍| 国产老妇伦熟女老妇高清| 久久精品人妻少妇| 特大巨黑吊av在线直播| 成人国产av品久久久| 国产亚洲91精品色在线| kizo精华| 在线观看免费高清a一片| 日本免费在线观看一区| 亚洲精品第二区| 成人黄色视频免费在线看| 成年av动漫网址| 久久精品熟女亚洲av麻豆精品| 国产av不卡久久| 欧美老熟妇乱子伦牲交| 国产精品久久久久久av不卡| 国产免费一级a男人的天堂| 麻豆乱淫一区二区| 一区二区av电影网| 青春草视频在线免费观看| 97在线视频观看| 国产精品99久久久久久久久| 能在线免费看毛片的网站| 国产高潮美女av| 性色av一级| 亚洲美女视频黄频| 久久精品熟女亚洲av麻豆精品| 99九九线精品视频在线观看视频| 亚洲成色77777| 五月天丁香电影| 亚洲精品,欧美精品| 内地一区二区视频在线| 久久久国产一区二区| 亚洲国产精品999| 欧美成人a在线观看| 午夜福利在线观看免费完整高清在| 亚洲精品第二区| 黄片wwwwww| 久久久久久久亚洲中文字幕| 国产亚洲最大av| 黄色一级大片看看| 在线 av 中文字幕| av播播在线观看一区| 直男gayav资源| 国产av不卡久久| 特级一级黄色大片| 亚洲综合精品二区| 男女边摸边吃奶| 中文字幕制服av| 大码成人一级视频| 欧美激情国产日韩精品一区| 欧美xxxx黑人xx丫x性爽| av网站免费在线观看视频| 男人狂女人下面高潮的视频| 国产精品国产三级国产av玫瑰| 人人妻人人看人人澡| 亚洲国产精品成人综合色| 大话2 男鬼变身卡| 色视频www国产| 免费播放大片免费观看视频在线观看| 国产色婷婷99| 最近中文字幕高清免费大全6| 日本免费在线观看一区| 韩国高清视频一区二区三区| 国产精品一及| kizo精华| 丰满少妇做爰视频| 婷婷色麻豆天堂久久| 国产久久久一区二区三区| 国产黄a三级三级三级人| 免费看光身美女| 毛片女人毛片| 国产午夜精品久久久久久一区二区三区| 久久精品国产a三级三级三级| 亚洲欧洲日产国产| 大香蕉久久网| 水蜜桃什么品种好| 精品国产乱码久久久久久小说| 日韩视频在线欧美| 视频中文字幕在线观看| 制服丝袜香蕉在线| 91精品国产九色| av又黄又爽大尺度在线免费看| 尾随美女入室| 婷婷色综合www| 一区二区三区精品91| 中文字幕人妻熟人妻熟丝袜美| 精品一区在线观看国产| 国产精品一区二区三区四区免费观看| 国产精品久久久久久精品电影小说 | 菩萨蛮人人尽说江南好唐韦庄| 80岁老熟妇乱子伦牲交| 大片免费播放器 马上看| 久久精品人妻少妇| 直男gayav资源| 国产精品爽爽va在线观看网站| 黄色日韩在线| 日韩一区二区三区影片| 色吧在线观看| 中文资源天堂在线| 婷婷色综合www| 亚洲自拍偷在线| 国产日韩欧美在线精品| av国产精品久久久久影院| 两个人的视频大全免费| 午夜精品国产一区二区电影 | 国产色爽女视频免费观看| 久久人人爽av亚洲精品天堂 | 日本猛色少妇xxxxx猛交久久| 精品酒店卫生间| 国产成人aa在线观看| 国产精品秋霞免费鲁丝片| 久久久欧美国产精品| 丰满少妇做爰视频| av播播在线观看一区| 国产日韩欧美在线精品| 国产av国产精品国产| 边亲边吃奶的免费视频| 九九爱精品视频在线观看| 精品一区在线观看国产| 九九在线视频观看精品| 成人二区视频| 国产成人福利小说| 黄色日韩在线| 国产亚洲一区二区精品| xxx大片免费视频| 久久久久国产网址| 黄色视频在线播放观看不卡| 精品久久久久久电影网| 国产毛片在线视频| 久久久久网色| 日韩不卡一区二区三区视频在线| 18+在线观看网站| 婷婷色综合www| 五月伊人婷婷丁香| 午夜视频国产福利| 久久久久九九精品影院| 国产精品福利在线免费观看| 夜夜爽夜夜爽视频| 免费观看无遮挡的男女| 亚洲av中文av极速乱| 最近最新中文字幕大全电影3| 久久久久久久精品精品| 22中文网久久字幕| 精品人妻一区二区三区麻豆| 超碰av人人做人人爽久久| 大话2 男鬼变身卡| 亚洲美女搞黄在线观看| 久久影院123| 欧美xxⅹ黑人| 熟妇人妻不卡中文字幕| 欧美性感艳星| 国产高潮美女av| 伊人久久国产一区二区| 天美传媒精品一区二区| 午夜福利网站1000一区二区三区| 久久精品国产亚洲av天美| 欧美一区二区亚洲| 欧美精品一区二区大全| 精品人妻一区二区三区麻豆| 嫩草影院新地址| 男女边摸边吃奶| 国产精品麻豆人妻色哟哟久久| 69av精品久久久久久| 蜜臀久久99精品久久宅男| 午夜精品一区二区三区免费看| 一二三四中文在线观看免费高清| 日产精品乱码卡一卡2卡三| 国产精品麻豆人妻色哟哟久久| 久久6这里有精品| 国产精品国产三级专区第一集| 日本一本二区三区精品| 卡戴珊不雅视频在线播放| 亚洲真实伦在线观看| 联通29元200g的流量卡| 丰满少妇做爰视频| 3wmmmm亚洲av在线观看| 欧美日本视频| 亚洲欧美中文字幕日韩二区|