范麗先， 葉圓慧， 尹靜波

(1. 上海大學(xué) 管理學(xué)院， 上海 200444； 2. 上海交通大學(xué) 船舶海洋與建筑工程學(xué)院， 上海 200240)

據(jù)聯(lián)合國貿(mào)易和發(fā)展會議(United Nations Conference on Trade And Development, UNCTAD)[1]統(tǒng)計，當(dāng)前全球海上貿(mào)易運(yùn)輸量占國際貿(mào)易運(yùn)輸總量的80%左右。我國作為出口貿(mào)易大國，海上貿(mào)易運(yùn)輸量占我國國際貿(mào)易運(yùn)輸總量的比例更高達(dá)90%以上。[2]隨著船舶建造技術(shù)和管理水平不斷提高，海上運(yùn)輸能以較低的成本完成更加安全、可靠和高效的運(yùn)輸任務(wù)，從而有效加快世界經(jīng)濟(jì)的全球化進(jìn)程。[3]

但是，海運(yùn)事業(yè)的蓬勃發(fā)展也給海上交通安全帶來很大隱患。根據(jù)國際海事組織(International Marithime Organization, IMO) 官方網(wǎng)站統(tǒng)計，從1986年初至2015年底，前10 a共發(fā)生海上事故1 257起，中間10 a共發(fā)生4 394起，近10 a下降到3 648起。海上交通事故發(fā)生率雖然遠(yuǎn)不及陸路交通事故發(fā)生率，但一旦發(fā)生，就會造成巨大的人員傷亡、財產(chǎn)損失和環(huán)境污染。[4]

海上事故發(fā)生后付出的沉重代價使人們不得不從中汲取教訓(xùn)，以盡量避免事故再次發(fā)生。在眾多研究中，貝葉斯網(wǎng)絡(luò)(Bayesion Network, BN)技術(shù)因具有系統(tǒng)的成本效益分析能力而廣泛應(yīng)用于對風(fēng)險控制手段的預(yù)測和評估中。[5]HNNINEN等[6-7]和 BANDA等[8]用貝葉斯模型分析海上安全管理與影響變量之間的關(guān)系，通過研究這些變量之間的關(guān)聯(lián)和變量與整體安全管理水平之間的關(guān)聯(lián)為安全管理提供決策支持。周國華等[9]用貝葉斯模型分析影響京滬高速鐵路項(xiàng)目質(zhì)量和風(fēng)險的關(guān)鍵因素，以進(jìn)行項(xiàng)目風(fēng)險管理。這些研究通常采用專家的意見建立貝葉斯網(wǎng)絡(luò)模型，具有較強(qiáng)的主觀性。為降低主觀性對貝葉斯模型的影響，不少學(xué)者采用客觀的船舶事故數(shù)據(jù)分析風(fēng)險概率，例如：YIN等[10]采用船齡、船型、船旗及船級社等指標(biāo)，通過Logit 模型擬合海上交通事故發(fā)生的概率，并在此基礎(chǔ)上建立一套海上交通安全指標(biāo)系統(tǒng)；LI等[11]用一種創(chuàng)新的方式將Logistic回歸與貝葉斯網(wǎng)絡(luò)相結(jié)合應(yīng)用到海上風(fēng)險分析中，研究各變量在不同狀態(tài)下對船舶損失和維修保養(yǎng)成本的影響。這種將客觀的條件概率與貝葉斯網(wǎng)絡(luò)相結(jié)合的方法可使航運(yùn)事故及其成本效益分析更加準(zhǔn)確、客觀，但只選取幾個基本的變量，僅將事故及系統(tǒng)因素狀態(tài)劃分為正常和失效2種狀態(tài)。本文對IMO和巴黎諒解備忘錄(Memoradum of Understanding, MoU)組織公布的船舶事故及事故影響程度進(jìn)行詳細(xì)分析，并將其影響因素(如檢查次數(shù))擴(kuò)展至多態(tài)。

1 數(shù)據(jù)來源及研究方法

1.1 數(shù)據(jù)來源

本文構(gòu)造貝葉斯模型的數(shù)據(jù)來自于以下4個方面。

1) 船舶基本信息：包括近12萬艘船舶(100總噸以上)的200多項(xiàng)基本信息(Lloyd’s全球船舶數(shù)據(jù)庫)。

2) 海上事故信息：IMO網(wǎng)站公布的2006—2015年3 648項(xiàng)海上事故報告。

3) 海上安全檢查報告：巴黎MoU自2006—2015年發(fā)布的有關(guān)港口國監(jiān)控(Port State Control, PSC)年度報告，包含20余萬條檢查結(jié)果。

4) 世界經(jīng)濟(jì)發(fā)展數(shù)據(jù)(World Bank及Clarkson數(shù)據(jù)庫)：所有國家近10 a的經(jīng)濟(jì)、貿(mào)易、人口及資源等數(shù)據(jù)。

1.2 貝葉斯網(wǎng)絡(luò)模型

貝葉斯網(wǎng)絡(luò)是一種概率圖模型，其網(wǎng)絡(luò)拓樸結(jié)構(gòu)是一個有向無環(huán)圖,由節(jié)點(diǎn)和帶方向的箭頭組成，節(jié)點(diǎn)表示隨機(jī)變量，帶方向的箭頭表示因果依賴關(guān)系。

根據(jù)文獻(xiàn)[10]、文獻(xiàn)[12]及巴黎MoU網(wǎng)站相關(guān)信息選擇6個變量，分別為船舶尺度(Size)、船齡(Age)、船旗(Flag)、船型(Ship Type)、PSC檢查次數(shù)(Inspection)和滯留情況(Detention)。前5個變量直接影響事故發(fā)生率(Accident)，PSC檢查次數(shù)直接影響滯留情況。圖1為船舶事故的貝葉斯模型。船舶尺度分為大船和小船；船齡分為新船和舊船；船旗分開放式和封閉式；PSC檢查次數(shù)沿襲CARIOU等[13]的分類方法，分為未被檢查過、檢查次數(shù)在1～10次和大于10次等3種；滯留情況分為曾滯留和未滯留過2種。船舶類型分為干散貨船(Bulker)、集裝箱船(Container)、油船(Tanker)、客船(Passenger)和一般貨船(General Cargo)。此外，YIN等[10]發(fā)現(xiàn)發(fā)達(dá)國家的事故率比發(fā)展中國家更低，因此對發(fā)達(dá)國家和發(fā)展中國家航運(yùn)安全因素分別進(jìn)行研究。

1.3 貝葉斯網(wǎng)絡(luò)模型的先驗(yàn)概率

船舶事故分為嚴(yán)重事故和較輕事故，因此對嚴(yán)重事故船舶和較輕事故船舶分別進(jìn)行建模分析。

利用建立的數(shù)據(jù)庫計算圖1中貝葉斯網(wǎng)絡(luò)模型的各先驗(yàn)概率，其中船舶尺度、船齡、船舶類型和船旗等4個變量均為外生自變量，采用各狀態(tài)所占比例作為其先驗(yàn)概率(見圖2和圖3)。

注：事故發(fā)生率、PSC檢查次數(shù)和滯留情況各概率值為其各狀態(tài)下條件概率期望值

圖3 發(fā)展中國家船舶事故先驗(yàn)概率

1.4 基于Logistic回歸模型的條件概率

為估計圖2和圖3中的貝葉斯模型，需進(jìn)一步計算各影響因素不同狀態(tài)下的船舶事故、檢查次數(shù)和滯留情況的條件概率。采用Logistic回歸方法估計各因素對事故、檢查和滯留情況的影響系數(shù)，進(jìn)而

對每艘船舶發(fā)生事故、被PSC檢查和滯留的概率進(jìn)行預(yù)測，最后根據(jù)圖1中各影響因素的不同狀態(tài)計算條件概率。

表1為各因素對事故發(fā)生率、PSC檢查次數(shù)和滯留情況的Logistic回歸分析結(jié)果，其中：模型1是各因素對事故的影響系數(shù)；模型2是各因素對檢查情況的回歸系數(shù)；模型3是各因素對船舶滯留情況的影響系數(shù)。3個模型的極大似然比均高度顯著，表明模型的整體擬合度均較高。由于國家類型、船旗和船舶類型均為類別變量，因此表1中該3個變量以啞變量形式代入回歸，即Develop代表發(fā)達(dá)國家，Open代表開放式登記旗船，Bulker，Container，Tanker，Passenger，General Cargo分別代表5種主要類型船舶啞變量，未放入回歸中的對比組為其他類型船舶。

表1 各因素對事故發(fā)生率、PSC檢查次數(shù)和滯留情況的Logistic回歸分析結(jié)果

2 結(jié)果分析

2.1 各因素對嚴(yán)重事故的影響及敏感性分析

模型結(jié)果顯示，各因素對嚴(yán)重事故和較輕事故的影響均類似，因此僅列出對嚴(yán)重事故影響的結(jié)果分析。表2～表7中：最后一列為不同狀態(tài)下后驗(yàn)概率的變化率；最后一行為發(fā)達(dá)國家對比發(fā)展中國家的變化率。

表2為船舶尺度對航運(yùn)安全影響的敏感性分析。最后一列為不同船舶尺度下事故后驗(yàn)概率的變化率，可見發(fā)達(dá)國家大船事故發(fā)生概率比小船高14.52%，而發(fā)展中國家大船事故發(fā)生概率比小船高8.38%。最后一行的變化率表明，無論是大船還是小船，發(fā)達(dá)國家船舶事故發(fā)生概率均比發(fā)展中國家高60%左右。

表3為船旗對航運(yùn)安全影響的敏感性分析。由后驗(yàn)概率變化率可知，開放式旗船事故發(fā)生概率平均比封閉式旗船高70%左右。

表2 船舶尺度對航運(yùn)安全影響的敏感性分析 %

表3 船旗對航運(yùn)安全影響的敏感性分析 %

表4為船齡對航運(yùn)安全影響的敏感性分析。由表4可知：舊船引發(fā)事故的概率平均比新船高40%～50%。

表4 船齡對航運(yùn)安全影響的敏感性分析 %

表5為船舶類型對航運(yùn)安全影響的敏感性分析。由表5可知：發(fā)達(dá)國家的船舶中客船事故發(fā)生概率最高，為5.7%；其次是集裝箱船，為5.11%；其后是一般貨船、油船和干散貨船。發(fā)展中國家各類型船舶發(fā)生事故的概率均低于發(fā)達(dá)國家，但相對排序類似。

表5 船舶類型對航運(yùn)安全影響的敏感性分析 %

表6為檢查次數(shù)對航運(yùn)安全影響的敏感性分析。由表6可知，發(fā)生事故概率較高的船舶累計PSC檢查的次數(shù)也較多。該結(jié)果與預(yù)期相反，當(dāng)前各大MoU均建立有一套較為系統(tǒng)、科學(xué)的船舶監(jiān)控機(jī)制，可滯留抽查出的問題船舶，并責(zé)令其修補(bǔ)缺陷。一般認(rèn)為被港口國監(jiān)控檢查次數(shù)越多，越有助于其提高船舶質(zhì)量，從而降低事故發(fā)生概率。[13]

表6 檢查次數(shù)對航運(yùn)安全影響的敏感性分析 %

該結(jié)果表明，當(dāng)前港口PSC監(jiān)控對篩選易發(fā)事故船舶非常有效，能準(zhǔn)確篩選出問題船舶并進(jìn)行抽查，但因被檢查船舶存在的嚴(yán)重問題未得到全面改善，并未消除其后續(xù)事故發(fā)生的可能性。

表7為檢查次數(shù)對滯留情況影響的敏感性分析。由表7可知：若船舶在港口國被檢查超過10次，則滯留的概率比檢查次數(shù)少的船舶高出很多；而發(fā)展中國家船舶滯留的概率高于發(fā)達(dá)國家。這充分說明被多次檢查并滯留后并未有效敦促船舶所有人提高船舶質(zhì)量。

2.2 模型整體分析

表8為模型整體有效性分析。由表8可知，當(dāng)船舶為大船且懸掛的是開放式船旗時，發(fā)達(dá)國家和發(fā)展中國家發(fā)生嚴(yán)重事故的概率增加到4.91%及3.05%。依此類推，發(fā)達(dá)國家和發(fā)展中國家在各因素最差狀態(tài)下發(fā)生事故的概率分別高達(dá)21.36%及9.68%。這說明船舶尺度較大、懸掛開放式船旗、船齡越大、被港口檢查次數(shù)較多的客船或集裝箱船，發(fā)生事故的概率較高。

與各因素對較輕事故的影響相比，各因素對嚴(yán)重事故發(fā)生概率的影響均較高，這在一定程度上說明盡管船舶發(fā)生事故的概率相對較低，但其發(fā)生的事故一般都較為嚴(yán)重，這也是各傳統(tǒng)航運(yùn)國、沿岸港口國及國際海事組織對航運(yùn)安全較為重視的原因。因此，該研究結(jié)果可幫助各相關(guān)國際、國內(nèi)機(jī)構(gòu)完善航運(yùn)政策，有效降低航運(yùn)事故發(fā)生的概率。

表7 檢查次數(shù)對滯留情況影響的敏感性分析 %

表8 模型整體有效性分析 %

3 結(jié)束語

本文將貝葉斯網(wǎng)絡(luò)與Logistic回歸相結(jié)合，研究各類因素對船舶事故的影響，挖掘當(dāng)前海上交通安全管理中存在的問題。分析發(fā)現(xiàn)，不論是發(fā)達(dá)國家還是發(fā)展中國家，易發(fā)生事故的船舶均已被港口PSC檢查多次，滯留概率較高的船舶同樣被檢查過多次。這充分說明當(dāng)前各MoU制定的船舶監(jiān)控制度具有一定的有效性。然而，因未能對高風(fēng)險船舶實(shí)施進(jìn)一步的控制，致使事故發(fā)生概率未能從根本上降低。船舶類型、船齡和船舶尺度對事故發(fā)生概率也有顯著影響。結(jié)合研究結(jié)果，為加強(qiáng)海上航行安全，提出以下建議：

1) 加強(qiáng)各MoU之間在船舶檢查數(shù)據(jù)上的溝通，建立一個完善的船舶監(jiān)控數(shù)據(jù)庫，以更好地對經(jīng)停的累積檢查次數(shù)過多的船舶進(jìn)行限制。建議對船舶累積檢查次數(shù)進(jìn)行分類，輕者限制船舶航行的某些權(quán)限(比如載重等)，重者禁止在某些海域航行。

2) 對不同類型的船舶進(jìn)行不同程度的檢查。不同類型船舶運(yùn)營條件的嚴(yán)格程度不同，船舶發(fā)生事故的概率也有所差別。

3) 敦促完善港口國PSC監(jiān)控的處罰力度，有效促進(jìn)船舶質(zhì)量的提升。由于船舶涉及到的相關(guān)國家較多，各相關(guān)國家的海事組織對船舶的監(jiān)控不力，停經(jīng)港口要求修補(bǔ)缺陷或滯留的處罰并未對其造成實(shí)質(zhì)性影響，使得缺陷船舶沒有動力通過增加投入來改善船舶質(zhì)量?？梢姡ㄟ^IMO協(xié)調(diào)各國家建立一套完善的監(jiān)控、處罰機(jī)制，可更加有效地促進(jìn)全球航運(yùn)安全的提升。

[1] UNCTAD. Review of Maritime Transport[R]. Geneva: UNCTAD, 2015.

[2] 呂靖, 高天航. 海上通道關(guān)鍵節(jié)點(diǎn)安全保障效率研究[J]. 中國軟科學(xué), 2015(10):1-8.

[3] LUO M. Emission Reduction in International Shipping—The Hidden Side Effects[J]. Maritime Policy & Management, 2013,40(7):694-708.

[4] MANSELL JNK. Development of International Law for Flag State Responsibility[Z] Berlin, Heidelberg: Springer-Verlag Berlin Heidelberg, 2009: 53-69.

[8] BANDA OAV, HNNINEN M, GOERLANDT F, et al. Bayesian Networks as a Decision Making Tool to Plan and Assess Maritime Safety Management Indicators[C]∥Probabilistic Safety Assessment and Management PSAM 12. Honolulu, Hawaii,2014.

[9] 周國華, 彭波. 基于貝葉斯網(wǎng)絡(luò)的建設(shè)項(xiàng)目質(zhì)量管理風(fēng)險因素分析——以京滬高速鐵路建設(shè)項(xiàng)目為例[J]. 中國軟科學(xué), 2009(9):99-106.

[10] YIN JB, LI KX, FAN LX. A New Safety Index for Maritime Transportation [J]. Transportation Research Part A: Policy and Practice, 2014,66:75-87.

[11] LI K, YIN J, BANG H, et al. Bayesian Network with Quantitative Input for Maritime Risk Analysis[J]. Transportmetrica A: Transport Science, 2014,2:89-118.

[12] LADAN M, HNNINEN M. Data Sources for Quantitative Marine Traffic Accident modeling[M]. Helsinki: School of Engineering, 2012.

[13] CARIOU P, WOLFF F. Identifying Substandard Vessels Through Port State Control Inspections: A New Methodology for Concentrated Inspection Campaigns[J]. Marine Policy, 2015,60:27-39.