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

    考慮碳排放的兩階段選址-路徑問題及其算法

    2023-11-03 03:43:06湯希峰
    西南交通大學(xué)學(xué)報 2023年5期
    關(guān)鍵詞:物流園區(qū)算例貨車

    湯希峰 ,何 杰 ,張 浩

    (1.河海大學(xué)土木與交通學(xué)院,江蘇 南京 210098;2.東南大學(xué)交通學(xué)院,江蘇 南京 211189)

    城市是生產(chǎn)、生活的主要集中地,也是物流活動最為集中的地方.兩層的物流網(wǎng)絡(luò)作為一種可行的城市物流解決方案正被國內(nèi)外許多大中城市所采用[1-4],其核心問題之一是兩階段選址-路徑問題 (2ELRP).經(jīng)典的2E-LRP 可以描述為:重型或中型貨車負責(zé)將貨物從一個位于城郊的物流園區(qū)運送至若干位于城市中心區(qū)周邊的配送中心,貨物在這些配送中心經(jīng)過重新配裝后再由輕型或微型貨車送達終端客戶,其目標(biāo)是尋求最優(yōu)的配送中心選址以及第一階段和第二階段貨車的最優(yōu)路徑以實現(xiàn)總的物流成本最小化[5].

    眾所周知,貨車相較于客車容易導(dǎo)致交通事故和交通擁堵,更是污染物和碳排放的主要貢獻者.以我國為例,2019 年全國貨車保有量雖然只占汽車總量的10.71%,貨車CO、HC、NOx、PM 排放量卻占到了汽車排放總量的29.7%、26.3%、83.5%、90.1%[6].隨著生態(tài)環(huán)境保護問題日漸受到重視,如何優(yōu)化2 層網(wǎng)絡(luò)結(jié)構(gòu)的物流運作以減少貨車的碳排放引起了國內(nèi)外學(xué)者越來越多的研究興趣.Govindan 等[7]圍繞生鮮配送,建立了以物流成本最小化和環(huán)境影響最小化為目標(biāo)的2E-LRP 多目標(biāo)優(yōu)化模型,并設(shè)計了多目標(biāo)粒子群優(yōu)化算法(MOPSO)和改進的多目標(biāo)可變鄰域搜索算法(AMOVNS)2 種求解算法;胡大偉等[8]針對燃油車和電動車不同的行駛與排放特性,分別建立了燃油車和電動車的2 階段開放式選址-路徑問題模型,并提出了一種改進的模擬退火算法;劉成清等[9]以經(jīng)典的2 階段開放式選址-路徑問題模型為基礎(chǔ),通過考慮車輛的燃油消耗和CO2排放以及車輛路徑選擇的靈活性,建立了基于路徑靈活性的兩階段開放式低碳選址-路徑問題模型,并運用Cplex進行了求解;Pitakaso 等[10]以燃油消耗最小化為目標(biāo),建立了考慮碳排放的綠色兩階段選址路徑問題(G2ELRP)優(yōu)化模型,并設(shè)計了相應(yīng)的可變鄰域策略的自適應(yīng)搜索(VANSAS)算法.

    一方面,鑒于車輛的碳排放受到車型、車速、載重、行駛距離以及駕駛行為等諸多因素的影響[11-12],現(xiàn)有文獻大多基于車輛的燃油消耗進行碳排放的計算,所需參數(shù)較多,實用性不強;另一方面,現(xiàn)有的用于求解各種2E-LRP 的算法多適用于求解小規(guī)模問題 (客戶數(shù)量n≤100, 候選配送中心數(shù)量m<10),求解大規(guī)模問題 (n>100,m≥10) 時計算耗時較長,很難實際應(yīng)用.為此,本文基于一種簡單實用的以排放因子為主要參數(shù)的碳排放計算方法,建立了考慮碳排放的2E-LRP 優(yōu)化模型,并針對所研究問題的計算復(fù)雜性,設(shè)計了一種可快速求解大規(guī)模問題的兩階段混合算法.

    1 模型構(gòu)建

    1.1 問題描述

    基于圖論,考慮碳排放的2E-LRP 可以描述為G={V,E},V=V0∪Vs∪Vc,V0為物流園區(qū),Vs= (1,2,…,m)為配送中心候選集,Vc= (m+1,m+2,…,m+n)為客戶集;E=E1∪E2,E1= {(i,j):i,j∈V0∪Vs;i≠j}為物流園區(qū)與配送中心之間的路線集,E2= {(i,j):i,j∈Vs∪Vc;i,j?Vs×Vs;i≠j} 為配送中心與客戶以及客戶與客戶之間的路線集;H為可供物流園區(qū)使用的重型或中型貨車的集合,每臺重型或中型貨車具有相同的載重能力,記為QH,滿載和空載時的碳排放系數(shù)分別為εvfH、εveH(單位: kg/km, 下同);K為可供配送中心共用的輕型或微型貨車的集合,每臺輕型或微型貨車也具有相同的載重能力,記為QK(QK

    圖1 2E-LRP 示意Fig.1 An illustration of the 2E-LRP

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

    若用qijh和qijk分別為貨車h(h∈H)在路線 (i,j)∈E1、貨車k(k∈K)在路線(i,j)∈E2上行駛時的載重量;決策變量xijh=1 表示貨車h(h∈H) 經(jīng)過路線 (i,j)∈E1,否則xijh=0;yijk=1 表示貨車k(k∈K)經(jīng)過路線 (i,j)∈E2,否則yijk=0;zi=1 表示配送中心i(i∈Vs) 被選用,否則zi=0;uij=1 表示客戶j由配送中心i提供服務(wù),否則uij=0;則以碳排放最小化為目標(biāo)的2E-LRP 數(shù)學(xué)模型可表示為

    式 (1) 為目標(biāo)函數(shù),表示從物流園區(qū)到配送中心以及配送中心到客戶2 個階段的貨車碳排放量之和最小,與經(jīng)典2E-LRP 模型的目標(biāo)函數(shù) (關(guān)于車輛通行成本或車輛行駛距離的函數(shù)) 不同,是一個關(guān)于車輛載重和車輛行駛距離的函數(shù);式(2)、(3) 表示2 個階段中某一貨車要么不被使用,要么至多被使用一次;式(4)、(5) 保證2 個階段中貨車路線的連續(xù)性,同時保證貨車完成任務(wù)后返回原先出發(fā)的物流園區(qū)或配送中心;式(6) 表示客戶只能由一個配送中心提供服務(wù);式(7) 表示客戶若由某一配送中心提供服務(wù),則必被從該配送中心出發(fā)的貨車訪問一次;式(8) 表示某一配送中心要么不被選用,若被選用,所服務(wù)的客戶總需求量不能超過其作業(yè)能力;式(9)、(10) 表示2 個階段中貨車服務(wù)的客戶需求均不能超過其載重能力;式(11)、(12) 為2 個階段中路線的完整性約束;式(13)、(14)表示2 個階段中貨車載重量在相鄰路段上的數(shù)量關(guān)系;式(15) ~(18)為決策變量的取值約束.

    2 算法設(shè)計

    由于設(shè)施選址問題和車輛路徑問題都屬于NPhard 問題,所以本文研究的考慮碳排放的2E-LRP是更為復(fù)雜的NP-hard 問題,求解難度很大,尤其是隨著問題規(guī)模的增大,計算耗時呈指數(shù)增加.鑒于此,本文設(shè)計了一種可快速求解大規(guī)模問題的2 階段混合算法(TSHA ),算法流程如圖2 所示.第1 階段,通過將2E-LRP 轉(zhuǎn)化成不考慮車輛路徑的2 階段設(shè)施選址問題 (2E-FLP),直接求解配送中心選址方案和客戶分配方案;第2 階段,基于得到的配送中心選址和客戶分配方案,物流園區(qū)到被選用的配送中心以及每個配送中心到所分配客戶的車輛路徑問題就變成了獨立且較易求解的VRP 問題,進而加以解決.

    圖2 兩階段混合算法的算法流程Fig.2 Flowchart of the TSHA

    2.1 第1 階段算法

    由于迭代算法在求解過程中非常低效,為提高所建模型的求解效率,TSHA 在第1 階段先將考慮碳排放的2E-LRP 轉(zhuǎn)化成不考慮車輛路徑的2EFLP,模型如下:

    式 (19) 為目標(biāo)函數(shù),表示客戶到配送中心以及配送中心到物流園區(qū)的以需求量為權(quán)重的距離之和最小,約束條件對應(yīng)原問題模型中的式 (6)、(8)、式(15)、(16).不難看出,上述2E-FLP 模型為線性規(guī)劃模型,可直接運用Cplex 求解器進行求解,從而得到配送中心選址和客戶分配方案.

    2.2 第2 階段算法

    基于得到的配送中心選址和客戶分配方案,物流園區(qū)到被選用的配送中心以及每個被選用配送中心到所分配客戶的車輛路徑問題都變成了獨立的VRP 問題.為求解以碳排放最小化為目標(biāo)的VRP 問題,TSHA 在第2 階段采用了一種改進的蟻群算法,如圖3 所示.

    圖3 改進蟻群算法的偽代碼Fig.3 Pseudocode of the improved ant colony algorithm

    與經(jīng)典的蟻群算法相同[13],改進的蟻群算法也包含2 個關(guān)鍵步驟:1) 主要負責(zé)路線的構(gòu)建,即基于當(dāng)前節(jié)點逐步選擇下一個將要訪問的節(jié)點;2) 主要負責(zé)信息素濃度的更新,具體包括局部信息素更新和全局信息素更新.

    從式 (1) 碳排放的計算式可以看出,貨車的碳排放與其載重量和行駛距離密切相關(guān).為此,改進的蟻群算法采用了一種新的路線構(gòu)建規(guī)則,如式 (20) 、式(21) 所示.

    式中:Ω為尚未被訪問的節(jié)點集,Ω∈Vs或Vc;τij為當(dāng)前節(jié)點i與未到訪節(jié)點j之間路線上的信息素濃度;ηij為節(jié)點i與節(jié)點j之間距離的倒數(shù);α、β為重要度系數(shù);q0為給定的常數(shù);q為 [0, 1] 上的隨機數(shù).

    當(dāng)隨機產(chǎn)生的q≤q0時,選擇使取值最大的j作為下一個將要訪問的節(jié)點;當(dāng)q>q0時,則根據(jù)的值所占的比例選擇下一個將要訪問的節(jié)點j,在TSHA 算法中采用輪盤賭的方式進行選擇.

    每只螞蟻都按照上述的偽比例選擇規(guī)則構(gòu)建路線,當(dāng)螞蟻構(gòu)建出一條路線(即達到貨車的載重限制時,返回原先出發(fā)的物流園區(qū)或配送中心),然后重新出發(fā)構(gòu)建下一條路線,如此反復(fù),直至遍訪所有客戶.當(dāng)螞蟻構(gòu)建出一個完整的路線方案后,對該路線方案中的每條路線進行局部信息素更新,第t+ 1 次信息素濃度為

    式中:ρ(0 ≤ρ≤1) 為信息素的揮發(fā)系數(shù);τ0為初始的信息素濃度,取值為1/(mEnns),Enns為基于最近鄰搜索算法得到的貨車碳排放量.

    當(dāng)蟻群中每只螞蟻均構(gòu)建出一個完整的路線方案后,從中選擇一個碳排放量最小的路線方案作為當(dāng)次最優(yōu)的路線方案,為在當(dāng)次最優(yōu)路線方案和當(dāng)前最優(yōu)路線之間取得平衡,亦即鼓勵隨后的螞蟻在當(dāng)次最優(yōu)路線和當(dāng)前最優(yōu)路線附近的解空間繼續(xù)搜索,改進的蟻群算法采用Ting 等[14]推薦的策略對當(dāng)次最優(yōu)和當(dāng)前最優(yōu)路線方案中的每條路線進行全局信息素更新,如式(23)、(24)所示.

    式中:EIW為當(dāng)次最差路線方案產(chǎn)生的碳排放.

    如此,通過若干次循環(huán)后即可得到所有以碳排放最小化為目標(biāo)的VRP 的最優(yōu)路線方案.考慮到由物流園區(qū)提供服務(wù)的配送中心數(shù)量較少而由每個配送中心提供服務(wù)的客戶數(shù)量相對較多,對得到的配送中心到客戶的VRP 最優(yōu)路線方案,本文運用3 種局部搜索算法對其再進行改進:cross_relocation 算法將從某個配送中心出發(fā)路線上的客戶逐一插入從其他配送中心出發(fā)的路線中,若得到的路線可行且碳排放量小于原路線,則用改進后的路線替代原路線;cross_swap 算法將從不同配送中心出發(fā)的兩條路線上的客戶兩兩交換并重新生成新的路線,若新路線可行且碳排放量較小,則替代原路線;最后,再對每一條路線運用3-opt 算法進行改進.至此,改進后的路線方案即可作為最終的最優(yōu)路線方案.

    3 算例測試

    由于目前還沒有可用于測試2E-LRP 碳排放的標(biāo)準(zhǔn)算例集,所以本文選取了Prodhon 標(biāo)準(zhǔn)算例集[15]中全部6 個最大規(guī)模的算例稍加修改,用以測試本文所提出的模型和算法.Prodhon 算例集中6 個最大規(guī)模的算例均包括200 個客戶和10 個候選配送中心.算例的修改包括去掉原算例中的配送中心選址費用、增加不同類型貨車的碳排放系數(shù),其中重型或中型貨車滿載和空載時的碳排放系數(shù)分別為2.392、1.638 kg/km,輕型或微型貨車滿載和空載時的碳排放系數(shù)分別為1.096、0.772 kg/km,其他數(shù)據(jù)保持不變.

    設(shè)計的TSHA 算法采用MATLAB 2018b 編程并嵌入Cplex V12.7 求解器,改進的蟻群算法控制參數(shù)包括:求解物流園區(qū)到被選用的配送中心VRP 問題時,循環(huán)次數(shù)設(shè)定為選用的配送中心個數(shù)的1/2,蟻群中螞蟻的數(shù)量與選用配送中心個數(shù)相同;求解配送中心到被分配客戶的VRP 問題時循環(huán)次數(shù)均設(shè)定為被分配客戶數(shù)量的1/2,蟻群中螞蟻的數(shù)量均等于被分配的客戶的數(shù)量;重要度系數(shù)α= 2,β= 1;給定常數(shù)q0= 0.5;揮發(fā)系數(shù)ρ= 0.2.針對每個算例,獨立運行程序20 次,并分別記下計算結(jié)果和時間.

    3.1 對TSHA 算法思想的驗證

    為驗證TSHA 的算法思想,本文先編寫了基于同樣算法思想的TSHA-Ⅱ算法用以求解Prodhon 算例集中6 個最大規(guī)模的原算例.TSHA-Ⅱ與TSHA略微不同的是:1) 用物流成本代替碳排放作為目標(biāo)函數(shù);2) 改進的蟻群算法中,路線構(gòu)建不考慮客戶的貨運需求量,即式 (20) 、(21) 中去掉啟發(fā)因子dj.TSHA-Ⅱ算法的求解結(jié)果與目前文獻中兩種最新算法的求解結(jié)果對比如表1 所示.表中:B為目前文獻所給出的最優(yōu)解;Cmin為基于某種算法得到的最優(yōu)解;Tmin為某種算法得到最優(yōu)解所需的時間;G為某個算法求得的最優(yōu)解與B之間的差異,即G=(Cmin-B)/B.

    表1 TSHA-II 算法與LNS-2E、Two-stage 算法求解2E-LRP 原算例的結(jié)果對比Tab.1 Comparison of the results obtained by LNS-2E, two-stage algorithm, and TSHA-II in solving original 2E-LRP instances

    通過對比可以看出:TSHA-Ⅱ算法求解大規(guī)模的2E-LRP 算例時,求解質(zhì)量稍遜于Breunig 等[16]提出的LNS-2E 算法,但明顯優(yōu)于Dai 等[17]提出的Twostage 算法;求解效率雖稍遜于后者,但大大優(yōu)于前者.這表明,TSHA-Ⅱ算法能夠在求解質(zhì)量與求解效率之間取得較好的平衡,可以在較短的時間內(nèi)得到較高質(zhì)量的解,TSHA 的算法思想可行有效.

    3.2 對考慮碳排放的2E-LRP 算例的求解

    算法思想通過驗證后,TSHA 再被用于求解考慮碳排放的2E-LRP 大規(guī)模算例,計算結(jié)果如表2所示.表中:Eavg為20 次運算得到的平均碳排放量;Emax、Emin分別為20 次運算得到的最大和最小碳排放量;Tavg為20 次運算的平均耗時.

    表2 TSHA 算法求解考慮碳排放的2E-LRP 算例的結(jié)果Tab.2 Results obtained by TSHA in solving 2E-LRP instances considering carbon emissions

    從表2 可以看出,TSHA 在求解考慮碳排放的2E-LRP 大規(guī)模算例時,求解結(jié)果的最大值、最小值和平均值以及計算時長的最小值和平均值之間相差很小,表現(xiàn)穩(wěn)定.可見,TSHA 可以作為求解考慮碳排放的2E-LRP 大規(guī)模算例的一種有效算法.

    4 結(jié) 論

    1) 相較于經(jīng)典的以物流成本最小化為目標(biāo)的2E-LRP 模型,本文基于一種簡單實用的以排放因子為主要參數(shù)的碳排放計算方法,建立了以碳排放最小化為目標(biāo)的2E-LRP 優(yōu)化模型.

    2) 對Prodhon 標(biāo)準(zhǔn)算例集中全部6 個最大規(guī)模算例的測試,表明TSHA-Ⅱ算法比現(xiàn)有算法更具兼顧求解結(jié)果和求解效率的優(yōu)勢,也證明了TSHA算法思想的可行性和有效性.

    3) TSHA 算法在求解考慮碳排放的2E-LRP 大規(guī)模算例時表現(xiàn)得高效而且穩(wěn)定,能夠滿足實際應(yīng)用場景的計算要求.

    4) 本文提出的模型和算法有助于物流企業(yè)通過優(yōu)化城市物流運作實現(xiàn)車輛的減排目標(biāo),服務(wù)國家“雙碳”戰(zhàn)略.

    猜你喜歡
    物流園區(qū)算例貨車
    貨車
    幼兒畫刊(2023年12期)2024-01-15 07:06:14
    貨車也便捷之ETC新時代!——看高速公路貨車ETC如何實現(xiàn)
    學(xué)與玩(2017年6期)2017-02-16 07:07:24
    物流園區(qū)出入口規(guī)劃設(shè)計及其優(yōu)化
    治超新規(guī)實施在即 深究貨車非法改裝亂象
    專用汽車(2016年9期)2016-03-01 04:16:52
    基于振蕩能量的低頻振蕩分析與振蕩源定位(二)振蕩源定位方法與算例
    物流園區(qū)的突圍之路
    中國儲運(2015年8期)2015-11-21 06:49:34
    互補問題算例分析
    一張圖帶你讀懂第四次全國物流園區(qū)(基地)調(diào)查報告 看看全國物流園區(qū)都有哪些“新”變化
    基于CYMDIST的配電網(wǎng)運行優(yōu)化技術(shù)及算例分析
    岛国毛片在线播放| 一级毛片我不卡| 亚洲av欧美aⅴ国产| 亚洲精品成人av观看孕妇| 亚洲精品国产色婷婷电影| 久久精品国产亚洲av高清一级| 韩国精品一区二区三区| 好男人视频免费观看在线| 欧美精品一区二区免费开放| 欧美黑人精品巨大| 欧美 亚洲 国产 日韩一| 精品午夜福利在线看| 巨乳人妻的诱惑在线观看| 男的添女的下面高潮视频| 精品国产一区二区三区久久久樱花| 国产成人一区二区在线| 久久久精品国产亚洲av高清涩受| 亚洲成人国产一区在线观看 | 多毛熟女@视频| 老鸭窝网址在线观看| 一级爰片在线观看| 国产老妇伦熟女老妇高清| 精品一区在线观看国产| 日韩av在线免费看完整版不卡| 国产淫语在线视频| 欧美日韩精品网址| 国产精品三级大全| 亚洲熟女精品中文字幕| 日韩,欧美,国产一区二区三区| 女的被弄到高潮叫床怎么办| 男男h啪啪无遮挡| 久久久久精品久久久久真实原创| 亚洲成人av在线免费| 人成视频在线观看免费观看| 最近中文字幕2019免费版| 亚洲精品国产区一区二| 成人影院久久| 免费观看人在逋| 国产野战对白在线观看| svipshipincom国产片| 国产精品一区二区在线观看99| 日日爽夜夜爽网站| 波野结衣二区三区在线| 亚洲国产精品国产精品| 老司机亚洲免费影院| 亚洲国产欧美在线一区| 亚洲精品自拍成人| 久久av网站| 亚洲国产欧美日韩在线播放| 国产精品一区二区在线不卡| 又大又黄又爽视频免费| 亚洲av日韩精品久久久久久密 | 波多野结衣一区麻豆| 亚洲精品一二三| 啦啦啦啦在线视频资源| 少妇猛男粗大的猛烈进出视频| 国产熟女午夜一区二区三区| 看免费成人av毛片| 久久人人97超碰香蕉20202| 精品国产国语对白av| 午夜激情久久久久久久| 亚洲精品久久久久久婷婷小说| 黄频高清免费视频| 国产成人一区二区在线| 纯流量卡能插随身wifi吗| 免费日韩欧美在线观看| 中文字幕最新亚洲高清| 操出白浆在线播放| 久久久久久久久久久久大奶| 成年人免费黄色播放视频| 亚洲欧美色中文字幕在线| 一本大道久久a久久精品| 国产成人啪精品午夜网站| av福利片在线| 满18在线观看网站| 久久久久精品人妻al黑| 国产亚洲精品第一综合不卡| 97精品久久久久久久久久精品| 欧美日韩成人在线一区二区| 久久精品久久久久久久性| 色婷婷av一区二区三区视频| 日本91视频免费播放| 日韩av免费高清视频| 搡老岳熟女国产| 欧美日韩视频高清一区二区三区二| 一个人免费看片子| 桃花免费在线播放| 伊人久久国产一区二区| 日韩精品有码人妻一区| 国产精品久久久久久久久免| 婷婷色综合大香蕉| 成人亚洲精品一区在线观看| 精品视频人人做人人爽| 黄色 视频免费看| 一本久久精品| 各种免费的搞黄视频| 亚洲精品日韩在线中文字幕| 亚洲,欧美,日韩| 新久久久久国产一级毛片| 免费看av在线观看网站| 午夜福利影视在线免费观看| 久久久久久久久免费视频了| 菩萨蛮人人尽说江南好唐韦庄| 欧美精品一区二区大全| 高清视频免费观看一区二区| 国产在线一区二区三区精| 无限看片的www在线观看| 精品一区在线观看国产| 我的亚洲天堂| 久久 成人 亚洲| 亚洲精品乱久久久久久| 亚洲七黄色美女视频| 熟女av电影| 亚洲伊人色综图| 亚洲精品第二区| 国精品久久久久久国模美| 亚洲综合色网址| 香蕉国产在线看| 欧美xxⅹ黑人| 欧美最新免费一区二区三区| 国产精品久久久人人做人人爽| 蜜桃国产av成人99| av电影中文网址| 精品第一国产精品| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久人妻精品电影 | 精品酒店卫生间| 中文字幕制服av| 曰老女人黄片| 国产成人免费无遮挡视频| 一级毛片 在线播放| 一区二区三区精品91| 久久99一区二区三区| 啦啦啦 在线观看视频| 国产亚洲一区二区精品| 精品一区二区三区四区五区乱码 | 欧美黑人精品巨大| 两个人看的免费小视频| 大陆偷拍与自拍| 无限看片的www在线观看| 日韩一本色道免费dvd| 在线 av 中文字幕| 午夜日本视频在线| 亚洲精品中文字幕在线视频| 日韩一卡2卡3卡4卡2021年| 美女午夜性视频免费| av一本久久久久| 男的添女的下面高潮视频| videosex国产| 少妇被粗大的猛进出69影院| 如何舔出高潮| 成年人午夜在线观看视频| 久久韩国三级中文字幕| 又大又黄又爽视频免费| 国产精品一国产av| 人体艺术视频欧美日本| 黄频高清免费视频| 乱人伦中国视频| 亚洲熟女毛片儿| 亚洲熟女毛片儿| 免费看av在线观看网站| 99热全是精品| 99热全是精品| 亚洲成人一二三区av| 久久精品久久久久久久性| 亚洲成人一二三区av| 欧美亚洲日本最大视频资源| 在线免费观看不下载黄p国产| 天天躁日日躁夜夜躁夜夜| 国产免费一区二区三区四区乱码| 99久久99久久久精品蜜桃| av卡一久久| 亚洲 欧美一区二区三区| 国产精品一区二区在线不卡| 中文字幕最新亚洲高清| 9色porny在线观看| 色视频在线一区二区三区| 亚洲精品av麻豆狂野| 亚洲,一卡二卡三卡| 毛片一级片免费看久久久久| 人妻人人澡人人爽人人| 夜夜骑夜夜射夜夜干| 午夜91福利影院| 日韩免费高清中文字幕av| 黄网站色视频无遮挡免费观看| 亚洲三区欧美一区| 啦啦啦啦在线视频资源| 美女视频免费永久观看网站| 亚洲精品中文字幕在线视频| 18禁动态无遮挡网站| 18禁观看日本| 熟女av电影| 老汉色av国产亚洲站长工具| 国产成人精品在线电影| 精品酒店卫生间| 肉色欧美久久久久久久蜜桃| 最近最新中文字幕大全免费视频 | av视频免费观看在线观看| 亚洲精品一二三| 国产在线视频一区二区| 人人妻人人澡人人爽人人夜夜| 精品国产一区二区三区久久久樱花| 免费高清在线观看日韩| 狠狠婷婷综合久久久久久88av| 国产成人精品无人区| 伊人亚洲综合成人网| 91国产中文字幕| 久久精品久久久久久噜噜老黄| 一区二区av电影网| 国产成人欧美在线观看 | 精品人妻一区二区三区麻豆| 久久精品国产a三级三级三级| 欧美精品高潮呻吟av久久| 国产精品香港三级国产av潘金莲 | 中文欧美无线码| 夫妻午夜视频| 女性被躁到高潮视频| 国产精品香港三级国产av潘金莲 | 深夜精品福利| 在线看a的网站| 色婷婷av一区二区三区视频| 久久ye,这里只有精品| 十八禁人妻一区二区| 日韩熟女老妇一区二区性免费视频| 嫩草影视91久久| 国产亚洲精品第一综合不卡| 美女中出高潮动态图| 少妇精品久久久久久久| avwww免费| 人妻 亚洲 视频| 一区二区日韩欧美中文字幕| 国产一区有黄有色的免费视频| 女人高潮潮喷娇喘18禁视频| 新久久久久国产一级毛片| 国产精品久久久av美女十八| 精品视频人人做人人爽| 蜜桃在线观看..| 国产毛片在线视频| 国产一区二区三区综合在线观看| 国产av码专区亚洲av| 欧美xxⅹ黑人| 免费高清在线观看视频在线观看| 日韩伦理黄色片| 久久久久视频综合| 亚洲免费av在线视频| 高清av免费在线| 日本av手机在线免费观看| 日日爽夜夜爽网站| 国产免费视频播放在线视频| 交换朋友夫妻互换小说| 免费人妻精品一区二区三区视频| 成人国产麻豆网| 男女边吃奶边做爰视频| 国产99久久九九免费精品| 欧美国产精品va在线观看不卡| 国产精品久久久久久精品电影小说| 青春草国产在线视频| 亚洲国产精品国产精品| 久热这里只有精品99| 国产成人精品福利久久| 人体艺术视频欧美日本| 一级,二级,三级黄色视频| 久久精品久久久久久噜噜老黄| 国产成人啪精品午夜网站| av卡一久久| 侵犯人妻中文字幕一二三四区| 欧美 日韩 精品 国产| 免费少妇av软件| 99国产精品免费福利视频| 男女免费视频国产| 亚洲天堂av无毛| av网站在线播放免费| 国精品久久久久久国模美| 欧美精品亚洲一区二区| 国产欧美日韩一区二区三区在线| 亚洲第一av免费看| 美女高潮到喷水免费观看| 母亲3免费完整高清在线观看| 香蕉国产在线看| 精品国产乱码久久久久久小说| 爱豆传媒免费全集在线观看| 人人妻人人添人人爽欧美一区卜| 国产精品熟女久久久久浪| 一区二区三区激情视频| 亚洲精品中文字幕在线视频| 新久久久久国产一级毛片| 宅男免费午夜| 国产福利在线免费观看视频| 婷婷色综合大香蕉| 精品第一国产精品| 97人妻天天添夜夜摸| 巨乳人妻的诱惑在线观看| 国产一区二区三区综合在线观看| 黄频高清免费视频| 99精品久久久久人妻精品| 亚洲av福利一区| 一区二区三区四区激情视频| 免费黄色在线免费观看| av在线老鸭窝| 久久鲁丝午夜福利片| 欧美xxⅹ黑人| 免费在线观看黄色视频的| 一级毛片电影观看| av电影中文网址| 飞空精品影院首页| avwww免费| 51午夜福利影视在线观看| 一区在线观看完整版| 男女边吃奶边做爰视频| 亚洲欧美一区二区三区久久| 久久精品久久久久久久性| 日韩 欧美 亚洲 中文字幕| 国产精品一区二区在线观看99| 18在线观看网站| 色婷婷av一区二区三区视频| 亚洲 欧美一区二区三区| 国产黄色视频一区二区在线观看| 高清在线视频一区二区三区| 亚洲精品美女久久久久99蜜臀 | 欧美日韩亚洲综合一区二区三区_| 成人手机av| 欧美中文综合在线视频| 99久久综合免费| 国产无遮挡羞羞视频在线观看| 国产又爽黄色视频| 久久久国产精品麻豆| 日本av免费视频播放| 激情五月婷婷亚洲| 欧美 亚洲 国产 日韩一| 欧美日韩亚洲高清精品| 国产伦理片在线播放av一区| 亚洲欧美一区二区三区黑人| 69精品国产乱码久久久| 国产熟女午夜一区二区三区| 日韩制服骚丝袜av| 色婷婷av一区二区三区视频| www.精华液| 我要看黄色一级片免费的| 久久久久久久国产电影| 成年动漫av网址| 肉色欧美久久久久久久蜜桃| 日韩不卡一区二区三区视频在线| 久久人人爽av亚洲精品天堂| 老司机影院成人| 97在线人人人人妻| 九草在线视频观看| 十八禁人妻一区二区| av网站在线播放免费| 又黄又粗又硬又大视频| 亚洲图色成人| 女性生殖器流出的白浆| 街头女战士在线观看网站| 亚洲欧洲日产国产| 在现免费观看毛片| 国产亚洲一区二区精品| 久热这里只有精品99| 久久免费观看电影| 国产精品人妻久久久影院| 欧美日韩视频高清一区二区三区二| 老司机在亚洲福利影院| 高清在线视频一区二区三区| 国产成人午夜福利电影在线观看| 香蕉丝袜av| 少妇人妻 视频| 啦啦啦中文免费视频观看日本| 大码成人一级视频| 丁香六月欧美| 韩国高清视频一区二区三区| 亚洲视频免费观看视频| 制服丝袜香蕉在线| 国产成人免费观看mmmm| 国精品久久久久久国模美| 国产亚洲一区二区精品| 欧美变态另类bdsm刘玥| 欧美在线一区亚洲| 女人精品久久久久毛片| 毛片一级片免费看久久久久| 蜜桃国产av成人99| 18禁国产床啪视频网站| 日韩av在线免费看完整版不卡| 亚洲综合色网址| 免费日韩欧美在线观看| 97在线人人人人妻| 成年人免费黄色播放视频| 亚洲人成电影观看| 欧美精品av麻豆av| 欧美日韩av久久| 国产又爽黄色视频| 在线看a的网站| av网站在线播放免费| 国产精品久久久久成人av| 亚洲国产成人一精品久久久| 国产精品.久久久| 男人操女人黄网站| 免费观看性生交大片5| 妹子高潮喷水视频| 欧美激情 高清一区二区三区| a级毛片在线看网站| 国产亚洲欧美精品永久| 一本大道久久a久久精品| av在线老鸭窝| 久久久国产精品麻豆| 在线观看免费午夜福利视频| 啦啦啦在线免费观看视频4| 在线亚洲精品国产二区图片欧美| 尾随美女入室| 麻豆av在线久日| 国产精品人妻久久久影院| 一本久久精品| 97在线人人人人妻| 国产成人精品在线电影| 国产成人91sexporn| 国产黄色视频一区二区在线观看| 国产亚洲av片在线观看秒播厂| 久久久久精品人妻al黑| 一区在线观看完整版| 人人妻人人添人人爽欧美一区卜| 777久久人妻少妇嫩草av网站| 日本欧美视频一区| 最近中文字幕2019免费版| 最近最新中文字幕免费大全7| 啦啦啦视频在线资源免费观看| 国产精品久久久久久精品古装| 熟女av电影| 日本色播在线视频| 9色porny在线观看| 丝袜在线中文字幕| 久久精品人人爽人人爽视色| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区乱码不卡18| 国产日韩欧美在线精品| 电影成人av| 色视频在线一区二区三区| 自线自在国产av| 大话2 男鬼变身卡| √禁漫天堂资源中文www| 婷婷成人精品国产| 免费高清在线观看日韩| www.精华液| 亚洲精品在线美女| 久久女婷五月综合色啪小说| 2018国产大陆天天弄谢| 久久久精品区二区三区| 久久精品人人爽人人爽视色| 中国三级夫妇交换| 最近2019中文字幕mv第一页| 欧美黄色片欧美黄色片| 成人国产麻豆网| 国产老妇伦熟女老妇高清| 欧美黑人欧美精品刺激| 亚洲国产av新网站| av女优亚洲男人天堂| 国产精品二区激情视频| 亚洲免费av在线视频| 97精品久久久久久久久久精品| av有码第一页| 大片电影免费在线观看免费| 天堂俺去俺来也www色官网| 亚洲成色77777| 熟女少妇亚洲综合色aaa.| 欧美亚洲日本最大视频资源| 欧美激情高清一区二区三区 | 成年av动漫网址| 中文字幕另类日韩欧美亚洲嫩草| 国产片内射在线| 青春草国产在线视频| 欧美最新免费一区二区三区| 国产精品三级大全| 在线观看www视频免费| 午夜福利在线免费观看网站| 夫妻午夜视频| 韩国精品一区二区三区| 最近最新中文字幕大全免费视频 | 欧美在线一区亚洲| 欧美精品一区二区免费开放| 毛片一级片免费看久久久久| 国产精品香港三级国产av潘金莲 | 亚洲精品自拍成人| 91精品国产国语对白视频| 日韩 欧美 亚洲 中文字幕| 如日韩欧美国产精品一区二区三区| 日韩电影二区| 国产精品偷伦视频观看了| netflix在线观看网站| 赤兔流量卡办理| 精品一区二区免费观看| 国产亚洲av高清不卡| 99精国产麻豆久久婷婷| 国产视频首页在线观看| 男人舔女人的私密视频| 国产精品 国内视频| 91精品伊人久久大香线蕉| 亚洲国产毛片av蜜桃av| 亚洲美女黄色视频免费看| 一本大道久久a久久精品| 欧美黑人欧美精品刺激| 亚洲一级一片aⅴ在线观看| 一边摸一边做爽爽视频免费| 天天添夜夜摸| 午夜久久久在线观看| 久热爱精品视频在线9| 秋霞伦理黄片| 女性生殖器流出的白浆| 欧美日韩成人在线一区二区| 午夜激情久久久久久久| 毛片一级片免费看久久久久| 麻豆精品久久久久久蜜桃| 欧美av亚洲av综合av国产av | 黄片播放在线免费| 国产深夜福利视频在线观看| 成人毛片60女人毛片免费| 观看美女的网站| 精品一品国产午夜福利视频| 国产欧美日韩综合在线一区二区| 国产精品免费视频内射| 亚洲av电影在线进入| 国产片内射在线| 天天躁夜夜躁狠狠躁躁| 午夜日本视频在线| 亚洲国产最新在线播放| 久久天堂一区二区三区四区| 精品国产乱码久久久久久小说| 丝袜脚勾引网站| 日韩熟女老妇一区二区性免费视频| 欧美中文综合在线视频| 日韩制服丝袜自拍偷拍| 十八禁人妻一区二区| 日本欧美国产在线视频| 99国产精品免费福利视频| av又黄又爽大尺度在线免费看| 婷婷成人精品国产| 狠狠婷婷综合久久久久久88av| 一区二区三区四区激情视频| 另类亚洲欧美激情| 最近最新中文字幕免费大全7| 亚洲精品成人av观看孕妇| 51午夜福利影视在线观看| 中文字幕高清在线视频| 亚洲精品美女久久av网站| 青春草亚洲视频在线观看| 精品一区二区三卡| 少妇人妻 视频| 热re99久久精品国产66热6| 高清黄色对白视频在线免费看| 老汉色av国产亚洲站长工具| 哪个播放器可以免费观看大片| 狂野欧美激情性bbbbbb| 国产av码专区亚洲av| 日本av手机在线免费观看| 久久婷婷青草| 女人精品久久久久毛片| 悠悠久久av| 丝袜人妻中文字幕| 老司机亚洲免费影院| 亚洲国产日韩一区二区| 成年女人毛片免费观看观看9 | 美女视频免费永久观看网站| 亚洲精品aⅴ在线观看| 9191精品国产免费久久| 亚洲,欧美精品.| 自拍欧美九色日韩亚洲蝌蚪91| 丝瓜视频免费看黄片| 国产精品偷伦视频观看了| 老司机在亚洲福利影院| 制服丝袜香蕉在线| 咕卡用的链子| 黄色一级大片看看| 丰满迷人的少妇在线观看| 国产av一区二区精品久久| 色吧在线观看| 在线亚洲精品国产二区图片欧美| 香蕉国产在线看| 一边摸一边做爽爽视频免费| 久久久久久久久久久久大奶| 日本欧美国产在线视频| 人体艺术视频欧美日本| 国产成人a∨麻豆精品| 国产探花极品一区二区| 精品一区二区三区av网在线观看 | 亚洲综合精品二区| 老司机影院毛片| 日日爽夜夜爽网站| 不卡视频在线观看欧美| 搡老乐熟女国产| 韩国精品一区二区三区| 男女下面插进去视频免费观看| 99九九在线精品视频| 1024香蕉在线观看| 亚洲欧美成人综合另类久久久| 国产精品一二三区在线看| 这个男人来自地球电影免费观看 | 九草在线视频观看| 国产99久久九九免费精品| 国产一区二区三区综合在线观看| 天美传媒精品一区二区| 久久久久精品人妻al黑| 久久99热这里只频精品6学生| 亚洲在久久综合| 久久人人爽av亚洲精品天堂| 人人妻人人爽人人添夜夜欢视频| 成人18禁高潮啪啪吃奶动态图| 18禁国产床啪视频网站| 欧美久久黑人一区二区| 国产探花极品一区二区| 久久97久久精品| 亚洲av电影在线进入| 亚洲国产最新在线播放| 香蕉丝袜av| 免费日韩欧美在线观看| 午夜福利网站1000一区二区三区| 国产精品久久久久久精品电影小说| 亚洲av欧美aⅴ国产| 狂野欧美激情性bbbbbb| 欧美亚洲 丝袜 人妻 在线| 在线观看www视频免费|