• <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
    日本成人三级电影网站| 18禁观看日本| 婷婷六月久久综合丁香| 国内毛片毛片毛片毛片毛片| 美女 人体艺术 gogo| 免费在线观看日本一区| 国产一区二区三区视频了| 亚洲欧美日韩无卡精品| 老司机在亚洲福利影院| 亚洲国产欧美网| 国产精品 欧美亚洲| 99精品久久久久人妻精品| 国产精品 国内视频| 90打野战视频偷拍视频| 欧美日韩中文字幕国产精品一区二区三区| 美女黄网站色视频| 久久久成人免费电影| 国产一区在线观看成人免费| 最新中文字幕久久久久 | 亚洲美女视频黄频| 欧美日本亚洲视频在线播放| 亚洲国产精品999在线| 亚洲国产看品久久| 亚洲美女黄片视频| 午夜福利成人在线免费观看| 高清毛片免费观看视频网站| 中国美女看黄片| 欧美日韩一级在线毛片| 久久人妻av系列| 小说图片视频综合网站| 中文字幕人妻丝袜一区二区| a在线观看视频网站| 亚洲国产精品久久男人天堂| xxx96com| 法律面前人人平等表现在哪些方面| 亚洲人成伊人成综合网2020| 欧美黄色淫秽网站| 国产免费男女视频| 亚洲成人精品中文字幕电影| а√天堂www在线а√下载| 两人在一起打扑克的视频| 丁香六月欧美| 国产极品精品免费视频能看的| 国产亚洲精品一区二区www| 法律面前人人平等表现在哪些方面| 欧美黄色淫秽网站| 99久久99久久久精品蜜桃| 18禁观看日本| 国产伦人伦偷精品视频| 法律面前人人平等表现在哪些方面| 色老头精品视频在线观看| 日韩欧美精品v在线| 国产不卡一卡二| 免费看a级黄色片| 日韩中文字幕欧美一区二区| 中文字幕人成人乱码亚洲影| 亚洲欧美一区二区三区黑人| 久久伊人香网站| 网址你懂的国产日韩在线| 亚洲美女视频黄频| 1024手机看黄色片| 日韩欧美国产在线观看| 午夜精品在线福利| 最近最新免费中文字幕在线| 欧美日韩福利视频一区二区| 国产激情偷乱视频一区二区| 亚洲av免费在线观看| 两个人视频免费观看高清| 国产69精品久久久久777片 | 网址你懂的国产日韩在线| 亚洲无线观看免费| 免费观看精品视频网站| 国产高清videossex| 精品国内亚洲2022精品成人| 美女被艹到高潮喷水动态| 757午夜福利合集在线观看| 国产精品98久久久久久宅男小说| 人妻夜夜爽99麻豆av| 2021天堂中文幕一二区在线观| avwww免费| 国产黄a三级三级三级人| 一个人免费在线观看电影 | 在线国产一区二区在线| 国产伦精品一区二区三区四那| 黄色成人免费大全| 99热只有精品国产| 久久久久性生活片| 日本熟妇午夜| 国产主播在线观看一区二区| 久久久国产成人免费| 免费一级毛片在线播放高清视频| 熟妇人妻久久中文字幕3abv| 国产午夜精品论理片| 在线国产一区二区在线| 悠悠久久av| 天堂√8在线中文| 最近在线观看免费完整版| 人妻丰满熟妇av一区二区三区| 午夜视频精品福利| 免费观看人在逋| 熟女少妇亚洲综合色aaa.| 日日夜夜操网爽| 亚洲欧美日韩高清在线视频| 美女扒开内裤让男人捅视频| 欧美高清成人免费视频www| 成人国产一区最新在线观看| 亚洲av第一区精品v没综合| www.精华液| 在线a可以看的网站| 久久中文看片网| 久久久色成人| 久久人妻av系列| 国产精品永久免费网站| 级片在线观看| 亚洲欧美精品综合久久99| 国产精品 欧美亚洲| 级片在线观看| 成年版毛片免费区| 亚洲美女黄片视频| 久久天躁狠狠躁夜夜2o2o| 色尼玛亚洲综合影院| 国产在线精品亚洲第一网站| 91老司机精品| 淫妇啪啪啪对白视频| 亚洲在线自拍视频| 国产激情偷乱视频一区二区| 国产精品98久久久久久宅男小说| 老熟妇仑乱视频hdxx| 国产免费av片在线观看野外av| 99久久精品热视频| 亚洲国产欧洲综合997久久,| 90打野战视频偷拍视频| 久久久久久久精品吃奶| 午夜福利视频1000在线观看| 成年免费大片在线观看| 日韩大尺度精品在线看网址| 好男人电影高清在线观看| www.999成人在线观看| 五月玫瑰六月丁香| 美女高潮喷水抽搐中文字幕| 午夜福利免费观看在线| 偷拍熟女少妇极品色| 给我免费播放毛片高清在线观看| 狂野欧美激情性xxxx| 国产av麻豆久久久久久久| 全区人妻精品视频| 国产成人福利小说| 亚洲aⅴ乱码一区二区在线播放| 天天添夜夜摸| 欧美日韩乱码在线| 夜夜爽天天搞| 国产精品美女特级片免费视频播放器 | 观看免费一级毛片| 亚洲最大成人中文| 日本精品一区二区三区蜜桃| 亚洲国产精品成人综合色| 亚洲精华国产精华精| 老鸭窝网址在线观看| 欧美黑人巨大hd| 日本撒尿小便嘘嘘汇集6| 久久人妻av系列| 成人高潮视频无遮挡免费网站| 亚洲精品一区av在线观看| av视频在线观看入口| 欧美极品一区二区三区四区| 国产人伦9x9x在线观看| 日韩大尺度精品在线看网址| 美女被艹到高潮喷水动态| 91字幕亚洲| 在线观看美女被高潮喷水网站 | 免费人成视频x8x8入口观看| 看免费av毛片| 国产蜜桃级精品一区二区三区| 久久天躁狠狠躁夜夜2o2o| 美女cb高潮喷水在线观看 | 黑人巨大精品欧美一区二区mp4| 国产精品久久久久久久电影 | 欧美日韩黄片免| 一本综合久久免费| 免费无遮挡裸体视频| 动漫黄色视频在线观看| 国产一区二区三区视频了| 亚洲欧美日韩无卡精品| 久久久久久大精品| 午夜精品在线福利| 国产又黄又爽又无遮挡在线| 波多野结衣高清无吗| 91久久精品国产一区二区成人 | 精品国产三级普通话版| 久久精品影院6| 国产精品女同一区二区软件 | 99久久精品一区二区三区| 午夜两性在线视频| 狠狠狠狠99中文字幕| 欧美中文日本在线观看视频| 色综合欧美亚洲国产小说| 熟女少妇亚洲综合色aaa.| 中文字幕久久专区| 国产精品影院久久| 久久精品91无色码中文字幕| 特级一级黄色大片| 国产精品日韩av在线免费观看| 校园春色视频在线观看| 国内精品美女久久久久久| 国产精品女同一区二区软件 | 成人亚洲精品av一区二区| 久久人妻av系列| 伊人久久大香线蕉亚洲五| 久久亚洲真实| 亚洲国产日韩欧美精品在线观看 | 又大又爽又粗| 无限看片的www在线观看| 99热只有精品国产| 午夜激情欧美在线| 婷婷丁香在线五月| 久久国产精品人妻蜜桃| 国产精品99久久久久久久久| 色吧在线观看| 成人鲁丝片一二三区免费| 熟女电影av网| 欧美激情在线99| 欧美最黄视频在线播放免费| 色吧在线观看| 日韩欧美在线二视频| www.自偷自拍.com| 亚洲av电影在线进入| a级毛片在线看网站| 99riav亚洲国产免费| 女警被强在线播放| 在线观看舔阴道视频| 久久这里只有精品中国| 桃色一区二区三区在线观看| 成年女人看的毛片在线观看| 99国产精品一区二区蜜桃av| 成人永久免费在线观看视频| 久久香蕉国产精品| 亚洲精品乱码久久久v下载方式 | 啦啦啦观看免费观看视频高清| 91九色精品人成在线观看| 亚洲中文字幕日韩| 国产极品精品免费视频能看的| 国产三级在线视频| 黄色片一级片一级黄色片| 麻豆成人av在线观看| 中文字幕熟女人妻在线| 国产aⅴ精品一区二区三区波| 超碰成人久久| av天堂中文字幕网| 18禁观看日本| 国产精品香港三级国产av潘金莲| 叶爱在线成人免费视频播放| 欧美中文日本在线观看视频| 日韩精品青青久久久久久| 一本久久中文字幕| 欧美日韩精品网址| 亚洲欧美日韩卡通动漫| 国产伦一二天堂av在线观看| a级毛片a级免费在线| 美女扒开内裤让男人捅视频| 亚洲一区二区三区不卡视频| 一区二区三区高清视频在线| 老熟妇仑乱视频hdxx| 久久香蕉国产精品| 国产精品一区二区免费欧美| 国产高清激情床上av| 中文字幕精品亚洲无线码一区| 亚洲国产高清在线一区二区三| 叶爱在线成人免费视频播放| 欧美3d第一页| 51午夜福利影视在线观看| 99国产综合亚洲精品| 又爽又黄无遮挡网站| 级片在线观看| 亚洲人成网站在线播放欧美日韩| www.999成人在线观看| а√天堂www在线а√下载| 精品无人区乱码1区二区| 亚洲av片天天在线观看| 国产成人精品久久二区二区免费| 最好的美女福利视频网| 最近最新中文字幕大全免费视频| 国产成人av激情在线播放| 两个人的视频大全免费| 国产一区二区激情短视频| 亚洲成a人片在线一区二区| 国产伦精品一区二区三区四那| 日韩欧美精品v在线| 欧美激情在线99| 中文字幕久久专区| 国产三级中文精品| 成人无遮挡网站| 人妻丰满熟妇av一区二区三区| 国产精品久久久久久亚洲av鲁大| 亚洲性夜色夜夜综合| 国产成人欧美在线观看| 免费无遮挡裸体视频| 国产三级黄色录像| 亚洲avbb在线观看| 男人舔奶头视频| 久久久久久久午夜电影| 少妇的丰满在线观看| 亚洲精品在线美女| 午夜亚洲福利在线播放| 男女那种视频在线观看| 嫩草影视91久久| 国产视频一区二区在线看| 一进一出好大好爽视频| 一进一出抽搐动态| 国产精品98久久久久久宅男小说| 给我免费播放毛片高清在线观看| 久久久久国产一级毛片高清牌| 国产免费男女视频| 无遮挡黄片免费观看| 性色av乱码一区二区三区2| 欧美日韩亚洲国产一区二区在线观看| 精品不卡国产一区二区三区| 欧美一区二区国产精品久久精品| 日本撒尿小便嘘嘘汇集6| 日韩精品中文字幕看吧| 1024香蕉在线观看| 成年人黄色毛片网站| 亚洲 欧美一区二区三区| 最好的美女福利视频网| 午夜福利在线观看吧| 黑人欧美特级aaaaaa片| 国产高清视频在线播放一区| 黄色视频,在线免费观看| 桃色一区二区三区在线观看| 久久久久久九九精品二区国产| 欧美成狂野欧美在线观看| 99视频精品全部免费 在线 | 91九色精品人成在线观看| 久久久久久久午夜电影| 色综合婷婷激情| 成年女人永久免费观看视频| 成人精品一区二区免费| 国产一区二区三区在线臀色熟女| 看黄色毛片网站| 亚洲精品美女久久av网站| 男女那种视频在线观看| 十八禁网站免费在线| 亚洲天堂国产精品一区在线| 97超视频在线观看视频| 亚洲在线观看片| 美女黄网站色视频| 精品国产乱子伦一区二区三区| 91麻豆精品激情在线观看国产| 熟女少妇亚洲综合色aaa.| www日本黄色视频网| 久久久国产成人精品二区| av天堂在线播放| 国内精品久久久久久久电影| 久久欧美精品欧美久久欧美| 日本撒尿小便嘘嘘汇集6| 人人妻,人人澡人人爽秒播| 三级毛片av免费| 1000部很黄的大片| 国产麻豆成人av免费视频| 午夜福利高清视频| 国产日本99.免费观看| 男人舔女人的私密视频| 国产视频内射| 久久国产乱子伦精品免费另类| 国产69精品久久久久777片 | 日本一本二区三区精品| 免费高清视频大片| 国产精品电影一区二区三区| 午夜福利免费观看在线| 欧美日韩国产亚洲二区| 岛国在线免费视频观看| 精品福利观看| 一边摸一边抽搐一进一小说| 国产99白浆流出| 男人和女人高潮做爰伦理| 身体一侧抽搐| 美女 人体艺术 gogo| 国产精品久久久久久亚洲av鲁大| 欧美黄色淫秽网站| 色噜噜av男人的天堂激情| 99热这里只有精品一区 | 国产单亲对白刺激| 巨乳人妻的诱惑在线观看| 老司机午夜福利在线观看视频| 精华霜和精华液先用哪个| 精品日产1卡2卡| 丁香欧美五月| 中出人妻视频一区二区| 久久婷婷人人爽人人干人人爱| 国产激情久久老熟女| avwww免费| 国产久久久一区二区三区| 在线永久观看黄色视频| 夜夜躁狠狠躁天天躁| 国产黄a三级三级三级人| 一级毛片高清免费大全| 国产精品自产拍在线观看55亚洲| 天天一区二区日本电影三级| 国产精品久久视频播放| 亚洲一区二区三区色噜噜| 99精品久久久久人妻精品| 亚洲 国产 在线| 久9热在线精品视频| 18禁观看日本| 欧美3d第一页| 草草在线视频免费看| 国产aⅴ精品一区二区三区波| 亚洲,欧美精品.| 怎么达到女性高潮| 亚洲av日韩精品久久久久久密| 91在线观看av| 亚洲aⅴ乱码一区二区在线播放| 黄色成人免费大全| 午夜福利视频1000在线观看| 日韩免费av在线播放| 国内精品久久久久久久电影| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产精品久久男人天堂| 国产精品,欧美在线| 午夜久久久久精精品| 欧美黄色片欧美黄色片| 一级毛片高清免费大全| 18禁美女被吸乳视频| 色老头精品视频在线观看| 欧美xxxx黑人xx丫x性爽| 国产亚洲av高清不卡| tocl精华| 日韩欧美国产一区二区入口| 欧美成人免费av一区二区三区| 最新中文字幕久久久久 | 最新中文字幕久久久久 | 男人舔女人下体高潮全视频| www.熟女人妻精品国产| 亚洲熟妇熟女久久| 亚洲国产日韩欧美精品在线观看 | 美女扒开内裤让男人捅视频| 亚洲av美国av| 99国产精品99久久久久| 91在线精品国自产拍蜜月 | 18禁国产床啪视频网站| 免费一级毛片在线播放高清视频| 国产伦人伦偷精品视频| 女同久久另类99精品国产91| 熟女人妻精品中文字幕| 久久久国产欧美日韩av| av黄色大香蕉| 变态另类丝袜制服| 国产精品99久久久久久久久| 一边摸一边抽搐一进一小说| 日韩免费av在线播放| 首页视频小说图片口味搜索| 丰满的人妻完整版| 国产精品野战在线观看| 男女之事视频高清在线观看| 熟女电影av网| 在线免费观看不下载黄p国产 | 夜夜躁狠狠躁天天躁| 久久久久久大精品| 国产精品一区二区三区四区免费观看 | 精华霜和精华液先用哪个| 亚洲精品久久国产高清桃花| 日韩免费av在线播放| 欧美丝袜亚洲另类 | 精品一区二区三区av网在线观看| 一个人免费在线观看的高清视频| 欧美日韩瑟瑟在线播放| 精品国产三级普通话版| 怎么达到女性高潮| 香蕉丝袜av| 黄色成人免费大全| 操出白浆在线播放| 国产一区在线观看成人免费| 国产在线精品亚洲第一网站| 国产精品自产拍在线观看55亚洲| 国产精品98久久久久久宅男小说| 久久性视频一级片| 久久欧美精品欧美久久欧美| 法律面前人人平等表现在哪些方面| 男女床上黄色一级片免费看| 国产午夜精品论理片| 亚洲精华国产精华精| 1000部很黄的大片| 狂野欧美激情性xxxx| 久久天堂一区二区三区四区| 小蜜桃在线观看免费完整版高清| 亚洲aⅴ乱码一区二区在线播放| 动漫黄色视频在线观看| 色视频www国产| 国产伦精品一区二区三区视频9 | 欧美日韩瑟瑟在线播放| 亚洲成av人片免费观看| 亚洲av成人不卡在线观看播放网| 法律面前人人平等表现在哪些方面| 免费高清视频大片| 不卡av一区二区三区| 桃红色精品国产亚洲av| 男人舔奶头视频| 精品久久久久久成人av| 97超视频在线观看视频| 亚洲无线观看免费| 嫁个100分男人电影在线观看| www.自偷自拍.com| 亚洲精品乱码久久久v下载方式 | 很黄的视频免费| 热99在线观看视频| 精品国内亚洲2022精品成人| 亚洲七黄色美女视频| 国产午夜福利久久久久久| 夜夜夜夜夜久久久久| 草草在线视频免费看| 一级毛片精品| 最近最新中文字幕大全免费视频| av在线蜜桃| 在线看三级毛片| 色综合亚洲欧美另类图片| 成人精品一区二区免费| 日本 欧美在线| 日韩欧美 国产精品| 91麻豆av在线| 国产av不卡久久| 在线十欧美十亚洲十日本专区| 欧美成狂野欧美在线观看| 亚洲自拍偷在线| xxx96com| 一夜夜www| 身体一侧抽搐| 成人三级做爰电影| 免费高清视频大片| 国产精品一区二区三区四区免费观看 | 三级国产精品片| 草草在线视频免费看| 国产在线男女| videossex国产| 汤姆久久久久久久影院中文字幕 | 秋霞伦理黄片| 天天躁日日操中文字幕| 国产精品久久久久久久电影| 免费看av在线观看网站| 国产高清视频在线观看网站| 亚洲最大成人中文| 午夜爱爱视频在线播放| 亚洲美女视频黄频| 久久久久久久久中文| 国模一区二区三区四区视频| 久久99精品国语久久久| videos熟女内射| 嫩草影院精品99| 国产高清不卡午夜福利| 国产中年淑女户外野战色| 国产精品久久久久久精品电影| 91精品一卡2卡3卡4卡| 久久久久九九精品影院| 国产精品福利在线免费观看| 一级爰片在线观看| 午夜亚洲福利在线播放| 纵有疾风起免费观看全集完整版 | 亚洲在线观看片| 亚洲国产欧美人成| 一级黄色大片毛片| 久久精品久久久久久久性| av又黄又爽大尺度在线免费看 | 成人特级av手机在线观看| 在线观看66精品国产| 青春草国产在线视频| 爱豆传媒免费全集在线观看| 亚洲最大成人中文| 精品人妻熟女av久视频| 亚洲美女视频黄频| 久久亚洲国产成人精品v| 禁无遮挡网站| 黄片wwwwww| 欧美变态另类bdsm刘玥| 久久久久九九精品影院| 久久国产乱子免费精品| 日韩视频在线欧美| 男女视频在线观看网站免费| 蜜臀久久99精品久久宅男| 日本一本二区三区精品| 久久6这里有精品| 国内少妇人妻偷人精品xxx网站| 精品久久久久久成人av| 嫩草影院精品99| 成年女人看的毛片在线观看| 老司机福利观看| 亚洲精品乱久久久久久| www.av在线官网国产| 乱系列少妇在线播放| 久久久久久久久大av| 日韩国内少妇激情av| 熟女人妻精品中文字幕| 国产成人免费观看mmmm| 天天一区二区日本电影三级| 国产精品久久久久久精品电影小说 | 99热这里只有是精品在线观看| 国产欧美另类精品又又久久亚洲欧美| 日韩国内少妇激情av| 国产成人aa在线观看| 日韩大片免费观看网站 | 99九九线精品视频在线观看视频| 国内精品美女久久久久久| 午夜福利成人在线免费观看| 国产精品爽爽va在线观看网站| 中文亚洲av片在线观看爽| 久久精品国产99精品国产亚洲性色| 麻豆国产97在线/欧美| 最近中文字幕2019免费版| 国产单亲对白刺激| 久久精品综合一区二区三区| 国产精品久久久久久久电影| 成人亚洲精品av一区二区| 嫩草影院精品99| 人妻夜夜爽99麻豆av| 国产乱人偷精品视频| 久久99热这里只频精品6学生 | 日韩亚洲欧美综合| 免费播放大片免费观看视频在线观看 |