廖詩(shī)管，楊冬，白茜文，翁金賢,1b

(1.上海海事大學(xué)，a.海洋科學(xué)與工程學(xué)院，b.交通運(yùn)輸學(xué)院，上海201306；2.香港理工大學(xué)，物流及航運(yùn)學(xué)系，香港999077；3.清華大學(xué)，工業(yè)工程系，北京100084)

0 引言

港口裝卸效率是港口競(jìng)爭(zhēng)力的重要指標(biāo)，對(duì)于縮短船舶在港停泊時(shí)間、加快船舶周轉(zhuǎn)、降低運(yùn)輸成本起著至關(guān)重要的作用，直接影響著船公司對(duì)航線掛靠港的選擇[1]。評(píng)估港口的裝卸效率需要考慮許多因素的影響，包含可獲取因素(岸吊數(shù)量等)、難以獲取因素(港口的操作和管理水平等)和隨機(jī)因素(天氣和岸吊機(jī)械故障等)。難以獲取因素和隨機(jī)因素為評(píng)估港口實(shí)際裝卸效率帶來(lái)巨大挑戰(zhàn)，使港口間的效率衡量及比較沒有唯一的方法可遵循。為準(zhǔn)確估計(jì)港口的裝卸效率值，比較不同港口裝卸效率的差異，亟需提出一種標(biāo)準(zhǔn)化的港口裝卸效率值估計(jì)方法，以更客觀地評(píng)估和比較復(fù)雜多因素下的港口裝卸效率，為港口提升管理水平和船公司決策提供參考。

較為典型的港口效率評(píng)估方法是Charnes 等[2]率先提出的數(shù)據(jù)包絡(luò)分析法(DEA)，該方法得到了廣泛的運(yùn)用與發(fā)展。例如：李丹等[3]采用改進(jìn)的三階段DEA 模型對(duì)我國(guó)沿海集裝箱碼頭2008-2012年的運(yùn)營(yíng)效率進(jìn)行測(cè)度分析；Cheon 等[4]構(gòu)建了基于DEA 模型的Malmquist 生產(chǎn)力指數(shù)(MPI)模型，用以評(píng)估世界集裝箱港口效率變化情況；Wu等[5]應(yīng)用DEA 技術(shù)估算了集裝箱港口效率，他們認(rèn)為難以觀測(cè)的非有形因素會(huì)影響集裝箱港口的效率。DEA更多的是考慮宏觀指標(biāo)(資本、勞動(dòng)和土地等)對(duì)港口效率變化的影響，對(duì)微觀運(yùn)作指標(biāo)(港口的裝卸設(shè)備和運(yùn)營(yíng)條件等)缺乏考慮。陳龍彪等[6]估算了反映港口生產(chǎn)力的指標(biāo)，利用了岸吊設(shè)備對(duì)碼頭作業(yè)效率進(jìn)行評(píng)估。Zhong 等[7]使用海事大數(shù)據(jù)，提出一種識(shí)別船舶裝卸作業(yè)事件的算法，估算了石油港口的裝卸效率、吞吐量等5個(gè)港口性能指標(biāo)。虞哲人[8]分析了新型岸吊設(shè)備對(duì)提高碼頭作業(yè)效率的作用，并預(yù)測(cè)其應(yīng)用前景。此類方法雖然考慮了港口的微觀運(yùn)作數(shù)據(jù)，但依賴較為詳細(xì)的碼頭裝卸設(shè)備的數(shù)據(jù)資料，故在應(yīng)用擴(kuò)展性上有局限。

AIS提供全球船舶實(shí)時(shí)動(dòng)態(tài)信息，被廣泛認(rèn)為是未來(lái)改變航運(yùn)業(yè)運(yùn)營(yíng)模式的大數(shù)據(jù)來(lái)源。近年來(lái)，基于AIS利用大數(shù)據(jù)技術(shù)發(fā)展出來(lái)的實(shí)際應(yīng)用越來(lái)越廣，評(píng)估港口競(jìng)爭(zhēng)力指標(biāo)被認(rèn)為是其未來(lái)的發(fā)展方向之一[9]。本文利用GMap可視化技術(shù)和港口的歷史裝卸數(shù)據(jù)，提出一種基于船舶大數(shù)據(jù)的港口裝卸效率值的估計(jì)方法，以評(píng)估多因素和微觀運(yùn)作數(shù)據(jù)不足的復(fù)雜條件下的港口績(jī)效指標(biāo)。該指標(biāo)所具備的實(shí)用價(jià)值，可為港口管理提供標(biāo)準(zhǔn)，同時(shí)用于估算和實(shí)時(shí)監(jiān)控港口的吞吐量，最終提升港口數(shù)字化管理水平。本文以2017年全球集裝箱吞吐量前四的港口作為算例，證明該方法的有效性。

1 港口裝卸效率值計(jì)算模型

1.1 港口裝卸效率值計(jì)算公式

《海港總體設(shè)計(jì)規(guī)范》[10](簡(jiǎn)稱《規(guī)范》)是我國(guó)根據(jù)近十多年來(lái)沿海港口總體設(shè)計(jì)經(jīng)驗(yàn)及專家知識(shí)制定而成。《規(guī)范》的主要內(nèi)容包括港址選擇、設(shè)計(jì)基礎(chǔ)條件、裝卸工藝等。其中，裝卸工藝所在章節(jié)中提出了港口裝卸效率值的估計(jì)方法。

根據(jù)《規(guī)范》，集裝箱船設(shè)計(jì)船時(shí)效率pi表示為船舶i?？扛劭跁r(shí)，單位時(shí)間內(nèi)的裝卸量，即

式中：Eport為港口的裝卸效率值；ni為滿足船舶i裝卸所需配備的岸吊臺(tái)數(shù)；qi為船舶i所配備的岸吊臺(tái)時(shí)效率基準(zhǔn)值(TEU·h-1)；ki為船舶i所配備的岸吊同時(shí)作業(yè)率；ri為船舶i的作業(yè)倒箱率。

《規(guī)范》中給定的港口裝卸效率值Eport定義為兩個(gè)因素的乘積，即集裝箱標(biāo)準(zhǔn)箱折算系數(shù)和新型裝卸設(shè)備(可吊雙箱的集裝箱岸吊等)對(duì)裝卸效率的提升系數(shù)，單位為無(wú)量綱。在實(shí)際生產(chǎn)過程中，不確定因素也可能影響港口的裝卸效率(不良天氣，設(shè)備故障等)。因此，《規(guī)范》中也僅給出裝卸效率值的經(jīng)驗(yàn)取值范圍，并沒有具體的計(jì)算方法。

理論上，港口實(shí)際吞吐量是所有靠港船舶的設(shè)計(jì)船時(shí)效率與靠泊時(shí)長(zhǎng)的乘積，即

式中：Qa為港口實(shí)際吞吐量(TEU)；ti為船舶i的靠泊時(shí)長(zhǎng)(h)；I為所有靠港船舶的集合。

船舶自動(dòng)識(shí)別系統(tǒng)(AIS)實(shí)時(shí)采集船舶在港口中的位置信息[6]。港口操作的隨機(jī)因素和難以觀測(cè)因素可通過歷史吞吐量與靠泊時(shí)間長(zhǎng)短觀測(cè)而得。AIS提供了一個(gè)提取靠泊時(shí)間的平臺(tái)，可通過設(shè)計(jì)算法識(shí)別船舶靠離泊時(shí)刻，估計(jì)靠泊時(shí)長(zhǎng)(詳見第1.2 節(jié))。根據(jù)《規(guī)范》提供的ni,qi,ki,ri的取值范圍，本文提出其唯一的取值方法(詳見第2.3 節(jié))，從而估計(jì)港口的裝卸效率值，即

本文提出的港口裝卸效率值估計(jì)方法，彌補(bǔ)了現(xiàn)有文獻(xiàn)對(duì)不可知或不確定因素考慮不足的情形，把《規(guī)范》中效率的范圍取值進(jìn)一步估算為確切的港口裝卸效率，從而使實(shí)時(shí)大規(guī)模港口吞吐量的估計(jì)和監(jiān)控成為可能。

1.2 靠泊時(shí)長(zhǎng)估計(jì)方法

船舶靠泊時(shí)長(zhǎng)一般可理解為船舶離泊時(shí)刻與船舶靠泊時(shí)刻的時(shí)間差，準(zhǔn)確地識(shí)別船舶靠離泊時(shí)刻是估計(jì)港口裝卸效率值的關(guān)鍵。本文通過兩個(gè)步驟估計(jì)船舶的靠泊時(shí)長(zhǎng)：步驟1，根據(jù)衛(wèi)星地圖和專業(yè)網(wǎng)站實(shí)時(shí)地圖信息，識(shí)別港口泊位位置以確定船舶靠泊區(qū)域；步驟2，輸入AIS數(shù)據(jù)中船舶位置信息，利用GMap可視化技術(shù)確定船舶是否進(jìn)入泊位區(qū)域以確定靠離泊時(shí)刻。圖1為船舶靠泊時(shí)長(zhǎng)估計(jì)方法的框架。

圖1 船舶靠泊時(shí)長(zhǎng)估計(jì)方法框架Fig.1 Methodology framework for ship berthing time estimation

如圖2所示，由于集裝箱碼頭分布著大量岸吊及堆場(chǎng)，從衛(wèi)星圖像中能夠明顯地觀察出這些特征，初步識(shí)別泊位區(qū)域。船訊網(wǎng)、勞氏船級(jí)社等專業(yè)網(wǎng)站提供了實(shí)時(shí)船舶位置信息，通過觀察是否有船舶?？吭诔醪阶R(shí)別的泊位區(qū)域，可校正所識(shí)別泊位區(qū)域的準(zhǔn)確性。本文定義的泊位區(qū)域?qū)挾纫圆次痪€沿水域方向100 m劃分[11]。

圖2 泊位區(qū)域確定方法Fig.2 Method for berth area identification

圖3為將原始AIS數(shù)據(jù)去噪后提取的實(shí)時(shí)船舶運(yùn)動(dòng)軌跡。通過GMap可視化技術(shù)，將標(biāo)定的泊位區(qū)域與實(shí)時(shí)船舶運(yùn)動(dòng)軌跡結(jié)合，確定船舶靠離泊時(shí)刻。例如圖3中的船舶MMSI 編號(hào)為477335600，在2017年8月26日12:56:56 識(shí)別為靠泊時(shí)刻，在2017年8月26日20:05:43識(shí)別為離泊時(shí)刻，靠泊時(shí)長(zhǎng)約為7 h。

圖3 基于GMap可視化技術(shù)的船舶靠離泊時(shí)刻識(shí)別Fig.3 Method for GMap-visualization-technology-based ship berthing and unberthing time identification

2 數(shù)據(jù)收集與分析

2.1 集裝箱吞吐量統(tǒng)計(jì)數(shù)據(jù)

規(guī)模以上(例如年吞吐量達(dá)到200萬(wàn)TEU)港口統(tǒng)計(jì)集裝箱吞吐量一般以月度數(shù)據(jù)為主，本文基于中國(guó)交通運(yùn)輸部和新加坡海事及港務(wù)管理局官方統(tǒng)計(jì)數(shù)據(jù)，梳理出2017年全球集裝箱吞吐量排名前四的港口：上海港、新加坡港、深圳港和寧波-舟山港的月度集裝箱吞吐量數(shù)據(jù)，即Qa的值，結(jié)果如表1所示。

表1 2017年4個(gè)港口集裝箱吞吐量月度數(shù)據(jù)Table 1 Monthly container throughput data of four ports in 2017

2.2 AIS數(shù)據(jù)

基于2017年全球AIS 數(shù)據(jù)，時(shí)間頻率為10～15 s，共計(jì)超2 TB的數(shù)據(jù)容量。收集的數(shù)據(jù)具有較高的數(shù)據(jù)精度，為識(shí)別船舶的靠離泊時(shí)刻提供可靠的數(shù)據(jù)保障，也能夠保證計(jì)算和驗(yàn)證所提出方法的有效性。

共收集到全球6320 艘有MMSI(船舶編碼)的集裝箱船數(shù)據(jù)，經(jīng)匹配AIS數(shù)據(jù)獲得4787艘船舶運(yùn)動(dòng)軌跡數(shù)據(jù)，每月數(shù)據(jù)規(guī)模達(dá)到約5 GB。匹配結(jié)果少于總集裝箱船數(shù)量，可能原因?yàn)椋?/p>

(1)某些船舶已經(jīng)報(bào)廢；

(2)部分船舶經(jīng)改造成非集裝箱船，造成MMSI改變；

(3)部分集裝箱船在一段時(shí)間內(nèi)并不進(jìn)行運(yùn)輸作業(yè)；

(4)2017年之后新造的船舶。

2.3 計(jì)算參數(shù)取值

依據(jù)《規(guī)范》，ni,qi,ki,ri的取值范圍僅與船舶噸級(jí)(GDWT)或載箱量(CTEU)有關(guān)，如表2和表3所示。其中，ni和qi的取值為整數(shù)，ki的取值隨載箱量的增大而減小，ri的取值隨船舶載箱量的增大而增大。

表2 不同噸級(jí)的集裝箱船岸吊配備數(shù)量Table 2 Number of quay cranes for container ships by different tonnages

表3 不同載箱量的集裝箱船岸吊臺(tái)時(shí)效率、同時(shí)作業(yè)率及倒箱率Table 3 Loading and unloading efficiency of quay cranes,simultaneous operation rate and turnover rate of container ships with different TEU capacity

進(jìn)一步確定4個(gè)參數(shù)唯一取值的方法，即

例如：依據(jù)表2，當(dāng)船舶噸級(jí)在[27500,45000)時(shí)，ni的取值按船舶噸級(jí)范圍的中值分為兩組，即ni=2,GDWT∈(2 7500,36250);ni=3,GDWT∈[3 6250,45000)。依據(jù)表3，當(dāng)船舶載箱量在(1900,5650]時(shí)，qi的取值根據(jù)船舶載箱量范圍按組距等量分為6 組，即其中表示向下取整符號(hào)；對(duì)應(yīng)ki的取值采用兩點(diǎn)之間的一次方程確定唯一取值，即ki=-0.00003CTEU+0.9507,CTEU∈(1 900,5650 ]；ri的取值按給定范圍的均值做近似估計(jì)，即ri=0.035,CTEU∈(1 900,9500 ]。

3 計(jì)算結(jié)果分析

3.1 集裝箱船靠泊時(shí)長(zhǎng)計(jì)算

根據(jù)1.2節(jié)提出的方法和所匹配的4787艘AIS數(shù)據(jù)，估計(jì)2017年上海港、新加坡港、深圳港和寧波-舟山港集裝箱船靠泊時(shí)長(zhǎng)的月度數(shù)據(jù)，利用@RISK軟件對(duì)1月的單船靠泊時(shí)長(zhǎng)進(jìn)行統(tǒng)計(jì)分析，結(jié)果如圖4所示。上海港單船靠泊時(shí)長(zhǎng)呈Person分布，平均靠泊時(shí)長(zhǎng)為14.56 h；新加坡港單船靠泊時(shí)長(zhǎng)呈LogLogistic 分布，平均靠泊時(shí)長(zhǎng)為15.58 h；深圳港單船靠泊時(shí)長(zhǎng)呈LogLogistic分布，平均靠泊時(shí)長(zhǎng)為13.89 h；寧波-舟山港單船靠泊時(shí)長(zhǎng)呈InvGauss分布，平均靠泊時(shí)長(zhǎng)為17.25 h。本文關(guān)于船舶靠泊時(shí)間分布規(guī)律與現(xiàn)有文獻(xiàn)[11]研究結(jié)果相似。

圖4 1月的靠泊時(shí)長(zhǎng)分布結(jié)果Fig.4 Distribution of calculated berthing time in January

表4為4個(gè)港口1月靠泊船舶的大小分布。結(jié)果表明，上海港和新加坡港靠泊船舶大小分布類似；深圳港和寧波-舟山港相對(duì)來(lái)說大型船舶靠泊比例較高。一般地，大型船舶的單位集裝箱操作時(shí)間較長(zhǎng)。

表4 1月份靠泊船舶的大小分布Table 4 Size distribution of ships berthing in January

3.2 港口裝卸效率值計(jì)算

根據(jù)式(3)和第2.3 節(jié)給出的參數(shù)取值方法，在求得船舶靠泊時(shí)長(zhǎng)的情況下，對(duì)4個(gè)港口2017年上半年的裝卸效率值進(jìn)行計(jì)算，結(jié)果如圖5所示。

圖5 2017年上半年港口裝卸效率值估計(jì)結(jié)果Fig.5 Estimated handling efficiency of ports in first half of 2017

研究結(jié)果表明：

(1)2017年上半年，上海港、新加坡港、深圳港和寧波-舟山港的裝卸效率月度均值分別為2.86、1.87、2.17 和2.10。在靠泊船舶的大小分布相似的情況下，上海港的裝卸效率高于新加坡港，深圳港的裝卸效率高于寧波-舟山港。根據(jù)2015年ITF/OECD報(bào)告結(jié)果，新加坡港每船每小時(shí)平均集裝箱移動(dòng)量?jī)H為66 TEU，而上海港為167TEU，深圳港為133 TEU[12]。本文研究結(jié)論與該結(jié)果一致。

(2)各港口月度裝卸效率值起伏變化不明顯，說明港口裝卸效率在一定時(shí)間內(nèi)不會(huì)有較大變化。一般地，港口引進(jìn)新型設(shè)備或智能化技術(shù)能夠明顯提高港口裝卸效率[8]。因此，有理由認(rèn)為港口裝卸效率在一段時(shí)間內(nèi)相對(duì)穩(wěn)定，較少出現(xiàn)突變特征。

3.3 算例驗(yàn)證

為驗(yàn)證所提出方法的準(zhǔn)確性，利用所估計(jì)的上半年裝卸效率值對(duì)4 個(gè)港口下半年的吞吐量進(jìn)行估計(jì)。根據(jù)式(2)，港口吞吐量的估計(jì)值可表示為

式中：為估計(jì)的港口裝卸效率值。假設(shè)港口下半年裝卸效率值維持上半年月度均值水平，估計(jì)下半年港口吞吐量與實(shí)際吞吐量的誤差，結(jié)果如圖6所示。其中，誤差百分比為正數(shù)表示估計(jì)吞吐量大于實(shí)際吞吐量，反之亦然。

圖6 2017年下半年港口的估計(jì)吞吐量與實(shí)際吞吐量誤差結(jié)果Fig.6 Estimation error of ports in second half of year of 2017

研究結(jié)果表明：

(1)上海港、新加坡港、深圳港和寧波-舟山港下半年月度吞吐量估計(jì)誤差絕對(duì)值的均值分別為2.77%、2.06%、2.93%和2.46%。

(2)4 個(gè)港口月度吞吐量的估計(jì)均值準(zhǔn)確率在97%以上，故認(rèn)為本文所提出的估計(jì)港口裝卸效率值的方法具有較高的準(zhǔn)確性，能夠在一定程度上反映港口確切的裝卸效率。

(3)基于AIS 實(shí)時(shí)數(shù)據(jù)，利用該方法可對(duì)全球集裝箱港口吞吐量進(jìn)行實(shí)時(shí)測(cè)算和監(jiān)控。

4 結(jié)論

本文基于高頻AIS船舶大數(shù)據(jù)，提出一種基于圖像技術(shù)的港口裝卸效率值的計(jì)算方法，用于評(píng)估多因素和微觀運(yùn)作數(shù)據(jù)不足的復(fù)雜條件下的港口實(shí)際裝卸效率。通過對(duì)2017年全球集裝箱吞吐量排名前四的港口裝卸效率值進(jìn)行估算發(fā)現(xiàn)，上海港的裝卸效率高于新加坡港，深圳港的裝卸效率高于寧波-舟山港，利用估算的效率值估算4 個(gè)港口的吞吐量，準(zhǔn)確率可達(dá)到97%以上，這表明此方法能夠準(zhǔn)確地對(duì)不同港口間裝卸效率進(jìn)行估算和比較。本文提出的效率參數(shù)計(jì)算方法，可用于估算全球集裝箱港口或碼頭的實(shí)時(shí)吞吐量，極大的減少了人工統(tǒng)計(jì)產(chǎn)生的成本和誤差，使市場(chǎng)數(shù)據(jù)變得透明，不僅為港口提升管理水平和船公司選擇掛靠港口的決策提供理論參考，也有助于提升航運(yùn)港口業(yè)數(shù)字化的科研和發(fā)展。