• <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)用研究
    99久国产av精品国产电影| 亚洲av男天堂| 亚洲av福利一区| 黄色毛片三级朝国网站| 狠狠精品人妻久久久久久综合| 欧美变态另类bdsm刘玥| 国产免费一级a男人的天堂| 亚洲精品视频女| 久久精品国产亚洲av涩爱| 亚洲美女视频黄频| 九九在线视频观看精品| 成人午夜精彩视频在线观看| 伊人久久国产一区二区| 国产亚洲一区二区精品| 99热这里只有精品一区| 内地一区二区视频在线| 国产一区有黄有色的免费视频| 国产淫语在线视频| 肉色欧美久久久久久久蜜桃| 大片电影免费在线观看免费| 亚洲欧美一区二区三区国产| 少妇被粗大的猛进出69影院 | 国产日韩欧美在线精品| 男女免费视频国产| 国产精品一国产av| 一级,二级,三级黄色视频| 亚洲成人av在线免费| 久久av网站| 欧美成人午夜免费资源| 欧美激情 高清一区二区三区| 久久女婷五月综合色啪小说| 菩萨蛮人人尽说江南好唐韦庄| 国精品久久久久久国模美| 啦啦啦啦在线视频资源| 国产毛片在线视频| 日产精品乱码卡一卡2卡三| √禁漫天堂资源中文www| 老司机影院成人| 亚洲国产毛片av蜜桃av| 国产成人一区二区在线| 国产色婷婷99| 女的被弄到高潮叫床怎么办| 老司机影院成人| 欧美日韩国产mv在线观看视频| 特大巨黑吊av在线直播| 十八禁网站网址无遮挡| 国产精品秋霞免费鲁丝片| 中文字幕av电影在线播放| 婷婷成人精品国产| videosex国产| 国产精品国产三级国产av玫瑰| videossex国产| 精品99又大又爽又粗少妇毛片| 免费看光身美女| 免费人妻精品一区二区三区视频| 免费观看性生交大片5| 在线观看www视频免费| 久久精品国产亚洲av天美| 国产精品久久久久久精品古装| 男男h啪啪无遮挡| 久久午夜综合久久蜜桃| 亚洲欧洲精品一区二区精品久久久 | 久久久国产精品麻豆| 五月伊人婷婷丁香| videossex国产| 中文字幕人妻熟人妻熟丝袜美| 亚洲人成77777在线视频| 亚洲美女视频黄频| 国产在视频线精品| 在现免费观看毛片| 欧美日韩综合久久久久久| 一个人免费看片子| 日韩视频在线欧美| www.av在线官网国产| 伊人亚洲综合成人网| 草草在线视频免费看| av又黄又爽大尺度在线免费看| 精品久久久久久电影网| 日韩欧美精品免费久久| av在线播放精品| 一边亲一边摸免费视频| 久久国产亚洲av麻豆专区| 美女视频免费永久观看网站| 美女脱内裤让男人舔精品视频| 国产成人免费观看mmmm| 欧美xxxx性猛交bbbb| 久久久国产欧美日韩av| 国产精品久久久久久久电影| 97在线视频观看| 少妇人妻精品综合一区二区| 国产在线视频一区二区| www.色视频.com| 日韩免费高清中文字幕av| 国产色爽女视频免费观看| 天堂中文最新版在线下载| 777米奇影视久久| 精品一区二区三区视频在线| 蜜桃在线观看..| 街头女战士在线观看网站| 黄片播放在线免费| 国产亚洲精品第一综合不卡 | 亚洲精品久久久久久婷婷小说| 在线观看国产h片| 日本黄色片子视频| 老司机影院毛片| 一级毛片 在线播放| 亚洲欧美清纯卡通| 久久毛片免费看一区二区三区| 最新的欧美精品一区二区| 成人影院久久| 国产深夜福利视频在线观看| 人成视频在线观看免费观看| 黄色怎么调成土黄色| 中文字幕免费在线视频6| 国产成人精品在线电影| 午夜久久久在线观看| 一边摸一边做爽爽视频免费| 波野结衣二区三区在线| 国产成人精品婷婷| 蜜桃国产av成人99| 肉色欧美久久久久久久蜜桃| 青青草视频在线视频观看| 特大巨黑吊av在线直播| 交换朋友夫妻互换小说| 日日摸夜夜添夜夜爱| 精品99又大又爽又粗少妇毛片| 国产有黄有色有爽视频| 丰满迷人的少妇在线观看| 欧美成人精品欧美一级黄| 久久女婷五月综合色啪小说| 久久久精品94久久精品| 亚洲国产欧美日韩在线播放| 国产成人精品久久久久久| 国产亚洲欧美精品永久| 另类精品久久| 国产精品久久久久久久电影| 国产精品女同一区二区软件| 日日啪夜夜爽| 22中文网久久字幕| 亚洲精品av麻豆狂野| 高清在线视频一区二区三区| 日韩av不卡免费在线播放| 9色porny在线观看| 乱码一卡2卡4卡精品| 成人毛片60女人毛片免费| 日韩强制内射视频| 乱人伦中国视频| 制服丝袜香蕉在线| 免费大片18禁| 在线观看一区二区三区激情| 亚洲精品,欧美精品| 中文字幕av电影在线播放| 免费黄频网站在线观看国产| 久久99精品国语久久久| av卡一久久| 日本欧美视频一区| 五月伊人婷婷丁香| 午夜激情av网站| 人妻制服诱惑在线中文字幕| 七月丁香在线播放| 九色亚洲精品在线播放| 老女人水多毛片| 91精品三级在线观看| 国产高清不卡午夜福利| 亚洲精品久久久久久婷婷小说| 日本av手机在线免费观看| 乱人伦中国视频| 天天操日日干夜夜撸| 中国国产av一级| 久久久久久久国产电影| 久久久久久久久久久久大奶| av卡一久久| 中文字幕久久专区| 欧美人与性动交α欧美精品济南到 | 国产成人精品在线电影| 国产视频首页在线观看| 自线自在国产av| 制服丝袜香蕉在线| 久久97久久精品| 十八禁高潮呻吟视频| 制服丝袜香蕉在线| 中文字幕av电影在线播放| 欧美精品高潮呻吟av久久| 亚洲激情五月婷婷啪啪| 能在线免费看毛片的网站| 成年女人在线观看亚洲视频| 26uuu在线亚洲综合色| 婷婷色综合www| 免费看不卡的av| 日韩中字成人| 大香蕉97超碰在线| 在线播放无遮挡| 美女主播在线视频| av女优亚洲男人天堂| 尾随美女入室| 亚洲综合色网址| 三级国产精品片| 午夜免费观看性视频| 美女脱内裤让男人舔精品视频| 亚洲国产精品专区欧美| 欧美日韩视频高清一区二区三区二| 一本大道久久a久久精品| 在线亚洲精品国产二区图片欧美 | 中国国产av一级| 三上悠亚av全集在线观看| 亚洲av日韩在线播放| 视频中文字幕在线观看| 免费少妇av软件| 国产一区二区三区av在线| 午夜91福利影院| 久久婷婷青草| 亚洲欧洲国产日韩| 国产日韩欧美视频二区| 免费观看在线日韩| 天美传媒精品一区二区| 欧美97在线视频| 黄片无遮挡物在线观看| av女优亚洲男人天堂| 日韩,欧美,国产一区二区三区| 亚洲人成77777在线视频| 97在线视频观看| 乱人伦中国视频| 汤姆久久久久久久影院中文字幕| 日本欧美国产在线视频| 久久久久久久久久久免费av| 亚洲情色 制服丝袜| 午夜影院在线不卡| 草草在线视频免费看| 亚洲欧美日韩另类电影网站| 国产精品一区二区在线观看99| 久久ye,这里只有精品| 亚洲精品一区蜜桃| 日本色播在线视频| 秋霞在线观看毛片| 伊人亚洲综合成人网| 99久久综合免费| 国产又色又爽无遮挡免| 大码成人一级视频| 国产一级毛片在线| 精品人妻偷拍中文字幕| 我要看黄色一级片免费的| 人人妻人人爽人人添夜夜欢视频| av一本久久久久| 亚洲少妇的诱惑av| 天天躁夜夜躁狠狠久久av| 王馨瑶露胸无遮挡在线观看| 国产精品人妻久久久久久| 中国三级夫妇交换| 免费黄网站久久成人精品| 国产免费一级a男人的天堂| 青春草亚洲视频在线观看| 国产乱来视频区| 欧美成人午夜免费资源| 九草在线视频观看| 国产午夜精品久久久久久一区二区三区| 日本猛色少妇xxxxx猛交久久| 91精品伊人久久大香线蕉| 亚洲精品乱码久久久v下载方式| 国产精品一区二区三区四区免费观看| 国产免费一级a男人的天堂| 日韩av免费高清视频| 久热久热在线精品观看| 国产精品国产三级专区第一集| 一本久久精品| 黑丝袜美女国产一区| 观看美女的网站| 亚洲精品久久成人aⅴ小说 | 午夜激情福利司机影院| 亚洲av日韩在线播放| 亚洲高清免费不卡视频| 国产欧美亚洲国产| 久久午夜综合久久蜜桃| 免费黄网站久久成人精品| 狠狠婷婷综合久久久久久88av| 国产熟女午夜一区二区三区 | 永久免费av网站大全| 激情五月婷婷亚洲| 如何舔出高潮| 亚洲av二区三区四区| 日韩中文字幕视频在线看片| 一边摸一边做爽爽视频免费| 久久鲁丝午夜福利片| 少妇人妻久久综合中文| 99久久精品国产国产毛片| 中国美白少妇内射xxxbb| 久久精品久久久久久噜噜老黄| 亚洲精品日本国产第一区| av在线老鸭窝| 18在线观看网站| 久久国产亚洲av麻豆专区| 日韩电影二区| 大香蕉久久成人网| 啦啦啦视频在线资源免费观看| 成年av动漫网址| 国产探花极品一区二区| 精品久久久久久久久亚洲| 国产成人免费无遮挡视频| 成年av动漫网址| 精品熟女少妇av免费看| 免费av中文字幕在线| 亚洲成人一二三区av| 亚洲精品亚洲一区二区| 国产精品99久久99久久久不卡 | 三级国产精品片| 亚洲不卡免费看| 国产成人freesex在线| 一区在线观看完整版| 国产成人a∨麻豆精品| 女的被弄到高潮叫床怎么办| 如日韩欧美国产精品一区二区三区 | 纯流量卡能插随身wifi吗| 国产精品偷伦视频观看了| 丝袜脚勾引网站| 熟女电影av网| 我要看黄色一级片免费的| 波野结衣二区三区在线| videossex国产| 五月开心婷婷网| 天堂8中文在线网| 国产不卡av网站在线观看| 国产亚洲欧美精品永久| 久久婷婷青草| 欧美国产精品一级二级三级| 一区二区日韩欧美中文字幕 | 天美传媒精品一区二区| 熟女av电影| 亚洲欧洲日产国产| 日韩欧美一区视频在线观看| xxx大片免费视频| 熟妇人妻不卡中文字幕| 美女主播在线视频| 成年av动漫网址| 午夜av观看不卡| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品一区二区三区在线| 一本一本综合久久| 91久久精品电影网| 曰老女人黄片| 欧美精品人与动牲交sv欧美| 人妻少妇偷人精品九色| 大话2 男鬼变身卡| av女优亚洲男人天堂| 满18在线观看网站| 中文乱码字字幕精品一区二区三区| 国产有黄有色有爽视频| 亚洲四区av| 国产av码专区亚洲av| 欧美日韩国产mv在线观看视频| 人人妻人人澡人人爽人人夜夜| 国产av国产精品国产| 伦精品一区二区三区| 国产高清国产精品国产三级| 少妇被粗大的猛进出69影院 | 久久97久久精品| 久久毛片免费看一区二区三区| 欧美国产精品一级二级三级| 亚洲人成网站在线观看播放| √禁漫天堂资源中文www| av专区在线播放| 一级a做视频免费观看| 一级毛片aaaaaa免费看小| 精品午夜福利在线看| 少妇精品久久久久久久| 久久人人爽人人爽人人片va| 80岁老熟妇乱子伦牲交| 日韩强制内射视频| 亚洲久久久国产精品| 日本av手机在线免费观看| 91精品三级在线观看| 精品少妇久久久久久888优播| 国产探花极品一区二区| 在线观看免费高清a一片| 免费大片黄手机在线观看| 国精品久久久久久国模美| 不卡视频在线观看欧美| 免费人妻精品一区二区三区视频| 99热这里只有是精品在线观看| 国产黄频视频在线观看| 日本黄色片子视频| 精品人妻偷拍中文字幕| 久久久久国产精品人妻一区二区| 汤姆久久久久久久影院中文字幕| 亚洲精品456在线播放app| 观看av在线不卡| 亚洲成人一二三区av| 国产 一区精品| 晚上一个人看的免费电影| 欧美日韩视频精品一区| 晚上一个人看的免费电影| 观看av在线不卡| 免费看av在线观看网站| 又大又黄又爽视频免费| 国产亚洲av片在线观看秒播厂| av免费在线看不卡| 免费少妇av软件| 亚洲内射少妇av| 赤兔流量卡办理| 成人漫画全彩无遮挡| 欧美+日韩+精品| 校园人妻丝袜中文字幕| 久久精品人人爽人人爽视色| 成人无遮挡网站| 欧美最新免费一区二区三区| 最近手机中文字幕大全| 人妻夜夜爽99麻豆av| 国产精品嫩草影院av在线观看| 欧美精品人与动牲交sv欧美| 久久久精品区二区三区| 3wmmmm亚洲av在线观看| 18禁观看日本| 国产精品秋霞免费鲁丝片| 精品亚洲乱码少妇综合久久| 爱豆传媒免费全集在线观看| 午夜免费观看性视频| 国产精品一区二区在线不卡| 精品人妻熟女毛片av久久网站| 久久青草综合色| 精品酒店卫生间| 嫩草影院入口| 高清黄色对白视频在线免费看| 九色亚洲精品在线播放| 欧美日韩视频精品一区| 亚洲美女搞黄在线观看| 色吧在线观看| 亚洲丝袜综合中文字幕| 久久ye,这里只有精品| 久久亚洲国产成人精品v| 久久青草综合色| 亚洲人成77777在线视频| 高清欧美精品videossex| 三级国产精品欧美在线观看| 99热全是精品| 69精品国产乱码久久久| 丝袜在线中文字幕| 亚洲欧洲国产日韩| 久久久国产欧美日韩av| 丰满乱子伦码专区| 91精品三级在线观看| av卡一久久| 肉色欧美久久久久久久蜜桃| 少妇 在线观看| 中文欧美无线码| 成年av动漫网址| 亚洲精品乱久久久久久| 国产成人午夜福利电影在线观看| 亚洲欧美日韩卡通动漫| 久久人妻熟女aⅴ| 狂野欧美白嫩少妇大欣赏| 日本av免费视频播放| 老熟女久久久| 国产精品免费大片| 十八禁网站网址无遮挡| 国产在视频线精品| 亚洲精品视频女| 国产精品久久久久久久久免| 亚洲综合色惰| 考比视频在线观看| 国产黄片视频在线免费观看| 全区人妻精品视频| 久久人人爽人人片av| 一区在线观看完整版| 97在线人人人人妻| 亚洲内射少妇av| 高清视频免费观看一区二区| 亚洲国产最新在线播放| 男女啪啪激烈高潮av片| 纯流量卡能插随身wifi吗| 成人国语在线视频| 人妻 亚洲 视频| 精品亚洲乱码少妇综合久久| 中文欧美无线码| 亚洲高清免费不卡视频| 黄色毛片三级朝国网站| 80岁老熟妇乱子伦牲交| 国国产精品蜜臀av免费| 久久久a久久爽久久v久久| 婷婷色综合大香蕉| 蜜桃国产av成人99| 国产熟女午夜一区二区三区 | av有码第一页| 在线观看美女被高潮喷水网站| 久久久精品免费免费高清| 国产伦精品一区二区三区视频9| 亚洲成人一二三区av| 91午夜精品亚洲一区二区三区| 插阴视频在线观看视频| 人妻一区二区av| 久久精品国产自在天天线| 国产精品蜜桃在线观看| 久久人妻熟女aⅴ| 美女cb高潮喷水在线观看| 看十八女毛片水多多多| 国产成人av激情在线播放 | 22中文网久久字幕| 女性被躁到高潮视频| 午夜福利在线观看免费完整高清在| 丝袜美足系列| videosex国产| 99九九在线精品视频| 夫妻性生交免费视频一级片| 久久 成人 亚洲| 久久久久久久久久人人人人人人| 午夜激情av网站| av福利片在线| 国产色婷婷99| 国产高清有码在线观看视频| 日本-黄色视频高清免费观看| 天堂俺去俺来也www色官网| 极品人妻少妇av视频| 考比视频在线观看| 精品少妇内射三级| 亚洲精品乱码久久久久久按摩| 在线亚洲精品国产二区图片欧美 | 欧美人与善性xxx| 国产在视频线精品| 久久99蜜桃精品久久| 国产亚洲av片在线观看秒播厂| 成人影院久久| 九九在线视频观看精品| 性色av一级| 成人国语在线视频| 麻豆成人av视频| 我要看黄色一级片免费的| 国产精品偷伦视频观看了| 91成人精品电影| 国产视频首页在线观看| 国产男女超爽视频在线观看| 人妻少妇偷人精品九色| 亚洲av成人精品一二三区| 人妻 亚洲 视频| 毛片一级片免费看久久久久| 亚洲精品视频女| 亚洲精品一区蜜桃| 中国国产av一级| 国产有黄有色有爽视频| 另类精品久久| 自拍欧美九色日韩亚洲蝌蚪91| 伦精品一区二区三区| 99久国产av精品国产电影| 日本色播在线视频| 大又大粗又爽又黄少妇毛片口| 亚洲,一卡二卡三卡| 中文天堂在线官网| 久久精品国产亚洲av天美| av一本久久久久| 大陆偷拍与自拍| 一级黄片播放器| 精品一区二区免费观看| a级片在线免费高清观看视频| 国模一区二区三区四区视频| freevideosex欧美| 日本91视频免费播放| 国产成人av激情在线播放 | 丰满少妇做爰视频| 汤姆久久久久久久影院中文字幕| 久久久a久久爽久久v久久| 欧美日韩视频高清一区二区三区二| 欧美xxⅹ黑人| 欧美成人精品欧美一级黄| 日本猛色少妇xxxxx猛交久久| 看十八女毛片水多多多| 99热6这里只有精品| 成人国产av品久久久| 成人漫画全彩无遮挡| 狠狠婷婷综合久久久久久88av| 欧美亚洲日本最大视频资源| 国产亚洲欧美精品永久| 女性生殖器流出的白浆| 99热全是精品| 国产一区亚洲一区在线观看| 伊人久久精品亚洲午夜| 男人爽女人下面视频在线观看| 国产欧美日韩综合在线一区二区| 99久久精品一区二区三区| 免费日韩欧美在线观看| 免费人妻精品一区二区三区视频| 午夜影院在线不卡| 91精品国产国语对白视频| 亚洲av国产av综合av卡| 精品亚洲成国产av| 搡老乐熟女国产| 全区人妻精品视频| 晚上一个人看的免费电影| 黑人猛操日本美女一级片| 久久女婷五月综合色啪小说| 天堂俺去俺来也www色官网| 中文精品一卡2卡3卡4更新| 黄色配什么色好看| 日韩 亚洲 欧美在线| 我的女老师完整版在线观看| 不卡视频在线观看欧美| 亚洲人与动物交配视频| 18禁裸乳无遮挡动漫免费视频| 精品少妇黑人巨大在线播放| 美女国产高潮福利片在线看| 亚洲国产av新网站| 日韩成人伦理影院| 午夜老司机福利剧场| 日本av免费视频播放| 日韩一区二区视频免费看| 一区二区三区免费毛片| 91aial.com中文字幕在线观看| 久久婷婷青草| av电影中文网址| 日韩中字成人| 亚洲精品一二三| 国产片特级美女逼逼视频| 亚洲国产av影院在线观看| 国产精品99久久久久久久久| 国产精品99久久99久久久不卡 | 亚洲国产精品999| 男人爽女人下面视频在线观看| 亚洲伊人久久精品综合| 亚洲在久久综合| 熟女人妻精品中文字幕|