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

    復(fù)雜水網(wǎng)中船舶航行路徑多目標(biāo)規(guī)劃模型

    2017-10-30 08:43:10周春輝陳銘章何義才徐周華文元橋
    中國(guó)航海 2017年2期
    關(guān)鍵詞:水網(wǎng)航道航行

    周春輝,陳銘章,何義才,徐周華,文元橋

    (1.武漢理工大學(xué) 航運(yùn)學(xué)院,武漢 430063;2.湖北省內(nèi)河航運(yùn)技術(shù)重點(diǎn)實(shí)驗(yàn)室,武漢 430079;3.武漢大學(xué) 測(cè)繪學(xué)院,武漢 430070)

    復(fù)雜水網(wǎng)中船舶航行路徑多目標(biāo)規(guī)劃模型

    周春輝1,2,陳銘章1,2,何義才1,3,徐周華1,文元橋1

    (1.武漢理工大學(xué) 航運(yùn)學(xué)院,武漢 430063;2.湖北省內(nèi)河航運(yùn)技術(shù)重點(diǎn)實(shí)驗(yàn)室,武漢 430079;3.武漢大學(xué) 測(cè)繪學(xué)院,武漢 430070)

    在深入研究船舶交通流理論的基礎(chǔ)上,結(jié)合內(nèi)河復(fù)雜水網(wǎng)獨(dú)有的特性,構(gòu)建復(fù)雜水網(wǎng)條件下的船舶航行路徑多目標(biāo)規(guī)劃模型;運(yùn)用層次分析法獲取多目標(biāo)規(guī)劃模型各子目標(biāo)的權(quán)重;利用線性加權(quán)法和Floyd算法求解得到綜合最優(yōu)航行路徑。該模型不僅綜合考慮了影響航行路徑選擇的各種因素,而且對(duì)其進(jìn)行了合理的量化,并通過算例證明了模型的合理性和有效性。

    復(fù)雜水網(wǎng);船舶交通流;最優(yōu)航行路徑;多目標(biāo)規(guī)劃;線性加權(quán)法;Floyd算法

    近年來,隨著水上通航需求不斷上升,國(guó)內(nèi)部分通航條件受限的河段(尤其是江蘇、浙江一帶的水網(wǎng)區(qū)域)擁擠堵塞狀況日益嚴(yán)重。無論是從船舶自主航行路徑選擇的角度還是從內(nèi)河航運(yùn)管理部門調(diào)度需要的角度考慮,都迫切需要整個(gè)水網(wǎng)層面的航行路徑規(guī)劃方案。目前國(guó)內(nèi)外相關(guān)學(xué)者已在這方面進(jìn)行一些初步研究。

    李文歡[1]對(duì)集裝箱班輪運(yùn)輸?shù)暮骄€網(wǎng)絡(luò)優(yōu)化問題進(jìn)行深入研究,提出利用圖論中的旅行商問題模型及整數(shù)規(guī)劃等方法求解集裝箱船最優(yōu)航行路徑,并依據(jù)長(zhǎng)江中下游集裝箱班輪運(yùn)輸航線網(wǎng)絡(luò)的統(tǒng)計(jì)信息設(shè)計(jì)優(yōu)化算法,對(duì)模型進(jìn)行驗(yàn)證。劉建秋[2]從系統(tǒng)總收益最大化的角度出發(fā),綜合考慮模糊變量和隨機(jī)變量,建立不確定環(huán)境下的支線集裝箱班輪優(yōu)化網(wǎng)絡(luò)模型,并以長(zhǎng)三角某運(yùn)輸公司為例,采用蟻群算法求取網(wǎng)絡(luò)運(yùn)輸路徑。PFLIEGL[3]針對(duì)內(nèi)河水網(wǎng)中河道及船閘通航條件的動(dòng)態(tài)變化使調(diào)度變得困難,提出在澳大利亞Danube河上建立內(nèi)河信息服務(wù)系統(tǒng),以促進(jìn)水網(wǎng)中船舶的調(diào)度與管理。

    總的來說,目前水網(wǎng)層面上的船舶航行路徑選擇和船舶聯(lián)合調(diào)度尚處在研究探索階段,主要以信息共享為前提,一定程度上運(yùn)用預(yù)調(diào)度協(xié)調(diào)機(jī)制來提高航道的利用率,尚沒有實(shí)用的調(diào)度和路徑規(guī)劃模型,且相關(guān)研究主要集中在涉及船閘的船舶調(diào)度方面,很少涉及復(fù)雜水網(wǎng)條件下的船舶調(diào)度和航行路徑規(guī)劃。因此,研究復(fù)雜水網(wǎng)中船舶航行路徑多目標(biāo)規(guī)劃模型不僅對(duì)船舶自主航行路徑選擇的研究具有一定的指導(dǎo)意義,而且對(duì)內(nèi)河船舶管理調(diào)度部門管理調(diào)度船舶具有一定的參考價(jià)值。

    1 航行路徑選擇的多目標(biāo)規(guī)劃模型

    1.1問題分析

    影響船舶航行路徑選擇的因素可分為主因素觀和客觀因素2類,其中:主觀因素主要包括駕駛員的駕駛經(jīng)驗(yàn)、駕駛路徑選擇偏好及航行目的等;客觀因素主要包括航行總路程、總時(shí)間、擁擠程度及通行費(fèi)用等。因此,最優(yōu)航行路徑選擇是一個(gè)各種因素交互作用的多目標(biāo)規(guī)劃問題。[4]

    1.2目標(biāo)函數(shù)

    根據(jù)對(duì)最優(yōu)航行路徑選擇影響因素的分析,選擇總路程最小化、航道綜合時(shí)間阻抗最小化、擁擠程度最小化和通行費(fèi)用最小化作為模型子目標(biāo)。

    1.2.1總路程最小化

    (1)

    1.2.2航道綜合時(shí)間阻抗最小化

    與道路交通中的交通阻抗不同,內(nèi)河航道的交通時(shí)間阻抗由航道基本段和船閘時(shí)間阻抗2部分組成。根據(jù)水網(wǎng)交通分配路徑的選擇理論,水網(wǎng)時(shí)間阻抗函數(shù)模型[5]為

    (2)

    文獻(xiàn)[5]基于船舶交通流理論,結(jié)合交通流特性中的流量、速度和密度等3個(gè)參數(shù)提出3個(gè)參數(shù)間的關(guān)系模型,其中速度-流量的關(guān)系為

    (3)

    式(3)中:q為船舶流流量,艘/h;v為船舶流航行速度,km/h;vf為自由流速度,km/h;k為船舶流密度,艘/km;kf為自由流狀態(tài)下航道中的最大船舶密度,艘/km;kj為阻塞密度,艘/km;ρ為波速系數(shù),0≤ρ≤1。

    基于船舶交通流理論,假設(shè)航道段a上船舶的平均航行時(shí)間Ra(q)與航道長(zhǎng)度la及船舶航行的平均速度va(q)有關(guān),即

    (4)

    根據(jù)交通流理論,假設(shè)航道上船舶的到達(dá)分布符合泊松分布[6],λb為船閘b處船舶的到達(dá)率,E(T)為過閘服務(wù)時(shí)間,Var[T]為服務(wù)時(shí)間的方差,得到服務(wù)強(qiáng)度ρb=λb°E(T)。由Pollaczck-Khintchine公式可得M/G/1模型系統(tǒng)內(nèi)船閘b的平均排隊(duì)等待船舶數(shù)Nb為

    (5)

    (6)

    綜上所述,假設(shè)船舶傾向于選擇航行總阻抗最小的路徑,最優(yōu)航行路徑選擇的目標(biāo)函數(shù)為T2,則聯(lián)立式(2)~式(6)可得

    (7)

    1.2.3擁擠程度最小化

    擁擠程度S[7]可用航道中交通流的飽和度來衡量。假設(shè)航道段a上的設(shè)計(jì)通行能力為ca,實(shí)際交通流量為fa,則飽和度sa=fa/ca。假設(shè)船舶傾向于選擇從起點(diǎn)到終點(diǎn)擁擠程度最小的航行路徑,則最優(yōu)航行路徑選擇的目標(biāo)函數(shù)T3為

    (8)

    1.2.4通行費(fèi)用最小化

    船舶的通行費(fèi)用主要包括燃料費(fèi)和過閘費(fèi)2部分。由于船舶在內(nèi)河中航行時(shí)周圍環(huán)境穩(wěn)定,風(fēng)浪影響幾乎可以忽略不計(jì),因此在船舶載重不變的情況下可假設(shè)航行所消耗的燃料費(fèi)與航程成正比。由此,在通行費(fèi)用最小化的目標(biāo)下,最優(yōu)航行路徑選擇的目標(biāo)函數(shù)T4為

    (9)

    1.3航行路徑選擇的多目標(biāo)規(guī)劃模型

    聯(lián)立式(1)、式(7)~式(9)得最優(yōu)航行路徑選擇的多目標(biāo)規(guī)劃模型為

    (10)

    1.4模型求解

    對(duì)于多目標(biāo)規(guī)劃問題,可用線性加權(quán)法、約束法、效用系數(shù)法和評(píng)價(jià)函數(shù)法等方法進(jìn)行求解。根據(jù)所建立的多目標(biāo)規(guī)劃模型和各分目標(biāo)的特點(diǎn),采取線性加權(quán)法進(jìn)行求解??紤]到上述多目標(biāo)規(guī)劃模型的4個(gè)子目標(biāo)的量綱不同,采用極差化法對(duì)變量進(jìn)行無量綱化處理。[8]

    (11)

    (12)

    利用線性加權(quán)法對(duì)多目標(biāo)規(guī)劃模型進(jìn)行求解,得

    (13)

    式(13)中:λ1,λ2,λ3,λ4分別為4個(gè)子目標(biāo)的權(quán)重系數(shù),既可采用層次分析法求取,也可根據(jù)船舶對(duì)各個(gè)子目標(biāo)的偏好程度或船舶調(diào)度者的調(diào)度意愿來主觀確定。

    2 算例分析

    2.1層次分析法確定權(quán)重

    由于建立的多目標(biāo)規(guī)劃模型有4個(gè)子目標(biāo),而每個(gè)子目標(biāo)的重要性無論是對(duì)調(diào)度者還是對(duì)船舶駕駛員而言都是不同的,因此根據(jù)層次分析法[9]確定4個(gè)子目標(biāo)的權(quán)重系數(shù)λ1,λ2,λ3和λ4的值。

    2.1.1建立評(píng)價(jià)體系

    根據(jù)層析分析法評(píng)價(jià)體系構(gòu)建規(guī)則,構(gòu)建最優(yōu)航行路徑評(píng)價(jià)體系(見圖1)。

    圖1 航行路徑選擇評(píng)價(jià)體系

    2.1.2確定準(zhǔn)則層判斷矩陣

    利用表1所示的1~9標(biāo)度法對(duì)準(zhǔn)則層中的4個(gè)因素(總航程B1,綜合時(shí)間阻抗B2,擁擠程度B3,通行費(fèi)用B4)進(jìn)行成對(duì)比較。

    表1 1~9標(biāo)度法的意義

    2.1.3最大特征根與一致性檢驗(yàn)

    求解可得最大特征根λmax=4.12,特征向量(權(quán)重系數(shù))λi=(0.05,0.57,0.26,0.12)。一致性指標(biāo)計(jì)算式為

    (14)

    查找Saaty給出的平均隨機(jī)一致性指標(biāo)RI[9],計(jì)算一致性比例可得

    (15)

    由于0.043<0.1,因此判斷矩陣的一致性是可以接受的,權(quán)重系數(shù)λi=(0.05,0.57,0.26,0.12)。

    2.2確定鄰近點(diǎn)之間的“距離”

    由式(13)可得4個(gè)子目標(biāo)的取值對(duì)于船舶或船舶調(diào)度者來說均是越小越好,因此式(13)的取值可認(rèn)為是“距離”。根據(jù)Floyd算法[10]的求解條件,需先確定鄰近點(diǎn)之間的距離。在能獲得數(shù)據(jù)和確定λ1~λ4的情況下,鄰近點(diǎn)之間的“距離”容易求取。假設(shè)模型4個(gè)子目標(biāo)的值計(jì)入4個(gè)數(shù)組(t1,t2,t3,t4)中,則多目標(biāo)規(guī)劃模型的值(即鄰近點(diǎn)之間的“距離”)計(jì)入數(shù)組t中,t=λ1t1+λ2t2+λ3t3+λ4t4。

    2.3計(jì)算結(jié)果與分析

    假設(shè)兩點(diǎn)之間沒有直接連接或沒有連接,則兩點(diǎn)之間的距離定義為無窮大,通過確定t1,t2,t3,t4的值,并取λi=(0.05,0.57,0.26,0.12),得鄰近點(diǎn)之間的“距離”t(即賦權(quán)圖的權(quán)矩陣,見圖2)。

    圖2 水網(wǎng)示意及鄰近點(diǎn)之間的“距離”

    利用Floyd算法進(jìn)行求解,上述t即為賦權(quán)圖的權(quán)矩陣,得到水網(wǎng)中任意兩點(diǎn)最短路徑的長(zhǎng)度d(i,j)和第i個(gè)點(diǎn)到第j個(gè)點(diǎn)最短路徑中第i個(gè)點(diǎn)后繼點(diǎn)的編號(hào)path(i,j),見表2。

    表2 i點(diǎn)到j(luò)點(diǎn)的最短路徑中i點(diǎn)后繼點(diǎn)編號(hào)

    至此,即求出圖1網(wǎng)絡(luò)中任意兩點(diǎn)之間的最短路徑及其最短“距離”。例如,從表2中可查到點(diǎn)1到點(diǎn)8的最短路徑path(1,8)為1—4—7—8。但是,該結(jié)果不是唯一的,會(huì)隨著各種因素的變化而變化。例如:船舶的偏好或船舶調(diào)度者的意愿不同而使得λ1,λ2,λ3及λ4的取值不同;航道改道會(huì)導(dǎo)致總航程發(fā)生變化;航行時(shí)間及緊張程度不同會(huì)導(dǎo)致理解阻抗不同;隨著時(shí)間的推移,船舶密度發(fā)生變化會(huì)導(dǎo)致?lián)頂D程度不同;油價(jià)或過閘費(fèi)用的變化導(dǎo)致通行費(fèi)用上升。

    3 結(jié)束語

    針對(duì)最優(yōu)航行路徑選擇問題的多樣性特點(diǎn),將該問題轉(zhuǎn)化成一個(gè)具有4個(gè)子目標(biāo)的多目標(biāo)規(guī)劃模型,并利用線性加權(quán)法和Floyd算法進(jìn)行求解。通過分析,在獲取航道資料、過閘費(fèi)用、油價(jià)、水網(wǎng)中船舶的實(shí)時(shí)分布、船舶本身航行時(shí)間的緊迫度及航行目的等信息的條件下,可根據(jù)建立的多目標(biāo)規(guī)劃模型及其求解算法求得船舶的最優(yōu)航行路徑,實(shí)現(xiàn)復(fù)雜水網(wǎng)條件下船舶的自主航行路徑選擇和調(diào)度部門對(duì)船舶的宏觀調(diào)度。

    但是,運(yùn)用建立的復(fù)雜水網(wǎng)條件下船舶聯(lián)合調(diào)度多目標(biāo)規(guī)劃模型的前提是獲得內(nèi)河航運(yùn)調(diào)度部門提供的航道信息和船舶交通流量信息,除了考慮航行路徑路程、航行路徑阻抗、航行路徑擁擠程度及航行費(fèi)用等4項(xiàng)評(píng)價(jià)指標(biāo)之外,還可考慮服務(wù)水平、上水與下水等因素。此外,建立一個(gè)更加適合復(fù)雜、動(dòng)態(tài)變化的水網(wǎng)交通調(diào)度模型是進(jìn)一步研究的方向。

    [1] 李文歡.集裝箱班輪運(yùn)輸?shù)暮骄€網(wǎng)絡(luò)優(yōu)化研究[D].上海:上海海事大學(xué),2007.

    [2] 劉建秋.不確定環(huán)境下支線集裝箱班輪網(wǎng)絡(luò)優(yōu)化研究 [D].大連:大連海事大學(xué),2010.

    [3] PFLIEGL R.Innovative Application for Dynamic Navigational Support and Transport Management on Inland Waterways:Experience from a Research Project on the Danube River[J].Journal of Transportation Research Board,2001,1763:85-89.

    [4] 宗傳苓,李相勇,王英濤.出行前路徑選擇的多目標(biāo)規(guī)劃模型[J].交通運(yùn)輸系統(tǒng)工程與信息,2005,5(6):58-61.

    [5] 朱俊.基于船舶流的交通時(shí)間阻抗模型[J].武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版),2010,34(3):591-594.

    [6] 孫榮恒,李建平.排隊(duì)論基礎(chǔ)[M].北京:科學(xué)出版社,2002.

    [7] 達(dá)慶東,張國(guó)伍.交通擁擠定量分析方法[J].交通運(yùn)輸系統(tǒng)工程信息,2002,2(4):45-48.

    [8] 呂寶軍.機(jī)動(dòng)車出行路線選擇及軟件實(shí)現(xiàn)[D].天津:天津大學(xué),2012:45-48.

    [9] 司守奎,孫璽菁.數(shù)學(xué)建模算法與應(yīng)用[M].北京:國(guó)防工業(yè)出版社,2011:167-171.

    [10] 張秀蘭,林峰.數(shù)學(xué)建模與實(shí)驗(yàn)[M].北京:化學(xué)工業(yè)出版社,2013:101-108.

    Multi-ObjectiveProgrammingModelofShipSailingPathforComplexWaterNetwork

    ZHOUChunhui1,2,CHENMingzhang1,2,HEYicai1,3,XUZhouhua1,WENYuanqiao1

    (1.School of Navigation,Wuhan University of Technology,Wuhan 430063,China;2.Hubei Inland Shipping Technology Key Laboratory,Wuhan 430070,China;3.School of Geodesy and Geomatics,Wuhan University,Wuhan 430079,China)

    On the basis of the vessel traffic flow theory,taking the unique characteristics of the inland water net,the multi-objective programming model of ship sailing path is established under the condition of complex water network,and the AHP is used to get the weights of each sub-goal of the multi-objective programming model.By means of the linear weighting method and Floyd algorithm the overall optimal sailing path is decided.This model,taking the factors that affect the sailing path selection into account quantitatively,is verified by a numerical example.

    complex water network; ship traffic flow; optimal sailing path; multi-objective programming; linear weighting; Floyd algorithm

    U697.1

    A

    2017-01-11

    國(guó)家自然科學(xué)基金(51209166;51679180);湖北省交通廳科技計(jì)劃項(xiàng)目(2014-721-4-2);武漢理工大學(xué)自主創(chuàng)新研究基金(2017-HY-B1-07)

    周春輝(1978—),男,湖北松滋人,博士,副教授,從事交通安全仿真方面研究。E-mail:chunhui@whut.edu.cn

    陳銘章(1996—),男,湖北黃岡人,學(xué)士,從事水路交通運(yùn)輸研究。E-mail:984393738@qq.com

    1000-4653(2016)04-0078-05

    猜你喜歡
    水網(wǎng)航道航行
    水利部:到2035 年基本形成國(guó)家水網(wǎng)總體格局
    水網(wǎng)藻生長(zhǎng)及鈾對(duì)其生長(zhǎng)影響的模型構(gòu)建
    到慧骃國(guó)的航行
    ●山西加快推進(jìn)縣域小水網(wǎng)配套建設(shè)
    小舟在河上航行
    新航道
    留學(xué)(2017年5期)2017-03-29 03:03:54
    航行
    青年歌聲(2017年6期)2017-03-13 00:57:56
    山西省大水網(wǎng)工程建設(shè)如火如荼
    山西水利(2017年1期)2017-03-07 08:50:31
    我國(guó)首條40萬噸級(jí)航道正式開建
    水道港口(2014年1期)2014-04-27 14:14:38
    我國(guó)首條40萬噸超大型航道將迎來開工建設(shè)
    水道港口(2013年3期)2013-04-27 14:14:28
    有码 亚洲区| 国产精品嫩草影院av在线观看| 中文字幕熟女人妻在线| 亚洲av电影在线观看一区二区三区 | 国产极品天堂在线| 精品一区二区三区人妻视频| 神马国产精品三级电影在线观看| 日韩人妻高清精品专区| 久久精品国产亚洲av涩爱| av黄色大香蕉| 久久久国产成人免费| 91久久精品国产一区二区成人| 国产精品,欧美在线| 色网站视频免费| 欧美一级a爱片免费观看看| 日韩 亚洲 欧美在线| 99国产精品一区二区蜜桃av| 亚洲av中文av极速乱| 国产又色又爽无遮挡免| 啦啦啦啦在线视频资源| 岛国在线免费视频观看| 久久久久久久久大av| av在线蜜桃| 97超碰精品成人国产| 老司机福利观看| 国语对白做爰xxxⅹ性视频网站| 国产黄a三级三级三级人| 亚洲国产欧美人成| 舔av片在线| 久久6这里有精品| 日本免费a在线| 欧美zozozo另类| 免费观看的影片在线观看| 日韩一区二区三区影片| 欧美另类亚洲清纯唯美| 免费在线观看成人毛片| 97热精品久久久久久| 午夜亚洲福利在线播放| 在线免费观看的www视频| 高清日韩中文字幕在线| 国产私拍福利视频在线观看| 久久这里有精品视频免费| 国产在线男女| 国产激情偷乱视频一区二区| kizo精华| 欧美不卡视频在线免费观看| 免费av不卡在线播放| 国产乱人偷精品视频| 国产中年淑女户外野战色| 成人午夜高清在线视频| 亚洲一级一片aⅴ在线观看| 国产免费一级a男人的天堂| 日韩视频在线欧美| 国产亚洲一区二区精品| 久久精品夜色国产| 插阴视频在线观看视频| 国产精品野战在线观看| 国产欧美日韩精品一区二区| 美女xxoo啪啪120秒动态图| 特级一级黄色大片| 国产又色又爽无遮挡免| 久久精品人妻少妇| 99热精品在线国产| 伊人久久精品亚洲午夜| 国产三级中文精品| 亚洲美女搞黄在线观看| 色噜噜av男人的天堂激情| 国产 一区 欧美 日韩| 在线观看美女被高潮喷水网站| 国产老妇女一区| 中文字幕制服av| 人人妻人人澡欧美一区二区| 性插视频无遮挡在线免费观看| 成人亚洲欧美一区二区av| 女人久久www免费人成看片 | 小说图片视频综合网站| 日韩av在线免费看完整版不卡| 七月丁香在线播放| 偷拍熟女少妇极品色| 国产v大片淫在线免费观看| 又爽又黄无遮挡网站| 真实男女啪啪啪动态图| videos熟女内射| 伊人久久精品亚洲午夜| 黑人高潮一二区| 久久人妻av系列| 91av网一区二区| 欧美色视频一区免费| www.av在线官网国产| 日日撸夜夜添| 国产探花在线观看一区二区| 国产 一区精品| 建设人人有责人人尽责人人享有的 | 网址你懂的国产日韩在线| 有码 亚洲区| 国产午夜精品一二区理论片| 国产亚洲午夜精品一区二区久久 | 欧美zozozo另类| 免费一级毛片在线播放高清视频| 精品免费久久久久久久清纯| 黄色一级大片看看| 成人av在线播放网站| 日本黄大片高清| av线在线观看网站| av免费在线看不卡| 国产精品久久久久久久久免| 99久久中文字幕三级久久日本| 亚洲人成网站在线播| 久久精品国产自在天天线| 小蜜桃在线观看免费完整版高清| 91精品伊人久久大香线蕉| 久久亚洲精品不卡| 国产毛片a区久久久久| 欧美成人一区二区免费高清观看| 看片在线看免费视频| 亚洲国产欧洲综合997久久,| 国产毛片a区久久久久| 国产一级毛片在线| 久久久精品94久久精品| 在线观看一区二区三区| 国产午夜福利久久久久久| 亚洲人成网站在线观看播放| 欧美成人精品欧美一级黄| 日韩av在线大香蕉| 亚洲av男天堂| 美女大奶头视频| 国产黄片美女视频| 国产精品电影一区二区三区| 一区二区三区四区激情视频| 午夜免费男女啪啪视频观看| 69av精品久久久久久| 亚洲精品国产成人久久av| 国产精品,欧美在线| 久久久久久九九精品二区国产| 熟妇人妻久久中文字幕3abv| 真实男女啪啪啪动态图| 久久精品久久久久久久性| 五月玫瑰六月丁香| 国产午夜精品久久久久久一区二区三区| 欧美不卡视频在线免费观看| 久久精品熟女亚洲av麻豆精品 | 九九爱精品视频在线观看| 亚洲国产精品成人久久小说| 国产男人的电影天堂91| 七月丁香在线播放| 婷婷色综合大香蕉| av播播在线观看一区| 女人久久www免费人成看片 | 精品酒店卫生间| 国产免费男女视频| 久久久久免费精品人妻一区二区| 成年av动漫网址| 久久精品国产99精品国产亚洲性色| 毛片一级片免费看久久久久| 丰满人妻一区二区三区视频av| 成人毛片60女人毛片免费| 九九在线视频观看精品| 精品熟女少妇av免费看| 日韩一区二区视频免费看| 国产真实乱freesex| 国产欧美另类精品又又久久亚洲欧美| 久久久久久九九精品二区国产| 91精品一卡2卡3卡4卡| 男女那种视频在线观看| 少妇熟女欧美另类| 国产真实乱freesex| 我要搜黄色片| 麻豆国产97在线/欧美| 成人毛片60女人毛片免费| 波多野结衣巨乳人妻| 最近2019中文字幕mv第一页| 在线观看美女被高潮喷水网站| 久久精品影院6| 婷婷六月久久综合丁香| 22中文网久久字幕| 高清午夜精品一区二区三区| 看免费成人av毛片| 色尼玛亚洲综合影院| 水蜜桃什么品种好| 搡老妇女老女人老熟妇| 日日撸夜夜添| 九草在线视频观看| 又黄又爽又刺激的免费视频.| 日本一二三区视频观看| 日本三级黄在线观看| 国产激情偷乱视频一区二区| 建设人人有责人人尽责人人享有的 | 亚洲欧美清纯卡通| 综合色丁香网| 国内精品美女久久久久久| 白带黄色成豆腐渣| 三级经典国产精品| 国产探花极品一区二区| 成人综合一区亚洲| 日本三级黄在线观看| 老司机影院成人| 亚洲va在线va天堂va国产| 水蜜桃什么品种好| 国模一区二区三区四区视频| 99热6这里只有精品| 国产高清三级在线| 欧美人与善性xxx| 免费av不卡在线播放| 精品人妻偷拍中文字幕| 黄片wwwwww| 变态另类丝袜制服| 卡戴珊不雅视频在线播放| 亚洲va在线va天堂va国产| 大香蕉97超碰在线| 久久人人爽人人爽人人片va| 波多野结衣高清无吗| 亚洲第一区二区三区不卡| 亚洲五月天丁香| 美女黄网站色视频| av又黄又爽大尺度在线免费看 | 国产又黄又爽又无遮挡在线| 国产成人91sexporn| 欧美性猛交黑人性爽| 丝袜美腿在线中文| 久久国内精品自在自线图片| 99热这里只有是精品在线观看| 国产69精品久久久久777片| 韩国高清视频一区二区三区| 国产精品电影一区二区三区| 国产精品综合久久久久久久免费| 国产av不卡久久| 午夜亚洲福利在线播放| 99久国产av精品国产电影| 欧美丝袜亚洲另类| 精品久久久久久久久久久久久| 国产又色又爽无遮挡免| 寂寞人妻少妇视频99o| 少妇被粗大猛烈的视频| 在线观看一区二区三区| 一个人免费在线观看电影| 久99久视频精品免费| 真实男女啪啪啪动态图| 少妇高潮的动态图| 亚洲av二区三区四区| 纵有疾风起免费观看全集完整版 | 一级二级三级毛片免费看| 日韩av在线大香蕉| 欧美日韩国产亚洲二区| 老司机影院成人| 亚洲va在线va天堂va国产| 久久精品久久精品一区二区三区| 国产精品女同一区二区软件| 六月丁香七月| 国产高清有码在线观看视频| 一二三四中文在线观看免费高清| 国产男人的电影天堂91| 午夜福利在线观看吧| 丝袜美腿在线中文| 亚洲欧美日韩卡通动漫| av国产久精品久网站免费入址| 国产精品一区www在线观看| 亚洲国产色片| 久久6这里有精品| www日本黄色视频网| 九色成人免费人妻av| 丰满乱子伦码专区| 免费大片18禁| av在线播放精品| 亚洲av电影在线观看一区二区三区 | 性色avwww在线观看| 国产精品日韩av在线免费观看| 国产不卡一卡二| 国产熟女欧美一区二区| 亚洲成av人片在线播放无| 亚洲欧美成人精品一区二区| 欧美激情在线99| 激情 狠狠 欧美| 在线免费观看的www视频| 亚洲精品,欧美精品| 身体一侧抽搐| 国产亚洲精品av在线| 国产久久久一区二区三区| 热99re8久久精品国产| 岛国在线免费视频观看| 舔av片在线| 婷婷六月久久综合丁香| 免费看a级黄色片| 成人欧美大片| av在线播放精品| 亚洲精品乱久久久久久| 成人三级黄色视频| 国产 一区精品| 99热网站在线观看| 国产精品一及| 青春草亚洲视频在线观看| 日韩av在线大香蕉| 男插女下体视频免费在线播放| 免费在线观看成人毛片| 国产麻豆成人av免费视频| 女人被狂操c到高潮| 桃色一区二区三区在线观看| 色尼玛亚洲综合影院| 特大巨黑吊av在线直播| 在线免费观看不下载黄p国产| 美女高潮的动态| 国产精品人妻久久久久久| 中文字幕av在线有码专区| 久久鲁丝午夜福利片| 亚洲国产日韩欧美精品在线观看| 中国美白少妇内射xxxbb| 亚洲成av人片在线播放无| av天堂中文字幕网| 亚洲第一区二区三区不卡| 日本黄大片高清| 91午夜精品亚洲一区二区三区| 午夜日本视频在线| 3wmmmm亚洲av在线观看| 婷婷六月久久综合丁香| 欧美区成人在线视频| 一级黄色大片毛片| 久久热精品热| 精华霜和精华液先用哪个| www日本黄色视频网| 国产成人免费观看mmmm| 舔av片在线| 好男人在线观看高清免费视频| 免费看美女性在线毛片视频| 国产精品精品国产色婷婷| 久久久精品欧美日韩精品| 国产黄色小视频在线观看| 久久国内精品自在自线图片| 久99久视频精品免费| 蜜桃亚洲精品一区二区三区| 最近最新中文字幕大全电影3| av黄色大香蕉| 晚上一个人看的免费电影| 黄色欧美视频在线观看| www.av在线官网国产| 青春草视频在线免费观看| www.av在线官网国产| 日韩人妻高清精品专区| 久久欧美精品欧美久久欧美| 国产乱人视频| 身体一侧抽搐| 国产精品综合久久久久久久免费| 日本午夜av视频| 亚洲av福利一区| 黑人高潮一二区| 亚洲av一区综合| 搡老妇女老女人老熟妇| 好男人在线观看高清免费视频| 国产单亲对白刺激| 啦啦啦观看免费观看视频高清| 日本黄色视频三级网站网址| 国产一区二区在线观看日韩| av在线播放精品| 一二三四中文在线观看免费高清| 卡戴珊不雅视频在线播放| 最近的中文字幕免费完整| 国产伦精品一区二区三区四那| 国产成人一区二区在线| 亚洲激情五月婷婷啪啪| 亚洲熟妇中文字幕五十中出| 亚洲国产最新在线播放| 99久久成人亚洲精品观看| 日本-黄色视频高清免费观看| 国产精品av视频在线免费观看| 两个人视频免费观看高清| 桃色一区二区三区在线观看| 女人久久www免费人成看片 | 成年女人永久免费观看视频| 国产精品一二三区在线看| 熟妇人妻久久中文字幕3abv| 只有这里有精品99| 国产av不卡久久| 桃色一区二区三区在线观看| 国产成人freesex在线| 亚洲欧美成人综合另类久久久 | 久久精品熟女亚洲av麻豆精品 | 乱码一卡2卡4卡精品| 国产一级毛片七仙女欲春2| 久久久久久大精品| 乱系列少妇在线播放| 老女人水多毛片| 精品不卡国产一区二区三区| 全区人妻精品视频| 神马国产精品三级电影在线观看| 床上黄色一级片| 男人和女人高潮做爰伦理| 亚洲人成网站高清观看| 久久精品91蜜桃| 国产精品,欧美在线| 午夜日本视频在线| 国产中年淑女户外野战色| 最近的中文字幕免费完整| 久久久国产成人精品二区| 丝袜喷水一区| 日韩在线高清观看一区二区三区| 人人妻人人澡人人爽人人夜夜 | 日韩人妻高清精品专区| 三级男女做爰猛烈吃奶摸视频| 欧美激情在线99| 日本黄色片子视频| 欧美日韩一区二区视频在线观看视频在线 | 一本一本综合久久| 国产又黄又爽又无遮挡在线| 插逼视频在线观看| eeuss影院久久| 日韩av在线大香蕉| 国产精品永久免费网站| 国产白丝娇喘喷水9色精品| 国产精品av视频在线免费观看| 亚洲高清免费不卡视频| 又爽又黄a免费视频| 熟女电影av网| 亚洲18禁久久av| 乱系列少妇在线播放| 中文资源天堂在线| 国产黄片美女视频| 久久久精品94久久精品| 女的被弄到高潮叫床怎么办| 国产三级中文精品| 插阴视频在线观看视频| 人人妻人人看人人澡| 日本wwww免费看| 简卡轻食公司| 色综合色国产| 亚洲av福利一区| 久久久国产成人免费| 韩国高清视频一区二区三区| 高清在线视频一区二区三区 | 日本一本二区三区精品| 国产精品野战在线观看| 成人欧美大片| 国产伦理片在线播放av一区| 午夜福利网站1000一区二区三区| 久久热精品热| 色吧在线观看| a级毛片免费高清观看在线播放| 一边摸一边抽搐一进一小说| 亚洲美女搞黄在线观看| 白带黄色成豆腐渣| 久久欧美精品欧美久久欧美| 亚洲熟妇中文字幕五十中出| 精品一区二区三区视频在线| 日本三级黄在线观看| 成人一区二区视频在线观看| 熟妇人妻久久中文字幕3abv| 日韩欧美 国产精品| 美女大奶头视频| 亚洲av男天堂| 一个人免费在线观看电影| 18禁在线播放成人免费| 91久久精品国产一区二区成人| 热99re8久久精品国产| 中文在线观看免费www的网站| 伦理电影大哥的女人| 亚洲自拍偷在线| 国产精品熟女久久久久浪| 亚洲色图av天堂| 成人特级av手机在线观看| 99久久成人亚洲精品观看| 亚洲欧美精品专区久久| 午夜福利在线观看吧| 又黄又爽又刺激的免费视频.| 一区二区三区高清视频在线| 亚洲av熟女| 一级毛片aaaaaa免费看小| 国产乱人视频| 中文字幕亚洲精品专区| 欧美日本视频| 22中文网久久字幕| 狠狠狠狠99中文字幕| 观看免费一级毛片| 日日啪夜夜撸| 亚洲婷婷狠狠爱综合网| 中文字幕久久专区| 97人妻精品一区二区三区麻豆| 一级毛片我不卡| 小蜜桃在线观看免费完整版高清| 久久精品国产自在天天线| 精品欧美国产一区二区三| 欧美日本视频| 七月丁香在线播放| 日韩一区二区三区影片| 三级男女做爰猛烈吃奶摸视频| 国产精品三级大全| 久久久久久久国产电影| 边亲边吃奶的免费视频| 狠狠狠狠99中文字幕| 精品久久久噜噜| 亚洲乱码一区二区免费版| 日韩亚洲欧美综合| 黄色配什么色好看| 久久人人爽人人爽人人片va| 国产一级毛片在线| 午夜福利高清视频| 亚洲va在线va天堂va国产| 国语自产精品视频在线第100页| 免费看美女性在线毛片视频| 日韩一区二区视频免费看| 美女黄网站色视频| 午夜福利成人在线免费观看| 欧美又色又爽又黄视频| 禁无遮挡网站| www.av在线官网国产| 97在线视频观看| 久久人妻av系列| 97超视频在线观看视频| 国产久久久一区二区三区| 村上凉子中文字幕在线| 水蜜桃什么品种好| 久久精品91蜜桃| av国产久精品久网站免费入址| 内射极品少妇av片p| 成年女人永久免费观看视频| 午夜久久久久精精品| 黄片wwwwww| 国产精品国产三级国产专区5o | 在线观看美女被高潮喷水网站| 欧美97在线视频| 搡女人真爽免费视频火全软件| 久久久国产成人精品二区| 日日啪夜夜撸| 国产精品一区www在线观看| 麻豆久久精品国产亚洲av| 好男人视频免费观看在线| 亚洲国产日韩欧美精品在线观看| 国产成人免费观看mmmm| 嘟嘟电影网在线观看| 日日啪夜夜撸| 日韩高清综合在线| 精品国内亚洲2022精品成人| 91久久精品电影网| 我要看日韩黄色一级片| 晚上一个人看的免费电影| 欧美区成人在线视频| 男女边吃奶边做爰视频| 国产亚洲一区二区精品| h日本视频在线播放| 日韩大片免费观看网站 | 六月丁香七月| 国产大屁股一区二区在线视频| 成年女人看的毛片在线观看| 国产精品一区二区三区四区久久| 国产美女午夜福利| 色5月婷婷丁香| 国产一区二区在线观看日韩| 国产精品久久久久久av不卡| 国产大屁股一区二区在线视频| 1000部很黄的大片| av在线天堂中文字幕| 精品99又大又爽又粗少妇毛片| 久久久久久久久久久免费av| 国产私拍福利视频在线观看| 国产精品久久久久久精品电影小说 | 啦啦啦啦在线视频资源| 91午夜精品亚洲一区二区三区| 欧美色视频一区免费| ponron亚洲| 亚洲成人久久爱视频| 国产爱豆传媒在线观看| 99久久无色码亚洲精品果冻| 亚洲经典国产精华液单| 亚洲高清免费不卡视频| 国产精品.久久久| 丝袜美腿在线中文| 熟妇人妻久久中文字幕3abv| 一级毛片久久久久久久久女| 中文资源天堂在线| 免费看美女性在线毛片视频| av在线亚洲专区| 最近视频中文字幕2019在线8| kizo精华| 一级二级三级毛片免费看| 中文字幕久久专区| 国产av不卡久久| 99热精品在线国产| 国产黄a三级三级三级人| 黄色欧美视频在线观看| 国产国拍精品亚洲av在线观看| 老司机福利观看| 精华霜和精华液先用哪个| 男人狂女人下面高潮的视频| 日本色播在线视频| 一卡2卡三卡四卡精品乱码亚洲| 久久精品久久久久久噜噜老黄 | 久久精品人妻少妇| 色视频www国产| 91av网一区二区| 亚洲美女搞黄在线观看| 日韩制服骚丝袜av| 国产av一区在线观看免费| 久久久久久伊人网av| 免费黄网站久久成人精品| 日韩精品有码人妻一区| 高清毛片免费看| 最近最新中文字幕免费大全7| 久久久精品大字幕| 中文字幕熟女人妻在线| 一个人看视频在线观看www免费| 色播亚洲综合网| av.在线天堂| 1024手机看黄色片| 三级国产精品欧美在线观看| 2021少妇久久久久久久久久久| 我的老师免费观看完整版| 精品人妻熟女av久视频| 汤姆久久久久久久影院中文字幕 | 超碰av人人做人人爽久久| 国产探花极品一区二区| 一区二区三区免费毛片| 国内精品美女久久久久久| 精品无人区乱码1区二区| 2022亚洲国产成人精品| 久久精品国产99精品国产亚洲性色| 深夜a级毛片| 男人狂女人下面高潮的视频| 久久久色成人| 中文字幕av成人在线电影|