• <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
    能在线免费看毛片的网站| 国产视频首页在线观看| 欧美区成人在线视频| kizo精华| 国产亚洲精品av在线| 亚洲精品一二三| 国产精品.久久久| 观看免费一级毛片| 人人妻人人澡人人爽人人夜夜 | 免费观看a级毛片全部| 国产亚洲91精品色在线| 国内精品宾馆在线| 青春草视频在线免费观看| 女人久久www免费人成看片| 国产精品福利在线免费观看| 又粗又硬又长又爽又黄的视频| 能在线免费看毛片的网站| 女人被狂操c到高潮| 国产在线男女| 日韩人妻高清精品专区| 亚洲国产欧美人成| 久久人人爽人人片av| 联通29元200g的流量卡| 边亲边吃奶的免费视频| 男女边吃奶边做爰视频| 秋霞在线观看毛片| 国产激情偷乱视频一区二区| 日韩av不卡免费在线播放| 中文精品一卡2卡3卡4更新| 一区二区三区免费毛片| 精品久久久久久久久久久久久| 国产成人精品福利久久| 九九久久精品国产亚洲av麻豆| 国产视频首页在线观看| 秋霞在线观看毛片| 久久精品国产亚洲av涩爱| 亚洲av中文av极速乱| 一边亲一边摸免费视频| 一本一本综合久久| 日韩一区二区视频免费看| 国产成人精品久久久久久| 免费电影在线观看免费观看| 男插女下体视频免费在线播放| 色尼玛亚洲综合影院| 看免费成人av毛片| 久久99热这里只有精品18| 亚洲国产精品成人综合色| 99久国产av精品| 亚洲,欧美,日韩| 久久99精品国语久久久| 久久久精品94久久精品| 欧美xxxx性猛交bbbb| 国产淫语在线视频| 22中文网久久字幕| 亚洲av一区综合| 日本黄色片子视频| 国产精品一区二区三区四区免费观看| 日韩制服骚丝袜av| 国产片特级美女逼逼视频| 一级爰片在线观看| 久久热精品热| 日本免费a在线| 久久精品久久久久久久性| 欧美高清成人免费视频www| 22中文网久久字幕| 久久99热这里只频精品6学生| 欧美3d第一页| 三级经典国产精品| 久久久久精品久久久久真实原创| 午夜福利在线观看吧| av线在线观看网站| 国内精品美女久久久久久| 亚洲18禁久久av| 99久国产av精品国产电影| 十八禁国产超污无遮挡网站| 国产高潮美女av| 亚洲内射少妇av| 欧美xxxx性猛交bbbb| 亚洲精品一二三| 亚洲欧美精品专区久久| 国产免费视频播放在线视频 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99久久精品国产国产毛片| 亚洲av在线观看美女高潮| 欧美最新免费一区二区三区| 亚洲精品乱久久久久久| 简卡轻食公司| 亚洲欧美日韩东京热| 男人舔女人下体高潮全视频| 亚洲va在线va天堂va国产| 免费黄频网站在线观看国产| 天天躁夜夜躁狠狠久久av| 丰满乱子伦码专区| 欧美成人精品欧美一级黄| 免费大片黄手机在线观看| 亚洲成人精品中文字幕电影| 大片免费播放器 马上看| 伊人久久国产一区二区| 午夜福利网站1000一区二区三区| 18禁在线播放成人免费| 高清在线视频一区二区三区| 男的添女的下面高潮视频| 亚洲欧洲日产国产| 国产69精品久久久久777片| 欧美+日韩+精品| av国产免费在线观看| 国产免费视频播放在线视频 | a级毛片免费高清观看在线播放| 国产麻豆成人av免费视频| 久久久欧美国产精品| 免费看a级黄色片| 少妇裸体淫交视频免费看高清| 91精品国产九色| 特级一级黄色大片| 国产有黄有色有爽视频| 久久久久久久国产电影| 高清欧美精品videossex| 简卡轻食公司| 精品久久久久久成人av| 国产在线男女| 你懂的网址亚洲精品在线观看| av播播在线观看一区| 亚州av有码| 最近中文字幕高清免费大全6| 久久久久精品性色| 国产一级毛片在线| 欧美97在线视频| 成人亚洲精品一区在线观看 | 久久精品久久精品一区二区三区| 精品久久久久久久久久久久久| 国产成人精品久久久久久| 激情五月婷婷亚洲| 亚洲精品乱码久久久v下载方式| 免费观看av网站的网址| 久久99热6这里只有精品| 国产午夜精品一二区理论片| 黄片wwwwww| 国产精品久久视频播放| 亚洲人成网站在线观看播放| 亚洲精品日韩av片在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲在线自拍视频| 街头女战士在线观看网站| 成人午夜高清在线视频| 亚洲在线自拍视频| 精华霜和精华液先用哪个| 欧美成人a在线观看| 免费少妇av软件| 男女边摸边吃奶| 综合色丁香网| 中文在线观看免费www的网站| 我的老师免费观看完整版| 成人漫画全彩无遮挡| 中文字幕制服av| 插逼视频在线观看| 深爱激情五月婷婷| 又黄又爽又刺激的免费视频.| 久久久亚洲精品成人影院| 99九九线精品视频在线观看视频| 色尼玛亚洲综合影院| 中文天堂在线官网| 欧美zozozo另类| 久久久a久久爽久久v久久| 久久久亚洲精品成人影院| 国产黄片视频在线免费观看| 精品一区二区三卡| 日韩不卡一区二区三区视频在线| 成人av在线播放网站| 国产伦理片在线播放av一区| 高清日韩中文字幕在线| 天堂俺去俺来也www色官网 | 日日摸夜夜添夜夜爱| 中文字幕人妻熟人妻熟丝袜美| 1000部很黄的大片| 日韩av在线大香蕉| kizo精华| 国产精品一区二区三区四区久久| 国产乱来视频区| 人人妻人人澡人人爽人人夜夜 | 国产亚洲av片在线观看秒播厂 | 好男人在线观看高清免费视频| 国产亚洲精品av在线| 精品一区二区免费观看| 国产亚洲一区二区精品| 纵有疾风起免费观看全集完整版 | 日韩一区二区视频免费看| av国产免费在线观看| 毛片女人毛片| 亚洲最大成人中文| 男女边吃奶边做爰视频| av黄色大香蕉| 亚洲欧美日韩无卡精品| 日本一二三区视频观看| 亚洲精华国产精华液的使用体验| 最近最新中文字幕免费大全7| 97精品久久久久久久久久精品| 亚洲最大成人中文| 91久久精品国产一区二区三区| 中文字幕制服av| 永久免费av网站大全| 亚洲综合色惰| 成人午夜高清在线视频| 亚洲怡红院男人天堂| 大陆偷拍与自拍| 午夜福利高清视频| 丰满乱子伦码专区| 国产亚洲一区二区精品| 天天一区二区日本电影三级| 日韩中字成人| 中文精品一卡2卡3卡4更新| 亚洲国产精品成人久久小说| 别揉我奶头 嗯啊视频| 黄色一级大片看看| 亚洲成人精品中文字幕电影| 最近视频中文字幕2019在线8| 在线免费观看不下载黄p国产| 成人一区二区视频在线观看| 精品熟女少妇av免费看| eeuss影院久久| 成人漫画全彩无遮挡| 欧美xxxx性猛交bbbb| 精品国内亚洲2022精品成人| 人人妻人人澡人人爽人人夜夜 | 中文欧美无线码| 日本-黄色视频高清免费观看| 麻豆国产97在线/欧美| 欧美日韩国产mv在线观看视频 | 亚洲av男天堂| 国产黄片视频在线免费观看| 午夜精品国产一区二区电影 | 高清毛片免费看| 成年女人在线观看亚洲视频 | 全区人妻精品视频| 亚洲欧洲国产日韩| 狂野欧美激情性xxxx在线观看| 国产片特级美女逼逼视频| 日韩视频在线欧美| 国产免费一级a男人的天堂| 午夜激情久久久久久久| 成年人午夜在线观看视频 | 91午夜精品亚洲一区二区三区| 97超视频在线观看视频| 插阴视频在线观看视频| 国产大屁股一区二区在线视频| av天堂中文字幕网| 菩萨蛮人人尽说江南好唐韦庄| 欧美成人一区二区免费高清观看| 免费看日本二区| 国产在线男女| 日产精品乱码卡一卡2卡三| 国产视频首页在线观看| 午夜激情福利司机影院| 久久精品国产亚洲av涩爱| 国产精品久久久久久久电影| 日日撸夜夜添| 婷婷色综合www| 蜜桃久久精品国产亚洲av| 色网站视频免费| 国产精品av视频在线免费观看| 欧美zozozo另类| 麻豆av噜噜一区二区三区| 精品久久久久久久人妻蜜臀av| 22中文网久久字幕| 亚洲丝袜综合中文字幕| 十八禁网站网址无遮挡 | 高清av免费在线| 日韩 亚洲 欧美在线| 联通29元200g的流量卡| 三级国产精品欧美在线观看| 国产成人91sexporn| 一区二区三区四区激情视频| 色尼玛亚洲综合影院| 国产免费福利视频在线观看| 国产精品一区二区三区四区久久| 国语对白做爰xxxⅹ性视频网站| 日本三级黄在线观看| 成年女人看的毛片在线观看| 一级毛片 在线播放| 国产麻豆成人av免费视频| 少妇的逼水好多| 久久国产乱子免费精品| 久久99热这里只频精品6学生| 狠狠精品人妻久久久久久综合| 亚洲成色77777| 亚洲一级一片aⅴ在线观看| 美女内射精品一级片tv| 午夜精品国产一区二区电影 | 欧美成人精品欧美一级黄| 91狼人影院| av免费观看日本| 久久精品久久精品一区二区三区| 男女那种视频在线观看| 大话2 男鬼变身卡| 国产一区亚洲一区在线观看| 2021少妇久久久久久久久久久| 波多野结衣巨乳人妻| 91aial.com中文字幕在线观看| 一级二级三级毛片免费看| 久久草成人影院| 在线天堂最新版资源| 亚洲av男天堂| 大片免费播放器 马上看| 丝瓜视频免费看黄片| 国产毛片a区久久久久| 日韩精品有码人妻一区| 在现免费观看毛片| 国产日韩欧美在线精品| 中文欧美无线码| 亚洲精品影视一区二区三区av| 亚洲久久久久久中文字幕| 色综合色国产| 国产精品国产三级国产av玫瑰| 国产乱来视频区| 99re6热这里在线精品视频| 美女国产视频在线观看| 日韩 亚洲 欧美在线| 免费大片黄手机在线观看| 久久久久精品久久久久真实原创| 久久久久九九精品影院| 亚洲天堂国产精品一区在线| 国产高清不卡午夜福利| 97在线视频观看| 午夜福利成人在线免费观看| 国内精品一区二区在线观看| 久久久成人免费电影| 在现免费观看毛片| a级一级毛片免费在线观看| 亚洲乱码一区二区免费版| 少妇猛男粗大的猛烈进出视频 | 免费黄频网站在线观看国产| 国产av在哪里看| 国产综合精华液| 99热这里只有是精品50| 一级毛片黄色毛片免费观看视频| 欧美精品国产亚洲| 日韩亚洲欧美综合| 久久久久免费精品人妻一区二区| www.色视频.com| 中文乱码字字幕精品一区二区三区 | 又黄又爽又刺激的免费视频.| 91午夜精品亚洲一区二区三区| 有码 亚洲区| 女人被狂操c到高潮| 国产成人一区二区在线| 99热网站在线观看| 色视频www国产| 欧美+日韩+精品| 777米奇影视久久| 啦啦啦啦在线视频资源| 人妻制服诱惑在线中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久国产一区二区| 秋霞在线观看毛片| 亚洲,欧美,日韩| 高清av免费在线| 女的被弄到高潮叫床怎么办| 身体一侧抽搐| 国国产精品蜜臀av免费| 乱系列少妇在线播放| 伊人久久精品亚洲午夜| 亚洲国产欧美在线一区| 国产淫片久久久久久久久| 少妇熟女aⅴ在线视频| 亚洲内射少妇av| 色尼玛亚洲综合影院| 只有这里有精品99| 国产精品嫩草影院av在线观看| 日本熟妇午夜| 久久久精品94久久精品| 国产乱人偷精品视频| 国产午夜精品久久久久久一区二区三区| 亚洲va在线va天堂va国产| 亚洲在久久综合| 国产成年人精品一区二区| 最近的中文字幕免费完整| 国产伦一二天堂av在线观看| 国内揄拍国产精品人妻在线| 熟女人妻精品中文字幕| 一个人看视频在线观看www免费| 精品午夜福利在线看| 99久国产av精品| 男女边摸边吃奶| 午夜免费激情av| 91狼人影院| 色吧在线观看| 51国产日韩欧美| www.av在线官网国产| 黑人高潮一二区| 免费在线观看成人毛片| 白带黄色成豆腐渣| or卡值多少钱| 久久久久久久久久久丰满| 天天躁夜夜躁狠狠久久av| 卡戴珊不雅视频在线播放| 国产淫语在线视频| 搡女人真爽免费视频火全软件| 大香蕉97超碰在线| 久久久久久久久久成人| 在线天堂最新版资源| 狂野欧美激情性xxxx在线观看| 国产 一区 欧美 日韩| av国产免费在线观看| 国产黄色视频一区二区在线观看| 国产精品一区二区在线观看99 | 成人综合一区亚洲| 日本猛色少妇xxxxx猛交久久| 一级毛片aaaaaa免费看小| 久久久国产一区二区| a级一级毛片免费在线观看| or卡值多少钱| 免费大片黄手机在线观看| av卡一久久| 熟女电影av网| 国产精品一区www在线观看| 日韩人妻高清精品专区| 麻豆av噜噜一区二区三区| 欧美日韩精品成人综合77777| 午夜福利成人在线免费观看| 国产精品久久久久久精品电影| 真实男女啪啪啪动态图| 国产精品日韩av在线免费观看| 午夜免费男女啪啪视频观看| 视频中文字幕在线观看| 国产精品.久久久| 最后的刺客免费高清国语| 国产免费福利视频在线观看| 亚洲精品第二区| 久久99精品国语久久久| 久久久久久久国产电影| 亚洲精品乱久久久久久| 亚洲人成网站在线观看播放| 免费观看性生交大片5| 亚洲经典国产精华液单| 日韩亚洲欧美综合| 免费看a级黄色片| 特大巨黑吊av在线直播| 成人国产麻豆网| av在线老鸭窝| 夫妻午夜视频| 亚洲欧美日韩东京热| 好男人视频免费观看在线| 免费观看av网站的网址| 非洲黑人性xxxx精品又粗又长| 国产 一区精品| 色尼玛亚洲综合影院| 亚洲18禁久久av| 午夜免费男女啪啪视频观看| 成人午夜高清在线视频| 成人欧美大片| 狂野欧美激情性xxxx在线观看| 亚洲,欧美,日韩| av卡一久久| 亚洲欧美一区二区三区黑人 | 日韩制服骚丝袜av| 联通29元200g的流量卡| 青青草视频在线视频观看| 亚洲色图av天堂| 久久久久九九精品影院| 亚洲精品自拍成人| 欧美成人精品欧美一级黄| 99热网站在线观看| 80岁老熟妇乱子伦牲交| 日本av手机在线免费观看| av福利片在线观看| 久99久视频精品免费| 欧美日韩在线观看h| 久久久久久久久久久丰满| 欧美精品国产亚洲| 国产精品99久久久久久久久| 亚洲国产最新在线播放| 91午夜精品亚洲一区二区三区| 久久草成人影院| 亚洲av男天堂| 免费看不卡的av| 啦啦啦中文免费视频观看日本| 欧美区成人在线视频| 国产美女午夜福利| 久久精品国产自在天天线| 久久久久久久亚洲中文字幕| 18禁在线播放成人免费| 最近中文字幕2019免费版| 乱系列少妇在线播放| 免费观看a级毛片全部| 国产亚洲av片在线观看秒播厂 | 国产成人a区在线观看| 国产精品日韩av在线免费观看| 伊人久久精品亚洲午夜| 大香蕉97超碰在线| 亚洲欧美精品专区久久| 在线观看人妻少妇| 听说在线观看完整版免费高清| 网址你懂的国产日韩在线| 国产亚洲av嫩草精品影院| 国产极品天堂在线| 高清毛片免费看| 亚洲av福利一区| 亚洲精品一二三| 六月丁香七月| 男女啪啪激烈高潮av片| 内地一区二区视频在线| 色哟哟·www| 久久精品国产亚洲av涩爱| 亚洲va在线va天堂va国产| 婷婷色麻豆天堂久久| 亚洲精品国产av蜜桃| 国产亚洲91精品色在线| 免费看光身美女| 国产在视频线精品| 2021天堂中文幕一二区在线观| 久久久久久久久久久丰满| 久久久久久国产a免费观看| 欧美激情久久久久久爽电影| 男人爽女人下面视频在线观看| 国产三级在线视频| 在线观看人妻少妇| 欧美一区二区亚洲| 全区人妻精品视频| 亚洲一区高清亚洲精品| 日韩三级伦理在线观看| 99re6热这里在线精品视频| 美女cb高潮喷水在线观看| 人人妻人人澡欧美一区二区| 国产乱人偷精品视频| 综合色av麻豆| 欧美日韩综合久久久久久| 国产成人aa在线观看| 亚洲成人中文字幕在线播放| 亚洲av成人精品一区久久| 国产综合精华液| 国产午夜福利久久久久久| 亚洲一区高清亚洲精品| 国产成年人精品一区二区| 99re6热这里在线精品视频| 亚洲最大成人手机在线| 国产精品福利在线免费观看| 男女那种视频在线观看| 日韩成人av中文字幕在线观看| 高清午夜精品一区二区三区| 99热这里只有是精品在线观看| 国产精品一及| 久久久久久久久久成人| 在线 av 中文字幕| 在线观看一区二区三区| 午夜视频国产福利| 97热精品久久久久久| 亚洲精品乱码久久久久久按摩| av线在线观看网站| 亚洲最大成人av| 国产男女超爽视频在线观看| 三级国产精品片| 在线免费观看的www视频| 久久久久久久久久成人| 好男人在线观看高清免费视频| 欧美日韩精品成人综合77777| 久久精品熟女亚洲av麻豆精品 | 在线播放无遮挡| av国产久精品久网站免费入址| 亚洲精品自拍成人| videos熟女内射| 国产精品精品国产色婷婷| 欧美区成人在线视频| 啦啦啦啦在线视频资源| 亚洲熟妇中文字幕五十中出| 黄色一级大片看看| 免费黄色在线免费观看| 精品一区二区免费观看| 国产成人一区二区在线| 中文乱码字字幕精品一区二区三区 | 在线免费十八禁| 成人国产麻豆网| 亚洲精华国产精华液的使用体验| 亚洲欧美清纯卡通| 婷婷色综合大香蕉| 日韩一区二区视频免费看| 18禁在线无遮挡免费观看视频| av在线观看视频网站免费| 久久久久久久亚洲中文字幕| 久久97久久精品| 欧美精品一区二区大全| 亚洲国产精品国产精品| 久久综合国产亚洲精品| 成人毛片60女人毛片免费| 欧美成人午夜免费资源| 一级二级三级毛片免费看| 大话2 男鬼变身卡| 亚洲美女搞黄在线观看| 男女国产视频网站| 免费观看av网站的网址| 中文字幕免费在线视频6| 激情五月婷婷亚洲| 欧美三级亚洲精品| 极品教师在线视频| 国产精品三级大全| 国产av码专区亚洲av| 一级av片app| 2018国产大陆天天弄谢| 久久久久久久久久人人人人人人| 成人鲁丝片一二三区免费| 亚洲av成人av| 国产精品1区2区在线观看.| 麻豆av噜噜一区二区三区| 国产av不卡久久| 国产精品.久久久| 99热这里只有是精品在线观看| 亚洲av成人av| 老司机影院毛片| 色视频www国产| 国产黄a三级三级三级人| 成年版毛片免费区| 久久精品国产亚洲网站| 亚洲av成人av| 国产高清不卡午夜福利| 嘟嘟电影网在线观看| 国产亚洲午夜精品一区二区久久 | 一级毛片 在线播放|