• <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
    av欧美777| 亚洲人成网站在线播放欧美日韩| 国产成+人综合+亚洲专区| 亚洲av成人一区二区三| 脱女人内裤的视频| 亚洲欧美日韩高清在线视频| 最新在线观看一区二区三区| 久99久视频精品免费| 别揉我奶头~嗯~啊~动态视频| 国产亚洲av嫩草精品影院| 男插女下体视频免费在线播放| 欧美成人性av电影在线观看| 欧美一区二区国产精品久久精品 | 亚洲精品久久成人aⅴ小说| 欧美在线黄色| 午夜a级毛片| 精品久久久久久,| 日本五十路高清| 亚洲专区中文字幕在线| 久久欧美精品欧美久久欧美| 久久久国产欧美日韩av| 国产精品久久久av美女十八| www日本黄色视频网| 色噜噜av男人的天堂激情| 天堂动漫精品| 在线观看一区二区三区| 国产私拍福利视频在线观看| 久久天堂一区二区三区四区| АⅤ资源中文在线天堂| a级毛片a级免费在线| 天天添夜夜摸| 国产成人精品久久二区二区91| 18禁国产床啪视频网站| 成人国产综合亚洲| 男女那种视频在线观看| 久久久久久亚洲精品国产蜜桃av| 最近在线观看免费完整版| 一夜夜www| 亚洲精品粉嫩美女一区| www日本黄色视频网| 熟女电影av网| 亚洲午夜理论影院| 午夜成年电影在线免费观看| av福利片在线| 免费人成视频x8x8入口观看| 亚洲性夜色夜夜综合| 天天添夜夜摸| 亚洲熟妇中文字幕五十中出| 亚洲欧美日韩无卡精品| 日韩免费av在线播放| 国产精品亚洲一级av第二区| 国产亚洲精品一区二区www| 国产又黄又爽又无遮挡在线| 久久精品国产99精品国产亚洲性色| 嫁个100分男人电影在线观看| 99精品欧美一区二区三区四区| 最好的美女福利视频网| 久久国产精品影院| 韩国av一区二区三区四区| 国内精品久久久久精免费| 亚洲 欧美一区二区三区| 亚洲av电影不卡..在线观看| 美女 人体艺术 gogo| 精品第一国产精品| 嫁个100分男人电影在线观看| 亚洲av日韩精品久久久久久密| 色综合亚洲欧美另类图片| 神马国产精品三级电影在线观看 | www.999成人在线观看| 一a级毛片在线观看| 女人高潮潮喷娇喘18禁视频| 欧美日韩黄片免| 757午夜福利合集在线观看| 免费观看精品视频网站| 日韩欧美在线二视频| 亚洲人与动物交配视频| 午夜福利在线观看吧| 舔av片在线| 色综合亚洲欧美另类图片| 国产亚洲精品第一综合不卡| 香蕉久久夜色| 免费看日本二区| 久久亚洲精品不卡| 欧美+亚洲+日韩+国产| 中文字幕最新亚洲高清| 久久久国产成人免费| tocl精华| 精品免费久久久久久久清纯| 国产私拍福利视频在线观看| 亚洲欧美日韩无卡精品| 午夜两性在线视频| 亚洲一区二区三区不卡视频| 精品日产1卡2卡| 啦啦啦免费观看视频1| 一级黄色大片毛片| 欧美高清成人免费视频www| 午夜成年电影在线免费观看| 麻豆国产av国片精品| 国产一区二区激情短视频| 国产不卡一卡二| 欧美一区二区精品小视频在线| 国产免费男女视频| 一级毛片精品| 亚洲av电影不卡..在线观看| 久久香蕉激情| 两性夫妻黄色片| 久久久久久大精品| av国产免费在线观看| 免费电影在线观看免费观看| 国产精品爽爽va在线观看网站| 国产亚洲精品第一综合不卡| ponron亚洲| 国产高清视频在线播放一区| 特级一级黄色大片| 亚洲一区高清亚洲精品| 别揉我奶头~嗯~啊~动态视频| 国产真实乱freesex| 国产熟女午夜一区二区三区| 久久久久九九精品影院| 日韩大码丰满熟妇| 两性夫妻黄色片| 亚洲av美国av| 成人精品一区二区免费| 国产av麻豆久久久久久久| 亚洲av成人精品一区久久| 91大片在线观看| 国产在线精品亚洲第一网站| 国产aⅴ精品一区二区三区波| 国产人伦9x9x在线观看| 精品久久久久久久久久久久久| 精品一区二区三区四区五区乱码| 1024视频免费在线观看| 日本a在线网址| 一个人免费在线观看电影 | 777久久人妻少妇嫩草av网站| 亚洲成av人片在线播放无| 亚洲中文日韩欧美视频| xxx96com| a级毛片在线看网站| 香蕉丝袜av| 亚洲成av人片在线播放无| 久久99热这里只有精品18| 日韩欧美精品v在线| 亚洲真实伦在线观看| 伊人久久大香线蕉亚洲五| 日韩中文字幕欧美一区二区| 小说图片视频综合网站| 夜夜看夜夜爽夜夜摸| 亚洲中文字幕一区二区三区有码在线看 | 久久久久久久久免费视频了| 欧洲精品卡2卡3卡4卡5卡区| 一卡2卡三卡四卡精品乱码亚洲| 精品日产1卡2卡| 性欧美人与动物交配| 丰满人妻熟妇乱又伦精品不卡| 欧美一级毛片孕妇| 亚洲avbb在线观看| 正在播放国产对白刺激| 日韩大尺度精品在线看网址| 中文字幕av在线有码专区| 波多野结衣巨乳人妻| 中文字幕av在线有码专区| 色综合婷婷激情| 91九色精品人成在线观看| 亚洲国产精品合色在线| 精品福利观看| 精品国产亚洲在线| ponron亚洲| 美女大奶头视频| 欧美黑人精品巨大| 99热只有精品国产| 一二三四社区在线视频社区8| 国产精品日韩av在线免费观看| 国产精品av久久久久免费| 岛国在线免费视频观看| 久久精品91无色码中文字幕| 成人手机av| 欧美色视频一区免费| 午夜激情福利司机影院| 中出人妻视频一区二区| 级片在线观看| 国产av麻豆久久久久久久| 亚洲最大成人中文| 国产午夜福利久久久久久| 天天躁夜夜躁狠狠躁躁| 村上凉子中文字幕在线| 国产成人系列免费观看| 国产精品久久久人人做人人爽| 国产欧美日韩一区二区三| 老汉色av国产亚洲站长工具| 亚洲男人天堂网一区| 久久久水蜜桃国产精品网| 久久国产乱子伦精品免费另类| 人人妻人人澡欧美一区二区| 国产精品乱码一区二三区的特点| av中文乱码字幕在线| 19禁男女啪啪无遮挡网站| 99在线人妻在线中文字幕| 午夜福利免费观看在线| 香蕉丝袜av| 色精品久久人妻99蜜桃| 亚洲五月婷婷丁香| 久9热在线精品视频| 亚洲成a人片在线一区二区| 亚洲第一欧美日韩一区二区三区| 黄色丝袜av网址大全| 99精品在免费线老司机午夜| 欧美乱色亚洲激情| 亚洲人成网站在线播放欧美日韩| 好男人电影高清在线观看| 啦啦啦免费观看视频1| 五月玫瑰六月丁香| 最近视频中文字幕2019在线8| cao死你这个sao货| 99riav亚洲国产免费| 国产真实乱freesex| 男女下面进入的视频免费午夜| 天堂影院成人在线观看| 此物有八面人人有两片| 香蕉av资源在线| 成人av在线播放网站| 欧美成人一区二区免费高清观看 | 少妇裸体淫交视频免费看高清 | 欧美色欧美亚洲另类二区| √禁漫天堂资源中文www| 99久久99久久久精品蜜桃| 亚洲人成网站在线播放欧美日韩| 国产av在哪里看| 国产亚洲精品一区二区www| 91麻豆av在线| 国产片内射在线| 人妻丰满熟妇av一区二区三区| 久久精品亚洲精品国产色婷小说| 午夜激情av网站| 国产午夜福利久久久久久| 国产一区二区在线观看日韩 | 欧美黑人巨大hd| 三级毛片av免费| 国产高清videossex| 麻豆一二三区av精品| 久久久久免费精品人妻一区二区| 在线观看66精品国产| 亚洲av成人不卡在线观看播放网| 亚洲精品久久成人aⅴ小说| 国产亚洲欧美在线一区二区| 最近最新中文字幕大全免费视频| 无人区码免费观看不卡| 亚洲国产中文字幕在线视频| 一级毛片高清免费大全| 女同久久另类99精品国产91| 男女午夜视频在线观看| 国产精品一区二区三区四区久久| 精品国产乱码久久久久久男人| 亚洲欧美日韩高清专用| 欧美中文综合在线视频| 国产成人精品久久二区二区免费| 国产精品一区二区精品视频观看| 日韩免费av在线播放| 免费av毛片视频| 91成年电影在线观看| 久久久久国内视频| 精品欧美一区二区三区在线| 人妻丰满熟妇av一区二区三区| 少妇熟女aⅴ在线视频| 此物有八面人人有两片| 在线免费观看的www视频| 久久精品国产亚洲av高清一级| 三级毛片av免费| 久久久精品大字幕| 婷婷精品国产亚洲av| 国产成人系列免费观看| 88av欧美| 母亲3免费完整高清在线观看| 中文字幕精品亚洲无线码一区| 国产又黄又爽又无遮挡在线| 国产精品一区二区三区四区免费观看 | 中文亚洲av片在线观看爽| 不卡一级毛片| 草草在线视频免费看| 叶爱在线成人免费视频播放| 黄片大片在线免费观看| 特大巨黑吊av在线直播| 亚洲av成人一区二区三| 日韩中文字幕欧美一区二区| 一a级毛片在线观看| 90打野战视频偷拍视频| 91国产中文字幕| 欧美中文综合在线视频| 91在线观看av| 国产精品香港三级国产av潘金莲| 免费人成视频x8x8入口观看| 桃色一区二区三区在线观看| 啦啦啦韩国在线观看视频| 三级男女做爰猛烈吃奶摸视频| 无遮挡黄片免费观看| 国产一区二区在线av高清观看| 18禁裸乳无遮挡免费网站照片| 美女扒开内裤让男人捅视频| 中文资源天堂在线| 精品电影一区二区在线| 国产精品,欧美在线| av视频在线观看入口| e午夜精品久久久久久久| 成人精品一区二区免费| 欧美日韩黄片免| 宅男免费午夜| 日本 欧美在线| 色av中文字幕| 又黄又爽又免费观看的视频| 亚洲成人精品中文字幕电影| 亚洲国产欧美人成| 亚洲色图 男人天堂 中文字幕| 正在播放国产对白刺激| 免费人成视频x8x8入口观看| 精品无人区乱码1区二区| 免费在线观看日本一区| 两个人的视频大全免费| 久久精品夜夜夜夜夜久久蜜豆 | 麻豆国产av国片精品| 国产高清videossex| 老熟妇乱子伦视频在线观看| 国产成人系列免费观看| 成在线人永久免费视频| 国产麻豆成人av免费视频| 久久久精品大字幕| 女同久久另类99精品国产91| 男人舔女人下体高潮全视频| 国内毛片毛片毛片毛片毛片| 狂野欧美白嫩少妇大欣赏| 国内精品久久久久精免费| 麻豆成人av在线观看| 巨乳人妻的诱惑在线观看| 制服丝袜大香蕉在线| 在线看三级毛片| 在线观看舔阴道视频| 又粗又爽又猛毛片免费看| 日韩三级视频一区二区三区| 欧美色视频一区免费| 亚洲av成人精品一区久久| 又紧又爽又黄一区二区| 色综合亚洲欧美另类图片| 久久天躁狠狠躁夜夜2o2o| 亚洲国产中文字幕在线视频| 国内精品一区二区在线观看| 亚洲乱码一区二区免费版| 成在线人永久免费视频| 国产区一区二久久| 搞女人的毛片| 亚洲国产精品999在线| 九九热线精品视视频播放| 无人区码免费观看不卡| 成人一区二区视频在线观看| 老司机靠b影院| 天堂√8在线中文| aaaaa片日本免费| 人人妻人人看人人澡| 老汉色∧v一级毛片| 日韩三级视频一区二区三区| 麻豆国产av国片精品| 亚洲人成网站高清观看| 久久欧美精品欧美久久欧美| 久久欧美精品欧美久久欧美| 欧美性猛交黑人性爽| 操出白浆在线播放| 丁香六月欧美| 日本免费一区二区三区高清不卡| 久久久水蜜桃国产精品网| 精品熟女少妇八av免费久了| 日日夜夜操网爽| 午夜福利高清视频| 亚洲精品粉嫩美女一区| 狂野欧美白嫩少妇大欣赏| 91在线观看av| 日韩精品青青久久久久久| 91成年电影在线观看| 在线看三级毛片| 岛国在线观看网站| 国产免费av片在线观看野外av| 免费搜索国产男女视频| 蜜桃久久精品国产亚洲av| 国产乱人伦免费视频| 草草在线视频免费看| 亚洲五月婷婷丁香| 夜夜躁狠狠躁天天躁| 老汉色∧v一级毛片| 久久久久久国产a免费观看| 久久这里只有精品19| 亚洲欧美日韩东京热| 在线视频色国产色| 亚洲七黄色美女视频| 国产单亲对白刺激| 日本熟妇午夜| 亚洲片人在线观看| 精品久久久久久久毛片微露脸| 欧美日韩一级在线毛片| 午夜日韩欧美国产| 久久久国产欧美日韩av| 日韩欧美免费精品| 很黄的视频免费| 国产免费男女视频| 国产高清激情床上av| 亚洲aⅴ乱码一区二区在线播放 | 老司机靠b影院| 床上黄色一级片| 97超级碰碰碰精品色视频在线观看| 91在线观看av| 村上凉子中文字幕在线| 99久久无色码亚洲精品果冻| 成人18禁高潮啪啪吃奶动态图| 精品国产亚洲在线| 久久久久久国产a免费观看| 欧美乱妇无乱码| 看免费av毛片| 久久人人精品亚洲av| 在线观看午夜福利视频| 啦啦啦韩国在线观看视频| av国产免费在线观看| 黄频高清免费视频| 国产视频内射| 国产亚洲精品一区二区www| 亚洲乱码一区二区免费版| 色播亚洲综合网| 天堂影院成人在线观看| 熟女电影av网| 亚洲,欧美精品.| 久久久精品国产亚洲av高清涩受| 女人爽到高潮嗷嗷叫在线视频| 国产精品av视频在线免费观看| 日韩大码丰满熟妇| 免费看日本二区| 亚洲专区国产一区二区| 日韩欧美精品v在线| 51午夜福利影视在线观看| 久久伊人香网站| 老司机靠b影院| 欧美性猛交黑人性爽| 国产亚洲欧美98| 国产精品1区2区在线观看.| 久久精品影院6| xxxwww97欧美| 欧美最黄视频在线播放免费| 亚洲欧美日韩无卡精品| 十八禁网站免费在线| 校园春色视频在线观看| 久久久久亚洲av毛片大全| 最近在线观看免费完整版| 国产黄片美女视频| 亚洲九九香蕉| 91老司机精品| 一a级毛片在线观看| 女人被狂操c到高潮| 人人妻,人人澡人人爽秒播| 麻豆国产av国片精品| 国产亚洲精品久久久久久毛片| 国产高清有码在线观看视频 | 蜜桃久久精品国产亚洲av| 婷婷六月久久综合丁香| 中文字幕久久专区| 香蕉丝袜av| 国产精品一区二区精品视频观看| 好男人在线观看高清免费视频| 欧美又色又爽又黄视频| 成人亚洲精品av一区二区| 国产精品免费一区二区三区在线| 一级毛片高清免费大全| 国产成人精品久久二区二区91| 中文亚洲av片在线观看爽| 国产熟女午夜一区二区三区| 日韩成人在线观看一区二区三区| 国产人伦9x9x在线观看| 91成年电影在线观看| 757午夜福利合集在线观看| 又大又爽又粗| 日韩欧美国产一区二区入口| 国产激情久久老熟女| 亚洲熟妇中文字幕五十中出| 国产伦一二天堂av在线观看| a在线观看视频网站| 国产高清视频在线播放一区| 国产v大片淫在线免费观看| 色综合亚洲欧美另类图片| 怎么达到女性高潮| 一个人免费在线观看电影 | 18禁国产床啪视频网站| 精品久久久久久久久久久久久| 亚洲熟妇熟女久久| 亚洲最大成人中文| 男女视频在线观看网站免费 | 成年女人毛片免费观看观看9| 国产精品久久久人人做人人爽| 日日爽夜夜爽网站| 男女之事视频高清在线观看| 精品国产亚洲在线| 欧美丝袜亚洲另类 | 中国美女看黄片| 一进一出抽搐gif免费好疼| 国产单亲对白刺激| 国产91精品成人一区二区三区| 精品高清国产在线一区| 母亲3免费完整高清在线观看| 成人av一区二区三区在线看| 国产成人欧美在线观看| 美女黄网站色视频| 好看av亚洲va欧美ⅴa在| 欧美成人一区二区免费高清观看 | 91成年电影在线观看| 精品福利观看| bbb黄色大片| 午夜a级毛片| 啦啦啦韩国在线观看视频| 男女之事视频高清在线观看| 久久久久久免费高清国产稀缺| 夜夜夜夜夜久久久久| 欧美午夜高清在线| 韩国av一区二区三区四区| 一边摸一边抽搐一进一小说| 亚洲精品在线观看二区| 中文字幕最新亚洲高清| 91成年电影在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产成人aa在线观看| 午夜a级毛片| 日日摸夜夜添夜夜添小说| 人成视频在线观看免费观看| 99久久久亚洲精品蜜臀av| 国产v大片淫在线免费观看| 日本精品一区二区三区蜜桃| 日韩国内少妇激情av| 一本一本综合久久| 巨乳人妻的诱惑在线观看| 人妻丰满熟妇av一区二区三区| 搡老妇女老女人老熟妇| 人妻夜夜爽99麻豆av| 91麻豆av在线| 黄色视频不卡| 国产黄片美女视频| 中文字幕精品亚洲无线码一区| 两个人免费观看高清视频| 人人妻人人看人人澡| 中国美女看黄片| av天堂在线播放| 国产av一区在线观看免费| 露出奶头的视频| 精品少妇一区二区三区视频日本电影| 日本撒尿小便嘘嘘汇集6| 国产私拍福利视频在线观看| 在线观看舔阴道视频| 黄片大片在线免费观看| 亚洲国产精品sss在线观看| 欧美一级a爱片免费观看看 | 在线观看www视频免费| www国产在线视频色| 国产欧美日韩精品亚洲av| 天天一区二区日本电影三级| 国产真实乱freesex| 日本 av在线| 一级毛片高清免费大全| 又大又爽又粗| 精品久久久久久,| 国内精品久久久久精免费| 国产成人精品久久二区二区91| 午夜福利高清视频| 丰满人妻熟妇乱又伦精品不卡| or卡值多少钱| 动漫黄色视频在线观看| 特大巨黑吊av在线直播| 久久久久久久午夜电影| 国产视频一区二区在线看| 97超级碰碰碰精品色视频在线观看| 亚洲,欧美精品.| 天天一区二区日本电影三级| 欧美日韩福利视频一区二区| 99久久精品热视频| 国产午夜精品论理片| 国产av又大| av欧美777| 亚洲av成人一区二区三| 禁无遮挡网站| 色噜噜av男人的天堂激情| 中文字幕高清在线视频| 在线观看午夜福利视频| 在线视频色国产色| 成人欧美大片| 又大又爽又粗| 一级a爱片免费观看的视频| 一级毛片精品| 日韩大码丰满熟妇| 18禁黄网站禁片免费观看直播| 国产1区2区3区精品| www.www免费av| 97人妻精品一区二区三区麻豆| 亚洲精品av麻豆狂野| 亚洲精品久久国产高清桃花| 欧美av亚洲av综合av国产av| 国产av麻豆久久久久久久| 免费看十八禁软件| 亚洲国产精品999在线| 国产精品野战在线观看| 国产成人影院久久av| av国产免费在线观看| 人人妻人人看人人澡| 两性夫妻黄色片| 欧美黑人精品巨大| 国产午夜精品论理片| 在线观看日韩欧美| 亚洲激情在线av| 757午夜福利合集在线观看| 非洲黑人性xxxx精品又粗又长| 国产精品一及| 国产免费男女视频| 一进一出抽搐动态| 岛国视频午夜一区免费看| 久久久久国产精品人妻aⅴ院| x7x7x7水蜜桃| 天堂√8在线中文|