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

    基于駕駛實(shí)踐的無人船智能避碰決策方法

    2021-03-08 02:49:28丁志國張新宇王程博黎泉安蘭軒
    中國艦船研究 2021年1期
    關(guān)鍵詞:船舶智能模型

    丁志國,張新宇*,王程博,黎泉,安蘭軒

    1 大連海事大學(xué) 航海動(dòng)態(tài)仿真與控制交通行業(yè)重點(diǎn)實(shí)驗(yàn)室,遼寧 大連 116026

    2 大連海事大學(xué) 航海學(xué)院,遼寧 大連 116026

    0 引 言

    隨著新一代人工智能技術(shù)和無人駕駛技術(shù)在航海領(lǐng)域的興起,海上運(yùn)輸正朝著智能化方向快速發(fā)展。無人駕駛船舶作為海上智能交通的載體,已成為國內(nèi)外航運(yùn)界關(guān)注的焦點(diǎn)。近年來,無人駕駛車輛的智能避碰決策已取得較好成果,國內(nèi)外許多專家學(xué)者積極借鑒無人駕駛車輛的經(jīng)驗(yàn),結(jié)合航海背景,試圖解決船舶的智能避碰決策問題。

    無人駕駛船舶自主航行包括感知、理解、決策、控制4 個(gè)階段,而決策階段是無人駕駛船舶安全航行及避障的關(guān)鍵,因此無人駕駛船舶的智能避碰決策成為國內(nèi)外學(xué)者的研究重點(diǎn)。無人駕駛船舶智能避碰決策問題的核心在于能否在海上環(huán)境復(fù)雜、多變的情況下快速做出符合《國際海上避碰規(guī)則》(以下簡稱《規(guī)則》)以及海員通常做法的安全而有效的決策。實(shí)際上,船舶智能避碰決策要解決的實(shí)際問題與船舶駕駛?cè)藛T在避碰實(shí)際中要解決的問題是一致的,區(qū)別在于是用人工方式解決還是由計(jì)算機(jī)自主做出決策,而目前基于駕駛實(shí)踐視角的無人船智能避碰決策研究較少,本文擬從駕駛實(shí)踐的視角開展無人船的智能避碰決策研究。

    基于駕駛實(shí)踐,無人駕駛船舶智能避碰決策包含航行態(tài)勢(shì)感知、航行態(tài)勢(shì)理解、碰撞危險(xiǎn)度決策、避碰行動(dòng)時(shí)機(jī)決策、避碰行動(dòng)有效性決策和復(fù)航?jīng)Q策等過程。

    航行態(tài)勢(shì)感知系統(tǒng)是借助多種船載傳感器設(shè)備主動(dòng)對(duì)船舶內(nèi)、外部信息以及航行環(huán)境信息進(jìn)行多模態(tài)感知。

    航行態(tài)勢(shì)理解是通過對(duì)感知系統(tǒng)獲取的多源異構(gòu)、多模態(tài)信息進(jìn)行分類與理解,將無人駕駛船舶的航行態(tài)勢(shì)進(jìn)行結(jié)構(gòu)化表述,并結(jié)合《規(guī)則》與海員的通常做法對(duì)船舶會(huì)遇場(chǎng)景進(jìn)行量化劃分。目前,對(duì)于會(huì)遇場(chǎng)景的劃分標(biāo)準(zhǔn)不一。陳國權(quán)[1]根據(jù)相對(duì)運(yùn)動(dòng)幾何規(guī)律分析總結(jié)了目標(biāo)交會(huì)特征(TEC),通過目標(biāo)船相對(duì)運(yùn)動(dòng)航向、船速比的大小以及目標(biāo)船與本船的航向差取值范圍歸納出了15 種TEC,并針對(duì)不同的TEC 給出了對(duì)應(yīng)的避碰方案;Perera等[2]將本船碰撞區(qū)域劃分為8 個(gè)區(qū),將目標(biāo)船的相對(duì)運(yùn)動(dòng)航向劃分為8 個(gè)區(qū),然后按照船舶相對(duì)距離范圍、速比條件將會(huì)遇情況劃分為了144 種,其中有效情況為63 種。會(huì)遇場(chǎng)景量化程度越高,對(duì)避碰決策越有利,但隨著引入的影響因子的增加,系統(tǒng)復(fù)雜性越高,遍歷需要的時(shí)間也越長。本文將借鑒無人駕駛車輛場(chǎng)景建模方法,通過引入本體論模型,建立航行態(tài)勢(shì)本體模型,并根據(jù)《規(guī)則》及駕駛實(shí)踐對(duì)會(huì)遇場(chǎng)景進(jìn)行知識(shí)表述,實(shí)現(xiàn)在線推理,給出會(huì)遇場(chǎng)景、本船的避讓責(zé)任以及避讓行動(dòng)方式,為智能避碰決策提供決策基礎(chǔ),提高決策效率。

    碰撞危險(xiǎn)度是采取避碰行動(dòng)的依據(jù)。Kearon等[3]首次提出由最近會(huì)遇距離DCPA和到達(dá)最近會(huì)遇點(diǎn)的時(shí)間TCPA加權(quán)的方法確定碰撞危險(xiǎn)度,其只考慮到了2 個(gè)因素對(duì)碰撞危險(xiǎn)度的影響,可能會(huì)對(duì)實(shí)際情況產(chǎn)生誤判;文獻(xiàn)[4-6]均是采用神經(jīng)網(wǎng)絡(luò)的方法來確定船舶碰撞危險(xiǎn)度,其將DCPA和TCPA、舷角以及航向航速作為系統(tǒng)輸入,具有良好的預(yù)測(cè)效果;文獻(xiàn)[7-9]利用模糊理論研究了碰撞危險(xiǎn)度;劉冬冬等[7]在綜合考慮DCPA和TCPA、來船方位、距離和船速比的基礎(chǔ)上,構(gòu)建了碰撞危險(xiǎn)度模糊評(píng)價(jià)模型,并運(yùn)用四元船舶領(lǐng)域(QSD)對(duì)DCPA的隸屬度函數(shù)進(jìn)行了改進(jìn),可為避碰決策提供準(zhǔn)確依據(jù);胥文等[8]在考慮以上5 個(gè)影響因素的基礎(chǔ)上,結(jié)合船舶操縱性能、能見度情況以及航行水域情況對(duì)模糊隸屬度函數(shù)進(jìn)行修正,提高了碰撞危險(xiǎn)度研究的準(zhǔn)確性。本文將在前人研究的基礎(chǔ)上,改進(jìn)碰撞危險(xiǎn)度隸屬度函數(shù),并根據(jù)航行駕駛實(shí)踐,重新設(shè)計(jì)碰撞危險(xiǎn)度合成模型。

    無人駕駛船舶完成航行態(tài)勢(shì)感知與理解以及碰撞危險(xiǎn)決策之后,需要根據(jù)駕駛實(shí)踐做出合理的避碰決策,一些學(xué)者遷移應(yīng)用機(jī)器人避碰路徑規(guī)劃研究中的勢(shì)場(chǎng)理論來解決船舶的避碰問題[10-12],雖然能夠在理論上實(shí)現(xiàn)障礙物的避碰,但避碰幅度以及避碰行動(dòng)方式受引力場(chǎng)及斥力場(chǎng)的影響,不能很好地契合航海實(shí)踐要求;一些學(xué)者認(rèn)為,船舶避碰時(shí)機(jī)與碰撞危險(xiǎn)度CRI直接相關(guān),并以CRI為評(píng)價(jià)指標(biāo)來求取避碰策略[13-14],但CRI受航行環(huán)境、能見度、兩船間距、兩船方位、DCPA和TCPA以及船速比等多個(gè)因素影響,且各因素之間存在相互制約關(guān)系,不同的研究者對(duì)不同的影響因素進(jìn)行分析,CRI的計(jì)算結(jié)果存在一定的差異,且均沒有完備地考慮所有影響,直接將CRI作為船舶避碰的決定性指標(biāo)會(huì)存在較大誤差;一些學(xué)者結(jié)合深度強(qiáng)化學(xué)習(xí)算法對(duì)船舶避碰行為決策展開了研究[15-16],其通過構(gòu)建適用于無人駕駛船舶的深度強(qiáng)化學(xué)習(xí)模型(DRL),通過強(qiáng)化訓(xùn)練得到了基本的避碰模型,但該方法需要大量的訓(xùn)練數(shù)據(jù),且模型結(jié)構(gòu)對(duì)數(shù)據(jù)的依賴性較強(qiáng),算法訓(xùn)練耗時(shí)長,對(duì)于未知環(huán)境下的避碰無法達(dá)到實(shí)時(shí)性的要求;李麗娜等[17]提出了擬人化智能避碰決策方法(PIDVCA),通過模仿經(jīng)驗(yàn)豐富駕駛員(避碰專家)的避碰決策,利用機(jī)器學(xué)習(xí)構(gòu)建動(dòng)態(tài)避碰規(guī)則庫,做出了符合《規(guī)則》和海員通常做法的合理決策,避碰效果較好,但存在的問題就是不能完備地考慮船舶在航行過程中的所有會(huì)遇情況,且隨著會(huì)遇情況的增加,規(guī)則庫規(guī)模逐步擴(kuò)大,遍歷需要更長的時(shí)間。

    在船舶避碰過程中如何合理、恰當(dāng)?shù)卮_定避碰時(shí)機(jī)、避碰行動(dòng)方式以及避碰幅度非常重要,《規(guī)則》第十六條規(guī)定:“需給他船讓路的船舶,應(yīng)盡可能及早地采取大幅度的行動(dòng),寬裕地讓清他船”。國內(nèi)外專家學(xué)者采用幾何避碰原理、速度障礙法、勢(shì)場(chǎng)理論、模糊數(shù)學(xué)方法、神經(jīng)網(wǎng)絡(luò)、智能優(yōu)化算法、人工智能及專家系統(tǒng)等方法建立了船舶避碰決策模型,這些方法對(duì)于避碰決策的發(fā)展有一定的理論指導(dǎo)意義,但共性問題是不能夠很好地貼近航海實(shí)際。

    目前,基于駕駛實(shí)踐視角的無人船智能避碰決策研究較少,鑒于此,本文擬從駕駛實(shí)踐的角度討論船舶避碰決策問題。首先,對(duì)無人駕駛船舶智能避碰決策的合理性和實(shí)時(shí)性以及唯一性問題進(jìn)行論述,然后建立智能避碰決策過程的航行態(tài)勢(shì)本體模型、多元復(fù)合避碰危險(xiǎn)度評(píng)估模型以及BOP 智能避碰決策模型這3 個(gè)模型,最后,設(shè)計(jì)仿真實(shí)驗(yàn),并在典型會(huì)遇場(chǎng)景下驗(yàn)證方法的有效性。

    1 問題分析

    從有人駕駛船舶駕駛實(shí)踐的角度分析無人駕駛船舶在做避碰決策時(shí)應(yīng)主要考慮的問題有:

    1) 決策的實(shí)時(shí)性。

    船舶在海上航行時(shí),受風(fēng)、浪、流、能見度、水深等因素的約束,在避碰決策時(shí)需要綜合予以考慮。在復(fù)雜、多變的環(huán)境中快速做出響應(yīng)是對(duì)避碰決策算法的要求,如果決策算法不能保證實(shí)時(shí)性或不能在可接受的時(shí)間范圍內(nèi)做出決策,將嚴(yán)重威脅船舶的航行安全。

    2) 決策的合理性。

    船舶在做自主避碰決策時(shí),需要最大限度地實(shí)現(xiàn)擬人化。有人駕駛船舶在操縱船舶的過程中,駕駛員會(huì)結(jié)合周圍的航行環(huán)境、航行規(guī)則及本船的操縱性能,將人、船環(huán)境視為整體進(jìn)行考慮,主動(dòng)做出綜合性的、有經(jīng)驗(yàn)的決策,無人駕駛船舶自主決策同樣也需要根據(jù)感知系統(tǒng)傳來的多源異構(gòu)信息進(jìn)行綜合處理分析,得出最優(yōu)的避碰策略,從而實(shí)現(xiàn)擬人化且最優(yōu)的決策。

    3) 決策的唯一性。

    船舶駕駛員在進(jìn)行避碰決策的過程中,是通過評(píng)判碰撞危險(xiǎn)程度,根據(jù)駕駛航行經(jīng)驗(yàn)給出一個(gè)準(zhǔn)確的決策信息(舵角或者航速)來控制船舶避碰,這就要求無人駕駛船舶在執(zhí)行決策算法的時(shí)候不能出現(xiàn)多個(gè)解,只能輸出唯一的、確定的、可執(zhí)行的決策指令。

    本文在進(jìn)行船舶智能避碰決策方法研究的同時(shí),充分考慮了以上問題,在此基礎(chǔ)上,構(gòu)建無人駕駛船舶航行態(tài)勢(shì)本體概念模型,并結(jié)合《規(guī)則》及海員的通常做法對(duì)船舶會(huì)遇場(chǎng)景進(jìn)行量化處理,建立船舶復(fù)合碰撞危險(xiǎn)度評(píng)估模型,從駕駛實(shí)踐的角度出發(fā),建立駕駛員視角(BOP)智能避碰決策模型,以船舶避碰總路徑最短為目標(biāo)函數(shù)求解最優(yōu)避碰策略,實(shí)現(xiàn)船舶智能航行及自主避碰。

    2 模型構(gòu)建

    無人駕駛船舶自主航行需要實(shí)時(shí)與航行環(huán)境交互,獲取航行狀態(tài)信息,通過學(xué)習(xí)算法做出行為決策,其系統(tǒng)架構(gòu)如圖1 所示。輸入船舶當(dāng)前航行狀態(tài)、周圍礙航物信息,信息傳入行為決策層,經(jīng)一系列解算后輸出船舶的決策信息并交由控制執(zhí)行層完成對(duì)船舶的運(yùn)動(dòng)控制。圖中, λ ,? 分別為船舶的經(jīng)度和緯度;Ψ為船舶航向;u為船舶前進(jìn)速度;v為船舶橫移速度;r為船舶轉(zhuǎn)艏角速度;ΔΨ為決策輸出的航向改變量;ΔV為決策輸出的航速改變量; δ為船舶舵角;V為船舶航速。

    2.1 無人駕駛船舶航行態(tài)勢(shì)本體模型

    準(zhǔn)確的無人駕駛船舶航行態(tài)勢(shì)理解是行為決策的基礎(chǔ),避碰場(chǎng)景理解主要是基于對(duì)《規(guī)則》以及海員通常做法的理解,劃分典型的會(huì)遇場(chǎng)景,為避碰決策提供先決性條件。

    圖1 無人駕駛船舶自主航行系統(tǒng)架構(gòu)圖Fig. 1 architecture diagram of autonomous navigation system of unmanned ship

    《規(guī)則》中給出了對(duì)遇、交叉相遇、追越3 種會(huì)遇場(chǎng)景,如圖2 所示。但是在實(shí)際的船舶航行過程中對(duì)于不同的會(huì)遇情況,船舶有不同的航路權(quán),故需要對(duì)會(huì)遇場(chǎng)景進(jìn)行詳細(xì)劃分。

    圖2 3 種典型會(huì)遇場(chǎng)景示意圖Fig. 2 Schematic diagram of three typical encounter scenarios

    目前,有關(guān)無人船的航行態(tài)勢(shì)理解與建模的研究較少,無統(tǒng)一化的標(biāo)準(zhǔn)定義?;诤唵巍癐F···THEN···”推理規(guī)則的避碰知識(shí)庫,難以描述各種會(huì)遇形式和避碰階段的避碰方案。本文遷移應(yīng)用無人車的場(chǎng)景建模方法[18],利用本體模型的強(qiáng)大信息表述能力,可邏輯化表述并量化航行態(tài)勢(shì),結(jié)合駕駛實(shí)踐,對(duì)先驗(yàn)航行知識(shí)進(jìn)行描述,提高避碰決策效率。

    航行態(tài)勢(shì)理解是通過對(duì)感知系統(tǒng)獲取的多源異構(gòu)信息進(jìn)行分類與理解,將無人駕駛船舶的航行態(tài)勢(shì)進(jìn)行結(jié)構(gòu)化表述。首先,根據(jù)“本體”理論將無人駕駛船舶的航行態(tài)勢(shì)分為實(shí)體和屬性。實(shí)體類代表著不同本質(zhì)的客觀存在物,屬性類用于描述對(duì)象實(shí)體以及實(shí)體之間的關(guān)系屬性。本文將屬性類又分為實(shí)體屬性和關(guān)系屬性,其中實(shí)體屬性包括實(shí)體的位置、航向航速,關(guān)系屬性包括位置關(guān)系、航向關(guān)系、方位關(guān)系和速度關(guān)系。

    感知系統(tǒng)獲取的信息包含靜態(tài)信息和動(dòng)態(tài)信息。靜態(tài)信息一般指靜態(tài)障礙物實(shí)體,包括具有結(jié)構(gòu)化的海圖實(shí)體以及其他靜態(tài)障礙物實(shí)體,其中海圖實(shí)體包括具有結(jié)構(gòu)化特征的航道、導(dǎo)助航標(biāo)志等,靜態(tài)障礙物實(shí)體是指靜止的非結(jié)構(gòu)化實(shí)體,包括錨泊的船舶、礁石等。動(dòng)態(tài)信息一般指航行環(huán)境中的動(dòng)態(tài)障礙物實(shí)體和環(huán)境實(shí)體。

    通過對(duì)感知系統(tǒng)獲取的多源異構(gòu)信息進(jìn)行實(shí)體與屬性描述,建立無人駕駛船舶航行態(tài)勢(shì)本體模型,然后在此基礎(chǔ)上結(jié)合《規(guī)則》及良好的船藝要求,將會(huì)遇場(chǎng)景進(jìn)行量化劃分,并基于prolog 語言對(duì)會(huì)遇場(chǎng)景進(jìn)行表述。

    根據(jù)《規(guī)則》及駕駛實(shí)踐,本文以目標(biāo)船的方位關(guān)系劃分了6 個(gè)會(huì)遇區(qū)域,分別為船艏、船艉、左舷正橫前、左舷正橫后、右舷正橫前及右舷正橫后。在每個(gè)會(huì)遇區(qū)域,又根據(jù)目標(biāo)船與本船的航向關(guān)系以及航速關(guān)系,將會(huì)遇場(chǎng)景(ES)劃分為12 種,然后基于船舶航行態(tài)勢(shì)本體模型及prolog語言建立航行態(tài)勢(shì)理解與劃分表,如表1 所示。表中:HO 為對(duì)遇場(chǎng)景;DHO為危險(xiǎn)對(duì)遇場(chǎng)景;CR為交叉相遇場(chǎng)景;OT 為追越場(chǎng)景;q為目標(biāo)船的舷角;Ψt為目標(biāo)船航向;Ψ0為本船航向;Vt為目標(biāo)船航速;V0為本船航速。

    本文針對(duì)海上環(huán)境復(fù)雜多變的特點(diǎn),結(jié)合駕駛實(shí)踐,提出基于航行態(tài)勢(shì)本體模型的場(chǎng)景劃分方法,對(duì)航行環(huán)境降維處理,實(shí)現(xiàn)了無人駕駛船舶對(duì)航行態(tài)勢(shì)的理解與劃分,為智能避碰決策提供基礎(chǔ),并且在一定程度上提高了決策效率。

    表1 航行態(tài)勢(shì)理解與劃分表Table 1 Navigation situation understanding and division

    2.2 多元復(fù)合避碰危險(xiǎn)度評(píng)估模型

    《規(guī)則》沒有給出碰撞危險(xiǎn)的定義,但其中很多條款是以構(gòu)成碰撞危險(xiǎn)為前提的,故需建立一定的模型將其量化。

    在航海實(shí)踐中,駕駛員在判斷是否存在碰撞危險(xiǎn)的同時(shí),也就對(duì)存在碰撞危險(xiǎn)的緊迫程度有了直觀的感受,以及定性與定量的分析,一般習(xí)慣于使用DCPA、兩船間距及其變化、TCPA來判斷是否存在碰撞危險(xiǎn)以及碰撞危險(xiǎn)的程度。但影響船舶碰撞危險(xiǎn)度大小的因素很多,且相互制約,故本文在前人研究的基礎(chǔ)上,利用模糊理論建立了DCPA和TCPA、兩船距離、相對(duì)方位以及船速比這5 個(gè)影響因素的隸屬度函數(shù),并依據(jù)駕駛實(shí)踐構(gòu)建了多元復(fù)合碰撞危險(xiǎn)度評(píng)估模型。

    由文獻(xiàn)[8]可知,船舶間距離越小,船舶間碰撞危險(xiǎn)越大,且從避讓難度上看,距離越小,避讓效果越差。任意時(shí)刻船舶間距離的隸屬度函數(shù)ud為:

    船速比越大,碰撞危險(xiǎn)度越大,任意時(shí)刻船速比的隸屬度函數(shù)uK為:

    DCPA對(duì)碰撞危險(xiǎn)度的影響較大,且DCPA的值越小,碰撞危險(xiǎn)程度越高。任意時(shí)刻DCPA的隸屬度函數(shù)uDCPA為:

    式中,d2為最小安全會(huì)遇距離。為簡化計(jì)算,取d2=2 n mile,考慮到船舶尺度的影響,在DCPA很小時(shí)也會(huì)發(fā)生碰撞,故取兩船船長的平均值,一般近似取為0.1 n mile。

    TCPA是針對(duì)碰撞危險(xiǎn)程度的評(píng)價(jià)指標(biāo)。TCPA值越小,碰撞危險(xiǎn)程度越高,任意時(shí)刻TCPA的隸屬度函數(shù)uTCPA為:

    原模型對(duì)于CRI的計(jì)算采用

    根據(jù)駕駛實(shí)踐規(guī)律,若DCPA大于安全會(huì)遇距離,盡管TCPA很小,船舶仍無碰撞危險(xiǎn),故改進(jìn)的船舶碰撞危險(xiǎn)度為

    2.3 BOP 智能避碰決策模型

    在前面幾個(gè)模型的基礎(chǔ)上,船舶智能避碰決策主要是解決船舶的最佳避讓時(shí)機(jī)及避讓幅度問題。由于船舶采取避碰行動(dòng)的時(shí)機(jī)和幅度受《規(guī)則》的約束以及駕駛航行經(jīng)驗(yàn)的影響,不同駕駛員對(duì)于不同的會(huì)遇情況采取的行動(dòng)不一致,在不同的時(shí)機(jī)采取避碰措施會(huì)產(chǎn)生不同的避讓效果,如圖3 所示。

    圖3 中:O點(diǎn)為本船當(dāng)前船位;A點(diǎn)為目標(biāo)船當(dāng)前船位;R1,R2,R3為船舶的不同避讓決策時(shí)機(jī);P1,P2,P3為不同的最近會(huì)遇點(diǎn); ?ψi(i=1,2,3)為對(duì)應(yīng)于不同避讓決策時(shí)機(jī)所采取的不同避讓行動(dòng)幅度,即航向變化量;A1,A2,A3為目標(biāo)船轉(zhuǎn)向Δψ1,Δψ2,Δψ3的船舶位置;TC為無人駕駛船舶的真航向;AiO為兩船相對(duì)航向。

    圖3 船舶轉(zhuǎn)向避碰行動(dòng)方式示意圖Fig. 3 Schematic diagram of ship steering and collision avoidance action mode

    本文將從駕駛員思維決策的角度出發(fā),通過分析本船與目標(biāo)船之間的相對(duì)運(yùn)動(dòng)幾何關(guān)系,在保證船舶能在安全會(huì)遇距離駛過讓清的基礎(chǔ)上建立BOP 智能避碰決策模型,以船舶避碰總路徑最短為目標(biāo)函數(shù),船舶操縱性能及《規(guī)則》為約束條件,量化避碰轉(zhuǎn)向區(qū)間,求解最優(yōu)的避碰策略。BOP智能避碰決策模型流程圖如圖4 所示。

    當(dāng)船舶存在碰撞危險(xiǎn)時(shí),需要根據(jù)船舶安全會(huì)遇距離,通過幾何解算出避讓幅度。在不同的位置采取避碰措施將導(dǎo)致船舶避讓路徑長度不同,為求取不同轉(zhuǎn)向幅度下的船舶避碰最優(yōu)路徑,建立以下目標(biāo)函數(shù):

    式中:xk,yk分別為船舶當(dāng)前的船位;xk+1,yk+1分別為船舶下一步的船位。

    因避碰決策時(shí)機(jī)與避碰路徑的對(duì)應(yīng)關(guān)系不易求解出解析式,故采用避碰行動(dòng)距離量化模型,來求解每個(gè)避碰距離下的轉(zhuǎn)向幅度以及避碰路徑。根據(jù)駕駛實(shí)踐,對(duì)于正橫前的來船,一般在兩船相距4~6 n mile 時(shí)采取避碰行動(dòng);對(duì)于正橫后的來船,在兩船相距3~4 n mile 時(shí)采取避碰行動(dòng)。本文對(duì)避碰轉(zhuǎn)向區(qū)間進(jìn)行量化,按步幅為0.1 n mile進(jìn)行取值,建立避讓決策時(shí)機(jī) (R1,R2,···,Rn),避讓行動(dòng)幅度一維向量為 (?ψ1,?ψ2,···,?ψn)。

    在計(jì)算避碰路徑時(shí),考慮到國際海上避碰規(guī)則要求以及船舶操縱性限制,需滿足以下約束條件:

    圖4 BOP 智能避碰決策模型流程圖Fig. 4 Flow chart of BOP intelligent collision avoidance decisionmaking model

    在駕駛實(shí)踐中,要求駕駛員對(duì)船舶性能有充分的了解,在做決策時(shí),要考慮船舶的轉(zhuǎn)舵時(shí)間以及旋回特性,鑒于此,本文考慮采用Abkowitz提出的三自由度整體型模型[10]:

    式中:X,Y為外力在x,y軸上的分量;m為船舶質(zhì)量;Izz為船體繞通過重心的鉛直軸的質(zhì)量慣性矩(轉(zhuǎn)動(dòng)慣量);N為作用于船舶的外力繞通過重心的鉛直軸的回轉(zhuǎn)力矩。

    針對(duì)水動(dòng)力導(dǎo)數(shù)的估算方法不作為本文的重點(diǎn),故此處不再贅述。本文的船舶操縱數(shù)學(xué)模型選取ITTC 典型的Mariner 標(biāo)準(zhǔn)船型作為研究對(duì)象,其船型數(shù)據(jù)如表2 所示,利用文獻(xiàn)[21]中的無量綱水動(dòng)力數(shù)據(jù)表,即可求得船舶的運(yùn)動(dòng)模型。

    表2 Mariner 標(biāo)準(zhǔn)船型數(shù)據(jù)Table 2 Coefficients in mathematical model of ship maneuvering motion of Mariner

    為貼合航海實(shí)際,在設(shè)計(jì)船舶轉(zhuǎn)向控制器時(shí)要考慮動(dòng)舵合理、航向變化快速平穩(wěn)、調(diào)節(jié)時(shí)間短,故本文采用PD 改向控制器模型。δ 的控制率為Cδ=?Kp(Ψ?Ψdesired)+Kdr,其中Ψ為船舶的當(dāng)前航向,Ψdesired為決策航向,對(duì)船舶PD 控制模型進(jìn)行參數(shù)整定,比例增益Kp= 0.3,微分增益Kd=1。

    船舶在避讓操作之后,需進(jìn)行復(fù)航操縱。在研究船舶復(fù)航時(shí)機(jī)的過程中,對(duì)于避碰行動(dòng)后是恢復(fù)至船舶初始運(yùn)動(dòng)狀態(tài)還是恢復(fù)到原航線上一直留有爭議。但就避碰而言,只需考慮恢復(fù)初始運(yùn)動(dòng)狀態(tài)即可,恢復(fù)原航線是航行自動(dòng)化系統(tǒng)的功能和任務(wù)。文獻(xiàn)[22]給出了常用的4 種恢復(fù)航向的時(shí)機(jī)模型:

    1) 當(dāng)目標(biāo)船駛過本船正橫處時(shí)復(fù)航;

    2) 當(dāng)讓路船復(fù)航后,保證船舶間的DCPA大于最小安全會(huì)遇距離;

    3) 當(dāng)目標(biāo)船駛過最近會(huì)遇距離點(diǎn)時(shí)復(fù)航;

    4) 當(dāng)目標(biāo)船駛過本船船艏后復(fù)航。

    無論哪種復(fù)航?jīng)Q策方式,都要保證兩船之間有一定的安全會(huì)遇距離。本文根據(jù)幾何避碰原理,解算了船舶預(yù)測(cè)復(fù)航時(shí)機(jī),圖5 所示為船舶復(fù)航時(shí)機(jī)求解示意圖。

    圖5 所示為兩船會(huì)遇,本船船首向上情況下,本船改向?Ψ避讓目標(biāo)船,R為避讓決策時(shí)機(jī),G為預(yù)測(cè)復(fù)航點(diǎn),RML1 為船舶避碰操縱前的相對(duì)運(yùn)動(dòng)路徑,RML2 為船舶改向后的相對(duì)運(yùn)動(dòng)路徑,RML1'為RML1 的平行路徑且與船舶的最小安全會(huì)遇距離圓相切,G點(diǎn)為RML2 與RML1'的交點(diǎn),船舶從R到G航行所用時(shí)間即為船舶的復(fù)航時(shí)機(jī)。計(jì)算出船舶復(fù)航時(shí)間點(diǎn)以及船舶復(fù)航時(shí)應(yīng)走航向,利用PD 航向控制器控制船舶完成復(fù)航操縱。

    圖5 船舶復(fù)航時(shí)機(jī)求解示意圖Fig. 5 Schematic diagram of solving the timing of re-sailing

    3 仿真試驗(yàn)及分析

    本節(jié)以案例仿真實(shí)驗(yàn)來驗(yàn)證基于駕駛實(shí)踐的無人船智能避碰決策研究方法的有效性。在仿真實(shí)驗(yàn)中,本船采用ITTC 典型的Mariner 標(biāo)準(zhǔn)船型參數(shù),借助Matlab 2016a,Simulink 仿真平臺(tái)與python3.6 設(shè)計(jì)仿真對(duì)比實(shí)驗(yàn),并與未考慮船舶操縱運(yùn)動(dòng)模型且利用探測(cè)圓(半徑取4 n mile)確定避碰時(shí)機(jī)的傳統(tǒng)避碰算法進(jìn)行比較,在典型會(huì)遇場(chǎng)景下驗(yàn)證算法的適用性及優(yōu)越性。無人駕駛船舶與動(dòng)態(tài)目標(biāo)船舶的初始航行參數(shù)如表3 所示,仿真步長取6 s。

    表3 船舶初始航行參數(shù)設(shè)置Table 3 Ship's initial navigation state data

    仿真目標(biāo)為無人駕駛船舶能夠趨向目標(biāo)點(diǎn)并避讓動(dòng)態(tài)目標(biāo)船舶。當(dāng)環(huán)境中不存在障礙物或障礙物不在安全會(huì)遇距離內(nèi)時(shí),無人駕駛船舶保持航向或者處于復(fù)航階段,調(diào)整方向朝著目標(biāo)點(diǎn)前進(jìn);當(dāng)障礙物出現(xiàn)在安全會(huì)遇距離內(nèi)時(shí),無人駕駛船舶利用本文算法及傳統(tǒng)避碰算法進(jìn)行實(shí)時(shí)避障;當(dāng)無人駕駛船舶接近目標(biāo)點(diǎn)時(shí),仿真結(jié)束。

    在避碰仿真過程中,3 艘動(dòng)態(tài)目標(biāo)船舶保向保速航行并先后與無人駕駛船舶構(gòu)成交叉相遇局面、對(duì)遇局面以及追越局面。圖6 所示為截取的船舶避碰仿真過程關(guān)鍵片段。圖中,軌跡交叉部分為無人駕駛船舶避讓后穿過動(dòng)態(tài)目標(biāo)船舶后方時(shí)的歷史軌跡。

    對(duì)于動(dòng)態(tài)目標(biāo)船1,在T=1 000 s 時(shí),無人駕駛船舶判斷目標(biāo)船舶進(jìn)入到本船的決策領(lǐng)域,分析目標(biāo)船的航向及方位信息,經(jīng)過航向態(tài)勢(shì)理解,決策為正橫后交叉會(huì)遇態(tài)勢(shì),本船為讓路船,采取右轉(zhuǎn)的避讓行動(dòng)方式,通過基于駕駛員視角的智能避碰決策方法求解船舶的最優(yōu)避碰幅度。在T=2 000 s 時(shí),無人駕駛船舶與目標(biāo)船2 構(gòu)成對(duì)遇局面,本船為讓路船,采取右轉(zhuǎn)的避讓行動(dòng)方式,通過算法求解船舶的最優(yōu)避碰幅度。在T=4 500 s 時(shí),無人駕駛船舶與目標(biāo)船構(gòu)成追越局面,本船為追越船,目標(biāo)船為讓路船,由于目標(biāo)點(diǎn)在右側(cè),所以選擇從右舷追越?;趥鹘y(tǒng)避碰模型決策得到的避讓時(shí)機(jī)固定,且均比本文算法早,轉(zhuǎn)向幅度較小,船舶路徑不光滑,船舶避碰總路徑為32.952 n mile,而基于駕駛員視角的智能避碰決策方法求解的船舶避碰總路徑為25.022 n mile,較大程度縮短了航程,提升了避讓性能。

    通過仿真可知,無人駕駛船舶在與目標(biāo)船舶構(gòu)成碰撞危險(xiǎn)時(shí),能夠有效做出避碰決策,并且通過智能避碰決策算法能夠給出最優(yōu)的避碰路徑。

    圖7 所示為無人駕駛船舶在避碰過程中的航向變化以及舵角變化情況。從圖中可知,航向變化快速平穩(wěn),調(diào)節(jié)時(shí)間短,震蕩較小,打舵平滑,動(dòng)舵合理,符合航海實(shí)踐。

    由仿真實(shí)驗(yàn)結(jié)果可知,無人駕駛船舶在不同的會(huì)遇態(tài)勢(shì)下均能安全避讓目標(biāo)船舶,并根據(jù)不同的會(huì)遇態(tài)勢(shì)生成最優(yōu)避碰策略,驗(yàn)證了本文所提算法的有效性。

    圖6 船舶避碰仿真過程圖Fig. 6 Ship collision avoidance simulation process diagram

    圖7 船舶避碰過程中航向舵角變化曲線Fig. 7 Course and rudder angle variation curve during collision avoidance

    4 結(jié) 語

    通過對(duì)無人駕駛船舶避碰決策所面臨問題的分析,提出了一種基于駕駛實(shí)踐的智能避碰決策方法,可使無人駕駛船舶根據(jù)感知系統(tǒng)傳來的信息,實(shí)時(shí)與環(huán)境交互,做出符合《規(guī)則》的避碰決策和路徑規(guī)劃。在典型會(huì)遇場(chǎng)景下進(jìn)行的仿真實(shí)驗(yàn)表明,該模型能夠在未知環(huán)境下完成自主避碰,并解算出最優(yōu)路徑,證明了算法的有效性。在接下來的研究中,還需要對(duì)模型不斷完善,考慮更為復(fù)雜的環(huán)境擾動(dòng),加入對(duì)靜、動(dòng)態(tài)障礙物領(lǐng)域模型,更加貼合航海實(shí)際。

    猜你喜歡
    船舶智能模型
    一半模型
    計(jì)算流體力學(xué)在船舶操縱運(yùn)動(dòng)仿真中的應(yīng)用
    《船舶》2022 年度征訂啟事
    船舶(2021年4期)2021-09-07 17:32:22
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    船舶!請(qǐng)加速
    智能前沿
    文苑(2018年23期)2018-12-14 01:06:06
    智能前沿
    文苑(2018年19期)2018-11-09 01:30:14
    智能前沿
    文苑(2018年17期)2018-11-09 01:29:26
    智能前沿
    文苑(2018年21期)2018-11-09 01:22:32
    亚洲精品久久国产高清桃花| 亚洲人与动物交配视频| 色综合亚洲欧美另类图片| 色哟哟哟哟哟哟| 亚洲色图av天堂| 精品久久久久久久人妻蜜臀av| 久久精品国产亚洲av天美| 久久精品国产亚洲av天美| 国产精品99久久久久久久久| 国产高清视频在线播放一区| 丰满人妻一区二区三区视频av| 色播亚洲综合网| 麻豆av噜噜一区二区三区| 黄色欧美视频在线观看| 欧美潮喷喷水| 免费观看精品视频网站| 欧美xxxx性猛交bbbb| 1024手机看黄色片| 夜夜爽天天搞| 麻豆国产97在线/欧美| 美女免费视频网站| 97人妻精品一区二区三区麻豆| 日本免费一区二区三区高清不卡| ponron亚洲| 久久久久久久久久黄片| 男人和女人高潮做爰伦理| 天堂影院成人在线观看| 国产极品精品免费视频能看的| 色综合亚洲欧美另类图片| 免费不卡的大黄色大毛片视频在线观看 | 日韩欧美精品免费久久| 国产精品日韩av在线免费观看| 午夜福利欧美成人| 少妇被粗大猛烈的视频| 日本免费a在线| 成年免费大片在线观看| 免费看光身美女| 亚洲乱码一区二区免费版| 色综合站精品国产| av女优亚洲男人天堂| 老师上课跳d突然被开到最大视频| 少妇的逼水好多| 国产精品一区二区免费欧美| 婷婷丁香在线五月| 精品免费久久久久久久清纯| 亚洲美女黄片视频| 亚洲经典国产精华液单| 伦精品一区二区三区| 人人妻人人看人人澡| 国产免费男女视频| 黄片wwwwww| 日韩欧美免费精品| 久久午夜福利片| 男人和女人高潮做爰伦理| 亚洲五月天丁香| 日韩欧美精品v在线| 亚洲精品一区av在线观看| 国产爱豆传媒在线观看| 色哟哟·www| 久9热在线精品视频| 国产又黄又爽又无遮挡在线| 在线a可以看的网站| 国产午夜精品久久久久久一区二区三区 | 欧美性猛交╳xxx乱大交人| avwww免费| 在线观看av片永久免费下载| 欧美区成人在线视频| 黄色配什么色好看| 欧美激情久久久久久爽电影| 网址你懂的国产日韩在线| 亚洲av中文av极速乱 | 麻豆一二三区av精品| 一进一出抽搐gif免费好疼| 国产欧美日韩精品亚洲av| 久久久色成人| 亚洲精华国产精华液的使用体验 | aaaaa片日本免费| av女优亚洲男人天堂| 亚洲精品国产成人久久av| 黄色女人牲交| 亚洲成人中文字幕在线播放| 97碰自拍视频| 美女免费视频网站| 我要看日韩黄色一级片| 极品教师在线视频| 国产aⅴ精品一区二区三区波| 免费人成在线观看视频色| 又黄又爽又刺激的免费视频.| 亚洲一区高清亚洲精品| 成人性生交大片免费视频hd| 美女 人体艺术 gogo| 亚洲人与动物交配视频| 国产成人av教育| 内地一区二区视频在线| 永久网站在线| av在线天堂中文字幕| 欧美3d第一页| 又紧又爽又黄一区二区| 男女边吃奶边做爰视频| 国产单亲对白刺激| 一个人看视频在线观看www免费| 一级黄色大片毛片| 不卡一级毛片| 亚洲专区国产一区二区| 国产伦一二天堂av在线观看| 免费人成在线观看视频色| 国产 一区 欧美 日韩| 午夜久久久久精精品| 成人国产综合亚洲| 国产精品国产三级国产av玫瑰| 美女cb高潮喷水在线观看| 老师上课跳d突然被开到最大视频| 欧美日本视频| 日本成人三级电影网站| 不卡一级毛片| 久久热精品热| 在线a可以看的网站| 天天一区二区日本电影三级| 亚洲精华国产精华液的使用体验 | 国产主播在线观看一区二区| 亚洲国产欧洲综合997久久,| 欧美丝袜亚洲另类 | 欧美成人a在线观看| 老师上课跳d突然被开到最大视频| 身体一侧抽搐| а√天堂www在线а√下载| 国产久久久一区二区三区| 亚洲成av人片在线播放无| 日本黄色视频三级网站网址| 国产极品精品免费视频能看的| 欧美日韩国产亚洲二区| 亚洲av中文字字幕乱码综合| 日韩一区二区视频免费看| 日韩欧美在线乱码| 亚洲专区国产一区二区| 色综合站精品国产| 国内精品久久久久久久电影| 三级男女做爰猛烈吃奶摸视频| 日韩欧美一区二区三区在线观看| 亚洲狠狠婷婷综合久久图片| 啪啪无遮挡十八禁网站| 亚洲黑人精品在线| 免费搜索国产男女视频| 精品久久久久久久末码| 高清毛片免费观看视频网站| 乱系列少妇在线播放| 成熟少妇高潮喷水视频| 亚洲在线自拍视频| 国产精品一区二区三区四区免费观看 | 在线观看66精品国产| 天堂动漫精品| 婷婷丁香在线五月| 国内少妇人妻偷人精品xxx网站| 嫩草影院精品99| 一本一本综合久久| 内地一区二区视频在线| 亚洲在线自拍视频| 久久精品国产亚洲av天美| 在线播放国产精品三级| 又黄又爽又刺激的免费视频.| 两性午夜刺激爽爽歪歪视频在线观看| 熟女电影av网| 国产乱人视频| 亚洲18禁久久av| 夜夜爽天天搞| 极品教师在线视频| 国产午夜精品论理片| 欧美日韩瑟瑟在线播放| 久久久久免费精品人妻一区二区| 国产av麻豆久久久久久久| 国内精品一区二区在线观看| 国产熟女欧美一区二区| av中文乱码字幕在线| 精品乱码久久久久久99久播| av女优亚洲男人天堂| 色综合婷婷激情| 亚洲成人中文字幕在线播放| 精品一区二区三区视频在线| 精品福利观看| 国产成人福利小说| 欧美国产日韩亚洲一区| 久久热精品热| 我的女老师完整版在线观看| 一区二区三区免费毛片| 97超级碰碰碰精品色视频在线观看| 热99re8久久精品国产| 久久久国产成人精品二区| 国产综合懂色| 午夜免费男女啪啪视频观看 | 欧美一级a爱片免费观看看| 一区二区三区高清视频在线| 久久6这里有精品| 级片在线观看| 一个人看视频在线观看www免费| 免费在线观看影片大全网站| 97人妻精品一区二区三区麻豆| 欧美黑人欧美精品刺激| 最近视频中文字幕2019在线8| 久久6这里有精品| 国内精品宾馆在线| 蜜桃久久精品国产亚洲av| 日韩在线高清观看一区二区三区 | 老熟妇仑乱视频hdxx| 日韩欧美一区二区三区在线观看| 熟女人妻精品中文字幕| 免费av不卡在线播放| 久久久久久久午夜电影| 国产v大片淫在线免费观看| 国产伦一二天堂av在线观看| 亚洲精品亚洲一区二区| 国产精品国产三级国产av玫瑰| 亚洲熟妇中文字幕五十中出| 亚洲最大成人av| 少妇猛男粗大的猛烈进出视频 | 婷婷丁香在线五月| 波野结衣二区三区在线| 综合色av麻豆| 天堂影院成人在线观看| 亚洲精品日韩av片在线观看| 无人区码免费观看不卡| 日本黄大片高清| 91在线精品国自产拍蜜月| 成人亚洲精品av一区二区| 黄色配什么色好看| 精品人妻偷拍中文字幕| 国产伦精品一区二区三区视频9| 色综合站精品国产| 天天躁日日操中文字幕| 欧美国产日韩亚洲一区| 国产精品不卡视频一区二区| 999久久久精品免费观看国产| 男女下面进入的视频免费午夜| 日韩强制内射视频| 能在线免费观看的黄片| 女的被弄到高潮叫床怎么办 | 亚洲av日韩精品久久久久久密| 欧美另类亚洲清纯唯美| 国产真实伦视频高清在线观看 | 亚洲精品在线观看二区| 乱人视频在线观看| 麻豆精品久久久久久蜜桃| 老女人水多毛片| 看片在线看免费视频| 国产伦在线观看视频一区| 亚洲色图av天堂| 18禁黄网站禁片免费观看直播| 日韩欧美免费精品| 一级黄片播放器| 国产精品不卡视频一区二区| 亚洲精品亚洲一区二区| 少妇猛男粗大的猛烈进出视频 | 午夜亚洲福利在线播放| 午夜免费男女啪啪视频观看 | 日本一二三区视频观看| 特级一级黄色大片| 日韩中字成人| 色播亚洲综合网| 在线播放国产精品三级| 欧美激情国产日韩精品一区| 在线观看免费视频日本深夜| 久久精品国产亚洲av香蕉五月| 美女xxoo啪啪120秒动态图| av女优亚洲男人天堂| 国产色爽女视频免费观看| 亚洲中文字幕日韩| 国产一区二区激情短视频| 国内久久婷婷六月综合欲色啪| 午夜福利在线在线| 在现免费观看毛片| a级一级毛片免费在线观看| 欧美又色又爽又黄视频| 女生性感内裤真人,穿戴方法视频| 99精品久久久久人妻精品| 亚洲四区av| 精品久久久久久成人av| 免费看光身美女| 欧美日本亚洲视频在线播放| 丰满乱子伦码专区| 简卡轻食公司| 免费av观看视频| 色综合亚洲欧美另类图片| 干丝袜人妻中文字幕| 国内少妇人妻偷人精品xxx网站| 亚洲精品久久国产高清桃花| 国产精品1区2区在线观看.| 深爱激情五月婷婷| 尾随美女入室| 最近视频中文字幕2019在线8| 国产综合懂色| 又粗又爽又猛毛片免费看| 一级毛片久久久久久久久女| 免费观看的影片在线观看| 久久久国产成人精品二区| 免费黄网站久久成人精品| 久久久久久久久中文| 我的老师免费观看完整版| 国产精品,欧美在线| 18禁裸乳无遮挡免费网站照片| 国产色婷婷99| 97超视频在线观看视频| 国内久久婷婷六月综合欲色啪| 亚洲av中文av极速乱 | 亚洲人成网站高清观看| 97碰自拍视频| 亚洲国产欧美人成| 毛片一级片免费看久久久久 | 久久99热6这里只有精品| 亚洲av中文字字幕乱码综合| 日本与韩国留学比较| 欧美性猛交黑人性爽| 午夜激情福利司机影院| 欧美精品啪啪一区二区三区| 久久久久久伊人网av| 人妻久久中文字幕网| 亚洲乱码一区二区免费版| 狂野欧美白嫩少妇大欣赏| 999久久久精品免费观看国产| 精品久久久久久久末码| 国产精品国产高清国产av| or卡值多少钱| 级片在线观看| 国产精品久久久久久亚洲av鲁大| 一个人免费在线观看电影| 久久精品国产鲁丝片午夜精品 | 国产爱豆传媒在线观看| 日韩欧美三级三区| 亚洲国产色片| 看黄色毛片网站| av国产免费在线观看| 久久久久国内视频| 在线国产一区二区在线| 我的女老师完整版在线观看| 美女高潮的动态| 99视频精品全部免费 在线| 国产男人的电影天堂91| 亚洲av免费高清在线观看| 乱码一卡2卡4卡精品| 亚洲成av人片在线播放无| 国产爱豆传媒在线观看| 欧美一区二区亚洲| 成年女人毛片免费观看观看9| 国产毛片a区久久久久| 九九久久精品国产亚洲av麻豆| 高清毛片免费观看视频网站| 国产伦精品一区二区三区四那| 免费不卡的大黄色大毛片视频在线观看 | 久久精品国产99精品国产亚洲性色| 欧美一区二区精品小视频在线| 亚洲av免费高清在线观看| 少妇熟女aⅴ在线视频| 国产亚洲91精品色在线| 久久久色成人| 亚洲第一电影网av| 国产探花极品一区二区| 亚洲国产精品sss在线观看| 69av精品久久久久久| 亚洲图色成人| 精品一区二区三区人妻视频| or卡值多少钱| 蜜桃久久精品国产亚洲av| 亚洲av中文av极速乱 | 精品久久久久久久久久免费视频| 男女啪啪激烈高潮av片| 99国产精品一区二区蜜桃av| 精品久久久久久久久久免费视频| 美女大奶头视频| 欧美丝袜亚洲另类 | 蜜桃亚洲精品一区二区三区| 一级a爱片免费观看的视频| 成年女人毛片免费观看观看9| 国产真实伦视频高清在线观看 | 色吧在线观看| 天堂网av新在线| 国产亚洲91精品色在线| 亚洲av中文字字幕乱码综合| 乱人视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 不卡视频在线观看欧美| 色播亚洲综合网| 一个人看的www免费观看视频| 色播亚洲综合网| 亚洲经典国产精华液单| 嫁个100分男人电影在线观看| 亚洲成人免费电影在线观看| 99在线视频只有这里精品首页| 99久久久亚洲精品蜜臀av| 国产人妻一区二区三区在| 日本免费一区二区三区高清不卡| aaaaa片日本免费| 日本五十路高清| 欧美丝袜亚洲另类 | 国产大屁股一区二区在线视频| 日本撒尿小便嘘嘘汇集6| 国产av在哪里看| 美女高潮的动态| 久久国产乱子免费精品| 九九爱精品视频在线观看| 99久久精品热视频| 精品久久久久久久久久免费视频| 有码 亚洲区| 麻豆一二三区av精品| 国产一区二区三区在线臀色熟女| 搞女人的毛片| 嫩草影院入口| 国产极品精品免费视频能看的| 精品一区二区三区视频在线| 人人妻,人人澡人人爽秒播| 很黄的视频免费| 国产午夜精品久久久久久一区二区三区 | 国产亚洲av嫩草精品影院| 色哟哟哟哟哟哟| 久99久视频精品免费| 神马国产精品三级电影在线观看| 国产精品伦人一区二区| 九九久久精品国产亚洲av麻豆| 国产大屁股一区二区在线视频| 欧美黑人欧美精品刺激| 午夜免费男女啪啪视频观看 | 欧美3d第一页| 久久亚洲精品不卡| 国产成人a区在线观看| 国产国拍精品亚洲av在线观看| 成人精品一区二区免费| 国产精品精品国产色婷婷| 日本熟妇午夜| 18禁裸乳无遮挡免费网站照片| 午夜福利在线观看免费完整高清在 | 在线免费观看的www视频| 亚洲人成网站高清观看| 亚洲精品国产成人久久av| 国产伦在线观看视频一区| 99国产精品一区二区蜜桃av| 夜夜看夜夜爽夜夜摸| 麻豆国产av国片精品| 又粗又爽又猛毛片免费看| 国产伦精品一区二区三区视频9| 久久久久性生活片| 午夜福利欧美成人| 日韩欧美国产在线观看| 国产一区二区激情短视频| 尾随美女入室| 男女视频在线观看网站免费| 欧美三级亚洲精品| 中文字幕精品亚洲无线码一区| 国产真实伦视频高清在线观看 | 天堂√8在线中文| 午夜日韩欧美国产| 69av精品久久久久久| 亚洲人成网站在线播放欧美日韩| 熟女人妻精品中文字幕| 国产一区二区三区av在线 | 免费在线观看成人毛片| av国产免费在线观看| 国产av一区在线观看免费| 天堂动漫精品| netflix在线观看网站| 亚洲欧美日韩卡通动漫| 精品一区二区三区视频在线观看免费| 一进一出抽搐动态| 搡老岳熟女国产| 中文字幕人妻熟人妻熟丝袜美| 久久人人爽人人爽人人片va| 亚洲国产日韩欧美精品在线观看| 女人十人毛片免费观看3o分钟| 国产视频一区二区在线看| 国产淫片久久久久久久久| 国产av在哪里看| 欧美极品一区二区三区四区| 亚洲av免费在线观看| 看黄色毛片网站| 欧美成人性av电影在线观看| 日本-黄色视频高清免费观看| 无人区码免费观看不卡| 久久精品国产自在天天线| 午夜久久久久精精品| 国产一区二区激情短视频| 亚洲,欧美,日韩| 午夜福利在线观看吧| 1000部很黄的大片| 男人狂女人下面高潮的视频| 国产精品免费一区二区三区在线| 天堂av国产一区二区熟女人妻| 国产精品女同一区二区软件 | 午夜爱爱视频在线播放| 一级a爱片免费观看的视频| 国产综合懂色| 久久久精品欧美日韩精品| 亚洲av五月六月丁香网| 一本一本综合久久| 国产真实伦视频高清在线观看 | 琪琪午夜伦伦电影理论片6080| 国产探花极品一区二区| 久久国内精品自在自线图片| 神马国产精品三级电影在线观看| 亚洲国产精品久久男人天堂| 亚洲熟妇熟女久久| 日韩,欧美,国产一区二区三区 | 久久久午夜欧美精品| 欧美成人免费av一区二区三区| 在线免费观看的www视频| 免费av观看视频| 99久久无色码亚洲精品果冻| 91久久精品电影网| 精品日产1卡2卡| 亚洲国产欧洲综合997久久,| 特级一级黄色大片| 久久国内精品自在自线图片| 午夜福利18| 国产大屁股一区二区在线视频| or卡值多少钱| 99久久九九国产精品国产免费| 亚洲国产精品成人综合色| 99久久九九国产精品国产免费| 毛片女人毛片| 村上凉子中文字幕在线| 99久国产av精品| 少妇被粗大猛烈的视频| 一本精品99久久精品77| 乱人视频在线观看| 国产熟女欧美一区二区| 国产亚洲欧美98| 一进一出好大好爽视频| 在线国产一区二区在线| 一区二区三区免费毛片| 亚洲精品在线观看二区| 级片在线观看| 极品教师在线免费播放| h日本视频在线播放| 亚洲成人久久爱视频| avwww免费| 日日摸夜夜添夜夜添小说| 欧美zozozo另类| 成人av一区二区三区在线看| 老司机福利观看| 床上黄色一级片| 国产精品一区二区三区四区免费观看 | 三级毛片av免费| 亚洲精品日韩av片在线观看| 欧美成人免费av一区二区三区| 久久精品久久久久久噜噜老黄 | 国产 一区 欧美 日韩| 美女黄网站色视频| 亚洲在线自拍视频| 国产黄a三级三级三级人| 成人精品一区二区免费| 俄罗斯特黄特色一大片| 少妇熟女aⅴ在线视频| 亚洲五月天丁香| 99视频精品全部免费 在线| 亚洲成人久久爱视频| 搡女人真爽免费视频火全软件 | 亚洲国产精品合色在线| 人妻少妇偷人精品九色| 最新中文字幕久久久久| 小蜜桃在线观看免费完整版高清| 国模一区二区三区四区视频| 人人妻人人澡欧美一区二区| 真实男女啪啪啪动态图| 日本a在线网址| 不卡一级毛片| 国产v大片淫在线免费观看| АⅤ资源中文在线天堂| 亚洲av五月六月丁香网| 中文字幕精品亚洲无线码一区| 免费av毛片视频| 韩国av一区二区三区四区| 国产亚洲精品久久久com| 日本色播在线视频| 啦啦啦韩国在线观看视频| 午夜亚洲福利在线播放| 国产高潮美女av| 亚洲国产精品久久男人天堂| 听说在线观看完整版免费高清| 一区二区三区激情视频| 91麻豆av在线| 国国产精品蜜臀av免费| 在线a可以看的网站| 麻豆av噜噜一区二区三区| 天堂动漫精品| 白带黄色成豆腐渣| 国内精品久久久久久久电影| 亚洲av免费在线观看| 欧美精品啪啪一区二区三区| 一进一出抽搐动态| 国产免费男女视频| 在线免费十八禁| 99久久精品热视频| 欧美日韩国产亚洲二区| 最近在线观看免费完整版| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕av在线有码专区| 在线观看午夜福利视频| 国产三级在线视频| 免费观看精品视频网站| 亚洲精华国产精华精| 久久人人精品亚洲av| 啦啦啦韩国在线观看视频| 最近最新免费中文字幕在线| 村上凉子中文字幕在线| 啦啦啦韩国在线观看视频| 欧美xxxx黑人xx丫x性爽| 亚洲成人精品中文字幕电影| 大又大粗又爽又黄少妇毛片口| 日本 欧美在线| 亚洲美女黄片视频| 欧美潮喷喷水| 欧美成人一区二区免费高清观看| 欧美又色又爽又黄视频| 精品久久久噜噜| 日本 欧美在线| 精品久久久久久,| av专区在线播放| 最后的刺客免费高清国语| 精华霜和精华液先用哪个|