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

    基于O2O的大型零售企業(yè)城市配送網(wǎng)絡(luò)優(yōu)化研究

    2017-10-13 03:25:33趙泉午趙軍平
    中國管理科學(xué) 2017年9期
    關(guān)鍵詞:分配距離物流

    趙泉午,趙軍平,林 婭

    (1.重慶大學(xué)經(jīng)濟(jì)與工商管理學(xué)院,重慶 400030;2.重慶大學(xué)現(xiàn)代物流重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 400030)

    基于O2O的大型零售企業(yè)城市配送網(wǎng)絡(luò)優(yōu)化研究

    趙泉午1,2,趙軍平1,2,林 婭2

    (1.重慶大學(xué)經(jīng)濟(jì)與工商管理學(xué)院,重慶 400030;2.重慶大學(xué)現(xiàn)代物流重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 400030)

    本文研究了O2O轉(zhuǎn)型背景下大型零售企業(yè)城市配送網(wǎng)絡(luò)優(yōu)化面臨的中轉(zhuǎn)中心選址及末端需求點(diǎn)分配問題,構(gòu)建了考慮配送中心到末端需求點(diǎn)近似配送距離的中轉(zhuǎn)中心選址及末端需求點(diǎn)分配聯(lián)合優(yōu)化模型,設(shè)計(jì)了集成遺傳算法和禁忌搜索算法的混合算法求解模型,通過混合算法與CPLEX的對比證明了本文算法的有效性。以蘇寧重慶主城核心區(qū)域的城市配送網(wǎng)絡(luò)優(yōu)化為實(shí)例進(jìn)行測算,給出了蘇寧重慶主城核心區(qū)域中轉(zhuǎn)中心的選址及需求點(diǎn)分配的聯(lián)合優(yōu)化方案。聯(lián)合優(yōu)化方案表明城市物流“最后一公里”末端配送成本占城市物流總成本的比例超過60%;車型是影響城市物流總成本的關(guān)鍵因素,選擇與中轉(zhuǎn)中心容量接近的車型能夠顯著降低城市物流總成本,減少運(yùn)輸距離;合理的中轉(zhuǎn)中心選擇能提高短距離配送比例和減少配送總距離。

    O2O轉(zhuǎn)型;選址-分配;“最后一公里”;遺傳算法;禁忌搜索算法

    1 引言

    布局O2O模式是國內(nèi)大型零售企業(yè)應(yīng)對電商對沖擊的重要途徑。2014年國內(nèi)連鎖百強(qiáng)有75家開展網(wǎng)絡(luò)零售,其中23%的門店提供自提服務(wù)。門店是實(shí)施O2O的關(guān)鍵環(huán)節(jié),目前國內(nèi)大型零售企業(yè)的O2O轉(zhuǎn)型模式分為三種。第一種是門店僅提供商品展示及消費(fèi)體驗(yàn)服務(wù),不涉及任何線上訂單的物流功能。第二種是門店提供線上訂單的物流中轉(zhuǎn)服務(wù),以少數(shù)門店作為中轉(zhuǎn)中心,為周邊末端需求點(diǎn)提供配送服務(wù)。第三種是門店在前兩種模式的基礎(chǔ)上,提供客戶自提、商品包裝等服務(wù),實(shí)現(xiàn)了線上線下庫存及服務(wù)的共享融合。目前國內(nèi)大型零售企業(yè)(蘇寧、國美等)O2O布局以前兩種模式為主。城市內(nèi)部核心區(qū)域需求密度大,需求量集中,中轉(zhuǎn)中心選址關(guān)系到城市物流系統(tǒng)運(yùn)作效率和成本,各企業(yè)紛紛優(yōu)化布局城市核心區(qū)域的中轉(zhuǎn)中心,完善“最后一公里”配送。

    中轉(zhuǎn)中心選址問題屬于物流設(shè)施選址問題,大量文獻(xiàn)進(jìn)行了研究。Aikens等[1]研究了現(xiàn)有多種選址模型,認(rèn)為模型的選擇應(yīng)根據(jù)問題特征而定。Shen等[2]研究了包含一個(gè)供應(yīng)商,多個(gè)需求不確定零售商的選址及庫存聯(lián)合優(yōu)化模型,建立非線性規(guī)劃模型進(jìn)行求解。Tsao等[3]綜合考慮區(qū)域配送中心選址、需求點(diǎn)分配及區(qū)域配送中心庫存策略使物流總成本最小,論文采用非線性規(guī)劃求解。Taniguchi等[4]應(yīng)用雙層規(guī)劃研究了公共物流站點(diǎn)選址問題,采用遺傳算法得出最優(yōu)解。Wang Yong等[5]采用遺傳算法和粒子群混合算法研究了兩層物流網(wǎng)絡(luò)設(shè)施選址及分配問題,并以中國貴陽為例對算法進(jìn)行檢驗(yàn),其設(shè)計(jì)的算法適用于優(yōu)化小規(guī)模的雙層物流系統(tǒng)。Diabat等[6]研究了單個(gè)產(chǎn)品多層次庫存選址聯(lián)合優(yōu)化問題,并采用拉格朗日松弛多階段算法進(jìn)行求解。Amir等[7]研究了價(jià)格影響需求和庫存容量對物流設(shè)施選址的影響,通過數(shù)據(jù)模擬得出所建模型適用于優(yōu)化價(jià)格、選址和庫存聯(lián)合優(yōu)化。Jouzdani等[8]通過構(gòu)建混合整數(shù)規(guī)劃模型研究了交通擁堵及需求不確定下奶制品供應(yīng)鏈規(guī)劃問題。唐金環(huán)[9-12]研究了碳排放對企業(yè)物流設(shè)施選址的影響,企業(yè)可通過設(shè)施選址、車輛選擇、路徑設(shè)計(jì)等方式以較小代價(jià)降低碳排放。O2O物流設(shè)施選址與B2C物流設(shè)施選址具有相似特征,針對B2C物流設(shè)施選址問題,劉必爭等[13]認(rèn)為B2C車輛配送業(yè)務(wù)需考慮車輛的回訪特性,并建立帶軟時(shí)間窗的選址-路徑集成模型,結(jié)合遺傳算法和退火算法求解。周翔等[14]考慮顧客滿意度最大和運(yùn)輸成本最低兩個(gè)目標(biāo),對B2C模式下城市配送物流網(wǎng)絡(luò)中末端節(jié)點(diǎn)數(shù)量、配送中心和末端節(jié)點(diǎn)的選址進(jìn)行研究,構(gòu)建了兩階段優(yōu)化模型并求解。蔣忠中等[15]建立選址-分配模型,設(shè)計(jì)基于嵌入表上作業(yè)法的遺傳算法求解。Xu Bin等[16]立足于提升物流系統(tǒng)效率和降低運(yùn)輸成本研究了B2C兩級物流分撥體系,采用Lingo軟件求解證明模型的有效性。Torabi等[17]研究了電商物流庫存設(shè)施選址和庫存轉(zhuǎn)運(yùn)問題,建立了混合整數(shù)規(guī)劃模型,通過數(shù)據(jù)仿真說明了庫存轉(zhuǎn)運(yùn)對于電商企業(yè)的重要性。Yu等[18]研究了引入第三方物流時(shí)考慮車輛不返回的設(shè)施選址及路徑問題,采用模擬退火算法進(jìn)行求解,通過算例證明算法有效性。文獻(xiàn)[19-22]研究了模糊需求下物流設(shè)施選址問題,集成優(yōu)化設(shè)施選址、分配及路徑規(guī)劃。

    現(xiàn)有文獻(xiàn)側(cè)重優(yōu)化算法設(shè)計(jì),以隨機(jī)生成中小規(guī)模算例驗(yàn)證算法有效性,缺乏明確應(yīng)用背景,并且末端需求點(diǎn)規(guī)模設(shè)置鮮有超過300個(gè)。本文以蘇寧等國內(nèi)大型零售企業(yè)O2O轉(zhuǎn)型實(shí)踐為背景,綜合考慮中轉(zhuǎn)中心運(yùn)營成本、車輛往返配送成本和車型選擇等因素,構(gòu)建整數(shù)規(guī)劃模型對具有容量限制的城市配送中轉(zhuǎn)中心選址、末端需求點(diǎn)分配和末端配送問題進(jìn)行研究。本文的創(chuàng)新處體現(xiàn)在如下幾點(diǎn):一是基于蘇寧等國內(nèi)大型零售企業(yè)O2O轉(zhuǎn)型實(shí)踐為背景研究城市配送中轉(zhuǎn)中心選址及末端需求點(diǎn)分配問題。二是針對問題特性及模型特點(diǎn),論文綜合了遺傳算法和禁忌搜索算法的優(yōu)點(diǎn),基于啟發(fā)式規(guī)則設(shè)計(jì)了集成混合算法,并通過混合算法與CPLEX軟件求解結(jié)果的對比,驗(yàn)證本文混合算法的有效性。三是采用Daganzo[24]提出的距離估算方法測算中轉(zhuǎn)中心到末端需求點(diǎn)的配送距離,該測算方法使得配送距離更加接近車輛日常配送距離。四是本文以蘇寧重慶主城核心區(qū)域?yàn)閷?shí)例進(jìn)行測算,給出了不同情景下選址及末端需求點(diǎn)分配的具體方案,所得結(jié)論契合實(shí)際,具有較高的決策參考價(jià)值。

    2 問題及模型

    2.1問題描述

    本文研究國內(nèi)大型零售企業(yè)實(shí)施O2O轉(zhuǎn)型面臨的城市配送網(wǎng)絡(luò)優(yōu)化設(shè)計(jì)問題,也就是從城市配送中心到消費(fèi)者過程中面臨的中轉(zhuǎn)中心選址和末端需求點(diǎn)分配問題,中轉(zhuǎn)中心的備選地址為城市區(qū)域內(nèi)全部零售門店。城市配送中心承擔(dān)庫存、包裝、分揀、運(yùn)輸甚至流通加工等功能。中轉(zhuǎn)中心承擔(dān)訂單的中轉(zhuǎn)分撥功能。末端需求點(diǎn)包括自提點(diǎn)、公共快遞取送點(diǎn)、自提柜等形態(tài)。城市配送中心到中轉(zhuǎn)中心一般采用輕型卡車運(yùn)輸,中轉(zhuǎn)中心到末端需求點(diǎn)采用三輪車、摩托車、面包車及電動車等車型進(jìn)行配送。末端需求點(diǎn)同時(shí)回收客戶的退換貨,將退換貨送至中轉(zhuǎn)中心,由中轉(zhuǎn)中心轉(zhuǎn)送至城市配送中心。城市配送物流總成本包括從城市配送中心到中轉(zhuǎn)中心的運(yùn)輸成本、中轉(zhuǎn)中心運(yùn)營成本和中轉(zhuǎn)中心到末端需求點(diǎn)的配送成本。

    2.2變量及參數(shù)定義

    為方便刻畫大型零售企業(yè)城市配送網(wǎng)絡(luò)優(yōu)化中面臨的中轉(zhuǎn)中心選址以及末端需求點(diǎn)分配問題,定義如下變量及參數(shù)。

    I:備選中轉(zhuǎn)中心集合,I={1,2,…,m};

    J:末端需求點(diǎn)集合,J={1,2,…,n};

    Oi:備選中轉(zhuǎn)中心i每天運(yùn)營固定成本;

    F(·):中轉(zhuǎn)中心i每天包裹處理成本,包括裝卸、分撥等成本,與包裹處理量有關(guān);

    f1,f2:配送中心到中轉(zhuǎn)中心、中轉(zhuǎn)中心到需求點(diǎn)運(yùn)輸(配送)車輛固定成本;

    c1:配送中心到中轉(zhuǎn)中心貨車單位距離運(yùn)輸成本;

    di:配送中心與中轉(zhuǎn)中心i之間的距離;

    Mi:備選中轉(zhuǎn)中心i日最大中轉(zhuǎn)量;

    vc1,vc2:配送中心到中轉(zhuǎn)中心、中轉(zhuǎn)中心到末端需求點(diǎn)車輛容量;

    c2,c3: 中轉(zhuǎn)中心到末端需求點(diǎn)的配送車輛每天人工成本、配送車輛單位運(yùn)輸成本;

    rj,qj:末端需求點(diǎn)j平均每天配送量、退貨量;

    Li:中轉(zhuǎn)中心i到末端需求點(diǎn)的配送車輛總行駛距離;

    zij:當(dāng)j客戶需求點(diǎn)由中轉(zhuǎn)中心i服務(wù)時(shí)為1,否則為0;

    yi:是否選中備選地點(diǎn)i作為中轉(zhuǎn)中心,選中為1,否則為0。

    2.3聯(lián)合優(yōu)化模型

    為更好構(gòu)建中轉(zhuǎn)中心選址及末端需求點(diǎn)分配的聯(lián)合優(yōu)化模型,本文做如下三個(gè)重要假定。

    假定1:以末端需求點(diǎn)的日平均配送量為基準(zhǔn)進(jìn)行優(yōu)化

    末端需求點(diǎn)的配送量是變化的,本文研究側(cè)重戰(zhàn)略層面的中轉(zhuǎn)中心選址問題,因此以末端需求點(diǎn)的日平均配送量為基準(zhǔn)進(jìn)行優(yōu)化。每天具體配送量的變化可以在每天運(yùn)作層面的路徑及調(diào)度優(yōu)化中實(shí)現(xiàn)。

    假定2:中轉(zhuǎn)中心到多個(gè)末端需求點(diǎn)配送距離采用近似算法

    為了使配送距離近似測算更接近日常狀況,采用連續(xù)逼近法測算配送車輛的行駛距離。Eilon[23]通過仿真得出當(dāng)末端點(diǎn)均勻分布在正方形區(qū)域,中轉(zhuǎn)中心i到所服務(wù)末端需求點(diǎn)的配送距離Li與所服務(wù)末端需求點(diǎn)個(gè)數(shù)N之間的關(guān)系可以近似表示為:

    假定3:中轉(zhuǎn)中心運(yùn)營成本存在規(guī)模效應(yīng),采用階梯函數(shù)刻畫運(yùn)營成本

    F(x)表示中轉(zhuǎn)中心的包裹處理費(fèi)用,通常隨著中轉(zhuǎn)中心中轉(zhuǎn)量的增長。由于存在規(guī)模效應(yīng),運(yùn)營成本不能認(rèn)為是簡單的線性函數(shù)[20]。成本函數(shù)可表示為:

    x表示i中轉(zhuǎn)中心包裹處理量,Ni和Mi是中轉(zhuǎn)中心的臨界參數(shù)、最大中轉(zhuǎn)量,Ej0為系數(shù),θ中轉(zhuǎn)中心規(guī)模效應(yīng)系數(shù),θ∈(0,1)。

    基于以上三個(gè)重要假定,構(gòu)建如下聯(lián)合優(yōu)化模型:

    (1)

    s.t.

    ∑izij=1,?j∈J

    (2)

    ∑jzijrj≤Mi,?i∈I

    (3)

    rj>qj,?j∈J

    (4)

    zij≤yi,?i∈I,?j∈J

    (5)

    ∑jzij≥yi,?i∈I,?j∈J

    (6)

    zij,yi∈{0,1},?i∈I,?j∈J

    (7)

    式(1)為最小化物流總成本,其中包括從配送中心到中轉(zhuǎn)中心運(yùn)輸成本,中轉(zhuǎn)中心運(yùn)營成本和包裹處理成本,中轉(zhuǎn)中心到末端需求點(diǎn)的配送成本。約束條件(2)表明末端需求點(diǎn)有且只有一個(gè)中轉(zhuǎn)中心為其提供服務(wù),約束條件(3)表明中轉(zhuǎn)中心日中轉(zhuǎn)量不得大于中轉(zhuǎn)中心容量,約束條件(4)表明末端需求點(diǎn)退換量小于需求量,保證車輛動態(tài)負(fù)載不超過其容量,約束條件(5)(6)(7)為決策變量之間的關(guān)系及變量屬性。

    3 模型求解

    對于模型的求解,當(dāng)前學(xué)者采用智能算法與啟發(fā)式算法相結(jié)合的方式。國外文獻(xiàn)方面,Tsao等[3]采用連續(xù)逼近算法求解,Wang Yong等[5]采用擴(kuò)展粒子群算法與遺傳算法相結(jié)合的混合算法,Diabat、Amir等[6-7,12]分別采用基于拉格朗日松弛的啟發(fā)式算法和遺傳算法,Torabi等[17]采用Bender分解法等。國內(nèi)方面,唐金環(huán)等[9]采用基于遺傳算法的細(xì)菌覓食和粒子群理論(BFA-PSO),楊珺等[10]采用禁忌搜索算法,周翔等[14]采用基于最小生成樹聚類算法,關(guān)菲等[19]采用多目標(biāo)粒子群算法求解模型。本文提出了基于末端需求點(diǎn)配送量及配送距離的相對大小的啟發(fā)式規(guī)則分配末端需求點(diǎn),綜合利用遺傳算法、禁忌搜索算法和啟發(fā)式算法的優(yōu)點(diǎn),設(shè)計(jì)混合優(yōu)化算法對模型進(jìn)行求解。

    3.1遺傳算法和禁忌搜索混合優(yōu)化算法

    本文使用遺傳算法生成中轉(zhuǎn)中心開放方案,充分利用遺傳算法搜索范圍廣的特點(diǎn),由于禁忌搜索算法對初始解有較強(qiáng)的依賴性,基于啟發(fā)式規(guī)則生成末端需求點(diǎn)初始分配方案,采用禁忌算法優(yōu)化末端需求點(diǎn)的分配方案?;旌纤惴ǖ木唧w步驟如下:

    步驟1隨機(jī)產(chǎn)生一定數(shù)量的選址方案作為初始種群,計(jì)算當(dāng)前種群中全部開放中轉(zhuǎn)中心總?cè)萘俊H艨側(cè)萘啃∮谀┒诵枨簏c(diǎn)的總需求,增選任意尚未開放的備選中轉(zhuǎn)中心,直到總?cè)萘看笥诳傂枨蟆?/p>

    步驟2對種群中個(gè)體執(zhí)行以下操作:

    步驟2.1按3.2所示方法生成末端需求點(diǎn)初始分配方案。

    步驟 2.2按3.3所示方法生成優(yōu)化分配方案。

    步驟 2.3若達(dá)到終止條件,轉(zhuǎn)步驟3,否則轉(zhuǎn)步驟2.2。

    步驟3記錄當(dāng)前最優(yōu)解和對應(yīng)目標(biāo)值。

    步驟4采用適應(yīng)度函數(shù)計(jì)算當(dāng)前種群個(gè)體的適應(yīng)度。

    步驟5對種群進(jìn)行選擇、交叉和變異操作,并將新的最優(yōu)解加入下一代種群中,同時(shí)更新種群代數(shù)。

    步驟6當(dāng)?shù)螖?shù)大于算法設(shè)定的最大代數(shù)時(shí),輸出優(yōu)化結(jié)果。否則轉(zhuǎn)步驟2。

    3.2基于啟發(fā)式規(guī)則的初始分配算法

    采用末端需求點(diǎn)配送量及配送距離比值由大到小的啟發(fā)式規(guī)則對末端需求點(diǎn)進(jìn)行分配:

    步驟1測算開放的中轉(zhuǎn)中心與末端需求點(diǎn)之間的距離dij。

    步驟3取矩陣Q中kij的最大值ki′j′,ki′j′表示末端需求點(diǎn)j′預(yù)分配中轉(zhuǎn)中心i′,檢查i′中轉(zhuǎn)中心剩余容量能否滿足末端需求點(diǎn)j.當(dāng)滿足時(shí),將末端需求點(diǎn)j′分配給中轉(zhuǎn)中心i′,同時(shí)刪除該末端需求點(diǎn)j′所在列;當(dāng)不能滿足時(shí),ki′j′點(diǎn)置零。

    步驟4重復(fù)執(zhí)行步驟3直到將全部末端需求點(diǎn)完全分配給中轉(zhuǎn)中心,生成初始分配方案。

    3.3基于禁忌搜索算法的優(yōu)化分配算法

    定義“移動”、“交換”兩個(gè)操作。以A、B為例,“移動”指將中轉(zhuǎn)中心A服務(wù)的一個(gè)或幾個(gè)末端需求點(diǎn)改由B中轉(zhuǎn)中心服務(wù);“交換”指屬于A服務(wù)的若干個(gè)末端需求與屬于B服務(wù)的若干末端需求點(diǎn)進(jìn)行交換。

    應(yīng)用禁忌搜索算法對分配方案進(jìn)行優(yōu)化,具體步驟如下:

    步驟1生成選中中轉(zhuǎn)中心與末端需求點(diǎn)距離表,1A表示中轉(zhuǎn)中心A服務(wù)需求點(diǎn)1,如表1所示:

    步驟2“移動”操作。隨機(jī)選擇兩個(gè)選中的中轉(zhuǎn)中心A和B,計(jì)算A中轉(zhuǎn)中心服務(wù)的末端需求點(diǎn)距離B中轉(zhuǎn)中心的距離,將距離最小的需求點(diǎn)改由B中轉(zhuǎn)中心服務(wù),判斷是否滿足B中轉(zhuǎn)中心容量約束;若不滿足容量約束,放棄當(dāng)前操作;若滿足執(zhí)行“移動”操作。

    步驟3“交換”操作。隨機(jī)選擇兩個(gè)選中的中轉(zhuǎn)中心A和B,分別計(jì)算A中轉(zhuǎn)中心服務(wù)的末端需求點(diǎn)距離B中轉(zhuǎn)中心的最近距離,以及B中轉(zhuǎn)中心服務(wù)的末端需求點(diǎn)距離A中轉(zhuǎn)中心的最近距離。交換中轉(zhuǎn)中心A和B中距離對方最近需求點(diǎn),同時(shí)判斷末端總需求是否超過容量。若未超過,執(zhí)行“交換”操作;否則放棄“交換”操作。

    步驟4當(dāng)候選解滿足藐視準(zhǔn)則時(shí),用滿足藐視準(zhǔn)則的最佳狀態(tài)替代當(dāng)前狀態(tài)成為新的當(dāng)前解,并用與原狀態(tài)對應(yīng)的禁忌對象替換最早進(jìn)入禁忌表的禁忌對象,并轉(zhuǎn)步驟6;否則轉(zhuǎn)步驟5。

    步驟5判斷候選解對應(yīng)的各對象的禁忌屬性,選擇候選解集中非禁忌對象對應(yīng)的最佳狀態(tài)為新的當(dāng)前解,同時(shí)用與之對應(yīng)的禁忌對象替換最早進(jìn)入禁忌表的禁忌對象。

    步驟6判斷終止條件,是則停止,否則轉(zhuǎn)步驟2。

    4 算法有效性分析

    CPLEX是IBM開發(fā)的求解規(guī)劃問題的商業(yè)軟件之一。能最少用戶干預(yù)、精確地解決問題,已有學(xué)者在研究解決較大規(guī)模物流與供應(yīng)鏈問題時(shí)證明CPLEX在解決規(guī)模較小問題時(shí)具有好的求解精確度,但在求解問題規(guī)模較大時(shí),CPLEX效率不高[27]。本文采用CPLEX求解結(jié)果作為驗(yàn)證混合算法求解有效性的基準(zhǔn)。

    構(gòu)建一組算例驗(yàn)證上述模型和算法的可行性和有效性。算例的參數(shù)采用5.1實(shí)際測算參數(shù),各項(xiàng)數(shù)據(jù)在對應(yīng)區(qū)間范圍內(nèi)隨機(jī)生成。表1中規(guī)模4*10指4個(gè)備選中轉(zhuǎn)中心和10個(gè)末端需求點(diǎn),其他情形以此類推。本文采用Matlab R14編程語言實(shí)現(xiàn)上述混合算法,運(yùn)行平臺為:

    AMD,A8-7650K,Radeon,R7,10Compute Cores 4C+6G@3.3GHZ,Win8 PC。

    針對問題特性,混合算法在求解模型過程中采用啟發(fā)式規(guī)則進(jìn)行人工干預(yù),生成一個(gè)較優(yōu)質(zhì)初始解,極大提高了算法后續(xù)搜索效率。模型求解配送距離方式較為復(fù)雜, CPLEX在求解時(shí)間顯著大于混合算法。混合算法執(zhí)行結(jié)果與CPLEX模型求解結(jié)果相對比,兩者效果相當(dāng),混合算法效率高于CPLEX,也說明了CPLEX在針對小規(guī)模問題時(shí)表現(xiàn)出更優(yōu)異的性能,當(dāng)問題規(guī)模較大情況較為復(fù)雜時(shí)混合算法表現(xiàn)更加優(yōu)異。

    表1 CPLEX求解與混合算法性能比較

    5 蘇寧重慶實(shí)例分析

    5.1模型測算

    本文以蘇寧重慶主城核心區(qū)域的城市配送網(wǎng)絡(luò)為實(shí)例,各項(xiàng)數(shù)據(jù)均來源于蘇寧重慶實(shí)地調(diào)研。蘇寧在重慶擁有1個(gè)城市配送中心,承擔(dān)線上訂單的分撥轉(zhuǎn)運(yùn)和線下商品的庫存、分揀、包裝和城市配送等職能。主城核心區(qū)域15家門店(表2),為實(shí)現(xiàn)線上線下城市配送的融合,在15家門店中選擇若干家作為線上包裹中轉(zhuǎn)中心。選中的門店作為實(shí)體門店展示銷售商品的同時(shí),承擔(dān)周邊一定范圍內(nèi)線上包裹的中轉(zhuǎn)和末端需求點(diǎn)配送服務(wù)。本文將重慶主城核心區(qū)域15家門店作為備選中轉(zhuǎn)中心,末端需求點(diǎn)劃分根據(jù)重慶主城核心區(qū)域的人口分布、消費(fèi)水平及交通狀況等,將重慶主城核心區(qū)域劃分為419個(gè)末端需求點(diǎn),各個(gè)末端需求點(diǎn)的線上包裹配送需求量依據(jù)實(shí)地調(diào)研數(shù)據(jù)和經(jīng)驗(yàn)估算。

    表2 蘇寧重慶主城區(qū)核心區(qū)域門店信息

    從配送中心到中轉(zhuǎn)中心采用直送方式,以福田歐馬可1.49T廂式輕卡為例,車輛購買費(fèi)用10.23萬元,車輛購置稅8744元,折舊年限5年,保險(xiǎn)稅5203元/年,駕駛員工資6000元/月,以此計(jì)算車輛平均行駛固定成本為272元/天,車輛行駛變動費(fèi)用3元/公里(包括過橋過路費(fèi)、燃油費(fèi)、輪胎損耗費(fèi)等)。額定載重1.49T運(yùn)輸車包裹輸送量約1000個(gè)/車·次。中轉(zhuǎn)中心配置面積50-70平方米,日處理能力800-1400個(gè)包裹,由于門店處于不同城市區(qū)域,租金略有不同,租金區(qū)間大致為每天3.5-5元/平方米,中轉(zhuǎn)中心除租金外固定開支1000元/天,包括了人員工資/水電費(fèi)和設(shè)備使用費(fèi)等。中轉(zhuǎn)中心到末端需求點(diǎn)的配送車輛多采用面包車或摩托車,容量平均在80個(gè)包裹,車輛固定成本與人工成本合計(jì)為210元/天,配送車輛變動成本3.5元/公里。根據(jù)前文設(shè)計(jì)的混合算法,結(jié)合相關(guān)文獻(xiàn)參數(shù),本文設(shè)定遺傳算法種群規(guī)模100,交叉概率0.5,變異概率0.1,禁忌表長度20,禁忌搜索400,迭代次數(shù)20,對實(shí)例進(jìn)行求解。將以上參數(shù)、測量數(shù)據(jù)帶入模型,編程運(yùn)行15次,剔除5個(gè)相同結(jié)果,得出10個(gè)不同方案(表4)。

    5.2結(jié)果討論

    5.1采用車型為福田歐馬可1.49T廂式輕卡,最大容納包裹1000個(gè),由于中轉(zhuǎn)中心處理量介于800-1400之間,車型改變會對中轉(zhuǎn)中心選擇和末端需求點(diǎn)分配產(chǎn)生影響。本節(jié)以慶鈴五十鈴1.25T廂式輕卡,福田歐馬可2.0T廂式輕卡為例進(jìn)行測算。表3為3種車型參數(shù)對比。

    表3 車型參數(shù)對比

    表5是基于1.25T運(yùn)輸車型得出的10個(gè)中轉(zhuǎn)中心開放方案,其中最大成本30548元,最小成本29290元,平均成本29878.2元,平均成本較5.1上升2.3%,平均運(yùn)輸距離204.83千米,增加32.7%,平均配送距離654.2千米,增加2.9%。表6是基于2.0T運(yùn)輸車型得出的10個(gè)中轉(zhuǎn)中心開放方案,最大成本29099元,最小成本27796元,平均成本28507.2元,較5.1平均成本下降2.43%,平均運(yùn)輸距離139.64千米,降低9.5%,平均配送距離597.58千米,下降6.0%。由于通常狀況下中轉(zhuǎn)中心包裹日處理量處于800-1400,因此采用實(shí)際容量小于800車型會增加運(yùn)輸車數(shù)量,增加總運(yùn)輸距離,增加固定成本和運(yùn)輸成本,大于1400車型會增加車輛固定成本。

    根據(jù)表3、表5和表6的中轉(zhuǎn)中心選址及末端需求點(diǎn)分配優(yōu)化方案,得出城市物流總成本中,配送中心到中轉(zhuǎn)中心的運(yùn)輸成本、中轉(zhuǎn)中心的運(yùn)營成本和中轉(zhuǎn)中心到末端需求點(diǎn)的配送成本占比分別約為6.7%、33.3%、60%?!白詈笠还铩蹦┒伺渌统杀菊颊麄€(gè)城市配送成本的60%,同時(shí)影響客戶體驗(yàn),是城市配送物流的關(guān)鍵環(huán)節(jié)。

    表4 中轉(zhuǎn)中心開放方案

    注1:南坪1對應(yīng)南坪店,南坪2對應(yīng)南坪江南大道店,南坪3對應(yīng)南坪西路店;楊家坪1對應(yīng)楊建設(shè)廣場店,楊家坪2對應(yīng)楊家坪店。

    2:運(yùn)輸距離指配送中心到所有開放中轉(zhuǎn)中心的總距離;配送距離是所有開放中轉(zhuǎn)中心到其服務(wù)末端需求點(diǎn)的配送總距離。下同。

    表5 基于慶鈴五十鈴1.25T廂式輕卡的中轉(zhuǎn)中心開放方案

    表6 基于福田歐馬可2.0T廂式輕卡的中轉(zhuǎn)中心開放方案

    表3、表5和表6方案中,蘇寧重慶主城核心區(qū)域的城市物流年總成本最大約為1150.02萬元,最小約為1014.55萬元。2014年蘇寧重慶線上線下總營業(yè)額54億元,主城占比約為65%,主城核心區(qū)域線上銷售占比約為20%,線上物流費(fèi)用占線上營業(yè)額約2.5%。主城核心區(qū)需求占主城需求量超過75%,以此計(jì)算2014年線上部分城市物流總成本約為1755萬元,該項(xiàng)費(fèi)用包含了城市配送中心的運(yùn)營成本,其占比約10%,因此蘇寧重慶主城核心區(qū)域城市物流總成本(不包括配送中心運(yùn)營成本)約1184.6萬元。本文優(yōu)化最大成本、最小成本相比現(xiàn)有狀況分別降低2.92%,14.36%。根據(jù)當(dāng)前測算每單包裹從配送中心到客戶的城市配送成本約為4.6元。按照當(dāng)前重慶市內(nèi)電商同城配送次日送達(dá)的時(shí)效要求條件下,外包物流訂單一般1kg以下6元,因此對于蘇寧重慶來說主城核心區(qū)域乃至主城區(qū)采用自有物流配送在提升效率的同時(shí)也有效節(jié)約成本。

    6 結(jié)語

    本文以蘇寧等國內(nèi)大型零售企業(yè)O2O轉(zhuǎn)型實(shí)踐為背景,對具有容量限制的城市配送中轉(zhuǎn)中心選址、末端需求點(diǎn)分配和末端配送問題進(jìn)行研究。不同于動態(tài)車輛路徑問題,本文考慮城市核心區(qū)域“最后一公里”配送問題,在配送距離的估算上采用了由Daganzo[24]提出的距離估算方式。構(gòu)建了中轉(zhuǎn)中心選址及末端需求點(diǎn)分配的聯(lián)合優(yōu)化模型。設(shè)計(jì)了基于啟發(fā)式規(guī)則的遺傳算法和禁忌搜索算法的混合算法,并通過混合算法和CPLEX的對比驗(yàn)證本文算法的有效性。以蘇寧重慶主城核心區(qū)域的O2O實(shí)踐為基礎(chǔ),運(yùn)用本文設(shè)計(jì)的算法給出了具體的中轉(zhuǎn)中心選址及末端需求點(diǎn)分配優(yōu)化方案,并給出了各個(gè)成本、運(yùn)輸距離和配送距離的具體數(shù)值,并與蘇寧重慶主城核心區(qū)域的實(shí)際運(yùn)作進(jìn)行對比,結(jié)論表明在不考慮其他因素情況下,本文所得方案城市物流總成本相比現(xiàn)有狀況最少降低2.92% ,最多降低14.36%。論文測算了三種不同車型對中轉(zhuǎn)中心選址方案產(chǎn)生的影響,得出選擇與中轉(zhuǎn)中心容量接近車型能夠有效降低物流成本,減少運(yùn)輸、配送距離。另外得出運(yùn)輸、運(yùn)營、配送費(fèi)用分別占到物流成本的比例約為6.7%,33.3%,60%等結(jié)論。城市交通限行及碳排放等問題都會對決策產(chǎn)生影響,考慮這些因素所得結(jié)論更具有參考價(jià)值,后續(xù)研究可以針對這些問題進(jìn)行探討。

    [1] Aikens C H.Facility location models for distribution Planning[J].European Journal of Operational Research,1985,22(3):263-279.

    [2] Shen Z J M,Coullard C,Daskin M S.A joint location-Inventory model[J].Transportation Science,2003,37(1):40-55.

    [3] Tsao Y C,Mangotra D,Lu J C,et al.A continuous approximation approach for the integrated facility-inventory allocation problem[J].European Journal of Operational Research,2012,22(2):216-228.

    [4] Taniguchi E,Noritake M,Yamada T,et al.Optimal size and location planning of public logistics terminals[J].Transportation Research,1999,35(3):207-222.

    [5] Wang Yong,Ma Xiaolei,Xu Maozeng,et al.Two-echelon logistics distribution region partitioning problem based on a hybrid particle swarm optimization-genetic algorithm[J].Expert Systems with Applications,2015,42(12):5019-5031.

    [6] Diabat A,Batta?a O,Nazzal D.An improved Lagrangian relaxation-based heuristic for a joint location-inventory problem[J].Computers &Operations Research,2015,61:170-178.

    [7] Ahmadi-Javid A,Hoseinpour P.A Location-inventory-pricing model in a supply chain distribution network with price-sensitive demands and inventory-capacity constraints[J].Transportation Research Part E:Logistics and Transportation Review,2015,82:238-255.

    [8] Jouzdani J,Sadjadi S J,Fathian M.Dynamic dairy facility location and supply chain planning under traffic congestion and demand uncertainty: A case of Tehran[J].Applied Mathematical Modeling,2013,37(18-19) :8467-8483.

    [9] 唐金環(huán),戢守峰,朱寶琳.考慮碳配額差值的選址-路徑-庫存集成問題優(yōu)化模型與算法[J].中國管理科學(xué),2014,22(9):114-122.

    [10] 楊珺,盧巍.低碳政策下多容量等級選址與配送問題研究[J].中國管理科學(xué),2014,22(5):51-60.

    [11] 趙泉午,楊茜.考慮CO2排放量的城市專業(yè)物流中心選址研究[J].中國管理科學(xué),2014,22(7):124-130.

    [12] Diabat A,Al-Salem M.An integrated supply chain problem with environmental considerations[J].International Journal of Production Economics,2015,164:330-338.

    [13] 劉必爭,毛超.電子商務(wù)下的配送中心選址問題及其優(yōu)化[J].系統(tǒng)工程,2008,26(10):17-20.

    [14] 周翔,許茂增,呂奇光.B2C模式下配送中心與末端節(jié)點(diǎn)的兩階段布局優(yōu)化模型[J].計(jì)算機(jī)集成制造系統(tǒng),2014,20(12):3140-3149.

    [15] 蔣忠中,汪定偉.B2C電子商務(wù)中配送中心選址問題及其優(yōu)化[J].控制與決策,2005,20(10):1125-1128.

    [16] Xu Bin,Wei Chunli,Shi Zhanjiang,et al.Study on location-selection of B2C e-commerce logistics distribution center[C]//Proceedings of 2009 International Forum on Computer Science-Technology and Applications,Chongqing,China,December 25-27,2009.

    [17] Torabi S A,Hassini E,Jeihoonian M.Fulfillment source allocation,inventory transshipment,and customer order transfer in e-tailing[J].Transportation Research Part E:Logistics and Transportation Review,2015,79:128-144.

    [18] Yu V F,Lin S Y.A simulated annealing heuristic for the open location-routing problem[J].Computers &Industrial Engineering,2015,62:184-196.

    [19] 關(guān)菲,張強(qiáng).模糊多目標(biāo)物流配送中心選址模型及其求解算法[J].中國管理科學(xué),2013,21(S1):57-62.

    [20] Lau H C W,Jiang Zhongzhong,Ip W H,et al.A credibility-based fuzzy location model with Huriwicz criteria for the desingn of distribution systems in B2C e-commerce[J].Computers &Industrial Engineering,2010,59:873-886.

    [21] 張曉楠,范厚明,李劍鋒.B2C物流配送網(wǎng)絡(luò)雙目標(biāo)模糊選址模型與算法[J].系統(tǒng)工程理論與實(shí)踐,2015,35(5):1202-1213.

    [22] Nadizadeh A,Nasab HH.Solving the dynamic capacitated location-routing problem with fuzzy demands by hybrid heuristic algorithm[J].European Journal of Operational Research,2014,238(2):458-470.

    [23] Eilon S,Waston-Gandy C D T,Christofides N.Distribution management: Mathematical modelling and practical analysis [M].New York:Hafner,1971.

    [24] Daganzo C F.The distance traveled to visit N points with a maximum of C stops per vehicle: An Analytic Model and an Application [J].Transportation Science,1984,18(4):331-350.

    [25] Figliozzi M A.Planning approximations to the average length of vehicles routing problems with varying customer demands and routing constraints [J].Transportation Research Record: Journal of the Transportation Research Board,2008,2089:1-8.

    [26] Davis B A,Figliozzi M A.A methodology to evaluate the competitiveness of electric delivery trucks[J].Transportation Research Part E:Logistics and Transportation Review,2013 49:8-23.

    [27] 周愉峰,馬祖軍,王恪銘.應(yīng)急物資儲備庫的可靠性P-中位選址模型[J].管理評論,2015,27(5):198-208.

    Abstract: In this article,joint location-assignment problem is studied for designing an Urban Logistics Network (ULN) with multiple intermediate depots (IDs) and terminals of ULN.The key decisions are where to locate the intermediate depots (IDs) and how to assign terminals to IDs such that the total network cost is minimized.A large-scale static and deterministic integer programming model is presented solving a joint location-assignment problem of large chain retailers.To solve this model,a more efficient hybrid algorithm integrated with Genetic algorithm and Tabu search algorithm is put forward based on heuristic searching criteria.Our methodology is illustrated with the Urban Logistics Network from a leading Chinese retailer (Suning) in Chongqing.Numerical analysis suggests that optimal solutions can reduce the total network cost between 2.92%~14.36% by comparison with the current ULN .Delivery costs of “Last Mile” account for total urban logistics network costs more than 60%.Vehicle type is a deterministic factor to total urban logistics network costs.The effects of changing parameter values on the optimal solutions are also studied and some management implications are pointed out.

    Keywords: online-offline channel integration;iocation-assignment problem;“l(fā)ast mile ";Genetic algorithm;tabu search algorithm

    A City Logistics Network Optimization Model for Large Chain Retailers under Online-Offline Channel Integration

    ZHAOQuan-wu1,2,ZHAOJun-ping1,2,LINYa2

    (School of Economics and Business Administration,Chongqing University,Chongqing 400030,China;Chongqing Key Laboratory of Logistics,Chongqing University,Chongqing 400030,China)

    U121;O221

    A

    1003-207(2017)09-0159-09

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

    2016-01-29;

    2016-10-14

    重慶市社會科學(xué)規(guī)劃一般項(xiàng)目(2015YBGL109);國家自然科學(xué)基金資助項(xiàng)目(71002070);國家社會科學(xué)基金重點(diǎn)項(xiàng)目(14AGL023);國家科技支撐計(jì)劃資助項(xiàng)目(2015BAH46F01)

    趙泉午(1976-),男(漢族),河南方城人,重慶大學(xué)經(jīng)濟(jì)與工商管理學(xué)院教授,博士,研究方向:現(xiàn)代物流與供應(yīng)鏈管理、營銷渠道及零售管理,E-mail:zhaoquanwu@126.com.

    猜你喜歡
    分配距離物流
    應(yīng)答器THR和TFFR分配及SIL等級探討
    本刊重點(diǎn)關(guān)注的物流展會
    遺產(chǎn)的分配
    一種分配十分不均的財(cái)富
    “智”造更長物流生態(tài)鏈
    汽車觀察(2018年12期)2018-12-26 01:05:44
    績效考核分配的實(shí)踐與思考
    每次失敗都會距離成功更近一步
    山東青年(2016年3期)2016-02-28 14:25:55
    基于低碳物流的公路運(yùn)輸優(yōu)化
    愛的距離
    母子健康(2015年1期)2015-02-28 11:21:33
    距離美
    絲綢之路(2014年21期)2014-12-31 00:00:00
    校园人妻丝袜中文字幕| 国产永久视频网站| 一个人看视频在线观看www免费| 欧美极品一区二区三区四区| 五月天丁香电影| 久久久精品欧美日韩精品| 国产三级在线视频| 国产熟女欧美一区二区| 久久这里有精品视频免费| 少妇人妻精品综合一区二区| 免费观看性生交大片5| 男女下面进入的视频免费午夜| 99久久九九国产精品国产免费| 中文字幕久久专区| 水蜜桃什么品种好| 一二三四中文在线观看免费高清| 99热网站在线观看| 国产黄色免费在线视频| 一个人看的www免费观看视频| 精品久久久久久久人妻蜜臀av| 久久久久久久国产电影| 国产精品久久久久久精品电影小说 | 最近中文字幕2019免费版| 国产成人a区在线观看| 韩国av在线不卡| 色播亚洲综合网| 最新中文字幕久久久久| 日韩欧美 国产精品| 色5月婷婷丁香| 午夜精品国产一区二区电影 | freevideosex欧美| 网址你懂的国产日韩在线| 久久久精品免费免费高清| 最近最新中文字幕大全电影3| 午夜激情福利司机影院| 久久午夜福利片| 久久久久久久午夜电影| 久久国产乱子免费精品| 日韩av在线免费看完整版不卡| 欧美激情久久久久久爽电影| av在线播放精品| 亚洲国产最新在线播放| 国语对白做爰xxxⅹ性视频网站| 色播亚洲综合网| 欧美bdsm另类| 夜夜看夜夜爽夜夜摸| 可以在线观看毛片的网站| 亚洲人成网站在线播| 成人二区视频| 晚上一个人看的免费电影| 亚洲精品中文字幕在线视频 | 亚洲经典国产精华液单| 久久这里有精品视频免费| 在线观看人妻少妇| 久久精品国产亚洲av天美| 国产麻豆成人av免费视频| 日韩亚洲欧美综合| 精品久久久精品久久久| 搡老妇女老女人老熟妇| 亚洲国产欧美在线一区| 国产精品一及| 少妇人妻一区二区三区视频| 男人舔奶头视频| 亚洲精品aⅴ在线观看| 波野结衣二区三区在线| 尾随美女入室| 国产毛片a区久久久久| 内射极品少妇av片p| 国产成人精品福利久久| 国内精品一区二区在线观看| 欧美日本视频| 人人妻人人看人人澡| 中文在线观看免费www的网站| 蜜桃久久精品国产亚洲av| 建设人人有责人人尽责人人享有的 | 国产av不卡久久| 色综合亚洲欧美另类图片| 日韩不卡一区二区三区视频在线| 2022亚洲国产成人精品| 日本黄大片高清| av福利片在线观看| 日韩,欧美,国产一区二区三区| 国产伦在线观看视频一区| 国产精品人妻久久久久久| 秋霞伦理黄片| av专区在线播放| 成人综合一区亚洲| 国产亚洲av片在线观看秒播厂 | 亚洲欧美成人精品一区二区| 免费少妇av软件| 嘟嘟电影网在线观看| 国产黄色视频一区二区在线观看| 欧美xxxx黑人xx丫x性爽| 日本免费a在线| a级毛片免费高清观看在线播放| 亚洲国产欧美在线一区| 日韩制服骚丝袜av| 精华霜和精华液先用哪个| 亚洲精品久久久久久婷婷小说| 一区二区三区高清视频在线| 午夜福利在线观看免费完整高清在| 色综合色国产| 黄色配什么色好看| 麻豆成人av视频| 80岁老熟妇乱子伦牲交| 99热这里只有是精品50| 天堂影院成人在线观看| 亚洲图色成人| 精品久久久久久成人av| 啦啦啦啦在线视频资源| 超碰97精品在线观看| 久久99热6这里只有精品| 男人爽女人下面视频在线观看| 国模一区二区三区四区视频| 色播亚洲综合网| 亚洲婷婷狠狠爱综合网| 一区二区三区免费毛片| 国产在视频线精品| 国模一区二区三区四区视频| 成人无遮挡网站| 最近手机中文字幕大全| 在线 av 中文字幕| 午夜爱爱视频在线播放| 中文乱码字字幕精品一区二区三区 | 麻豆精品久久久久久蜜桃| 亚洲精品乱码久久久v下载方式| 日本一本二区三区精品| 日韩欧美一区视频在线观看 | 身体一侧抽搐| 日韩欧美精品v在线| 26uuu在线亚洲综合色| 久久精品国产自在天天线| 少妇被粗大猛烈的视频| 久久久久久伊人网av| 高清av免费在线| or卡值多少钱| 一级二级三级毛片免费看| 久久精品综合一区二区三区| 日韩成人av中文字幕在线观看| 两个人视频免费观看高清| 欧美日韩一区二区视频在线观看视频在线 | 亚洲18禁久久av| 五月玫瑰六月丁香| 国产 亚洲一区二区三区 | 国产 一区精品| 少妇人妻一区二区三区视频| 黄片wwwwww| 婷婷六月久久综合丁香| 最新中文字幕久久久久| 亚洲欧美成人综合另类久久久| 精品久久久久久久久av| 插阴视频在线观看视频| 亚洲四区av| 简卡轻食公司| 在线免费十八禁| 国产人妻一区二区三区在| 在线免费十八禁| 欧美xxxx性猛交bbbb| 非洲黑人性xxxx精品又粗又长| 99热6这里只有精品| 国产精品美女特级片免费视频播放器| 免费看日本二区| 国产久久久一区二区三区| 高清在线视频一区二区三区| 亚洲在线观看片| 精品一区二区三卡| 精品国内亚洲2022精品成人| 午夜免费男女啪啪视频观看| 国产在线男女| 欧美激情国产日韩精品一区| 国产白丝娇喘喷水9色精品| 国产精品.久久久| 国产高潮美女av| 人人妻人人澡欧美一区二区| 国产一区二区在线观看日韩| 免费播放大片免费观看视频在线观看| 日韩伦理黄色片| 精品一区二区免费观看| 久久久久久久久久黄片| 少妇的逼好多水| 又粗又硬又长又爽又黄的视频| 色吧在线观看| 久久人人爽人人爽人人片va| 亚洲伊人久久精品综合| 久久精品夜夜夜夜夜久久蜜豆| 久久精品久久久久久噜噜老黄| www.色视频.com| 免费播放大片免费观看视频在线观看| 亚洲av男天堂| 亚洲精品乱久久久久久| 一级毛片 在线播放| 国国产精品蜜臀av免费| 激情五月婷婷亚洲| 国产精品精品国产色婷婷| 一个人观看的视频www高清免费观看| 六月丁香七月| 一本一本综合久久| 久久精品国产亚洲网站| 女的被弄到高潮叫床怎么办| 老女人水多毛片| 国产亚洲av片在线观看秒播厂 | 大香蕉久久网| 九草在线视频观看| 精品一区二区三区人妻视频| 亚洲真实伦在线观看| 国产精品无大码| 高清日韩中文字幕在线| 日韩av免费高清视频| 日本爱情动作片www.在线观看| 国产精品久久视频播放| 国产av国产精品国产| 日韩亚洲欧美综合| 国产久久久一区二区三区| 日本熟妇午夜| freevideosex欧美| 热99在线观看视频| 九色成人免费人妻av| 免费观看性生交大片5| 一本久久精品| 日韩一区二区三区影片| 中文在线观看免费www的网站| 久久久精品免费免费高清| 欧美不卡视频在线免费观看| 亚洲av男天堂| av卡一久久| 亚洲图色成人| 看十八女毛片水多多多| 精品一区二区免费观看| 欧美zozozo另类| 啦啦啦啦在线视频资源| 看黄色毛片网站| 国产有黄有色有爽视频| 深夜a级毛片| 日韩强制内射视频| 大香蕉97超碰在线| 高清毛片免费看| 色综合色国产| 免费高清在线观看视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 乱人视频在线观看| 国产人妻一区二区三区在| 成人美女网站在线观看视频| 91狼人影院| 亚洲激情五月婷婷啪啪| 国产精品一区二区三区四区久久| 一个人免费在线观看电影| 又爽又黄a免费视频| 小蜜桃在线观看免费完整版高清| 天天一区二区日本电影三级| 国产熟女欧美一区二区| 97人妻精品一区二区三区麻豆| 精品久久久久久久久av| 美女cb高潮喷水在线观看| 国产精品国产三级国产专区5o| 国产精品蜜桃在线观看| 成人欧美大片| av福利片在线观看| 亚洲成人精品中文字幕电影| 99久久中文字幕三级久久日本| av在线天堂中文字幕| 亚洲国产最新在线播放| 欧美3d第一页| 国产免费又黄又爽又色| 一个人看视频在线观看www免费| 国产午夜福利久久久久久| 亚洲成人中文字幕在线播放| 久久久国产一区二区| 少妇猛男粗大的猛烈进出视频 | 天堂√8在线中文| 国产黄色小视频在线观看| 久久久久网色| 别揉我奶头 嗯啊视频| 91精品一卡2卡3卡4卡| 亚洲成色77777| 亚洲精品,欧美精品| 亚洲一级一片aⅴ在线观看| 狂野欧美激情性xxxx在线观看| 久久久久久国产a免费观看| 亚洲av电影在线观看一区二区三区 | 中文字幕免费在线视频6| 成人综合一区亚洲| 精品不卡国产一区二区三区| 国产精品国产三级专区第一集| 在线免费观看的www视频| 超碰av人人做人人爽久久| 午夜精品在线福利| 熟妇人妻不卡中文字幕| 人人妻人人看人人澡| 精品国产一区二区三区久久久樱花 | 人体艺术视频欧美日本| 国内揄拍国产精品人妻在线| 国产黄色小视频在线观看| 建设人人有责人人尽责人人享有的 | 欧美 日韩 精品 国产| 亚洲人成网站在线播| 国产精品人妻久久久久久| 国产免费福利视频在线观看| 精品久久久久久久久亚洲| 一级二级三级毛片免费看| 少妇熟女欧美另类| 亚洲人成网站在线观看播放| freevideosex欧美| 国产亚洲午夜精品一区二区久久 | 一个人观看的视频www高清免费观看| 美女高潮的动态| 国产成人午夜福利电影在线观看| 午夜福利网站1000一区二区三区| 亚洲国产精品专区欧美| 综合色丁香网| 成年免费大片在线观看| 成人午夜精彩视频在线观看| 国产精品一区二区性色av| 国内揄拍国产精品人妻在线| 欧美激情在线99| 亚洲av不卡在线观看| 国产探花极品一区二区| 中文欧美无线码| 超碰av人人做人人爽久久| 韩国av在线不卡| 中国美白少妇内射xxxbb| 亚洲av国产av综合av卡| 亚洲精品国产成人久久av| 免费黄色在线免费观看| 国产探花在线观看一区二区| 国产一级毛片七仙女欲春2| 99热网站在线观看| 能在线免费看毛片的网站| 亚洲在线自拍视频| 国产精品.久久久| av国产免费在线观看| 久久久精品欧美日韩精品| 久久精品夜色国产| 看黄色毛片网站| 欧美bdsm另类| 免费观看无遮挡的男女| 男女那种视频在线观看| 天堂中文最新版在线下载 | 欧美xxⅹ黑人| or卡值多少钱| 九九爱精品视频在线观看| 国产有黄有色有爽视频| 丰满乱子伦码专区| 99热这里只有精品一区| 少妇熟女aⅴ在线视频| 日本wwww免费看| 能在线免费看毛片的网站| 乱码一卡2卡4卡精品| 人体艺术视频欧美日本| 丰满人妻一区二区三区视频av| 18禁在线无遮挡免费观看视频| 国产在线一区二区三区精| 波野结衣二区三区在线| 亚洲国产欧美在线一区| 精品久久久精品久久久| 晚上一个人看的免费电影| 内射极品少妇av片p| 亚洲精品日韩av片在线观看| 啦啦啦中文免费视频观看日本| 免费看av在线观看网站| 亚洲精品乱码久久久v下载方式| 可以在线观看毛片的网站| 日韩av不卡免费在线播放| 欧美日韩视频高清一区二区三区二| 狂野欧美白嫩少妇大欣赏| 少妇的逼水好多| 日韩av在线免费看完整版不卡| av在线天堂中文字幕| 国产69精品久久久久777片| 亚洲,欧美,日韩| 久久精品久久久久久噜噜老黄| 卡戴珊不雅视频在线播放| 高清毛片免费看| 在线a可以看的网站| 99热网站在线观看| 狂野欧美白嫩少妇大欣赏| 中文字幕av在线有码专区| 少妇人妻一区二区三区视频| 日日干狠狠操夜夜爽| 美女高潮的动态| 成人高潮视频无遮挡免费网站| 热99在线观看视频| 禁无遮挡网站| 国产欧美日韩精品一区二区| 中文欧美无线码| 一个人免费在线观看电影| 亚洲精品第二区| 小蜜桃在线观看免费完整版高清| 亚洲精品视频女| av在线天堂中文字幕| 晚上一个人看的免费电影| 成年版毛片免费区| 亚州av有码| 非洲黑人性xxxx精品又粗又长| 十八禁网站网址无遮挡 | 伊人久久精品亚洲午夜| 又大又黄又爽视频免费| 国产亚洲精品久久久com| 国产欧美日韩精品一区二区| 国产三级在线视频| 亚洲婷婷狠狠爱综合网| 国产日韩欧美在线精品| 日本与韩国留学比较| 亚洲欧美日韩卡通动漫| 久久久久久久久久久丰满| 成年女人看的毛片在线观看| 日韩av免费高清视频| 狂野欧美白嫩少妇大欣赏| 久久久久免费精品人妻一区二区| 日韩伦理黄色片| 婷婷色综合www| a级一级毛片免费在线观看| 街头女战士在线观看网站| 免费观看无遮挡的男女| 十八禁国产超污无遮挡网站| 亚洲精品aⅴ在线观看| 搡老乐熟女国产| 亚洲国产成人一精品久久久| av天堂中文字幕网| 国产男人的电影天堂91| 一级黄片播放器| 精品少妇黑人巨大在线播放| 亚洲美女视频黄频| 免费观看的影片在线观看| 久久久精品免费免费高清| 亚洲欧美精品专区久久| 亚洲婷婷狠狠爱综合网| 免费少妇av软件| 亚洲成人精品中文字幕电影| 午夜福利视频精品| 国产成人a区在线观看| 亚洲精品亚洲一区二区| 69人妻影院| 青春草视频在线免费观看| 国产精品福利在线免费观看| 久久午夜福利片| 亚洲精华国产精华液的使用体验| 五月伊人婷婷丁香| 最后的刺客免费高清国语| 秋霞伦理黄片| 哪个播放器可以免费观看大片| 97热精品久久久久久| 久久久久久久亚洲中文字幕| 男人舔奶头视频| 国产免费视频播放在线视频 | 国产一区二区在线观看日韩| 亚洲精品国产av成人精品| 日韩制服骚丝袜av| 国产精品熟女久久久久浪| 亚洲人成网站在线观看播放| 少妇被粗大猛烈的视频| 99久国产av精品国产电影| 日韩成人伦理影院| 久久99热这里只有精品18| 人妻少妇偷人精品九色| 黄色一级大片看看| 国产男人的电影天堂91| 免费看av在线观看网站| 国产精品久久久久久久电影| 日韩强制内射视频| 日日啪夜夜爽| 国产黄色免费在线视频| 天天躁日日操中文字幕| 免费人成在线观看视频色| 91av网一区二区| 亚洲国产色片| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美一区二区三区国产| 日韩国内少妇激情av| 国产成人精品婷婷| 欧美日韩亚洲高清精品| av在线蜜桃| 一区二区三区四区激情视频| 久久久久久久亚洲中文字幕| 街头女战士在线观看网站| 美女cb高潮喷水在线观看| 中国美白少妇内射xxxbb| 亚洲一级一片aⅴ在线观看| 蜜桃亚洲精品一区二区三区| 欧美xxⅹ黑人| 最近视频中文字幕2019在线8| 乱码一卡2卡4卡精品| 国产午夜精品久久久久久一区二区三区| 大香蕉97超碰在线| 国产精品一区二区三区四区免费观看| 精品久久国产蜜桃| 中文在线观看免费www的网站| 日韩国内少妇激情av| av播播在线观看一区| 精品一区二区三区视频在线| 人妻制服诱惑在线中文字幕| 男女边吃奶边做爰视频| 亚洲图色成人| 亚洲精品久久久久久婷婷小说| 插阴视频在线观看视频| 伦精品一区二区三区| 国产男人的电影天堂91| 别揉我奶头 嗯啊视频| 亚洲精品日韩av片在线观看| 欧美xxxx性猛交bbbb| 国产熟女欧美一区二区| 少妇裸体淫交视频免费看高清| 夜夜爽夜夜爽视频| 又大又黄又爽视频免费| 亚洲成人一二三区av| 中文字幕av成人在线电影| 国模一区二区三区四区视频| 伦精品一区二区三区| 精品久久久久久久久av| 别揉我奶头 嗯啊视频| 日产精品乱码卡一卡2卡三| 男人和女人高潮做爰伦理| 免费看光身美女| 国产精品日韩av在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 美女大奶头视频| 韩国av在线不卡| 国产一级毛片七仙女欲春2| 卡戴珊不雅视频在线播放| 亚洲欧美成人精品一区二区| 高清av免费在线| 99久久人妻综合| 国产一区二区三区av在线| 亚洲国产色片| 插逼视频在线观看| 69人妻影院| 国产午夜精品论理片| 精品久久国产蜜桃| 亚洲在线观看片| 最近最新中文字幕免费大全7| 国产v大片淫在线免费观看| 亚洲精品乱码久久久久久按摩| 久久久久九九精品影院| 欧美激情国产日韩精品一区| 日韩av免费高清视频| 亚洲人成网站高清观看| 国产精品一及| av线在线观看网站| 五月伊人婷婷丁香| 免费在线观看成人毛片| 高清午夜精品一区二区三区| 丝袜美腿在线中文| 久久精品国产亚洲av涩爱| 久久久精品欧美日韩精品| 免费人成在线观看视频色| 国产精品国产三级国产av玫瑰| 蜜桃久久精品国产亚洲av| 国产av在哪里看| 国产黄色视频一区二区在线观看| 天堂中文最新版在线下载 | 国产永久视频网站| 男人和女人高潮做爰伦理| 色网站视频免费| 美女内射精品一级片tv| 久久鲁丝午夜福利片| 一夜夜www| 成人高潮视频无遮挡免费网站| 嫩草影院新地址| 十八禁网站网址无遮挡 | 欧美日韩视频高清一区二区三区二| 赤兔流量卡办理| 床上黄色一级片| 一级爰片在线观看| 高清日韩中文字幕在线| 成年av动漫网址| 亚洲精品日韩av片在线观看| 嘟嘟电影网在线观看| 日本色播在线视频| 亚洲精品中文字幕在线视频 | 夫妻性生交免费视频一级片| 一二三四中文在线观看免费高清| 国产精品1区2区在线观看.| 精品久久久久久久末码| 三级男女做爰猛烈吃奶摸视频| 国产黄片视频在线免费观看| 亚洲av中文字字幕乱码综合| 在线免费观看的www视频| 中文字幕亚洲精品专区| 亚洲激情五月婷婷啪啪| 99久久精品国产国产毛片| 熟妇人妻久久中文字幕3abv| 观看美女的网站| 搞女人的毛片| 亚洲精品久久午夜乱码| 国产探花极品一区二区| 69人妻影院| 亚洲欧美日韩东京热| 国产av国产精品国产| 欧美激情久久久久久爽电影| 国产av国产精品国产| 少妇高潮的动态图| 天天躁日日操中文字幕| 九九在线视频观看精品| 婷婷色综合www| 身体一侧抽搐| 熟女人妻精品中文字幕| 成人毛片a级毛片在线播放| 国产片特级美女逼逼视频| 欧美精品一区二区大全| 国产麻豆成人av免费视频| 亚洲av国产av综合av卡| 青春草国产在线视频| 中文字幕av成人在线电影| 久久精品国产亚洲网站| 亚洲精品,欧美精品| 成人二区视频| 少妇丰满av| 天堂影院成人在线观看| 精品午夜福利在线看| 国产色婷婷99| 淫秽高清视频在线观看| 熟妇人妻久久中文字幕3abv| 天美传媒精品一区二区| 免费看光身美女| 免费观看无遮挡的男女| 国产精品一区二区在线观看99 |