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

    城市突發(fā)事件下的應(yīng)急物資配送路徑尋優(yōu)

    2016-08-25 06:16:42王付宇葉春明王濤安徽工業(yè)大學(xué)管理科學(xué)與工程學(xué)院安徽馬鞍山4303上海理工大學(xué)管理學(xué)院上海00093
    關(guān)鍵詞:救援車輛螢火蟲(chóng)物資

    王付宇,葉春明,王濤(.安徽工業(yè)大學(xué)管理科學(xué)與工程學(xué)院,安徽馬鞍山4303;.上海理工大學(xué)管理學(xué)院,上海00093)

    城市突發(fā)事件下的應(yīng)急物資配送路徑尋優(yōu)

    王付宇1,2,葉春明2,王濤1
    (1.安徽工業(yè)大學(xué)管理科學(xué)與工程學(xué)院,安徽馬鞍山243032;2.上海理工大學(xué)管理學(xué)院,上海200093)

    針對(duì)城市突發(fā)事件環(huán)境下的應(yīng)急物資配送車輛行駛路徑尋優(yōu)問(wèn)題,分析城市道路擁堵?tīng)顩r對(duì)于救援車輛行駛速度的影響。以車輛到達(dá)待救點(diǎn)的時(shí)間和行駛成本最少為目標(biāo),構(gòu)建雙層規(guī)劃路徑尋優(yōu)模型。為提高全局搜索和局部尋優(yōu)能力,采用一種改進(jìn)的兩階段式螢火蟲(chóng)算法。算例結(jié)果表明所提出的模型和算法可有效解決城市突發(fā)事件下的應(yīng)急物資配送問(wèn)題。

    應(yīng)急救援;路徑尋優(yōu);雙層理論;螢火蟲(chóng)算法

    近年來(lái),除了地震、洪澇等大規(guī)模自然災(zāi)害頻發(fā),城市內(nèi)部公共突發(fā)性事件也層出不窮,其帶來(lái)的生命威脅與財(cái)產(chǎn)損失不可忽視。此類城市范圍內(nèi)突發(fā)性事件的救援工作對(duì)于救援響應(yīng)時(shí)間和救援時(shí)間周期的緊迫性更加嚴(yán)格。然而城市區(qū)域內(nèi)復(fù)雜的路網(wǎng)和時(shí)變的路況對(duì)于救援車輛行駛速度有著明顯的影響。因此,為了確保城市發(fā)生各種突發(fā)事件后城市應(yīng)急救援中心能夠及時(shí)響應(yīng)、選擇最佳救援路線、實(shí)施高效救援,面向城市交通復(fù)雜路況下的應(yīng)急物資配送路徑尋優(yōu)已成為目前亟待解決的問(wèn)題[1-2]。

    目前,國(guó)內(nèi)外學(xué)者就應(yīng)急救援物資配送路徑問(wèn)題作了大量研究。Dantzig等[3]于1959年提出車輛路徑問(wèn)題,并指出該問(wèn)題本質(zhì)上屬于運(yùn)籌學(xué)與組合優(yōu)化領(lǐng)域。Balcik等[4]針對(duì)災(zāi)害后不確定信息下的救援和物資配送進(jìn)行研究,以總成本最小為目標(biāo),建立了從供應(yīng)點(diǎn)到資源需求點(diǎn)的多物品網(wǎng)絡(luò)流模型。Ceyhun[5]研究基于救援中心屬于覆蓋模型環(huán)境下的車輛路徑尋優(yōu)問(wèn)題,為使得每臺(tái)救援車輛能夠高效地服務(wù)更多的待救點(diǎn),建立了一種多目標(biāo)車輛路徑尋優(yōu)模型。Pretolani[6]分析某城市路段行駛時(shí)間的結(jié)果分布規(guī)律,總結(jié)出該城市路段行駛時(shí)間呈正態(tài)分布,該結(jié)果對(duì)應(yīng)急救援車輛行駛速度的判定有著很大啟發(fā)作用。馬祖軍等[7-8]對(duì)洪災(zāi)地震等自然災(zāi)害和城市公共突發(fā)事件下的應(yīng)急救援中心的選址及車輛路徑尋優(yōu)問(wèn)題進(jìn)行研究,提出應(yīng)急物資配送中心定位與配送車輛路徑安排的聯(lián)合決策問(wèn)題,對(duì)此建立了模糊多目標(biāo)定位及路徑尋優(yōu)模型。喬茹等[9]基于改進(jìn)蟻群算法研究移動(dòng)機(jī)器人全局路徑規(guī)劃問(wèn)題。姜金貴等[10]針對(duì)城市內(nèi)澇突發(fā)情況下,引用連通系數(shù)和暢通系數(shù)來(lái)計(jì)算車行速度,從而反映出道路積水阻礙交通等影響因素,建立相應(yīng)的救援路徑尋優(yōu)模型,以提高救援效率降低災(zāi)害損失。在求解方法上,Purezza等[11]運(yùn)用禁忌搜索算法求解車輛優(yōu)化調(diào)度的問(wèn)題;Torki等[12]提出用自組織的神經(jīng)網(wǎng)絡(luò)算法來(lái)求解車輛路線優(yōu)化問(wèn)題;宋曉宇等[13]用遺傳算法求解了灰色規(guī)劃模型,通過(guò)測(cè)試種群數(shù)目、交叉率和變異率3種控制參數(shù)值以提高算法性能。

    上述有關(guān)應(yīng)急救援車輛路徑尋優(yōu)研究中,僅僅采用1-0決策判斷方法處理道路是否能夠通行,而忽略了城市本身道路交通擁堵特征。其中采用的多目標(biāo)路徑尋優(yōu)模型在建立目標(biāo)函數(shù)時(shí),大多是采用權(quán)重分配的方法處理時(shí)間最短和成本最低多目標(biāo)之間的關(guān)系,而該方法所得出的結(jié)果容易忽略應(yīng)急救援過(guò)程中的及時(shí)性,即使時(shí)間的權(quán)重比較大,也會(huì)出現(xiàn)結(jié)果偏向于所耗時(shí)間和成本較小的可能。在求解方法上基本采取智能啟發(fā)式算法,對(duì)于群體智能算法的應(yīng)用相對(duì)不夠成熟?;诖?,本文針對(duì)城市車輛行駛速度特征進(jìn)行分析,提出道路擁擠條件下的雙層路徑尋優(yōu)模型,并構(gòu)建混沌螢火蟲(chóng)算法進(jìn)行求解。

    1 路徑尋優(yōu)模型

    1.1問(wèn)題描述

    傳統(tǒng)物流配送車輛路徑尋優(yōu)問(wèn)題,是在物流配送中心以及物資需求點(diǎn)坐標(biāo)位置、道路交通狀況、需求點(diǎn)物資需求量等因素已知的基礎(chǔ)上進(jìn)行的。而應(yīng)急救援車輛路徑規(guī)劃的過(guò)程中,救援需求點(diǎn)位置的隨機(jī)性、道路修復(fù)狀況、時(shí)變下的道路擁堵信息等不確定因素使得該情景下的車輛路徑尋優(yōu)問(wèn)題變得異常復(fù)雜。對(duì)應(yīng)急救援路徑選擇問(wèn)題的描述如下:選擇車輛行駛時(shí)間最短以及車輛運(yùn)輸成本最少作為目標(biāo)。系統(tǒng)中只有一個(gè)救援配送中心,且地理坐標(biāo)位置及所擁有的救援車型、數(shù)量、單車容量已知,救援中心派出車輛對(duì)多個(gè)待救點(diǎn)進(jìn)行服務(wù),完成物資配送后返回救援中心。

    車輛行駛過(guò)程中應(yīng)滿足以下幾個(gè)假設(shè)條件:

    1)道路節(jié)點(diǎn)之間有可行通路;

    2)每個(gè)待救點(diǎn)只能接受一輛車一次救援;

    3)每輛車在一次救援巡回過(guò)程中的載重不能超過(guò)其最大容量;

    4)車輛行駛速度小于等于其最大限制行駛速度。

    以城市交通路網(wǎng)作為建立應(yīng)急救援車輛最優(yōu)路徑模型基礎(chǔ),節(jié)點(diǎn)(待救點(diǎn))與節(jié)點(diǎn)之間的路線設(shè)定為一條弧段。應(yīng)急救援車輛路徑問(wèn)題可以用圖論形式表示出來(lái):記G(I,L)為賦權(quán)圖,I為頂點(diǎn)集,L為邊集,定義決策變量如下:

    其中:i,j表示待救點(diǎn)序號(hào),i,j=0,1,…,m;k代表派使的第k臺(tái)救援車輛,k=1,2,…,K;其中i,j=0,表示救援中心。式(1),(2)為決策變量,判斷救援車輛是否通過(guò)i和j節(jié)點(diǎn),以及節(jié)點(diǎn)i是否被第k輛救援車所提供服務(wù)。

    1.2數(shù)學(xué)模型

    雙層路徑尋優(yōu)模型不同于以往的多目標(biāo)路徑尋優(yōu)模型[14]。多目標(biāo)模型一般是以車輛行駛路徑最短、成本最低、時(shí)間最短等作為目標(biāo),并且諸多目標(biāo)之間以設(shè)定權(quán)值綜合考慮。而雙層路徑尋優(yōu)模型更加強(qiáng)調(diào)突發(fā)事件發(fā)生后對(duì)應(yīng)急物資配送的及時(shí)性,該模型以應(yīng)急救援物資配送車輛行駛總時(shí)間最短為主要優(yōu)化目標(biāo),在充分體現(xiàn)及時(shí)性的基礎(chǔ)上再以車輛行駛所耗總成本最低為輔助目標(biāo),從而使救援中心位置能夠及時(shí)服務(wù)于待救點(diǎn)傷員并且盡量降低救援車輛所耗成本。雙層路徑尋優(yōu)數(shù)學(xué)模型表達(dá)如下:

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

    其中:

    約束條件:

    其中:(x0,y0)表示救援中心位置;(xi,yi)表示待救點(diǎn)位置;di表示各待救點(diǎn)對(duì)該物資的需求量;D表示每輛救援車的最大容量;ck表示第k輛車,每行駛單位距離所消耗的成本;lij表示i與j之間距離;vij表示車輛行駛于i與j之間的速度;τij表示車輛行駛于i與j之間所耗時(shí)間;vˉ表示車輛的平均速度;vmax表示車輛的最大速度;qij表示i與j之間道路的暢通度;pij表示i與j之間道路的最大通行量。式(5)引入最大速度和平均速度,當(dāng)路徑暢通度與路容量以及平均速度的乘積小于最大速度時(shí),則實(shí)際速度為三者之積;反之,速度則控制在最大速度,不允許超速;式(3),(4)分別為應(yīng)急物資配送車輛行駛時(shí)間和所耗成本優(yōu)化函數(shù);式(7),(8)表示每個(gè)待救點(diǎn)均得到車輛救助,并且每個(gè)待救點(diǎn)只出現(xiàn)在一臺(tái)救援車輛的路徑上;式(9)保證每臺(tái)救援車輛歷經(jīng)的所有待救點(diǎn)的總需求量小于該車最大容量;式(10)~(12)保證每一臺(tái)車輛能夠形成可行回路。

    2 算法設(shè)計(jì)

    2.1混沌螢火蟲(chóng)算法思路

    在標(biāo)準(zhǔn)螢火蟲(chóng)算法的基礎(chǔ)上引入慣性權(quán)重,使算法的全局搜索和局部尋優(yōu)能力得到一定改善。但是,實(shí)驗(yàn)發(fā)現(xiàn),一旦達(dá)到局部極值點(diǎn)附近,算法收斂速率慢,且在極值點(diǎn)附近振蕩;對(duì)于一些多峰多極值點(diǎn)問(wèn)題,算法仍存在陷入局部最優(yōu)問(wèn)題?;煦邕\(yùn)動(dòng)具有遍歷性、隨機(jī)性等特點(diǎn),為此,本文在螢火蟲(chóng)算法的基礎(chǔ)上引入混沌思想,從而提高螢火蟲(chóng)種群的多樣性和尋優(yōu)的遍歷性,增加算法擺脫陷入局部極值點(diǎn)的能力。改進(jìn)后的混沌螢火蟲(chóng)算法是利用載波的方式,將混沌變量映射到螢火蟲(chóng)優(yōu)化算法的搜索空間,再利用混沌變量進(jìn)行搜索尋優(yōu)[16]。搜索的過(guò)程分為2個(gè)階段。

    1)粗略搜索

    為提高算法收斂速度,在迭代運(yùn)算初期,利用基于慣性權(quán)重的螢火蟲(chóng)算法進(jìn)行大范圍全局搜索,尋找整個(gè)搜索空間中的滿意極值點(diǎn)。

    2)細(xì)致搜索

    選取粗略搜索法所得適應(yīng)度高的10%至20%螢火蟲(chóng),結(jié)合混沌映射再進(jìn)行局部搜索,在粗略搜索的極值點(diǎn)附近尋找出全局下最優(yōu)解。

    2.2混沌螢火蟲(chóng)算法求解步驟

    改進(jìn)混沌螢火蟲(chóng)算法流程如下。

    1)系統(tǒng)初始化,生成螢火蟲(chóng)初始種群的規(guī)模、位置信息(Xa(t)),設(shè)置光強(qiáng)吸收系數(shù)(γ)、最大吸引度(β0)和步長(zhǎng)因子(α)等參數(shù)。

    2)計(jì)算群體中每個(gè)螢火蟲(chóng)的熒光亮度,

    其中I0表示螢火蟲(chóng)的最大熒光亮度,即自身r=0處。

    3)對(duì)所有螢火蟲(chóng)個(gè)體進(jìn)行鄰域搜索,此時(shí)的搜索半徑最大,并且計(jì)算出鄰域中所有候選螢火蟲(chóng)個(gè)體的吸引度,

    4)根據(jù)式(15),(16)更新螢火蟲(chóng)位置,最亮的螢火蟲(chóng)個(gè)體隨機(jī)移動(dòng),

    其中:t表示迭代次數(shù);tmax表示最大迭代次數(shù);ω(t)表示慣性權(quán)重;ωmax,ωmin表示最大、最小權(quán)重。

    5)在有限的迭代次數(shù)內(nèi),通過(guò)以上步驟可以得出一組螢火蟲(chóng)亮度數(shù)據(jù),從中選取螢火蟲(chóng)種群中適應(yīng)度高的20%個(gè)體,進(jìn)行混沌局部搜索。

    6)將螢火蟲(chóng)位置信息Xbt按下式映射為0到1之間的混沌參量C(Xbt),

    其中Xmax s和Xmin s分別為搜索空間中s維的上界和下界(假設(shè)螢火蟲(chóng)的位置空間為s維,當(dāng)s=1時(shí),Xmaxs=Xmax,Xmins=Xmin)。

    7)當(dāng)?shù)螖?shù)t=t+1時(shí),根據(jù)式(17)計(jì)算得到下步迭代的混沌參量,并將混沌參量)按式(18)映射轉(zhuǎn)化為Xt+1b,

    8)候選螢火蟲(chóng)個(gè)體適應(yīng)度值大于當(dāng)前螢火蟲(chóng)的適應(yīng)度值,則替換當(dāng)前螢火蟲(chóng)個(gè)體,實(shí)現(xiàn)螢火蟲(chóng)對(duì)較優(yōu)個(gè)體位置的移動(dòng)。

    9)按下式收縮搜索區(qū)域

    其中Xmaxlight.s表示當(dāng)前最亮螢火蟲(chóng)個(gè)體的第s維變量的值。

    10)若混沌搜索已到預(yù)先設(shè)定的精度或迭代次數(shù),則將新解作為算法的最終結(jié)果,否則令t=t+1并返回步驟6)。

    3 算例

    3.1算例背景

    某城市道路車輛行駛速度一般在20 km/h左右,但是在高峰期時(shí)段內(nèi),由于車流量過(guò)大導(dǎo)致車速受到嚴(yán)重影響。由該城市道路交通部門(mén)的車輛行駛監(jiān)控資料分析路網(wǎng)擁堵?tīng)顩r,道路暢通系數(shù)取值于0.5~1.0,道路交通容量取值范圍2~4,根據(jù)抽檢的車輛速度計(jì)算得該城市內(nèi)車輛行駛的最大速度為9.6 km/h、平均速度為3.6 km/h。另外,該城市大多數(shù)道路為東西和南北方向,為了方便理解在求解過(guò)程中直接取絕對(duì)值距離。

    取應(yīng)急救援中心周圍長(zhǎng)寬均為10 km區(qū)域?yàn)檠芯糠秶O(shè)待救點(diǎn)的坐標(biāo)為(xi,yi),根據(jù)該城市歷史突發(fā)事件相關(guān)數(shù)據(jù),救援中心(5,5)和待救點(diǎn)具體數(shù)據(jù)如表1所示。此外假設(shè),應(yīng)急物資配送車輛最大的容量為10 t,車輛行駛1 km所耗成本為5元。

    以該城市路網(wǎng)擁堵環(huán)境下的應(yīng)急物資配送問(wèn)題為研究對(duì)象,優(yōu)化目標(biāo)是在應(yīng)急物資配送車輛總行駛時(shí)間(T)可以接受的范圍內(nèi)([Tmin,Tmin(1%+2%)])使得車輛行駛成本最低。

    表1 需求點(diǎn)坐標(biāo)及需求量Tab.1 Demand Point CoordinateAnd Demand

    3.2算例求解

    1)參數(shù)設(shè)定

    在Matlab環(huán)境下建立相關(guān)數(shù)據(jù)庫(kù),其中包括救援中心及待救點(diǎn)坐標(biāo)(xi,yi)、各待救點(diǎn)的物資需求量(di)、車輛數(shù)目(K)、車輛k行駛單位距離消耗成本(ck)、單車最大容量(D)、車輛平均速度)、最大速度(vmax)、待救點(diǎn)間距離(lij)。另外,隨機(jī)產(chǎn)生21×21組道路暢通系數(shù)和通行容量(pij,qij)。該案例涉及的救援中心及待救點(diǎn)地理分布如圖1所示。螢火蟲(chóng)算法求解參數(shù)設(shè)定為螢火蟲(chóng)數(shù)m=10,光強(qiáng)吸收系數(shù)γ= 1.0,最大吸引度 β0=1.0,步長(zhǎng)因子α=0.2,迭代次數(shù)tmax=200,ωmax=1,ωmin=0.4。

    圖1 需求點(diǎn)分布Fig.1 Demand Point Distribution

    2)初始路徑方案

    根據(jù)最近鄰域法得出第一個(gè)車輛行駛路徑方案,共派出3輛應(yīng)急物資配送車輛,第1輛共載重為9.3 t,行駛22 km,路徑為0-2-3-4-5-20-7-6-0;第2輛共載重為9.0 t,行駛24.4 km,路徑為0-16-9-19-14-8-10-1-0;第3輛共載重為9.5 t,行駛38 km,路徑為0-15-12-13-11-17-18-0(簡(jiǎn)記①:0-2-3-4-5-20-7-6-0;0-16-9-19-14-8-10-1-0;0-15-12-13-11-17-18-0),如圖2所示。

    圖2 初始車輛行駛路徑方案Fig.2 Initial Vehicle Routing Scheme

    對(duì)第一個(gè)車輛行駛路徑方案進(jìn)行鄰域搜索,再產(chǎn)生9個(gè)初始方案,方法如下:

    (a)同一車輛,將某2個(gè)待救點(diǎn)位置進(jìn)行互換,共3種方案。

    每輛車歷經(jīng)的第2,4點(diǎn)交換;每輛車歷經(jīng)的第3,5點(diǎn)交換;每輛車歷經(jīng)的第4,7點(diǎn)交換。

    (b)兩個(gè)車輛路線,各選取一個(gè)待救點(diǎn)位置進(jìn)行互換,共3種方案。

    第1,2輛車所歷經(jīng)的第3點(diǎn)交換;第1,3輛車所歷經(jīng)的4點(diǎn)交換;第2,3輛車所歷經(jīng)的6點(diǎn)交換。

    (c)兩個(gè)車輛路線,取前者當(dāng)中某個(gè)待救點(diǎn),將其插入到后者車輛路徑當(dāng)中,共3種方案。

    將第1輛車經(jīng)過(guò)的第2待救點(diǎn),轉(zhuǎn)移到第2輛車經(jīng)過(guò)的第2待救點(diǎn);將第1輛車經(jīng)過(guò)的第2待救點(diǎn),轉(zhuǎn)移到第3輛車經(jīng)過(guò)的第2待救點(diǎn);將第2輛車經(jīng)過(guò)的第5待救點(diǎn),轉(zhuǎn)移到第3輛車經(jīng)過(guò)的第5待救點(diǎn)。

    3)全局尋優(yōu)

    將以上10個(gè)應(yīng)急物資配送車輛路徑方案及對(duì)應(yīng)的目標(biāo)函數(shù)值(由式(5),(6)計(jì)算得出)看作螢火蟲(chóng)初始種群,根據(jù)慣性權(quán)重螢火蟲(chóng)算法進(jìn)行求解得最優(yōu)的10個(gè)路徑方案及時(shí)間值。以①車輛行駛路徑方案為例,再次進(jìn)行鄰域搜索,比較當(dāng)前螢火蟲(chóng)亮度值和候選螢火蟲(chóng)亮度值,再根據(jù)位置更新公式實(shí)現(xiàn)螢火蟲(chóng)位置的優(yōu)化,經(jīng)過(guò)200次迭代,最終得出①方案對(duì)應(yīng)的最優(yōu)結(jié)果,記作①",同樣方法,得到另外9個(gè)最優(yōu)方案。全局搜索環(huán)境下車輛耗時(shí)的尋優(yōu)結(jié)果如圖3所示。

    圖3 全局搜索尋優(yōu)結(jié)果Fig.3 Global Search Optimization Results

    4)混沌序列

    從上述10個(gè)最優(yōu)方案當(dāng)中選取目標(biāo)值最小的前2個(gè)方案。路徑方案A(0-2-14-6-5-4-17-18-1-0;0-3-10-19-13-7-11-0;0-8-9-16-15-12-20-0),車輛總行駛時(shí)間為9.312 h,車輛總行駛距離為76.4 km,車輛總行駛所耗成本382元;路徑方案B(0-2-14-6-4-5-17-18-1-0;0-3-10-19-13-7-11-0;0-8-9-16-15-12-20-0),車輛總行駛時(shí)間為9.191 h,車輛總行駛距離為75.6 km,車輛總行駛所耗成本378元。

    將以上兩個(gè)結(jié)果對(duì)應(yīng)的螢火蟲(chóng)位置信息,按照式(17)映射到0到1之間的混沌參量,其中解為1維空間(j=1),xmin代表路徑方案A(0-2-14-6-5-4-17-18-1-0;0-3-10-19-13-7-11-0;0-8-9-16-15-12-20-0)作為搜索空間下限,xmax代表路徑方案B(0-2-14-6-4-5-17-18-1-0;0-3-10-19-13-7-11-0;0-89-16-15-12-20-0)作為搜索空間上限。

    5)局部尋優(yōu)

    在局部尋優(yōu)過(guò)程中先不考慮救援中心,把路徑方案看作20進(jìn)制的一個(gè)數(shù)字(且每個(gè)位置上的數(shù)字不重復(fù)),由此可根據(jù)式(17)來(lái)實(shí)現(xiàn)方案A和方案B的[0,1]映射。需要注意的是,在搜索的過(guò)程中,會(huì)產(chǎn)生不可行解,對(duì)于不可行解點(diǎn)進(jìn)行“+1”處理,直至找到該數(shù)值的最近可行解,然后通過(guò)式(18)進(jìn)行循環(huán)比較擇優(yōu),在限制的迭代次數(shù)內(nèi)得到最優(yōu)局部解點(diǎn)。局部搜索環(huán)境下車輛耗時(shí)的尋優(yōu)結(jié)果如圖4所示。

    圖4 局部搜索尋優(yōu)結(jié)果Fig.4 Local Search Optimization Curve

    3.3結(jié)果

    通過(guò)先后對(duì)解集進(jìn)行全局搜索及局部尋優(yōu),最后得到以車輛行駛時(shí)間最短為目標(biāo)的前三組最優(yōu)解,如表2所示。

    表2 路徑尋優(yōu)方案Tab.2 Route Optimization Scheme

    從上表可以看出車輛行駛路徑方案耗時(shí)最少為9.067 h,則救援時(shí)間可接受范圍是[9.067,9.248](由[Tmin,Tmin(1%+2%)]計(jì)算得出)。以上3個(gè)方案的車輛行駛時(shí)間值均在可接受范圍以內(nèi),則選取其中所耗成本最低的方案作為最優(yōu)方案,即方案2:共派出3輛救援車進(jìn)行應(yīng)急物資配送,總路線長(zhǎng)度為72.8 km,共耗時(shí)9.153 h,車輛行駛所耗成本為364元。3輛救援車的具體路徑分別為:0-14-6-2-4-5-17-18-1-0;0-3-10-19-13-11-7-0;0-8-9-16-15-12-20-0,如圖5所示。

    為進(jìn)一步證明提出的兩階段式混沌螢火蟲(chóng)算法在求解過(guò)程與求解結(jié)果中的優(yōu)越性,文中禁忌搜索算法對(duì)相同的案例進(jìn)行求解,2種算法求解結(jié)果對(duì)比見(jiàn)表3。由表3可看出,本文給出的混沌螢火蟲(chóng)算法尋優(yōu)能力強(qiáng)、車輛行駛耗時(shí)更短,費(fèi)用更少。

    圖5 車輛行駛路徑最優(yōu)方案Fig.5 Vehicle Travel Path Optimal Plan

    表3 兩種算法結(jié)果對(duì)比Tab.3 Comparison Of TwoAlgorithms

    4 結(jié) 語(yǔ)

    本文主要分析了城市上下班高峰期環(huán)境下道路車流量和道路通行容量對(duì)應(yīng)急物資配送車輛行駛速度的影響。在模型求解過(guò)程中,通過(guò)介入慣性權(quán)重和混沌算法有效地平衡了螢火蟲(chóng)算法的全局搜索能力和局部尋優(yōu)能力。算例結(jié)果表明,在求解突發(fā)事件情境下應(yīng)急救援車輛路徑優(yōu)化問(wèn)題時(shí),新模型及算法相對(duì)于就算法的求解結(jié)果更優(yōu),以耗時(shí)和費(fèi)用為目標(biāo)函數(shù)時(shí),總耗時(shí)更少,總費(fèi)用更少,所提出的路徑優(yōu)化模型及算法對(duì)于應(yīng)急救援工作的順利開(kāi)展有一定的指導(dǎo)性意義。

    不足之處在于選擇車輛路徑時(shí)假設(shè)所有節(jié)點(diǎn)之間均可以通行,并且直接取絕對(duì)值距離作為節(jié)點(diǎn)間距離,而城市實(shí)際路網(wǎng)錯(cuò)綜復(fù)雜、各節(jié)點(diǎn)間的道路方向多種多樣,因此模擬出城市實(shí)際路網(wǎng)作為路徑尋優(yōu)的依據(jù)、計(jì)算出精準(zhǔn)的節(jié)點(diǎn)間距離是應(yīng)急救援車輛路徑尋優(yōu)的下一步研究重點(diǎn)。

    [1]陳則輝,劉誠(chéng),呂品,等.不確定環(huán)境下應(yīng)急物資配送問(wèn)題研究[J].鐵道科學(xué)與工程學(xué)報(bào),2014(5):82-89.

    [2]李創(chuàng).國(guó)內(nèi)外應(yīng)急物流研究綜述[J].華東經(jīng)濟(jì)管理,2013,27(6):160-165.

    [3]DANTIZIG G,RAMSER J.The truck dispatching pmblem[J].Management Science,1959,6:80-91.

    [4]BALCIKB,BEAMONBM,SMILOWITZK.Lastmiledis-tributioninhumanitarianrelief[J].JournalofIntelligentTransportation System,2008,12(2):51-63.

    [5]ARAZ C,SELIM H,OZKARAHAN I.A fuzzy multi-objective covering-based vehicle location model for emergency services[J]. Computers&Operations Research,2007,34(3):705-726.

    [6]PRETOLANI D.A directed hypergraph model for random time-dependent shortest path[J].European Journal of Operational Research,2000,123(8):315-342.

    [7]馬祖軍,代穎,李雙琳.帶限制期的震后應(yīng)急物資配送模糊多目標(biāo)開(kāi)放式定位-路徑問(wèn)題[J].系統(tǒng)管理學(xué)報(bào),2014,23(5):658-667.

    [8]代穎,馬祖軍.應(yīng)急物流系統(tǒng)中的隨機(jī)定位-路徑問(wèn)題[J].系統(tǒng)管理學(xué)報(bào),2012,21(2):212-217.

    [9]喬茹,章小兵,趙光興.基于改進(jìn)蟻群算法的移動(dòng)機(jī)器人全局路徑規(guī)劃[J].安徽工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,26(1):77-80.

    [10]姜金貴,張鵬飛.基于改進(jìn)蟻群算法的城市內(nèi)澇救援路徑優(yōu)化[J].計(jì)算機(jī)應(yīng)用,2014,34(7):2103-2106.

    [11]PUREZA V M,F(xiàn)RANCA P M.Vehicle routing problems via tabu search metaheuritic[R].Centre De Recherche Sur Les Transports Publication,1991.

    [12]TORKIA,SOMHON S,ENKAWAT.Acompetitive neural network algorithm for solving vehicle routing problem[J].Computers &Industrial Engineering,1997,33(3):473-476.

    [13]宋曉宇,劉鋒,常春光.面向應(yīng)急物資調(diào)度的一種灰色規(guī)劃模型[J].計(jì)算機(jī)應(yīng)用研究,2010,27(4):1259-1262.

    [14]管小俊,王喜富,王翠華,等.基于競(jìng)爭(zhēng)的物流中心選址雙層規(guī)劃模型及算法研究[J].武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版),2009,33(5):956-959.

    [15]陳榮,李超群.基于ABC法和自適應(yīng)混合遺傳算法的倉(cāng)儲(chǔ)區(qū)域布局優(yōu)化策略[J].安徽工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,28(2):183-187.

    [16]劉長(zhǎng)平,葉春明.具有混沌搜索策略的螢火蟲(chóng)優(yōu)化算法[J].系統(tǒng)管理學(xué)報(bào),2013,22(4):538-543.

    [17]湯金春,王培珍.基于遺傳算法的鋼錠配切優(yōu)化[J].安徽工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,29(1):67-69.

    [18]郜振華,梅莉,祝遠(yuǎn)鑒.復(fù)合策略慣性權(quán)重的粒子群優(yōu)化算法[J].計(jì)算機(jī)應(yīng)用,2012,32(8):2216-2218.

    責(zé)任編輯:丁吉海

    Route Optimization of Emergency Material Distribution under Urban Emergency Situations

    WANG Fuyu1,2,YE Chunming2,WANG Tao1
    (1.School of management science and Engineering,Anhui University of Technology,Ma'anshan 243032,China;2.School of management,University of Shanghai for Science and Technology,Shanghai 200093,China)

    Aiming at the take of the vehicle routing optimization of emergency material distribution in the urban emergency environment,the influence of urban road congestion on the speed of the vehicle was analyzed.Take the least time and cost of the vehicle as object,a routing optimization model of bi-level programming is built.To improve the global search and local search capability,an improved two stage firefly algorithm was employed to solve the problem.The numerical example shows that the proposed model and algorithm can effectively solve the problem of emergency material distribution in urban emergency.

    emergency rescue;route optimization;two-level theory;firefly algorithm

    TP391

    Adoi:10.3969/j.issn.1671-7872.2016.02.016

    1671-7872(2016)02-0177-08

    2015-11-13

    國(guó)家自然科學(xué)基金項(xiàng)目(71271138);教育部人文社會(huì)科學(xué)青年基金項(xiàng)目(14YJC630119);住建部軟科學(xué)研究項(xiàng)目(2015-R2-057);安徽省高校人文社科研究重大項(xiàng)目(SK2014ZD016)

    王付宇(1977-),男,河南泌陽(yáng)人,副教授,主要研究方向?yàn)橹悄軆?yōu)化、工業(yè)工程。

    猜你喜歡
    救援車輛螢火蟲(chóng)物資
    被偷的救援物資
    螢火蟲(chóng)
    電力企業(yè)物資管理模式探討
    國(guó)務(wù)院辦公廳印發(fā)《關(guān)于國(guó)家綜合性消防救援車輛懸掛應(yīng)急救 援專用號(hào)牌有關(guān)事項(xiàng)的通知》
    螢火蟲(chóng)
    應(yīng)急救援車輛產(chǎn)品概況
    救援物資
    抱抱就不哭了
    夏天的螢火蟲(chóng)
    PKPM物資管理系統(tǒng)應(yīng)用實(shí)踐
    在线看a的网站| 欧美日韩亚洲高清精品| 久久久久精品国产欧美久久久 | 久久久精品区二区三区| av国产精品久久久久影院| 亚洲第一欧美日韩一区二区三区 | 日本91视频免费播放| 国产精品免费视频内射| 91精品伊人久久大香线蕉| 成年女人毛片免费观看观看9 | 青春草亚洲视频在线观看| 亚洲精品乱久久久久久| 欧美日韩福利视频一区二区| www.熟女人妻精品国产| 黄色视频在线播放观看不卡| 日韩制服丝袜自拍偷拍| 久久中文看片网| 国产福利在线免费观看视频| 麻豆国产av国片精品| 免费人妻精品一区二区三区视频| 免费观看av网站的网址| 精品人妻一区二区三区麻豆| 自线自在国产av| 亚洲中文字幕日韩| 一边摸一边抽搐一进一出视频| 丰满少妇做爰视频| 人妻一区二区av| 精品亚洲乱码少妇综合久久| 韩国高清视频一区二区三区| 欧美精品啪啪一区二区三区 | 两性午夜刺激爽爽歪歪视频在线观看 | 美女福利国产在线| 伦理电影免费视频| 9色porny在线观看| 欧美性长视频在线观看| 亚洲精品成人av观看孕妇| 嫩草影视91久久| 一区二区日韩欧美中文字幕| 亚洲人成77777在线视频| 欧美精品一区二区大全| 国产精品欧美亚洲77777| 国产1区2区3区精品| a 毛片基地| 亚洲 国产 在线| 久久久精品免费免费高清| 在线天堂中文资源库| 久久人妻福利社区极品人妻图片| 亚洲美女黄色视频免费看| 精品久久久久久电影网| 91精品国产国语对白视频| 欧美激情久久久久久爽电影 | 纯流量卡能插随身wifi吗| 久久久国产欧美日韩av| 精品人妻熟女毛片av久久网站| 欧美日韩成人在线一区二区| 免费久久久久久久精品成人欧美视频| 久久精品国产综合久久久| 日韩一区二区三区影片| h视频一区二区三区| 51午夜福利影视在线观看| 欧美人与性动交α欧美精品济南到| 熟女少妇亚洲综合色aaa.| 欧美人与性动交α欧美软件| av线在线观看网站| kizo精华| 美女大奶头黄色视频| 欧美一级毛片孕妇| 日日摸夜夜添夜夜添小说| 久久国产精品影院| 午夜福利一区二区在线看| 亚洲九九香蕉| 日本wwww免费看| 亚洲中文日韩欧美视频| 国产日韩欧美视频二区| 亚洲第一av免费看| 欧美黑人欧美精品刺激| 国产精品 国内视频| 国产成人免费观看mmmm| 我的亚洲天堂| 秋霞在线观看毛片| 欧美在线黄色| 日韩视频一区二区在线观看| 美女主播在线视频| 精品少妇内射三级| 老汉色∧v一级毛片| av视频免费观看在线观看| 精品国产一区二区三区久久久樱花| 精品国产乱码久久久久久小说| 国产一区二区 视频在线| av电影中文网址| 下体分泌物呈黄色| av电影中文网址| 色综合欧美亚洲国产小说| 99久久精品国产亚洲精品| 菩萨蛮人人尽说江南好唐韦庄| av福利片在线| 亚洲人成电影免费在线| 国产精品秋霞免费鲁丝片| 国产极品粉嫩免费观看在线| 亚洲激情五月婷婷啪啪| 亚洲精品中文字幕一二三四区 | 久久毛片免费看一区二区三区| 男女午夜视频在线观看| 精品福利观看| 超碰97精品在线观看| 他把我摸到了高潮在线观看 | 国产精品成人在线| 男女高潮啪啪啪动态图| a级毛片黄视频| 波多野结衣av一区二区av| 黑丝袜美女国产一区| 精品国产乱码久久久久久小说| 精品第一国产精品| 久久久久久免费高清国产稀缺| 久久精品国产亚洲av香蕉五月 | 久久久国产欧美日韩av| 国产福利在线免费观看视频| av又黄又爽大尺度在线免费看| 亚洲欧美一区二区三区久久| 777久久人妻少妇嫩草av网站| 亚洲专区国产一区二区| 另类精品久久| 老熟女久久久| www.自偷自拍.com| 国产一区有黄有色的免费视频| 91成人精品电影| 搡老岳熟女国产| 成人影院久久| 久久久久精品国产欧美久久久 | www.熟女人妻精品国产| 高清av免费在线| 欧美精品av麻豆av| 欧美性长视频在线观看| 成年动漫av网址| 免费人妻精品一区二区三区视频| 制服人妻中文乱码| 秋霞在线观看毛片| 国产91精品成人一区二区三区 | 精品久久蜜臀av无| 天天躁日日躁夜夜躁夜夜| 人人妻,人人澡人人爽秒播| 国产在线免费精品| 男人操女人黄网站| 高潮久久久久久久久久久不卡| 黑人巨大精品欧美一区二区mp4| av欧美777| 欧美精品一区二区大全| 欧美少妇被猛烈插入视频| 久久人人爽人人片av| 久久精品aⅴ一区二区三区四区| 男女免费视频国产| 精品一区二区三区av网在线观看 | 老司机在亚洲福利影院| a 毛片基地| 深夜精品福利| 男女免费视频国产| 国产男女内射视频| 国产亚洲一区二区精品| 久久亚洲国产成人精品v| 欧美精品人与动牲交sv欧美| 黄片播放在线免费| 婷婷成人精品国产| 在线精品无人区一区二区三| 一进一出抽搐动态| 丁香六月欧美| 日本五十路高清| 精品福利永久在线观看| svipshipincom国产片| 男人爽女人下面视频在线观看| 国产又色又爽无遮挡免| 最新的欧美精品一区二区| 大片电影免费在线观看免费| 黄网站色视频无遮挡免费观看| 蜜桃在线观看..| 久久久久国内视频| 午夜福利视频精品| 99国产极品粉嫩在线观看| 高清视频免费观看一区二区| 日韩有码中文字幕| 精品卡一卡二卡四卡免费| 国产一区有黄有色的免费视频| 久久精品成人免费网站| 人人澡人人妻人| 久久人人97超碰香蕉20202| www.av在线官网国产| 亚洲国产av新网站| 欧美日韩av久久| 久久天堂一区二区三区四区| 欧美97在线视频| 欧美黑人欧美精品刺激| 日韩制服骚丝袜av| 欧美精品啪啪一区二区三区 | 国产精品成人在线| 高潮久久久久久久久久久不卡| 韩国高清视频一区二区三区| 成人免费观看视频高清| 国产97色在线日韩免费| 亚洲专区字幕在线| 美女高潮喷水抽搐中文字幕| 免费观看a级毛片全部| 人人妻人人添人人爽欧美一区卜| 每晚都被弄得嗷嗷叫到高潮| 性高湖久久久久久久久免费观看| 免费不卡黄色视频| 黑人操中国人逼视频| 曰老女人黄片| 水蜜桃什么品种好| 国产亚洲av高清不卡| 欧美av亚洲av综合av国产av| 亚洲,欧美精品.| 视频区欧美日本亚洲| 成人免费观看视频高清| 亚洲精品国产精品久久久不卡| 纵有疾风起免费观看全集完整版| 国产精品偷伦视频观看了| 在线观看免费午夜福利视频| 精品熟女少妇八av免费久了| 午夜激情久久久久久久| 男女国产视频网站| 亚洲成人手机| 老熟妇乱子伦视频在线观看 | 精品国产一区二区三区四区第35| 国产精品一区二区在线观看99| 别揉我奶头~嗯~啊~动态视频 | 在线十欧美十亚洲十日本专区| 丝袜喷水一区| 99香蕉大伊视频| 午夜91福利影院| 亚洲国产看品久久| 精品熟女少妇八av免费久了| 亚洲va日本ⅴa欧美va伊人久久 | 久久香蕉激情| 亚洲国产欧美网| 国产精品麻豆人妻色哟哟久久| 国产精品偷伦视频观看了| 亚洲中文日韩欧美视频| 成年动漫av网址| 桃红色精品国产亚洲av| 欧美人与性动交α欧美软件| 久久av网站| 午夜老司机福利片| 精品国产超薄肉色丝袜足j| 国产一区二区激情短视频 | 在线观看免费视频网站a站| 这个男人来自地球电影免费观看| 亚洲中文日韩欧美视频| 黄片小视频在线播放| 亚洲欧美成人综合另类久久久| 国产亚洲精品一区二区www | 日韩 亚洲 欧美在线| 久久亚洲国产成人精品v| av网站免费在线观看视频| 精品欧美一区二区三区在线| 又紧又爽又黄一区二区| 女性生殖器流出的白浆| 大片电影免费在线观看免费| 亚洲精品中文字幕在线视频| 国产精品九九99| 欧美中文综合在线视频| 三上悠亚av全集在线观看| 日本精品一区二区三区蜜桃| 精品欧美一区二区三区在线| 黄网站色视频无遮挡免费观看| 十八禁网站免费在线| 91麻豆av在线| 免费日韩欧美在线观看| 欧美午夜高清在线| 亚洲精品日韩在线中文字幕| 91精品国产国语对白视频| 欧美精品高潮呻吟av久久| 国产一区二区三区av在线| 少妇猛男粗大的猛烈进出视频| 日本91视频免费播放| 亚洲熟女精品中文字幕| 一区二区三区乱码不卡18| 久久狼人影院| 成人黄色视频免费在线看| 久久 成人 亚洲| 99久久人妻综合| 亚洲av日韩在线播放| 精品少妇久久久久久888优播| 国产深夜福利视频在线观看| cao死你这个sao货| 国产99久久九九免费精品| 亚洲中文字幕日韩| 91成年电影在线观看| www.精华液| av免费在线观看网站| 中文字幕人妻丝袜一区二区| 亚洲精品乱久久久久久| 超碰97精品在线观看| 亚洲人成77777在线视频| 久久久久精品国产欧美久久久 | 午夜福利免费观看在线| av又黄又爽大尺度在线免费看| 亚洲中文字幕日韩| 俄罗斯特黄特色一大片| 久久久国产一区二区| 97人妻天天添夜夜摸| 国产精品二区激情视频| 黄片大片在线免费观看| 性少妇av在线| 亚洲精品在线美女| 一级片免费观看大全| 亚洲欧美日韩高清在线视频 | 欧美老熟妇乱子伦牲交| 国产国语露脸激情在线看| 久久精品亚洲熟妇少妇任你| 国产成人精品在线电影| 黑人猛操日本美女一级片| 999精品在线视频| 日本a在线网址| 老司机福利观看| 一区二区三区激情视频| 老汉色av国产亚洲站长工具| 人妻久久中文字幕网| 狂野欧美激情性bbbbbb| 国产成人精品在线电影| 黄片播放在线免费| 国产精品影院久久| 香蕉丝袜av| 啦啦啦在线免费观看视频4| 国产日韩欧美在线精品| 麻豆国产av国片精品| 亚洲午夜精品一区,二区,三区| 高清视频免费观看一区二区| 中文字幕色久视频| 久久久久久久国产电影| 欧美在线一区亚洲| 亚洲精品av麻豆狂野| 91麻豆精品激情在线观看国产 | 9191精品国产免费久久| 国产真人三级小视频在线观看| 在线看a的网站| 国产一区二区在线观看av| 99久久精品国产亚洲精品| 视频区图区小说| 精品欧美一区二区三区在线| netflix在线观看网站| 老熟妇仑乱视频hdxx| 欧美日韩中文字幕国产精品一区二区三区 | 免费一级毛片在线播放高清视频 | 国产亚洲精品一区二区www | 亚洲中文日韩欧美视频| 日韩精品免费视频一区二区三区| a在线观看视频网站| 色老头精品视频在线观看| 国产精品av久久久久免费| 91精品国产国语对白视频| 国产一区二区三区综合在线观看| 淫妇啪啪啪对白视频 | 免费观看人在逋| 看免费av毛片| 国产伦理片在线播放av一区| 黑人欧美特级aaaaaa片| 美女午夜性视频免费| 母亲3免费完整高清在线观看| 欧美 亚洲 国产 日韩一| 19禁男女啪啪无遮挡网站| 50天的宝宝边吃奶边哭怎么回事| 久久毛片免费看一区二区三区| 操美女的视频在线观看| 亚洲少妇的诱惑av| 黄色片一级片一级黄色片| 久久久久久久大尺度免费视频| 免费看十八禁软件| 在线 av 中文字幕| 国内毛片毛片毛片毛片毛片| 啦啦啦视频在线资源免费观看| 91精品三级在线观看| 婷婷丁香在线五月| 波多野结衣av一区二区av| 首页视频小说图片口味搜索| 国产在线一区二区三区精| 老鸭窝网址在线观看| 国产免费一区二区三区四区乱码| 亚洲精品美女久久av网站| 青青草视频在线视频观看| 国产淫语在线视频| 国产成人一区二区三区免费视频网站| 国产深夜福利视频在线观看| bbb黄色大片| 亚洲 国产 在线| 丝袜人妻中文字幕| 国产视频一区二区在线看| 99国产精品一区二区蜜桃av | 99热国产这里只有精品6| 亚洲国产成人一精品久久久| 这个男人来自地球电影免费观看| 精品福利观看| 国产免费福利视频在线观看| 老司机亚洲免费影院| 色94色欧美一区二区| 50天的宝宝边吃奶边哭怎么回事| 国产在线免费精品| 国精品久久久久久国模美| 免费在线观看完整版高清| 99久久精品国产亚洲精品| 91成人精品电影| 精品国产一区二区久久| 亚洲精品国产av蜜桃| 人人妻人人澡人人看| 丝袜在线中文字幕| 欧美成狂野欧美在线观看| 黑人巨大精品欧美一区二区mp4| 香蕉丝袜av| 免费日韩欧美在线观看| 精品免费久久久久久久清纯 | 在线av久久热| 不卡一级毛片| 成人手机av| 曰老女人黄片| 黄片大片在线免费观看| 一二三四在线观看免费中文在| 这个男人来自地球电影免费观看| 亚洲中文字幕日韩| 国产黄色免费在线视频| 女人被躁到高潮嗷嗷叫费观| 亚洲色图 男人天堂 中文字幕| 中文字幕制服av| 国产一区二区三区av在线| 色老头精品视频在线观看| 亚洲九九香蕉| 午夜91福利影院| 黄频高清免费视频| 在线十欧美十亚洲十日本专区| 电影成人av| 欧美性长视频在线观看| 一级片免费观看大全| www日本在线高清视频| 国产av精品麻豆| 亚洲自偷自拍图片 自拍| 精品国产一区二区三区久久久樱花| 日日摸夜夜添夜夜添小说| 蜜桃在线观看..| 精品少妇久久久久久888优播| 亚洲成人国产一区在线观看| 搡老乐熟女国产| 十八禁网站免费在线| 女性被躁到高潮视频| h视频一区二区三区| 十分钟在线观看高清视频www| 精品一区二区三卡| 美女扒开内裤让男人捅视频| 国产av一区二区精品久久| 狂野欧美激情性xxxx| 无遮挡黄片免费观看| www.精华液| 一级a爱视频在线免费观看| 国产在线免费精品| 亚洲成人手机| 亚洲一码二码三码区别大吗| 亚洲国产欧美网| 精品国产一区二区三区四区第35| 人人妻人人澡人人看| 精品高清国产在线一区| 黄色视频,在线免费观看| 久久国产亚洲av麻豆专区| 亚洲精品粉嫩美女一区| 久久 成人 亚洲| 99久久99久久久精品蜜桃| 在线观看一区二区三区激情| 大码成人一级视频| 国产xxxxx性猛交| 91大片在线观看| 国产一卡二卡三卡精品| 国产av国产精品国产| 在线天堂中文资源库| 法律面前人人平等表现在哪些方面 | 国产成人影院久久av| 欧美成人午夜精品| 别揉我奶头~嗯~啊~动态视频 | 亚洲一区二区三区欧美精品| 亚洲欧美日韩另类电影网站| 成年美女黄网站色视频大全免费| av免费在线观看网站| 亚洲欧美精品自产自拍| 男女国产视频网站| 日本五十路高清| 天天添夜夜摸| 国产精品 欧美亚洲| 亚洲专区字幕在线| 国产精品久久久人人做人人爽| 性色av乱码一区二区三区2| 亚洲欧洲精品一区二区精品久久久| 三上悠亚av全集在线观看| 天天添夜夜摸| 久久久久久久久免费视频了| www.自偷自拍.com| 亚洲熟女精品中文字幕| 久久久久网色| 国产三级黄色录像| 日本猛色少妇xxxxx猛交久久| 国产无遮挡羞羞视频在线观看| 一级,二级,三级黄色视频| 中文字幕高清在线视频| 十八禁网站免费在线| 久久综合国产亚洲精品| 国产av精品麻豆| 国产在视频线精品| 久久久久国内视频| 婷婷色av中文字幕| 69精品国产乱码久久久| 无遮挡黄片免费观看| 精品久久蜜臀av无| 午夜福利,免费看| av在线播放精品| 天堂中文最新版在线下载| 老司机亚洲免费影院| 一级片免费观看大全| 纯流量卡能插随身wifi吗| 精品一区二区三区四区五区乱码| 久久精品aⅴ一区二区三区四区| 91大片在线观看| a级片在线免费高清观看视频| 久久久久久久精品精品| 久久女婷五月综合色啪小说| 手机成人av网站| 超色免费av| 一区二区三区激情视频| 十分钟在线观看高清视频www| 美女午夜性视频免费| 黄色怎么调成土黄色| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲av美国av| 久久 成人 亚洲| 啪啪无遮挡十八禁网站| 久久人人爽人人片av| 日本vs欧美在线观看视频| 亚洲国产成人一精品久久久| 国产熟女午夜一区二区三区| 搡老岳熟女国产| 黄色片一级片一级黄色片| 丰满饥渴人妻一区二区三| 久久精品亚洲av国产电影网| 黄色 视频免费看| 最黄视频免费看| 涩涩av久久男人的天堂| 国产一区二区 视频在线| 岛国毛片在线播放| 欧美精品亚洲一区二区| 欧美精品一区二区大全| 下体分泌物呈黄色| 热99re8久久精品国产| 好男人电影高清在线观看| 国产男女内射视频| 秋霞在线观看毛片| 午夜老司机福利片| 蜜桃在线观看..| 一区二区三区四区激情视频| 99国产精品一区二区蜜桃av | 日本黄色日本黄色录像| 亚洲欧洲精品一区二区精品久久久| 丝袜喷水一区| 国产亚洲欧美精品永久| 国产三级黄色录像| 一本久久精品| tocl精华| 免费在线观看影片大全网站| 欧美黄色淫秽网站| av视频免费观看在线观看| 成人影院久久| 免费女性裸体啪啪无遮挡网站| 狂野欧美激情性bbbbbb| av片东京热男人的天堂| av天堂在线播放| 欧美久久黑人一区二区| 19禁男女啪啪无遮挡网站| 老司机在亚洲福利影院| 后天国语完整版免费观看| 国产野战对白在线观看| 夫妻午夜视频| 亚洲精品国产av蜜桃| 丰满人妻熟妇乱又伦精品不卡| 十八禁网站网址无遮挡| 9热在线视频观看99| 十八禁人妻一区二区| 亚洲精品一区蜜桃| 欧美一级毛片孕妇| 国产三级黄色录像| 精品一区二区三区av网在线观看 | 日韩制服骚丝袜av| 午夜福利视频在线观看免费| 中文字幕人妻熟女乱码| 亚洲欧美激情在线| 最近最新免费中文字幕在线| 97人妻天天添夜夜摸| 亚洲黑人精品在线| 欧美日韩av久久| 亚洲欧美一区二区三区黑人| 国产日韩一区二区三区精品不卡| 日韩电影二区| 日韩,欧美,国产一区二区三区| 永久免费av网站大全| 久久影院123| av在线app专区| 国产成人免费观看mmmm| 自线自在国产av| 欧美日韩视频精品一区| 秋霞在线观看毛片| 日韩大码丰满熟妇| 咕卡用的链子| 亚洲成人免费电影在线观看| 久久国产精品男人的天堂亚洲| 最新在线观看一区二区三区| 婷婷色av中文字幕| 一级,二级,三级黄色视频| 在线十欧美十亚洲十日本专区| 精品国产一区二区三区久久久樱花| 黄片播放在线免费| 女性生殖器流出的白浆| 色视频在线一区二区三区| 国产一区二区三区综合在线观看| 制服人妻中文乱码| 精品少妇黑人巨大在线播放|