• <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
    夜夜爽夜夜爽视频| 男插女下体视频免费在线播放| 18+在线观看网站| 久久久久久伊人网av| 亚洲久久久久久中文字幕| 三级国产精品欧美在线观看| 一级a做视频免费观看| 久久久久久久久久成人| 只有这里有精品99| 最近的中文字幕免费完整| 日产精品乱码卡一卡2卡三| 69av精品久久久久久| 日韩 亚洲 欧美在线| av在线观看视频网站免费| 免费黄网站久久成人精品| 热99在线观看视频| 五月天丁香电影| 男女啪啪激烈高潮av片| 精品久久久久久久久久久久久| 日日摸夜夜添夜夜爱| 国产真实伦视频高清在线观看| 精品久久久噜噜| a级毛色黄片| 日本熟妇午夜| 一级毛片 在线播放| 国产乱来视频区| 日韩av不卡免费在线播放| 一级毛片aaaaaa免费看小| 18禁动态无遮挡网站| 国产高清有码在线观看视频| 美女大奶头视频| 午夜免费男女啪啪视频观看| 亚洲国产最新在线播放| 亚洲成人中文字幕在线播放| 99久国产av精品| av国产久精品久网站免费入址| 午夜福利高清视频| 亚洲三级黄色毛片| 国产精品久久久久久av不卡| 麻豆久久精品国产亚洲av| 色5月婷婷丁香| 午夜福利成人在线免费观看| 2021少妇久久久久久久久久久| 亚洲精品色激情综合| 免费播放大片免费观看视频在线观看| h日本视频在线播放| 色5月婷婷丁香| 观看免费一级毛片| 久久99热6这里只有精品| 69av精品久久久久久| 舔av片在线| 18禁在线播放成人免费| 国产成人freesex在线| 韩国高清视频一区二区三区| av在线播放精品| 国产精品国产三级国产av玫瑰| 国产精品一区二区三区四区免费观看| 特大巨黑吊av在线直播| 国语对白做爰xxxⅹ性视频网站| 国产av不卡久久| 国产三级在线视频| 久久久久久久久久成人| 精品99又大又爽又粗少妇毛片| 国产精品久久久久久av不卡| 日韩av免费高清视频| 最新中文字幕久久久久| 久久久久免费精品人妻一区二区| 18禁动态无遮挡网站| 一本一本综合久久| 色播亚洲综合网| 日韩精品有码人妻一区| 欧美xxⅹ黑人| 亚洲人成网站在线观看播放| 一级av片app| 五月天丁香电影| 女人十人毛片免费观看3o分钟| 亚洲自偷自拍三级| 美女脱内裤让男人舔精品视频| 免费观看精品视频网站| 欧美xxxx性猛交bbbb| 大片免费播放器 马上看| 一边亲一边摸免费视频| 日日撸夜夜添| 亚洲一级一片aⅴ在线观看| av免费在线看不卡| 国产真实伦视频高清在线观看| 边亲边吃奶的免费视频| 日本黄大片高清| 国产精品综合久久久久久久免费| 伊人久久精品亚洲午夜| 一级爰片在线观看| 国产乱来视频区| 国产精品美女特级片免费视频播放器| 亚洲国产色片| 午夜福利在线观看吧| 大片免费播放器 马上看| 69av精品久久久久久| 久久久久久久久大av| 亚洲自拍偷在线| 麻豆久久精品国产亚洲av| 精品国产一区二区三区久久久樱花 | 老司机影院毛片| 日韩,欧美,国产一区二区三区| 又黄又爽又刺激的免费视频.| 嘟嘟电影网在线观看| 久久精品夜色国产| 午夜福利视频1000在线观看| 久久韩国三级中文字幕| 听说在线观看完整版免费高清| 国产高清国产精品国产三级 | 一区二区三区高清视频在线| 欧美不卡视频在线免费观看| 午夜免费男女啪啪视频观看| 免费无遮挡裸体视频| 久久久久久久大尺度免费视频| 男女视频在线观看网站免费| 免费看a级黄色片| 熟妇人妻久久中文字幕3abv| 欧美成人午夜免费资源| av女优亚洲男人天堂| 天美传媒精品一区二区| 国产高清国产精品国产三级 | 在线观看人妻少妇| 亚洲精品中文字幕在线视频 | 男人舔女人下体高潮全视频| 欧美日韩在线观看h| 神马国产精品三级电影在线观看| 91久久精品电影网| 少妇熟女aⅴ在线视频| 国产淫语在线视频| 免费看日本二区| 国产69精品久久久久777片| 麻豆精品久久久久久蜜桃| 日本黄色片子视频| 免费人成在线观看视频色| 尾随美女入室| 欧美日本视频| www.av在线官网国产| 少妇熟女aⅴ在线视频| 日韩,欧美,国产一区二区三区| 国产男女超爽视频在线观看| 午夜爱爱视频在线播放| 男人舔奶头视频| 纵有疾风起免费观看全集完整版 | 欧美精品国产亚洲| 久久久精品94久久精品| 性插视频无遮挡在线免费观看| 午夜激情欧美在线| 精品久久久久久成人av| 看黄色毛片网站| 不卡视频在线观看欧美| 国产黄片美女视频| 男女啪啪激烈高潮av片| 亚洲av男天堂| 久久99热这里只频精品6学生| 国产中年淑女户外野战色| 亚洲精华国产精华液的使用体验| 看非洲黑人一级黄片| 91精品伊人久久大香线蕉| 精品久久久久久成人av| 22中文网久久字幕| 欧美一区二区亚洲| 欧美高清性xxxxhd video| 免费av不卡在线播放| 欧美精品国产亚洲| 国产乱来视频区| av在线亚洲专区| 欧美日韩在线观看h| 国产成人a∨麻豆精品| 亚洲天堂国产精品一区在线| 亚洲av中文av极速乱| 国产视频首页在线观看| 亚洲av在线观看美女高潮| www.av在线官网国产| 高清日韩中文字幕在线| 国产精品一区www在线观看| 亚洲国产日韩欧美精品在线观看| 99久久精品热视频| 三级毛片av免费| 天天一区二区日本电影三级| 日韩欧美三级三区| 国产成人aa在线观看| av在线亚洲专区| 3wmmmm亚洲av在线观看| 春色校园在线视频观看| 亚洲在久久综合| av国产免费在线观看| 最后的刺客免费高清国语| av在线老鸭窝| 国产一区二区三区av在线| 成人午夜高清在线视频| 国产免费又黄又爽又色| 亚洲真实伦在线观看| 天天躁日日操中文字幕| 免费观看精品视频网站| .国产精品久久| 午夜福利在线在线| 亚洲av日韩在线播放| 国产av在哪里看| 国产伦精品一区二区三区视频9| 成人国产麻豆网| 有码 亚洲区| 精品久久久精品久久久| 亚洲精品中文字幕在线视频 | av专区在线播放| 久久久久久久久久久免费av| 亚洲精品自拍成人| 国产黄a三级三级三级人| 搡老乐熟女国产| 国产男人的电影天堂91| 国产成人福利小说| 美女黄网站色视频| 淫秽高清视频在线观看| 有码 亚洲区| 一夜夜www| 禁无遮挡网站| 非洲黑人性xxxx精品又粗又长| 国产91av在线免费观看| 久久精品国产亚洲av天美| 日日啪夜夜爽| 99热这里只有精品一区| 亚洲天堂国产精品一区在线| 色尼玛亚洲综合影院| 夜夜爽夜夜爽视频| 亚洲丝袜综合中文字幕| 少妇的逼水好多| 国产成人aa在线观看| 亚州av有码| 久久久久精品久久久久真实原创| 搡老乐熟女国产| 99九九线精品视频在线观看视频| 国产精品人妻久久久影院| 最新中文字幕久久久久| 国产三级在线视频| 男女边摸边吃奶| 汤姆久久久久久久影院中文字幕 | 亚洲最大成人中文| 午夜免费观看性视频| 人妻一区二区av| 夫妻性生交免费视频一级片| 国产有黄有色有爽视频| 日韩欧美 国产精品| 中文欧美无线码| 精品久久国产蜜桃| 一区二区三区免费毛片| 99热6这里只有精品| 高清在线视频一区二区三区| 波多野结衣巨乳人妻| 内射极品少妇av片p| 天堂中文最新版在线下载 | 日本黄色片子视频| 久久精品国产亚洲av天美| 免费av观看视频| 国产成人一区二区在线| 欧美xxⅹ黑人| 久久国产乱子免费精品| 最近最新中文字幕免费大全7| 国产一区二区亚洲精品在线观看| 青青草视频在线视频观看| 亚洲四区av| 国产视频内射| 亚洲一区高清亚洲精品| 亚洲成人一二三区av| 成年av动漫网址| 欧美一级a爱片免费观看看| 亚洲精品一二三| 网址你懂的国产日韩在线| 一级毛片黄色毛片免费观看视频| 国产真实伦视频高清在线观看| 亚洲自偷自拍三级| 69人妻影院| 亚洲成人一二三区av| 尾随美女入室| 亚洲精品久久久久久婷婷小说| 色综合站精品国产| 午夜爱爱视频在线播放| 午夜视频国产福利| 亚洲精品久久久久久婷婷小说| 久久久久久久午夜电影| 国产精品久久久久久精品电影小说 | 亚洲色图av天堂| 美女黄网站色视频| 国产乱人偷精品视频| 美女大奶头视频| 狂野欧美白嫩少妇大欣赏| 熟女电影av网| 国产精品久久视频播放| 免费大片黄手机在线观看| 国产高清国产精品国产三级 | 亚洲图色成人| 国内揄拍国产精品人妻在线| 亚洲一区高清亚洲精品| 自拍偷自拍亚洲精品老妇| 99久国产av精品国产电影| 99久久九九国产精品国产免费| 久久综合国产亚洲精品| 人妻一区二区av| 欧美另类一区| 亚洲av中文字字幕乱码综合| 国产亚洲午夜精品一区二区久久 | 亚洲电影在线观看av| 国产av国产精品国产| 日日摸夜夜添夜夜添av毛片| 国产精品日韩av在线免费观看| 国产精品麻豆人妻色哟哟久久 | 欧美97在线视频| 亚洲怡红院男人天堂| 国产精品久久久久久久久免| 亚洲成人精品中文字幕电影| 欧美丝袜亚洲另类| ponron亚洲| 最近手机中文字幕大全| 精品国产一区二区三区久久久樱花 | 免费无遮挡裸体视频| 在线观看美女被高潮喷水网站| 亚洲丝袜综合中文字幕| 日本猛色少妇xxxxx猛交久久| 日韩欧美国产在线观看| av免费在线看不卡| 最近最新中文字幕大全电影3| 男人舔女人下体高潮全视频| 国产91av在线免费观看| 男女下面进入的视频免费午夜| 亚洲精品中文字幕在线视频 | 中文乱码字字幕精品一区二区三区 | 我要看日韩黄色一级片| 少妇的逼水好多| 少妇的逼好多水| 国产 一区 欧美 日韩| 99久久精品热视频| 日日啪夜夜爽| 色吧在线观看| 亚洲天堂国产精品一区在线| 国产 亚洲一区二区三区 | 中文字幕亚洲精品专区| 免费观看a级毛片全部| 99热全是精品| 亚洲精品久久久久久婷婷小说| 水蜜桃什么品种好| 日韩精品青青久久久久久| 熟女人妻精品中文字幕| av网站免费在线观看视频 | 亚洲经典国产精华液单| 嫩草影院入口| 国产男女超爽视频在线观看| 一区二区三区乱码不卡18| 99re6热这里在线精品视频| 国产成人免费观看mmmm| 午夜福利成人在线免费观看| 精品国产三级普通话版| 国产精品国产三级国产专区5o| 国产三级在线视频| 国产麻豆成人av免费视频| 国产精品综合久久久久久久免费| 国产麻豆成人av免费视频| 国产精品综合久久久久久久免费| 欧美xxxx性猛交bbbb| 22中文网久久字幕| 少妇熟女欧美另类| av在线观看视频网站免费| 韩国av在线不卡| 午夜福利在线在线| 干丝袜人妻中文字幕| 日本与韩国留学比较| 国产精品一二三区在线看| 国产爱豆传媒在线观看| 亚洲人与动物交配视频| 只有这里有精品99| 久热久热在线精品观看| 亚洲欧洲国产日韩| 久久久久久久久久久丰满| 菩萨蛮人人尽说江南好唐韦庄| 国产91av在线免费观看| 美女国产视频在线观看| 观看免费一级毛片| 天堂√8在线中文| 99热这里只有是精品50| 国产黄a三级三级三级人| 国产欧美日韩精品一区二区| or卡值多少钱| 国产探花极品一区二区| 少妇猛男粗大的猛烈进出视频 | 秋霞伦理黄片| 精品熟女少妇av免费看| 久久久久久国产a免费观看| 国产麻豆成人av免费视频| 少妇熟女aⅴ在线视频| 又粗又硬又长又爽又黄的视频| 免费看日本二区| 黄色日韩在线| 亚洲国产av新网站| 亚洲国产最新在线播放| 大又大粗又爽又黄少妇毛片口| 天天躁夜夜躁狠狠久久av| 国产高清有码在线观看视频| 日本三级黄在线观看| 身体一侧抽搐| 国产男人的电影天堂91| 99热这里只有是精品在线观看| 国产成人精品福利久久| 亚洲怡红院男人天堂| 最后的刺客免费高清国语| 国产精品国产三级国产专区5o| 亚洲av中文字字幕乱码综合| 日日啪夜夜撸| 国产一区二区亚洲精品在线观看| 亚洲一区高清亚洲精品| 三级男女做爰猛烈吃奶摸视频| 久久久国产一区二区| 蜜桃久久精品国产亚洲av| 国内揄拍国产精品人妻在线| 国产 一区 欧美 日韩| xxx大片免费视频| 久久精品熟女亚洲av麻豆精品 | 久热久热在线精品观看| 精品久久久久久电影网| 国产色婷婷99| 嫩草影院精品99| 精品一区二区三卡| 亚洲成人中文字幕在线播放| 老司机影院成人| 麻豆成人av视频| 在线免费观看的www视频| 最近最新中文字幕大全电影3| 中文精品一卡2卡3卡4更新| 99久久精品一区二区三区| 久久久精品欧美日韩精品| 十八禁网站网址无遮挡 | 国产精品一区二区性色av| 成年人午夜在线观看视频 | 亚洲欧美清纯卡通| 在线观看av片永久免费下载| 干丝袜人妻中文字幕| 黄片wwwwww| 亚洲精品乱码久久久v下载方式| 又大又黄又爽视频免费| 久久久久久久亚洲中文字幕| 亚洲精品国产成人久久av| 特级一级黄色大片| videossex国产| 国产免费福利视频在线观看| 国产麻豆成人av免费视频| 亚洲欧美精品专区久久| 直男gayav资源| av又黄又爽大尺度在线免费看| 国产伦精品一区二区三区四那| av在线天堂中文字幕| 日韩国内少妇激情av| 欧美xxxx黑人xx丫x性爽| 免费观看av网站的网址| 秋霞在线观看毛片| 精品久久久久久久久亚洲| 日韩欧美国产在线观看| 久久久久久久午夜电影| 搞女人的毛片| 午夜福利在线观看吧| 天堂中文最新版在线下载 | 午夜福利成人在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久综合国产亚洲精品| 午夜福利视频1000在线观看| 国产淫语在线视频| 亚洲熟女精品中文字幕| 国产精品久久久久久av不卡| 国产淫语在线视频| 国内揄拍国产精品人妻在线| 久久人人爽人人片av| 欧美日韩综合久久久久久| 男女国产视频网站| 免费黄网站久久成人精品| 亚洲精品日韩在线中文字幕| 午夜福利视频1000在线观看| 久久久久久久久久人人人人人人| 天美传媒精品一区二区| 亚洲欧美成人精品一区二区| 国产av不卡久久| 极品少妇高潮喷水抽搐| 看非洲黑人一级黄片| 97精品久久久久久久久久精品| 国产高清不卡午夜福利| 搡老妇女老女人老熟妇| 天堂av国产一区二区熟女人妻| 又大又黄又爽视频免费| 麻豆精品久久久久久蜜桃| 我的老师免费观看完整版| 久久韩国三级中文字幕| 精品少妇黑人巨大在线播放| 男女国产视频网站| 国产在视频线在精品| 婷婷六月久久综合丁香| 中文字幕亚洲精品专区| 黄片无遮挡物在线观看| 在线免费十八禁| 国产成人精品久久久久久| 久久韩国三级中文字幕| 国产伦一二天堂av在线观看| 精品久久久久久久久av| 国产在视频线在精品| 97人妻精品一区二区三区麻豆| 国产片特级美女逼逼视频| av天堂中文字幕网| 午夜福利在线观看吧| 国产精品1区2区在线观看.| 免费观看无遮挡的男女| 黄色配什么色好看| 少妇裸体淫交视频免费看高清| 国产成人午夜福利电影在线观看| 色哟哟·www| 免费黄频网站在线观看国产| 久久精品夜夜夜夜夜久久蜜豆| 人妻系列 视频| 国产免费一级a男人的天堂| 少妇猛男粗大的猛烈进出视频 | 中文字幕亚洲精品专区| freevideosex欧美| 国产片特级美女逼逼视频| 亚洲av二区三区四区| 嫩草影院精品99| 青春草视频在线免费观看| av免费在线看不卡| 特大巨黑吊av在线直播| 亚洲av国产av综合av卡| 国产老妇伦熟女老妇高清| 国产精品一区二区在线观看99 | 男女视频在线观看网站免费| 国产色婷婷99| 99久久精品热视频| 女人被狂操c到高潮| 国产高潮美女av| 一区二区三区乱码不卡18| 岛国毛片在线播放| 听说在线观看完整版免费高清| 国产免费又黄又爽又色| 婷婷六月久久综合丁香| 国产一级毛片在线| 亚洲精品成人av观看孕妇| 精品国产露脸久久av麻豆 | 国内精品宾馆在线| 伦理电影大哥的女人| 在线观看美女被高潮喷水网站| 免费看不卡的av| 激情 狠狠 欧美| 欧美一级a爱片免费观看看| 欧美精品国产亚洲| 99热这里只有是精品50| 精品一区二区免费观看| 中文字幕亚洲精品专区| 亚洲成人精品中文字幕电影| 最近2019中文字幕mv第一页| 亚洲欧美精品自产自拍| 高清欧美精品videossex| 黄色日韩在线| 久久6这里有精品| 熟妇人妻久久中文字幕3abv| 日韩强制内射视频| 免费人成在线观看视频色| 97超视频在线观看视频| 乱码一卡2卡4卡精品| 少妇丰满av| 亚洲天堂国产精品一区在线| 国产一区亚洲一区在线观看| 国产成人精品一,二区| 久久久久久久午夜电影| 国产午夜精品论理片| 久久97久久精品| 亚洲色图av天堂| 一本久久精品| 精品熟女少妇av免费看| 在现免费观看毛片| 国产午夜福利久久久久久| av国产久精品久网站免费入址| 欧美激情久久久久久爽电影| 日日摸夜夜添夜夜爱| 成人亚洲精品一区在线观看 | 色视频www国产| 亚洲欧美日韩无卡精品| 男人狂女人下面高潮的视频| 自拍偷自拍亚洲精品老妇| 欧美丝袜亚洲另类| 2021天堂中文幕一二区在线观| 美女高潮的动态| 午夜日本视频在线| 在线免费观看不下载黄p国产| 国产黄色免费在线视频| 简卡轻食公司| 免费观看在线日韩| 一本一本综合久久| 91久久精品电影网| 少妇人妻精品综合一区二区| 日韩一本色道免费dvd| 一级黄片播放器| 亚洲av中文av极速乱| 久久久久网色| 1000部很黄的大片| 国产淫片久久久久久久久| 我的老师免费观看完整版| 成人性生交大片免费视频hd| 国产爱豆传媒在线观看| 国产乱来视频区| 亚洲国产精品sss在线观看| 99久国产av精品| 少妇丰满av| av线在线观看网站| 爱豆传媒免费全集在线观看| 欧美潮喷喷水| 亚洲最大成人av| 欧美高清成人免费视频www| 亚洲国产欧美人成| a级一级毛片免费在线观看| 美女xxoo啪啪120秒动态图| 亚洲,欧美,日韩| 69人妻影院| h日本视频在线播放| 春色校园在线视频观看| 日本黄大片高清|