海上落水人員漂移軌跡可視化預測研究

夏山宏 徐純潔 羅修波 謝方偉 冷江南 田祖織

（1.黑龍江省水上交通救援中心，黑龍江哈爾濱 1 500119；2.中國礦業(yè)大學 機電工程學院，江蘇徐州 221116；3.泰州市金海運船用設備有限責任公司，江蘇泰州 225300）

0 引 言

近年來，我國海洋經(jīng)濟持續(xù)發(fā)展壯大，海上運輸貿易蓬勃發(fā)展，海上航行船舶與日俱增，隨之而來的是海上險情逐年增加、海難事故頻發(fā)[1-2]。海洋事故發(fā)生后的海洋搜救以海上搜尋和海上救援為主，搜尋是海上救援的前提，快速、準確、動態(tài)地預測海難事故中落水目標及其漂移位置，將大大縮小搜救范圍、提高海上救援效率，因此，如何高效地開展海上搜尋成為海洋搜救的當務之急[3-4]。

針對海上搜尋，美國、英國和法國等西方發(fā)達國家較早地開展了相關研究，開發(fā)出了搜救目標漂移預測的計算機輔助決策系統(tǒng)，如POSSE、SAR IS 和MOTHY 等系統(tǒng)[5]。然而上述輔助決策系統(tǒng)在預測漂移軌跡時僅僅考慮了海風與海流的作用，預測結果與實際往往還存在較大差距。

為進一步提升海上落水目標漂移軌跡預測的精度，國內外學者對此開展了大量研究。劉同木等[6]提出了一種基于受力分析的落水人員漂移軌跡方法，并通過海上仿真人漂移實驗驗證了其可靠性；劉柯等[7]采用Leeway 模型和拉格朗日追蹤法對內河流道漂移物體的軌跡進行預測，提出了結合航道邊界特征的最終漂移位置預測方法；王茹意等[8]在利用蒙特卡羅法預測海上落水遇險人員漂移軌跡時，選取分段線性插值法來表述流場數(shù)據(jù)，獲得了更佳的預測效果；Zhu 等[9-10]基于大量的觀測數(shù)據(jù)和多種數(shù)值模擬方法對海上目標漂移運動模型和海上目標搜尋區(qū)域計算模型進行了系統(tǒng)研究，并提出了相應的模型改進和算法。

本文在上述研究基礎之上，開展了海上落水人員漂移軌跡可視化預測研究，基于海上落水人員受力數(shù)學模型分別使用解析法和蒙特卡洛模擬法仿真得到落水人員漂流軌跡，并進行對比；在Matlab AppDesigner 軟件上引入經(jīng)緯度信息和電子海圖搭建了海上落水人員漂移軌跡精準可視化系統(tǒng)，實現(xiàn)了更高效更準確地海上落水人員漂移預測。

1 海上落水人員漂移軌跡預測模型

海上落水人員在漂移的過程中，主要受到水上部分和水下部分兩方面的影響，其中，水上部分主要受到海流和波浪的影響，水下部分則主要受到海面風的影響，另外海上落水人員漂移還受到科氏力的影響，因此，海上落水人員的漂移速度方程可以表示為：

根據(jù)牛頓第二定律，海上落水人員漂移軌跡的運動方程也可表示為：

式中：m 為落水人員質量；V 為漂移運動速度；t 為時間；Fc為科氏力；Fwind為風作用力；Fcurrent為海流作用力；Fwave為海浪作用力。

以東西方向x 軸，南北方向為y 軸建立平面直角坐標系，其中正東方向為x 軸正方向，正北向為y 軸正方向。落水人員在水平和垂向速度分別設為u 和v，將海風作用力、海浪作用力、海流作用力和科氏力分別分解到水平和垂直兩個方向，則上式可變換為：

2 海上落水人員漂移軌跡預測仿真參數(shù)設置

海上落水人員漂移軌跡與自身質量、海風和海流相關，本文落水人員質量取70 kg，海風作用橫截面積取0.4 m2，其中風速取5.2 m/s，風向角取15°，海流作用橫截面積取1.0 m2，其中流速取0.5 m/s，流向角取10°，波高平均取值為1.6 m，浪向角取10°。由于海風海流處于實時變化之中，仿真模擬時設置一定范圍，具體參數(shù)見表1。

表1 模擬海上落水人員漂移軌跡具體參數(shù)

3 海上落水人員漂移軌跡預測對比

目前，海上落水人員漂移軌跡預測方法主要有解析法、蒙特卡洛模擬法、拉格朗日粒子追蹤法等，為了更好地驗證模型的有效性，本文選擇解析法和蒙特卡洛模擬法兩種方法進行海上落水人員漂移軌跡預測的相互印證，漂移軌跡預測和印證流程如圖1 所示。

圖1 漂移軌跡預測和對比流程

解析法是通過分析海上落水人員受力情況，得到落水人員漂流動力學方程，通過解析動力學方程計算出落水人員漂移軌跡的方法。從理論上講，只要能夠綜合考慮到所有受力情況，解析法得到的落水人員漂流軌跡是準確無誤的。但是由于海上環(huán)境復雜，落水人員自身和海上環(huán)境一直處于動態(tài)變化之中，在建立動力學方程時會忽略一些因素，因此會產(chǎn)生一定的誤差。本文考慮落水人員自身質量、海風作用力、海流作用力和科氏力綜合作用，基本吻合海上落水人員受力情況，在Matlab 軟件中通過求解公式（3）和公式（4）可以得到落水人員的漂移軌跡。

蒙特卡洛模擬法是隨著科學技術和電子計算機技術發(fā)展產(chǎn)生的數(shù)值模擬統(tǒng)計方法，通過概率統(tǒng)計的方法來計算某種事件發(fā)生的概率。海上落水人員漂移軌跡模擬中，通過擬合大量抽象粒子的位置來預測落水人員漂移曲線，仿真中抽象出的粒子數(shù)越多，仿真結果越接近實際值。蒙特卡洛模擬法以公式（1）和公式（2）為基礎，仿真模擬落水人員的漂流軌跡。

圖2 所示為解析法和蒙特卡洛法兩種方法得到的海上落水人員的漂移軌跡對比，紅色軌跡表示解析法得到的漂移軌跡，黑色軌跡表示蒙特卡洛模擬法得到的預測漂移軌跡。由圖2 的漂移軌跡對比圖可以看出，解析法和蒙特卡洛模擬法仿真得到的海上落水人員漂移軌跡初始階段保持一致。隨著漂移時間增加，兩者軌跡出現(xiàn)一定偏差，但是漂移軌跡基本保持吻合，兩種漂移軌跡預測方法實現(xiàn)了互相印證，模型的有效性得到了印證。但是表示落水人員位置時采用了直角坐標系南北坐標和東西坐標來反映落水人員位置，這與實際海上利用經(jīng)緯度信息搜尋人員不符合，為了更準確更直觀地反映出海上落水人員的漂移軌跡，更好地服務海上搜救，需對模型進行進一步的優(yōu)化，下一章節(jié)將以解析法為基礎進一步研究海上落水人員漂移軌跡精準可視化預測。

圖2 兩種方法漂移軌跡對比圖

4 MATLAB 模擬仿真與可視化預測

解析法建立的模型中引入直角坐標系來表示海上落水人員位置，這與海上實際搜救情況和北斗定位系統(tǒng)顯示的數(shù)據(jù)不符合。因此，在運用Matlab AppDesigner 軟件搭建的海上落水人員漂移軌跡可視化預測系統(tǒng)中引入了經(jīng)緯度信息，將直角坐標轉化為經(jīng)緯度坐標信息，為了更方便搜救人員定位和搜救指揮中心調度搜救力量，同時還引入電子海圖來顯示落水人員漂移軌跡。

海上落水人員漂移軌跡可視化預測系統(tǒng)主程序界面、漂移軌跡數(shù)據(jù)界面和海況數(shù)據(jù)界面三個板塊組成。其中主程序界面主要分為三個區(qū)域，如圖3所示：第一個區(qū)域具備參數(shù)錄入功能，救援人員可以根據(jù)實際海況數(shù)據(jù)和遇險人員具體信息將落水人員的關鍵參數(shù)、海況參數(shù)、落水時間和落水位置輸入系統(tǒng)；第二個區(qū)域具備查詢漂移軌跡和落水人員位置功能，救援人員可以根據(jù)落水人員的落水時長來獲取落水人員漂移各個時刻的的經(jīng)緯度位置信息，以便于快速定位救援地點；第三個區(qū)域具備漂移軌跡顯示功能，引入的電子海圖可以更清晰顯示遇險人員落水區(qū)域，更方便搜救指揮中心調度搜救力量，提高搜救成功率。

圖3 海上落水人員漂移軌跡可視化預測系統(tǒng)

圖4 所示為系統(tǒng)運行界面顯示，紅色曲線代表落水人員的漂移預測軌跡，根據(jù)系統(tǒng)顯示可以準確地獲取落水人員落水點、飄移終點和漂移軌跡曲線，也可以移動星標位置來查看某位置處的時間信息。

圖4 海上落水人員漂移軌跡可視化預測

同時，系統(tǒng)的漂移軌跡數(shù)據(jù)和海況數(shù)據(jù)兩個板塊可以實時記錄信息并存儲導出，為后續(xù)海上搜救行動總結經(jīng)驗和后續(xù)優(yōu)化提供支撐材料，其界面顯示結果如圖5 和圖6 所示。

圖5 漂移軌跡數(shù)據(jù)和海況數(shù)據(jù)記錄界面

圖6 海況數(shù)據(jù)和海況數(shù)據(jù)記錄界面

5 結 論

本文通過分析海上落水人員受力信息，建立了海上落水人員漂移軌跡數(shù)學模型，分別基于解析法和蒙特卡洛模擬法得到了海上落水人員的漂移軌跡，對比結果表明：兩種方法得到的漂移軌跡基本一致，驗證了模型的有效性。在Matlab AppDesigner 軟件上引入經(jīng)緯度坐標和電子海圖搭建了海上落水人員漂移軌跡可視化預測系統(tǒng)，可以實現(xiàn)海上落水人員漂移軌跡的精準可視化預測，大大提高海上搜救遇險落水人員和指揮中心調度海上救援力量的效率，研究內容和結果對海上搜救具備一定參考指導意義。