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

    基于空間模糊層次分析的航行環(huán)境風(fēng)險(xiǎn)評估

    2023-05-30 10:48:04杜思雨諶頌程亮朱大明李澤明
    貴州大學(xué)學(xué)報(bào)(自然科學(xué)版) 2023年3期
    關(guān)鍵詞:廈門港航行海域

    杜思雨 諶頌 程亮 朱大明 李澤明

    摘 要:隨著海上交通日益繁忙,保障船舶安全航行成為全球性重要議題。本文從航行自然環(huán)境和海上交通環(huán)境兩方面構(gòu)建航行環(huán)境風(fēng)險(xiǎn)評估模型,將地理空間技術(shù)和模糊層次分析法(fuzzy analytic hierarchy process,F(xiàn)AHP)結(jié)合,對廈門港及其周邊海域的航行環(huán)境風(fēng)險(xiǎn)評估,形成風(fēng)險(xiǎn)分級圖并對其空間分布特征進(jìn)行分析。結(jié)果表明:受自然環(huán)境影響,該區(qū)域航行環(huán)境風(fēng)險(xiǎn)大致呈“U”型分布,風(fēng)險(xiǎn)等級沿“SW—NE”逐漸降低;在交通環(huán)境影響下,風(fēng)險(xiǎn)等級≥8級的區(qū)域呈散點(diǎn)狀分布在距海岸線較遠(yuǎn)海域,尤其是兄弟嶼東部海域;綜合以上兩方面,廈門港及其周邊海域航行環(huán)境風(fēng)險(xiǎn)整體處于中等偏低風(fēng)險(xiǎn),但仍有約6.17%的區(qū)域風(fēng)險(xiǎn)等級≥7級,主要分布在“南澎列島—兄弟嶼”連線的延長線以南和廈門灣海域。

    關(guān)鍵詞:航行環(huán)境風(fēng)險(xiǎn);地理空間技術(shù);模糊層次分析法;風(fēng)險(xiǎn)空間分布

    中圖分類號:P208;P732.1

    文獻(xiàn)標(biāo)志碼:A

    文章編號 1000-5269(2023)03-0054-08

    DOI:10.15958/j.cnki.gdxbzrb.2023.03.08

    隨著海上運(yùn)輸在全球貨運(yùn)、客運(yùn)中占據(jù)著愈發(fā)重要的地位,船舶數(shù)量日漸增多,港口水域的運(yùn)輸壓力日益加重 [1],從而導(dǎo)致近幾年海上交通事故頻繁發(fā)生,造成了嚴(yán)重的生命威脅和財(cái)物損失[2]。因此,航行風(fēng)險(xiǎn)評估是海上交通運(yùn)輸領(lǐng)域的重點(diǎn)研究課題。

    航行風(fēng)險(xiǎn)評估主要從指標(biāo)選取和權(quán)重確定方法兩方面出發(fā)。根據(jù)評估的側(cè)重點(diǎn)不同,評估指標(biāo)的選取上也存在較大差異[3]。CARMEN等[4]主要針對航道環(huán)境對大型集裝箱船的航行風(fēng)險(xiǎn)進(jìn)行了評估。張新放等[5]通過分析船舶因素、管理因素和人為因素,建立了基于綜合安全評估的指標(biāo)體系。張青等[6]將重點(diǎn)放在氣象因素的研究,對在不良天氣下航行的船舶航行有著一定的指導(dǎo)意義。李振福等[7]從政治、經(jīng)濟(jì)、社會方面對北極航行風(fēng)險(xiǎn)進(jìn)行了評估。在評估方法的選取上,1998年,馬會等[8]最先將灰色聚類和灰色統(tǒng)計(jì)結(jié)合起來對航道危險(xiǎn)度進(jìn)行了綜合評估。陳堯等[9]將模糊數(shù)學(xué)引入層次分析法中,利用模糊層次分析法對仰光河航道風(fēng)險(xiǎn)進(jìn)行了綜合評價。王晨等[10]建立了基于逼近理想解排序法(technique for order preference by similarity to ideal solution,TOPSIS)的航道引航風(fēng)險(xiǎn)評估模型對長江江蘇段4段航道的引航環(huán)境風(fēng)險(xiǎn)進(jìn)行了評價。以上研究多由觀測和統(tǒng)計(jì)數(shù)據(jù)確定評估指標(biāo)參數(shù),缺乏一定的客觀性,并且尚未考慮航行風(fēng)險(xiǎn)等級的空間差異,評估結(jié)果大多呈現(xiàn)為抽象的分值,并不能直觀地反映研究區(qū)的風(fēng)險(xiǎn)等級空間分布情況。

    近年來,遙感、全球衛(wèi)星導(dǎo)航定位(global positioning system,GPS)等地理大數(shù)據(jù)的日益增多,給海上航行風(fēng)險(xiǎn)提供了新視角。李振福等[11]通過船舶地理數(shù)據(jù)提取船舶的自身屬性對“21世紀(jì)海上絲綢之路”的航行風(fēng)險(xiǎn)進(jìn)行了評估。甄榮等[12]通過地理大數(shù)據(jù)研究了一種用于海上交通監(jiān)控管理的船舶會遇計(jì)算方法。然而,已有文獻(xiàn)對這些數(shù)據(jù)的量化形式缺乏詳盡考慮。

    基于此,本文將地理空間技術(shù)和模糊層次分析法相結(jié)合,基于海洋環(huán)境和船舶自動識別系統(tǒng)(Automatic Identification System,AIS)地理大數(shù)據(jù), 從航行自然環(huán)境和海上交通環(huán)境兩個方面構(gòu)建航行環(huán)境風(fēng)險(xiǎn)評估模型,并以廈門港及其周邊海域?yàn)槔?,對其航行環(huán)境風(fēng)險(xiǎn)進(jìn)行評估,實(shí)現(xiàn)風(fēng)險(xiǎn)等級空間分布的可視化,為船舶安全航行提供直觀的參考依據(jù)。

    1 研究區(qū)和數(shù)據(jù)

    1.1 研究區(qū)概況

    廈門港位于中國福建省東南部,是中國東南國際航運(yùn)中心、全國綜合運(yùn)輸體系的重要樞紐,也是閩南、閩西的主要出??诎叮?3]。港口轄區(qū)包含廈門市(東渡、海滄、翔安)和漳州市(招銀、后石、石碼、古雷、東山、詔安)的九大港區(qū)。2021年,廈門港集裝箱吞吐量達(dá)1 205萬標(biāo)箱,成為全球第十三大港口。2022年,增加集裝箱班輪航線,總計(jì)達(dá)到162條,港口調(diào)度和交通壓力明顯增大,船舶航行所面臨的風(fēng)險(xiǎn)也愈來愈大[14]。據(jù)中華人民共和國廈門海事局(http://port.xm.gov.cn/)發(fā)布的報(bào)告統(tǒng)計(jì)得,自2020起,轄區(qū)內(nèi)發(fā)生不同程度的事故和險(xiǎn)情共計(jì)53次,造成了不同程度的損失。廈門港地處東南亞熱帶季風(fēng)區(qū),季節(jié)性特征較為典型;沿海多霧且多為平流霧,每年的2月至5月為海霧多發(fā)期,霧數(shù)日占全年的60%以上;降水主要集中在每年的 5月至10月,其中,6月降雨量最大[15]。本文將廈門港及其周邊海域作為研究區(qū),如圖1所示。

    1.2 數(shù)據(jù)概況及預(yù)處理

    本文所使用到的數(shù)據(jù)主要包括船舶航行數(shù)據(jù)、海洋環(huán)境數(shù)據(jù)兩類,數(shù)據(jù)相關(guān)信息見表1。

    AIS數(shù)據(jù)可以客觀準(zhǔn)確地獲取船舶的航行信息和軌跡點(diǎn)位的空間分布,但是由于AIS數(shù)據(jù)中部分屬性需要人工輸入,所以不可避免地會出現(xiàn)錯誤。因此,為了保證數(shù)據(jù)質(zhì)量,確保后續(xù)試驗(yàn)的準(zhǔn)確性和客觀性,從字段屬性值和空間分布角度對AIS進(jìn)行缺失值補(bǔ)充、錯誤值去除或更新 [16]。海洋環(huán)境數(shù)據(jù)為.NC格式的文件,需要進(jìn)行預(yù)處理。首先將.NC文件全部轉(zhuǎn)為.TIF文件后,把所有.TIF文件重采樣至15弧秒,并裁剪至研究區(qū)范圍。最后將所有數(shù)據(jù)統(tǒng)一轉(zhuǎn)換至WGS-84坐標(biāo)系。

    2 研究方法

    從航行自然環(huán)境和海上交通環(huán)境兩方面構(gòu)建評估指標(biāo)體系,用模糊層析分析法確定權(quán)重,利用線性加權(quán)和地理空間技術(shù)對廈門港及其周邊海域的航行環(huán)境風(fēng)險(xiǎn)進(jìn)行評估,如圖2所示。

    2.1 指標(biāo)體系構(gòu)建

    在遵循指標(biāo)的可行性、時效性和數(shù)據(jù)可獲取性的前提下,從航行自然環(huán)境、海上交通環(huán)境兩方面,選取大風(fēng)頻數(shù)、大浪頻數(shù)、霧、降雨量、水深、船舶速度、交通密度、危險(xiǎn)船舶、助航標(biāo)志 9項(xiàng)影響因素作為航行風(fēng)險(xiǎn)評估指標(biāo),見表2。

    2.2 指標(biāo)值計(jì)算

    對9個指標(biāo)進(jìn)行數(shù)據(jù)處理,得到各指標(biāo)量化數(shù)值,具體處理過程見表2。

    由于每個指標(biāo)都有各自的量綱和特征,需進(jìn)行無量綱化處理,得到范圍在1~10的標(biāo)準(zhǔn)化指標(biāo)值。其中,正向和逆向指標(biāo)的無量綱化處理公式見式(1)[18-19]。式中,X′為標(biāo)準(zhǔn)化后的指標(biāo)值,X為原始指標(biāo)值,Xmin和Xmax分別為指標(biāo)值中的最小值和最大值。

    2.3 權(quán)重確定

    模糊層次分析法,是在層次分析法(analytic hierarchy process,AHP)的基礎(chǔ)上引入模糊數(shù)學(xué),最大程度降低了決策過程中的不確定性帶來的影響[20]。

    本文中采用文獻(xiàn)[21]中的矩陣解算方法代替?zhèn)鹘y(tǒng)方法[22]:

    1)根據(jù)專家提供的評價數(shù)據(jù)構(gòu)造三角模糊判斷矩陣A=(aij)n×n,矩陣中的元素aij=(lij,mij,uij)為三角模糊數(shù)[23](triangular fuzzy number,TFNs),在每一個三角模糊數(shù)(l,m,u)中,下界l、中值m、上界u的取值均依據(jù)層次分析法中的1~9標(biāo)度法確定[24];

    2)對判斷矩陣依照公式(2)、(3)進(jìn)行一致性檢驗(yàn),當(dāng)R<0.1時,矩陣滿足一致性要求,可以對矩陣進(jìn)行計(jì)算,否則,應(yīng)對元素aij加以調(diào)整,直至通過一致性檢驗(yàn);

    式中,C表示偏離一致性指標(biāo),I為隨機(jī)一致性指標(biāo),λmax表示矩陣的最大特征根;

    3)構(gòu)建模糊評判因子矩陣E,并計(jì)算調(diào)整判斷矩陣Q,如式(4)。其中,mij為三角模糊中值;

    4)將矩陣Q按列轉(zhuǎn)化為對角線為1的判斷矩陣;

    5)計(jì)算矩陣內(nèi)每行元素的n次方根,并對該數(shù)值進(jìn)行歸一化處理,即得到每項(xiàng)指標(biāo)的權(quán)重值。

    2.4 風(fēng)險(xiǎn)量化分級

    通過線性加權(quán)得到航行環(huán)境的風(fēng)險(xiǎn)值,如式(5)。利用相等間隔法將航行環(huán)境風(fēng)險(xiǎn)劃分為10個等級,航行環(huán)境風(fēng)險(xiǎn)從1級到10級依次增高。

    式中,Dk為第k個柵格單元的航行環(huán)境風(fēng)險(xiǎn)值;uik為指標(biāo)i在第k個柵格單元中的指標(biāo)值;Wi為指標(biāo)i的權(quán)重值;n為指標(biāo)的總個數(shù)。

    3 評估結(jié)果與分析

    3.1 指標(biāo)計(jì)算結(jié)果

    圖3顯示了各指標(biāo)的空間分布特征。在航行自然環(huán)境方面,研究區(qū)遠(yuǎn)海水域極易受到大風(fēng)、大浪的影響,從而影響船舶的操縱性,給船舶航行帶來一定的風(fēng)險(xiǎn),而近岸海域海霧覆蓋面積和降雨量則高于遠(yuǎn)海水域,能見度不佳,會使船舶航行安全度大幅降低。在水深指標(biāo)中,研究區(qū)近岸和島嶼處水深較淺,船舶航行會存在擱淺或觸礁風(fēng)險(xiǎn)。在海上交通環(huán)境方面,研究區(qū)內(nèi)船舶的航行速度整體較低,遠(yuǎn)海區(qū)域船速略高于近海區(qū)域。在交通密度指標(biāo)中,研究區(qū)內(nèi)的船舶主要集中在航道上,其中,東北部的進(jìn)港航道和港口內(nèi)部的船舶數(shù)量較多,交通密度較大時,很容易發(fā)生碰撞事故。對于危險(xiǎn)船舶指標(biāo),研究區(qū)內(nèi)大型油輪和運(yùn)輸危險(xiǎn)物品的貨船整體集中在“西南-東北”的近岸航道上,其中,風(fēng)險(xiǎn)較高的區(qū)域呈散點(diǎn)狀分布在航道及航道以南。研究區(qū)內(nèi)的助航標(biāo)志主要沿航道走向分布,也會布設(shè)在入港處或者航道重要轉(zhuǎn)向點(diǎn)處,供船舶定位、確定航向等,在該指標(biāo)上,東南部的航行環(huán)境風(fēng)險(xiǎn)較高。

    3.2 權(quán)重計(jì)算

    為了保證權(quán)重的客觀性、合理性,邀請三位海上交通安全領(lǐng)域相關(guān)的專家及教授,對各層級指標(biāo)進(jìn)行兩兩比較打分。在所有判斷矩陣滿足一致性條件(0.017 8~0.068 9)的情況下,利用模糊層次分析法計(jì)算得出指標(biāo)權(quán)重。在航行自然環(huán)境中,水深的權(quán)重最高,為0.119 0。大浪頻數(shù)和大風(fēng)頻數(shù)的權(quán)重分別為0.096 7和0.094 3。霧和降雨量對船舶航行的影響則相對較小,權(quán)重緊隨大浪頻數(shù)和大風(fēng)頻數(shù)之后,分別為0.089 1和0.074 9。在航行交通環(huán)境中,交通密度、助航標(biāo)志和船舶速度的重要程度相當(dāng),分別為0.144 5、0.143 6和0.143 0。危險(xiǎn)船舶的權(quán)重為0.094 9。

    3.3 航行環(huán)境風(fēng)險(xiǎn)分析

    根據(jù)線性加權(quán)法分別得到自然環(huán)境、交通環(huán)境影響下的航行環(huán)境風(fēng)險(xiǎn)等級分布和綜合風(fēng)險(xiǎn)等級分布,如圖4、圖5所示。其中,受水深限制,大型船舶在水深不足9 m的海域航行時擱淺風(fēng)險(xiǎn)較大[25],通過疊置分析將水深小于9 m的區(qū)域劃定為危險(xiǎn)水深區(qū)(圖5灰色區(qū)域)。對不同風(fēng)險(xiǎn)等級的分布面積進(jìn)行占比統(tǒng)計(jì),如圖6所示。

    1)自然環(huán)境影響:由圖4(a)和圖6(a)可知,在自然環(huán)境的影響下,研究區(qū)的航行環(huán)境風(fēng)險(xiǎn)等級整體呈“U”型分布,“U”字拐點(diǎn)可延伸到南澎列島北側(cè)海域,且風(fēng)險(xiǎn)等級自“U”字拐點(diǎn)起沿“SW—NE”依次降低。風(fēng)險(xiǎn)等級≥7級的區(qū)域占總面積的13.27%,主要分布在研究區(qū)西南部、詔安灣、廈門灣和同安灣。其中,風(fēng)險(xiǎn)等級≥9級的區(qū)域僅占1.85%,為“南澎列島—兄弟嶼”連線以南和同安灣以西的海域,大風(fēng)、大浪頻數(shù)高和水深條件不佳是分別導(dǎo)致這兩個區(qū)域風(fēng)險(xiǎn)較高的主要原因。風(fēng)險(xiǎn)等級≤6級的區(qū)域分布范圍最廣,共占86.73%,主要分布在研究區(qū)中部海域。

    2)海上交通環(huán)境影響:由圖4(b)和圖6(b)可知,在交通環(huán)境的影響下,研究區(qū)航行環(huán)境風(fēng)險(xiǎn)等級整體≤6級。風(fēng)險(xiǎn)等級處于1~5級的區(qū)域共占93.34%,主要分布在研究區(qū)各個港灣和航道上。風(fēng)險(xiǎn)等級處于6~8級的區(qū)域主要受交通密度影響,分布在主航道、廈門灣、圍頭灣海域,占6.54%。兄弟嶼東部海域的風(fēng)險(xiǎn)等級整體處于5~7級,風(fēng)險(xiǎn)等級整體由兄弟嶼東部最遠(yuǎn)處向海岸線方向遞減,這種遞減的趨勢主要是由于助航標(biāo)志的分布趨于密集,且船舶速度較低。風(fēng)險(xiǎn)等級≥8級的區(qū)域面積極小,僅占0.38%,主要呈散點(diǎn)狀分布在距海岸線較遠(yuǎn)海域,尤其是兄弟嶼東部海域,危險(xiǎn)船舶占比較多且助航標(biāo)志設(shè)立不完善是導(dǎo)致該區(qū)域風(fēng)險(xiǎn)較大的主要原因。

    3)綜合兩方面影響:由圖5和圖6(c)分析得,在自然環(huán)境和交通環(huán)境的綜合影響下,研究區(qū)的航行環(huán)境風(fēng)險(xiǎn)整體較低,且風(fēng)險(xiǎn)等級的空間分布格局大致呈現(xiàn)“SW—NE”逐漸降低?!澳吓炝袓u—兄弟嶼”連線的延長線以南和廈門灣海域航行環(huán)境風(fēng)險(xiǎn)等級偏高,達(dá)到7級及以上,面積占總研究區(qū)的6.17%。造成該區(qū)域風(fēng)險(xiǎn)等級較高的主要原因分別是大風(fēng)、大浪頻數(shù)較高和交通密度較大。其中,風(fēng)險(xiǎn)等級為8~10級的海域占比僅為1.62%,主要分布在研究區(qū)南部。但受交通環(huán)境影響,研究區(qū)內(nèi)的交通情況整體較好且助航標(biāo)志配備較為完善,在“南澎列島—兄弟嶼”連線的延長線以北區(qū)域的航行環(huán)境風(fēng)險(xiǎn)等級整體降低,大多處于5級及以下,面積占比共計(jì)73.7%。其中,東山灣、浮頭灣處風(fēng)險(xiǎn)等級較低,而詔安灣、廈門灣、圍頭灣和同安灣周邊風(fēng)險(xiǎn)等級相對略高。值得注意的是,危險(xiǎn)水深區(qū)水深不足9 m,主要分布在近岸和島嶼附近海域,占總面積的12.17%,船舶在此區(qū)域內(nèi)航行時擱淺風(fēng)險(xiǎn)較大。

    4 結(jié)論

    本文將地理空間技術(shù)與模糊層次分析相結(jié)合,從航行自然環(huán)境和海上交通環(huán)境兩方面選取大風(fēng)頻數(shù)、大浪頻數(shù)、水深、船舶速度、交通密度等9個指標(biāo)建立船舶航行環(huán)境風(fēng)險(xiǎn)評估體系,對廈門港及其周邊海域的航行環(huán)境進(jìn)行風(fēng)險(xiǎn)評估,形成了綜合航行自然環(huán)境和海上交通環(huán)境兩方面的航行環(huán)境風(fēng)險(xiǎn)分級圖,并分析了其空間分布特征。研究結(jié)果表明:

    1)在自然環(huán)境影響下,廈門港及其周邊海域的航行環(huán)境風(fēng)險(xiǎn)等級整體呈“U”型分布,且由西南向東北依次降低。風(fēng)險(xiǎn)等級較高的區(qū)域主要位于 “南澎列島—兄弟嶼”連線以南和同安灣以西。

    2)在海上交通環(huán)境的影響下,廈門港及其周邊海域的航行環(huán)境風(fēng)險(xiǎn)等級整體≤6級,僅有0.38%的區(qū)域風(fēng)險(xiǎn)等級≥8級,主要呈散點(diǎn)狀分布在距海岸線較遠(yuǎn)區(qū)域,尤其是兄弟嶼東部海域。

    3)綜合自然環(huán)境和交通環(huán)境兩方面的影響,廈門港及其周邊海域的航行環(huán)境整體處于中等偏低風(fēng)險(xiǎn),“南澎列島—兄弟嶼”連線的延長線以南和廈門灣海域的航行環(huán)境風(fēng)險(xiǎn)呈現(xiàn)較高趨勢。風(fēng)險(xiǎn)等級的空間分布格局大致呈“SW—NE”逐漸降低。其中,水深不足9 m的危險(xiǎn)水深區(qū)擱淺風(fēng)險(xiǎn)較大,主要分布在近岸和島嶼附近海域。經(jīng)統(tǒng)計(jì)得,研究區(qū)內(nèi)航行環(huán)境風(fēng)險(xiǎn)等級≥7級的區(qū)域占總面積的6.17%,為偏高風(fēng)險(xiǎn)區(qū)。對于本研究識別出的高風(fēng)險(xiǎn)區(qū)域,建議完善引航標(biāo)志的配置,加強(qiáng)日常巡邏措施,從而保障該區(qū)域的船舶航行安全。

    本文著重對航行的環(huán)境風(fēng)險(xiǎn)進(jìn)行評估,并獲取了風(fēng)險(xiǎn)等級空間分布特征。在指標(biāo)中尚未考慮海域沿岸國家相關(guān)政策、人為因素等軟性指標(biāo),未來將在指標(biāo)體系上進(jìn)一步優(yōu)化評估模型。

    參考文獻(xiàn):

    [1]GERLITZ L M, MEYER C. Small and medium-sized ports in the TEN-T network and nexus of europe's twin transition: the way towards sustainable and digital port service ecosystems[J]. Sustainability, 2021, 13(8): 4386.

    [2] SABINE K, CHRISTIAAN H J. Evaluation of total risk exposure and insurance premiums in the maritime industry[J]. Transportation Research Part D, 2017, 54: 321-334.

    [3] LUO M, SHIN S H. Half-century research developments in maritime accidents: Future directions[J]. Accid Anal Prev, 2016, 123: 448-460.

    [4] CARMEN G, LEONARD D. Navgation safety assessment based on short term seakeeping limit criteria for a feeder type container ship of 1300 TEU[J]. Galati University Press, 2015, 38: 5-8.

    [5] 張新放, 關(guān)克平. 船舶航行安全的影響因素及安全性評估研究[J]. 中國水運(yùn), 2015, 15(11): 50-52.

    [6] 張青, 李建, 顧婷婷, 等. 江海聯(lián)運(yùn)船舶航行安全氣象評價指標(biāo)體系的建立與應(yīng)用[J]. 中國水運(yùn), 2020(8): 27-29.

    [7] 李振福, 史曉梅, 周玉濤. 北極航行外部環(huán)境風(fēng)險(xiǎn)評估[J]. 海洋經(jīng)濟(jì), 2022, 12(1): 11-20.

    [8] 馬會, 吳兆麟. 港口航道水域操船環(huán)境危險(xiǎn)度的綜合評價[J]. 大連海事大學(xué)學(xué)報(bào), 1998(3): 17-20.

    [9] 陳堯, 肖英杰. 基于模糊綜合評判方法的仰光河航道安全評估[J]. 上海海事大學(xué)學(xué)報(bào), 2007(1): 50-56.

    [10]王晨, 江福才, 馬全黨. 基于熵權(quán)TOPSIS模型的航道引航環(huán)境風(fēng)險(xiǎn)評價[J]. 安全與環(huán)境學(xué)報(bào), 2016, 16(3): 33-37.

    [11]李振福, 陳霄, 趙梓州. 基于AIS數(shù)據(jù)的“21世紀(jì)海上絲綢之路”航行安全風(fēng)險(xiǎn)評估[J]. 水利水電技術(shù), 2021, 52(5): 1-12.

    [12]甄榮, RIVEIRO M, 金永興. 一種基于DBSCAN的船舶會遇實(shí)時識別方法[J]. 上海海事大學(xué)學(xué)報(bào), 2018, 39(1): 1-5.

    [13]任偉. 廈門港加快融入“一帶一路”的路徑探究[J]. 中國集體經(jīng)濟(jì), 2022(1): 28-29.

    [14]李振福, 史曉梅, 周玉濤. 北極航行外部環(huán)境風(fēng)險(xiǎn)評估[J]. 海洋經(jīng)濟(jì), 2022, 12(4): 93-102.

    [15]陳毅. 基于IWRAP和PAWSA模型的廈門港通航環(huán)境分析[D]. 廈門: 集美大學(xué), 2019.

    [16]郭昊. 基于AIS數(shù)據(jù)的軌跡段聚類方法研究[D]. 南京: 南京信息工程大學(xué), 2021.

    [17]汪恒, 蘭培真. 區(qū)間直覺模糊集的港口水域船舶航行環(huán)境風(fēng)險(xiǎn)評價[J]. 中國航海, 2021, 44(2): 38-44,52.

    [18]MUHAMMAD A H, STUART P, CHRIS R, et al. Assessing tropical cyclone risks using geospatial techniques[J]. Applied Geography, 2018, 98: 22-23.

    [19]LI S P, GONG Q X, YANG S L. A sustainable, regional agricultural development measurement system based on dissipative structure theory and the entropy weight method: a case study in Chengdu, China[J]. Sustainability, 2019, 11(19): 5313.

    [20]FAWAZ A, AHMED A, ASIF I K, et al. Evaluating security of big data through fuzzy based decision-making technique[J]. Computer Systems Science and Engineering, 2023, 44(1): 859-872.

    [21]曹馨勻. 基于三角模糊層次分析法的重慶地區(qū)建筑低碳化評價指標(biāo)體系研究[D]. 重慶: 重慶大學(xué), 2014.

    [22]CHANG D Y. Applications of the extent analysis method on fuzzy AHP[J]. European Journal of Operational Research, 1996, 95(3): 649-655.

    [23]BHADRAN A, GIRISHBAI D, JESIYA N P, et al. A GIS based Fuzzy-AHP for delineating groundwater potential zones in tropical river basin, southern part of India[J]. Geosystems and Geoenvironment, 2022, 1(4): 10093.

    [24]李志強(qiáng), 徐斌, 王升, 等. 基于改進(jìn)的AHP-SPA的邊坡穩(wěn)定性評價及工程應(yīng)用[J]. 公路交通科技, 2022, 39(1): 56-64.

    [25]HUANG J C, NIEH C Y, KUO H C. Risk assessment of ships maneuvering in an approaching channel based on AIS data[J]. Ocean Engineering, 2019, 173: 399-414.

    (責(zé)任編輯:于慧梅)

    Abstract: As the maritime traffic gets increasingly busy, ensuring the safe navigation of ships has become an important global issue. A risk assessment model of navigation environment from natural and maritime traffic is proposed in this paper. To evaluate the navigation risk of Xiamen Port and its surrounding areas, this study combines geographic space technology and fuzzy analytic hierarchy process (FAHP), forms the risk classification map and analyzes its spatial distribution characteristics. The results show that: Affected by the natural environment, the navigational environmental risks in this area are roughly distributed in a "U" shape, and the risk level gradually decreases from southwest to northeast. Under the influence of the traffic environment, the areas with risk greater than 8 level are scattered in the seas far from the coastline, especially in the eastern sea area of Brother Island. Taking the above two aspects into account, the overall risk of the navigation environment in Xiamen Port and its surrounding sea areas is moderately low. However, nearly 6.17% of the sea areas of Xiamen Port and its surrounding areas have the navigation risk greater than level 7, which is mainly distributed in the south of the extension of the “Nan Peng Islands-Brother Island” line and the Xiamen Bay waters.

    Key words: navigation risk; geographic space technology; fuzzy analytic hierarchy process; spatial distribution of risk

    基金項(xiàng)目:國家自然科學(xué)基金資助項(xiàng)目(42001401)

    作者簡介:杜思雨(1997—),女,在讀碩士,研究方向:GIS應(yīng)用、海上大數(shù)據(jù)挖掘分析,E-mail:631519737@qq.com.

    *通訊作者:朱大明,E-mail:634617255@qq.com.

    猜你喜歡
    廈門港航行海域
    遺落海域
    中國寶玉石(2021年5期)2021-11-18 07:42:32
    到慧骃國的航行
    埕島海域海上獨(dú)立樁拆除方案探討
    廈門港航道通過能力建模及應(yīng)用
    中國航海(2019年2期)2019-07-24 08:26:46
    小舟在河上航行
    自貿(mào)區(qū)背景下廈門港競爭力分析及轉(zhuǎn)型研究
    物流科技(2017年5期)2017-07-06 09:04:43
    廣東省海域使用統(tǒng)計(jì)分析
    航行
    青年歌聲(2017年6期)2017-03-13 00:57:56
    廈門港同時靠泊兩艘20萬噸級集箱船
    海域評估技術(shù)指引
    97在线人人人人妻| 又大又爽又粗| 国产区一区二久久| 国产成人精品久久二区二区免费| 日韩中文字幕欧美一区二区| 亚洲,欧美精品.| 久久人妻福利社区极品人妻图片| 久久久久久亚洲精品国产蜜桃av| 久久久久久久精品精品| 一区二区三区激情视频| 精品一品国产午夜福利视频| 亚洲七黄色美女视频| 少妇人妻久久综合中文| 黄色毛片三级朝国网站| 少妇粗大呻吟视频| 亚洲久久久国产精品| 亚洲av国产av综合av卡| 中文字幕最新亚洲高清| 日韩中文字幕欧美一区二区| 亚洲久久久国产精品| 免费高清在线观看视频在线观看| 搡老乐熟女国产| 欧美激情久久久久久爽电影 | 极品少妇高潮喷水抽搐| 黄色视频,在线免费观看| 麻豆国产av国片精品| 精品第一国产精品| 日韩中文字幕欧美一区二区| 十八禁网站网址无遮挡| 久久香蕉激情| 脱女人内裤的视频| 在线亚洲精品国产二区图片欧美| 久久久国产精品麻豆| 麻豆国产av国片精品| 50天的宝宝边吃奶边哭怎么回事| 老汉色∧v一级毛片| 亚洲av美国av| 亚洲一卡2卡3卡4卡5卡精品中文| 丰满迷人的少妇在线观看| 少妇 在线观看| 日韩欧美一区视频在线观看| 午夜福利在线观看吧| 国产亚洲欧美精品永久| 一级a爱视频在线免费观看| 啦啦啦免费观看视频1| 欧美日韩亚洲综合一区二区三区_| 国产麻豆69| 日日夜夜操网爽| 蜜桃国产av成人99| 亚洲精品美女久久av网站| 两性夫妻黄色片| 成人18禁高潮啪啪吃奶动态图| 久久久久国产精品人妻一区二区| 美女高潮喷水抽搐中文字幕| 最黄视频免费看| 精品国产乱子伦一区二区三区 | 国产麻豆69| 亚洲中文日韩欧美视频| 无遮挡黄片免费观看| 欧美日韩亚洲综合一区二区三区_| 日韩一卡2卡3卡4卡2021年| 亚洲成人手机| 少妇精品久久久久久久| 99国产精品一区二区蜜桃av | 少妇人妻久久综合中文| 老司机亚洲免费影院| 天天影视国产精品| 男女下面插进去视频免费观看| 国产男人的电影天堂91| 秋霞在线观看毛片| 在线天堂中文资源库| 久久香蕉激情| av电影中文网址| 少妇精品久久久久久久| 亚洲激情五月婷婷啪啪| 久久久久久久久免费视频了| www.自偷自拍.com| 肉色欧美久久久久久久蜜桃| 亚洲欧洲日产国产| 一区二区三区激情视频| 午夜福利免费观看在线| 亚洲国产欧美在线一区| 国产亚洲av高清不卡| 久久中文看片网| 免费人妻精品一区二区三区视频| 国精品久久久久久国模美| 好男人电影高清在线观看| 日韩熟女老妇一区二区性免费视频| 国产极品粉嫩免费观看在线| 老司机深夜福利视频在线观看 | 两个人看的免费小视频| 成年av动漫网址| 国产一区二区三区在线臀色熟女 | 少妇被粗大的猛进出69影院| 国产欧美日韩一区二区精品| 丰满少妇做爰视频| 欧美97在线视频| 亚洲精华国产精华精| av天堂久久9| 另类精品久久| 精品久久久久久久毛片微露脸 | 老熟妇仑乱视频hdxx| 欧美变态另类bdsm刘玥| 最近中文字幕2019免费版| 亚洲天堂av无毛| 亚洲精品国产色婷婷电影| 欧美性长视频在线观看| 97在线人人人人妻| 亚洲中文av在线| 国产黄色免费在线视频| 免费久久久久久久精品成人欧美视频| 不卡av一区二区三区| 18禁观看日本| 欧美午夜高清在线| 精品国产乱码久久久久久男人| 欧美成狂野欧美在线观看| 成人三级做爰电影| 日韩中文字幕欧美一区二区| 色94色欧美一区二区| 久久免费观看电影| svipshipincom国产片| 秋霞在线观看毛片| 午夜福利免费观看在线| 午夜久久久在线观看| 少妇 在线观看| 美女高潮喷水抽搐中文字幕| 精品国产一区二区久久| 男人爽女人下面视频在线观看| 欧美黑人精品巨大| 在线观看免费午夜福利视频| 人妻人人澡人人爽人人| 国产免费av片在线观看野外av| 久久性视频一级片| 国产成+人综合+亚洲专区| 青春草视频在线免费观看| 免费看十八禁软件| 亚洲中文日韩欧美视频| 欧美变态另类bdsm刘玥| 精品国产一区二区久久| a 毛片基地| 欧美黄色片欧美黄色片| av片东京热男人的天堂| 午夜两性在线视频| 亚洲国产日韩一区二区| 色综合欧美亚洲国产小说| 国产精品1区2区在线观看. | 午夜视频精品福利| 久久人妻熟女aⅴ| 亚洲欧洲日产国产| 12—13女人毛片做爰片一| 最近中文字幕2019免费版| 2018国产大陆天天弄谢| 另类精品久久| 久久99一区二区三区| 精品一区在线观看国产| 十八禁高潮呻吟视频| 男女下面插进去视频免费观看| 久久影院123| 12—13女人毛片做爰片一| 99国产综合亚洲精品| 亚洲欧洲日产国产| 99国产极品粉嫩在线观看| 久久99一区二区三区| 精品第一国产精品| 免费少妇av软件| 视频在线观看一区二区三区| 亚洲精品av麻豆狂野| 亚洲成人免费av在线播放| 老鸭窝网址在线观看| 人人妻人人添人人爽欧美一区卜| 丝袜喷水一区| 午夜精品久久久久久毛片777| 99久久99久久久精品蜜桃| 久久亚洲精品不卡| 日本wwww免费看| 欧美性长视频在线观看| 大香蕉久久网| 成人三级做爰电影| 悠悠久久av| 肉色欧美久久久久久久蜜桃| 91老司机精品| 国产淫语在线视频| 亚洲天堂av无毛| 无限看片的www在线观看| 久久久欧美国产精品| 免费在线观看黄色视频的| 中文字幕精品免费在线观看视频| 亚洲男人天堂网一区| 老司机午夜十八禁免费视频| 97精品久久久久久久久久精品| 黄色视频不卡| av天堂久久9| 日本五十路高清| 中亚洲国语对白在线视频| 免费一级毛片在线播放高清视频 | 久久免费观看电影| 午夜福利视频精品| 欧美激情 高清一区二区三区| 一区二区av电影网| 午夜两性在线视频| 久久99一区二区三区| 日本av手机在线免费观看| 新久久久久国产一级毛片| 亚洲精品自拍成人| 男男h啪啪无遮挡| 日韩欧美一区视频在线观看| 成人国产av品久久久| 97在线人人人人妻| 亚洲天堂av无毛| 欧美精品啪啪一区二区三区 | 国产男女内射视频| 精品人妻1区二区| 亚洲国产欧美一区二区综合| av电影中文网址| 午夜福利,免费看| 在线av久久热| 国产欧美日韩一区二区三区在线| 麻豆av在线久日| 国产成人a∨麻豆精品| 久久亚洲精品不卡| 一本色道久久久久久精品综合| 秋霞在线观看毛片| 免费女性裸体啪啪无遮挡网站| 青春草亚洲视频在线观看| 日韩欧美国产一区二区入口| 久久性视频一级片| 久久久久国产精品人妻一区二区| 男女之事视频高清在线观看| 亚洲av片天天在线观看| 91精品三级在线观看| 国产男人的电影天堂91| 极品少妇高潮喷水抽搐| 99热国产这里只有精品6| 欧美日韩精品网址| 久久人妻福利社区极品人妻图片| 日日摸夜夜添夜夜添小说| 久久久精品免费免费高清| 精品国产一区二区久久| 老司机靠b影院| 天天操日日干夜夜撸| 中文字幕高清在线视频| 国产一区二区三区综合在线观看| 成人影院久久| 涩涩av久久男人的天堂| 男人操女人黄网站| 黑人巨大精品欧美一区二区mp4| 大陆偷拍与自拍| 久久人人爽av亚洲精品天堂| 久热这里只有精品99| 性少妇av在线| 狠狠精品人妻久久久久久综合| av超薄肉色丝袜交足视频| 黑人操中国人逼视频| 欧美老熟妇乱子伦牲交| 欧美另类亚洲清纯唯美| 亚洲欧美色中文字幕在线| 久9热在线精品视频| 岛国在线观看网站| 日韩欧美免费精品| 丝袜美足系列| 黄色视频不卡| 亚洲精品久久午夜乱码| 蜜桃在线观看..| 国产精品亚洲av一区麻豆| 另类精品久久| 大香蕉久久成人网| 久热爱精品视频在线9| 美女脱内裤让男人舔精品视频| 天堂8中文在线网| 美女中出高潮动态图| 18禁黄网站禁片午夜丰满| 97精品久久久久久久久久精品| 中文字幕另类日韩欧美亚洲嫩草| 精品国产国语对白av| 国产精品.久久久| 菩萨蛮人人尽说江南好唐韦庄| 久热爱精品视频在线9| 欧美日韩亚洲综合一区二区三区_| 精品少妇久久久久久888优播| 日本五十路高清| 91av网站免费观看| 国产免费视频播放在线视频| a级片在线免费高清观看视频| 十八禁人妻一区二区| 男男h啪啪无遮挡| 久久av网站| av网站在线播放免费| 亚洲av成人不卡在线观看播放网 | 乱人伦中国视频| 亚洲五月婷婷丁香| 国产视频一区二区在线看| 亚洲av日韩在线播放| 精品卡一卡二卡四卡免费| 欧美黄色片欧美黄色片| 精品久久蜜臀av无| 精品亚洲成国产av| 亚洲欧美成人综合另类久久久| 精品国产超薄肉色丝袜足j| 交换朋友夫妻互换小说| av不卡在线播放| av有码第一页| 看免费av毛片| 亚洲成人国产一区在线观看| 国产高清国产精品国产三级| 另类亚洲欧美激情| 国产精品 欧美亚洲| 天天影视国产精品| av又黄又爽大尺度在线免费看| av线在线观看网站| 国产欧美日韩综合在线一区二区| 悠悠久久av| 一区二区三区精品91| 97人妻天天添夜夜摸| 丝袜美腿诱惑在线| 91九色精品人成在线观看| 老汉色av国产亚洲站长工具| 日本av手机在线免费观看| 1024香蕉在线观看| 国产亚洲欧美精品永久| 美女中出高潮动态图| 国产成人免费观看mmmm| 人人妻,人人澡人人爽秒播| 国产欧美日韩一区二区精品| 老司机靠b影院| 一区二区三区精品91| 一二三四在线观看免费中文在| 满18在线观看网站| 美女国产高潮福利片在线看| 久久精品国产亚洲av高清一级| 午夜日韩欧美国产| 国产精品九九99| 精品第一国产精品| 高清视频免费观看一区二区| 正在播放国产对白刺激| 午夜福利乱码中文字幕| 免费av中文字幕在线| 国产又爽黄色视频| 欧美一级毛片孕妇| 日韩电影二区| 欧美另类亚洲清纯唯美| 国产主播在线观看一区二区| 黑人猛操日本美女一级片| 日本欧美视频一区| 亚洲第一欧美日韩一区二区三区 | 国产av一区二区精品久久| 高清av免费在线| 欧美精品高潮呻吟av久久| 亚洲av男天堂| 国产一区有黄有色的免费视频| 国产一区二区在线观看av| 久久人人97超碰香蕉20202| av福利片在线| 蜜桃国产av成人99| 午夜精品久久久久久毛片777| 超碰97精品在线观看| 最新的欧美精品一区二区| 九色亚洲精品在线播放| 亚洲黑人精品在线| 一级黄色大片毛片| 亚洲五月婷婷丁香| 老司机靠b影院| 亚洲熟女精品中文字幕| 久久久久久人人人人人| 中文字幕人妻熟女乱码| www.自偷自拍.com| 国产一区二区在线观看av| 日韩制服丝袜自拍偷拍| 亚洲国产av新网站| 亚洲成人免费av在线播放| 桃花免费在线播放| 国产精品久久久av美女十八| 国产精品亚洲av一区麻豆| 在线观看人妻少妇| 少妇的丰满在线观看| 一个人免费看片子| 亚洲综合色网址| 欧美av亚洲av综合av国产av| 亚洲国产欧美在线一区| 日韩欧美一区视频在线观看| 搡老乐熟女国产| 久久中文看片网| 两个人看的免费小视频| 欧美日本中文国产一区发布| 真人做人爱边吃奶动态| 国产精品av久久久久免费| 国产在视频线精品| 狠狠婷婷综合久久久久久88av| 一区二区三区乱码不卡18| 丝袜美腿诱惑在线| 欧美激情久久久久久爽电影 | 国产男女内射视频| 久久天堂一区二区三区四区| 我的亚洲天堂| 日韩熟女老妇一区二区性免费视频| 真人做人爱边吃奶动态| 欧美日韩亚洲综合一区二区三区_| av电影中文网址| 男女之事视频高清在线观看| 777久久人妻少妇嫩草av网站| cao死你这个sao货| 色老头精品视频在线观看| 天堂中文最新版在线下载| 91精品国产国语对白视频| 欧美黄色片欧美黄色片| 性少妇av在线| 99久久国产精品久久久| 亚洲伊人久久精品综合| 在线亚洲精品国产二区图片欧美| 欧美精品人与动牲交sv欧美| 亚洲欧美清纯卡通| 亚洲中文av在线| 国产高清国产精品国产三级| 欧美精品高潮呻吟av久久| 99热网站在线观看| 自线自在国产av| 免费高清在线观看视频在线观看| 久久久久网色| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品久久午夜乱码| kizo精华| 午夜久久久在线观看| 91成人精品电影| 男人操女人黄网站| 欧美精品亚洲一区二区| 国产免费现黄频在线看| 搡老乐熟女国产| 午夜福利视频在线观看免费| 国产色视频综合| 丝袜脚勾引网站| 精品久久久久久久毛片微露脸 | 欧美人与性动交α欧美软件| 国产一区二区三区综合在线观看| 欧美少妇被猛烈插入视频| 丝袜喷水一区| 一进一出抽搐动态| 日韩一区二区三区影片| 最近最新中文字幕大全免费视频| 热re99久久国产66热| 欧美老熟妇乱子伦牲交| 色综合欧美亚洲国产小说| 亚洲精品一区蜜桃| 亚洲色图 男人天堂 中文字幕| 一边摸一边抽搐一进一出视频| 中亚洲国语对白在线视频| 激情视频va一区二区三区| 欧美日韩福利视频一区二区| 精品人妻一区二区三区麻豆| 国产亚洲欧美精品永久| 亚洲男人天堂网一区| 啦啦啦视频在线资源免费观看| 亚洲人成77777在线视频| 亚洲va日本ⅴa欧美va伊人久久 | 女性被躁到高潮视频| 最近最新免费中文字幕在线| 久久精品熟女亚洲av麻豆精品| 中文字幕最新亚洲高清| 亚洲中文av在线| 中文字幕精品免费在线观看视频| 最近最新免费中文字幕在线| 12—13女人毛片做爰片一| 性少妇av在线| 精品人妻在线不人妻| 亚洲av电影在线观看一区二区三区| 亚洲精品粉嫩美女一区| 又黄又粗又硬又大视频| 亚洲一区二区三区欧美精品| 老司机靠b影院| 久久久久久久久久久久大奶| 欧美成人午夜精品| bbb黄色大片| 天天躁狠狠躁夜夜躁狠狠躁| 夫妻午夜视频| 久久久久精品人妻al黑| 高清欧美精品videossex| 亚洲精华国产精华精| 蜜桃国产av成人99| 十分钟在线观看高清视频www| 一区二区日韩欧美中文字幕| 天堂8中文在线网| 亚洲人成77777在线视频| 美女国产高潮福利片在线看| 亚洲熟女精品中文字幕| 无遮挡黄片免费观看| 久久久精品94久久精品| 日韩,欧美,国产一区二区三区| 丰满少妇做爰视频| 亚洲精品自拍成人| 自线自在国产av| 人妻久久中文字幕网| 亚洲,欧美精品.| 99久久综合免费| 满18在线观看网站| 蜜桃国产av成人99| 久久久久视频综合| 美女高潮喷水抽搐中文字幕| 女人爽到高潮嗷嗷叫在线视频| 人成视频在线观看免费观看| 欧美成人午夜精品| 日本撒尿小便嘘嘘汇集6| 亚洲第一av免费看| 午夜免费鲁丝| 国内毛片毛片毛片毛片毛片| 99热全是精品| 老司机午夜十八禁免费视频| 国产免费av片在线观看野外av| 日日爽夜夜爽网站| 永久免费av网站大全| 欧美日韩黄片免| 欧美黄色淫秽网站| 91麻豆精品激情在线观看国产 | 伊人亚洲综合成人网| 美女视频免费永久观看网站| 国产成人欧美在线观看 | 在线观看免费视频网站a站| 黄色a级毛片大全视频| 久久人妻福利社区极品人妻图片| 日韩一区二区三区影片| 国产极品粉嫩免费观看在线| avwww免费| 啦啦啦在线免费观看视频4| 99国产综合亚洲精品| 久久人人爽人人片av| 亚洲国产av影院在线观看| 久久精品国产亚洲av香蕉五月 | av天堂久久9| 亚洲人成77777在线视频| 亚洲欧美精品综合一区二区三区| 日本一区二区免费在线视频| 精品视频人人做人人爽| 久久 成人 亚洲| 亚洲黑人精品在线| 12—13女人毛片做爰片一| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av美国av| 99国产精品一区二区三区| 十八禁人妻一区二区| 新久久久久国产一级毛片| 丰满人妻熟妇乱又伦精品不卡| 99精品欧美一区二区三区四区| 中文字幕精品免费在线观看视频| 两性夫妻黄色片| 搡老岳熟女国产| 欧美大码av| 国产精品一区二区在线观看99| 久热爱精品视频在线9| 亚洲国产av新网站| 1024香蕉在线观看| 一本久久精品| 欧美精品一区二区免费开放| 日本欧美视频一区| 日日夜夜操网爽| 啪啪无遮挡十八禁网站| www.自偷自拍.com| 搡老乐熟女国产| 国产有黄有色有爽视频| 999久久久精品免费观看国产| 久久久久久亚洲精品国产蜜桃av| 人人妻人人爽人人添夜夜欢视频| 国产精品熟女久久久久浪| 99香蕉大伊视频| 精品一品国产午夜福利视频| 色播在线永久视频| 桃红色精品国产亚洲av| 国产一区有黄有色的免费视频| 久久久久久久精品精品| cao死你这个sao货| 自线自在国产av| 成人影院久久| 国产伦人伦偷精品视频| 超碰97精品在线观看| 成人三级做爰电影| 成人免费观看视频高清| 狂野欧美激情性bbbbbb| 亚洲av男天堂| 一二三四社区在线视频社区8| 不卡av一区二区三区| 精品熟女少妇八av免费久了| 考比视频在线观看| 美国免费a级毛片| 制服诱惑二区| 亚洲精品久久成人aⅴ小说| 亚洲av片天天在线观看| 久久性视频一级片| 美女中出高潮动态图| 午夜福利视频精品| 欧美亚洲 丝袜 人妻 在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲 欧美一区二区三区| 美女大奶头黄色视频| 亚洲少妇的诱惑av| 黄片播放在线免费| 天天操日日干夜夜撸| 国产在线一区二区三区精| 久久影院123| 在线看a的网站| 黑人巨大精品欧美一区二区蜜桃| 久久国产精品人妻蜜桃| 国产成人精品无人区| 乱人伦中国视频| 黄色a级毛片大全视频| av网站免费在线观看视频| 成人国产一区最新在线观看| 亚洲综合色网址| 国产伦理片在线播放av一区| 色播在线永久视频| 成年人黄色毛片网站| 久久亚洲精品不卡| 亚洲精品av麻豆狂野| 久久精品国产亚洲av高清一级| 国产精品一二三区在线看| 国产97色在线日韩免费| 中文字幕人妻丝袜一区二区| 午夜激情av网站| 首页视频小说图片口味搜索| 国产三级黄色录像|