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

    基于AIS實(shí)時(shí)數(shù)據(jù)的航道交通擁堵快速判定方法

    2018-07-09 13:06:28何正偉劉力榮王開(kāi)森
    中國(guó)航海 2018年2期
    關(guān)鍵詞:交通流服務(wù)水平航道

    何正偉, 劉力榮, 楊 帆, 王開(kāi)森

    (1.武漢理工大學(xué) 航運(yùn)學(xué)院, 武漢 430063;2.國(guó)家水運(yùn)安全工程技術(shù)研究中心, 武漢 430063;3.內(nèi)河航運(yùn)技術(shù)湖北省重點(diǎn)實(shí)驗(yàn)室, 武漢 430063)

    隨著大型船舶不斷增多,航道交通擁堵問(wèn)題時(shí)有發(fā)生,不僅會(huì)延誤船舶航行和進(jìn)出港,還會(huì)增加船舶交通事故發(fā)生的可能性。如何快速、準(zhǔn)確地判定航道交通狀態(tài)成為目前亟待解決的問(wèn)題。

    已有許多學(xué)者對(duì)航道和公路交通擁堵問(wèn)題進(jìn)行研究。劉予笑[1]分析船舶定位數(shù)據(jù),建立航道擁堵分級(jí)指標(biāo)模型,提出基于模糊綜合評(píng)價(jià)的航道擁堵識(shí)別方法;劉賽龍等[2]提出航道服務(wù)水平指標(biāo)和等級(jí),將其與船舶交通流理論相結(jié)合,建立內(nèi)河航道通過(guò)能力計(jì)算模型;劉明俊等[3]分析航道通過(guò)能力的影響因素,選取修正系數(shù),建立基于船舶流的航道通過(guò)能力計(jì)算模型;祝付玲[4]借鑒國(guó)外的公路通行能力手冊(cè)(Highway Capacity Manual,HCM),建立城市道路交通擁堵評(píng)價(jià)指標(biāo)體系;KNORR等[5]基于三相交通理論,提出通信檢測(cè)和控制擁堵技術(shù),使公路交通擁堵可觀、可控;LAKAS等[6]基于Geocast協(xié)議,制造車對(duì)車通信系統(tǒng),可有效判定公路交通狀況。

    當(dāng)前,公路交通擁堵判定方法已較為成熟,擁有科學(xué)的判定指標(biāo)體系,而航道交通擁堵判定大多基于航道通過(guò)能力來(lái)確定擁堵分級(jí)原則,屬于人工判定和事后判定,缺乏實(shí)時(shí)性,海事部門(mén)難以及時(shí)處理航道交通擁堵問(wèn)題。為解決該問(wèn)題,本文借鑒國(guó)內(nèi)外道路交通擁堵判定方法,結(jié)合交通流參數(shù)和船舶領(lǐng)域模型等,實(shí)時(shí)計(jì)算船舶自動(dòng)識(shí)別系統(tǒng)(Automatic Identification System,AIS)數(shù)據(jù),建立航道交通擁堵快速判定模型。

    1 航道交通擁堵判定原則和流程

    1.1 判定原則

    將航道交通狀態(tài)分為嚴(yán)重?fù)矶?、一般擁堵和暢通?種。嚴(yán)重?fù)矶率侵负降乐写昂芏?,船舶時(shí)動(dòng)時(shí)?;蛲耆?;一般擁堵是指航道中的船舶雖然較多,但船舶尚能低速航行;暢通是指航道中的船舶不多,船舶能按正常速度航行。

    國(guó)外相關(guān)研究[4]通過(guò)HCM確立V/C(交通量與通行能力的比值)、服務(wù)水平(公路運(yùn)行服務(wù)質(zhì)量)等公路交通擁堵評(píng)價(jià)指標(biāo)。在研究航道交通擁堵時(shí),可借鑒這些指標(biāo),并參照美國(guó)公路服務(wù)水平(見(jiàn)表1)和服務(wù)水平與速度的關(guān)系(見(jiàn)表2)。

    表1 美國(guó)公路服務(wù)水平的分類

    本文借鑒E等和F等服務(wù)水平。當(dāng)服務(wù)水平為E等時(shí),公路開(kāi)始出現(xiàn)交通擁堵,車輛行駛速度比自由流速度下降約70%;當(dāng)服務(wù)水平為F等時(shí),公路嚴(yán)重?fù)矶拢珻/V>1。同理,航道服務(wù)水平和交通狀態(tài)也遵循該規(guī)律。

    1.2 判定流程

    1) Spark實(shí)時(shí)數(shù)據(jù)處理平臺(tái)實(shí)時(shí)接收AIS在線數(shù)據(jù),對(duì)數(shù)據(jù)進(jìn)行解析、預(yù)處理、匹配和傳送等,然后自動(dòng)統(tǒng)計(jì)數(shù)據(jù)個(gè)數(shù)m1。

    2) 將經(jīng)緯度為負(fù)值或大大偏離航道的數(shù)據(jù)清除掉,剩下m2個(gè)數(shù)據(jù)。

    表2 服務(wù)水平與速度的關(guān)系 %

    3) 船舶數(shù)據(jù)集中分布的區(qū)域才是有效的,反映航道的走向。選取合適的經(jīng)緯度范圍作為有效區(qū)域,有效區(qū)域內(nèi)有m3個(gè)數(shù)據(jù)。

    4) 多數(shù)航道是彎曲的,對(duì)整段航道進(jìn)行精確計(jì)算不容易,可在其中選取某個(gè)經(jīng)緯度范圍內(nèi)的小段航道作為研究區(qū)域,該區(qū)域內(nèi)有m4個(gè)數(shù)據(jù)。

    5) 對(duì)m4個(gè)數(shù)據(jù)進(jìn)行分類,將MMSI相同的數(shù)據(jù)歸為一類,共分成m5類,即有m5艘船。

    6) 用Polynomial(多項(xiàng)式)函數(shù)將研究區(qū)域內(nèi)的數(shù)據(jù)擬合成曲線1,清除偏離主航道的數(shù)據(jù),剩下m6個(gè)數(shù)據(jù),再用Polynomial將其擬合成曲線2,曲線2可視為該段航道的主航道。

    7) 在曲線2上取(x1,y1)和(x2,y2)這2點(diǎn)作為該段航道的主航道起止點(diǎn),用經(jīng)緯度坐標(biāo)計(jì)算兩點(diǎn)間的距離,求出其長(zhǎng)度s0。

    8) 從步驟“5)”中已分類的數(shù)據(jù)中找出各船船長(zhǎng),結(jié)合s0和標(biāo)準(zhǔn)船舶艘次求出航道實(shí)際密度ρ1。

    9) 將各船船長(zhǎng)和標(biāo)準(zhǔn)船舶艘次結(jié)合起來(lái),借助船舶領(lǐng)域模型求出航道阻塞密度ρ0。

    11) 確定航道自由流速度u0。

    1.3 Spark實(shí)時(shí)數(shù)據(jù)處理平臺(tái)

    由于普通數(shù)據(jù)處理軟件難以快速處理大量AIS數(shù)據(jù),故采用Spark實(shí)時(shí)數(shù)據(jù)處理平臺(tái)。Spark是一種可處理大數(shù)據(jù)集合的低延遲集群分布式計(jì)算系統(tǒng)。[7-8]

    首先運(yùn)用Spark實(shí)時(shí)接收并解析AIS在線數(shù)據(jù),對(duì)數(shù)據(jù)進(jìn)行格式轉(zhuǎn)換、數(shù)據(jù)清洗、數(shù)據(jù)集成和剔除無(wú)效數(shù)據(jù)等預(yù)處理。隨后對(duì)數(shù)據(jù)進(jìn)行匹配,通過(guò)HDFS文件存儲(chǔ)系統(tǒng)和Hbase數(shù)據(jù)庫(kù)進(jìn)行存儲(chǔ),使數(shù)據(jù)能傳送到Spark實(shí)時(shí)數(shù)據(jù)處理平臺(tái)上。

    在該平臺(tái)上,先將AIS數(shù)據(jù)轉(zhuǎn)換成Stream數(shù)據(jù),用Spark Streaming將其按6 min/段分割成若干片段,用batch批量處理這些數(shù)據(jù),再采用Spark Engine高層算法處理,轉(zhuǎn)變成一批批結(jié)果流。

    2 構(gòu)建航道交通擁堵判定模型

    2.1 模型參數(shù)確定

    1) 平均速度。船舶交通流平均速度是指單位時(shí)間內(nèi)船舶平均對(duì)地速度,計(jì)算式為

    (1)

    式(1)中:u為船舶交通流平均速度;∑vi為船舶對(duì)地速度之和;n為船舶數(shù)量。

    2) 區(qū)域密度。船舶交通流區(qū)域密度[9]是指單位矩形區(qū)域內(nèi)存在的船舶數(shù)量,計(jì)算式為

    (2)

    式(2)中:k為船舶交通流區(qū)域密度;n為船舶數(shù)量;s為區(qū)域面積。

    3) 臨界密度。航道臨界密度[2]分為轉(zhuǎn)折密度和阻塞密度。轉(zhuǎn)折密度指當(dāng)船舶可安全暢行時(shí),航道能允許的最大交通流密度;阻塞密度指船舶密集到難以移動(dòng)時(shí)的交通流密度。本文根據(jù)阻塞密度提出一種新的航道交通擁堵判定方法。

    4) 線密度。航道線密度實(shí)際上是區(qū)域密度的一種特殊形式,指航道單位長(zhǎng)度上某個(gè)時(shí)間段內(nèi)的船舶數(shù)量,計(jì)算式為

    (3)

    式(3)中:k為航道線密度;n為船舶數(shù)量;l為航道長(zhǎng)度。

    2.2 標(biāo)準(zhǔn)船舶艘次計(jì)算

    船舶的大小會(huì)影響航道擁堵程度,故不能只考慮船舶艘次。將各船船長(zhǎng)換算成標(biāo)準(zhǔn)船長(zhǎng),設(shè)置相應(yīng)系數(shù),將船舶艘次轉(zhuǎn)換成標(biāo)準(zhǔn)船舶艘次[10](見(jiàn)表3)。標(biāo)準(zhǔn)船舶艘次的計(jì)算式為

    S=Ciμi

    (4)

    式(4)中:S為標(biāo)準(zhǔn)船舶艘次;Ci為同類船舶艘次;μi為每類船舶的轉(zhuǎn)換系數(shù)。

    表3 船舶按船長(zhǎng)分類后的轉(zhuǎn)換系數(shù)

    2.3 航道交通擁堵判定模型

    根據(jù)藤井船舶領(lǐng)域模型[11],在按船長(zhǎng)分類之后,單船的船舶領(lǐng)域長(zhǎng)軸計(jì)算式為

    li=k0Li

    (5)

    式(5)中:li為單船的船舶領(lǐng)域長(zhǎng)軸;Li為單船的船長(zhǎng);i的取值范圍為1~n2;n2為船舶數(shù)量;k0為船舶領(lǐng)域長(zhǎng)軸系數(shù)。

    若要判斷航道是否嚴(yán)重?fù)矶?,須用到臨界密度中的阻塞密度。阻塞密度等于臨界船舶數(shù)量除以航道長(zhǎng)度,是一種線密度。臨界船舶數(shù)量是指船舶可安全暢行時(shí)航道內(nèi)的最大標(biāo)準(zhǔn)船舶艘次。單船對(duì)應(yīng)的航道阻塞密度分量的計(jì)算式為

    (6)

    式(6)中:ρ0i為單船對(duì)應(yīng)的航道阻塞密度分量;n0i為單船對(duì)應(yīng)的臨界船舶數(shù)量分量;s0為航道長(zhǎng)度。

    臨界船舶數(shù)量等于航道長(zhǎng)度除以船舶領(lǐng)域長(zhǎng)軸,n0i的計(jì)算式為

    (7)

    將式(7)代入式(6)可得

    (8)

    在計(jì)算各船的平均阻塞密度時(shí),要考慮標(biāo)準(zhǔn)船舶艘次的影響,計(jì)算式為

    (9)

    式(9)中:Ci為單船對(duì)應(yīng)的船舶艘次,均等于1;ρ0為航道阻塞密度;k0為船舶領(lǐng)域長(zhǎng)軸系數(shù),一般取1.5。[2]

    判斷航道是否嚴(yán)重?fù)矶拢瑢?shí)際上是用航道實(shí)際密度與阻塞密度相比較。航道實(shí)際密度等于標(biāo)準(zhǔn)船舶艘次之和除以航道長(zhǎng)度,計(jì)算式為

    (10)

    式(10)中:ρ1為航道實(shí)際密度;n1為標(biāo)準(zhǔn)船舶艘次之和。

    借鑒公路交通擁堵判定方法,由表1可知,當(dāng)V/C>1時(shí),交通擁堵很嚴(yán)重。同理,若航道水域的總面積大于航道中所有船舶的領(lǐng)域之和,則航道處于嚴(yán)重?fù)矶聽(tīng)顟B(tài),結(jié)合式(3),可對(duì)ρ1和ρ0定義一個(gè)標(biāo)準(zhǔn):若ρ1≥ρ0,則航道處于嚴(yán)重?fù)矶聽(tīng)顟B(tài)。

    (11)

    vi對(duì)應(yīng)q個(gè)AIS點(diǎn),計(jì)算式為

    (12)

    式(12)中:vij為vi在q個(gè)點(diǎn)中對(duì)應(yīng)的第j個(gè)點(diǎn);j的取值范圍為1~q。

    自由流速度u0是指暢通狀態(tài)下的航道最高限速值,該值可定為8 kn。[12]

    綜上所述,判定航道交通狀態(tài)的標(biāo)準(zhǔn)為

    (13)

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

    3.1 長(zhǎng)江南通段航道數(shù)據(jù)試驗(yàn)

    該試驗(yàn)的數(shù)據(jù)來(lái)源于長(zhǎng)江南通段,該航段見(jiàn)圖2。試驗(yàn)于2017年2月6日10:00—10:05進(jìn)行,初始數(shù)據(jù)個(gè)數(shù)m1=2 822。

    選取東經(jīng)120.5°~121.25°、北緯31.5°~32.5°的范圍作為有效區(qū)域(見(jiàn)圖3)。

    選取東經(jīng)120.85°~120.9°、北緯31.9°~31.95°范圍內(nèi)的小段航道作為研究區(qū)域,其數(shù)據(jù)個(gè)數(shù)m4=240,分類后有24艘船。將數(shù)據(jù)擬合成曲線1(見(jiàn)圖4)。

    在圖4中,清除偏離主航道的數(shù)據(jù),將剩余數(shù)據(jù)擬合成曲線2(見(jiàn)圖5)。

    取曲線2上的(120.860,31.952)和(120.882,31.903)2點(diǎn)作為該段航道的主航道的起止點(diǎn),用式(14)求出其長(zhǎng)度s0=5.829 8 km。

    (14)

    式(14)中:R為地球半徑,約6 370 km;(x1,y1)和(x2,y2)為地球上任意2點(diǎn)。

    試驗(yàn)結(jié)果見(jiàn)表4第1行。將10:00—10:59分成10段,重復(fù)進(jìn)行試驗(yàn),結(jié)果見(jiàn)表4。

    表4 各時(shí)間間隔內(nèi)的試驗(yàn)結(jié)果

    將表4的擁堵?tīng)顟B(tài)判定結(jié)果與實(shí)際擁堵?tīng)顟B(tài)用折線圖表示,結(jié)果見(jiàn)圖6。

    3.2 湘江湘陰段航道數(shù)據(jù)試驗(yàn)

    該試驗(yàn)的數(shù)據(jù)來(lái)源于湘江湘陰段,該航段見(jiàn)圖7。試驗(yàn)于2017年2月18日17:00—17:05進(jìn)行,初始數(shù)據(jù)個(gè)數(shù)m1=301。

    選取東經(jīng)112.5°~113°、北緯28°~29°的范圍作為有效區(qū)域(見(jiàn)圖8)。

    選取東經(jīng)112.805°~112.815°、北緯28.57°~28.59°范圍內(nèi)的航道作為研究區(qū)域,其數(shù)據(jù)個(gè)數(shù)m4=77,分類后有10艘船。將數(shù)據(jù)擬合成曲線1(見(jiàn)圖9)。

    由于沒(méi)有偏離主航道的數(shù)據(jù),可取曲線1上的(112.805 2,28.571 0)和(112.811 8,28.590 0)作為該段航道的主航道的起止點(diǎn),用式(14)求出其長(zhǎng)度s0=2.208 5 km。

    試驗(yàn)結(jié)果見(jiàn)表5第1行。將17:00—17:47分成8段,重復(fù)進(jìn)行試驗(yàn),結(jié)果見(jiàn)表5。

    將表5的擁堵?tīng)顟B(tài)判定結(jié)果與實(shí)際擁堵?tīng)顟B(tài)用折線圖表示,結(jié)果見(jiàn)圖10。

    表5 各時(shí)間間隔內(nèi)的試驗(yàn)結(jié)果

    4 結(jié)束語(yǔ)

    本文借鑒國(guó)內(nèi)外道路交通擁堵判定方法,將航道交通狀態(tài)劃分為嚴(yán)重?fù)矶?、一般擁堵和暢通?種,運(yùn)用交通流參數(shù)、標(biāo)準(zhǔn)船舶艘次和船舶領(lǐng)域模型,根據(jù)參數(shù)計(jì)算公式,構(gòu)建航道交通擁堵快速判定模型。在Spark實(shí)時(shí)數(shù)據(jù)處理平臺(tái)上實(shí)時(shí)處理AIS數(shù)據(jù),運(yùn)用該模型分別對(duì)長(zhǎng)江南通段航道和湘江湘陰段航道進(jìn)行多組在線試驗(yàn)。

    試驗(yàn)結(jié)果表明:2段航道數(shù)據(jù)試驗(yàn)的判定準(zhǔn)確率均達(dá)到80%以上,準(zhǔn)確率較高;判定耗時(shí)均為1.50 s左右,判定速度較快;多次重復(fù)試驗(yàn)可使判定準(zhǔn)確率達(dá)到90%以上。

    本文提出的航道交通擁堵判定方法是一種全新的方法,考慮的參量與已有方法的參量有很大不同,思路也不一致,因此難以直接對(duì)比。本文的創(chuàng)新點(diǎn)和優(yōu)勢(shì)可概括為:

    1) 實(shí)時(shí)性高,利用詳細(xì)的AIS實(shí)時(shí)數(shù)據(jù),實(shí)時(shí)運(yùn)用判定模型,快速判定航道交通狀態(tài)。

    2) 利用計(jì)算機(jī)平臺(tái)進(jìn)行自動(dòng)化判定,大大降低了人為因素的影響。

    3) 判定準(zhǔn)確率高,判定速度快,具有可重復(fù)性、可驗(yàn)證性,明確了航道交通狀態(tài)的定量劃分,克服了以往人工判定和事后判定的缺陷,可為海事部門(mén)實(shí)時(shí)掌握航道交通狀態(tài)提供有效的技術(shù)支持。

    [1] 劉予笑. 基于模糊綜合評(píng)價(jià)的航道擁堵?tīng)顟B(tài)識(shí)別研究[D]. 杭州:浙江工商大學(xué), 2015.

    [2] 劉賽龍, 蔣璘暉. 內(nèi)河航道服務(wù)水平及通過(guò)能力研究[J]. 水運(yùn)工程, 2014(3):134-139.

    [3] 劉明俊, 萬(wàn)長(zhǎng)征. 航道通過(guò)能力影響因素的分析[J]. 船海工程, 2008,37(5): 116-118.

    [4] 祝付玲. 城市道路交通擁堵評(píng)價(jià)指標(biāo)體系研究[D]. 南京:東南大學(xué), 2006.

    [5] KNORR F, SCHRECKENBERG M. Influence of Inter-Vehicle Communication on Peak Hour Traffic Flow[J]. Physica A: Statistical Mechanics & Its Applications, 2012, 391(6): 2225-2231.

    [6] LAKAS A, CHAQFEH M. A Novel Method for Reducing Road Traffic Congestion Using Vehicular Communication[C]//International Wireless Communications & Mobile Computing Conference, 2010: 16-20.

    [7] 李文棟. 基于Spark的大數(shù)據(jù)挖掘技術(shù)的研究與實(shí)現(xiàn)[D]. 濟(jì)南:山東大學(xué), 2015.

    [8] 李天喜. 基于Spark Streaming的試驗(yàn)數(shù)據(jù)處理系統(tǒng)的研究與實(shí)現(xiàn)[D]. 西安:西安電子科技大學(xué), 2015.

    [9] 馮宏祥, 孔凡邨, 肖英杰,等. 基于AIS的元胞自動(dòng)機(jī)模型的船舶交通流特征參數(shù)分析[J]. 武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版), 2014(2):324-328.

    [10] 鄧順江, 劉明俊. 船舶領(lǐng)域模型的對(duì)比研究[J]. 中國(guó)水運(yùn), 2009, 9(6):7-9.

    [11] 丁法. 基于AIS數(shù)據(jù)的開(kāi)闊水域船舶領(lǐng)域分析研究[D]. 大連:大連海事大學(xué), 2016.

    [12] 徐春, 鄭中義, 吳兆麟. 受限航道中船舶的限速[J].大連海事大學(xué)學(xué)報(bào), 2002, 28(2):6-9.

    猜你喜歡
    交通流服務(wù)水平航道
    遂寧市:提升社保服務(wù)水平 夯實(shí)保障民生基礎(chǔ)
    加強(qiáng)圖書(shū)館管理 提高服務(wù)水平
    活力(2019年19期)2020-01-06 07:35:32
    提升糧食流通社會(huì)化服務(wù)水平的舉措構(gòu)思
    新航道
    留學(xué)(2017年5期)2017-03-29 03:03:54
    Vipersat升級(jí)版
    ——HeightsTM用高效率和智能化提升服務(wù)水平
    交通流隨機(jī)行為的研究進(jìn)展
    路內(nèi)停車對(duì)交通流延誤影響的定量分析
    具有負(fù)壓力的Aw-Rascle交通流的Riemann問(wèn)題
    我國(guó)首條40萬(wàn)噸級(jí)航道正式開(kāi)建
    水道港口(2014年1期)2014-04-27 14:14:38
    我國(guó)首條40萬(wàn)噸超大型航道將迎來(lái)開(kāi)工建設(shè)
    水道港口(2013年3期)2013-04-27 14:14:28
    av一本久久久久| 久久精品国产自在天天线| 亚洲综合色网址| 国精品久久久久久国模美| 久久久久精品性色| 精品一区二区三卡| 肉色欧美久久久久久久蜜桃| 人人妻人人澡人人看| 两个人看的免费小视频| 精品人妻在线不人妻| 90打野战视频偷拍视频| 一区在线观看完整版| 夫妻午夜视频| 免费日韩欧美在线观看| 久久久a久久爽久久v久久| 青春草视频在线免费观看| 免费大片黄手机在线观看| 一级,二级,三级黄色视频| 日韩三级伦理在线观看| 热re99久久国产66热| 美女视频免费永久观看网站| 宅男免费午夜| 母亲3免费完整高清在线观看 | 国产精品女同一区二区软件| 777久久人妻少妇嫩草av网站| 自线自在国产av| 欧美激情 高清一区二区三区| 香蕉丝袜av| 亚洲美女视频黄频| 久久精品久久精品一区二区三区| 天堂俺去俺来也www色官网| 成年女人在线观看亚洲视频| 母亲3免费完整高清在线观看 | 国产精品蜜桃在线观看| 亚洲中文av在线| 波多野结衣一区麻豆| 成人漫画全彩无遮挡| 久久久久网色| 亚洲伊人色综图| 国产欧美亚洲国产| 一区二区av电影网| 91精品伊人久久大香线蕉| 久热久热在线精品观看| 国产成人精品福利久久| 精品国产一区二区久久| 一区二区三区乱码不卡18| 中文字幕亚洲精品专区| 午夜福利在线免费观看网站| 视频区图区小说| 免费观看性生交大片5| 一区二区三区激情视频| 国产精品二区激情视频| 国产黄频视频在线观看| 精品第一国产精品| 青春草国产在线视频| av视频免费观看在线观看| 天天操日日干夜夜撸| 自拍欧美九色日韩亚洲蝌蚪91| 黄色怎么调成土黄色| 久久精品久久久久久久性| 日韩一区二区视频免费看| 叶爱在线成人免费视频播放| 国产男人的电影天堂91| 久久久久精品性色| 国产伦理片在线播放av一区| av在线观看视频网站免费| 一边亲一边摸免费视频| 777久久人妻少妇嫩草av网站| 国产一区有黄有色的免费视频| 波野结衣二区三区在线| 日本wwww免费看| 午夜福利视频精品| tube8黄色片| 免费观看av网站的网址| 777米奇影视久久| av不卡在线播放| 亚洲欧美色中文字幕在线| 免费观看av网站的网址| 男女高潮啪啪啪动态图| 啦啦啦啦在线视频资源| 国产在线免费精品| 日韩伦理黄色片| 久久国产亚洲av麻豆专区| 欧美成人午夜免费资源| 亚洲欧美精品自产自拍| 色哟哟·www| 欧美国产精品va在线观看不卡| 香蕉精品网在线| 黄色配什么色好看| 中文欧美无线码| 91久久精品国产一区二区三区| 人妻少妇偷人精品九色| 亚洲精品一二三| 中文欧美无线码| 天堂8中文在线网| 久热久热在线精品观看| av.在线天堂| 国产97色在线日韩免费| 一二三四在线观看免费中文在| 美女中出高潮动态图| 美女大奶头黄色视频| 国产亚洲最大av| 天堂俺去俺来也www色官网| 十八禁网站网址无遮挡| 色婷婷av一区二区三区视频| 水蜜桃什么品种好| 看免费成人av毛片| 亚洲伊人久久精品综合| 国产精品免费视频内射| 亚洲成av片中文字幕在线观看 | 亚洲精品日本国产第一区| 日韩欧美一区视频在线观看| 国产熟女午夜一区二区三区| av国产久精品久网站免费入址| 中文字幕制服av| 亚洲中文av在线| 国产97色在线日韩免费| 少妇熟女欧美另类| 性高湖久久久久久久久免费观看| 桃花免费在线播放| 亚洲精品成人av观看孕妇| 国产精品二区激情视频| 中文字幕色久视频| 少妇 在线观看| 高清视频免费观看一区二区| 亚洲av成人精品一二三区| 亚洲欧美一区二区三区久久| 成年女人在线观看亚洲视频| av福利片在线| 春色校园在线视频观看| 久久婷婷青草| av在线观看视频网站免费| 黑人欧美特级aaaaaa片| 777久久人妻少妇嫩草av网站| 性色av一级| 色视频在线一区二区三区| 国产精品久久久久久精品电影小说| 国产无遮挡羞羞视频在线观看| 国产精品香港三级国产av潘金莲 | 热re99久久国产66热| 婷婷色av中文字幕| 亚洲欧美成人精品一区二区| 美女高潮到喷水免费观看| 两性夫妻黄色片| www.精华液| 超碰成人久久| 欧美成人午夜免费资源| 成人亚洲精品一区在线观看| 国产欧美日韩综合在线一区二区| 多毛熟女@视频| 人妻少妇偷人精品九色| 伊人久久国产一区二区| 中文字幕人妻丝袜一区二区 | 国产成人精品久久二区二区91 | 亚洲精品久久成人aⅴ小说| 人人妻人人澡人人爽人人夜夜| 2022亚洲国产成人精品| 国产亚洲午夜精品一区二区久久| 国产精品女同一区二区软件| 国产精品久久久久成人av| 久热久热在线精品观看| 日韩视频在线欧美| 免费大片黄手机在线观看| 亚洲av国产av综合av卡| 国产精品人妻久久久影院| 人人妻人人添人人爽欧美一区卜| 中文字幕精品免费在线观看视频| 激情视频va一区二区三区| 精品一品国产午夜福利视频| 人体艺术视频欧美日本| 成人手机av| 男女高潮啪啪啪动态图| 少妇 在线观看| 大香蕉久久网| 两性夫妻黄色片| 中文字幕最新亚洲高清| 亚洲,一卡二卡三卡| 九色亚洲精品在线播放| 纯流量卡能插随身wifi吗| 国产成人午夜福利电影在线观看| 亚洲国产日韩一区二区| 国产精品偷伦视频观看了| 色哟哟·www| 人妻人人澡人人爽人人| 午夜精品国产一区二区电影| 亚洲精品一二三| 伊人亚洲综合成人网| 在线观看三级黄色| av国产久精品久网站免费入址| 国产精品国产三级国产专区5o| 久久久久久久亚洲中文字幕| 国产欧美亚洲国产| av片东京热男人的天堂| 亚洲av综合色区一区| 欧美激情高清一区二区三区 | 国语对白做爰xxxⅹ性视频网站| 日韩视频在线欧美| 欧美bdsm另类| 久久久久网色| 欧美人与善性xxx| 亚洲伊人久久精品综合| 少妇 在线观看| 久久女婷五月综合色啪小说| 美女xxoo啪啪120秒动态图| 熟妇人妻不卡中文字幕| 一区二区三区精品91| 99九九在线精品视频| 18+在线观看网站| 国产麻豆69| 少妇人妻精品综合一区二区| 叶爱在线成人免费视频播放| 一本色道久久久久久精品综合| 日韩三级伦理在线观看| 久久av网站| 久久久久久久久久久久大奶| 爱豆传媒免费全集在线观看| 国产精品女同一区二区软件| 水蜜桃什么品种好| 久久精品aⅴ一区二区三区四区 | 天堂中文最新版在线下载| 性少妇av在线| 黄色毛片三级朝国网站| 久久久精品区二区三区| 久热这里只有精品99| 亚洲av国产av综合av卡| 久久久久久久精品精品| 18禁观看日本| 91午夜精品亚洲一区二区三区| 久久久久久久久免费视频了| av视频免费观看在线观看| 国产成人91sexporn| 一二三四在线观看免费中文在| 在线天堂中文资源库| 亚洲内射少妇av| 成年av动漫网址| 午夜福利网站1000一区二区三区| av卡一久久| 精品视频人人做人人爽| 亚洲av男天堂| 免费黄色在线免费观看| 国产成人精品无人区| 热re99久久国产66热| 宅男免费午夜| 中文精品一卡2卡3卡4更新| 人妻系列 视频| 成人亚洲精品一区在线观看| 久久久久久久精品精品| 久久这里有精品视频免费| 日韩中字成人| 亚洲欧美中文字幕日韩二区| 成人亚洲精品一区在线观看| 在线观看美女被高潮喷水网站| 国产在线一区二区三区精| 男女下面插进去视频免费观看| 亚洲国产精品一区三区| 美女中出高潮动态图| 亚洲国产av新网站| 老汉色∧v一级毛片| 一级片免费观看大全| 性色avwww在线观看| 捣出白浆h1v1| 91久久精品国产一区二区三区| 中文字幕亚洲精品专区| 两性夫妻黄色片| 欧美日韩精品网址| 精品国产露脸久久av麻豆| 校园人妻丝袜中文字幕| 久久精品熟女亚洲av麻豆精品| 国产精品无大码| 大片免费播放器 马上看| 宅男免费午夜| 国产精品免费大片| 一级爰片在线观看| 久久久久精品久久久久真实原创| 亚洲成人手机| 国产精品久久久久久精品古装| 精品少妇久久久久久888优播| 亚洲婷婷狠狠爱综合网| 韩国精品一区二区三区| 国产有黄有色有爽视频| 免费不卡的大黄色大毛片视频在线观看| 婷婷色综合www| 亚洲欧美日韩另类电影网站| 涩涩av久久男人的天堂| 在线看a的网站| 欧美少妇被猛烈插入视频| 久久久久久免费高清国产稀缺| 午夜福利视频在线观看免费| 99国产综合亚洲精品| 亚洲av男天堂| 午夜福利视频精品| 精品人妻一区二区三区麻豆| 免费高清在线观看日韩| 大码成人一级视频| 亚洲精品成人av观看孕妇| 日本黄色日本黄色录像| 色婷婷av一区二区三区视频| 美女国产高潮福利片在线看| 欧美黄色片欧美黄色片| 精品一区二区三区四区五区乱码 | 日日啪夜夜爽| 肉色欧美久久久久久久蜜桃| 一级毛片电影观看| 国产成人欧美| 色94色欧美一区二区| 国产成人精品婷婷| 少妇人妻精品综合一区二区| 国产精品二区激情视频| 亚洲精品,欧美精品| 美女xxoo啪啪120秒动态图| 丰满饥渴人妻一区二区三| 老汉色av国产亚洲站长工具| 亚洲国产av影院在线观看| 在线 av 中文字幕| 深夜精品福利| 日韩伦理黄色片| 国产精品熟女久久久久浪| 1024香蕉在线观看| 少妇人妻 视频| 欧美国产精品va在线观看不卡| 成年动漫av网址| 观看av在线不卡| 欧美国产精品va在线观看不卡| 亚洲人成网站在线观看播放| 街头女战士在线观看网站| 日韩欧美精品免费久久| 天天躁狠狠躁夜夜躁狠狠躁| av在线观看视频网站免费| 少妇猛男粗大的猛烈进出视频| 国产xxxxx性猛交| 成人国语在线视频| 精品亚洲乱码少妇综合久久| 69精品国产乱码久久久| 如何舔出高潮| 亚洲精品久久久久久婷婷小说| 色婷婷久久久亚洲欧美| 国产福利在线免费观看视频| av女优亚洲男人天堂| 久久精品国产亚洲av高清一级| 亚洲国产欧美网| 岛国毛片在线播放| av网站在线播放免费| 999久久久国产精品视频| 亚洲中文av在线| 日韩三级伦理在线观看| 欧美日韩精品网址| 高清视频免费观看一区二区| 亚洲一级一片aⅴ在线观看| 2021少妇久久久久久久久久久| videosex国产| 最黄视频免费看| 欧美国产精品va在线观看不卡| 欧美老熟妇乱子伦牲交| 一级毛片 在线播放| 成人二区视频| 99国产综合亚洲精品| 美女xxoo啪啪120秒动态图| 久久精品国产a三级三级三级| 午夜福利网站1000一区二区三区| 久久久久精品久久久久真实原创| 日本av手机在线免费观看| 中文字幕制服av| 免费观看在线日韩| √禁漫天堂资源中文www| 欧美精品亚洲一区二区| 大香蕉久久网| 亚洲成国产人片在线观看| 少妇 在线观看| 免费高清在线观看日韩| 亚洲欧洲精品一区二区精品久久久 | 午夜福利乱码中文字幕| 少妇被粗大猛烈的视频| 我要看黄色一级片免费的| 国产伦理片在线播放av一区| 18禁裸乳无遮挡动漫免费视频| 在线免费观看不下载黄p国产| 国产成人精品福利久久| 欧美中文综合在线视频| 久久精品国产亚洲av天美| 中文天堂在线官网| 亚洲国产精品一区三区| 一区二区三区精品91| 国产免费视频播放在线视频| 久久国产精品男人的天堂亚洲| 欧美亚洲 丝袜 人妻 在线| av国产久精品久网站免费入址| 亚洲人成网站在线观看播放| 亚洲美女视频黄频| 如何舔出高潮| 岛国毛片在线播放| 亚洲精品国产一区二区精华液| 麻豆精品久久久久久蜜桃| 亚洲在久久综合| 成人国产av品久久久| 91aial.com中文字幕在线观看| 国语对白做爰xxxⅹ性视频网站| 国产精品国产av在线观看| 七月丁香在线播放| 黄色 视频免费看| 亚洲一区二区三区欧美精品| 亚洲精品国产av成人精品| 国产成人一区二区在线| 国产成人免费观看mmmm| 黄色怎么调成土黄色| 伊人久久大香线蕉亚洲五| 汤姆久久久久久久影院中文字幕| 日韩 亚洲 欧美在线| 精品福利永久在线观看| 一边摸一边做爽爽视频免费| 69精品国产乱码久久久| 97在线人人人人妻| 精品午夜福利在线看| 天天影视国产精品| 男女边吃奶边做爰视频| 午夜免费男女啪啪视频观看| 成人漫画全彩无遮挡| 国产伦理片在线播放av一区| 黄色毛片三级朝国网站| 一级毛片黄色毛片免费观看视频| 99久国产av精品国产电影| 午夜福利在线观看免费完整高清在| 欧美人与性动交α欧美软件| 欧美精品一区二区大全| 在线免费观看不下载黄p国产| 亚洲精品视频女| 另类精品久久| 少妇人妻 视频| 精品亚洲成国产av| 有码 亚洲区| 久久久欧美国产精品| 在线观看国产h片| 国产一区二区三区综合在线观看| 少妇猛男粗大的猛烈进出视频| 国产免费现黄频在线看| 免费日韩欧美在线观看| 国产精品国产三级专区第一集| 建设人人有责人人尽责人人享有的| 国产黄频视频在线观看| 丝瓜视频免费看黄片| 亚洲成国产人片在线观看| 日韩精品免费视频一区二区三区| 国产精品 欧美亚洲| 人人妻人人爽人人添夜夜欢视频| 国产av国产精品国产| 久久国产精品大桥未久av| 亚洲四区av| 国产成人精品久久二区二区91 | 麻豆精品久久久久久蜜桃| 91国产中文字幕| 男人舔女人的私密视频| 777米奇影视久久| 18+在线观看网站| 亚洲 欧美一区二区三区| 叶爱在线成人免费视频播放| 最近2019中文字幕mv第一页| 成人手机av| 韩国精品一区二区三区| 大片电影免费在线观看免费| 深夜精品福利| 性色avwww在线观看| 美女大奶头黄色视频| 免费看av在线观看网站| 人人妻人人澡人人看| 一级爰片在线观看| 精品第一国产精品| 老司机影院毛片| av不卡在线播放| 2022亚洲国产成人精品| 99国产综合亚洲精品| 侵犯人妻中文字幕一二三四区| 亚洲经典国产精华液单| 91精品国产国语对白视频| 丝袜美腿诱惑在线| 国产精品免费视频内射| 欧美日韩一级在线毛片| 欧美精品高潮呻吟av久久| 免费av中文字幕在线| 亚洲在久久综合| 少妇猛男粗大的猛烈进出视频| 宅男免费午夜| 国产精品三级大全| 亚洲精品日本国产第一区| 欧美成人精品欧美一级黄| 午夜福利网站1000一区二区三区| 亚洲国产最新在线播放| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av网站在线播放免费| 久久久国产欧美日韩av| 亚洲精品乱久久久久久| 一级爰片在线观看| √禁漫天堂资源中文www| 午夜免费观看性视频| 国产成人精品婷婷| 国产淫语在线视频| 男女高潮啪啪啪动态图| 中文字幕人妻丝袜一区二区 | av国产精品久久久久影院| 国产精品国产av在线观看| 国产男女内射视频| 国产又色又爽无遮挡免| 精品第一国产精品| 久久av网站| 精品一区在线观看国产| 久久久欧美国产精品| 国产男人的电影天堂91| av线在线观看网站| 一区二区av电影网| 99久久人妻综合| 如日韩欧美国产精品一区二区三区| 人妻人人澡人人爽人人| 黄色配什么色好看| 精品一区在线观看国产| 日韩欧美一区视频在线观看| 777米奇影视久久| 亚洲三级黄色毛片| 久久久久视频综合| 国产精品久久久av美女十八| 亚洲av福利一区| 美女中出高潮动态图| 色哟哟·www| 欧美97在线视频| 免费久久久久久久精品成人欧美视频| 老熟女久久久| 99热网站在线观看| 久久精品国产a三级三级三级| 26uuu在线亚洲综合色| 国产精品久久久久久久久免| 一区二区三区精品91| 婷婷成人精品国产| 精品国产乱码久久久久久小说| 热re99久久精品国产66热6| 国产 一区精品| 热re99久久精品国产66热6| 日本色播在线视频| 国产精品久久久久久精品古装| 91国产中文字幕| 久热久热在线精品观看| 免费在线观看黄色视频的| 极品少妇高潮喷水抽搐| 丝袜脚勾引网站| 中文天堂在线官网| 97精品久久久久久久久久精品| 高清黄色对白视频在线免费看| 黄频高清免费视频| 国产激情久久老熟女| 国产成人午夜福利电影在线观看| 亚洲av福利一区| 亚洲人成网站在线观看播放| 亚洲国产欧美在线一区| av在线老鸭窝| 自拍欧美九色日韩亚洲蝌蚪91| 久久国产精品男人的天堂亚洲| 久久久久视频综合| 日韩欧美一区视频在线观看| 搡老乐熟女国产| 精品少妇黑人巨大在线播放| 黄色毛片三级朝国网站| 精品国产一区二区三区久久久樱花| 亚洲婷婷狠狠爱综合网| 国产av精品麻豆| 在现免费观看毛片| 人成视频在线观看免费观看| 亚洲精品成人av观看孕妇| 久久99精品国语久久久| 久久久久久免费高清国产稀缺| 国产精品偷伦视频观看了| 青春草国产在线视频| 大香蕉久久成人网| 黄网站色视频无遮挡免费观看| 国产黄频视频在线观看| 久久久久久久国产电影| 国产精品秋霞免费鲁丝片| 日本91视频免费播放| www.精华液| 青春草视频在线免费观看| 日韩三级伦理在线观看| 国产日韩欧美在线精品| 各种免费的搞黄视频| 日韩欧美一区视频在线观看| 91久久精品国产一区二区三区| av视频免费观看在线观看| 久久ye,这里只有精品| 欧美老熟妇乱子伦牲交| 最近中文字幕高清免费大全6| 欧美日本中文国产一区发布| 亚洲,一卡二卡三卡| 午夜福利在线免费观看网站| 热99国产精品久久久久久7| 成人毛片a级毛片在线播放| 日韩 亚洲 欧美在线| 亚洲精品久久久久久婷婷小说| 激情视频va一区二区三区| 久久久久精品久久久久真实原创| 老熟女久久久| av网站在线播放免费| 久久亚洲国产成人精品v| 精品少妇久久久久久888优播| 黄色配什么色好看| 精品卡一卡二卡四卡免费| 久久人妻熟女aⅴ| 亚洲五月色婷婷综合| 最黄视频免费看| 欧美精品高潮呻吟av久久| 久久午夜综合久久蜜桃| 成人免费观看视频高清| 久久国内精品自在自线图片| 中文天堂在线官网| 丝袜在线中文字幕| 日日撸夜夜添| 日本欧美视频一区| 青春草亚洲视频在线观看| 岛国毛片在线播放| 亚洲国产欧美网| 成人毛片a级毛片在线播放| 久久久精品区二区三区|