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

    震后應急物資配送的模糊定位-路徑問題研究

    2016-12-28 02:04:34劉長石
    中國管理科學 2016年5期
    關鍵詞:需求量約束條件救災

    劉長石,彭 怡,寇 綱

    (1.電子科技大學經濟與管理學院,四川 成都 611731; 2.湖南商學院工商管理學院,湖南 長沙 410205;3.西南財經大學工商管理學院,四川 成都 610074;4.湖南商學院湖南省移動電子商務協(xié)調創(chuàng)新中心,湖南 長沙 410205;5.湖南商學院移動商務智能湖南省重點實驗室,湖南 長沙 410205)

    ?

    震后應急物資配送的模糊定位-路徑問題研究

    劉長石1,2,4,5,彭 怡1,寇 綱3

    (1.電子科技大學經濟與管理學院,四川 成都 611731; 2.湖南商學院工商管理學院,湖南 長沙 410205;3.西南財經大學工商管理學院,四川 成都 610074;4.湖南商學院湖南省移動電子商務協(xié)調創(chuàng)新中心,湖南 長沙 410205;5.湖南商學院移動商務智能湖南省重點實驗室,湖南 長沙 410205)

    從系統(tǒng)集成優(yōu)化的角度研究震后應急物資配送的一種新模糊定位-路徑問題(LRP),綜合考慮救災點所在地理位置和地形導致的應急車輛行駛時間的隨機性、救災點應急物資需求量的不確定性與應急物資配送的時間緊迫性,以應急物資總運達時間最短與總配送成本最小為目標,構建一個基于機會約束規(guī)劃的多目標模糊LRP優(yōu)化模型,并根據模型的特點設計了一種混合免疫遺傳算法予以求解。最后,通過算例驗證了本文方法能有效解決震后應急物資配送的模糊多目標LRP,實現(xiàn)了震后應急物流中心定位和應急車輛路徑規(guī)劃的聯(lián)合決策。

    應急物資配送;車輛隨機行駛時間;定位-路徑問題;混合免疫遺傳算法

    1 引言

    近年來地震災害頻發(fā),給人類帶來了巨大的人員傷亡和財產損失。例如2008年5月12日,中國四川省汶川縣發(fā)生里氏7.8級地震,約6.9萬人遇難;2010年1月12日,海地發(fā)生7.3級地震,超過27萬人死亡;2015年4月25日,尼泊爾發(fā)生8.1級地震,造成重大人員傷亡。為盡可能降低地震損失和人員傷亡,應急物資必須在震后盡可能短的時間內配送到各救災點。通常,震后應急物資配送首先考慮在有限的時間、空間和資源等各種約束條件下將應急物資從各級應急物資儲備庫運往設置在災區(qū)附近的各應急物流中心,然后再由各應急物流中心配送至各救災點[1]。因此,應急物流中心定位(Location Allocation Problem, LAP)與應急車輛路徑規(guī)劃(Vehicle Routing Problem, VRP)是震后應急物資配送的兩個關鍵問題。有部分學者采用分別解決LAP與VRP的方法求解LRP[2-3]。事實上,LAP與VRP兩者之間存在相互依賴、相互影響的關系。例如,應急物資配送中的不同LAP方案導致不一樣的VRP規(guī)劃,VRP規(guī)劃結果反過來也會影響LAP方案[4-6]。因此,必須綜合考慮應急物資配送中的LAP和VRP的聯(lián)合決策問題,即定位-路徑問題(Location Routing Problem, LRP)。

    當前,關于商業(yè)物資配送中的LRP研究成果比較多[7-8],但是關注應急物資配送中的LRP研究文獻相對較少。代穎等[9]考慮救援時間范圍模糊的LRP,構建了一個以救援時間滿意度之和最大、應急物流系統(tǒng)總成本最小為目標的多目標LRP優(yōu)化模型,并設計了一種混合遺傳算法予以求解。孫華麗等[10]考慮了不確定應急物資需求與路網通行能力風險,以系統(tǒng)總成本最小為目標,建立了一個應急物流LRP模型,并設計了一種改進遺傳算法求解該問題。Sun Huali等[11]采用雙層規(guī)劃理論,以救援時間滿意度最大與系統(tǒng)總成本最小為目標進行了應急物流中的LRP規(guī)劃。Han Yunjun等[12]考慮時間窗和道路容量限制的LRP,采用拉格朗日松弛法將其轉化為LAP與VRP分別予以解決。李雙林等[13]研究了震后初期應急物資配送的模糊多目標LRP,以應急物資總配送時間最短和受災點應急物資未滿足的總損失最小為目標建立了一個多目標LRP優(yōu)化模型,并設計了一種遺傳算法予以求解。Wang Haujun等[14]研究了震后開放式應急物資供應的LRP,構建了一個非線性整數(shù)LRP優(yōu)化模型,并設計了一種混合啟發(fā)式算法予以求解。最近還有馬祖軍等[15]研究了帶限制期的震后應急物資配送模糊多目標開放式LRP, Rath與Gutjahr[16]研究了災后應急物資配送中的LRP。

    已有研究大多假設救災點的應急物資需求量小于車輛容量、應急物流中心與救災點分屬兩個不同的集合、應急物流網絡始終正常連通、車輛在各路段上的行駛速度固定,甚少考慮到救災點所處的地理位置和地形對車輛行駛時間的影響,同時鮮少考慮災后的道路橋梁損毀與路網動態(tài)恢復等現(xiàn)實情況,難以體現(xiàn)震后應急物資配送的真實特性。在實際的震后應急物資配送過程中,重災區(qū)可能有部分救災點的應急物資需求量非常大;由于地震災害的破壞性,應急物流網絡可能存在不同損毀程度的部分道路,致使配送車輛通行不暢,導致行駛在損毀道路上的車輛的行駛時間難以確定;救災點所處的地理位置和地形、單位作業(yè)單元所需修復時間、道路性質及恢復難易程度等等因素對車輛行駛時間都具有非常大的影響[18-21]。因為與地震中心距離不同的道路的損毀情況是不同的,不同地形條件下的道路損毀后的搶修難度是有區(qū)別的。例如,距離地震中心越近,道路損毀情況越嚴重;山區(qū)道路損毀后的搶修難度要遠遠大于平原道路損毀后的搶修難度。

    基于以上分析,本文研究一種新的LRP:假設震后有n個救災點,在這n個救災點內選擇m個點作為沒有容量限制的應急物流中心,即M?N,并綜合考慮了救災點所處地理位置和地形導致的車輛行駛時間的隨機性、應急物資需求的不確定性、限制期、車輛多次往返配送、同時進行車輛“滿載直配”與“巡回配送”等限制條件,采用機會約束規(guī)劃理論并以應急物資總運達時間最短和總配送成本最小為目標,構建基于機會約束規(guī)劃的震后應急物資配送多目標模糊LRP優(yōu)化模型,并根據模型的特點設計一種混合免疫遺傳算法予以求解。

    2 多目標模糊LRP數(shù)學模型

    2.1 問題描述

    震后急需在盡可能短的時間、有限的空間和資源約束條件下將應急物資配送到各個救災點。本文假設地震中心的位置已知,有n個地理位置和地形已知的救災點,救災點需求量為一個三角模糊數(shù),在這n個救災點內選擇m個點作為沒有容量限制的應急物流中心(該中心必須服務屬于自己的那個救災點),各應急物流中心擁有多輛不同容量限制的車輛,車輛必須從應急物流中心出發(fā),完成任務后返回到出發(fā)點,要求在一定的時間期限內完成所有救災點的應急物資配送任務。決策問題:確定應急物流中心選址、并規(guī)劃各物流中心的服務其所分配救災點的車輛路徑,在規(guī)定時間期限內滿足各救災點的應急物資需求,使應急物資總運達時間最短、總配送成本最小。

    災區(qū)可能有部分救災點的應急物資需求量非常大,超過車輛容量。因此本文將救災點分為大需求量救災點與小需求量救災點(應急物資需求量大于車輛容量的救災點為大需求量救災點,否則為小需求量救災點)。關于大需求量救災點采用“需求分割”策略同時進行滿載直配與巡回配送。即把大需求量救災點的應急物資需求分為兩部分,先對該大需求量救災點滿足滿載直接配送的應急物資需求處理為一個虛擬大需求量救災點,再對該大需求量救災點剩余不足最小車輛容量的應急物資需求處理為一個虛擬小需求量救災點,和其他小需求量救災點一起進行巡回配送[15]。

    2.2 符號說明

    定義如下符號和變量:

    2.3 模型構建

    針對救災點的模糊需求量必須滿足車輛容量限制這一約束條件,本文引進決策者置信水平這一概念,采用概率方法將其處理為機會約束條件,構建震后應急物資配送的多目標LRP模型如下:

    (1)

    (2)

    (3)

    (4)

    (5)

    (6)

    (7)

    xijk=0,?(i,j)∈M,?k∈V

    (8)

    (9)

    (10)

    Tkj=Tki+tijkxijk,?(i,j)∈AN,k∈V

    (11)

    Tki≤Lti,?i∈AN,k∈V

    (12)

    xijk(Tkj-Tki)≥0, ?(i,j)∈AN,k∈V

    (13)

    yg={0,1}, ?g∈M

    (14)

    βgi={0,1},?g∈M,i∈AN

    (15)

    wk={0,1},?k∈V

    (16)

    xijk={0,1}, ?(i,j)∈AN,k∈V

    (17)

    目標函數(shù)(1)表示最小化應急物資總運達時間;目標函數(shù)(2)表示最小化應急物資總配送成本(物流中心固定成本、車輛滿載直配成本、車輛巡回配送成本、車輛使用固定成本)。約束條件(3)保證分配給巡回配送車輛的所有小需求量救災點的應急物資需求量之和不超過該車容量的可能性不低于預先確定的置信水平;約束條件(4)表示運到大需求量救災點的應急物資要滿足該救災點的需求,其中rem(di,Qk)表示該救災點需求量與應急車輛容量的余數(shù);約束條件(5)保證開設的應急物流中心有車發(fā)出;約束條件(6)保證未開設的應急物流中心不能發(fā)車;約束條件(7)為車輛路徑連續(xù)條件,即進入救災點的車輛必須從該救災點離開;約束條件(8)保證任意2個應急物流中心不直接相連;約束條件(9)保證應急車輛至多只能分配給一個應急物流中心;約束條件(10)保證每個救災點有且僅有一輛車進行巡回配送服務;約束條件(11)給出車輛到達救災點的時間計算方法;約束條件(12)表示救災點的時間窗約束;約束條件(13)保證車輛路線時間的順序性;約束條件(14)至(17)為整數(shù)變量約束。其中,ceil(fbgi/Qk)表示不小于fbgi/Qk的整數(shù)。

    2.4 模糊需求量約束條件處理

    由于約束條件(3)具有模糊參數(shù)di,無法直接計算,需要進行適當處理。本文綜合參考了代穎等[6,13,15,17]的模糊需求處理方法。令Ek表示應急車輛k配送完所分配救災點后的剩余容量,Ek=(Ek1,Ek2,Ek3),即:

    (18)

    (19)

    故約束條件(3)等價于

    Pos(Ek≥0)≥ρ,?k∈V

    (20)

    2.5 隨機行駛時間約束條件處理

    研究應急車輛行駛時間的文獻甚少。僅僅孫華麗等[10]考慮了應急救援過程中的超期風險、路段通行能力風險、路徑復雜性等因素下的車輛行駛時間;Shinozuka等[18]和Bertini等[19]都把道路容量和車輛數(shù)量作為影響車輛行駛的關鍵因素,但都只分析了道路正常狀態(tài)下的車輛行駛時間;蘇幼坡[20]根據災后道路需要修復的作業(yè)單元數(shù)、道路性質及恢復難易程度來預計車輛行駛時間;Peng Min等[21]對震后車輛行駛時間進行了模擬仿真。

    (21)

    3 混合免疫遺傳算法設計

    LRP的求解算法主要分為兩階段求解算法與整體求解算法。前者運行時間較短,后者運行時間相對稍長,但是求得的LRP規(guī)劃方案較優(yōu)[7-9]。免疫遺傳算法不僅保留了遺傳算法的搜索特性,且具有免疫算法的多機制求解多目標函數(shù)最優(yōu)解的自適應特性,能有效避免算法早熟與收斂于局部極值[22,23]。因此,本文基于整體求解的思路,設計了一個混合免疫遺傳算法求解LRP。

    3.1 算法設計

    3.1.1 初始抗體群產生

    將LRP的目標函數(shù)及約束條件對應于抗原,問題的解對應于抗體,按照預先設定的群體規(guī)模H,隨機生成一組抗體作為初始抗體群B。

    本文中的抗體由三個子串組成:子串1有n個基因位(n為救災點數(shù)目),每個基因位的取由1到n的自然數(shù)隨機排列生成且相互不能重復的值;子串2有k個基因位,每個基因位的取由1到m的自然數(shù)中隨機生成且相互不能重復的值(m為候選應急物流中心數(shù)目);子串3有n個基因位,每個基因位的取值由1到k的自然數(shù)中隨機選取(k為應急車輛數(shù)目)。其中,子串1表示路徑中各需求點的排列順序,子串2表示各應急車輛歸屬哪一個應急物流中心,子串3和子串1中的基因位一一對應,表示為各救災點服務的車輛編號。

    3.1.2 抗體適應度與親和力計算

    本文LRP模型的目標函數(shù)與所有約束條件均為線性函數(shù),可行域是凸集,根據文獻[15,25]可知,凸性多目標規(guī)劃問題可以用加權法處理。因此,抗體適應度函數(shù)可通過對目標函數(shù)值進行變換得到。

    (22)

    而且fitness(b)的值越大,抗體b與抗原的親和力G(b)越大,表明該抗體越優(yōu)。

    3.1.3 抗體促進與抑制

    在抗體進化過程中,選擇親和力較高而且濃度較低的抗體進行促進。同時,為確??贵w群體更新過程中抗體的多樣性、避免未成熟收斂,本文規(guī)定當抗體的濃度達到一定的閥值時隨機生成超濃度規(guī)模數(shù)量的若干新抗體替換舊抗體。根據抗體的編碼特點,本文定義抗體b濃度為:

    (23)

    其中,

    如果種群中一共有ε種不同的親和力,設抗體濃度閥值為ε/H。

    3.1.4 遺傳操作

    采用傳統(tǒng)遺傳算法的3基本操作方法,即選擇、交叉和變異,這里不再累述。

    3.1.5 免疫算子

    免疫算子一般包括提取疫苗、接種疫苗和免疫選擇等3個步驟。

    (1)提取疫苗。抽取當前k+1代的最優(yōu)個體群中的最優(yōu)個體,同時找到前k代保留下來的最優(yōu)個體群中的最優(yōu)個體,把這2條最優(yōu)個體的基因位的共同特點和有效信息作為疫苗。

    (2)接種疫苗。抽取一定概率的帶接種的抗體,進行疫苗接種。

    (3)免疫選擇。對已接種疫苗的個體進行檢測,若子代最高適應值的個體優(yōu)于父代,子代替代父代進入下一種群,否則采用賭輪方式選擇父代中的個體進入下一種群。

    3.1.6 終止條件

    設算法最大迭代次數(shù)為M, 當抗體群的進化代數(shù)達到了M時算法終止。

    3.2 算法復雜度分析

    本文算法采用三參數(shù)級聯(lián)編碼方式,第一子串長度為n,第二子串長度為m,第三子串長度為n,即每條抗體占用的內存空間為2n+m。如果增加b個救災點,抗體長度為2(n+b)+m+ω,0≤ω≤b,則每條抗體內存空間的增量為O(2b+ω)。又O(2b)≤O(2b+ω)≤O(3b),即隨著算例增大本文算法占用內存空間呈線性增加。

    本文算法的群體規(guī)模為H,最大迭代次數(shù)為M,問題規(guī)模為Q=n+n+k。算法計算時間復雜度分析如下:初始抗體群產生后,每次迭代都要經過適應度計算、抗體促進與抑制、遺傳操作、提取疫苗、接種疫苗和免疫選擇等計算?;旌厦庖哌z傳算法的時間復雜度為:

    Ο(H,M,n,n,k)=HΟ(2n+k)+M(Ο(Hn2k)+Ο(HlgQ)+Ο(HlgQ)+Ο(Hn2k)+Ο(Hn2kλ)+Ο(Hn2kλ))≈2HM(Ο(n2k))

    其中,λ表示進行疫苗接種抗體的概率。

    可見,本文算法的計算量與抗體群規(guī)模、迭代次數(shù)、需求點數(shù)量與車輛數(shù)量呈正比。

    4 算例分析

    以Solomon的VRPTW BENCHMARK PROBLEMS測試數(shù)據庫中的RC201算例數(shù)據為基礎,同時根據本文所研究LRP的特點對需求數(shù)據進行修改并補充相關數(shù)據,產生測試算例。具體方法如下:從RC201數(shù)據中隨機抽取26個坐標數(shù)據作為救災點的地理位置,救助點之間的距離采用直線距離表示,然后隨機產生各救災點的物資需求數(shù)量與地形。候選應急物流中心的是從這26個救災點內選取使目標函數(shù)值最小的3個點。設定各應急物流中心具有2種不同容量類型的車輛。各相關數(shù)據如表1與表2所示。

    相關變量設置如下:群體規(guī)模為40,記憶庫容量為10,最大迭代次數(shù)為500,交叉概率為0.5,變異概率為0.4,多樣性評價函數(shù)值為0.95,應急物流中心數(shù)量為3,固定使用成本為20000元/個,置信水平為0.9,時間目標函數(shù)值的權重w1=0.6,成本目標函數(shù)值權重w2=0.4,車輛容量分別為280件與350件,車型1的行駛速度為60km/h,車型2行駛速度為50km/h,平均單位行駛成本為10元/km,地震中心的坐標為(50,50),Tfp=3,α=5,χ=4,β=3,φ=2。

    采用Matlab R2013a對本文設計的混合免疫遺傳算法進行編程實現(xiàn),在AMD A8-4500m APU 1.90GHz、內存4G的微機上對測試算例進行了20次求解。程序運行時間比較短,計算效率比較高,最短計算時間為37.72秒,平均計算時間為38.74秒,最大計算時間為40.76秒。計算結果也比較穩(wěn)定:20次方案都是啟用了12臺容量為280件的應急車輛,目標函數(shù)值1的最優(yōu)解為3359.78分鐘,最差解為3431.46分鐘,前者與后者的差距只有2.09%;目標函數(shù)值2的最優(yōu)解為75093.1元,最差解為76130.6元,前者與后者的差距只有1.36%。選擇的3個應急物流中心坐標是:(2,45)、(65,20)、(45,65)。各應急物流中心的最佳車輛行駛路線如圖1所示。具體應急物流中心車輛路線安排見表3,車輛路線安排中的數(shù)字0表示應急物流中心,其余數(shù)字表示救災點序號。從圖1和表3可知,本文的算法有效實現(xiàn)了應急物流中心定位與車輛路徑規(guī)劃的聯(lián)合決策,合理進行了救災點的滿載直配與巡回配送,在規(guī)定時間內完成了應急物資配送任務。

    為進一步說明本文方法的有效性,引進阮俊虎等[24]的2個考核指標:平均等待時間、最大等待時間。平均等待時間表示應急物資到達救災點的平均時間,越短意味著應急物資(尤其是醫(yī)療器械與藥品)的效用越大。最大等待時間在一定程度上反映了應急物資到達救災點時間的差距,越短表示應急物資配送越公平[24]。經過計算,本文的混合免疫遺傳算法的平均等待時間為108.01分、最大等待時間為259.36分;Yi Wei[3]的蟻群算法的平均等待時間為110.79分、最大等待時間為282.52分。前者的2個考核指標比后者分別降低1.81%與8.19%。說明本文的混合免疫遺傳算法具有良好的性能。

    在抗體群規(guī)模不變的情況下,通過擴展算例1(表1)形成算例2-4,采用本文的混合免疫遺傳算法(HIGA)與Yi Wei[3]的蟻群算法(ACA)分別計算各個算例,計算結果如表4所示。表4中,n表示救災點數(shù)量,m表示物流中心數(shù)量,k表示車輛數(shù)量,計算結果為求解20次的均值。從表4可知,隨著問題規(guī)模的增大,程序運行時間會有所增加,與前面分析得出的混合免疫遺傳算法的時間復雜度2HM(Ο(n2k))一致,且HIGA的計算時間稍優(yōu)于ACA的計算時間。

    表1 救災點數(shù)據

    表2 應急車輛數(shù)據

    圖1 車輛路線安排示意圖

    表3 應急物流中心車輛路線安排

    表4 不同問題規(guī)模計算時間

    5 結語

    震后,為盡可能降低地震帶來的人員傷亡與財產損失,應急物資必須在盡可能短的時間內配送到各個救災點。研究震后應急物流系統(tǒng)、提高應急物資配送績效,具有非常大的現(xiàn)實價值。本文從系統(tǒng)集成的角度研究了一種新的LRP,并綜合考慮了救災點與地震中心的距離、救災點的地理位置與地形等因素對車輛行駛時間的影響,構建了一個以應急物資總運達時間最短與總配送成本最小為目標的多目標模糊LRP優(yōu)化模型,并根據模型的特點設計了一種混合免疫遺傳算法予以求解,通過算例驗證了本文模型的正確性與算法的穩(wěn)定性、有效性。

    本文構建的模型主要有以下適用條件:在震后應急救援初期,應急物資(尤其是醫(yī)療器械與藥品)需求不確定、即來即走,從各級應急物資儲備庫運達各應急物流中心后,馬上進行配送安排,救災點要求應急物資在規(guī)定時間內運達,物流中心的車輛與裝卸能力充足,應急物流網絡存在部分損毀道路;不適合于應急物流中心的應急物資有一定的存儲時間這種情況。本文只考慮了單品種應急物資配送的LRP,而現(xiàn)實中的應急物資配送可能涉及多種物資,因此多品種應急物資配送的LRP有待進一步研究。另外,在某些情況下,應急物流具有動態(tài)調度的特性,應急車輛路線動態(tài)優(yōu)化也將是下一個研究重點。

    [1] 鄭斌,馬祖軍,方濤.應急物流系統(tǒng)中的模糊多目標定位-路徑問題[J].系統(tǒng)工程, 2008,27(8):21-25.

    [2] Chang M S, Tseng Y L, Chen Jingwen.A scenario planning approach for the flood emergency logistics preparation problem under uncertainty[J].Transportation Research E, 2007,43(6): 737-754.

    [3] Yi Wei, Kumar A.Ant colony optimization for disaster relief operations[J].Transportation Research E, 2007,43(6): 660-672.

    [4] 曾敏剛,崔增收,余高輝.基于應急物流的減災系統(tǒng)LRP研究[J].中國管理科學, 2010, 18(2):75-80.

    [5] 代穎,馬祖軍.應急物流系統(tǒng)中的隨機定位-路徑問題[J].系統(tǒng)管理學報, 2012,21(2): 212-218.

    [6] 代穎,馬祖軍,朱道立,等.震后應急物資配送的模糊動態(tài)定位-路徑問題[J].管理科學學報, 2012,15(7): 212-218.

    [7] 汪壽陽,趙秋紅,夏國平.集成物流管理系統(tǒng)中的定位-運輸路線安排問題研究[J].管理科學學報, 2000,3(2): 69-75.

    [8] 林巖,胡祥培,王旭茵.物流系統(tǒng)優(yōu)化中的定位-運輸路線安排問題(LRP)研究評述[J].管理工程學報, 2004,18(4): 45-49.

    [9] 代穎,馬祖軍,鄭斌.突發(fā)公共事件應急系統(tǒng)中的模糊多目標定位-路徑問題研究[J].管理評論, 2010,22(1): 121-128.

    [10] 孫華麗,周戰(zhàn)杰,薛耀鋒.考慮路徑風險的不確定需求應急物流定位-路徑問題[J].上海交通大學學報, 2013,47(6): 962-966.

    [11] Sun Huali, Wang Xunqing, Xue Yaofeng.A bi-level programming model for a multi-facility location-routing problem in urban emergency system[M]//Zhang Liangchi,Zhang Chunliang.Lecture Notes in Electrical Engineering, 2011,111: 75-80.

    [12] Han Yunjun, Guan Xiahong, Shi Leyuan.Optimization based method for supply location selection and routing in large scale emergency material delivery[J].IEEE Transactions on Automation Science and Engineering, 2011,8(4): 683-693.

    [13] 李雙琳,馬祖軍,鄭斌,等.震后初期應急物資配送的模糊多目標選址-多式聯(lián)運問題[J].中國管理科學, 2013,21(2): 144-151.

    [14] Wang Haijun, Du Lijing, Ma Shihua.Multi-objective open location-routing model with split delivery for optimized relief distribution in post-earthquake [J].Transportation Research Part E, 2014,69: 160-179.

    [15] 馬祖軍,代穎,李雙琳.帶限制期的震后應急物資配送模糊多目標開放式定位-路徑問題[J].系統(tǒng)管理學報,2014,23(5): 658-667.

    [16] Rath S, Gutjahr W J.A math-heuristic for the warehouse location-routing problem in disaster relief[J].Computers & Operations Research, 2014,42(2): 25-39.

    [17] 張建勇,李軍,郭耀煌.模糊需求信息條件下的實時動態(tài)車輛調度問題研究[J].管理工程學報,2004,18(4): 69-72.

    [18] Shinozuka M, Murachi Y, Dong Xuejiang, et al.Effect of seismic retrofit of bridges on transportation networks[J].Earthquake Engineering and Engineering Vibration, 2003,2(2): 169-179.

    [19] Bertini R L, Hansen S, Byrd A,et al.Experience implementing a user service for archived intelligent transportation systems data[J].Transportation Research Record: Journal of the Transportation Research Board, 2005,1917(1): 90-99.

    [20] 蘇幼坡.城市生命線系統(tǒng)震后恢復過程優(yōu)化的研究[D].天津:天津大學,2001.

    [21] Peng Min, Peng Yi, Chen Hong.Post-seismic supply chain risk management: A system dynamics disruption analysis approach for inventory and logistics planning[J].Computers & Operations Research, 2014,42(1): 14-24.

    [22] Molla-Alizadeh-Zavardehi S, Hajiaghaei-Keshteli M, Tavakkoli-Moghaddam R.Solving a capacitated fixed-charge transportation problem by artificial immune and genetic algorithms with a prufer number representation [J].Expert Systems with Applications, 2011,38(8): 10462-10474.

    [23] Chui C Y, Kuo I T, Lin C H.Applying artificial immune system and ant algorithm in air-conditioner market segmentation[J].Expert Systems with Applications, 2009,36(3): 4437-4442.

    [24] 阮俊虎,王旭坪,楊挺.大規(guī)模災害中基于聚類的醫(yī)療物資聯(lián)合運送優(yōu)化[J].中國管理科學, 2014,22(10): 80-89.

    [25] Ehrgott M,Shao Lizhen,Schobel A.An approximation algorithm for convex multi-objective programming problems[J].Journal of Global Optimization, 2011,50(3): 397-416.

    [26] Leyffer S.A complementarity constraint formulation of convex multiobjective optimization problems[J].Informs Journal on Computing, 2009,21(2): 257-267.

    Research on Fuzzy Location-routing Problem in Post-earthquake Delivery of Relief Materials

    LIU Chang-shi1,2,4,5, PENG Yi1, KOU Gang3

    (1.School of Management and Economics, University of Electronic Science and Technology of China, Chengdu 611731, China;2.School of Management, Hunan University of Commerce, Changsha 410205,China;3.School of Business Administration, Southwestern University of Finance and Economics, Chengdu 610074, China;4.Mobile E-business Collaborative Innovation Center of Hunnan Province,Hunan University of Commerce,Changsha 410205, China; 5.Key Laborartory of Huna Province for Mobile Bsiness Intelligence,Hunan University of Commerce,Changsha 410205,China)

    In order to improve the efficiency of emergency logistics, a new location-routing problem (LRP) in post-earthquake delivery of relief materials was studied from the view point of integrated optimization.By considering the stochastic vehicle time caused by the location and topography of relief points, the urgent window constraints, the fuzzy demands of relief materials and the urgency of rescue time in post-earthquake, a chance-constrained programming model for the fuzzy LRP was developed, and the goal is to minimize the total time in delivering relief materials and the total costs.A hybrid immune genetic algorithm was proposed to solve the model.Finally,the feasibility and validity of the model and algorithm was demonstrated by a numerical example.The test data are obtained from the Solomon’s VRPTW BENCHMARK PROBLEMS.The results show that the proposed approaches effectively solve the fuzzy multi objective LRP, and realize joint decision of emergency logistics center location and emergency vehicle routing planning.

    delivery of relief materials;stochastic vehicle travel time;location-routing problem;hybrid immune genetic algorithm

    1003-207(2016)05-0111-08

    10.16381/j.cnki.issn1003-207x.2016.05.013

    2014-08-17;

    2015-10-10

    國家自然科學基金面上資助項目(71173028);國家自然科學基金資助項目(71222108);湖南省教育廳優(yōu)秀青年項目(15B131);湖南省科技廳重點研發(fā)計劃項目(2015ZK3049)

    簡介:劉長石(1975-),男(漢族),湖南邵陽人,電子科技大學經濟與管理學院博士究生,湖南商學院教師,研究方向:應急物流管理、供應鏈管理,E-mail:liuchangshi964@126.com.

    F252

    A

    猜你喜歡
    需求量約束條件救災
    航母救災
    基于一種改進AZSVPWM的滿調制度死區(qū)約束條件分析
    從數(shù)學角度看“彈性”
    防汛救災 應急先鋒沖在前
    滅火救災分秒不能耽誤
    水上消防(2020年3期)2020-07-25 02:36:26
    A literature review of research exploring the experiences of overseas nurses in the United Kingdom (2002–2017)
    全省齊心協(xié)力共打抗洪搶險救災攻堅戰(zhàn)
    河北水利(2016年8期)2017-01-13 12:39:00
    線性規(guī)劃的八大妙用
    2017年我國汽車軟管需求量將達6.4億m
    橡膠科技(2015年3期)2015-02-26 14:45:02
    基于BP神經網絡人均豬肉需求量預測
    成熟少妇高潮喷水视频| 很黄的视频免费| www.熟女人妻精品国产| 老司机在亚洲福利影院| 国产精品av视频在线免费观看| 99精品欧美一区二区三区四区| 成人午夜高清在线视频| 欧美日韩亚洲国产一区二区在线观看| 国产伦在线观看视频一区| 成人特级黄色片久久久久久久| 在线视频色国产色| 岛国在线免费视频观看| e午夜精品久久久久久久| 窝窝影院91人妻| 亚洲中文字幕一区二区三区有码在线看 | 中文字幕最新亚洲高清| 婷婷精品国产亚洲av| 美女免费视频网站| 国产成人精品久久二区二区91| 桃色一区二区三区在线观看| 一本大道久久a久久精品| 此物有八面人人有两片| 哪里可以看免费的av片| 操出白浆在线播放| 亚洲成av人片在线播放无| 午夜福利高清视频| 一个人免费在线观看的高清视频| 在线观看日韩欧美| 国产aⅴ精品一区二区三区波| 免费无遮挡裸体视频| 精华霜和精华液先用哪个| 天天添夜夜摸| 久久精品国产99精品国产亚洲性色| 亚洲成av人片在线播放无| 嫩草影视91久久| 麻豆国产97在线/欧美 | 在线观看www视频免费| 99久久精品热视频| 禁无遮挡网站| 中文字幕久久专区| 日韩欧美国产一区二区入口| 国产亚洲精品久久久久5区| 日本撒尿小便嘘嘘汇集6| 一级毛片女人18水好多| 欧美一区二区精品小视频在线| 国产高清视频在线观看网站| 亚洲人成电影免费在线| 精品久久久久久久人妻蜜臀av| 午夜视频精品福利| 国产精品爽爽va在线观看网站| 别揉我奶头~嗯~啊~动态视频| 91老司机精品| 亚洲 欧美 日韩 在线 免费| 麻豆国产97在线/欧美 | 久久精品人妻少妇| 精品国产乱码久久久久久男人| 欧美色欧美亚洲另类二区| 看黄色毛片网站| 老熟妇仑乱视频hdxx| 变态另类丝袜制服| 免费在线观看亚洲国产| 亚洲欧美日韩东京热| 一级a爱片免费观看的视频| 免费观看人在逋| 久久草成人影院| 久久精品亚洲精品国产色婷小说| 国产精品av视频在线免费观看| 亚洲精华国产精华精| www.精华液| 曰老女人黄片| 在线免费观看的www视频| 国产精华一区二区三区| 国产精品野战在线观看| 日韩欧美国产一区二区入口| 亚洲性夜色夜夜综合| aaaaa片日本免费| 91在线观看av| 女人高潮潮喷娇喘18禁视频| 国产成人精品久久二区二区91| 亚洲avbb在线观看| 亚洲欧美精品综合一区二区三区| 国产高清视频在线观看网站| aaaaa片日本免费| 欧美乱色亚洲激情| 久久这里只有精品中国| 香蕉国产在线看| xxx96com| 国产精品影院久久| 亚洲九九香蕉| 高清在线国产一区| 国产爱豆传媒在线观看 | 国产精品一区二区三区四区免费观看 | 精品乱码久久久久久99久播| 久久国产乱子伦精品免费另类| 国产一区在线观看成人免费| 美女高潮喷水抽搐中文字幕| 琪琪午夜伦伦电影理论片6080| 国产精品av久久久久免费| 看黄色毛片网站| 婷婷亚洲欧美| 少妇人妻一区二区三区视频| 国产精品免费一区二区三区在线| 久久久久久国产a免费观看| 国模一区二区三区四区视频 | 欧美 亚洲 国产 日韩一| 中文字幕人成人乱码亚洲影| 久久天躁狠狠躁夜夜2o2o| 91大片在线观看| 久久久久久免费高清国产稀缺| 国产av不卡久久| 久9热在线精品视频| 看片在线看免费视频| 女人高潮潮喷娇喘18禁视频| 日韩高清综合在线| 色综合婷婷激情| 2021天堂中文幕一二区在线观| 国产亚洲精品久久久久5区| 一卡2卡三卡四卡精品乱码亚洲| 免费无遮挡裸体视频| 国产私拍福利视频在线观看| 在线看三级毛片| 欧美性长视频在线观看| 深夜精品福利| 国产69精品久久久久777片 | 欧美黄色片欧美黄色片| 1024手机看黄色片| 日韩国内少妇激情av| 午夜a级毛片| 国产成人影院久久av| 人妻丰满熟妇av一区二区三区| 国产99白浆流出| 国内少妇人妻偷人精品xxx网站 | 哪里可以看免费的av片| 欧美日韩亚洲综合一区二区三区_| 久久久久九九精品影院| 国产成人欧美在线观看| 免费在线观看日本一区| 变态另类成人亚洲欧美熟女| 国产精品影院久久| 中文字幕人妻丝袜一区二区| 久久人妻福利社区极品人妻图片| tocl精华| 国产精品 欧美亚洲| 亚洲国产欧美网| 欧美三级亚洲精品| 亚洲美女视频黄频| 欧美日韩亚洲综合一区二区三区_| 久久久精品国产亚洲av高清涩受| 国产1区2区3区精品| 此物有八面人人有两片| 免费在线观看完整版高清| 1024视频免费在线观看| √禁漫天堂资源中文www| 欧美乱码精品一区二区三区| 成人国产综合亚洲| 国内少妇人妻偷人精品xxx网站 | 色噜噜av男人的天堂激情| 午夜免费观看网址| 18禁美女被吸乳视频| 色综合欧美亚洲国产小说| 一区二区三区国产精品乱码| 12—13女人毛片做爰片一| 丰满的人妻完整版| 色播亚洲综合网| 两人在一起打扑克的视频| 亚洲国产欧美网| 国产男靠女视频免费网站| 88av欧美| 欧美成人一区二区免费高清观看 | 一本精品99久久精品77| 中文资源天堂在线| 精品电影一区二区在线| 亚洲av美国av| 免费观看人在逋| 可以在线观看的亚洲视频| 国内久久婷婷六月综合欲色啪| 欧美 亚洲 国产 日韩一| 97碰自拍视频| 长腿黑丝高跟| 在线观看午夜福利视频| 99久久综合精品五月天人人| 久久精品国产综合久久久| 成人av一区二区三区在线看| 88av欧美| 国产精华一区二区三区| 国内久久婷婷六月综合欲色啪| 亚洲午夜精品一区,二区,三区| 国产精品久久久av美女十八| 一夜夜www| 成人国语在线视频| 久久久久久人人人人人| 最近最新免费中文字幕在线| 韩国av一区二区三区四区| 久久婷婷人人爽人人干人人爱| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品98久久久久久宅男小说| 亚洲中文av在线| 欧美另类亚洲清纯唯美| 国产男靠女视频免费网站| 久久精品91无色码中文字幕| 久久久久免费精品人妻一区二区| 老汉色av国产亚洲站长工具| 婷婷丁香在线五月| 午夜免费成人在线视频| 制服诱惑二区| 变态另类成人亚洲欧美熟女| 亚洲精品中文字幕一二三四区| 99热这里只有是精品50| 在线国产一区二区在线| 久久久久久久久久黄片| 国产精品久久电影中文字幕| 久久精品人妻少妇| 久久中文字幕一级| 黄色丝袜av网址大全| 亚洲av成人不卡在线观看播放网| 一个人观看的视频www高清免费观看 | 国产成年人精品一区二区| 国产不卡一卡二| 亚洲 欧美一区二区三区| 成人永久免费在线观看视频| www.自偷自拍.com| 久久久水蜜桃国产精品网| 国产精品久久久久久人妻精品电影| 老司机午夜十八禁免费视频| 亚洲欧美日韩高清在线视频| 亚洲成人久久性| 国产午夜精品久久久久久| 1024手机看黄色片| 脱女人内裤的视频| 午夜福利高清视频| 久久久久久大精品| 国产成人一区二区三区免费视频网站| 免费在线观看日本一区| 一进一出抽搐gif免费好疼| 变态另类丝袜制服| 国产成人aa在线观看| 看片在线看免费视频| 99国产极品粉嫩在线观看| 一个人免费在线观看的高清视频| 久久国产精品人妻蜜桃| 99re在线观看精品视频| 97超级碰碰碰精品色视频在线观看| 欧美性长视频在线观看| 国产一区二区三区在线臀色熟女| 国产成人av激情在线播放| 亚洲欧美日韩无卡精品| 久久久国产成人免费| 超碰成人久久| 露出奶头的视频| 国产高清videossex| 久久欧美精品欧美久久欧美| 最近在线观看免费完整版| 亚洲av五月六月丁香网| 亚洲欧美一区二区三区黑人| 女人爽到高潮嗷嗷叫在线视频| 深夜精品福利| 午夜久久久久精精品| 无遮挡黄片免费观看| 久久久国产成人免费| 国语自产精品视频在线第100页| 一个人免费在线观看的高清视频| 国产精品电影一区二区三区| 一进一出好大好爽视频| 亚洲国产精品久久男人天堂| 日韩 欧美 亚洲 中文字幕| 国产成人av教育| 国产精品免费视频内射| 亚洲av第一区精品v没综合| 日日摸夜夜添夜夜添小说| 久久天堂一区二区三区四区| 亚洲精品国产一区二区精华液| 亚洲 国产 在线| 18禁裸乳无遮挡免费网站照片| 日韩国内少妇激情av| 亚洲av日韩精品久久久久久密| 久久国产精品人妻蜜桃| 午夜免费观看网址| 国内久久婷婷六月综合欲色啪| 村上凉子中文字幕在线| xxxwww97欧美| 人人妻,人人澡人人爽秒播| tocl精华| 国产一级毛片七仙女欲春2| 老鸭窝网址在线观看| 床上黄色一级片| 亚洲第一欧美日韩一区二区三区| 午夜免费成人在线视频| 精品无人区乱码1区二区| 欧美日韩黄片免| 最近最新免费中文字幕在线| 免费无遮挡裸体视频| 日韩欧美免费精品| 91在线观看av| 黄色视频,在线免费观看| 午夜精品一区二区三区免费看| 国产免费av片在线观看野外av| 少妇的丰满在线观看| www.自偷自拍.com| 黄色片一级片一级黄色片| 成人特级黄色片久久久久久久| 国产精品久久久久久人妻精品电影| 国产午夜精品论理片| 国产精品日韩av在线免费观看| 精品第一国产精品| 日韩欧美三级三区| 级片在线观看| 久久天堂一区二区三区四区| 久久久久精品国产欧美久久久| 国产又色又爽无遮挡免费看| 国产v大片淫在线免费观看| 伊人久久大香线蕉亚洲五| 后天国语完整版免费观看| 黄色成人免费大全| 琪琪午夜伦伦电影理论片6080| 成人三级做爰电影| 好看av亚洲va欧美ⅴa在| 日韩欧美 国产精品| 国产欧美日韩精品亚洲av| www国产在线视频色| 全区人妻精品视频| 亚洲18禁久久av| 精品熟女少妇八av免费久了| 亚洲一区二区三区不卡视频| avwww免费| 免费看十八禁软件| 国产探花在线观看一区二区| 欧美日本亚洲视频在线播放| av超薄肉色丝袜交足视频| 夜夜躁狠狠躁天天躁| 神马国产精品三级电影在线观看 | 日本一本二区三区精品| 亚洲 国产 在线| 岛国在线免费视频观看| 最新美女视频免费是黄的| 国产三级中文精品| 免费人成视频x8x8入口观看| 亚洲av中文字字幕乱码综合| 久久久水蜜桃国产精品网| 50天的宝宝边吃奶边哭怎么回事| 十八禁网站免费在线| 日韩欧美国产在线观看| 午夜免费观看网址| 国产精品一区二区免费欧美| 波多野结衣高清作品| 欧美久久黑人一区二区| www.999成人在线观看| 日韩欧美在线二视频| 在线观看一区二区三区| 久久精品国产99精品国产亚洲性色| 国产精品98久久久久久宅男小说| 天天添夜夜摸| 香蕉久久夜色| 午夜亚洲福利在线播放| 欧美日韩瑟瑟在线播放| 亚洲aⅴ乱码一区二区在线播放 | 真人一进一出gif抽搐免费| 一边摸一边做爽爽视频免费| 青草久久国产| 999精品在线视频| 午夜免费激情av| 在线观看66精品国产| 国产精品久久久久久人妻精品电影| 久久久久国产精品人妻aⅴ院| 狠狠狠狠99中文字幕| www日本在线高清视频| 亚洲18禁久久av| 国产高清videossex| 日韩欧美一区二区三区在线观看| 国内久久婷婷六月综合欲色啪| 久久婷婷人人爽人人干人人爱| 日本三级黄在线观看| 国产精品av视频在线免费观看| www.精华液| 久久婷婷人人爽人人干人人爱| 久久伊人香网站| 少妇被粗大的猛进出69影院| 一本精品99久久精品77| 亚洲国产欧洲综合997久久,| 99精品在免费线老司机午夜| 一个人免费在线观看的高清视频| 无限看片的www在线观看| 久久中文字幕人妻熟女| 免费在线观看视频国产中文字幕亚洲| 国产伦人伦偷精品视频| 欧美精品啪啪一区二区三区| 国产精品久久久久久人妻精品电影| 久久久久国产精品人妻aⅴ院| 99国产精品一区二区蜜桃av| 亚洲精品一卡2卡三卡4卡5卡| 母亲3免费完整高清在线观看| 老司机在亚洲福利影院| 黄片大片在线免费观看| 色哟哟哟哟哟哟| 国产真实乱freesex| 午夜视频精品福利| 国产视频内射| 亚洲国产精品999在线| 老司机福利观看| 亚洲精品久久成人aⅴ小说| 国产精品久久电影中文字幕| 丰满人妻一区二区三区视频av | 久久午夜亚洲精品久久| 日韩欧美 国产精品| 国产成人av教育| 男女午夜视频在线观看| 亚洲avbb在线观看| 免费电影在线观看免费观看| 国产精品电影一区二区三区| 极品教师在线免费播放| 日韩精品中文字幕看吧| 18美女黄网站色大片免费观看| 全区人妻精品视频| 国产成人aa在线观看| 女人高潮潮喷娇喘18禁视频| 校园春色视频在线观看| 亚洲激情在线av| 欧美一级a爱片免费观看看 | 国产精品98久久久久久宅男小说| 给我免费播放毛片高清在线观看| 日韩欧美 国产精品| 久久久精品国产亚洲av高清涩受| 美女大奶头视频| 亚洲av电影不卡..在线观看| 免费电影在线观看免费观看| 亚洲精品av麻豆狂野| 精品久久久久久久毛片微露脸| 韩国av一区二区三区四区| 亚洲欧美一区二区三区黑人| 欧美成人午夜精品| 最近在线观看免费完整版| 亚洲全国av大片| 亚洲成人久久爱视频| 国产三级在线视频| 日韩欧美在线二视频| 免费在线观看完整版高清| 国产激情欧美一区二区| 黄色视频,在线免费观看| 不卡av一区二区三区| 久久亚洲真实| 身体一侧抽搐| 国产亚洲欧美98| 日韩精品免费视频一区二区三区| 国产aⅴ精品一区二区三区波| 久久精品影院6| 人人妻人人看人人澡| 亚洲精品中文字幕一二三四区| 黄色成人免费大全| 一级作爱视频免费观看| 日韩欧美精品v在线| www.精华液| 亚洲av美国av| 黑人操中国人逼视频| 亚洲一码二码三码区别大吗| 欧美成人午夜精品| 狂野欧美白嫩少妇大欣赏| 欧美日本亚洲视频在线播放| 日本在线视频免费播放| 变态另类丝袜制服| 欧美 亚洲 国产 日韩一| 色精品久久人妻99蜜桃| 久久久久久九九精品二区国产 | 他把我摸到了高潮在线观看| 丁香六月欧美| 日本免费a在线| 变态另类成人亚洲欧美熟女| 操出白浆在线播放| www日本黄色视频网| www.熟女人妻精品国产| 国产三级黄色录像| 欧美日韩瑟瑟在线播放| 香蕉国产在线看| 午夜免费成人在线视频| 久久伊人香网站| 老汉色av国产亚洲站长工具| 一本大道久久a久久精品| 麻豆av在线久日| 国产精品综合久久久久久久免费| 免费看日本二区| 精品国产美女av久久久久小说| 夜夜躁狠狠躁天天躁| 黄色视频,在线免费观看| 免费无遮挡裸体视频| 制服诱惑二区| 日韩欧美精品v在线| 伦理电影免费视频| 90打野战视频偷拍视频| 777久久人妻少妇嫩草av网站| 香蕉国产在线看| 最好的美女福利视频网| 国产精品久久电影中文字幕| 99在线人妻在线中文字幕| 欧美乱色亚洲激情| 日本一二三区视频观看| 国产男靠女视频免费网站| 亚洲国产精品久久男人天堂| 国产精品九九99| 成人高潮视频无遮挡免费网站| 在线观看午夜福利视频| 亚洲av成人不卡在线观看播放网| 午夜成年电影在线免费观看| 久久久久久国产a免费观看| 成人高潮视频无遮挡免费网站| 91国产中文字幕| 久久欧美精品欧美久久欧美| 色尼玛亚洲综合影院| 国产三级黄色录像| 亚洲色图av天堂| 99精品久久久久人妻精品| 少妇人妻一区二区三区视频| 久久久久久久久免费视频了| 少妇裸体淫交视频免费看高清 | 久久精品国产亚洲av高清一级| 日韩中文字幕欧美一区二区| 国产亚洲欧美98| 一进一出好大好爽视频| 亚洲午夜理论影院| 这个男人来自地球电影免费观看| 婷婷精品国产亚洲av在线| 成人18禁在线播放| 91成年电影在线观看| 国语自产精品视频在线第100页| 亚洲五月婷婷丁香| 老汉色∧v一级毛片| 女警被强在线播放| 亚洲av成人不卡在线观看播放网| 在线观看日韩欧美| 欧美+亚洲+日韩+国产| av福利片在线| 免费在线观看影片大全网站| 久久香蕉激情| 白带黄色成豆腐渣| 亚洲黑人精品在线| 欧美乱色亚洲激情| 国产成人精品久久二区二区91| www日本在线高清视频| av视频在线观看入口| 欧美激情久久久久久爽电影| 在线观看美女被高潮喷水网站 | 国产亚洲精品一区二区www| 欧美日韩亚洲综合一区二区三区_| 韩国av一区二区三区四区| 叶爱在线成人免费视频播放| 欧美日本视频| 久久人妻福利社区极品人妻图片| 久久久久免费精品人妻一区二区| 国产精品av久久久久免费| 国产午夜精品论理片| 男女之事视频高清在线观看| 亚洲va日本ⅴa欧美va伊人久久| 性色av乱码一区二区三区2| 老司机深夜福利视频在线观看| 18禁美女被吸乳视频| 两性午夜刺激爽爽歪歪视频在线观看 | 很黄的视频免费| 一级a爱片免费观看的视频| 一本精品99久久精品77| 亚洲中文av在线| 亚洲va日本ⅴa欧美va伊人久久| a级毛片在线看网站| 亚洲一区二区三区色噜噜| 婷婷丁香在线五月| 亚洲国产欧美人成| 小说图片视频综合网站| xxxwww97欧美| 亚洲18禁久久av| 亚洲一区中文字幕在线| 久99久视频精品免费| 最新在线观看一区二区三区| 男女下面进入的视频免费午夜| 波多野结衣巨乳人妻| 在线观看午夜福利视频| 久久久国产欧美日韩av| 精品国内亚洲2022精品成人| 欧美zozozo另类| 此物有八面人人有两片| 久久国产精品人妻蜜桃| 国产高清有码在线观看视频 | a在线观看视频网站| 在线视频色国产色| 国产av麻豆久久久久久久| 午夜福利在线观看吧| 欧美日本亚洲视频在线播放| 国内精品久久久久精免费| 国产高清有码在线观看视频 | 制服诱惑二区| 久久精品夜夜夜夜夜久久蜜豆 | 一二三四社区在线视频社区8| 色老头精品视频在线观看| 热99re8久久精品国产| 一进一出抽搐gif免费好疼| 法律面前人人平等表现在哪些方面| 一个人免费在线观看的高清视频| 日本黄色视频三级网站网址| 中出人妻视频一区二区| 黄色视频不卡| av福利片在线观看| 成年免费大片在线观看| 日本精品一区二区三区蜜桃| 老司机在亚洲福利影院| 成年免费大片在线观看| 波多野结衣高清作品| 岛国在线观看网站| 99热6这里只有精品| 黄片小视频在线播放| 欧美精品啪啪一区二区三区| 非洲黑人性xxxx精品又粗又长| 亚洲男人天堂网一区| 老熟妇乱子伦视频在线观看| av片东京热男人的天堂| 久久久水蜜桃国产精品网| 久久国产精品影院| 精品欧美一区二区三区在线| 午夜福利免费观看在线| 动漫黄色视频在线观看| 校园春色视频在线观看|