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

    考慮工作量平衡的餐飲垃圾多行程收運路線優(yōu)化

    2023-12-28 02:54:20張燕李子鑫劉進平
    關(guān)鍵詞:收運時段算子

    張燕,李子鑫,劉進平

    (大連海事大學,交通運輸工程學院,遼寧大連 116026)

    0 引言

    作為具有悠久飲食文化的國家,近年來隨著經(jīng)濟發(fā)展和城市化進程的加快,我國城市的各大餐館和企事業(yè)單位食堂產(chǎn)生的餐飲垃圾也隨之激增。為改善環(huán)境,更有效地利用資源,很多城市開始試點餐飲垃圾分類收運和處理。然而,高收運成本成為推進餐飲垃圾分類管理的一大障礙,如何以更低的成本完成餐飲垃圾的分類收運成為實施垃圾分類管理的重中之重[1]。同時,隨著智慧城市成為未來城市的發(fā)展理念,很多行業(yè)和企業(yè)嘗試采用優(yōu)化算法降本增效,但大多數(shù)算法只關(guān)注成本,沒有考慮工作人員的人本需求[2]。本文針對餐飲垃圾的收運特點提出相應的收運模型及算法,致力于在降低餐飲垃圾分類收運成本的同時,提升駕駛員的工作滿意度,對于提升城市廢物管理水平,推動城市的可持續(xù)發(fā)展具有重要意義。

    在以往的相關(guān)研究中,不少學者針對居民生活垃圾的收運路線優(yōu)化問題進行了研究。例如,LEI等[3]考慮生活垃圾的周期性收運路線優(yōu)化和車輛調(diào)度問題,建立多周期的車輛路徑問題模型,并采用兩階段算法求解問題;ASEFI 等[4]考慮多種廢物類型的收運問題,建立車隊規(guī)模和混合車輛路徑優(yōu)化問題模型,并對小規(guī)模問題進行求解;ALIAHMADI等[5]考慮運營成本和收運時間兩個目標,建立帶時間窗的多行程車輛路徑問題模型,并采用多目標優(yōu)化算法求解問題;閆芳等[6]結(jié)合智能垃圾桶提出動態(tài)生活垃圾收運優(yōu)化策略,通過粒子群算法獲得初始收運方案,然后,采用重置優(yōu)化算法實時調(diào)整收運路徑。同時,也有學者結(jié)合新的技術(shù)方法,例如,李明月等[7]考慮時間成本及環(huán)境成本,提出一種基于GIS技術(shù)與改進蟻群算法的垃圾收運路徑優(yōu)化方法;GLAESER 等[8]針對地下廢物收運系統(tǒng),建立以收運成本最小化為目標的垃圾箱分配和車輛路徑問題的數(shù)學模型,并設計變鄰域啟發(fā)式算法求解問題。

    與城市生活垃圾的收運不同,餐飲垃圾的收運路線優(yōu)化問題需要考慮以下因素:首先,由于各個餐飲企業(yè)及企事業(yè)單位食堂產(chǎn)生餐飲垃圾量和高峰時段不同,因此,需要針對不同的企業(yè)類型設置不同的收運頻率、收運時段和服務時間窗,以便均衡各個時段的垃圾收運量;第二,餐飲垃圾含水量大,難以壓縮,餐飲垃圾收運車在每個收運時段內(nèi)需要多次進出餐飲垃圾處理廠進行卸貨;第三,城市餐飲企業(yè)集中于商業(yè)區(qū)或者小吃街附近,因此,在設計優(yōu)化路線時需考慮實際道路的擁堵因素以及限行要求。

    因此,餐飲垃圾的收運路線優(yōu)化是一個考慮時間窗約束的多行程車輛路徑問題(Multi-Trip Vehicle Routing Problem with Time Windows,MTVRPTW),學者們針對該問題提出了不同的模型和算法。例如,F(xiàn)RANCOIS 等[9]建立以總行駛距離為目標函數(shù)的多行程車輛路徑數(shù)學模型,并設計帶有多行程算子的自適應大鄰域搜索算法;彭勇等[10]考慮帶時間窗的多行程交換箱甩掛路徑問題,并設計混合遺傳算法,以充分利用司機日常工作時間;PAN 等[11]以最小總運輸距離為目標,研究城市物流配送中時變路網(wǎng)下的多行程車輛路徑問題,并采用混合元啟發(fā)式算法求解問題。但是,現(xiàn)有研究大多考慮的都是單一車場,且大多是以總行駛距離或時間為目標函數(shù),缺少針對餐飲垃圾的收運特點所提出的相應模型和算法。

    本文首先建立具有時間窗約束的多行程餐飲垃圾收運問題的混合整數(shù)規(guī)劃模型,并根據(jù)問題特點設計混合自適應大鄰域搜索算法。在模型和算法中,考慮駕駛員的各項工作需求,例如,保障駕駛員的午休以及工作強度不超標,保證駕駛員之間的工作平衡;同時,還考慮城市實際路網(wǎng)的擁堵情況,得到既節(jié)省成本又降低碳排放,還能提高駕駛員工作滿意度的最優(yōu)路線。

    1 問題描述與模型建立

    1.1 問題描述

    本文考慮城市中的餐飲垃圾收運路線優(yōu)化問題,其中,餐飲垃圾的收運來源包括各種不同類型的餐館和企事業(yè)單位食堂。首先,根據(jù)各類型餐館和單位食堂的營業(yè)特點,設定兩種收運頻率:每天兩次或每天一次。對于夜市及早餐為主的餐飲企業(yè),每天上午(時段1)收運一次,對于不供應晚餐的企事業(yè)單位食堂等,每天下午(時段2)收運一次;對于全天營業(yè)的餐館及食堂,每天上午和下午各收運一次。

    餐飲垃圾的收運網(wǎng)絡由餐飲垃圾的收運點、車場和餐飲垃圾處理廠組成。每輛餐飲垃圾收運車在每個收運時段內(nèi)完成一個混合的多行程任務,如圖1 所示。每輛車收運任務的第一段行程為開放式行程,車輛從車場空車出發(fā),沿途訪問收運點之后,到達處理廠卸貨,之后的收運任務是以處理廠為中心的多個閉合式行程,最后一段是完成收運任務后空載返回車場的行程。

    圖1 混合多行程的餐飲垃圾收運路線示例Fig.1 Illustration of hybrid multi-trip food waste collection route

    駕駛員在上午時間段的收運任務完成后,返回車場午休后,再進行下午的收運。駕駛員每天的工作任務總時長不超過最長工作時間,并且各個駕駛員之間的工作量盡可能平衡。同時,本文考慮車輛在實際道路行駛過程中的擁堵因素,將車輛在任意兩點間的行駛時間計為正常行駛時間加上擁堵時間。假設餐飲垃圾的處理設施容量足夠大,每個收運點的餐飲垃圾不允許分批收運。該問題的目標是在保證駕駛員工作量不超標且平衡的情況下為每個時段設計合理的收運路線,從而最小化收運過程中的車輛固定成本和燃油成本。

    1.2 模型參數(shù)與變量

    (1)相關(guān)集合

    Na={1,2,…,n}——所有餐飲垃圾收運點的集合,其中,n為收運點數(shù)目;

    Ne={0,1,2,…,n,n+1}——包含車場、收運點和處理廠的擴展集合,其中,0為車場,n+1 為餐飲垃圾處理廠;

    Np——在時段p收運的所有收運點集合;

    V={1,2,…,Vmax}——收運車輛集合,其中,Vmax為收運車輛的總數(shù)量;

    P={1,2,…,T}——餐飲垃圾收運時段集合。

    (2)參數(shù)定義

    對于收運點(i,j)∈Ne,車輛k∈V,時段p∈P,定義以下參數(shù)。

    tˉij——收運車在收運點(i,j)間的正常行駛時間(min);

    Qmax——餐飲垃圾收運車的最大凈載量(kg);

    cV——收運車的單位固定運營成本(元·d-1·輛-1),包含車輛折舊和駕駛員薪酬;

    cS——收運車在擁堵和裝卸服務時的怠速燃油成本(元·min-1);

    cD——收運車正常行駛時的單位燃油成本(元·min-1);

    H——收運車在每個時段內(nèi)可以進行多行程收運的最大次數(shù);

    M——極大的正數(shù)。

    (3)決策變量

    1.3 數(shù)學模型

    根據(jù)以上參數(shù)和變量定義,本文所考慮的帶時間窗約束的混合多行程的餐飲垃圾收運問題(HMTVRPTW)的數(shù)學模型(模型M)為

    目標函數(shù)式(1)表示最小化所有時段的總收運成本,包括車輛的固定運營成本、車輛在交通擁堵以及裝卸作業(yè)時的怠速燃油成本和正常行駛時的燃油成本。

    式(2)~式(15)是關(guān)于混合多行程收運路線的約束。其中,式(2)表示第一段的開放行程需要從車場出發(fā);式(3)表示收運車的第一段開放行程中不可以從車場直接到達處理廠;式(4)表示最后一段的空載行程必須返回車場;式(5)和式(6)表示實際派出收運車的數(shù)量不超過最大數(shù)量;式(7)表示收運車需完成當前時段內(nèi)所有收運點的收運任務;式(8)表示收運車不訪問不屬于當前時段的收運點;式(9)表示收運車在每個收運點收運完成后,需離開;式(10)表示各收運點自身無路徑產(chǎn)生;式(11)表示收運車訪問收運點之后,需前往處理廠卸貨,不可以從收運點直接返回車場;式(12)表示每個時段內(nèi)收運車執(zhí)行多行程的最大數(shù)量限制;式(13)~式(15)表示收運車被使用時才會產(chǎn)生多行程路徑。

    式(16)~式(18)為收運車容量相關(guān)的約束。其中,式(16)表示各收運點的垃圾量需一次性收集完;式(17)表示收運車收運的垃圾不超過其最大凈載量;式(18)表示開放行程中,收運車需在車場空車出發(fā),并在到達處理廠后卸空垃圾。

    式(19)~式(26)為時間相關(guān)的約束。其中,式(19)表示加上擁堵時間后的收運車實際行駛時間;式(20)~式(23)表示收運車到達各點的訪問時間需滿足路線訪問順序;式(24)表示收運車的到達時間必須小于各收運點的最晚收運時間窗;式(25)和式(26)為工作平衡性約束。

    2 問題求解與算法設計

    本文建立的模型M為混合整數(shù)規(guī)劃模型,對于小規(guī)模問題,采用AMPL+GUROBI直接求解,得到最優(yōu)解。當問題規(guī)模較大時,模型M無法在有效時間內(nèi)求解。為滿足求解大規(guī)模實際問題的需要,本文設計了混合自適應大鄰域搜索算法(Hybrid Adaptive Large Neighborhood Search algorithm,HALNS),在短時間內(nèi)得到問題的近似最優(yōu)解。

    2.1 啟發(fā)式算法的設計及實現(xiàn)

    自適應大鄰域搜索算法(Adaptive Large Neighborhood Search,ALNS)是在大鄰域搜索算法(Large Neighborhood Search,LNS)的基礎上,加上算子權(quán)重的自適應調(diào)整,即在每次迭代時會根據(jù)當前解的質(zhì)量更新算子的權(quán)重,從而使算法在尋優(yōu)過程中能自主調(diào)整各個算子被選擇的概率。本文提出的HALNS 算法在ALNS 算法的基礎上,加入局部搜索算子,以提升算法的尋優(yōu)能力。HALNS 算法的整體流程如圖2所示。

    圖2 HALNS算法流程Fig.2 Flowchart of HALNS

    2.2 初始解的構(gòu)造及初始解

    根據(jù)問題多行程路線的特點,本文采用三維矩陣對解進行編碼,如圖3所示。所有車的路線組成一個完整的解,每輛車的路線由多個行程組成。為保證解的可行性,需要保證每輛車的行程總數(shù)不超過車輛最大行程數(shù),且每段行程的收運量均不超過收運車的最大容量。

    為獲得初始可行解,采用貪婪算法構(gòu)造初始路線,具體過程如下:首先,按節(jié)點編號順序選擇收運點,以加入空車的行程路線,當該行程中垃圾量達到收運車的最大載量Qmax時,新增下一段行程,并將剩余的收運點增加到新行程中,循環(huán)此過程直至該收運車的總行駛時間達到駕駛員的工作最晚時間;然后,選取新的收運車,并重復上述步驟;直至該時段所有收運點均完成。

    2.3 適應度函數(shù)

    為滿足收運時間窗和駕駛員工作平衡的約束,構(gòu)造帶懲罰函數(shù)的適應度函數(shù)fs,即

    式中:zs——當前解s的目標函數(shù);

    Es,Ls——當前解s中,違反式(25)和式(26)的工作時間;

    η——懲罰系數(shù)(經(jīng)實驗,取η=500),表示當前解s違反時間窗式(25)和式(26)時的單位時間懲罰成本。

    適應度函數(shù)fs的值越小,表明當前解的質(zhì)量越高。

    2.4 鄰域解的更新及算子設計

    根據(jù)問題特性,本文設計了3 種移除算子和2 種插入算子,以獲得新的鄰域解。為產(chǎn)生新的鄰域解,需要先從當前解中移除m個收運點,再將這些點插入到新的位置里。在每次迭代中,移除的收運點數(shù)目m從范圍中隨機產(chǎn)生,其中,|Np|表示當前時段收運點的總數(shù)量,λ1和λ2分別表示要移除的收運點占總收運點數(shù)的下界與上界比例值(0<λ1<λ2<1)。經(jīng)實驗,本文選取λ1=0.10,λ2=0.25。各個算子的規(guī)則如下。

    (1)移除算子

    ①R1隨機移除

    從當前解中隨機移除一個收運點。

    ②R2最差點移除

    移除使當前解適應度值增加最多的收運點。定義當前解s 中收運點i 的邊際成本,其中,表示從當前解s 中移除收運點i 后的適應度值,將最大的點移除。

    ③R3相關(guān)移除

    根據(jù)收運節(jié)點之間的相似度進行移除[12],即先隨機移除一個收運點i,再移除與收運點i相似度最高的收運點j。本文通過距離、交通狀況、時間窗要求及餐飲垃圾量這4方面衡量相似度,收運節(jié)點i和j的相似度R(i,j)為

    式中:dij——收運點i和j之間的距離;

    cij——收運點之間的交通狀況;

    |li-lj|——收運點之間最晚開始服務的時間差異;

    |qi-qj|——收運點之間的需求差異;

    ?、ρ、σ、?——影響因素的權(quán)重。

    R(i,j)越小,表示相似度最大。

    (2)插入算子

    ①INS1貪婪插入

    將從當前時段的收運路線中移除的收運點插入到使得當前解的適應度值增加最少的位置。

    ②INS2最大后悔值準則插入

    令fk表示當前解的適應度值,用fik表示當前解插入收運點i后第k低的適應度值,令ΔFk=fik-fk,當k≤k'時,滿足ΔFk≤ΔFk',定義后悔值ci=ΔF2-ΔF1,表示把收運點i插入到當前解的最優(yōu)位置和次優(yōu)位置的差值,將最大后悔值的收運點插入到當前解中。

    (3)局部搜索算子

    使用移除和插入算子產(chǎn)生新的鄰域解之后,采用局部搜索算子改善產(chǎn)生的鄰域解,改善后的解作為當前解。每次迭代的改進次數(shù)為ω次,每次使用的算子從以下5種隨機選擇。

    ①LS1轉(zhuǎn)移(Shift)

    隨機選擇一輛車,從該車的一段行程中隨機選擇一個收運點,移動到其他行程中。

    ②LS2交換(Swap)

    隨機選擇兩輛車,從已有行程中隨機選擇兩個收運點,交換兩收運點位置。

    ③LS3 2-opt搜索

    從某一行程中隨機選擇兩個收運點,翻轉(zhuǎn)兩點之間的收運點順序。

    ④LS4 2-opt*搜索

    隨機選擇兩輛車,從不同行程路徑中隨機選擇兩個收運點,交換該兩點之后的收運路徑。

    ⑤LS5 平衡性改進(Balance)

    對于工作結(jié)束時間超過駕駛員最晚工作時間的多行程路線,從中隨機選取一段行程,并將其插入到另外一條工作時間較短的多行程路線中,使各個路線的總工作時間盡可能平衡。

    對于經(jīng)過局部搜索改進的鄰域解,需要使用如下規(guī)則判斷是否將其接受為新的當前解。當新解優(yōu)于當前最好解時,則接受新解;否則,根據(jù)模擬退火接受準則,以概率接受新解,其中,f'為新解的適應度值,f為當前解的適應度值,D表示當前溫度,其初始溫度為D=D0,每次迭代后,D通過D=θD更新,0<θ <1,本文選取D0=30,θ=0.94。最后,若算法連續(xù)迭代φ次未能更新解,或達到最大迭代次數(shù)ψ時,算法終止。本文選取φ=30,ψ=300。

    (4)算子權(quán)重的更新

    在迭代下一代鄰域解之前,需要依據(jù)當前解的質(zhì)量更新移除和插入算子的權(quán)重。令po表示算子o被選擇的概率,表示在第t次迭代時算子o的權(quán)重,算子被選擇的概率以及算子權(quán)重的迭代方法為

    設定移除和插入算子每20代更新一次,其中,算子o的初始得分為10;算子更新參數(shù)α=0.7;表示算子o在第t次更新時,在近20代迭代中被選中的次數(shù);表示算子o在第t次更新時,在近20次迭代中得到的解的總評分。評分方式如下:若在迭代后得到的鄰域解snew不劣于當前的全局最好解sbest,得 分40;若sbest≤snew≤sold,則得分30;若snew>sold,但該解被接受,則得分20;否則,得分為0。

    3 算例驗證

    為驗證所建立的數(shù)學模型和提出算法的有效性,采用不同規(guī)模的算例,將GUROBI 求解數(shù)學模型M得到的解與HALNS和ALNS算法得到的解進行比較。

    3.1 算例設置

    本文從Solomon 標準算例庫中選取不同規(guī)模的車輛路徑問題算例,并對原數(shù)據(jù)進行以下調(diào)整,適應求解本文問題:(1)劃分兩個收運時段,將需求點劃分為3種收運方式(上午、下午及上下午),并調(diào)整需求點的時間窗和收運量;(2)原數(shù)據(jù)中僅含有車場,無處理廠,本文選取接近中心的位置虛擬為處理廠位置,并設定車輛在處理廠的卸料時間為5 min;(3)為適應多行程收運,縮減了車容量;(4)原數(shù)據(jù)中不考慮擁堵,故該部分將擁堵時間設置為0;(5)原數(shù)據(jù)中不含收運成本,本文設定每輛車的固定使用成本為200 元·時段-1,燃油成本為1 元·min-1。具體參數(shù)的調(diào)整情況如表1所示。

    表1 算例規(guī)模及相關(guān)參數(shù)Table 1 Instance setting and parameters

    3.2 算例驗證

    算例R1~R5的求解結(jié)果如表2所示,每個算例均采用3 種方法求解,即AMPL+GUROBI 9.5.2 直接求解模型、ALNS 算法和HALNS 算法。其中,ALNS算法采用與HALNS算法相同的初始解產(chǎn)生方案,以及相同的移除和插入算子,但沒有加入局部搜索算子。HALNS 算法中,根據(jù)算例規(guī)模的不同設定不同的局部搜索次數(shù)值(ω),其中,R1和R2中,ω=50;R3 中,ω=100;R4 中,ω=200;R5 中,ω=1000。HALNS 和ALNS 的算法均采用Python編程實現(xiàn)。

    表2 算例求解結(jié)果對比Table 2 Comparison of solutions obtained by different algorithms

    表2 的后兩列給出每個算例求解10 次后得到的平均求解時間和平均目標值??梢钥闯?,對于小規(guī)模算例R1~R3,雖然直接采用GUROBI求解模型M 需要較長時間,但是,可以得到問題的最優(yōu)解或近似最優(yōu)解。HALNS算法和ALNS算法的求解時間較短,但不是每次求解都可以得到最優(yōu)解。相比較之下,HALNS 算法得到最優(yōu)解的次數(shù)多于ALNS 算法,所以,HALNS 算法的平均目標值低于ALNS 算法。當算例規(guī)模增大時,GUROBI 無法在合理時間內(nèi)得到最優(yōu)解,但HALNS 和ALNS 算法都能在短時間內(nèi)得到近似最優(yōu)解,且HALNS 算法在解的質(zhì)量和求解時間上明顯優(yōu)于ALNS算法,所以,HALNS算法更適合求解大規(guī)模問題。綜上,算例求解結(jié)果驗證了本文模型的正確性以及HALNS算法的有效性。

    3.3 駕駛員工作平衡性檢驗

    算例R1~R5 在有平衡性約束和無平衡約束的情況下得到的結(jié)果對比如表3 所示。令ΔTB表示考慮平衡約束時駕駛員的工作時間差(最大工作時長減最小工作時長),ΔTN表示不考慮平衡約束時駕駛員的總工作時間差,當考慮平衡約束時,對駕駛員的工作量平衡性的改善程度為αB;同樣,令CB表示考慮平衡約束時的總成本,CN表示不考慮平衡約束時的路線總成本,考慮平衡約束后總成本的增加程度為αC,即

    表3 工作平衡性對比結(jié)果Table 3 Comparison of workload balancing

    表4 油耗參數(shù)定義與取值Table 4 Definition and value for parameters related to fuel consumption

    從表3 可以看出,當考慮工作量平衡性約束時,總成本雖有小幅增加(小于2%),但駕駛員的工作平衡性能得到大幅度改善,尤其當算例規(guī)模增大時,考慮平衡性約束的優(yōu)勢更為明顯。

    4 實例分析

    為進一步驗證本文的模型和算法在求解實際問題時的有效性,選取大連市中山區(qū)的餐飲垃圾收運網(wǎng)絡作為實際案例進行分析。在實例求解中,根據(jù)實際調(diào)研的結(jié)果以及相關(guān)文獻中的數(shù)據(jù),設定相關(guān)數(shù)據(jù)。

    4.1 實例數(shù)據(jù)

    (1)收運點及產(chǎn)生量

    首先,依據(jù)大連市目前餐飲垃圾的分類實施情況,收運來源僅考慮大中型餐飲企業(yè)和各企事業(yè)單位食堂。本文通過調(diào)研得到大連市中山區(qū)的228 個餐飲垃圾收運點,并將其位置信息導入到地理信息系統(tǒng)軟件(ARCGIS)中,如圖4所示。根據(jù)文獻[13]確定大中型餐飲企業(yè)餐飲垃圾的產(chǎn)生量均值分別為360 kg·d-1和191 kg·d-1,以均值上下浮動10%為取值范圍隨機產(chǎn)生不同餐飲企業(yè)的產(chǎn)生量。企事業(yè)單位食堂的餐飲垃圾產(chǎn)生量根據(jù)各企事業(yè)單位的在職人數(shù)乘以人均日餐飲垃圾產(chǎn)生量(0.23 kg·d-1·人-1)得到。同時,根據(jù)各餐飲企業(yè)的性質(zhì)和營業(yè)高峰期,確定上下午兩個收運時段和兩種收運頻率。其中,營業(yè)高峰期為晚上的餐飲企業(yè)在每天的時段1,即8:00-12:00收運;營業(yè)高峰期在早上和中午的餐飲企業(yè)在每天的時段2,即14:00-18:00收運;全天營業(yè)的餐飲企業(yè)在時段1和時段2各收運一次。

    圖4 大連市中山區(qū)餐飲垃圾收運點示例Fig.4 Illustration of food waste collection points in Zhongshan District,Dalian

    (2)收運車輛油耗及道路擁堵系數(shù)

    在實例求解時,本文根據(jù)實際調(diào)研的車輛數(shù)據(jù),結(jié)合BARTH 等[14]開發(fā)的Comprehensive Emission Model 計算收運車以不同速度行駛的油耗,其中,收運車輛以速度v行駛t段時間的油耗為

    另外,本文對速度和載重量對于油耗的影響進行靈敏度分析發(fā)現(xiàn),餐飲垃圾收運車的油耗量對速度的敏感性較高,而對載重量的敏感性較低。因此,本文按照收運車輛在城區(qū)和非城區(qū)的實際行駛路程采用兩種速度。根據(jù)文獻[5,15],并結(jié)合實際調(diào)研,設定收運車輛在城區(qū)的運行速度v為30 km·h-1,在非城區(qū)的行駛速度v為70 km·h-1。當收運車遇到擁堵路段或在收運點和垃圾處理廠進行作業(yè)時,收運車處在怠速狀態(tài),根據(jù)文獻[16]設定收運車輛怠速時油耗FC為2.20 L·h-1。收運車的凈載量設定為2.5 t,車輛固定成本為500元·時段-1(每個時段為4 h),柴油價格為7.2元·L-1。

    為確定主要道路的擁堵系數(shù),首先,通過高德地圖獲取大連市中山區(qū)的歷史交通擁堵數(shù)據(jù),并通過計算得到主干道路在各個時段的擁堵系數(shù),如表5 所示;然后,通過擁堵系數(shù)得到車輛在路線上的擁堵時間(擁堵時間為兩點距離乘擁堵系數(shù))。

    表5 常發(fā)生擁堵的主干路及擁堵系數(shù)Table 5 Congestion coefficients of some main roads

    4.2 結(jié)果分析

    采用HALNS算法求解實例得到的車輛路線如表6所示,其中,HALNS 算法每次迭代的局部搜索次數(shù)ω設定為1000,連續(xù)迭代未更新解次數(shù)φ=30,最大迭代次數(shù)ψ=200。

    表6 HALNS算法求解實例得到的收運路線Table 6 Food collection routes obtained by HALNS for real case

    HALNS算法與當前某餐飲垃圾智能收運系統(tǒng)中采用的優(yōu)化算法(BAU)求解結(jié)果的對比情況如表7所示。其中,BAU中采用兩階段方法,首先,將收運區(qū)域劃分為6 個子區(qū)域;然后,求解一個旅行商問題,得到各個區(qū)域的收運路線。同時,當前實際優(yōu)化算法設定每個收運點的收運頻率均為每天一次,也未考慮駕駛員的工作量平衡情況以及道路擁堵情況。

    表7 餐飲垃圾收運路線優(yōu)化結(jié)果對比Table 7 Comparison of food waste collection route solutions

    表7 的對比結(jié)果表明,本文提出的HALNS 算法能對收運路線進行更合理的安排,不僅可以降低14.3%的總收運成本,還能在少用1 輛車的情況下減少2%的總工作時間,使駕駛員的工作平衡性提升57.3%。另外,從環(huán)境角度,通過在模型中考慮擁堵時間,使求解結(jié)果能更有效地規(guī)避擁堵路段,使收運過程中的碳排放量降低12.7%。

    5 結(jié)論

    本文拓展了對于車輛路線優(yōu)化方面的研究,為城市餐飲垃圾的收運提供了降本增效的實際解決方案,對于提升城市管理水平,加快智慧城市的建設,促進經(jīng)濟、社會和環(huán)境的可持續(xù)發(fā)展具有重要意義。本文得到的主要結(jié)論如下:

    (1)所建立的模型M 和提出的算法HALNS 能在降低成本的同時提高駕駛員工作量的平衡性。首先,通過將收運點劃分到不同收運時段,均衡各個時段的垃圾收運量;然后,在優(yōu)化每個時段的收運路線時,考慮駕駛員的工作量平衡約束。優(yōu)化結(jié)果表明,在成本增加值不超過2%的情況下,可以將工作量的平衡性提高36%~98%,而且算例的規(guī)模越大,對于不平衡性的改善效果越明顯。

    (2)采用模型和算法能有效求解具有時間窗約束的多行程餐飲垃圾收運問題。對于小規(guī)模問題,直接求解所建立的混合整數(shù)規(guī)劃模型能得到問題的最優(yōu)解;對于中大規(guī)模算例,采用提出的HALNS算法能在較短時間內(nèi)得到問題的近似最優(yōu)解。

    (3)研究思路及算法適用于實際問題的求解與優(yōu)化。以大連市中山區(qū)的餐飲垃圾收運為例,考慮實際路網(wǎng)的道路擁堵情況和駕駛員的工作平衡性約束等現(xiàn)實因素時,相比當前實際應用的優(yōu)化算法,HALNS算法能在經(jīng)濟方面降低14.3%的總收運成本,環(huán)境方面降低12.7%的碳排放量,社會方面減少駕駛員總的工作時長,并將工作量的不平衡性改善57.3%。

    猜你喜歡
    收運時段算子
    基于物聯(lián)網(wǎng)的智慧垃圾收運系統(tǒng)分析
    2025年山西垃圾收運覆蓋90%以上自然村
    蘇州工業(yè)園區(qū)餐廚垃圾產(chǎn)生現(xiàn)狀及收運方案研究
    擬微分算子在Hp(ω)上的有界性
    各向異性次Laplace算子和擬p-次Laplace算子的Picone恒等式及其應用
    農(nóng)村生活垃圾收運員量化考核指標
    一類Markov模算子半群與相應的算子值Dirichlet型刻畫
    四個養(yǎng)生黃金時段,你抓住了嗎
    Roper-Suffridge延拓算子與Loewner鏈
    傍晚是交通事故高發(fā)時段
    国产精品九九99| 巨乳人妻的诱惑在线观看| 国产淫语在线视频| 俄罗斯特黄特色一大片| 欧美日韩一级在线毛片| 午夜免费鲁丝| 精品国产超薄肉色丝袜足j| 精品国产一区二区三区四区第35| 51午夜福利影视在线观看| 视频区图区小说| 久久国产精品人妻蜜桃| 亚洲欧美激情在线| kizo精华| 国产成人影院久久av| 国产精品亚洲av一区麻豆| 人人妻人人添人人爽欧美一区卜| 12—13女人毛片做爰片一| 久久国产精品男人的天堂亚洲| 久久国产精品大桥未久av| 免费观看a级毛片全部| 久久久久精品国产欧美久久久| 欧美人与性动交α欧美软件| 国产在线观看jvid| 亚洲九九香蕉| 久久久精品国产亚洲av高清涩受| av国产精品久久久久影院| 久久久精品区二区三区| 亚洲精品美女久久av网站| 丝袜人妻中文字幕| √禁漫天堂资源中文www| 99国产精品一区二区三区| 精品亚洲乱码少妇综合久久| 肉色欧美久久久久久久蜜桃| 精品欧美一区二区三区在线| 国产精品久久久久久精品电影小说| 丝袜喷水一区| 国产成人精品久久二区二区91| 桃花免费在线播放| 超碰成人久久| 国产又色又爽无遮挡免费看| 亚洲国产欧美一区二区综合| 欧美另类亚洲清纯唯美| 一二三四在线观看免费中文在| 巨乳人妻的诱惑在线观看| 香蕉丝袜av| 99热网站在线观看| 纵有疾风起免费观看全集完整版| 在线av久久热| 在线 av 中文字幕| 久久久欧美国产精品| 久久免费观看电影| 伦理电影免费视频| 不卡一级毛片| 女人久久www免费人成看片| 国产欧美日韩一区二区精品| 99国产精品99久久久久| 老熟妇仑乱视频hdxx| 亚洲精品美女久久久久99蜜臀| 欧美亚洲 丝袜 人妻 在线| 成人18禁在线播放| 女同久久另类99精品国产91| 精品人妻熟女毛片av久久网站| 汤姆久久久久久久影院中文字幕| 高潮久久久久久久久久久不卡| 国产亚洲一区二区精品| 99在线人妻在线中文字幕 | 汤姆久久久久久久影院中文字幕| 日本av手机在线免费观看| 动漫黄色视频在线观看| videosex国产| 久久精品国产99精品国产亚洲性色 | 十八禁网站免费在线| 在线天堂中文资源库| 一个人免费在线观看的高清视频| 久久免费观看电影| 亚洲人成伊人成综合网2020| 国产高清国产精品国产三级| 亚洲av片天天在线观看| 精品福利观看| 久久久久久久大尺度免费视频| 亚洲精品在线观看二区| 日韩精品免费视频一区二区三区| 亚洲久久久国产精品| 久久精品人人爽人人爽视色| a在线观看视频网站| 日本一区二区免费在线视频| 亚洲欧洲精品一区二区精品久久久| 久久青草综合色| 麻豆乱淫一区二区| 国产在线免费精品| 少妇的丰满在线观看| 色婷婷久久久亚洲欧美| kizo精华| 精品国产一区二区久久| 免费日韩欧美在线观看| 久久中文看片网| 久久性视频一级片| 99国产精品99久久久久| 久久久国产欧美日韩av| 99国产精品一区二区蜜桃av | 9热在线视频观看99| 美女午夜性视频免费| 国产三级黄色录像| 女人久久www免费人成看片| 性高湖久久久久久久久免费观看| 亚洲九九香蕉| 欧美人与性动交α欧美精品济南到| 一本色道久久久久久精品综合| 两个人看的免费小视频| 我要看黄色一级片免费的| 老司机靠b影院| 精品亚洲乱码少妇综合久久| 成年动漫av网址| 国产精品二区激情视频| 国产精品自产拍在线观看55亚洲 | 高清黄色对白视频在线免费看| 国产男女内射视频| 老司机在亚洲福利影院| 精品少妇黑人巨大在线播放| 桃花免费在线播放| 黄网站色视频无遮挡免费观看| 欧美日本中文国产一区发布| 久久精品国产99精品国产亚洲性色 | 免费黄频网站在线观看国产| 高清毛片免费观看视频网站 | 最黄视频免费看| 亚洲成人手机| 久久狼人影院| 久久性视频一级片| 免费av中文字幕在线| 久久午夜亚洲精品久久| 国产伦人伦偷精品视频| 免费黄频网站在线观看国产| 午夜福利一区二区在线看| 成年人免费黄色播放视频| 精品国产乱码久久久久久男人| 成年人黄色毛片网站| 成人三级做爰电影| netflix在线观看网站| 十八禁高潮呻吟视频| 男女下面插进去视频免费观看| 亚洲第一欧美日韩一区二区三区 | 美女扒开内裤让男人捅视频| 免费一级毛片在线播放高清视频 | 999久久久国产精品视频| 久久国产精品男人的天堂亚洲| 久久精品aⅴ一区二区三区四区| 国产不卡一卡二| 18禁国产床啪视频网站| 99在线人妻在线中文字幕 | 久久九九热精品免费| www.精华液| 咕卡用的链子| 国产不卡av网站在线观看| 免费观看a级毛片全部| 丰满少妇做爰视频| a级片在线免费高清观看视频| 亚洲一卡2卡3卡4卡5卡精品中文| 色老头精品视频在线观看| 午夜福利乱码中文字幕| 男人操女人黄网站| 曰老女人黄片| 久久精品人人爽人人爽视色| av网站在线播放免费| 黄网站色视频无遮挡免费观看| 日韩熟女老妇一区二区性免费视频| 老司机午夜福利在线观看视频 | 国产男女超爽视频在线观看| 中文欧美无线码| 国产麻豆69| 国产极品粉嫩免费观看在线| 法律面前人人平等表现在哪些方面| 久久影院123| 中文字幕精品免费在线观看视频| 欧美+亚洲+日韩+国产| 90打野战视频偷拍视频| 久久久精品94久久精品| 免费看a级黄色片| 亚洲国产欧美一区二区综合| 脱女人内裤的视频| 99国产精品一区二区三区| 最新的欧美精品一区二区| 人人妻人人澡人人看| av网站在线播放免费| 91精品国产国语对白视频| 麻豆成人av在线观看| 少妇精品久久久久久久| 90打野战视频偷拍视频| 免费观看a级毛片全部| 精品第一国产精品| 亚洲va日本ⅴa欧美va伊人久久| 精品人妻在线不人妻| 精品国产乱码久久久久久小说| 在线十欧美十亚洲十日本专区| 亚洲成av片中文字幕在线观看| 99久久精品国产亚洲精品| 又紧又爽又黄一区二区| 国产男女超爽视频在线观看| 日本av手机在线免费观看| 免费少妇av软件| 久久精品亚洲av国产电影网| 一本综合久久免费| 19禁男女啪啪无遮挡网站| 精品午夜福利视频在线观看一区 | 十八禁网站免费在线| 成年女人毛片免费观看观看9 | 97在线人人人人妻| 欧美日韩视频精品一区| 水蜜桃什么品种好| 肉色欧美久久久久久久蜜桃| 美女福利国产在线| 在线永久观看黄色视频| 一本色道久久久久久精品综合| 动漫黄色视频在线观看| 交换朋友夫妻互换小说| netflix在线观看网站| 1024香蕉在线观看| 亚洲精品国产一区二区精华液| 午夜精品国产一区二区电影| 999精品在线视频| 男女免费视频国产| 亚洲欧美一区二区三区久久| 亚洲成人国产一区在线观看| 亚洲精品国产一区二区精华液| 大型黄色视频在线免费观看| 国产成人影院久久av| 天天躁日日躁夜夜躁夜夜| 1024香蕉在线观看| 午夜久久久在线观看| 天天躁夜夜躁狠狠躁躁| 日本av免费视频播放| 99国产精品一区二区三区| 国产欧美日韩精品亚洲av| 老司机福利观看| 天天躁夜夜躁狠狠躁躁| 久久午夜亚洲精品久久| 亚洲五月色婷婷综合| 国产一区二区三区在线臀色熟女 | 熟女少妇亚洲综合色aaa.| 久久久久久久国产电影| 中亚洲国语对白在线视频| 王馨瑶露胸无遮挡在线观看| 欧美黄色片欧美黄色片| 高清在线国产一区| 巨乳人妻的诱惑在线观看| 国产精品美女特级片免费视频播放器 | 久久久水蜜桃国产精品网| 国产激情久久老熟女| 亚洲av成人不卡在线观看播放网| 亚洲熟女精品中文字幕| 岛国在线观看网站| 少妇裸体淫交视频免费看高清 | 日韩有码中文字幕| 亚洲中文字幕日韩| 男女边摸边吃奶| 国产在视频线精品| 一个人免费在线观看的高清视频| 窝窝影院91人妻| 他把我摸到了高潮在线观看 | 欧美老熟妇乱子伦牲交| 多毛熟女@视频| 极品少妇高潮喷水抽搐| 精品国产乱码久久久久久小说| 亚洲欧美精品综合一区二区三区| 男女下面插进去视频免费观看| 国产成人精品久久二区二区免费| 99精品欧美一区二区三区四区| 777米奇影视久久| 亚洲专区中文字幕在线| 午夜精品久久久久久毛片777| 国产精品98久久久久久宅男小说| 久久99一区二区三区| 免费在线观看日本一区| 午夜福利乱码中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 日本黄色视频三级网站网址 | 男女高潮啪啪啪动态图| 丁香欧美五月| 极品人妻少妇av视频| 精品国产亚洲在线| 欧美日韩成人在线一区二区| 国产精品免费视频内射| 免费女性裸体啪啪无遮挡网站| 国产免费视频播放在线视频| 亚洲久久久国产精品| svipshipincom国产片| 国产亚洲精品第一综合不卡| 99热网站在线观看| av天堂在线播放| 久久香蕉激情| 每晚都被弄得嗷嗷叫到高潮| 久久午夜亚洲精品久久| 国产无遮挡羞羞视频在线观看| 欧美精品一区二区大全| 国产在视频线精品| 久久人妻福利社区极品人妻图片| 久久精品国产综合久久久| 老汉色av国产亚洲站长工具| 国产人伦9x9x在线观看| 一进一出好大好爽视频| 久久久久精品国产欧美久久久| 啦啦啦免费观看视频1| 免费不卡黄色视频| 国产日韩欧美亚洲二区| 日韩一卡2卡3卡4卡2021年| 999精品在线视频| 丝袜人妻中文字幕| 亚洲成人免费电影在线观看| 黄片大片在线免费观看| 成人特级黄色片久久久久久久 | 成人永久免费在线观看视频 | 色老头精品视频在线观看| 叶爱在线成人免费视频播放| netflix在线观看网站| 亚洲精品国产精品久久久不卡| 黑人欧美特级aaaaaa片| 黄片小视频在线播放| 法律面前人人平等表现在哪些方面| 国产精品麻豆人妻色哟哟久久| 99久久国产精品久久久| 99久久人妻综合| 一进一出抽搐动态| 好男人电影高清在线观看| 欧美精品一区二区免费开放| 国产精品久久久久久精品古装| 国产成人影院久久av| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品中文字幕一二三四区 | 日韩大码丰满熟妇| av免费在线观看网站| 久久99一区二区三区| 亚洲自偷自拍图片 自拍| 国产一区二区三区在线臀色熟女 | 色视频在线一区二区三区| 他把我摸到了高潮在线观看 | 五月开心婷婷网| 亚洲精品久久午夜乱码| 十八禁网站网址无遮挡| 中文字幕av电影在线播放| 高清视频免费观看一区二区| 免费在线观看日本一区| 国产一区二区激情短视频| 精品国产乱子伦一区二区三区| 久久中文看片网| 中文字幕最新亚洲高清| 久久中文看片网| 国产精品一区二区在线不卡| 人妻一区二区av| 国产aⅴ精品一区二区三区波| 久久毛片免费看一区二区三区| 美女扒开内裤让男人捅视频| av视频免费观看在线观看| 久久久久国产一级毛片高清牌| 女人高潮潮喷娇喘18禁视频| 一本—道久久a久久精品蜜桃钙片| 精品亚洲成国产av| 一本色道久久久久久精品综合| 91成人精品电影| 老司机午夜十八禁免费视频| 国产日韩欧美亚洲二区| 老司机午夜福利在线观看视频 | 亚洲国产av影院在线观看| 韩国精品一区二区三区| 国产区一区二久久| 啦啦啦免费观看视频1| 精品欧美一区二区三区在线| 一区福利在线观看| av天堂久久9| 亚洲第一青青草原| 免费在线观看日本一区| 欧美日韩中文字幕国产精品一区二区三区 | 少妇粗大呻吟视频| 女警被强在线播放| 男女午夜视频在线观看| 国产精品亚洲av一区麻豆| 91成人精品电影| 18禁裸乳无遮挡动漫免费视频| 精品亚洲成国产av| 2018国产大陆天天弄谢| 国产不卡av网站在线观看| 搡老乐熟女国产| 精品一区二区三区av网在线观看 | 亚洲精品国产一区二区精华液| 国产精品香港三级国产av潘金莲| 丰满人妻熟妇乱又伦精品不卡| 在线永久观看黄色视频| 50天的宝宝边吃奶边哭怎么回事| 成年人免费黄色播放视频| 啦啦啦视频在线资源免费观看| 国产精品.久久久| 亚洲国产欧美网| 亚洲人成77777在线视频| 丰满少妇做爰视频| 国产成人精品无人区| 97人妻天天添夜夜摸| 中亚洲国语对白在线视频| 亚洲五月色婷婷综合| 大香蕉久久成人网| 免费在线观看影片大全网站| 精品久久久久久电影网| 天天影视国产精品| 国产片内射在线| 久久中文看片网| 高潮久久久久久久久久久不卡| 啦啦啦 在线观看视频| 淫妇啪啪啪对白视频| 一二三四社区在线视频社区8| 久久青草综合色| 久久久国产精品麻豆| 亚洲午夜理论影院| 日本a在线网址| 在线 av 中文字幕| 日韩大码丰满熟妇| 纯流量卡能插随身wifi吗| 午夜老司机福利片| tube8黄色片| 国产精品免费大片| 免费在线观看完整版高清| 淫妇啪啪啪对白视频| 国产精品九九99| 99九九在线精品视频| 黄频高清免费视频| 黑人巨大精品欧美一区二区mp4| 久久精品亚洲av国产电影网| 免费女性裸体啪啪无遮挡网站| 亚洲av成人一区二区三| 天堂中文最新版在线下载| 欧美激情久久久久久爽电影 | 亚洲精品国产区一区二| 精品久久久久久电影网| 天天影视国产精品| 一级毛片电影观看| 在线天堂中文资源库| 亚洲av成人一区二区三| 亚洲七黄色美女视频| 亚洲色图av天堂| 最新的欧美精品一区二区| 老汉色∧v一级毛片| 欧美日韩一级在线毛片| 色在线成人网| av一本久久久久| 一二三四社区在线视频社区8| 热re99久久国产66热| 免费高清在线观看日韩| 下体分泌物呈黄色| 制服诱惑二区| 老熟妇乱子伦视频在线观看| 亚洲av日韩精品久久久久久密| 极品人妻少妇av视频| 国产成人啪精品午夜网站| 飞空精品影院首页| 欧美精品人与动牲交sv欧美| 亚洲av片天天在线观看| av视频免费观看在线观看| 久久久久国产一级毛片高清牌| 午夜福利乱码中文字幕| 淫妇啪啪啪对白视频| 国产免费现黄频在线看| 成年人午夜在线观看视频| 亚洲人成电影观看| 韩国精品一区二区三区| 亚洲国产欧美在线一区| 亚洲第一av免费看| 热99re8久久精品国产| 在线 av 中文字幕| 国产伦人伦偷精品视频| 欧美人与性动交α欧美软件| 国产精品.久久久| 国产1区2区3区精品| 精品少妇黑人巨大在线播放| 最近最新中文字幕大全电影3 | 人人妻人人添人人爽欧美一区卜| 少妇被粗大的猛进出69影院| av一本久久久久| 亚洲国产看品久久| 久久久国产一区二区| 久久精品熟女亚洲av麻豆精品| 婷婷丁香在线五月| 18禁观看日本| 美女主播在线视频| 日本黄色视频三级网站网址 | 日韩大码丰满熟妇| a级毛片黄视频| 欧美日韩亚洲综合一区二区三区_| 日日爽夜夜爽网站| 欧美黑人欧美精品刺激| 国产片内射在线| 精品人妻在线不人妻| 亚洲avbb在线观看| 天天影视国产精品| 女警被强在线播放| 婷婷成人精品国产| 久久精品亚洲av国产电影网| 女同久久另类99精品国产91| 1024视频免费在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久亚洲真实| 亚洲五月色婷婷综合| 一进一出抽搐动态| 国产高清激情床上av| 日日夜夜操网爽| 男女免费视频国产| 老司机靠b影院| 可以免费在线观看a视频的电影网站| 蜜桃在线观看..| 亚洲av成人一区二区三| 免费看十八禁软件| 人人妻,人人澡人人爽秒播| 熟女少妇亚洲综合色aaa.| 一本一本久久a久久精品综合妖精| 王馨瑶露胸无遮挡在线观看| 亚洲精品自拍成人| 日韩免费av在线播放| 18禁裸乳无遮挡动漫免费视频| 宅男免费午夜| 国产在视频线精品| 国产色视频综合| 新久久久久国产一级毛片| 一级毛片电影观看| 国产男女超爽视频在线观看| 黑人巨大精品欧美一区二区mp4| 久久这里只有精品19| 午夜福利一区二区在线看| 午夜福利视频精品| 久久天躁狠狠躁夜夜2o2o| 韩国精品一区二区三区| 在线播放国产精品三级| 成在线人永久免费视频| 亚洲欧美激情在线| 首页视频小说图片口味搜索| 亚洲成av片中文字幕在线观看| 高清欧美精品videossex| 成人手机av| 欧美精品高潮呻吟av久久| 成人18禁高潮啪啪吃奶动态图| 俄罗斯特黄特色一大片| 精品国产亚洲在线| 久久久久国产一级毛片高清牌| 十八禁网站免费在线| 欧美精品啪啪一区二区三区| 99香蕉大伊视频| 久久狼人影院| 成人精品一区二区免费| 亚洲av第一区精品v没综合| 亚洲精品在线美女| 亚洲第一欧美日韩一区二区三区 | 中文字幕最新亚洲高清| 叶爱在线成人免费视频播放| 精品国产超薄肉色丝袜足j| 少妇的丰满在线观看| 一二三四社区在线视频社区8| 欧美精品啪啪一区二区三区| 最黄视频免费看| 欧美黑人精品巨大| 久久热在线av| 国产有黄有色有爽视频| 妹子高潮喷水视频| 大码成人一级视频| 日韩三级视频一区二区三区| 国产精品一区二区精品视频观看| 极品人妻少妇av视频| 日本欧美视频一区| 9色porny在线观看| 国产福利在线免费观看视频| 男人舔女人的私密视频| 久久ye,这里只有精品| av有码第一页| 久久国产精品影院| 色婷婷av一区二区三区视频| 国产精品免费一区二区三区在线 | 黄片播放在线免费| 免费观看人在逋| 老司机靠b影院| 欧美中文综合在线视频| 99九九在线精品视频| 久久99一区二区三区| 黑人巨大精品欧美一区二区mp4| 亚洲国产av新网站| 午夜福利,免费看| 欧美精品啪啪一区二区三区| 18禁黄网站禁片午夜丰满| 久久久欧美国产精品| 亚洲专区字幕在线| 亚洲天堂av无毛| 日韩欧美一区视频在线观看| 亚洲av片天天在线观看| 国产色视频综合| 久久精品熟女亚洲av麻豆精品| 色综合婷婷激情| 成人手机av| 精品视频人人做人人爽| 精品国产国语对白av| 欧美成狂野欧美在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲国产欧美网| 中文欧美无线码| 一区福利在线观看| 天天躁日日躁夜夜躁夜夜| 亚洲天堂av无毛| 考比视频在线观看| 无限看片的www在线观看| 最黄视频免费看| 蜜桃在线观看..| 在线十欧美十亚洲十日本专区| 黄色片一级片一级黄色片| 国产不卡av网站在线观看| 欧美性长视频在线观看| 黄色片一级片一级黄色片| 免费在线观看完整版高清| 欧美精品啪啪一区二区三区| 国产熟女午夜一区二区三区| 另类亚洲欧美激情| 中文字幕精品免费在线观看视频| 国产欧美日韩综合在线一区二区| 91麻豆精品激情在线观看国产 |