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

    面向海運統(tǒng)計的AIS大數(shù)據(jù)挖掘分析研究

    2022-01-18 11:49:08趙龍飛姜曉軼呂憧憬
    海洋科學(xué) 2021年12期
    關(guān)鍵詞:航次貨品泊位

    趙龍飛, 姜曉軼, 孫 苗, 曹 磊, 郭 雪, 呂憧憬

    面向海運統(tǒng)計的AIS大數(shù)據(jù)挖掘分析研究

    趙龍飛1, 2, 姜曉軼1, 2, 孫 苗1, 2, 曹 磊1, 2, 郭 雪1, 2, 呂憧憬1, 2

    (1. 國家海洋信息中心, 天津 300171; 2. 自然資源部海洋信息技術(shù)創(chuàng)新中心, 天津 300171)

    為解決傳統(tǒng)海運產(chǎn)業(yè)統(tǒng)計方法數(shù)據(jù)質(zhì)量不高、時效性差、統(tǒng)計產(chǎn)品陳舊短缺、公信力不足等的問題, 提出基于AIS大數(shù)據(jù)挖掘分析開展海運統(tǒng)計的方法。詳細(xì)闡述了基于AIS大數(shù)據(jù)的海運統(tǒng)計分析技術(shù)路線、大數(shù)據(jù)平臺技術(shù)架構(gòu), 以及電子圍欄分析、航行事件分析、航次分析和統(tǒng)計指標(biāo)生成等關(guān)鍵模型算法。以2019年3月至5月全球大宗貨品船舶的AIS數(shù)據(jù)應(yīng)用為例表明, 該方法可提供港口、海上通道和大宗貨品三方面的海運大數(shù)據(jù)統(tǒng)計指標(biāo), 為實現(xiàn)海運即時化、準(zhǔn)確化、精細(xì)化的統(tǒng)計分析與展現(xiàn)提供了新思路。

    海運統(tǒng)計; 大數(shù)據(jù); 海洋經(jīng)濟; AIS; 挖掘分析

    海運業(yè)是國際貿(mào)易的橋梁和紐帶[1]。海運統(tǒng)計是摸清海運貿(mào)易量、提供海運產(chǎn)業(yè)運行情況的重要工具, 對于發(fā)現(xiàn)海運產(chǎn)業(yè)存在的問題, 科學(xué)制定海運產(chǎn)業(yè)政策, 實施正確、有力的海洋經(jīng)濟宏觀調(diào)控措施, 從而推進海洋強國和海運強國建設(shè)等具有重要意義[2-3]。當(dāng)前, 作為政府統(tǒng)計中的部門統(tǒng)計之一, 海運統(tǒng)計仍遵循制度設(shè)計、調(diào)查采集、數(shù)據(jù)整理上報、數(shù)據(jù)發(fā)布與開發(fā)利用的流程開展, 面臨著數(shù)據(jù)質(zhì)量不高、時效性差、統(tǒng)計產(chǎn)品陳舊短缺、數(shù)據(jù)公信力不足等的問題[4]。大數(shù)據(jù)統(tǒng)計通過對產(chǎn)業(yè)運行過程中產(chǎn)生的樣本海量、實時、非結(jié)構(gòu)化等傳統(tǒng)統(tǒng)計調(diào)查無法獲取的數(shù)據(jù)進行挖掘分析, 如電商交易數(shù)據(jù)、位置數(shù)據(jù)等, 為從根本上解決傳統(tǒng)統(tǒng)計存在問題打開了新視角[5]。

    近幾年基于船舶自動識別系統(tǒng)(automatic identifi-cation system, AIS)大數(shù)據(jù)進行船舶航行、航道和港口作業(yè)統(tǒng)計以及經(jīng)濟運行態(tài)勢分析成為新的趨勢[6-13]。何建春[14]基于船舶AIS、RFID、視頻監(jiān)控等多傳感融合的全樣本水運統(tǒng)計模型與方法, 對內(nèi)河航道船舶流量流向、貨物流量流向、船舶航行里程、航道航行密度進行了統(tǒng)計分析與展現(xiàn); 劉滿娜[15]建立了基于AIS的港口監(jiān)控與分析系統(tǒng), 實現(xiàn)船舶進出港、船舶錨泊和系泊等的分析統(tǒng)計; 祝天虎[16]以AIS和大數(shù)據(jù)挖掘技術(shù)為基礎(chǔ), 設(shè)計了基于Hadoop的聚類算法和船舶航次預(yù)測算法, 實現(xiàn)了對鐵礦石貿(mào)易量的預(yù)測。航運電商平臺如船訊網(wǎng)、Hifleet船隊在線等將海量零散、動態(tài)的AIS等海運數(shù)據(jù)搜集、整理、建模后形成多維基礎(chǔ)大數(shù)據(jù)模型, 并提供了大宗貨品海運統(tǒng)計、港口運營效能、船舶運力分布及預(yù)測、集裝箱準(zhǔn)班率等分析產(chǎn)品和服務(wù)。

    宏觀層面的海運統(tǒng)計包括港口、對外貿(mào)易和大宗貨品等方面, 基于AIS開展海運統(tǒng)計分析不僅需要全球范圍的港口、泊位數(shù)據(jù), 還要接入、處理并分析所有船舶的海量全球?qū)崟r和歷史AIS數(shù)據(jù), 并能據(jù)此挖掘出各類型貨物的運量情況[17]。但從筆者搜索到的業(yè)界及學(xué)術(shù)界的公開資料及文獻來看, 對基于大數(shù)據(jù)平臺架構(gòu)的AIS接入、存儲、處理分析技術(shù)和宏觀海運統(tǒng)計挖掘分析模型構(gòu)建的整體技術(shù)方案仍描述不足。因此, 本文面向港口、對外貿(mào)易及大宗貨品的海運大數(shù)據(jù)統(tǒng)計分析需求, 給出了基于AIS開展海運統(tǒng)計的技術(shù)方法, 詳細(xì)闡述了AIS大數(shù)據(jù)平臺的構(gòu)建方案, 以及電子圍欄分析、航行事件挖掘、航次分析等關(guān)鍵模型算法, 并以歷史全球大宗貨品船舶的AIS數(shù)據(jù)為例進行了實證分析, 為變革海運統(tǒng)計模式, 實現(xiàn)海運即時化、準(zhǔn)確化、精細(xì)化的大數(shù)據(jù)統(tǒng)計分析與展現(xiàn)奠定技術(shù)基礎(chǔ)。

    1 方法概述

    1.1 AIS數(shù)據(jù)

    本文提出的海運大數(shù)據(jù)統(tǒng)計分析方法主要基于AIS數(shù)據(jù)。根據(jù)IMO《國際海上人命安全公約》規(guī)定, 總噸大于300 t的海運船舶都需要安裝AIS終端設(shè)備。船舶發(fā)送的AIS信息通??梢员桓浇陌痘蛘咝l(wèi)星AIS接收裝置接收, 基本可實現(xiàn)全球范圍的船舶追蹤。AIS數(shù)據(jù)主要包括動態(tài)信息和靜態(tài)信息, 其中動態(tài)信息主要包括船舶位置(經(jīng)緯度坐標(biāo))、航速、航向等, 靜態(tài)信息包括船舶載重、吃水、目的港和預(yù)到時間等。動態(tài)信息一般由與AIS相連接的物聯(lián)網(wǎng)設(shè)備自動提供, 這些信息在航行時通常幾秒鐘更新一次, 在停泊時幾分鐘更新1次, 且準(zhǔn)確性較高。靜態(tài)信息一般由人工輸入, 存在信息缺失或者不正確問題, 可靠性相對較差。

    1.2 技術(shù)路線

    基于船舶AIS開展海運統(tǒng)計的數(shù)據(jù)模型涉及AIS動態(tài)感知信息以及船舶檔案、港口、泊位、海上通道等基礎(chǔ)數(shù)據(jù)以及這些信息在業(yè)務(wù)上、邏輯上、時空上的相互關(guān)系。大數(shù)據(jù)統(tǒng)計模型的建立有賴于船舶流量流向、貨物流量流向、船舶進出港、途經(jīng)海上通道等信息的生成, 其實現(xiàn)的整體技術(shù)流程如圖1所示。

    圖1 基于AIS大數(shù)據(jù)的海運統(tǒng)計技術(shù)路線

    1) 數(shù)據(jù)資料收集: 包括全球范圍內(nèi)實時及歷史衛(wèi)星AIS數(shù)據(jù)、船舶檔案數(shù)據(jù)、全球港口及泊位數(shù)據(jù)等, 來源包括商業(yè)采購、英國勞氏海事數(shù)據(jù)庫、人工調(diào)研補充等。

    2) 大數(shù)據(jù)平臺及分析模型構(gòu)建: 利用分布式存儲和云計算技術(shù), 搭建AIS大數(shù)據(jù)存儲和計算平臺, 研發(fā)并集成電子圍欄分析、航行事件判斷、航次生成等挖掘分析模型算法, 滿足AIS實時數(shù)據(jù)接入、存儲和處理分析的要求。

    3) 數(shù)據(jù)清洗及數(shù)據(jù)倉庫構(gòu)建: 基于大數(shù)據(jù)平臺, 開展AIS數(shù)據(jù)的抽取轉(zhuǎn)換和清洗等, 處理結(jié)果以船舶航跡點和航跡線的形式形成數(shù)據(jù)倉庫, 滿足后續(xù)挖掘分析的需要。

    4) 航行事件挖掘: 通過對AIS航跡數(shù)據(jù)的挖掘分析, 將船舶航行過程劃分為預(yù)裝航行中、準(zhǔn)備裝貨、正在裝貨、裝貨完成、卸貨航行中、準(zhǔn)備卸貨、正在卸貨、卸貨完成8種航行狀態(tài)。

    5) 航次分析: 分析船舶事件, 形成船舶航行從裝貨到卸貨的完整航次周期, 分析結(jié)果包含船舶航行航次的裝貨港口、卸貨港口、途經(jīng)海上通道、貨品類型和載貨量等信息。

    6) 統(tǒng)計指標(biāo)生成: 在船舶航次信息基礎(chǔ)上, 進行港口船舶進出港統(tǒng)計、海上通道船舶流量統(tǒng)計, 以及大宗貨品裝卸貨統(tǒng)計、進出口統(tǒng)計等, 并形成多維度的大數(shù)據(jù)統(tǒng)計指標(biāo)。

    7) 結(jié)果分析及可視化展示: 分析評估航次信息和成果指標(biāo)的準(zhǔn)確度, 迭代優(yōu)化挖掘分析模型, 同時開發(fā)前后臺功能, 實現(xiàn)基于GIS的分析結(jié)果關(guān)聯(lián)性交互查詢展示。

    2 大數(shù)據(jù)平臺架構(gòu)

    AIS大數(shù)據(jù)平臺的建設(shè)是實現(xiàn)海運大數(shù)據(jù)統(tǒng)計的技術(shù)基礎(chǔ), 其目標(biāo)是基于實時與歷史AIS數(shù)據(jù), 結(jié)合船舶檔案、港口資料、調(diào)研數(shù)據(jù)等挖掘知識庫, 輔以在線加工, 通過不斷的數(shù)據(jù)挖掘和經(jīng)驗積累, 形成特定的基礎(chǔ)數(shù)據(jù)庫和數(shù)據(jù)分析模型。最終自動化定時輸出港口、海上通道及貨品運量等統(tǒng)計指標(biāo), 以及各類商品、船舶相關(guān)的動態(tài)信息, 并支持開放式編輯運算、可視化交互展示和報表導(dǎo)出。

    2.1 平臺總體架構(gòu)

    秉承“高內(nèi)聚低耦合”的設(shè)計理念, AIS大數(shù)據(jù)平臺包括基礎(chǔ)設(shè)施層、數(shù)據(jù)資源層、數(shù)據(jù)分析層和服務(wù)封裝層四個部分, 如圖2所示。

    圖2 AIS大數(shù)據(jù)平臺總體架構(gòu)圖

    基礎(chǔ)設(shè)施層: 由網(wǎng)絡(luò)、主機、存儲、安全等設(shè)備組成, 為大數(shù)據(jù)平臺的存儲和計算提供基礎(chǔ)軟硬件支撐, 也為基于大數(shù)據(jù)平臺的各類應(yīng)用系統(tǒng)運行提供基礎(chǔ)支撐環(huán)境。數(shù)據(jù)資源層: 實現(xiàn)包括實時AIS數(shù)據(jù)和歷史AIS數(shù)據(jù)的采集、抽取轉(zhuǎn)換和清洗, 并基于大數(shù)據(jù)存儲架構(gòu)滿足各類基礎(chǔ)數(shù)據(jù)、過程數(shù)據(jù)和成果數(shù)據(jù)存儲管理、查詢檢索等的需求。數(shù)據(jù)分析層: 集成大數(shù)據(jù)分析模型算法, 基于數(shù)據(jù)存儲層中的各類數(shù)據(jù), 通過挖掘分析, 生成船舶航行事件、航次、航跡等, 并最終形成各類航運統(tǒng)計指標(biāo)。服務(wù)封裝層: 負(fù)責(zé)根據(jù)不同需求對各種大數(shù)據(jù)處理和分析功能進行封裝并對外提供服務(wù)。

    2.2 平臺數(shù)據(jù)架構(gòu)

    基于大數(shù)據(jù)平臺總體架構(gòu), 結(jié)合AIS海運統(tǒng)計技術(shù)路線, 實現(xiàn)對大數(shù)據(jù)平臺的數(shù)據(jù)流程設(shè)計, 如圖3所示, 包括數(shù)據(jù)源、數(shù)據(jù)接入處理、數(shù)據(jù)存儲管理、數(shù)據(jù)計算分析和數(shù)據(jù)展示應(yīng)用5個部分。

    1) 數(shù)據(jù)源: 包括實時AIS數(shù)據(jù)、歷史AIS數(shù)據(jù), 以及船舶檔案、港口、泊位、海上通道、電子海圖等基礎(chǔ)支撐數(shù)據(jù), 其中實時AIS數(shù)據(jù)是指實時產(chǎn)生并在較短時間間隔內(nèi)(一般不超過10 min)動態(tài)增加的流式AIS大數(shù)據(jù), 歷史AIS數(shù)據(jù)指過去某個時間段內(nèi)所有船舶發(fā)送的AIS數(shù)據(jù), 以文本文件的形式存儲, 數(shù)據(jù)格式為標(biāo)準(zhǔn)AIS原始報文。

    圖 3 AIS大數(shù)據(jù)平臺數(shù)據(jù)架構(gòu)圖

    2) 數(shù)據(jù)接入處理: 對實時AIS數(shù)據(jù), 首先搭建流數(shù)據(jù)處理框架, 通過TCP協(xié)議接入外部實時AIS數(shù)據(jù)源, 并統(tǒng)一交由Kafka消息中間件進行存儲(Kafka實現(xiàn)高吞吐的同時可以保證準(zhǔn)確的發(fā)送給Flink流數(shù)據(jù)處理引擎, 避免因Flink出現(xiàn)故障導(dǎo)致數(shù)據(jù)丟失的問題)并分發(fā)給Flink進行多機多線程處理。Flink流數(shù)據(jù)處理引擎首先需要完成AIS數(shù)據(jù)的清洗工作, 包括去除無效數(shù)據(jù)與重復(fù)數(shù)據(jù), 整合位置、吃水、目的地等, 并與上一條AIS數(shù)據(jù)對比補全航行相關(guān)數(shù)據(jù), 同時, Flink還將處理后的AIS數(shù)據(jù)存儲到AIS日志文件中, 并由Flume日志采集器定時獲取AIS日志文件并合并到AIS月日志文件中; 對歷史AIS數(shù)據(jù), 采用ETL工具Kettle執(zhí)行批量任務(wù)抽取報文文件, 并進行清洗、去重、去無效等批操作, 形成統(tǒng)一格式的待處理文件, 同樣將處理結(jié)果存儲到AIS日志文件中。

    3) 數(shù)據(jù)存儲管理: 對Flink處理后的實時AIS數(shù)據(jù), 一方面以“MMSI+日期”為Key存儲到Redis內(nèi)存數(shù)據(jù)庫中, 供Spark離線任務(wù)進行航行事件生成和進出海上通道判斷等實時性分析工作; 另一方面, 在空間內(nèi)存數(shù)據(jù)庫H2GIS中記錄船舶最近一次的位置, 滿足后期實時顯示船舶當(dāng)前位置的需求。對清洗后形成AIS日志文件, 通過Flume自動導(dǎo)入到Hive數(shù)據(jù)倉庫中, 方便后期問題回溯, 并為數(shù)據(jù)計算分析提供支撐?;贖Base存儲通過數(shù)據(jù)計算分析形成的航跡線、航行事件、航次信息和數(shù)據(jù)集市。船舶檔案、港口、泊位等基礎(chǔ)支撐數(shù)據(jù), 整體數(shù)據(jù)量不大, 且在進行數(shù)據(jù)挖掘分析時需頻繁調(diào)用, 因此基于關(guān)系型數(shù)據(jù)庫MySQL進行存儲。

    4) 數(shù)據(jù)計算分析: 首先利用Spark事件分析任務(wù)批量從Redis內(nèi)存數(shù)據(jù)庫及Hive中提取AIS航跡點數(shù)據(jù), 結(jié)合MySQL中的船舶檔案、港口、泊位等基礎(chǔ)支撐數(shù)據(jù), 進行分析形成航行事件數(shù)據(jù)并存儲到HBase中, 然后通過Spark航次分析任務(wù)從HBase中獲取待分析的航行事件并分析形成航次數(shù)據(jù)并存儲在HBase中。數(shù)據(jù)集市是最終呈現(xiàn)的統(tǒng)計指標(biāo)表與維度表, 包括港口統(tǒng)計、通道統(tǒng)計、貨品統(tǒng)計等方面, 其多維數(shù)據(jù)統(tǒng)計查詢通過分析型數(shù)據(jù)倉庫Kylin提供的接口實現(xiàn), Hive和HBase服務(wù)于Kylin, 為Kylin提供數(shù)據(jù), 并用HBase存放Kylin生成的Cube。利用Geomesa時空大數(shù)據(jù)處理引擎, 通過讀取存儲于HBase和Hive中的AIS數(shù)據(jù), 建立時空索引并基于Spark進行計算分析, 進而實現(xiàn)單船軌跡回放、多船航跡熱力圖、航跡線擬合等空間大數(shù)據(jù)分析, 并具備OGC標(biāo)準(zhǔn)服務(wù)接口數(shù)據(jù)訪問能力, 滿足后期數(shù)據(jù)展示應(yīng)用的需求。

    5) 數(shù)據(jù)展示應(yīng)用: 以數(shù)據(jù)存儲和分析計算服務(wù)為支撐, 實現(xiàn)分析結(jié)果的統(tǒng)計圖表、統(tǒng)計報表和地圖可視化展示。本文基于UReport2報表引擎和ECharts可視化庫, 讀取存儲在HBase中的航次數(shù)據(jù)和統(tǒng)計指標(biāo)數(shù)據(jù), 實現(xiàn)報表化服務(wù)和圖表可視化服務(wù)。以H2GIS中的實時AIS點位數(shù)據(jù)及HBase中的航跡線數(shù)據(jù)為基礎(chǔ), 結(jié)合Geomesa的空間分析服務(wù), 利用Geoserver地圖服務(wù)器發(fā)布形成標(biāo)準(zhǔn)地圖服務(wù), 基于OpenLayers地圖庫, 實現(xiàn)海量船舶AIS航跡點及航跡線的地圖可視分析服務(wù)。

    3 關(guān)鍵模型算法

    根據(jù)港口、泊位等信息, 對海量AIS軌跡數(shù)據(jù)進行分割, 解析出船舶的所有航次, 是實現(xiàn)海運統(tǒng)計的關(guān)鍵。基于AIS軌跡數(shù)據(jù)對船舶航行中所在位置的航行狀態(tài)判斷, 即航行事件分析, 是航次生成的基礎(chǔ)。而航行事件分析有賴于AIS軌跡點與港口及泊位的位置關(guān)系分析, 即需通過電子圍欄技術(shù)實現(xiàn)。

    3.1 電子圍欄分析

    電子圍欄分析的核心內(nèi)容包括港口、泊位的電子圍欄構(gòu)建以及與船舶AIS軌跡點的相交計算。電子圍欄構(gòu)建是基于港口和泊位數(shù)據(jù)中的點位位置信息, 生成其多邊形范圍, 包括計算生成和人工調(diào)研繪制兩種方法, 計算生成是以港口和泊位的點位位置為圓心, 以特定值或以港口(泊位)與周邊港口(泊位)距離的最小值為半徑畫圓; 人工調(diào)研繪制是對重點區(qū)域港口和泊位通過實地調(diào)研或基于影像圖的方法手動繪制電子圍欄, 以獲得更為準(zhǔn)確的港口和泊位的空間范圍。在構(gòu)建的電子圍欄基礎(chǔ)上, 通過空間相交分析計算, 判斷AIS軌跡點是否在港口或泊位中, 詳細(xì)方法在此不再贅述。

    3.2 航行事件分析

    航行事件分析的主要流程是: 首先判斷船舶AIS軌跡點數(shù)據(jù)的航行狀態(tài)和航速, 然后與港口、泊位的電子圍欄進行相交計算, 最后根據(jù)港口、泊位的裝卸貨標(biāo)識及船舶吃水值對船舶裝卸貨狀態(tài)及貨品類型進行判斷, 形成包含8種航行狀態(tài)的船舶航行事件并按時間升序排列。

    3.3 航次分析

    航次分析通過分析按時間升序排列的航行事件, 形成航行周期并對應(yīng)到不同的運輸航次。需要注意的是, 在一個航次周期中, 可能出現(xiàn)“多次裝貨”及“多次卸貨”的情況, 需分別根據(jù)“裝貨停止”和“卸貨停止”后下一事件的狀態(tài)進行判別, 若“裝貨停止”事件后再次出現(xiàn)“預(yù)裝航行中”事件, 則該航次為“多次裝貨”; 若“卸貨停止”事件后再次出現(xiàn)“卸貨航行中”事件, 則該航次為“多次卸貨”, 詳細(xì)算法流程見圖4。

    圖 4 航次分析算法流程圖

    航次中的裝貨港口和卸貨港口可根據(jù)航行事件判斷中經(jīng)電子圍欄分析得出的泊位的所屬港口確定; 途經(jīng)海上通道通過航次航跡線與海上通道進行相交計算分析判斷得到; 各航次的載貨類型可以依據(jù)裝卸貨泊位或港口的貨品類型屬性確定, 對應(yīng)航次的載貨量根據(jù)航次裝貨和卸貨后的吃水差進行估算得出。

    3.4 統(tǒng)計指標(biāo)生成

    在航次分析形成的航次數(shù)據(jù)基礎(chǔ)上, 生成包含港口、海上通道、大宗貨品三方面的統(tǒng)計指標(biāo), 具體流程如下:

    1) 港口統(tǒng)計

    (1) 篩選出航次數(shù)據(jù)的裝貨港或卸貨港為需統(tǒng)計港口的航次;

    (2) 篩選出進離港時間在預(yù)先設(shè)定的統(tǒng)計周期(如月、季、半年、年等)的航次;

    (3) 根據(jù)航次數(shù)據(jù)的裝貨港、卸貨港、途經(jīng)海上通道、載貨類型、載貨量等信息, 并結(jié)合船舶檔案數(shù)據(jù), 形成包含進離港船舶累計載重及艘次、大宗貨品運量、途經(jīng)海上通道、進口國/出口國等多維度的港口統(tǒng)計指標(biāo)。

    2) 海上通道統(tǒng)計

    (1) 篩選出航次數(shù)據(jù)途經(jīng)海上通道包含需統(tǒng)計海上通道的航次;

    (2) 篩選出途經(jīng)該海上通道的時間在預(yù)先設(shè)定的統(tǒng)計周期(如月、季、半年、年等)的航次;

    (3) 根據(jù)航次數(shù)據(jù)的裝貨港及國家、卸貨港及國家、載貨類型、載貨量等信息, 并結(jié)合船舶檔案數(shù)據(jù), 形成包含進出海上通道船舶累計載重及艘次、大宗貨品運量、各港口/國家途經(jīng)該通道船舶累計載重及艘次等多維度的港口統(tǒng)計指標(biāo)。

    3) 大宗貨品統(tǒng)計

    (1) 對某一港口或海上通道的貨品統(tǒng)計, 可按港口統(tǒng)計或海上通道統(tǒng)計方法, 形成包括運輸大宗貨品船舶的累計載重及艘次、貨運量、來源/去向港口及國家等多維統(tǒng)計指標(biāo);

    (2) 對多個港口乃至某個國家的貨品統(tǒng)計, 則需在港口統(tǒng)計的基礎(chǔ)上進行匯總統(tǒng)計, 形成包括運輸船舶的累計載重及艘次、貨運量、進口/出口國、途經(jīng)海上通道等多維統(tǒng)計指標(biāo)。

    4 實驗與分析

    4.1 實驗數(shù)據(jù)

    考慮到承擔(dān)大宗貨品海運的船舶船型、裝卸貨港口、泊位都較為固定, 與本文提出的航行事件及航次分析算法適用場景一致, 因此, 以全球大宗貨品船舶的AIS數(shù)據(jù)挖掘分析為例, 驗證本文方法的可行性。如表1及圖5所示, 本文實驗采用的數(shù)據(jù)包括承擔(dān)大宗貨品運輸?shù)娜? 005個港口, 41 222個泊位, 25 497艘船舶檔案, 以及這些船舶在2019年3月至5月產(chǎn)生的衛(wèi)星AIS數(shù)據(jù)(來源于美國的ORBCOMM公司), 約32億條。

    表 1 實驗數(shù)據(jù)介紹

    圖 5 實驗數(shù)據(jù)地圖可視化

    4.2 實驗結(jié)果與分析

    按本文所述方法并采用6臺服務(wù)器作為計算節(jié)點搭建了大數(shù)據(jù)平臺, 其中, 主節(jié)點(Master)1個, 從節(jié)點(Slave)5個。構(gòu)建了港口和泊位的電子圍欄(如圖5b所示), 經(jīng)過對AIS數(shù)據(jù)的抽取轉(zhuǎn)換、清洗和挖掘分析(其中數(shù)據(jù)抽取與清洗耗時約17 h, 挖掘分析耗時約33 h), 共得到航次數(shù)據(jù)84 827條, 按船舶船型各航次分布情況如圖6所示。

    圖 6 按船舶類型各航次分布情況

    在航次數(shù)據(jù)基礎(chǔ)上, 形成了包含中國進出口船舶過主要海上通道情況、沿海港口進離港船舶情況, 以及煤炭、鐵礦石和原油3類大宗貨品進口情況的統(tǒng)計指標(biāo)。圖7和圖8是2019年3—5月大宗貨品船舶中國進出口所經(jīng)過主要海上通道的航次及載重統(tǒng)計情況, 其中載重是根據(jù)船舶檔案中的“滿載載重”累加計算得到。

    圖 7 中國進口船舶經(jīng)過海上通道統(tǒng)計

    對2019年3—5月中國沿海港口的進離港情況進行統(tǒng)計, 得出進離港船舶數(shù)量排前十的港口, 統(tǒng)計結(jié)果如圖9如示, 從圖中可以看出, 寧波舟山港無論是從艘次還是載重看, 都排第一位, 這與中國大宗貨品海運港口運量的排名基本一致, 進一步證明了一定程度上航次分析得出的“艘次”及“載重”大數(shù)據(jù)指標(biāo)可以從趨勢分析的層面代替?zhèn)鹘y(tǒng)統(tǒng)計的港口“吞吐量”指標(biāo)。

    圖 8 中國出口船舶經(jīng)過海上通道統(tǒng)計

    圖 9 沿海港口進離港船舶統(tǒng)計

    通過大宗貨品統(tǒng)計, 分析得出中國2019年3—5月的原油、鐵礦石和煤炭3類大宗貨品的進口量(表2),經(jīng)與海關(guān)公布的統(tǒng)計數(shù)據(jù)(網(wǎng)址: www.customs.gov.cn)進行對比發(fā)現(xiàn), 3類貨品的大數(shù)據(jù)統(tǒng)計值均小于海關(guān)統(tǒng)計數(shù)據(jù)(約小3%~20%), 且2019年3月數(shù)據(jù)與海關(guān)數(shù)據(jù)差距更大, 2019年4—5月數(shù)據(jù)與海關(guān)數(shù)據(jù)差距逐步減小, 經(jīng)分析, 應(yīng)是缺乏2019年3月之前的AIS數(shù)據(jù), 對航次的挖掘不夠完整造成的。但大數(shù)據(jù)統(tǒng)計的結(jié)果與官方統(tǒng)計的趨勢一致, 且統(tǒng)計結(jié)果的準(zhǔn)確性會隨著基礎(chǔ)數(shù)據(jù)的完善及算法的優(yōu)化不斷提高, 這也進一步從另一方面證明了本文提出的方法能滿足實際需要。

    表 2 中國大宗貨品進口量統(tǒng)計(×104 t)

    為進一步分析各港口、各貨品類型船舶的來源和流向, 基于大數(shù)據(jù)平臺對船舶AIS數(shù)據(jù)按航次進行了航跡線可視化, 部分樣例圖如圖10所示, 從航跡線可以明顯的看出中國進口鐵礦石的主要來源地包括澳大利亞、巴西、印度、加拿大等; 而通過航跡熱力圖可以反映大宗貨品海運中不同沿海港口或海上通道的繁忙程度。

    圖 10 AIS大數(shù)據(jù)可視化

    5 結(jié)論與展望

    本文面向港口、海上通道及大宗貨品海運大數(shù)據(jù)統(tǒng)計分析需求, 給出了基于AIS開展海運統(tǒng)計的技術(shù)方法、AIS大數(shù)據(jù)平臺的技術(shù)架構(gòu), 以及電子圍欄分析、航行事件挖掘、航次分析等關(guān)鍵模型算法?;?019年3—5月全球大宗貨品船舶的AIS數(shù)據(jù), 以及船舶檔案、港口及泊位數(shù)據(jù), 對本文提出的方法進行了實證分析。結(jié)果表明本文方法可以實現(xiàn)全球范圍內(nèi)海量船舶AIS數(shù)據(jù)的接入、清洗、處理及挖掘分析, 并在生成的航次數(shù)據(jù)基礎(chǔ)上, 實現(xiàn)了大宗貨品海運的港口統(tǒng)計、海上通道統(tǒng)計及貨品運量統(tǒng)計, 形成的大數(shù)據(jù)統(tǒng)計指標(biāo)一定程度上可以作為傳統(tǒng)統(tǒng)計方法及指標(biāo)的補充; 結(jié)合大數(shù)據(jù)可視化方法, 對航跡線基于地圖的展示還可以直觀反應(yīng)船舶、貨品的來源和流向, 為變革海運統(tǒng)計模式, 實現(xiàn)海運即時化、準(zhǔn)確化、精細(xì)化的大數(shù)據(jù)統(tǒng)計分析與展現(xiàn)提供了新思路。

    港口、泊位等基礎(chǔ)數(shù)據(jù)的完善是分析結(jié)果準(zhǔn)確性的保障, 但實際情況中, 數(shù)據(jù)的獲取及更新相當(dāng)困難, 因此, 引入機器學(xué)習(xí)、神經(jīng)網(wǎng)絡(luò)、關(guān)聯(lián)規(guī)則等方法, 基于已有航次數(shù)據(jù)和AIS船舶軌跡對港口、泊位等的位置及屬性數(shù)據(jù)進行補充、更新和完善是下一階段的研究方向; 此外, 由于模型和數(shù)據(jù)的限制, 本文并沒有考慮集裝箱船、漁船等船型船舶的海運特點及其AIS信息挖掘分析方法, 這也是在今后的工作中需要研究的重點。

    [1] 蘇興國. 歐盟海運統(tǒng)計簡介[J]. 中國統(tǒng)計, 2015(10): 30-31.

    SU Xingguo. Brief introduction of EU maritime statistics[J]. China Statistics, 2015(10): 30-31.

    [2] 王月恒, 王晶, 劉欣俠. 大數(shù)據(jù)時代新型對外貨物貿(mào)易統(tǒng)計體系研究[J]. 中國市場, 2017(19): 101, 103.

    Wang Yueheng, Wang Jing, Liu Xinxia. Research on the new statistical system of foreign trade in goods in the era of big data[J]. China Market, 2017(19): 101, 103.

    [3] 王杰, 唐順梅, 于海. 海運大國與海運強國辨析[J]. 大連海事大學(xué)學(xué)報(社會科學(xué)版), 2009(3): 44-47, 57.

    Wang Jie, Tang Shunmei, Yu Hai. Differentiation and analysis of maritime power and big maritime country[J]. Journal of Dalian Maritime University (Social Sciences Edition), 2009(3): 44-47, 57.

    [4] 徐叢春. 大數(shù)據(jù)時代海洋統(tǒng)計改革的思考[J]. 海洋經(jīng)濟, 2016, 6(6): 54-60.

    Xu Congchun. Reflections on the reform of marine statistics in the era of big data[J]. Marine Economy, 2016, 6(6): 54-60.

    [5] 趙銳. 我國海洋經(jīng)濟統(tǒng)計存在的問題及完善途徑分析[J]. 中國統(tǒng)計, 2013(2): 48-50.

    Zhao Rui. Analysis on the problems and improvement ways of China’s marine economic statistics[J]. China Statistics, 2013(2): 48-50.

    [6] 田璐, 張敖木翰, 張翼, 等. 船舶AIS大數(shù)據(jù)資源管理及分析應(yīng)用架構(gòu)設(shè)計[J]. 交通運輸研究, 2019, 5(5): 31-40.

    Tian Lu, Zhang Aomuhan, Zhang Yi, et al. Ship AIS big data management, analysis and application architecture design[J]. Transport Research, 2019, 5(5): 31-40.

    [7] 李紅祥, 方遜. 基于AIS的船舶交通流量統(tǒng)計方法研究[J]. 武漢理工大學(xué)學(xué)報(交通科學(xué)與工程版), 2011, 35(4): 853-857.

    Li Hongxiang, Fang Xun. Statistical method research of vessel traffic flow based on AIS[J]. Journal of Wuhan University of Technology (Transportation Science & Engineering), 2011, 35(4): 853-857.

    [8] 徐明強. 基于AIS數(shù)據(jù)的港口船舶營運效率評估系統(tǒng)[J]. 中國航海, 2018, 41(3): 119-122.

    Xu Mingqiang. Evaluation system of operation efficiency of ships in port based on AIS data[J]. Navigation of China, 2018, 41(3): 119-122.

    [9] 朱倩茹, 廖程浩, 王龍, 等. 基于AIS數(shù)據(jù)的精細(xì)化船舶排放清單方法[J]. 中國環(huán)境科學(xué), 2017, 37(12): 4493-4500.

    Zhu Qianru, Liao Chenghao, Wang Long, et al. Application of fine vessel emission inventory compilation method based on AIS data[J]. China Environmental Science, 2017, 37(12): 4493-4500.

    [10] Jalkanen J P, Johansson L, Kukkonen J. A comprehensive inventory of ship traffic exhaust emissions in the European sea areas in 2011[J]. Atmospheric Chemistry and Physics, 2016, 16: 71-84.

    [11] Maragkogianni A, Papaefthimious S. Eva-lua-ting the social cost of cruise ships air emissions in major ports of Greece[J]. Transportation Research-Part D, 2015, 36: 10-17.

    [12] Lonati G, Cernuschi S, Sidi S. Air quality impact assessment of at-berth ship emissions: Case-study for the project of a new freight port[J]. Science of the Total Environment, 2010, 409: 192-200.

    [13] 張建雄. 基于AIS的港口船舶動態(tài)統(tǒng)計系統(tǒng)[J]. 中國水運, 2015(2): 56-57.

    Zhang Jianxiong. Port ship dynamic statistics system based on AIS[J]. China Shipping, 2015(2): 56-57.

    [14] 何建春. 試述基于多傳感信息融合的內(nèi)河水運統(tǒng)計應(yīng)用研究[J]. 中國水運, 2017, 561(10): 39-40.

    He Jianchun. Application research on inland water transport statistics based on multi-sensor information fusion[J]. China Water Transport, 2017, 561(10): 39-40.

    [15] 劉滿娜. 基于AIS的港口監(jiān)控與分析系統(tǒng)的設(shè)計與實現(xiàn)[D]. 北京: 北京郵電大學(xué), 2014.

    Liu Manna. Design and implemrnt of port monitoring and analysis system based on AIS[D]. Beijing: Beijing University of Posts and Telecommunications, 2014.

    [16] 祝天虎. 全球鐵礦石貿(mào)易與航運分析系統(tǒng)設(shè)計與實現(xiàn)[D]. 北京: 北京郵電大學(xué), 2019.

    Zhu Tianhu. Design and implementation of global iron ore trade and shipping analysis system[D]. Beijing: Beijing University of Posts and Telecommunications, 2019.

    [17] 朱吉雙, 邱偉維, 菅文濤. 水運大數(shù)據(jù)開發(fā)與應(yīng)用研究——以船舶自動識別系統(tǒng)(AIS)為例[J]. 綜合運輸, 2020, 42(4): 33-37.

    Zhu Jishuang, Qiu Weiwei, JIan Wentao. Development and application of big data in the waterborne transportation—A case study of the ship automatic identification system[J]. China Transportation Review, 2020, 42(4): 33-37.

    AIS big data mining for maritime statistics

    ZHAO Long-fei1, 2, JIANG Xiao-yi1, 2, SUN Miao1, 2, CAO Lei1, 2, GUO Xue1, 2, Lü Chong-jing1, 2

    (1. National Marine Data and Information Service, Tianjin 300171, China; 2. Marine Information Technology Innovation Center of the Ministry of Natural Resources, Tianjin 300171, China)

    To address the problems of poor data quality, poor timeliness, outdated statistical products, and the unreliable traditional statistical methods of the maritime industry, a method of maritime statistics based on automatic identification system (AIS) big data mining analysis is proposed in this paper. The technical route of a maritime statistical analysis based on AIS big data, the technical architecture of big data platform, and the key model algorithm of electronic fence analysis, navigation event analysis, voyage analysis, and statistical index generation are elaborated in detail. With the application of AIS data of global bulk cargo ships from March to May in 2019 as an example, this method can provide big data statistical indicators of port, sea passage, and bulk cargo, as well as provide a new idea for realizing real-time, accurate, and refined statistical analysis and presentation of maritime transportation.

    maritime transport statistics; big data; marine economy; AIS; mining analysis

    Dec. 16, 2020

    TP399

    A

    1000-3096(2021)12-0055-10

    10.11759/hykx20201216003

    2020-12-16;

    2021-03-02

    國家重點研發(fā)計劃(2017YFC1405300)

    [National Key R & D plan, No. 2017YFC1405300]

    趙龍飛(1991—), 男, 山東鄆城人, 助理研究員, 碩士, 主要從事海洋大數(shù)據(jù)挖掘分析與海洋信息化研究, E-mail: zhaolongfei220@163.com; 曹磊(1987—),通信作者, 男, 工程師, 碩士, E-mail: 544341126@qq.com

    (本文編輯: 趙衛(wèi)紅)

    猜你喜歡
    航次貨品泊位
    多出/入口倉庫的貨位優(yōu)化研究
    化學(xué)品船適裝貨品的新要求及實船應(yīng)用
    中國船檢(2021年11期)2021-12-04 14:02:26
    湄洲灣港斗尾港區(qū)部分泊位竣工驗收
    水道港口(2016年3期)2016-04-07 13:50:11
    基于排隊論的區(qū)域路內(nèi)停車最優(yōu)泊位占用率研究
    考慮疲勞和工作負(fù)荷的人工揀選貨品排程研究
    Anti-ageing effects of a new Dimethylaminoethanol-based formulation on DGalactose induced skin ageing model of rat
    我國集裝箱航運企業(yè)實施作業(yè)成本管理法面臨的困難及解決方案
    集裝箱化(2014年10期)2014-10-31 18:26:46
    OBM型服裝企業(yè)電子商務(wù)貨品管理問題分析
    天津今年將新建對外開放泊位97個
    水道港口(2013年2期)2013-04-27 14:14:21
    制服诱惑二区| 热re99久久国产66热| 狂野欧美激情性bbbbbb| 一级爰片在线观看| 亚洲中文av在线| 最新的欧美精品一区二区| 一级毛片 在线播放| 啦啦啦中文免费视频观看日本| 国产激情久久老熟女| 国语对白做爰xxxⅹ性视频网站| www日本在线高清视频| 久久精品国产a三级三级三级| 久久精品熟女亚洲av麻豆精品| 人人妻人人澡人人爽人人夜夜| 日产精品乱码卡一卡2卡三| 免费播放大片免费观看视频在线观看| 少妇被粗大猛烈的视频| 欧美xxⅹ黑人| 亚洲精品国产av蜜桃| 免费黄频网站在线观看国产| 青春草视频在线免费观看| 国产毛片在线视频| 少妇高潮的动态图| 久久精品人人爽人人爽视色| 国产精品久久久久久久电影| 人成视频在线观看免费观看| 亚洲一区二区三区欧美精品| 日本与韩国留学比较| 精品一区二区三区四区五区乱码 | 自线自在国产av| 九九爱精品视频在线观看| 亚洲伊人色综图| 国产xxxxx性猛交| 欧美激情 高清一区二区三区| 国产在线视频一区二区| 免费av不卡在线播放| 国产福利在线免费观看视频| 2022亚洲国产成人精品| 美女主播在线视频| 美女国产视频在线观看| 日日爽夜夜爽网站| 亚洲美女搞黄在线观看| 亚洲久久久国产精品| 欧美亚洲 丝袜 人妻 在线| 久久久精品94久久精品| 黄网站色视频无遮挡免费观看| 亚洲丝袜综合中文字幕| 精品少妇黑人巨大在线播放| 五月伊人婷婷丁香| 亚洲色图综合在线观看| 乱人伦中国视频| 成人免费观看视频高清| 亚洲av男天堂| 成年人午夜在线观看视频| 国产色爽女视频免费观看| 国产有黄有色有爽视频| 少妇的逼好多水| 免费日韩欧美在线观看| 久久久精品免费免费高清| 亚洲av欧美aⅴ国产| 国产欧美另类精品又又久久亚洲欧美| 国产成人a∨麻豆精品| 亚洲国产成人一精品久久久| 在线天堂最新版资源| 99久久中文字幕三级久久日本| 狂野欧美激情性xxxx在线观看| 国产精品熟女久久久久浪| 久久久国产欧美日韩av| 永久免费av网站大全| 欧美精品国产亚洲| 水蜜桃什么品种好| 成人午夜精彩视频在线观看| 街头女战士在线观看网站| 成人18禁高潮啪啪吃奶动态图| 老司机亚洲免费影院| 午夜福利视频在线观看免费| 成人毛片a级毛片在线播放| 大码成人一级视频| 国产精品久久久久久精品电影小说| 国产高清三级在线| 寂寞人妻少妇视频99o| 日韩大片免费观看网站| 久久国内精品自在自线图片| 久久久精品免费免费高清| 大话2 男鬼变身卡| 亚洲欧洲精品一区二区精品久久久 | 日韩欧美一区视频在线观看| 国产精品偷伦视频观看了| 久久影院123| 婷婷成人精品国产| 夫妻性生交免费视频一级片| 精品久久久久久电影网| 下体分泌物呈黄色| 欧美人与善性xxx| 制服诱惑二区| videossex国产| 久久久a久久爽久久v久久| 欧美xxⅹ黑人| 多毛熟女@视频| 三级国产精品片| 久久av网站| 亚洲熟女精品中文字幕| 一本久久精品| 精品视频人人做人人爽| 国产精品国产三级专区第一集| 狂野欧美激情性xxxx在线观看| 51国产日韩欧美| 国产又色又爽无遮挡免| 最新中文字幕久久久久| 色婷婷久久久亚洲欧美| 亚洲国产精品999| 欧美日韩视频精品一区| 18禁在线无遮挡免费观看视频| 亚洲国产欧美日韩在线播放| 国产精品久久久久久久电影| 欧美精品av麻豆av| 激情视频va一区二区三区| 免费大片黄手机在线观看| 日本wwww免费看| 美国免费a级毛片| 久久狼人影院| 亚洲天堂av无毛| 美女福利国产在线| 美女中出高潮动态图| 亚洲一区二区三区欧美精品| 亚洲精品久久午夜乱码| 国产高清三级在线| 亚洲美女黄色视频免费看| 国精品久久久久久国模美| 精品亚洲乱码少妇综合久久| 欧美精品人与动牲交sv欧美| 免费观看性生交大片5| 日日啪夜夜爽| 亚洲内射少妇av| 在现免费观看毛片| 日韩成人伦理影院| 最黄视频免费看| 一级毛片我不卡| h视频一区二区三区| 国产极品天堂在线| 午夜福利影视在线免费观看| 色5月婷婷丁香| 看非洲黑人一级黄片| 精品卡一卡二卡四卡免费| 国产色爽女视频免费观看| 国产精品 国内视频| www.熟女人妻精品国产 | 欧美日韩成人在线一区二区| 免费久久久久久久精品成人欧美视频 | 亚洲国产最新在线播放| 国产色婷婷99| 亚洲天堂av无毛| 丰满少妇做爰视频| 欧美老熟妇乱子伦牲交| 中文字幕人妻熟女乱码| 国产成人aa在线观看| 伦精品一区二区三区| 亚洲欧洲国产日韩| 狂野欧美激情性bbbbbb| 亚洲国产日韩一区二区| 国产亚洲欧美精品永久| 少妇熟女欧美另类| 最近中文字幕2019免费版| 国产欧美另类精品又又久久亚洲欧美| 国产 精品1| 捣出白浆h1v1| 免费久久久久久久精品成人欧美视频 | 久久狼人影院| 在线观看免费日韩欧美大片| 春色校园在线视频观看| 高清视频免费观看一区二区| 少妇的丰满在线观看| 在线免费观看不下载黄p国产| 热99久久久久精品小说推荐| 人人妻人人澡人人看| 欧美 亚洲 国产 日韩一| 午夜免费男女啪啪视频观看| 亚洲av.av天堂| 久久久久精品性色| 亚洲中文av在线| 又粗又硬又长又爽又黄的视频| 免费日韩欧美在线观看| videosex国产| 久久午夜综合久久蜜桃| 十八禁高潮呻吟视频| 天天躁夜夜躁狠狠久久av| 国产综合精华液| 男女午夜视频在线观看 | 超碰97精品在线观看| 99国产综合亚洲精品| √禁漫天堂资源中文www| 黄色 视频免费看| 全区人妻精品视频| 啦啦啦中文免费视频观看日本| 97在线人人人人妻| 制服诱惑二区| 曰老女人黄片| 国产有黄有色有爽视频| 最近中文字幕高清免费大全6| 高清毛片免费看| 最近最新中文字幕大全免费视频 | 国产在线免费精品| 国产成人精品婷婷| 高清在线视频一区二区三区| 中文欧美无线码| 亚洲婷婷狠狠爱综合网| 久久亚洲国产成人精品v| av片东京热男人的天堂| 国产成人欧美| 2022亚洲国产成人精品| 精品一区二区三区视频在线| 精品国产一区二区三区久久久樱花| 高清av免费在线| 纯流量卡能插随身wifi吗| 成年美女黄网站色视频大全免费| 日韩成人av中文字幕在线观看| 国产深夜福利视频在线观看| 一级片免费观看大全| 91在线精品国自产拍蜜月| 久久国产精品男人的天堂亚洲 | 伊人久久国产一区二区| 欧美亚洲 丝袜 人妻 在线| 性色avwww在线观看| 精品酒店卫生间| 久久久久视频综合| 乱人伦中国视频| 国产极品天堂在线| 久久久久精品性色| 免费观看在线日韩| 嫩草影院入口| 中国美白少妇内射xxxbb| 美女国产视频在线观看| 三上悠亚av全集在线观看| 母亲3免费完整高清在线观看 | 日日摸夜夜添夜夜爱| 日韩视频在线欧美| 美女国产视频在线观看| 久久久久久久国产电影| 亚洲av国产av综合av卡| 五月天丁香电影| 欧美激情极品国产一区二区三区 | 亚洲丝袜综合中文字幕| 桃花免费在线播放| 黄网站色视频无遮挡免费观看| 十八禁网站网址无遮挡| 国产片内射在线| 久久久国产欧美日韩av| 国产麻豆69| 韩国高清视频一区二区三区| 免费观看av网站的网址| 精品卡一卡二卡四卡免费| 国产精品国产av在线观看| 亚洲国产av新网站| 丝袜在线中文字幕| 男女边摸边吃奶| 国产免费一区二区三区四区乱码| 欧美激情 高清一区二区三区| 美女国产视频在线观看| 成年av动漫网址| 亚洲一区二区三区欧美精品| 老司机亚洲免费影院| 看十八女毛片水多多多| 青春草亚洲视频在线观看| 国产亚洲av片在线观看秒播厂| 在线看a的网站| 国产1区2区3区精品| 欧美变态另类bdsm刘玥| 精品午夜福利在线看| 免费看av在线观看网站| 天天躁夜夜躁狠狠久久av| 咕卡用的链子| 日韩熟女老妇一区二区性免费视频| 多毛熟女@视频| 精品亚洲成国产av| 亚洲美女视频黄频| a级毛色黄片| 在线看a的网站| 亚洲成人手机| 亚洲国产日韩一区二区| 亚洲国产精品一区三区| 成人午夜精彩视频在线观看| 免费看av在线观看网站| 国产免费福利视频在线观看| 卡戴珊不雅视频在线播放| 尾随美女入室| 亚洲欧洲精品一区二区精品久久久 | 国产精品秋霞免费鲁丝片| 人妻系列 视频| 永久网站在线| 最后的刺客免费高清国语| 黄色配什么色好看| 久久人人爽av亚洲精品天堂| 亚洲人成网站在线观看播放| 国产精品免费大片| 成年动漫av网址| 国产成人免费观看mmmm| 久久久久久久大尺度免费视频| 寂寞人妻少妇视频99o| 美女内射精品一级片tv| 国产精品.久久久| 国产综合精华液| 久久毛片免费看一区二区三区| 日韩制服骚丝袜av| 成年av动漫网址| 午夜av观看不卡| 国产成人aa在线观看| 国产成人精品久久久久久| 99九九在线精品视频| 亚洲国产日韩一区二区| 久久97久久精品| 亚洲国产精品999| 欧美少妇被猛烈插入视频| 国产xxxxx性猛交| 一级毛片 在线播放| 自线自在国产av| 免费观看a级毛片全部| 天堂俺去俺来也www色官网| 午夜福利网站1000一区二区三区| 丰满少妇做爰视频| 插逼视频在线观看| 国产av国产精品国产| 蜜桃在线观看..| 一级a做视频免费观看| 免费少妇av软件| 国产精品秋霞免费鲁丝片| 国产精品蜜桃在线观看| 午夜av观看不卡| 99热这里只有是精品在线观看| 最近最新中文字幕大全免费视频 | 久久久亚洲精品成人影院| 春色校园在线视频观看| 国产福利在线免费观看视频| 亚洲欧美清纯卡通| 99国产综合亚洲精品| 丝袜美足系列| 免费av不卡在线播放| 国产免费视频播放在线视频| 色94色欧美一区二区| av国产精品久久久久影院| 日韩成人伦理影院| 天堂8中文在线网| 亚洲精品一二三| 国产激情久久老熟女| 成人综合一区亚洲| 午夜激情久久久久久久| 精品一区二区三卡| 一边亲一边摸免费视频| 成年女人在线观看亚洲视频| 色哟哟·www| 日韩电影二区| 欧美少妇被猛烈插入视频| 午夜免费观看性视频| 国产男人的电影天堂91| 亚洲丝袜综合中文字幕| 精品一区在线观看国产| 热99久久久久精品小说推荐| 亚洲av电影在线观看一区二区三区| 久久婷婷青草| 亚洲精品国产色婷婷电影| 亚洲人成网站在线观看播放| 国产精品一区二区在线不卡| a级片在线免费高清观看视频| 中文精品一卡2卡3卡4更新| 午夜福利视频在线观看免费| 欧美精品国产亚洲| 午夜免费观看性视频| 精品久久国产蜜桃| 亚洲精品一二三| 搡老乐熟女国产| 少妇被粗大的猛进出69影院 | 久久精品国产a三级三级三级| 亚洲天堂av无毛| 国产一级毛片在线| 少妇人妻 视频| 国产爽快片一区二区三区| 婷婷色综合大香蕉| 久久久久久人人人人人| 极品人妻少妇av视频| 伊人久久国产一区二区| 久久久久视频综合| 在线观看免费视频网站a站| 国产男人的电影天堂91| 国产av一区二区精品久久| 男女边吃奶边做爰视频| 亚洲,欧美精品.| 国产日韩一区二区三区精品不卡| 午夜视频国产福利| 狂野欧美激情性xxxx在线观看| 久热这里只有精品99| 久久人人爽人人爽人人片va| 精品视频人人做人人爽| 国产片内射在线| 宅男免费午夜| 亚洲色图综合在线观看| 亚洲成人av在线免费| 日韩免费高清中文字幕av| 成年美女黄网站色视频大全免费| 性色avwww在线观看| 中文字幕另类日韩欧美亚洲嫩草| 午夜91福利影院| 欧美日韩综合久久久久久| 人妻人人澡人人爽人人| 成人毛片60女人毛片免费| 在线 av 中文字幕| 日韩一本色道免费dvd| 欧美成人午夜精品| 成人黄色视频免费在线看| 婷婷成人精品国产| 精品一区二区三区四区五区乱码 | 欧美日韩视频精品一区| 欧美成人午夜免费资源| 99久久中文字幕三级久久日本| av黄色大香蕉| xxx大片免费视频| 成年女人在线观看亚洲视频| 亚洲国产精品一区三区| 大香蕉久久网| 久久韩国三级中文字幕| 精品国产露脸久久av麻豆| 春色校园在线视频观看| 天天躁夜夜躁狠狠久久av| 成人亚洲精品一区在线观看| 久久久久久久精品精品| 色94色欧美一区二区| 午夜福利影视在线免费观看| 美女大奶头黄色视频| 亚洲欧美日韩另类电影网站| 国产熟女午夜一区二区三区| av在线app专区| 亚洲精品色激情综合| 男女无遮挡免费网站观看| 男女边吃奶边做爰视频| 欧美激情国产日韩精品一区| 好男人视频免费观看在线| 天堂中文最新版在线下载| 人妻少妇偷人精品九色| 国产一区二区三区av在线| 久久人人97超碰香蕉20202| 女人精品久久久久毛片| 久久狼人影院| 9色porny在线观看| 亚洲中文av在线| videossex国产| 在线看a的网站| 99热全是精品| 视频区图区小说| 少妇精品久久久久久久| 精品一区在线观看国产| 九草在线视频观看| av电影中文网址| 亚洲av欧美aⅴ国产| 久久99热6这里只有精品| 久久精品国产自在天天线| 老熟女久久久| 国产男女内射视频| 欧美激情国产日韩精品一区| 成人无遮挡网站| 成年美女黄网站色视频大全免费| 亚洲国产色片| 亚洲精品,欧美精品| 狠狠精品人妻久久久久久综合| 最新中文字幕久久久久| 天天影视国产精品| 精品第一国产精品| 国产男女超爽视频在线观看| 97精品久久久久久久久久精品| 大话2 男鬼变身卡| 建设人人有责人人尽责人人享有的| 91aial.com中文字幕在线观看| 免费看av在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 久久久国产欧美日韩av| 有码 亚洲区| 宅男免费午夜| 亚洲成人手机| av网站免费在线观看视频| a级毛色黄片| 精品国产一区二区久久| 日韩精品免费视频一区二区三区 | 国产xxxxx性猛交| 亚洲国产精品一区三区| 欧美最新免费一区二区三区| 色婷婷久久久亚洲欧美| 国产精品无大码| 韩国精品一区二区三区 | 亚洲欧美成人精品一区二区| 啦啦啦视频在线资源免费观看| 大陆偷拍与自拍| 中文乱码字字幕精品一区二区三区| 一级,二级,三级黄色视频| 晚上一个人看的免费电影| 成人国产av品久久久| 亚洲国产欧美日韩在线播放| 一本色道久久久久久精品综合| 午夜福利乱码中文字幕| 久久青草综合色| 人体艺术视频欧美日本| 少妇人妻精品综合一区二区| 日韩不卡一区二区三区视频在线| 丁香六月天网| 成人18禁高潮啪啪吃奶动态图| 成年人免费黄色播放视频| av电影中文网址| freevideosex欧美| 亚洲国产精品国产精品| 九色成人免费人妻av| 国产在线视频一区二区| 亚洲av免费高清在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲图色成人| 日本猛色少妇xxxxx猛交久久| 午夜福利,免费看| 最后的刺客免费高清国语| 男女边摸边吃奶| 久久人人97超碰香蕉20202| 日韩av免费高清视频| 成年人免费黄色播放视频| 色婷婷久久久亚洲欧美| av免费观看日本| 天天影视国产精品| 亚洲情色 制服丝袜| 老司机影院毛片| 熟女人妻精品中文字幕| 久久精品国产自在天天线| 青春草视频在线免费观看| 久久青草综合色| 在线看a的网站| 丝袜在线中文字幕| 久久久国产一区二区| 99re6热这里在线精品视频| 国语对白做爰xxxⅹ性视频网站| 边亲边吃奶的免费视频| 18禁观看日本| 男女高潮啪啪啪动态图| 夫妻性生交免费视频一级片| 最近的中文字幕免费完整| 少妇人妻 视频| 成人午夜精彩视频在线观看| 18禁在线无遮挡免费观看视频| 免费女性裸体啪啪无遮挡网站| 一级黄片播放器| www.色视频.com| av视频免费观看在线观看| 嫩草影院入口| 蜜臀久久99精品久久宅男| 乱码一卡2卡4卡精品| 搡老乐熟女国产| 精品国产露脸久久av麻豆| 久久久精品区二区三区| 日韩中文字幕视频在线看片| 亚洲久久久国产精品| 国产老妇伦熟女老妇高清| av免费观看日本| 久久婷婷青草| 亚洲婷婷狠狠爱综合网| 蜜臀久久99精品久久宅男| 久久久国产欧美日韩av| 久久久久久久精品精品| 中文天堂在线官网| av线在线观看网站| 成人漫画全彩无遮挡| 丰满迷人的少妇在线观看| 美女脱内裤让男人舔精品视频| 侵犯人妻中文字幕一二三四区| 久久久久精品久久久久真实原创| 日本欧美视频一区| 91在线精品国自产拍蜜月| 人成视频在线观看免费观看| 高清毛片免费看| 久久影院123| 天堂中文最新版在线下载| 久久久久久人人人人人| 欧美少妇被猛烈插入视频| 人体艺术视频欧美日本| 亚洲精品456在线播放app| 午夜激情久久久久久久| 久久精品人人爽人人爽视色| 日韩一区二区三区影片| 女人精品久久久久毛片| 十八禁网站网址无遮挡| 黄色视频在线播放观看不卡| 日本vs欧美在线观看视频| 国产午夜精品一二区理论片| 人人澡人人妻人| 波野结衣二区三区在线| 国产精品99久久99久久久不卡 | 伦精品一区二区三区| 色婷婷久久久亚洲欧美| 亚洲av成人精品一二三区| 熟女av电影| 午夜老司机福利剧场| 成人黄色视频免费在线看| av在线老鸭窝| 亚洲精品av麻豆狂野| 亚洲一区二区三区欧美精品| 久久久a久久爽久久v久久| 久久97久久精品| 婷婷色麻豆天堂久久| 国产乱人偷精品视频| 日韩成人av中文字幕在线观看| 大话2 男鬼变身卡| 久久精品国产亚洲av天美| 两个人看的免费小视频| 最近中文字幕2019免费版| 香蕉精品网在线| 国产视频首页在线观看| 一级毛片黄色毛片免费观看视频| 天堂俺去俺来也www色官网| 高清av免费在线| 久久女婷五月综合色啪小说| 亚洲成人一二三区av| 国产精品欧美亚洲77777| 欧美日韩一区二区视频在线观看视频在线| 人妻少妇偷人精品九色| av片东京热男人的天堂|