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

    基于航道與泊位資源協(xié)調(diào)調(diào)度的船舶交通優(yōu)化

    2015-11-29 02:58:35張新宇郭子堅(jiān)王金濤大連海事大學(xué)航海動(dòng)態(tài)仿真和控制交通行業(yè)重點(diǎn)實(shí)驗(yàn)室遼寧大連606大連理工大學(xué)水利工程學(xué)院遼寧大連603
    中國航海 2015年3期
    關(guān)鍵詞:港口船舶模型

    張新宇, 郭子堅(jiān), 王金濤, 林 俊(. 大連海事大學(xué) 航海動(dòng)態(tài)仿真和控制交通行業(yè)重點(diǎn)實(shí)驗(yàn)室, 遼寧 大連 606; . 大連理工大學(xué) 水利工程學(xué)院, 遼寧 大連 603)

    基于航道與泊位資源協(xié)調(diào)調(diào)度的船舶交通優(yōu)化

    張新宇1,2, 郭子堅(jiān)2, 王金濤1, 林 俊1
    (1. 大連海事大學(xué) 航海動(dòng)態(tài)仿真和控制交通行業(yè)重點(diǎn)實(shí)驗(yàn)室, 遼寧 大連 116026; 2. 大連理工大學(xué) 水利工程學(xué)院, 遼寧 大連 116023)

    將港口航道與泊位資源協(xié)調(diào)利用,可最大限度地提高港口運(yùn)營效率。對此,以所有船舶總在港時(shí)間和等待時(shí)間最少為目標(biāo),建立基于單向航道的多目標(biāo)船舶調(diào)度優(yōu)化模型。根據(jù)港口不同區(qū)域的交通流特征,建立初始化約束、流量轉(zhuǎn)換約束、時(shí)隙分配約束和泊位沖突消解約束等多個(gè)約束模型。設(shè)計(jì)多目標(biāo)遺傳算法進(jìn)行求解,并設(shè)計(jì)有針對性的模擬場景進(jìn)行驗(yàn)證;以港口某一繁忙時(shí)段的20艘船舶為例,進(jìn)行調(diào)度試驗(yàn)。最終得到8個(gè)Pareto最優(yōu)解和2個(gè)目標(biāo)的最優(yōu)解:總調(diào)度時(shí)間為2.618 6 h,添加懲罰后總等待時(shí)間為23.012 4 h。結(jié)果表明:該模型及算法給出的調(diào)度方案能有效提高船舶的調(diào)度效率。

    水路運(yùn)輸; 單向航道; 船舶調(diào)度優(yōu)化; 多目標(biāo)遺傳算法; 時(shí)隙分配; 流量轉(zhuǎn)換

    海上運(yùn)輸是目前國際貨物運(yùn)輸中最主要的方式,其服務(wù)質(zhì)量很大程度上受港口對船舶進(jìn)出港服務(wù)效率的影響。目前關(guān)于船舶進(jìn)出港調(diào)度的研究主要側(cè)重于航道流量和泊位分配中的單一方面,而對于兩者協(xié)調(diào)調(diào)度的研究相對較少。

    在航道通過能力方面:邵俊崗等[1]運(yùn)用排隊(duì)論對港口通過能力進(jìn)行了綜合評估;李子強(qiáng)[2]通過分析交通流,對航道通過能力進(jìn)行了研究;代君等[3]運(yùn)用船舶領(lǐng)域模型對港口受限航道通過能力的計(jì)算方法進(jìn)行了探討。這些研究針對的是靜態(tài)模型,不能反映動(dòng)態(tài)交通情況。此外,張瑋等[4-6]分別針對蘇北大運(yùn)河、長江三峽和葛洲壩的單向航段進(jìn)行了規(guī)劃研究,研究重點(diǎn)是通航閘口管制,無須考慮航道通過時(shí)間,并不適用于單向航道的船舶調(diào)度問題。

    在航道使用效率方面:徐國裕等[7]根據(jù)船載自動(dòng)識別系統(tǒng)(Automatic Identification System,AIS)提供的與船舶航行安全相關(guān)的資料,以泊位遠(yuǎn)近、船舶大小、船舶類型和船舶吃水等因素為權(quán)重,提出了單向水道船舶進(jìn)出港最佳排序模式;LLOYD等[8-9]對保持安全水域空間、構(gòu)建整合系統(tǒng)和減少排序沖突給予了充分關(guān)注;OZGECAN等[10]采用多種排隊(duì)模型分析了伊斯坦布爾海峽的狹長水道內(nèi)船舶航行的風(fēng)險(xiǎn)和效率,并針對水道中突發(fā)的交通阻塞問題進(jìn)行了多種中斷排隊(duì)優(yōu)化。PULUGURTHAL等[11]從航道流量角度進(jìn)行了船舶調(diào)度優(yōu)化,側(cè)重于航道的運(yùn)輸效率,忽略了船舶進(jìn)入泊位的作業(yè)效率。

    在模型求解方面:李平等[12]將遺傳算法和混合優(yōu)化策略應(yīng)用到了求解集裝箱港口的泊位調(diào)度問題中;KAP等[13]以船舶因延誤離港而導(dǎo)致的罰金與因不合理停泊導(dǎo)致的額外搬運(yùn)集裝箱費(fèi)用之和最小為目標(biāo),提出了一種基于模擬退火算法的泊位分配模型;CHUANG等[14]設(shè)計(jì)了模糊遺傳算法求解航線調(diào)度優(yōu)化模型,可為船舶航線調(diào)度問題提供輔助支持。

    綜上所述,目前關(guān)于船舶進(jìn)出港調(diào)度的研究主要側(cè)重于航道容量和泊位分配中的單一方面,在實(shí)際的港口調(diào)度中,船舶進(jìn)出港管理效率較低。對此,以航道為切入點(diǎn),協(xié)調(diào)泊位資源,構(gòu)建多目標(biāo)港口船舶調(diào)度優(yōu)化模型,設(shè)計(jì)適于求解的多目標(biāo)遺傳算法,并設(shè)計(jì)模擬數(shù)據(jù)進(jìn)行驗(yàn)證。

    1 問題描述

    船舶調(diào)度的目的是在遵循航道管理規(guī)則和保證港口水域通航安全的條件下,使船舶盡快、盡量多地通過航道,減少船舶延誤。對于單向航道,在港口繁忙的情況下,港口船舶調(diào)度的核心問題是泊位和航道資源的優(yōu)化配置。交通量的大小受船舶進(jìn)出港初始狀態(tài)、同向船舶的安排協(xié)調(diào)、進(jìn)出流量的安排控制以及泊位沖突消解等因素的制約。因此,港口船舶調(diào)度問題是典型的NP-Hard (Non-deterministic Polynomial-time Hard)問題。

    2 模型建立

    2.1概念解釋

    1) 單向航道:同一時(shí)間只允許同一個(gè)方向通航的航道。

    2) 時(shí)隙分配:這里的一個(gè)時(shí)隙指從一艘船舶開始進(jìn)出港到下一艘船舶開始進(jìn)出港的時(shí)間間隔。

    3) 流量轉(zhuǎn)換:船舶進(jìn)港流量和出港流量的協(xié)調(diào)。

    2.2模型假設(shè)

    1) 港口已做好船舶靠泊作業(yè)計(jì)劃,船舶交通管理中心負(fù)責(zé)港口船舶進(jìn)出港調(diào)度。

    2) 天氣對港口船舶調(diào)度沒有影響,拖船和引航員充足且能準(zhǔn)時(shí)到位。

    2.3船舶調(diào)度模型建立

    2.3.1目標(biāo)函數(shù)的建立

    港口方和船方在船舶調(diào)度過程中的利益訴求不同。單純從時(shí)間的維度考慮,港口方希望所有船舶完成調(diào)度的時(shí)間最短,而船方則希望各自船舶在港的等待時(shí)間最短。這里根據(jù)兩利益方關(guān)心的問題,構(gòu)建多目標(biāo)目標(biāo)函數(shù)Min[F,W]

    (1)

    Min(bi-t0i)

    (2)

    式(1)~式(2)中:F為所有船舶完成進(jìn)出港的最少總時(shí)間;W為修正后所有船舶等候進(jìn)出港的最少時(shí)間總和;bi為第i艘船舶進(jìn)/出港開始時(shí)刻;fi為第i艘船舶進(jìn)/出港結(jié)束時(shí)刻;t0i為船舶i進(jìn)入報(bào)告線的時(shí)刻;s為報(bào)告線與航道入口間的距離;si為理論上第i艘船舶的泊位與航道入口間的距離;vi為理論上第i艘船舶從進(jìn)入航道到抵達(dá)泊位處的平均速度;b0為上一時(shí)間段內(nèi)最后一艘船舶進(jìn)港(出港)開始時(shí)刻;f0為上一時(shí)間段內(nèi)最后一艘船舶進(jìn)港(出港)結(jié)束時(shí)刻;h0為本階段第一艘船舶與上一階段最后一艘船舶的安全時(shí)間間隔。

    式(1)表示本階段所有船舶完成進(jìn)出港所使用的總時(shí)間最少;式(2)表示所有船舶等候進(jìn)出港時(shí)間總和最少。這里的船舶等候進(jìn)出港時(shí)間指船舶實(shí)際使用的進(jìn)出港時(shí)間與理論上船舶最快進(jìn)出港時(shí)間的差值,即船舶在港浪費(fèi)的時(shí)間成本。

    2.3.2船舶調(diào)度模型約束的建立

    IOi=(0,1)i=1,2,3,…,n

    (3)

    Rij=(0,1)i=1,2,3,…,n

    (4)

    Bij=(0,1)i=1,2,3,…,n

    (5)

    Berthim=(0,1)i=1,2,3,…,n

    (6)

    (7)

    (bi-b0-h-Bijh1)×[1-(IOi-IO0)2]≥0

    i=1,2,3,…,n

    (8)

    (fi-f0-h-Bijh1)×[1-(IOi-IO0)2]≥0

    i=1,2,3,…,n

    (9)

    (10)

    fi-bi-pi≥0i=1,2,3,…,n

    (11)

    (12)

    (13)

    (IOi-IOj)2≥0i

    (14)

    (bi-bj+MRij-h-Bijh1)×[1-

    (IOi-IOj)2]≥0i

    (15)

    [-bi+bj+M(1-Rij)-h-Bijh1]×[1-

    (IOi-IOj)2]≥0i

    (16)

    (fi-fj+MRij-h-Bijh1)×

    [1-(IOi-IOj)2]≥0i

    (17)

    [-fi+fj+M(1-Rij)-h-Bijh]×[1-

    (IOi-IOj)2]≥0i

    (18)

    (1-Berthin)+M(1-IOi)>0t=bi

    (19)

    式(7)為船舶i開始調(diào)度時(shí)刻必須晚于船舶i預(yù)計(jì)到達(dá)航道入口的時(shí)間; 式(8)~式(10)為保證模型的連續(xù)性構(gòu)建的初始化約束; 式(11)~式(14)為流量轉(zhuǎn)換約束,主要實(shí)現(xiàn)進(jìn)出港船舶的流量在時(shí)間上的轉(zhuǎn)換; 式(15)~式(18)為時(shí)隙分配約束,用于協(xié)調(diào)先后調(diào)度的2艘船舶同向航行時(shí)的安全問題,主要考慮安全時(shí)間間隔、船舶的先后次序和泊位遠(yuǎn)近;式(19)為泊位沖突消解約束,保證任意一艘進(jìn)港船舶調(diào)度開始時(shí)刻,其將要停靠的泊位是空閑的。

    3 算法設(shè)計(jì)

    遺傳算法具有較好的并行計(jì)算能力,是在多目標(biāo)優(yōu)化領(lǐng)域應(yīng)用最成功的算法之一。這里針對港口船舶調(diào)度設(shè)計(jì)遺傳算法進(jìn)行求解。為模型設(shè)計(jì)的多目標(biāo)遺傳算法流程圖見圖1。

    4 實(shí)例驗(yàn)證

    使用MATLAB 2010b設(shè)計(jì)港口船舶調(diào)度多目標(biāo)遺傳算法程序,基于MATLAB GUI實(shí)現(xiàn)計(jì)算結(jié)果的可視化。對某港口某天某時(shí)段的船舶進(jìn)行調(diào)度試驗(yàn),模擬20艘船舶(14進(jìn)6出)準(zhǔn)備時(shí)刻的間隔服從λ=20的負(fù)指數(shù)分布規(guī)律的場景。準(zhǔn)備時(shí)刻是指只要泊位和航道條件允許,即可調(diào)度船舶的時(shí)刻;航行距離是指航道起點(diǎn)與泊位間的距離。模擬數(shù)據(jù)的船舶信息見表1。

    圖1 多目標(biāo)遺傳算法流程圖

    假設(shè)上一階段最后一艘船舶的進(jìn)出港開始時(shí)刻為07:00。表2給出了上一階段最后一艘船舶的相關(guān)信息。

    表2 上一階段最后一艘船舶信息

    泊位是碼頭集團(tuán)按照申請預(yù)先分配完成的,且泊位的條件都滿足要??看暗囊?。表3給出了本階段調(diào)度開始時(shí)刻泊位的空閑和占用狀態(tài),模擬數(shù)據(jù)假設(shè)共有18個(gè)泊位。

    表3 泊位開始時(shí)刻的空閑和占用狀態(tài)

    將上述數(shù)據(jù)讀入程序進(jìn)行運(yùn)算,參數(shù)為:種群大小200,代溝0.9,交叉概率0.95,變異概率0.05,最大遺傳代數(shù)200。運(yùn)行結(jié)果見圖2??煽闯龉灿?個(gè)Pareto最優(yōu)解,20艘船舶調(diào)度完成所用總時(shí)間為2.618 6 h,添加懲罰后的總等待時(shí)間為23.012 4 h;種群在不斷尋優(yōu),并在200代以內(nèi)收斂。

    8個(gè)Pareto最優(yōu)解見表4。在總調(diào)度時(shí)間的方向上,最優(yōu)解為2.618 6 h,此時(shí)總等待時(shí)間的適應(yīng)度值為24.776 7 h。在總等待時(shí)間方向上,最優(yōu)解為23.012 4 h,此時(shí)總調(diào)度時(shí)間為2.698 5 h。這符合多目標(biāo)遺傳算法結(jié)果的特點(diǎn):當(dāng)一個(gè)目標(biāo)改進(jìn)時(shí),另一個(gè)目標(biāo)會(huì)相應(yīng)地衰退。由于船舶準(zhǔn)備調(diào)度時(shí)刻比較集中,因此雖然能增加船舶可調(diào)度時(shí)間段,但同時(shí)也會(huì)造成船舶等待時(shí)間延長。

    對于單向航道,航道上的船舶不能對遇,在進(jìn)行流量轉(zhuǎn)換前必須保證航道上沒有船舶,進(jìn)而會(huì)消耗較多的時(shí)間。因此,理論上對于較長的單向航道而言,流量轉(zhuǎn)換的次數(shù)越少,所有船舶完成調(diào)度所用時(shí)間就越短。8個(gè)Pareto最優(yōu)解的調(diào)度方案都只進(jìn)行3次流量轉(zhuǎn)換,大大提高了船舶的調(diào)度效率。

    每個(gè)Pareto最優(yōu)解代表一個(gè)調(diào)度方案,經(jīng)檢驗(yàn),8個(gè)Pareto最優(yōu)解均符合以上約束,可保證船舶的安全。在實(shí)際應(yīng)用中,若港口追求總調(diào)度時(shí)間最短,可選用第1目標(biāo)值較小的調(diào)度方案Pareto1;若兼顧每艘船舶的等待時(shí)間,可選用第2目標(biāo)值較小的調(diào)度方案Pareto8。

    表4 Pareto 最優(yōu)解

    5 結(jié)束語

    研究了基于單向航道和泊位資源協(xié)調(diào)調(diào)度的船舶交通管理問題,建立了以所有船舶總在港時(shí)間和等待時(shí)間最少為目標(biāo)的基于單向航道的多目標(biāo)船舶調(diào)度優(yōu)化模型。使用MATLAB 2010b編寫了港口船舶調(diào)度多程序,基于MATLAB GUI實(shí)現(xiàn)了計(jì)算結(jié)果的可視化。經(jīng)驗(yàn)證,該模型能大大提高船舶進(jìn)出港作業(yè)的效率。

    [1] 邵俊崗, 許小兵,王煜,等. 洋山港區(qū)運(yùn)營階段的港口通過能力[J]. 上海海事大學(xué)學(xué)報(bào), 2008, 29(4): 25-28.

    [2] 李子強(qiáng). 基于交通流分析的航道通過能力研究[J]. 山東交通學(xué)院學(xué)報(bào),2009,17(4):32-35.

    [3] 代君, 王當(dāng)利, 劉克中. 基于船舶領(lǐng)域模型的港口受限航道通過能力計(jì)算方法[J]. 武漢理工大學(xué)學(xué)報(bào):交通科學(xué)與工程版, 2009, 33(4): 679-682.

    [4] 張瑋, 廖鵬,吳玲莉,等. 船閘通過能力主要影響因素[J].交通運(yùn)輸工程學(xué)報(bào), 2004, 4(3): 108-110.

    [5] 胡洋. 三峽工程五級船閘通過能力初步分析和提高樞紐綜合通過能力的對策[J]. 水運(yùn)工程, 2004(10): 74-77.

    [6] 王小平, 齊歡, 肖恒輝, 等. 基于串聯(lián)排隊(duì)網(wǎng)絡(luò)的三峽葛洲壩水利樞紐聯(lián)合調(diào)度模型[J].交通運(yùn)輸工程學(xué)報(bào), 2006, 6(3): 82-86.

    [7] 徐國裕, 郭涂城, 吳兆麟. 單向水道船舶進(jìn)出港最佳排序模式[J]. 大連海事大學(xué)學(xué)報(bào), 2008, 34(4):150-153.

    [8] LLOYD D N.Integrating VTS with Port Operation[Z]. Hongkong: Proceedings of VTS 2004,2004: 44-52.

    [9] MEINE J,KRAMER M.Creation of a Regionee VTS with Data Fusion and Integration of Vessel Data and Movement Details Over the Area of Several Independent Ports, Coastal And Inland Waterways[Z]. Hongkong: Proceedings of VTS 2004,2004: 235-248.

    [10] OZGECAN S U T. Performance Modeling and Risk Analysis of Transit Vessel Traffic in the Istanbul Strait: Studies on Queues with Multiple Types of Interruptions [D]. New Brunswick: Rutgers the State University of New Jersey, 2011.

    [11] PULUGURTHAL S S, NAMBISAN S S. Using Genetic Algorithms to Evaluate Aircraft Ground Holding Policy Under Static Conditions[J]. Journal of Transportation Engineering, 2011,127(5): 433-441.

    [12] 李平, 孫俊清, 韓梅 .泊位調(diào)度問題的GATS混合優(yōu)化策略[J] .天津理工大學(xué)學(xué)報(bào), 2009, 22(4): 58-61.

    [13] KAP H K, KYUNG C. Moon Beth Scheduling by Simulated Annealing[J]. Transportation Research Part B, 2012(37): 541-560.

    [14] CHUANG T N, LIN C T, KUNG J Y,etal. Planning the Route of Container Vessels: A Fuzzy Genetic Approach[J]. Expert Systems with Applications, 2010(37): 2948-2956.

    OptimizationofVesselSchedulingBasedonCoordinationofOne-WayChannelandBerths

    ZHANGXinyu1,2,GUOZijian2,WANGJintao1,LINJun1
    (1. Key Laboratory of Maritime Dynamic Simulation and Control of Ministry of Transportation, Dalian Maritime University, Dalian 116026, China; 2. School of Hydraulic Engineering, Dalian University of Technology, Dalian 116023, China)

    Careful coordinating the channel and berths can maximize the efficiency of a port. With the objective that total sum of time of all vessels spending in the harbor and waiting for available channel is minimized, a multi-objective optimization model of vessel scheduling based on one-way channel is built. Initialization constraints, traffic conversion constraints, slot allocation constraints and berths conflict resolution models are defined respectively based on traffic flow characteristics in different regions of the port. A multi-objective genetic algorithm is designed to solve the proposed model, and a practical scenario of 20 vessels at a peak period is designed to validate the reasonability and effectiveness of the model. The simulation results show that there are 8 Pareto optimal scheduling scheme options as the optimal solution, in which the total of time is 2.6186 hours and the waiting time is 23.0124 hours. The results of experiment suggest that the proposed scheduling model is of great value to improve the efficiency of vessel traffic.

    waterway transportation; one-way channel; vessel scheduling optimization; multi-objective genetic algorithm; slot allocation; traffic flow transition

    2015-05-11

    國家自然科學(xué)基金(51309043);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(3132014202); 交通部應(yīng)用基礎(chǔ)研究項(xiàng)目(2014329225020);中國博士后科學(xué)基金(2014M551095);遼寧省自然科學(xué)基金(2014025005);遼寧省高校杰出青年學(xué)者成長計(jì)劃(LJQ2014052)

    張新宇(1978—),男,吉林長春人,副教授,博士生,主要研究方向?yàn)榻煌ㄐ畔⒐こ碳翱刂啤-mail:41544391@qq.com

    1000-4653(2015)03-0033-04

    U691+.3

    A

    猜你喜歡
    港口船舶模型
    一半模型
    聚焦港口國際化
    金橋(2022年10期)2022-10-11 03:29:46
    中國港口,屹立東方
    金橋(2022年10期)2022-10-11 03:29:22
    計(jì)算流體力學(xué)在船舶操縱運(yùn)動(dòng)仿真中的應(yīng)用
    《船舶》2022 年度征訂啟事
    船舶(2021年4期)2021-09-07 17:32:22
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    船舶!請加速
    港口上的笑臉
    惠東港口
    海洋世界(2016年12期)2017-01-03 11:33:00
    久久狼人影院| www.熟女人妻精品国产| 久久久久久久国产电影| 国产一区二区激情短视频| 热re99久久国产66热| 成年女人毛片免费观看观看9 | 亚洲精品久久午夜乱码| 一a级毛片在线观看| 777久久人妻少妇嫩草av网站| 免费黄频网站在线观看国产| 午夜日韩欧美国产| 狠狠狠狠99中文字幕| 黑人巨大精品欧美一区二区mp4| 欧美精品人与动牲交sv欧美| 中文字幕色久视频| 欧美日韩亚洲综合一区二区三区_| 啦啦啦免费观看视频1| 又紧又爽又黄一区二区| 久久久久久免费高清国产稀缺| 亚洲成人免费av在线播放| 在线观看一区二区三区激情| 一进一出好大好爽视频| 露出奶头的视频| 90打野战视频偷拍视频| 日韩视频一区二区在线观看| av天堂在线播放| 国产男靠女视频免费网站| 韩国av一区二区三区四区| 91在线观看av| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久人人做人人爽| 午夜视频精品福利| 免费日韩欧美在线观看| 一级毛片女人18水好多| 色尼玛亚洲综合影院| 日本vs欧美在线观看视频| 成人手机av| 国产精品久久久久久人妻精品电影| 十八禁人妻一区二区| 免费观看a级毛片全部| 免费在线观看亚洲国产| 黑人猛操日本美女一级片| 真人做人爱边吃奶动态| 又黄又爽又免费观看的视频| 999久久久精品免费观看国产| 中文欧美无线码| 黄片小视频在线播放| 国产精品国产av在线观看| 亚洲av成人av| 久久久国产精品麻豆| 欧美精品人与动牲交sv欧美| 桃红色精品国产亚洲av| 久久久国产成人免费| 日本一区二区免费在线视频| 极品人妻少妇av视频| 亚洲免费av在线视频| 欧美日韩一级在线毛片| 99香蕉大伊视频| 免费在线观看亚洲国产| 美女扒开内裤让男人捅视频| 人妻丰满熟妇av一区二区三区 | 国产精品永久免费网站| 色婷婷av一区二区三区视频| 欧美色视频一区免费| av网站在线播放免费| 成人三级做爰电影| 高清在线国产一区| 中文字幕最新亚洲高清| 我的亚洲天堂| 国产在视频线精品| 中文亚洲av片在线观看爽 | 久久香蕉激情| 婷婷丁香在线五月| 黄片播放在线免费| 热99国产精品久久久久久7| 悠悠久久av| 亚洲欧美色中文字幕在线| 亚洲av美国av| 中文字幕人妻丝袜制服| 国产精华一区二区三区| 久久久久视频综合| 精品久久久久久久久久免费视频 | 久久热在线av| 超碰成人久久| 欧美日韩亚洲高清精品| 日韩 欧美 亚洲 中文字幕| 亚洲伊人色综图| 变态另类成人亚洲欧美熟女 | 午夜精品久久久久久毛片777| 国产欧美日韩一区二区精品| 青草久久国产| 免费看a级黄色片| 亚洲欧洲精品一区二区精品久久久| 日日夜夜操网爽| 在线观看66精品国产| 国产精品免费视频内射| 在线观看www视频免费| av不卡在线播放| 成人永久免费在线观看视频| 香蕉久久夜色| 天堂中文最新版在线下载| 91精品三级在线观看| 成年动漫av网址| 免费av中文字幕在线| 国产亚洲欧美98| 久久久久国内视频| 宅男免费午夜| 午夜福利免费观看在线| 又大又爽又粗| 波多野结衣一区麻豆| 国产成人av教育| 熟女少妇亚洲综合色aaa.| 人妻 亚洲 视频| 老熟女久久久| 精品一品国产午夜福利视频| 十八禁高潮呻吟视频| 男男h啪啪无遮挡| 99国产精品一区二区三区| 天堂俺去俺来也www色官网| 久久国产精品男人的天堂亚洲| 日韩三级视频一区二区三区| 交换朋友夫妻互换小说| 大香蕉久久成人网| 国产亚洲欧美98| 亚洲精品国产一区二区精华液| 国产人伦9x9x在线观看| 久久精品亚洲熟妇少妇任你| ponron亚洲| 欧美黄色片欧美黄色片| 亚洲av欧美aⅴ国产| 亚洲成国产人片在线观看| 真人做人爱边吃奶动态| 欧美日韩黄片免| 精品国产一区二区三区四区第35| 极品教师在线免费播放| 久久久国产成人免费| 国产人伦9x9x在线观看| 色老头精品视频在线观看| 日韩成人在线观看一区二区三区| 国产精品.久久久| 国产精品成人在线| www.精华液| 久久精品亚洲熟妇少妇任你| 精品人妻1区二区| 成年人免费黄色播放视频| 最近最新免费中文字幕在线| 亚洲中文字幕日韩| 日韩欧美免费精品| 久久久久久久久免费视频了| 村上凉子中文字幕在线| 成人亚洲精品一区在线观看| 国产一区二区激情短视频| 免费一级毛片在线播放高清视频 | 国产一区二区三区视频了| 精品国产国语对白av| 老熟妇仑乱视频hdxx| 天堂动漫精品| 国产色视频综合| 精品午夜福利视频在线观看一区| av天堂久久9| 91老司机精品| 精品卡一卡二卡四卡免费| 欧美最黄视频在线播放免费 | 日本wwww免费看| 两人在一起打扑克的视频| 国产成人av教育| 在线观看免费高清a一片| 999精品在线视频| 国产深夜福利视频在线观看| 在线免费观看的www视频| 91精品三级在线观看| av线在线观看网站| 亚洲熟女毛片儿| 美国免费a级毛片| 下体分泌物呈黄色| 国产精品久久久久成人av| 麻豆乱淫一区二区| 国产精品国产高清国产av | 女警被强在线播放| а√天堂www在线а√下载 | 天天躁日日躁夜夜躁夜夜| 精品久久蜜臀av无| 女人高潮潮喷娇喘18禁视频| 精品少妇一区二区三区视频日本电影| 免费黄频网站在线观看国产| 成年人免费黄色播放视频| 亚洲人成电影免费在线| 国产成人影院久久av| 亚洲欧美一区二区三区黑人| 18禁裸乳无遮挡免费网站照片 | 大型黄色视频在线免费观看| 建设人人有责人人尽责人人享有的| 欧美av亚洲av综合av国产av| 亚洲,欧美精品.| 国产精品九九99| 在线看a的网站| 在线观看免费日韩欧美大片| 国产精品欧美亚洲77777| 成人国语在线视频| 国产精品成人在线| 麻豆av在线久日| 一级作爱视频免费观看| 69av精品久久久久久| a级毛片在线看网站| 久久中文字幕一级| 国产av又大| 一级毛片女人18水好多| 精品国产超薄肉色丝袜足j| 国产欧美日韩精品亚洲av| 国产国语露脸激情在线看| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产毛片av蜜桃av| 免费高清在线观看日韩| 亚洲人成伊人成综合网2020| 不卡av一区二区三区| 精品国产一区二区三区四区第35| 香蕉丝袜av| 亚洲精品国产一区二区精华液| 九色亚洲精品在线播放| 欧美人与性动交α欧美精品济南到| 国产真人三级小视频在线观看| av欧美777| 美女高潮到喷水免费观看| 男女之事视频高清在线观看| 超色免费av| 亚洲人成电影免费在线| 亚洲av日韩在线播放| 精品熟女少妇八av免费久了| 大码成人一级视频| 成人18禁在线播放| 国产亚洲精品一区二区www | 激情视频va一区二区三区| 悠悠久久av| 国产不卡av网站在线观看| 黄色视频,在线免费观看| 在线观看日韩欧美| 成人18禁高潮啪啪吃奶动态图| 亚洲色图av天堂| 国产真人三级小视频在线观看| 亚洲精品中文字幕一二三四区| 丰满人妻熟妇乱又伦精品不卡| 激情在线观看视频在线高清 | 久久九九热精品免费| 亚洲精品久久午夜乱码| 亚洲五月色婷婷综合| 精品国产一区二区三区四区第35| 人人妻人人澡人人爽人人夜夜| 在线观看日韩欧美| 免费在线观看日本一区| 久久久久久久国产电影| 亚洲av成人av| 性少妇av在线| 精品久久久久久电影网| 精品国产乱子伦一区二区三区| 老司机靠b影院| 国产精品二区激情视频| 中国美女看黄片| 欧美国产精品一级二级三级| 精品卡一卡二卡四卡免费| 十八禁高潮呻吟视频| 午夜福利免费观看在线| 俄罗斯特黄特色一大片| 久久精品国产亚洲av高清一级| 久久精品国产99精品国产亚洲性色 | 99香蕉大伊视频| 午夜免费鲁丝| 国产一区二区三区视频了| 丰满的人妻完整版| а√天堂www在线а√下载 | 12—13女人毛片做爰片一| 99热只有精品国产| 黄色视频,在线免费观看| 精品国产美女av久久久久小说| 国产精品久久久久久人妻精品电影| 91麻豆av在线| 精品福利观看| 欧美精品高潮呻吟av久久| 久久人妻福利社区极品人妻图片| 亚洲欧美激情综合另类| 亚洲一区二区三区不卡视频| 精品亚洲成a人片在线观看| av有码第一页| 久久精品国产99精品国产亚洲性色 | 午夜两性在线视频| 精品久久久久久,| 午夜福利,免费看| 亚洲欧美激情综合另类| 国产乱人伦免费视频| 80岁老熟妇乱子伦牲交| 亚洲五月婷婷丁香| 91老司机精品| 老熟妇仑乱视频hdxx| 91成年电影在线观看| 在线av久久热| 又紧又爽又黄一区二区| 天堂√8在线中文| av片东京热男人的天堂| 色老头精品视频在线观看| 亚洲男人天堂网一区| 无人区码免费观看不卡| 超碰成人久久| 精品国产美女av久久久久小说| 欧美日韩福利视频一区二区| 在线观看免费视频日本深夜| 99在线人妻在线中文字幕 | av天堂在线播放| 久久亚洲真实| 一级毛片女人18水好多| 天天躁狠狠躁夜夜躁狠狠躁| 色精品久久人妻99蜜桃| 99精品在免费线老司机午夜| 91成人精品电影| 天天操日日干夜夜撸| 色尼玛亚洲综合影院| 少妇粗大呻吟视频| 99精国产麻豆久久婷婷| 50天的宝宝边吃奶边哭怎么回事| 国产野战对白在线观看| 精品视频人人做人人爽| 欧美精品一区二区免费开放| 欧美精品啪啪一区二区三区| 国产欧美日韩一区二区三区在线| 国产有黄有色有爽视频| 免费av中文字幕在线| 久久人人97超碰香蕉20202| 黄频高清免费视频| 国内久久婷婷六月综合欲色啪| 欧美 日韩 精品 国产| 国内毛片毛片毛片毛片毛片| 国产精品 欧美亚洲| 飞空精品影院首页| 中文字幕另类日韩欧美亚洲嫩草| 777米奇影视久久| 少妇 在线观看| 国产免费现黄频在线看| 国产单亲对白刺激| 一区二区三区激情视频| 成在线人永久免费视频| 最新在线观看一区二区三区| 色老头精品视频在线观看| 国产欧美日韩精品亚洲av| 欧美丝袜亚洲另类 | 夫妻午夜视频| 成年人午夜在线观看视频| 国产又色又爽无遮挡免费看| 超色免费av| 看黄色毛片网站| 国产一区在线观看成人免费| 桃红色精品国产亚洲av| 色综合婷婷激情| 久久精品国产亚洲av香蕉五月 | 超碰97精品在线观看| 亚洲九九香蕉| 在线播放国产精品三级| 黑人巨大精品欧美一区二区蜜桃| 亚洲av熟女| 亚洲精品久久午夜乱码| 久久国产乱子伦精品免费另类| 桃红色精品国产亚洲av| 成人永久免费在线观看视频| a级片在线免费高清观看视频| 欧美日韩一级在线毛片| 精品国内亚洲2022精品成人 | 亚洲精品粉嫩美女一区| 精品无人区乱码1区二区| 麻豆成人av在线观看| 91av网站免费观看| 精品少妇久久久久久888优播| 他把我摸到了高潮在线观看| x7x7x7水蜜桃| 一区二区三区激情视频| 亚洲欧美一区二区三区黑人| 一二三四在线观看免费中文在| 国产精品二区激情视频| av网站在线播放免费| 久久久精品区二区三区| 大香蕉久久网| 亚洲国产中文字幕在线视频| 国产高清视频在线播放一区| 中文字幕人妻丝袜制服| 大码成人一级视频| 国产精品久久久久久精品古装| 成人亚洲精品一区在线观看| e午夜精品久久久久久久| 女人爽到高潮嗷嗷叫在线视频| 在线视频色国产色| 国产av又大| a级毛片黄视频| 色94色欧美一区二区| 免费av中文字幕在线| av天堂在线播放| 黄片大片在线免费观看| 日韩免费高清中文字幕av| 亚洲伊人色综图| 国产野战对白在线观看| 久久亚洲真实| 怎么达到女性高潮| 国产一卡二卡三卡精品| 法律面前人人平等表现在哪些方面| 村上凉子中文字幕在线| 亚洲av第一区精品v没综合| 女人精品久久久久毛片| 老司机深夜福利视频在线观看| 日本欧美视频一区| 在线观看日韩欧美| 色播在线永久视频| 中文亚洲av片在线观看爽 | 高清av免费在线| 亚洲熟女毛片儿| 捣出白浆h1v1| 欧美人与性动交α欧美精品济南到| 美女福利国产在线| 国产成人欧美| 国产精品.久久久| 欧美国产精品一级二级三级| 亚洲五月色婷婷综合| 日本精品一区二区三区蜜桃| 亚洲人成电影免费在线| 亚洲一区二区三区不卡视频| 国产蜜桃级精品一区二区三区 | 日韩中文字幕欧美一区二区| 久久久久久久午夜电影 | 久久久国产成人免费| 国产精品电影一区二区三区 | 亚洲一区高清亚洲精品| 国产无遮挡羞羞视频在线观看| 亚洲av日韩在线播放| 黑人欧美特级aaaaaa片| 亚洲精品中文字幕在线视频| 久久亚洲真实| 最新美女视频免费是黄的| 一区二区日韩欧美中文字幕| 18禁裸乳无遮挡免费网站照片 | 女性被躁到高潮视频| 欧美日韩中文字幕国产精品一区二区三区 | 午夜成年电影在线免费观看| 国产精品99久久99久久久不卡| 亚洲精品一二三| 久久香蕉激情| 免费在线观看视频国产中文字幕亚洲| 天堂动漫精品| 在线观看免费日韩欧美大片| 久久人妻av系列| 国产精品久久久av美女十八| 国产一区二区三区综合在线观看| 在线av久久热| 91字幕亚洲| 精品福利永久在线观看| 亚洲成人手机| 人人妻,人人澡人人爽秒播| 性色av乱码一区二区三区2| 大片电影免费在线观看免费| 亚洲精品一卡2卡三卡4卡5卡| av有码第一页| 成人三级做爰电影| 国产精品99久久99久久久不卡| 国产一区二区激情短视频| 女人精品久久久久毛片| 亚洲专区中文字幕在线| 免费观看a级毛片全部| 午夜视频精品福利| 欧美乱妇无乱码| 精品少妇一区二区三区视频日本电影| 国产精品久久视频播放| av免费在线观看网站| 一级,二级,三级黄色视频| 亚洲欧美日韩另类电影网站| 欧美精品人与动牲交sv欧美| 在线观看免费高清a一片| 视频区欧美日本亚洲| 天堂√8在线中文| 丝瓜视频免费看黄片| 久久99一区二区三区| 亚洲欧美日韩高清在线视频| 精品亚洲成a人片在线观看| a级毛片在线看网站| 美女福利国产在线| 日韩欧美一区视频在线观看| 午夜精品久久久久久毛片777| av免费在线观看网站| 中亚洲国语对白在线视频| 国产激情欧美一区二区| 性少妇av在线| 久久久久久人人人人人| 51午夜福利影视在线观看| www.熟女人妻精品国产| 老熟妇仑乱视频hdxx| 久久香蕉激情| 久久久久久人人人人人| 精品熟女少妇八av免费久了| 一进一出抽搐gif免费好疼 | 无人区码免费观看不卡| 伊人久久大香线蕉亚洲五| 国产国语露脸激情在线看| 日韩免费av在线播放| 日本一区二区免费在线视频| 变态另类成人亚洲欧美熟女 | 欧美精品av麻豆av| 男人舔女人的私密视频| 久久久久久久久久久久大奶| 欧美最黄视频在线播放免费 | 在线观看舔阴道视频| 欧美黄色片欧美黄色片| √禁漫天堂资源中文www| 成人国产一区最新在线观看| 18禁裸乳无遮挡免费网站照片 | 国产在视频线精品| 国产伦人伦偷精品视频| 9191精品国产免费久久| 久久精品人人爽人人爽视色| 精品高清国产在线一区| 岛国毛片在线播放| 国产99久久九九免费精品| 亚洲视频免费观看视频| 亚洲精品国产一区二区精华液| 亚洲欧美色中文字幕在线| 女警被强在线播放| 免费在线观看影片大全网站| 亚洲国产欧美一区二区综合| 国产精品乱码一区二三区的特点 | 一区二区三区国产精品乱码| 久久精品国产清高在天天线| 99精品久久久久人妻精品| 国产视频一区二区在线看| 99国产综合亚洲精品| 一级毛片高清免费大全| 极品教师在线免费播放| 亚洲专区字幕在线| 国产成人精品无人区| 黄色毛片三级朝国网站| 成人免费观看视频高清| 久久人妻av系列| 亚洲av成人av| 亚洲午夜理论影院| а√天堂www在线а√下载 | 午夜福利乱码中文字幕| 精品午夜福利视频在线观看一区| 精品一区二区三卡| 国内久久婷婷六月综合欲色啪| 国产精品二区激情视频| 丝袜美足系列| 黄色视频,在线免费观看| 女性被躁到高潮视频| 一级,二级,三级黄色视频| 美女国产高潮福利片在线看| 精品免费久久久久久久清纯 | 国内毛片毛片毛片毛片毛片| 欧美人与性动交α欧美软件| 99国产精品免费福利视频| 精品少妇久久久久久888优播| 久久久久久久午夜电影 | 我的亚洲天堂| 国产精品九九99| 欧美最黄视频在线播放免费 | 久热这里只有精品99| 丰满人妻熟妇乱又伦精品不卡| 国产精品免费一区二区三区在线 | 欧美av亚洲av综合av国产av| 午夜成年电影在线免费观看| 久久精品国产清高在天天线| 亚洲精品久久午夜乱码| 在线天堂中文资源库| 国产aⅴ精品一区二区三区波| 一区福利在线观看| 国产aⅴ精品一区二区三区波| 亚洲熟妇熟女久久| av电影中文网址| 欧美日韩乱码在线| 亚洲五月色婷婷综合| 亚洲国产精品合色在线| 黄片大片在线免费观看| 男人舔女人的私密视频| 精品久久蜜臀av无| 国产免费av片在线观看野外av| 精品国产乱子伦一区二区三区| 亚洲成av片中文字幕在线观看| 亚洲午夜理论影院| 一级片免费观看大全| 国产精品免费视频内射| 99精品久久久久人妻精品| 女人被躁到高潮嗷嗷叫费观| 国产亚洲欧美98| 1024香蕉在线观看| 看片在线看免费视频| 久久99一区二区三区| 黄色 视频免费看| 国产97色在线日韩免费| 欧美日韩亚洲高清精品| 久久久久久久久久久久大奶| 亚洲成人手机| 91成年电影在线观看| 窝窝影院91人妻| 亚洲av成人av| 久久青草综合色| 久久精品人人爽人人爽视色| 99国产精品99久久久久| 天天添夜夜摸| 久久精品亚洲熟妇少妇任你| 少妇裸体淫交视频免费看高清 | 亚洲国产欧美日韩在线播放| 亚洲色图 男人天堂 中文字幕| 精品久久久精品久久久| 精品免费久久久久久久清纯 | 欧美精品一区二区免费开放| 精品亚洲成a人片在线观看| av一本久久久久| 精品人妻熟女毛片av久久网站| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品自拍成人| 精品电影一区二区在线| 国产亚洲欧美精品永久| 精品人妻熟女毛片av久久网站| 99国产精品一区二区三区| 青草久久国产| 国产视频一区二区在线看|