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

    集裝箱港口岸橋動態(tài)分配與貝位調(diào)度研究

    2016-05-25 00:37:00梁承姬王丹丹
    關(guān)鍵詞:離港進港集裝箱

    梁承姬,王丹丹

    (上海海事大學 物流研究中心,上海 201306)

    集裝箱港口岸橋動態(tài)分配與貝位調(diào)度研究

    梁承姬,王丹丹

    (上海海事大學 物流研究中心,上海 201306)

    為提高集裝箱港口岸橋的利用率和對船舶的服務(wù)水平,針對多船各貝位的任務(wù)量,考慮船舶優(yōu)先度、可作業(yè)岸橋數(shù)量、岸橋沖突和安全距離等因素,以船舶的在港時間最短為目標,建立了多船岸橋動態(tài)分配與貝位調(diào)度混合整數(shù)規(guī)劃模型,通過不同規(guī)模的算例,分別用C#調(diào)用Cplex和C#對遺傳算法編程求解。結(jié)果表明:該模型能有效地解決多船的岸橋分配與調(diào)度問題,且遺傳算法求解時間更短、求解質(zhì)量更高,進而驗證了該模型和算法的有效性。

    交通運輸工程;岸橋動態(tài)分配;岸橋貝位調(diào)度;岸橋沖突;安全距離;遺傳算法

    0 引 言

    集裝箱港口作為現(xiàn)代物流的重要組成部分,在全球物流系統(tǒng)中起著重要的作用。隨著競爭的加劇,各個港口都面臨著提高服務(wù)質(zhì)量、降低服務(wù)成本、增加吞吐量的壓力。要想緩解這些壓力,必須有效利用港口的各種資源,如泊位、岸橋、場橋、堆存、集卡等。岸橋是港口的昂貴資源,也是稀缺資源,它的有效利用對于港口的運作效率至關(guān)重要。岸橋資源計劃分為岸橋分配計劃(Quay Crane Allocation Problem, QCAP)和岸橋調(diào)度計劃(Quay Crane Scheduling Problem, QCSP),兩者緊密聯(lián)系,同時考慮岸橋分配和調(diào)度問題可以更加有效的利用岸橋資源。集裝箱港口的岸橋分配與調(diào)度問題(Quay Crane Allocation and Scheduling Problem, QCASP)研究如何為集裝箱船舶分配岸橋以及安排每臺岸橋的任務(wù)操作順序,使得在滿足實際約束的情況下盡早完成每艘船舶的裝卸操作。

    對于單船的QCSP,K.H.Kim等[1]提出了分支定界算法和貪婪隨機自適應(yīng)搜索算法。L.Moccia等[2]在K.H.Kim[1]的基礎(chǔ)上,設(shè)計了分支切割算法,改進了解的質(zhì)量。M.Sammarra等[3]把單船的QCSP分解為兩個子問題,在此基礎(chǔ)上設(shè)計了一種禁忌搜索算法。C.Bierirth等[4]設(shè)計了一個基于分支定界的單向調(diào)度啟發(fā)式算法,算例結(jié)果表明該算法能夠在更短時間內(nèi)產(chǎn)生比以往文獻質(zhì)量更好的解。楊明珠[5]也對該問題進行了研究,給出了一種改進的啟發(fā)式算法。F.Meisel等[6]對目前已有的單船岸橋調(diào)度優(yōu)化模型和算法進行了分析和評價。

    對于多船的QCASP,很多學者都是聯(lián)合泊位資源一起考慮的。R.Tavakkoli-Moghaddam[7]擴展了K.H.Kim[1]提出的單船岸橋調(diào)度模型,構(gòu)建了同時考慮一組船舶的混合整數(shù)規(guī)劃(Mixed Integer Programming ,MIP)模型,并設(shè)計了遺傳算法求解該問題。Liang Chengji等[8]和Zhang Canrong等[9]研究了船舶的靠泊位置、作業(yè)時間以及分配給每艘船的岸橋數(shù)量,前者建立了作業(yè)時間、等待時間及延遲時間之和最小的優(yōu)化模型,后者建立了成本最小的優(yōu)化模型,考慮產(chǎn)生成本的因素是船舶靠泊位置、停泊時間和完成時間。Lee Der-horng等[10]建立了求船舶在港時間最小的MIP模型,并證明了該模型是完全NP問題。梁承姬等[11]研究了考慮泊位的多船QCAP,建立了泊位分配MIP模型和岸橋調(diào)度MIP模型,并用Gurobi軟件進行求解。A.Imai等[12]研究了泊位和岸橋的協(xié)調(diào)分配問題,并運用遺傳算法進行求解。

    以上研究大多是將岸橋分配和岸橋調(diào)度分開考慮的,即使考慮岸橋分配與調(diào)度的協(xié)同,也只是考慮單船的貝位調(diào)度問題。為此,筆者將岸橋分配和調(diào)度問題結(jié)合起來,針對多船各個貝位的任務(wù)量動態(tài)地考慮岸橋分配和調(diào)度問題,提出了相應(yīng)的優(yōu)化模型和求解算法,以尋求港口岸邊操作效率最高的岸橋分配和調(diào)度方案。

    1 問題描述

    船舶靠泊之后由岸橋進行集裝箱的裝卸操作,完成裝卸操作之后離港,如圖1。在該過程中,對于某一個確定的計劃時段,如船舶的進港時間、最晚離港時間、最大裝卸時間、船舶上每個貝位的裝卸任務(wù)量都是已知的,且假定船舶按照進港時間在岸邊依次從左到右靠泊。在為船舶分配岸橋以及安排每臺岸橋的任務(wù)操作順序時,不僅要考慮船舶的進港時間、最晚離港時間和船舶優(yōu)先度,還要考慮可作業(yè)岸橋數(shù)量、岸橋沖突以及岸橋安全距離等因素。筆者研究的是對某一個計劃期內(nèi)(通常一個計劃期為12 h),對進港的所有船舶進行動態(tài)的岸橋分配和調(diào)度,在保證船舶按時離港的前提下最小化每艘船舶的在港時間。

    圖1 集裝箱港口岸邊示意Fig.1 The illustration of quay side at container terminal

    2 模型的建立

    2.1 模型假設(shè)

    岸橋在其服務(wù)的船舶沒有全部完成任務(wù)之前可以移動到另外一艘船舶;每臺岸橋擁有相同的工作能力,且可以雙向移動;所有的岸橋都位于一條水平軌道上,并且將其從左到右依次編號;考慮岸橋之間的安全距離,不考慮岸橋的移動時間;按進港先后順序為船舶依次編號,并將所有船舶的所有貝位依次編號;根據(jù)船舶進港時間生成船舶優(yōu)先度序列,即先進港的船舶優(yōu)先度高,其中同時進港的船舶根據(jù)其離港時間產(chǎn)生優(yōu)先度,即先離港的船舶優(yōu)先度高;根據(jù)所有船舶的總?cè)蝿?wù)量安排可作業(yè)的岸橋數(shù)量;每個貝位的任務(wù)量由整數(shù)時間表示,單位為h;每個貝位的任務(wù)可以由多臺岸橋在不同時段完成。

    2.2 集合、參數(shù)及決策變量的定義

    集合:K為船舶集合,由k表示,k代表船舶序號;Jk為船舶k的貝位集合,由j表示;Q為岸橋集合,由q表示,q代表岸橋序號;T為時間集合,由t表示。

    參數(shù):lkj為船舶k的貝位j的位置序號;m為相鄰岸橋所在貝位序號的最小差值,這里定義為最小安全距離;pkj為船舶k的貝位j的任務(wù)量;Dk為船舶k的最大裝卸時間;wk為船舶k的優(yōu)先度。

    Lk為船舶k的進港時間與第一個進港船舶的進港時間之差,這里定義為相距時間。

    決策變量:Xkjqt為若岸橋q在第t時段分配給了船舶k的第j貝位則為1,否則為0;fkjt為若船舶k的第j貝位在第t時段已全部完成任務(wù)則為1,否則為0;rkt為若船舶k在第t時段已全部完成任務(wù)則為1,否則為0。

    2.3 目標函數(shù)和約束條件

    (1)

    (2)

    (3)

    (4)

    ?q∈Q,?t∈T

    (5)

    (6)

    ?k∈K,?j∈Jk,?t∈T

    (7)

    fkjt≥rkt, ?k∈K,?j∈Jk,?t∈T

    (8)

    (9)

    ?k∈K,?j∈Jk

    (10)

    Xkjqt,fkjt,rkt∈{0,1},

    ?k∈K,?j∈Jk,?q∈Q,?t∈T

    (11)

    目標函數(shù)(1)將最小化每艘船舶的在港時間之和,轉(zhuǎn)換為最大化一個計劃期內(nèi)船舶離港之后的時間之和;約束條件(2)保證在任意一個時段一臺岸橋最多只能停放在一個貝位上;約束條件(3)保證在任意一個時段任意一個貝位上最多只能停放一臺岸橋;約束條件(4)保證在任意一個時段正在工作的岸橋總數(shù)量不能超過可用的岸橋數(shù)量;約束條件(5)保證在岸橋從左到右依次排序的情況下,相鄰岸橋之間保持至少m貝位的安全距離;約束條件(6)保證所有的集裝箱任務(wù)在計劃時段內(nèi)都必須完成;約束條件(7)定義了每一個貝位上的工作完成標志,當t時段船舶k的第j貝位上的任務(wù)全部完成時fkjt的值為1,此時在某一貝位上的所有岸橋的時間段總和要大于該貝位的任務(wù)量;約束條件(8)保證某一船舶的所有貝位上的任務(wù)全部完成時才表示該船舶操作完成;約束條件(9)確保所有船舶不超過其最晚離港時間離港;約束條件(10)確保任意船舶的任意貝位是在船舶進港之后開始操作的;約束條件(11)確保所有的決策變量都是0~1變量。

    3 遺傳算法的實現(xiàn)

    遺傳算法(Genetic Algorithm,GA)是一種高度并行、隨機和自適應(yīng)的優(yōu)化算法,其編碼技術(shù)與遺傳操作較簡單,優(yōu)化不受限制性條件的約束,其已經(jīng)在各個領(lǐng)域有了廣泛的應(yīng)用。由于問題的復(fù)雜性,考慮使用遺傳算法來求解。

    3.1 染色體編碼

    圖2 染色體編碼示意Fig.2 The illustration of chromosome coding

    3.2 生成初始種群

    生成初始種群有兩種方法。一種是隨機生成法,該方法有利于保持初始種群的多樣性;另一種是啟發(fā)式方法,它雖然可以生成質(zhì)量較高的初始個體,但可能會導(dǎo)致局部最優(yōu)解。筆者采用隨機生成法生成初始種群,并且只允許可行解成為初始種群的一部分,具體規(guī)則如下:

    2)每行中的整數(shù)順序隨機排列,但要保證每艘船舶的各個貝位要在進港之后開始作業(yè)以及在最晚離港時間之前結(jié)束作業(yè);

    3)每列中的整數(shù)按降序排列,即保證岸橋不跨越交叉,且每列中相鄰整數(shù)的差值不能小于安全距離。

    3.3 適應(yīng)度計算和父代選擇

    根據(jù)目標函數(shù),目標函數(shù)值大的給予高適應(yīng)度值。采用輪盤賭方式產(chǎn)生與種群相同數(shù)量的染色體為父代。

    3.4 遺傳操作

    3.4.1 交叉操作

    采用修補交叉,它基于均勻交叉。即首先生成一個二進制字符串矩陣,其中1和0的個數(shù)相等。如果在該二進制字符串矩陣中的值為1則子代相應(yīng)位置的基因值取父代1相應(yīng)位置的基因值,然后根據(jù)子代已有的基因值劃去父代2的一部分基因值,再用父代2剩余的基因值依次填滿子代剩余的基因位置。如圖3為交叉操作的過程。

    圖3 交叉操作示意Fig.3 The illustration of crossover operation

    3.4.2 變異操作

    采用交換變異,其變異操作過程如圖4。

    圖4 變異操作示意Fig.4 The illustration of mutation operation

    3.5 結(jié)束規(guī)則

    通過以上規(guī)則產(chǎn)生的新的子代將會替代舊種群中適應(yīng)度值最小的染色體,該過程中將保持種群規(guī)模,且迭代到最大迭代次數(shù)時終止迭代。

    4 算例與結(jié)果分析

    為了更好地理解和驗證筆者所提出的模型和算法的有效性,編制了不同規(guī)模的算例。設(shè)定遺傳算法的交叉概率是0.9,變異概率是0.2,最大迭代次數(shù)是500次,種群規(guī)模是150。運行環(huán)境為Intel(R) Core(TM) i5處理器、內(nèi)存4G、硬盤500G的個人計算機,使用C#_2010對遺傳算法編程并計算。

    4.1 兩艘船舶的輸入數(shù)據(jù)與結(jié)果分析

    船舶數(shù)量為2,可用的岸橋數(shù)量為3,安全距離為1,其它輸入數(shù)據(jù)見表1和表2。計算結(jié)果如圖5,可以看出該結(jié)果既滿足模型中船舶的約束也滿足岸橋的約束,且能夠使每艘船舶在最晚離港時間之前盡早離港。圖6所示是將結(jié)果用岸橋和時間的坐標顯示,坐標里的數(shù)字代表貝位編號。

    表1 船舶任務(wù)量數(shù)據(jù)Table 1 The task data of every vessel

    表2 船舶時間數(shù)據(jù)Table 2 The time data of every vessel

    圖5 兩艘船舶的貝位-時間甘特圖Fig.5 Bay & time Gantt chart of double vessels

    圖6 兩艘船舶的岸橋-時間甘特圖Fig.6 Quay crane & time Gantt chart of double vessels

    4.2 5艘船舶的輸入數(shù)據(jù)與結(jié)果分析

    船舶數(shù)量為5,可用的岸橋數(shù)量為7,安全距離為1,其它輸入數(shù)據(jù)見表3和表4。計算結(jié)果如圖7所示,可以看出該結(jié)果既滿足模型中船舶的約束也滿足岸橋的約束,且能夠使每艘船舶在最晚離港時間之前盡早離港。圖8是將結(jié)果用岸橋和時間的坐標顯示,坐標里的數(shù)字代表貝位編號。

    表3 船舶任務(wù)量數(shù)據(jù)Table 3 The task data of every vessel

    表4 船舶時間數(shù)據(jù)Table 4 The time data of every vessel

    圖7 5艘船舶的貝位-時間甘特圖Fig.7 Bay & time Gantt chart of five vessels

    圖8 5艘船舶的岸橋-時間甘特圖Fig.8 Quay crane & time Gantt chart of five vessels

    4.3 不同規(guī)模的算例結(jié)果比較

    為了驗證筆者設(shè)計的遺傳算法求解該模型的有效性,使用C#_2010編程并調(diào)用Cplex求解作為對比。問題規(guī)模從小規(guī)模擴展到實際問題的規(guī)模,其中每艘船舶的貝位數(shù)以及每個貝位上的任務(wù)量是隨機生成的。船舶的時間信息和可用的岸橋數(shù)量是根據(jù)問題的規(guī)模已知的。表5是不同算例規(guī)模的情形下兩種方法的計算結(jié)果。對比可以看出:從能否求解來講,在所有算例中遺傳算法均求得了最優(yōu)解,而當問題規(guī)模擴大到一定程度時Cplex無法求解;在計算時間方面,盡管在較小規(guī)模時遺傳算法求解的時間稍長,但隨著問題規(guī)模的擴大,其計算時間遠遠小于Cplex的求解時間,如圖9;在求解質(zhì)量方面,小規(guī)模問題時兩者的求解結(jié)果一致,但隨著問題規(guī)模的擴大,遺傳算法的目標函數(shù)值(OFV)越來越優(yōu)于Cplex。總之,遺傳算法能夠在合理的時間內(nèi)求得最優(yōu)解,相比Cplex求解時間較短,求解質(zhì)量較高,體現(xiàn)了其有效性,且對于該問題的大規(guī)模問題,遺傳算法更加有效。

    圖9 兩種方法求解時間對比Fig.9 Comparison of computing time by two methods

    表5 Cplex和GA求解不同規(guī)模的算例結(jié)果Table 5 The results of Cplex and GA by solving examples of different scales

    5 結(jié) 語

    對集裝箱港口的岸橋分配與調(diào)度問題進行了研究。考慮了船舶和岸橋的實際因素,以船舶在港時間最短為目標建立了多船岸橋動態(tài)分配與貝位調(diào)度混合整數(shù)規(guī)劃模型。設(shè)計了遺傳算法進行求解,結(jié)果表明該模型可以有效解決岸橋的分配與調(diào)度問題。最后,通過不同規(guī)模的算例,分別用Cplex與遺傳算法進行求解。結(jié)果表明遺傳算法的求解時間更短,求解質(zhì)量更高,隨著問題規(guī)模的擴大其優(yōu)越性更加明顯。因此筆者提出的優(yōu)化方法,可為解決大規(guī)模的多船情形下岸橋的統(tǒng)一優(yōu)化問題提供科學合理的解決方案和決策依據(jù)。需提出的是,由于泊位計劃影響岸橋的分配和調(diào)度,所以未來也可以考慮把泊位計劃加入QCASP中,進而更加有效的優(yōu)化集裝箱港口的岸邊操作。

    [1] KIM K H,PARK Y M.A crane scheduling method for port container terminals[J].EuropeanJournalofOperationalResearch,2004,156(3):752-768.

    [2] MOCCIA L,CORDEAU J F,GAUDIOSO M ,et al.A branch-and-cut algorithm for the quay crane scheduling problem in a container terminal[J].NavalResearchLogistics,2006,53(1):45-59.

    [3] SAMMARRA M,CORDEAU J F,LAPORTE G,et al.A tabu search heuristic for the quay crane scheduling problem[J].JournalofScheduling,2007,10(4/5):327-336.

    [4] BIERWIRTH C,MEISEL F.A fast heuristic for quay crane scheduling with interference constraints[J].JournalofScheduling,2009,12(4): 345-360.

    [5] 楊明珠.單船裝卸作業(yè)的岸橋調(diào)度[J].計算機工程與應(yīng)用,2011,47(10): 224-228. YANG Mingzhu.Quay crane scheduling method for single vessel[J].ComputerEngineeringandApplications,2011,47(10):224-228.

    [6] MEISEL F,BIERWIRTH C.A unified approach for the evaluation of quay crane scheduling models and algorithms[J].Computers&OperationsResearch,2011,38(3):683-693.

    [7] TAVAKKOLI-MOGHADDAM R,MAKUIA A,SALAHI S ,et al .An efficient algorithm for solving a new mathematical model for a quay crane scheduling problem in container ports[J].Computers&IndustrialEngineering,2009,56(1):241-248.

    [8] LIANG Chengji,HUANG Youfang,YANG Yang.A quay crane dynamic scheduling problem by hybrid evolutionary algorithm for berth allocation planning[J].Computers&IndustrialEngineering,2009,56(3):1021-1028.

    [9] ZHANG Canrong,ZHENG Li,ZHANG Zhihai,et al.The allocation of berths and quay cranes by using a sub-gradient optimization technique[J].Computers&IndustrialEngineering,2010,58(1):40-50.

    [10] LEE Der-horng,WANG Huiqiu.Integrated discrete berth allocation and quay crane scheduling in port container terminals[J].EngineeringOptimization,2010,42(8):747-761.

    [11] 趙坤強,韓曉龍,梁承姬.連續(xù)泊位下集裝箱港口泊位與橋吊協(xié)同調(diào)度優(yōu)化研究[J].武漢理工大學學報,2011,33(11):60-65. ZHAO Kunqiang,HAN Xiaolong,LIANG Chengji.Berth and quay crane collaborative optimization research at continuous terminals[J].JournalofWuhanUniversityofTechnology,2011,33(11):60-65.

    [12] IMAI A,CHEN H C,NISHIMURA E,et al.The simultaneous berth and quay crane allocation problem[J].TransportationResearchPartE:LogisticsandTransportationReview,2008,44(5):900-920.

    Quay Crane Dynamic Allocation and Scheduling with Bay Research at Container Terminal

    LIANG Chengji,WANG Dandan

    (Logistics Research Center, Shanghai Maritime University, Shanghai 201306, P. R. China)

    In order to improve the quay crane’s efficiency and service level for vessels at container terminal, aiming at the heavy workload in every bay of multi-vessels, by considering such factors as ship priority, quantity of available quay crane, quay crane collision and safe distance, a mixed integer programming model on dynamic allocation of quay crane for vessels and scheduling and deployment of bay was establizhed in an attempt to minimize vessels’ staying time at container terminal. Through different scales of examples, the results obtained from Cplex and genetic algorithm programmed by C# respectively showed that the proposed model can effectively solve the problems of quay crane allocation and scheduling of multi-vessels and genetic algorithm has shorter computing time and higher quality solution thus verifying the effectiveness of this model and algorithm.

    transportation engineering; quay crane dynamic allocation; quay crane scheduling with bay; quay crane collision; safe distance; genetic algorithm

    10.3969/j.issn.1674-0696.2016.02.34

    2014-10-27;

    2014-12-29

    國家自然科學基金項目(71471110,71301101);中國博士后科學基金項目(2014M550084)

    梁承姬(1970—),女(朝鮮族),吉林龍井人,教授,博士,主要從事集裝箱港口物流方面的研究。E-mail:liangcj@shmtu.edu.cn。

    王丹丹(1989—),女,山西長治人,碩士研究生,主要從事集裝箱港口設(shè)備資源優(yōu)化方面的研究。E-mail:783964457@qq.com。

    U694

    A

    1674-0696(2016)02-169-05

    猜你喜歡
    離港進港集裝箱
    美軍一架C-130J正在投放集裝箱
    軍事文摘(2023年5期)2023-03-27 09:13:10
    離港航班延誤成本研究
    價值工程(2022年15期)2022-07-13 05:37:08
    長三角機場群運行相關(guān)性分析
    大型滿載油輪使用魚山作業(yè)區(qū)南部進港航道航行方法探討
    航海(2020年6期)2020-12-23 06:53:37
    虛實之間——集裝箱衍生出的空間折疊
    船舶進靠浙能臺二電煤炭碼頭風險的研究
    珠江水運(2019年9期)2019-06-02 17:00:04
    旺角暴亂嫌犯被禁止離港
    我家住在集裝箱
    中國公路(2017年8期)2017-07-21 14:26:20
    一種新型自卸式污泥集裝箱罐
    專用汽車(2015年2期)2015-03-01 04:06:52
    寧波梅山港區(qū)進港航道將于今年建成
    水道港口(2014年1期)2014-04-27 14:14:35
    久久精品aⅴ一区二区三区四区| 午夜免费观看网址| 91成年电影在线观看| 久久久久国内视频| √禁漫天堂资源中文www| 欧美乱码精品一区二区三区| 免费在线观看影片大全网站| 久久婷婷人人爽人人干人人爱| 中文字幕熟女人妻在线| 国内精品久久久久久久电影| 啦啦啦观看免费观看视频高清| 一级黄色大片毛片| 两个人免费观看高清视频| 久久久久久国产a免费观看| 午夜影院日韩av| 欧美在线一区亚洲| 久久久久久久久久黄片| 精品人妻1区二区| ponron亚洲| 女人高潮潮喷娇喘18禁视频| av在线播放免费不卡| av在线播放免费不卡| 国产成年人精品一区二区| 亚洲欧美日韩无卡精品| 久久中文字幕一级| 国产精品亚洲av一区麻豆| 中出人妻视频一区二区| 日日爽夜夜爽网站| 999久久久国产精品视频| 国产高清视频在线观看网站| 别揉我奶头~嗯~啊~动态视频| 亚洲欧洲精品一区二区精品久久久| 亚洲欧美日韩东京热| av天堂在线播放| 久久香蕉国产精品| 久久这里只有精品19| 精品国产乱子伦一区二区三区| 国产亚洲av高清不卡| 日韩 欧美 亚洲 中文字幕| 免费看美女性在线毛片视频| 午夜老司机福利片| 黄色毛片三级朝国网站| 亚洲美女黄片视频| 最近最新中文字幕大全免费视频| 欧美日韩亚洲国产一区二区在线观看| 国产精品一区二区精品视频观看| 久久天堂一区二区三区四区| 18美女黄网站色大片免费观看| 欧美另类亚洲清纯唯美| 女同久久另类99精品国产91| 亚洲精品色激情综合| 国模一区二区三区四区视频 | 日韩欧美国产在线观看| 国产精品亚洲一级av第二区| 岛国视频午夜一区免费看| 免费看十八禁软件| 国产成人av激情在线播放| 亚洲成人免费电影在线观看| 久久精品影院6| 久久香蕉激情| 免费观看精品视频网站| 日本免费一区二区三区高清不卡| 极品教师在线免费播放| 熟女少妇亚洲综合色aaa.| 成人国产综合亚洲| 国产欧美日韩一区二区三| 亚洲精华国产精华精| 国产成人啪精品午夜网站| 午夜激情av网站| 婷婷六月久久综合丁香| 91国产中文字幕| 国产在线精品亚洲第一网站| 好男人电影高清在线观看| 国产一区二区激情短视频| 精品熟女少妇八av免费久了| 特级一级黄色大片| 日韩 欧美 亚洲 中文字幕| 大型av网站在线播放| 蜜桃久久精品国产亚洲av| 曰老女人黄片| 在线免费观看的www视频| 久久性视频一级片| 一级毛片高清免费大全| 丁香欧美五月| а√天堂www在线а√下载| 久久久久国产一级毛片高清牌| 成人国产综合亚洲| 久久天堂一区二区三区四区| 亚洲成人久久爱视频| 成人欧美大片| 日韩欧美精品v在线| 一个人免费在线观看电影 | 亚洲va日本ⅴa欧美va伊人久久| 巨乳人妻的诱惑在线观看| 亚洲国产精品成人综合色| 国内久久婷婷六月综合欲色啪| 精品少妇一区二区三区视频日本电影| 深夜精品福利| 国内精品久久久久久久电影| 午夜激情av网站| 亚洲激情在线av| 桃红色精品国产亚洲av| 看片在线看免费视频| 伊人久久大香线蕉亚洲五| 午夜福利免费观看在线| 91麻豆av在线| 欧美日韩黄片免| 久久香蕉国产精品| 人妻久久中文字幕网| 亚洲一区中文字幕在线| 日韩欧美一区二区三区在线观看| 在线观看日韩欧美| 看黄色毛片网站| 在线a可以看的网站| 亚洲人成网站高清观看| 一本久久中文字幕| 久久久精品国产亚洲av高清涩受| 亚洲avbb在线观看| 久久久久久亚洲精品国产蜜桃av| 午夜精品久久久久久毛片777| 99热6这里只有精品| 桃色一区二区三区在线观看| 久久久久精品国产欧美久久久| 夜夜躁狠狠躁天天躁| 99久久99久久久精品蜜桃| 99国产精品一区二区三区| 动漫黄色视频在线观看| 亚洲精品一区av在线观看| 宅男免费午夜| 亚洲av第一区精品v没综合| 香蕉国产在线看| 亚洲国产精品久久男人天堂| 免费一级毛片在线播放高清视频| 亚洲av五月六月丁香网| 人妻夜夜爽99麻豆av| 日韩三级视频一区二区三区| 亚洲第一欧美日韩一区二区三区| 久久精品91蜜桃| 久久久久久国产a免费观看| 最新在线观看一区二区三区| 国产乱人伦免费视频| 丰满人妻一区二区三区视频av | 99re在线观看精品视频| 亚洲一卡2卡3卡4卡5卡精品中文| 在线十欧美十亚洲十日本专区| 国产蜜桃级精品一区二区三区| 欧美成狂野欧美在线观看| 中文字幕av在线有码专区| 男女午夜视频在线观看| 1024视频免费在线观看| 真人做人爱边吃奶动态| 国产真人三级小视频在线观看| 俺也久久电影网| 看片在线看免费视频| 两性夫妻黄色片| 中出人妻视频一区二区| 在线观看www视频免费| 色综合站精品国产| 久久伊人香网站| 999久久久精品免费观看国产| a级毛片a级免费在线| 性欧美人与动物交配| 一二三四在线观看免费中文在| 深夜精品福利| 很黄的视频免费| 亚洲av五月六月丁香网| 久久久国产精品麻豆| 色av中文字幕| 黄色 视频免费看| 天天添夜夜摸| 级片在线观看| 一区二区三区高清视频在线| 国产精品久久久久久人妻精品电影| 伦理电影免费视频| 午夜视频精品福利| 国产成人精品无人区| 午夜久久久久精精品| 亚洲最大成人中文| 91麻豆精品激情在线观看国产| 在线观看日韩欧美| 国产黄片美女视频| 丝袜美腿诱惑在线| 日韩大尺度精品在线看网址| 男女下面进入的视频免费午夜| 亚洲欧美日韩高清专用| 老司机深夜福利视频在线观看| 丰满的人妻完整版| 久久久国产成人精品二区| 国内精品一区二区在线观看| 久久中文看片网| 欧美zozozo另类| 国产亚洲av嫩草精品影院| 欧美国产日韩亚洲一区| 好看av亚洲va欧美ⅴa在| 日本黄大片高清| 欧美在线黄色| 久久精品综合一区二区三区| 国产精华一区二区三区| 久久草成人影院| 亚洲中文日韩欧美视频| 一进一出抽搐gif免费好疼| av欧美777| 大型av网站在线播放| 亚洲一区高清亚洲精品| 中文字幕高清在线视频| 身体一侧抽搐| 91麻豆精品激情在线观看国产| 国产精品 国内视频| 99久久精品国产亚洲精品| 欧美乱码精品一区二区三区| 欧美+亚洲+日韩+国产| 麻豆国产97在线/欧美 | 亚洲精品在线美女| 亚洲九九香蕉| 亚洲av成人精品一区久久| 免费看a级黄色片| 国产精品久久视频播放| 久久精品亚洲精品国产色婷小说| 床上黄色一级片| 亚洲男人的天堂狠狠| 日韩欧美 国产精品| 日韩av在线大香蕉| 久久久久免费精品人妻一区二区| 美女黄网站色视频| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美国产一区二区入口| 真人一进一出gif抽搐免费| 欧美日韩瑟瑟在线播放| 国产一区二区三区在线臀色熟女| 一级毛片女人18水好多| 免费av毛片视频| 老司机午夜十八禁免费视频| 12—13女人毛片做爰片一| 一个人观看的视频www高清免费观看 | 亚洲国产精品成人综合色| 天堂av国产一区二区熟女人妻 | 99久久综合精品五月天人人| 97碰自拍视频| 欧美日韩一级在线毛片| 国产精品久久久久久久电影 | 色综合站精品国产| 亚洲欧美日韩高清在线视频| 国产91精品成人一区二区三区| 桃色一区二区三区在线观看| 伊人久久大香线蕉亚洲五| 18禁裸乳无遮挡免费网站照片| 韩国av一区二区三区四区| 中文字幕人妻丝袜一区二区| 日韩欧美在线二视频| 亚洲精品一区av在线观看| 成年人黄色毛片网站| 又黄又粗又硬又大视频| 欧美性猛交黑人性爽| 91大片在线观看| 五月伊人婷婷丁香| 性欧美人与动物交配| 国产爱豆传媒在线观看 | 国产爱豆传媒在线观看 | 久久欧美精品欧美久久欧美| 精品熟女少妇八av免费久了| 不卡一级毛片| 嫁个100分男人电影在线观看| 少妇熟女aⅴ在线视频| 九九热线精品视视频播放| 成人18禁在线播放| 久久精品夜夜夜夜夜久久蜜豆 | 国产99白浆流出| 国内毛片毛片毛片毛片毛片| 一区二区三区国产精品乱码| 99热只有精品国产| 欧美国产日韩亚洲一区| 中文字幕人成人乱码亚洲影| 99久久国产精品久久久| 欧美三级亚洲精品| 精品久久久久久久末码| 国产精品国产高清国产av| 国产三级黄色录像| 欧美久久黑人一区二区| 久久婷婷人人爽人人干人人爱| 两人在一起打扑克的视频| 亚洲中文字幕一区二区三区有码在线看 | 色播亚洲综合网| 亚洲av片天天在线观看| 国产一区二区在线av高清观看| 亚洲男人的天堂狠狠| 精品人妻1区二区| tocl精华| 夜夜看夜夜爽夜夜摸| 亚洲一码二码三码区别大吗| 久久久久免费精品人妻一区二区| 久久久久久国产a免费观看| 91大片在线观看| 国产一区二区在线观看日韩 | 俄罗斯特黄特色一大片| 国产激情久久老熟女| 窝窝影院91人妻| x7x7x7水蜜桃| 怎么达到女性高潮| 亚洲免费av在线视频| 久久久国产成人精品二区| 午夜日韩欧美国产| 99国产极品粉嫩在线观看| 免费看美女性在线毛片视频| 99久久精品国产亚洲精品| 国产亚洲精品av在线| 免费看a级黄色片| 嫩草影院精品99| 欧美3d第一页| 亚洲欧美日韩东京热| 99re在线观看精品视频| 久久天躁狠狠躁夜夜2o2o| 国产高清videossex| 在线观看www视频免费| 19禁男女啪啪无遮挡网站| 一进一出抽搐动态| 特级一级黄色大片| 在线永久观看黄色视频| 真人一进一出gif抽搐免费| 国产精品1区2区在线观看.| 午夜亚洲福利在线播放| 亚洲专区字幕在线| 色综合站精品国产| 真人做人爱边吃奶动态| 国产精品日韩av在线免费观看| 男女那种视频在线观看| 欧美黄色片欧美黄色片| 亚洲在线自拍视频| 国产成年人精品一区二区| 久久伊人香网站| 男人舔女人下体高潮全视频| 婷婷丁香在线五月| 欧美在线黄色| www.精华液| av在线播放免费不卡| 国产v大片淫在线免费观看| av天堂在线播放| 非洲黑人性xxxx精品又粗又长| 老熟妇仑乱视频hdxx| 好男人在线观看高清免费视频| 亚洲电影在线观看av| 色精品久久人妻99蜜桃| 亚洲人成网站在线播放欧美日韩| 男人的好看免费观看在线视频 | 99久久无色码亚洲精品果冻| 婷婷精品国产亚洲av| 亚洲黑人精品在线| 国产黄片美女视频| cao死你这个sao货| 欧美不卡视频在线免费观看 | 搞女人的毛片| 国产在线精品亚洲第一网站| 日韩精品中文字幕看吧| 亚洲熟妇熟女久久| 亚洲五月婷婷丁香| 欧美 亚洲 国产 日韩一| 免费高清视频大片| 国产v大片淫在线免费观看| 婷婷精品国产亚洲av| 久久久水蜜桃国产精品网| 欧美一级毛片孕妇| 欧美日韩黄片免| 69av精品久久久久久| 午夜成年电影在线免费观看| ponron亚洲| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品综合一区在线观看 | 一级作爱视频免费观看| 免费在线观看视频国产中文字幕亚洲| 中文字幕人妻丝袜一区二区| 亚洲人与动物交配视频| 中文资源天堂在线| 欧美黑人精品巨大| 欧美色视频一区免费| 免费在线观看视频国产中文字幕亚洲| 宅男免费午夜| 亚洲人成电影免费在线| 天天一区二区日本电影三级| 国产精品av视频在线免费观看| 国产成人一区二区三区免费视频网站| www日本在线高清视频| 欧美黄色片欧美黄色片| 欧美色欧美亚洲另类二区| 老司机在亚洲福利影院| 色综合站精品国产| 国产精品永久免费网站| 欧美一区二区精品小视频在线| 欧美在线一区亚洲| 亚洲国产欧美网| 黄色女人牲交| 蜜桃久久精品国产亚洲av| 草草在线视频免费看| 黄色丝袜av网址大全| 国产精品久久久久久精品电影| 天堂影院成人在线观看| 亚洲男人天堂网一区| 人妻丰满熟妇av一区二区三区| 亚洲熟妇中文字幕五十中出| 波多野结衣巨乳人妻| 嫩草影院精品99| 一本大道久久a久久精品| 国产高清有码在线观看视频 | 国产精品爽爽va在线观看网站| 一级毛片高清免费大全| 制服丝袜大香蕉在线| 欧美精品啪啪一区二区三区| 欧美成人一区二区免费高清观看 | 亚洲男人天堂网一区| 啦啦啦观看免费观看视频高清| 国产一区二区在线av高清观看| 亚洲中文日韩欧美视频| 精品免费久久久久久久清纯| 国产精品98久久久久久宅男小说| 青草久久国产| 国产三级中文精品| 国产精品亚洲av一区麻豆| 中文字幕人成人乱码亚洲影| 中文字幕最新亚洲高清| 床上黄色一级片| 亚洲国产欧美网| 国产亚洲精品久久久久5区| 免费在线观看视频国产中文字幕亚洲| 中亚洲国语对白在线视频| 夜夜看夜夜爽夜夜摸| 99在线视频只有这里精品首页| 精品电影一区二区在线| 国产单亲对白刺激| 国产精品久久久人人做人人爽| 男女那种视频在线观看| 九九热线精品视视频播放| 亚洲一区二区三区不卡视频| 老鸭窝网址在线观看| 熟女电影av网| 三级毛片av免费| 在线视频色国产色| 级片在线观看| 亚洲中文av在线| 国产伦一二天堂av在线观看| 日本 欧美在线| 精品国产亚洲在线| 很黄的视频免费| 一级a爱片免费观看的视频| 色尼玛亚洲综合影院| 国产精品免费一区二区三区在线| 真人做人爱边吃奶动态| 天天添夜夜摸| 欧美精品亚洲一区二区| 青草久久国产| 国产av在哪里看| 精品日产1卡2卡| 午夜免费观看网址| 波多野结衣巨乳人妻| 亚洲成人久久性| 热99re8久久精品国产| 亚洲成人免费电影在线观看| 色精品久久人妻99蜜桃| 久久亚洲真实| 男女那种视频在线观看| avwww免费| 夜夜夜夜夜久久久久| 黄色 视频免费看| a在线观看视频网站| 又黄又爽又免费观看的视频| 中文字幕高清在线视频| 波多野结衣高清作品| 91大片在线观看| 少妇被粗大的猛进出69影院| 少妇的丰满在线观看| 欧美成人午夜精品| www.自偷自拍.com| 亚洲性夜色夜夜综合| 欧美久久黑人一区二区| 国内揄拍国产精品人妻在线| 国产区一区二久久| 亚洲精品在线观看二区| 欧美色欧美亚洲另类二区| 999久久久国产精品视频| 在线观看免费午夜福利视频| 午夜两性在线视频| 国产成人精品久久二区二区91| 在线国产一区二区在线| 国产欧美日韩精品亚洲av| 亚洲国产日韩欧美精品在线观看 | 亚洲乱码一区二区免费版| 欧美中文综合在线视频| 国产高清视频在线观看网站| 成人永久免费在线观看视频| 日韩大码丰满熟妇| 亚洲av成人不卡在线观看播放网| 色av中文字幕| 老汉色av国产亚洲站长工具| 国产精品亚洲美女久久久| 国模一区二区三区四区视频 | 久久久精品欧美日韩精品| 亚洲国产精品合色在线| 香蕉国产在线看| 亚洲乱码一区二区免费版| 欧美中文综合在线视频| 国产一区二区激情短视频| 无限看片的www在线观看| 一本久久中文字幕| 男女午夜视频在线观看| 亚洲天堂国产精品一区在线| 天天添夜夜摸| 女同久久另类99精品国产91| 免费在线观看成人毛片| 日韩欧美一区二区三区在线观看| 一区二区三区国产精品乱码| 国产精品 国内视频| 久久久久久九九精品二区国产 | www.999成人在线观看| 91在线观看av| 国产成人影院久久av| 亚洲全国av大片| 精品久久久久久成人av| 国产一区二区激情短视频| 亚洲一区中文字幕在线| 在线观看日韩欧美| 国产精品香港三级国产av潘金莲| 99精品在免费线老司机午夜| 国产成人精品久久二区二区91| 999精品在线视频| 久久久久久亚洲精品国产蜜桃av| 日韩欧美三级三区| 亚洲自偷自拍图片 自拍| 十八禁人妻一区二区| 日韩精品青青久久久久久| 亚洲国产精品久久男人天堂| 久久久国产精品麻豆| 国产精品1区2区在线观看.| 黄色片一级片一级黄色片| 亚洲国产看品久久| 日日爽夜夜爽网站| 人成视频在线观看免费观看| 大型黄色视频在线免费观看| 国产97色在线日韩免费| 男女午夜视频在线观看| 国产精品,欧美在线| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美三级亚洲精品| 国产av又大| 久久香蕉精品热| 嫁个100分男人电影在线观看| 成人特级黄色片久久久久久久| 免费看十八禁软件| 国产精品久久久久久久电影 | 亚洲中文字幕一区二区三区有码在线看 | 午夜福利在线观看吧| 欧美乱码精品一区二区三区| 最近视频中文字幕2019在线8| 久久久久久久久免费视频了| 深夜精品福利| 淫妇啪啪啪对白视频| 露出奶头的视频| 香蕉久久夜色| 欧美在线一区亚洲| 亚洲电影在线观看av| 一进一出抽搐动态| 精品久久久久久,| 久久精品国产99精品国产亚洲性色| 巨乳人妻的诱惑在线观看| 久久伊人香网站| 亚洲,欧美精品.| 成人欧美大片| av有码第一页| 人人妻人人看人人澡| 亚洲av第一区精品v没综合| 老司机在亚洲福利影院| 久久久久九九精品影院| 一区二区三区高清视频在线| 99久久精品国产亚洲精品| 精华霜和精华液先用哪个| 在线播放国产精品三级| 99精品欧美一区二区三区四区| 国产1区2区3区精品| 一个人免费在线观看电影 | 亚洲一区二区三区不卡视频| 国产一区二区三区在线臀色熟女| 9191精品国产免费久久| 淫秽高清视频在线观看| 全区人妻精品视频| 国产成人欧美在线观看| 亚洲精品粉嫩美女一区| 国产精品一区二区精品视频观看| av片东京热男人的天堂| 精品乱码久久久久久99久播| 欧美日韩瑟瑟在线播放| 午夜老司机福利片| 免费在线观看完整版高清| 少妇熟女aⅴ在线视频| 日本一区二区免费在线视频| 亚洲激情在线av| 日本熟妇午夜| 午夜老司机福利片| 亚洲男人天堂网一区| 欧美不卡视频在线免费观看 | 久久99热这里只有精品18| 欧美日韩国产亚洲二区| 别揉我奶头~嗯~啊~动态视频| 日韩欧美一区二区三区在线观看| 欧美乱色亚洲激情| 欧美人与性动交α欧美精品济南到| 久久精品国产亚洲av香蕉五月| 亚洲一区二区三区不卡视频| 别揉我奶头~嗯~啊~动态视频| 国产蜜桃级精品一区二区三区| 无限看片的www在线观看| 欧美日韩福利视频一区二区| 国产精品,欧美在线| 岛国在线观看网站| aaaaa片日本免费| 天堂√8在线中文| 一级黄色大片毛片| 99久久综合精品五月天人人| 一个人免费在线观看的高清视频| 999久久久精品免费观看国产|