李世友

(湖南信息學院 電子科學與工程學院，長沙 410151)

船舶自動識別系統(tǒng)(Automatic Identification System, AIS)[1-2]采用自組織時分多址(Self-organizing Time Division Multiple Access，SOTDMA)通信協(xié)議[2]，自動地播發(fā)出本船的靜態(tài)、動態(tài)和航次信息，達到規(guī)避船舶碰撞、領航調(diào)度和航運管理等航行決策的目的, 主要用于對海上航行的船舶進行動態(tài)跟蹤，實現(xiàn)船舶標識、位置、航向、航速等航行信息的自動交換，用于領航及避撞[3]。

星載AIS系統(tǒng)[4-7]有著全天時、全天候、遠距離、大范圍、多參數(shù)等特點，以及衛(wèi)星探測的合法性、低軌道衛(wèi)星發(fā)射不受國籍限制等優(yōu)點，因此星載AIS在海洋監(jiān)管領域具有獨到的優(yōu)勢[6-8]。衛(wèi)星AIS技術彌補了岸基AIS系統(tǒng)在地域覆蓋上的局限性，拓展數(shù)據(jù)接收范圍，進而實現(xiàn)船舶全球監(jiān)控。星載AIS系統(tǒng)信息在國民經(jīng)濟建設中可用于掌握全球船舶動態(tài)、維護航行安全、分析全球經(jīng)濟態(tài)勢等，在軍事上可用于監(jiān)視全球海域艦船目標，特別是在與其他手段配合的基礎上，可以提高定位精度、準確判定目標屬性，并生成全球海洋目標態(tài)勢圖。目前，美國、挪威、德國、日本等多個國家和組織發(fā)射了自己的AIS衛(wèi)星[9]；我國也積極開展AIS相關技術的研究，并于2012年發(fā)射了第一顆AIS海事衛(wèi)星“天拓一號”[10]，隨后2015年發(fā)射“天拓三號”[11-12]、2017年發(fā)射“和德一號”[13]以及2018年和2020年發(fā)射的“海洋一號C”和“海洋一號D”[14]等多顆衛(wèi)星上搭載了AIS載荷，極大地豐富了我國的星載AIS數(shù)據(jù)庫。

近些年來，隨著衛(wèi)星AIS技術的快速發(fā)展，其應用領域越來越廣泛，并且逐漸持續(xù)朝向縱深層面發(fā)展。它可以提高船舶監(jiān)測能力，加強船舶交通監(jiān)管、打擊海上恐怖主義、協(xié)助海上搜尋救援。除上述基礎應用外，可將衛(wèi)星AIS數(shù)據(jù)與岸基AIS數(shù)據(jù)進行融合，有利于全面和完整地收集船舶AIS數(shù)據(jù)，實現(xiàn)對船舶的全程跟蹤與監(jiān)控。同時，可使電子海圖系統(tǒng)上顯示的船舶航跡等信息更加連貫，融合數(shù)據(jù)在利用過程中的效率和效益都將大幅提高。

隨著星載AIS數(shù)據(jù)庫的日益豐富，國內(nèi)在用AIS數(shù)據(jù)特別是星載AIS數(shù)據(jù)的研究方面開展了很多相關工作。例如，2014年姜洪濤和楊曉峰開發(fā)了一套可以在PC 平臺上運行的船舶實時監(jiān)控系統(tǒng)[15]；2015年肖瀟等利用軌跡聚類算法分析AIS 信息，以獲取船舶典型運動軌跡[16]；2016年李立春提出一種基于船舶自動識別系統(tǒng)(AIS)時態(tài)數(shù)據(jù)庫的船舶實時監(jiān)控和歷史軌跡查詢解決方案,有助于提高水上交通管理的效率并改善水上交通的安全[17]；2017～2018年李世友等利用天拓三號星載AIS的數(shù)據(jù)分析了全球船舶的檢測概率并得到了全球船舶的靜態(tài)分布[18]、分析了長程AIS數(shù)據(jù)的全球船舶分布和發(fā)展趨勢[19]，為下一代星載AIS的發(fā)展提供了必要的理論支持；2018年周田瑞等基于AIS數(shù)據(jù)獲取船舶航行典型運動軌跡，結(jié)合船舶歷史航速、航跡向等數(shù)據(jù)，預測船舶到港時間[20]；2019年劉奧博利用AIS衛(wèi)星獲取的船舶軌跡數(shù)據(jù)對2015年6月～2016年5月南極半島地區(qū)的人類活動進行分析[21]。

然而利用AIS數(shù)據(jù)對面向海面船舶越界方面的研究還少有涉及。本文結(jié)合星載AIS數(shù)據(jù)和海圖設計一款基于星載AIS的海面船舶越界預警系統(tǒng)，并用國內(nèi)自主研制的天拓三號[11-12]星載AIS的實測數(shù)據(jù)實現(xiàn)對系統(tǒng)的驗證。本系統(tǒng)有如下功能：(1) 船舶動態(tài)顯示自身以及附近船舶在海圖上的位置；(2) 結(jié)合海圖實現(xiàn)對船舶靠近領海、毗鄰區(qū)、專屬經(jīng)濟區(qū)域、國際公共海域等分界線進行預警。本文設計的船舶越界預警系統(tǒng)能極大地提高我國漁民海上捕撈作業(yè)的安全性，具有較好的實用性。

1 電子海圖與海洋分界線

1.1 地球表面任意兩點距離

精確計算地球表面上任意兩地理坐標點之間的距離對本系統(tǒng)具有重要的意義。常用的計算方法有球面余弦公式法、Haversine公式法等。根據(jù)球面余弦推導得到的Great-circle distance公式需要用到大量余弦函數(shù)的計算，因此在兩點之間距離過近時余弦函數(shù)會造成比較大的誤差；Haversine公式采用了正弦函數(shù)，因此在距離很小的情況下也能保持足夠的有效數(shù)字。本文采用Haversine公式來計算地球表面任意兩點距離

式中：d為待計算的兩點之間的距離；R為地球的半徑，取其平均值為6 371 km；φ1、φ2為兩點的緯度；Δλ則為兩點經(jīng)度差值。本節(jié)以北京到廣州兩地距離的計算為例驗證。北京和廣州的經(jīng)緯度坐標分別為(116.4°E, 39.9°N)和(113.15°E, 23.06°N)。利用Haversine公式得到了兩地相距1 889.175 km，該結(jié)果與高德地圖上的距離誤差在1 km以內(nèi)。

1.2 海洋分界線

注：其中A1～A4為插值點。圖1 船舶位置P到邊界線 AB的距離計算示意圖Fig.1 The calculation diagram of the distance between the ship position P and the boundary line AB

根據(jù)《聯(lián)合國海洋法公約》規(guī)定每一個沿海國家擁有的海域由內(nèi)海、領海、毗連區(qū)、專屬經(jīng)濟區(qū)四部分組成。根據(jù)公約規(guī)定領?；c以內(nèi)為內(nèi)海；領海區(qū)域是根據(jù)領?；c確定的，一般領海外部界線是根據(jù)領?；c線往外推12 n mile。領海相鄰的區(qū)域為毗連區(qū)(連接區(qū)、專門管制區(qū)等)，該區(qū)域是從領?；€往外推移24 n mile而劃分出的一個區(qū)域。而專屬經(jīng)濟區(qū)域則是不超過領?；c200 n mile；超過這個距離的區(qū)域?qū)儆诠！?/p>

由于現(xiàn)有的電子海圖上很少有對海洋分界線的劃分，本文通過各國已有的基點畫出基線，并通過基線推算出領海界線(距離基線12 n mile)、毗連區(qū)界線(距離基線24 n mile)、專屬經(jīng)濟區(qū)域(距離基線200 n mile)、公海(專屬經(jīng)濟區(qū)域之外)，然后把推算出的領?；€結(jié)合電子海圖顯示出來。本文設計的越界預警系統(tǒng)主要服務于中國沿海地區(qū)，因此本文首先利用中國政府公布的領?；c數(shù)據(jù)利用1.1節(jié)介紹的Haversine公式進行距離外推(外推算法略)從而得到海洋分界線，分別得到領海界線、毗連區(qū)界線和專屬經(jīng)濟區(qū)界線。

2 船舶越界判斷

2.1 船舶到邊界線距離

根據(jù)點到直線距離的計算方法，圖1給出某船舶P所在位置到邊界線AB最短距離的計算。實際計算過程，對AB之間進行數(shù)據(jù)插值的方式，使得AB之間的數(shù)據(jù)點足夠密集，然后利用Haversine公式來計算地球表面任意兩點距離的方法，計算得到PA、PA1、PA2、PA3、PA4和PB的距離，然后取其最小值作為P到邊界線AB的距離。

圖2 某船舶靠近邊界線越界判斷示意圖Fig.2 The judgment diagram of a ship′s crossing near the boundary line

2.2 船舶越界判斷

圖2給出某船舶在靠近某海洋界線的航行過程中的越界判斷示意圖。黑色船形為船體，點線為某邊界線，5個黑色箭頭代表5種可能的航行方向，分別用A, B, C, D 和E表示，對應的五個虛箭頭代表下一時刻的航行方向。當前時刻，船舶距離邊界線的最短距離為L0，下一時刻船舶距離邊界線的最短距離為L1。圖2中L1(a),L1(b),L1(c),L1(d)和L1(e)分別為5個可能航行方向上的下一時刻船舶所在位置與邊界線的最近距離。L1

3 船舶越界預警系統(tǒng)

3.1 系統(tǒng)要求

圖3 船舶客戶端的船舶越界預警 系統(tǒng)界面Fig.3 The interface of the ship boundary crossing warning system (For Vessel Client)

本文設計的越界預警系統(tǒng)要求在海圖上對海域進行符合國際法規(guī)范的海域劃分并顯示海域邊界線，實時解析接收到的星載AIS的數(shù)據(jù)并將船舶自身所在海域一定范圍之內(nèi)的船舶和軍艦的船舶標識、位置分布、動向以及船舶靠近各國海域邊境線警告?？紤]到便利使用，對客戶端的船舶越界預警系統(tǒng)的設計要求主要包括如下幾個方面的因素：(1) 數(shù)據(jù)選擇：實時接入衛(wèi)星數(shù)據(jù)或者用歷史數(shù)據(jù)模擬；(2) 海圖選擇；(3) 可選擇是否使用分界線顯示，或者顯示哪些分界線；(4) 越界距離選擇1 n mile、2 n mile，或根據(jù)實際需要輸入；(5) 預警信息提醒方式設置：界面上警示或者/語音警示。

3.2 總體設計

本文結(jié)合JAVA和MATLAB設計海面船舶越界預警系統(tǒng)。JAVA語言主要具有跨平臺、具有安全性、簡單等特點。MATLAB主要用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計算的高級技術計算機語言和交互式環(huán)境，尤其因為具備完善的數(shù)據(jù)可視化功能而備受軟件開發(fā)者的青睞。本文基于JAVA設計平臺對海面船舶越界預警系統(tǒng)進行界面設計，利用MATLAB實現(xiàn)數(shù)據(jù)計算和海圖表現(xiàn)，然后打包成jar包供Java調(diào)用。Java調(diào)用MATLAB中要注意以下幾點：(1) MATLAB的版本必須是2006b+(包括2006b或更高版本)，因為只有在這些版本中才有MATLAB Builder for Java(也叫Java Builder)；(2) 運行機器上必須裝有jre 并且其版本要和MATLAB自帶的jre版本保持一致(位數(shù)也要保持一致)；(3) 驗證系統(tǒng)jdk和jre。

圖4 海面船舶越界預警系統(tǒng) 工作流程圖Fig.4 The working flowing chart of the boundary crossing early warning system

設計的系統(tǒng)分為船舶界面端和監(jiān)管界面端，分別對應海面船舶用戶和監(jiān)管單位用戶，本文介紹的是面向海面船舶用戶的船舶客戶端。進入海面船舶客戶端的船舶越界預警系統(tǒng)界面如圖3所示。在該系統(tǒng)界面中用戶可以進行參數(shù)設置，包括數(shù)據(jù)選擇、海圖選擇、分界線顯示設置、越界警告距離設置、是否進行區(qū)域動態(tài)顯示以及語音播報設置等。數(shù)據(jù)選擇包括實時衛(wèi)星數(shù)據(jù)和歷史數(shù)據(jù)；其中，實時衛(wèi)星數(shù)據(jù)是默認選項。海圖選擇包括ECDIS海圖、SENC海圖等；需要說明的是，由于這些海圖獲取途徑比較困難，在本文中用MATLAB中默認的世界地圖作為演示。分界線顯示包括領?；€顯示、領海線顯示、毗鄰區(qū)線顯示、專屬經(jīng)濟區(qū)分界線顯示等，可以根據(jù)本船航行的位置進行設置。越界參數(shù)設置主要是警告距離設置，用戶可在界面上進行選擇1 n mile和2 n mile，也可以根據(jù)實際船舶速度大小進行自行輸入；一般選擇的越界警告距離是1 n mile。如果用戶選擇區(qū)域動態(tài)顯示，則在實際顯示中，除了用戶自身船舶之外，在區(qū)域內(nèi)的其他船舶也會實時顯示；否則只顯示本船。系統(tǒng)默認為區(qū)域動態(tài)顯示。在本系統(tǒng)中，語音播報暫未實現(xiàn)實際的功能，在界面設計的時候考慮到了語音播報。

3.3 軟件設計與工作流程

船舶越界預警系統(tǒng)的工作流程圖如圖4所示。系統(tǒng)啟動之后進行初始化，給出船舶自身的MMSI號，然后自動導入實時衛(wèi)星數(shù)據(jù)并在界面右側(cè)進行動態(tài)顯示。動態(tài)顯示系統(tǒng)通過對自身船舶經(jīng)緯度進行獲取，并默認設置海圖的初始顯示范圍為以本船為中心經(jīng)緯度正負20°的范圍。從接入數(shù)據(jù)到動態(tài)顯示的過程包括數(shù)據(jù)解析和導入海圖等過程。系統(tǒng)默認采用ECDIS海圖。天拓三號衛(wèi)星的數(shù)據(jù)解析過程參見文獻[11]。

圖5 海面船舶越界預警流程Fig.5 The boundary crossing early warning working flow

系統(tǒng)在工作過程中接受客戶端參數(shù)設置。如果船舶用戶導入歷史數(shù)據(jù)，則顯示一個月之內(nèi)本船的歷史航跡。船舶動態(tài)或者歷史航跡使用的海圖為默認的ECDIS海圖，用戶可以根據(jù)需要設置為其他的海圖。船舶用戶可根據(jù)自己的需要顯示分界線。例如，本船目前的工作地點在東海附近的專屬經(jīng)濟區(qū)之內(nèi)進行漁業(yè)作業(yè)，可顯示毗鄰區(qū)和專屬經(jīng)濟區(qū)的分界線，以便清楚自身的位置?？傊?，船舶根據(jù)自身的需要動態(tài)顯示在電子海圖上，形成一種直觀的表達方式，達到區(qū)域船舶動態(tài)的實時更新，對于船舶駕駛員能更好地了解船舶周圍情況，及時做出應對措施。

船舶在海面作業(yè)時如果靠近分界線，特別是靠近專屬經(jīng)濟區(qū)的分界線，則很容易發(fā)生一些國際漁業(yè)糾紛，甚至導致被鄰國扣押。因此，越界預警非常重要。本系統(tǒng)最重要的組成部分是越界預警模塊。該模塊通過星載AIS船舶反饋的數(shù)據(jù)在電子海圖上顯示，能夠在海圖上顯示出本船舶的基本信息(船艏向、航行速度、船舶位置等)，并且可以通過給定的航線判斷船舶是否偏離航向；也能更加方便本船獲悉所處位置以及船舶距離領海、專屬經(jīng)濟區(qū)域等界線的距離，能夠有足夠的時間對船舶進行航向調(diào)整。

本系統(tǒng)的越界預警模塊中默認的預警距離是1 n mile；船舶用于根據(jù)實際需要選擇2 n mile或者自行設置越界預警距離。船舶越界預警模塊的預警流程見圖5所示。首先，船舶獲取越界設置中越界距離D0并獲取邊界線的數(shù)據(jù)；然后，根據(jù)1.1節(jié)中的Haversine公式計算本船位置與邊界線的距離，取其最小值L0，并根據(jù)當前位置和航向計算下一時刻的位置，在此基礎之上計算下一個時刻與邊界線距離的最小值(L1)。根據(jù)D0、L0和L1三個值來判斷本船是否即將越界。越界預警條件為：L0

4 基于星載AIS的船舶越界系統(tǒng)驗證

4.1 星載AIS系統(tǒng)與數(shù)據(jù)

船舶越界預警系統(tǒng)主要通過船載AIS設備實時向衛(wèi)星發(fā)送AIS信息，經(jīng)過衛(wèi)星星載AIS接收設備實時獲取本船和周邊海域相關船只的AIS信息進行數(shù)據(jù)分析，進而得到相關的預警信息。本系統(tǒng)主要依賴“天拓三號”衛(wèi)星星載AIS進行數(shù)據(jù)接收和傳輸，所有數(shù)據(jù)均使用的是“天拓三號”衛(wèi)星所接收到的船舶信息[11-12]。星載AIS接收到的每一條報文信息主要有發(fā)送時間、數(shù)據(jù)消息識別碼、用戶識別碼、航向狀態(tài)、對地航速、對地航向、經(jīng)緯度、通信狀態(tài)等信息等幾部分組成。特別說明每一條報文的發(fā)送時間由前面的時間信標(精確到分鐘)和報文中給出的秒值共同合成[12]。

圖6 2016年3月13日MMSI為503****的船舶穿越領海 基線進入領海并將越過領海分界線示意圖Fig.6 Schematic diagram of ship crossing the baseline and entering the territorial water (2016-03-13)

4.2 越界預警案例一

本節(jié)根據(jù)“天拓三號”衛(wèi)星星載AIS接收到報文信息解析結(jié)果進行越界預警案例篩選。選擇的基本原則是，某船舶在某海域的分界線(領海線、毗鄰區(qū)界線或者專屬經(jīng)濟分界線)附近存在至少3個AIS偵收數(shù)據(jù)點。圖6為船舶號MMSI為503****的一艘漁船在2016年3月13日的一段航行動態(tài)圖。需要注意的是此時在參數(shù)設置的時候沒有進行區(qū)域動態(tài)顯示。本文不從地圖規(guī)范方面考慮，因此在圖6中將具體的地圖數(shù)據(jù)隱去，經(jīng)緯度坐標為相對經(jīng)緯度。方點線、十字點線和星點線分別代表領?；€、領海分界線和毗鄰區(qū)分界線。該漁船在2016年3月13日從福建附近海域穿越領?；€進入領海，并即將穿越領海分界線進入毗鄰區(qū)。該漁船在8:40時距離領海分界線2 n mile，此時收到越界預警：本船距離領海分界線2 n mile，預計6.5 min后進入毗鄰區(qū)。

如果本系統(tǒng)運行在監(jiān)管界面端，從監(jiān)管的角度來看，如該類船舶屬于只能在領海航行的船舶應該對該船警告并提醒該船舶返航，如該船舶具備在領海外區(qū)域航行的資格，則提醒該船舶已經(jīng)靠近領海區(qū)域。

4.3 其他越界預警案例

本小節(jié)簡要闡述另外2個船舶越界案例。2018年5月13日08:50 UT左右天拓三號衛(wèi)星星載AIS偵收到在中國沿海航行的MMSI號為503****的船舶動態(tài)。該船在2018年5月13日8:48時分距離領?；€連接線約3 km，再繼續(xù)往前航行時，在距離界線2 n mile時船舶給出進入領海的越界通知。由于該船舶航行速度較快，該船舶于8:53時候已經(jīng)航行在領海基線之外。

2017年10月14日06:52UT左右天拓三號衛(wèi)星監(jiān)測到MMSI號為11111****的船舶的航行動態(tài)。該船在東海海域位于中國專屬經(jīng)濟區(qū)和日本專屬經(jīng)濟區(qū)的邊界線附近，距離邊界線僅僅為1 n mile，且其航行方向為向中國專屬經(jīng)濟區(qū)一側(cè)。在本文的模擬應用系統(tǒng)中，該船舶在東海海域作業(yè)，在距離專屬經(jīng)濟區(qū)界線2 n mile即將越過中國專屬經(jīng)濟區(qū)邊界時收到即將越界進入鄰國專屬經(jīng)濟區(qū)的即將越界預警信息，并且，在靠近邊界線時不斷收到預警。在即將越過界線并收到越界預警之后，立即調(diào)整航向返回中國專屬經(jīng)濟區(qū)，從而避免越界。

5 結(jié)論

本文結(jié)合星載AIS數(shù)據(jù)和海圖設計一款基于星載AIS的海面船舶越界預警系統(tǒng)，并用國內(nèi)自主研制的天拓三號星載AIS的實測數(shù)據(jù)實現(xiàn)對系統(tǒng)的驗證。該系統(tǒng)有如下功能：

(1) 船舶動態(tài)顯示自身以及附近船舶在海圖上的位置；

(2) 結(jié)合海圖實現(xiàn)對船舶靠近領海、專屬經(jīng)濟區(qū)域、國際公共海域等分界線進行預警。本越界預警系統(tǒng)在海面船舶預警中具備以下幾點意義：

①減少漁船發(fā)生越界捕撈的概率；

②為漁船作業(yè)提供安全保障；

③提高海事管理機構對于海洋監(jiān)控的實時性。

該系統(tǒng)的應用能極大地提高海面船舶的航行安全，尤其是我國漁民海上捕撈作業(yè)，具有極大的實用性。下一步可以安裝在海洋或漁政相關監(jiān)管部門作為系統(tǒng)監(jiān)控之用。監(jiān)管單位只需要將海圖范圍設定在自己的監(jiān)管區(qū)域之內(nèi)，利用實時接收的星載AIS信號解析得到本系統(tǒng)內(nèi)的相關船舶動態(tài)圖，然后對可能的越界行為發(fā)出預警即可。

本系統(tǒng)目前只設計了本船的越界預警功能。未來根據(jù)需要，可擴展對友鄰船舶的越界預警功能和對區(qū)域船舶的越界預警等功能?！坝燕彺浇珙A警”可對一同出海的友鄰船舶是否越界提出預警提示，以便在特殊情況下能夠及時對友鄰船舶進行必要的援助；“區(qū)域船舶越界預警”主要是為了方便了解船主所關心的區(qū)域船舶是否有越界行為，也可以方便海上管理機構對于每艘船舶的管理與調(diào)度。