摘 要 考慮到影響船舶吃水的主要因素：船舶對波浪的響應(yīng)、橫傾引起的吃水變化、淺水中船舶的下沉和首傾以及航道底部情況分布特征是非常復(fù)雜和隨機的，鑒于對動態(tài)吃水引起可航水深不足而導(dǎo)致擱淺的研究相對偏少，本文主要圍繞船舶動態(tài)吃水采用蒙特卡洛方法建立船舶擱淺概率計算模型，為船舶安全通過受限水域提供依據(jù)，具有科學(xué)和現(xiàn)實應(yīng)用的雙重意義。

關(guān)鍵詞 動態(tài)吃水 蒙特卡洛方法 擱淺概率

中圖分類號：U661 文獻標(biāo)識碼：A

The Probability of Ship Grounding Based on Monte Carlo Method

LIU Chao

（Anhui Communications Vocational Technical College， Hefei， Anhui 230051）

Abstract Considering the main factors affecting the ship draft： response to wave， ship heeling induced draught， sinking ships in shallow water and the first bend and channel bottom case distribution is very complex and random， in view of the dynamic draft caused by lack of navigable depth and led to the stranding of relatively low， the main focus of this article the ship dynamic draft a probability model using Monte Carlo method， provide the basis for the safety of the ship in restricted waters， which has dual significance of science and reality application.

Key words Dynamic draft； Monte Carlo method； probability of ship grounding

隨著船舶的大型化和高速化，船舶擱淺的問題日益凸顯。擱淺事故常帶來船體結(jié)構(gòu)損壞、環(huán)境污染甚至人員傷亡等嚴(yán)重后果，擱淺是個非常復(fù)雜的科學(xué)與工程技術(shù)問題，但其發(fā)生的必要條件比較簡單，即航行水域水深相對船舶吃水不足?？紤]到影響船舶吃水的主要因素：船舶對波浪的響應(yīng)、橫傾引起的吃水變化、淺水中船舶的下沉和首傾以及航道底部情況分布特征是非常復(fù)雜和隨機的，鑒于對動態(tài)吃水引起可航水深不足而導(dǎo)致擱淺的研究相對偏少，本文主要圍繞船舶動態(tài)吃水采用蒙特卡洛方法建立船舶擱淺概率計算模型，為船舶安全通過受限水域提供依據(jù)，具有科學(xué)和現(xiàn)實應(yīng)用的雙重意義。

1 國內(nèi)外研究現(xiàn)狀

自上世紀(jì)50年代末以來國內(nèi)外學(xué)者對船舶擱淺進行了一系列研究。這些研究成果大概可歸納為三類：第一類主要從船舶實際操作入手開展研究，如采用實船事故調(diào)查、事件樹、模糊事件樹、決策樹、貝葉斯網(wǎng)絡(luò)等方法，探索擱淺的原因和機理，對船舶偏離航道、機械失效、救援失效、惡劣海況等導(dǎo)致擱淺的因素進行分析，定性或定量地描述船舶擱淺風(fēng)險的大小；①第二類主要從船體結(jié)構(gòu)力學(xué)入手研究，采用縮尺或?qū)嵆叨饶Ｐ驮囼?、有限元分析、簡化?shù)值模擬和解析經(jīng)驗公式等方法，研究擱淺事故中船舶力學(xué)問題，分析船底結(jié)構(gòu)損傷、變形受力，為改善船體結(jié)構(gòu)提供參考；②第三類則側(cè)重于航道中船舶通航條件等相關(guān)數(shù)據(jù)的采集和分析，以充分利用航道水深，減小船舶擱淺風(fēng)險。黃蘊和③針對重載船舶在進出浙江海門港乘潮作業(yè)問題，結(jié)合港外航道潮汐變化規(guī)律，研討了現(xiàn)有航道水深條件下乘潮作業(yè)的合理性。胡勤友④綜合潮汐數(shù)據(jù)、電子海圖、船位和航向信息等開發(fā)了船用擱淺預(yù)警系統(tǒng)。

蒙特卡洛（Monte Carlo）仿真是一種研究船舶擱淺概率較新的方法。它在考慮影響船舶吃水各因素的基礎(chǔ)上，根據(jù)其分布規(guī)律產(chǎn)生大量隨機數(shù)，最后統(tǒng)計得到擱淺的頻率，進而得出概率。N.W.A. Schoenmakers采用該方法對開普敦港（Cape Town）進出港船舶個擱淺概率進行了研究，波浪處理采用了swan（近岸海浪數(shù)值）模型，船舶對波浪的響應(yīng)采用了簡化的經(jīng)驗公式，而且沒有考慮船舶在波浪中的橫傾和縱傾問題。Gucma采用了類似的方法，更偏向于經(jīng)驗型，相關(guān)數(shù)據(jù)的獲取采用了模擬器與實船實驗相結(jié)合的方法，另外又考慮了可航的泥底、各種測量誤差等因素，求出擱淺的頻率。但受限于試驗次數(shù)和實驗條件的設(shè)置，誤差等參數(shù)選取進行了大量簡化。為此，本文擬根據(jù)蒙特卡洛仿真的思想，結(jié)合CFD（計算流體動力學(xué)）計算結(jié)果，開展設(shè)計代表船舶在航道中擱淺的概率計算。

2 蒙特卡洛方法概述

蒙特卡羅（Monte Carlo）方法，又稱計算機隨機模擬方法，是一種基于“隨機數(shù)”的計算方法。Monte Carlo 方法以概率統(tǒng)計理論為基礎(chǔ)，以隨機抽樣（隨機變量的抽樣）為手段，在很多方面有重要的應(yīng)用。它的優(yōu)點表現(xiàn)在三個方面：方法和程序的結(jié)構(gòu)簡單，易分析、易理解；收斂的概率性和收斂速度與問題的維數(shù)無關(guān)，很好地避免了維數(shù)問題；受問題條件限制的影響較小，很好地提高可行性。

3 建立船舶擱淺概率計算模型

影響船舶動態(tài)吃水的主要因素包括：潮汐、船舶對波浪的響應(yīng)、淺水中船舶下坐、船舶操縱運動所致浮態(tài)變化以及航道底部情況等?？紤]到船舶動態(tài)吃水采用蒙特卡洛方法建立船舶擱淺概率計算模型。

3.1 獲取影響船舶吃水各因素的分布特征函數(shù)

影響船舶動態(tài)吃水的因素如：潮汐、船舶對波浪的響應(yīng)、淺水中船舶下坐、船舶操縱運動所致浮態(tài)變化以及航道底部情況等。具有很強的隨機性和不確定性，但也有其分布規(guī)律。因此我們可統(tǒng)計其分布規(guī)律，構(gòu)建概率密度函數(shù)再基于蒙特卡洛產(chǎn)生大量隨機數(shù)，進行多次仿真。完成船舶對波浪的響應(yīng)，可以采用簡化的波浪模型，簡化航道邊界條件，選取合理波浪響應(yīng)模型。對于橫傾引起的吃水變化，可以通過統(tǒng)計記錄船舶橫傾角的分布規(guī)律，然后根據(jù)給定船型的基本數(shù)據(jù)計算得到橫傾引起的吃水變化的規(guī)律。至于淺水中船舶的下沉和首傾，從水動力學(xué)角度使用三維面元法可以獲取船舶在淺水航道的船體下沉和首傾規(guī)律。通過航道水深基本數(shù)據(jù)，采用計算機擬合的方法，來描述航道底部情況。

3.2 建立船舶擱淺概率計算模型

通過上一步研究得來的船舶對波浪的響應(yīng)、橫傾引起的吃水變化、淺水中船舶的下沉和首傾以及航道底部情況的分布特征函數(shù)，采用蒙特卡洛方法產(chǎn)生一系列隨機數(shù)，根據(jù)下列公式求取給定航道受限水深的船舶擱淺概率P（Z<0）。

Z=d+ Htide -（D+r+Bf+h+S）（d：航道深度（m）；Htide：潮高（m）；D：設(shè)計船型吃水（m）；r：船舶對波浪的垂向響應(yīng)（m）；Bf：航道底部情況（m）；h：橫傾引起的吃水變化（m）；s：淺水引起的下坐（squat）（m）；）

因此根據(jù)蒙特卡洛方法求取船舶擱淺概率，即P（Z<0），可以通過研究上式中各因素的概率密度函數(shù)，產(chǎn)生大量的隨機數(shù)序列，最后統(tǒng)計船舶擱淺的頻率以代替擱淺概率。

3.3 對擱淺發(fā)生的后果進行量化評估，建立風(fēng)險模型

風(fēng)險（Risk）要從危險發(fā)生的概率（P， Probability of danger）和危險發(fā)生的后果（C， Consequence of danger）兩個方面同時確定，即：Risk=P*C。第二步已經(jīng)求取了擱淺發(fā)生的概率，因此接下來要用統(tǒng)計的方法，確定擱淺發(fā)生的后果，從而建立風(fēng)險模型，描述風(fēng)險隨航道水深變化的分布規(guī)律。

4 結(jié)束語

本文通過對影響船舶吃水的主要因素：船舶對波浪的響應(yīng)、橫傾引起的吃水變化、淺水中船舶的下沉和首傾以及航道底部情況的分布規(guī)律進行研究，考慮到影響因素的復(fù)雜性和隨機性，使用蒙特卡洛方法構(gòu)建出船舶擱淺概率計算模型，是具有科學(xué)和現(xiàn)實應(yīng)用的雙重意義。

注釋

① 孫小林.長江航道新治：5年內(nèi)投1500億建成黃金水道[J].21世紀(jì)經(jīng)濟報道，2011.

② 李文超.港口航道通過能力研究綜述[J].民晉科技，2012.11.

③ 黃蘊和，金丕信.海門港外航道通航條件研究[J].中國航海，1996.1：77-83.

④ 龔安祥，胡勤友，徐鐵.SAGA：一種新的船舶擱淺預(yù)警模型[J].大連海事大學(xué)學(xué)報，2007.28（1）：106-110.