汪士寒，張麗娜

(1.寧波財經(jīng)學院國際經(jīng)濟貿(mào)易學院，浙江寧波 315175；2.浙江海洋大學港航與交通運輸工程學院，浙江舟山 316022)

科學合理的碼頭泊位數(shù)量布局是提升港口競爭力的一個重要方面，泊位數(shù)過剩會造成岸線資源浪費和費用成本上升，泊位數(shù)不足會導(dǎo)致到港船舶靠泊等待時間過長使港口失去市場吸引力。鼠浪湖碼頭占有全國可靠泊40萬t級礦石船的7個泊位中的2個，是世界第一大港寧波舟山港的鐵礦石集疏運系統(tǒng)中最重要的碼頭，也是國內(nèi)最大的鐵礦石中轉(zhuǎn)碼頭和國家最重要的鐵礦石儲備中轉(zhuǎn)基地之一。但隨著鐵礦石吞吐量的日益增長，碼頭泊位數(shù)量配置問題凸顯。筆者根據(jù)碼頭吞吐量的近期實際和遠期預(yù)測，利用排隊論模型，在岸線資源容量、航道通航能力上限之下，以港、船雙方發(fā)生和損失的綜合費用最小為目標，對鼠浪湖碼頭泊位數(shù)量布局優(yōu)化問題進行研究，以期為鼠浪湖鐵礦石碼頭的規(guī)劃發(fā)展提供科學依據(jù)。

1 研究思路與分布檢驗

1.1 研究原理與思路

船舶到港分布是一種典型的隨機分布，泊位是服務(wù)機構(gòu)中的服務(wù)臺，船舶和泊位構(gòu)成了一個排隊系統(tǒng)。船舶到港后排隊靠泊接受服務(wù)屬于排隊論中描述的排隊現(xiàn)象的一種，確定碼頭的最優(yōu)泊位數(shù)適用排隊模型。

港口排隊服務(wù)系統(tǒng)中，不同的輸入過程、排隊規(guī)則和服務(wù)機構(gòu)適用不同的模型，符合港口服務(wù)系統(tǒng)的模型主要有M/M/S和M/Ek/S。在等待制的排隊規(guī)則下，當船舶到港和占泊服務(wù)時間服從泊松分布和負指數(shù)分布時，適用M/M/S模型；當船舶到港和服務(wù)時間服從泊松分布和K階愛爾朗分布時，適用M/Ek/S模型[1-2]。

等待制規(guī)則：不考慮因其他因素使船舶中途退出的情況，船舶排隊隊長無限制，泊位對于待泊船舶采用先到先服務(wù)原則，泊位一有空閑，船舶立即前去?？?，并且在港的每艘船舶早晚都會服務(wù)到。

港口系統(tǒng)的排隊規(guī)則屬于等待制，因此本文在判別鼠浪湖碼頭船舶到港與占泊時間的分布屬性的基礎(chǔ)上選用適用的模型，根據(jù)碼頭吞吐量的近期實際和遠期預(yù)測，在綜合費用最小的目標和資源限制的約束下，得到鼠浪湖碼頭泊位近期與遠期的最佳數(shù)量配置。

1.2 碼頭設(shè)施與資源

1.2.1 碼頭泊位與基礎(chǔ)設(shè)施

寧波舟山港鼠浪湖碼頭2016年建成投入營運，目前在營泊位5個，其中2個30～40萬t級，1個10萬t級和2個5萬t級(詳見表1)。碼頭現(xiàn)有5臺卸船機和2臺裝船機，卸船機額定能力3 000 t·h-1，裝船機額定能力為 7 500 t·h-1。

表1 鼠浪湖碼頭泊位情況Tab.1 Berth situation at Shulanghu wharf

1.2.2 航道現(xiàn)狀與通航能力上限

鼠浪湖碼頭位于蛇移門水道東部，其所用航道有南、北向主航道及其支航道，航道水深25.1 m，可滿足進出本碼頭船舶安全航行要求(詳見表2)。

表2 鼠浪湖碼頭蛇移門航道情況一覽表Tab.2 A list of the status quo of Sheyimen waterway at Shulanghu wharf

1.2.3 岸線資源上限

根據(jù)舟山市港航管理局《舟山港域碼頭設(shè)施與資源報告》，鼠浪湖島擁有岸線16.03 km，其中水深25 m以上岸線1.8 km，可建(30～40)萬t級泊位4～5個；水深15 m以上岸線1.9 km，可建10萬t級左右泊位5～6個；另有5萬t級以下泊位8個左右的建設(shè)空間。

1.3 分布檢驗

多年來國內(nèi)外大量資料和學者研究證實船舶到港和服務(wù)時間服從泊松分布和負指數(shù)分布[3-4]。為嚴密起見，這里根據(jù)船舶的到港信息用χ2檢驗來驗證該港口船舶到港和占泊服務(wù)時間服從泊松分布和負指數(shù)分布。

式中，g 為樣本總體 X 實數(shù)軸上所取得等距區(qū)間數(shù)(-∞，t1)，(t1，t2)，……(tg-1，+∞)；n 為到達港口的船舶數(shù)；fi表示頻數(shù)，即樣本中 x1，x2，x3，…，xn 發(fā)生在第 i段區(qū)間的個數(shù)(i=1，2，3，…，g)，其對應(yīng)的頻率是 fi/n；pi是頻率fi/n將發(fā)生在第i段區(qū)間上相應(yīng)的概率，其對應(yīng)的理論頻數(shù)為Fi=npi[5]。

用n(t)表示在時間間隔t內(nèi)到達船舶的數(shù)量，則當船舶到達符合泊松分布時，在時間間隔t內(nèi)到達n(t)艘船舶的概率P可以表示為：

式中：λ表示平均到達率，即單位時間間隔內(nèi)平均到達的船舶數(shù)。μ表示單位時間間隔內(nèi)單個泊位完成裝卸服務(wù)的船舶數(shù)，即泊位的平均服務(wù)率。

鼠浪湖碼頭是以鐵礦石運輸為主的裝卸中轉(zhuǎn)碼頭，實際可利用裝卸作業(yè)的泊位有5個。2018年鼠浪湖碼頭到港的鐵礦石船舶共1 224艘，船舶平均裝卸量40 866 t·艘-1。

據(jù)碼頭所在地港航氣象統(tǒng)計，鼠浪湖碼頭建成以來雖然遇到的臺風等天氣較多，但無大破壞力的極其惡劣天氣，港口營運天數(shù)為365 d，這樣船舶日平均到達率λ=1 224/365=3.353 4艘·d-1。

根據(jù)鼠浪湖碼頭統(tǒng)計數(shù)據(jù)計算得出單個泊位的平均裝卸一艘次船的時間為1.17 d，則μ=1/1.17=0.85艘次·d-1。統(tǒng)計1 d內(nèi)實際到達船舶的經(jīng)驗頻數(shù)fi(i=0，1，2，3…10)，由此計算出1 d內(nèi)船舶到達的經(jīng)驗頻率。根據(jù)公式(2)求得符合泊松分布一天內(nèi)到港船舶 1，2，3，…，9，10 艘的理論頻率 pi(i=1，2，3….10)，再求出相應(yīng)的理論頻數(shù)Fi。計算結(jié)果如表3所示。

表3 各個參數(shù)計算結(jié)果Tab.3 Calculation results of each parameter

根據(jù)表3數(shù)據(jù)運用χ2檢驗來判斷鼠浪湖碼頭船舶到港是否服從泊松分布。

(1)建立原假設(shè)H0：船舶到港服從泊松分布。

(2)將表中數(shù)據(jù)代入公式(1)(g=9組，將理論頻數(shù)≤5的合并為一組)，求得χ2=2.70

自由度 v=9-2=7，取 α=0.05，查 χ2分布表得 χ20.05=14.067

(3)得出結(jié)果χ20.05≥χ2，可知原假設(shè)成立，鼠浪湖碼頭船舶到港服從泊松分布。且實際到港的船舶相繼達到的時間間隔相互獨立，占用泊位服務(wù)時間服從負指數(shù)分布。

2 基本模型

分布檢驗證明輸入過程和服務(wù)時間服從泊松分布與負指數(shù)分布，適用M/M/S模型。

2.1 排隊模型

設(shè)泊位數(shù)為S，在平均到達率λ不變時，當?shù)竭_船數(shù)n＜S時，系統(tǒng)的平均服務(wù)效率μ會隨著到達船數(shù)的增加而發(fā)生改變μ、2μ、3μ……(S-1)μ，當?shù)竭_船數(shù)n≥S時，最多有S艘船舶接受服務(wù)，n-S艘船舶待泊，平均服務(wù)效率保持臨界值Sμ不變[6]。

圖1 M/M/S排隊系統(tǒng)的狀態(tài)轉(zhuǎn)移圖Fig.1 State transition diagram of M/M/S queuing system

如圖1所示，從狀態(tài)(1)轉(zhuǎn)移為狀態(tài)(0)表示當系統(tǒng)中只有一艘船舶接受服務(wù)并離開的轉(zhuǎn)移率為μP1，從狀態(tài)(2)轉(zhuǎn)移到狀態(tài)(1)表示系統(tǒng)中有兩艘船在泊位接受服務(wù)，其中一艘船接受服務(wù)完畢離開的轉(zhuǎn)移率為2μP2，以此類推得到穩(wěn)態(tài)狀態(tài)下系統(tǒng)的方程組：

用遞推法求解方程得穩(wěn)態(tài)狀態(tài)下無船在港的概率P0[7]

《菜根譚》中有言：冷眼觀人，冷耳聽語，冷情當感，冷心思理。處世“冷”一點，遇事“靜”一點，才不會趕走好運氣。

有n艘船在港的概率Pn·S

2.2 費用模型

本研究優(yōu)化目標涉及相關(guān)費用，推算費用模型如下：

式中，Tb為泊位年總費用；Ts為船舶在港年總費用；T為船、港年總費用；cb為泊位日均營運費用；cs為船舶日均在港待泊費用；Cb為泊位年閑置費用；Cs為船舶年待泊損失費用；ρs為泊位數(shù)為S的泊位利用率；N為365 d。

3 優(yōu)化計算

3.1 近期優(yōu)化

鼠浪湖碼頭當前有關(guān)數(shù)據(jù)：單個泊位的平均裝卸時間是1.17 d·艘-1；單個服務(wù)臺的平均服務(wù)率μ=0.85；日均到達率λ=3.35艘·d-1；cb可近似估算為投建泊位的成本和港內(nèi)動力照明費用、碼頭設(shè)施修理費用、設(shè)施折舊費用等部分的花費均攤到單個泊位的費用總和[8]，經(jīng)綜合加權(quán)取平均值得cb=21.20萬元；根據(jù)實際到港船型的平均情況調(diào)研與計算得cs=26.48萬元。應(yīng)用前述的模型，計算得到相關(guān)數(shù)據(jù)如表4和表5所示。

表4 M/M/S模型最優(yōu)泊位數(shù)計算結(jié)果Tab.4 Optimal berth number calculation results by M/M/S model

表5 不同泊位數(shù)產(chǎn)生的泊位閑置損失費和船舶待泊損失費用計算結(jié)果Tab.5 Calculation results of berth idle loss and berth waiting loss incurred by different

由表4和表5可以看出泊位數(shù)越多，Ts和Cs越低，但Tb和Cb越高，因此，碼頭在規(guī)劃泊位時不能盲目建設(shè)，要科學規(guī)劃泊位數(shù)。

由表4和表5可知，當泊位數(shù)S=6時，發(fā)生在港、船雙方的綜合費用T的值最小，雙方的綜合損失費用也最小。因此，當前鼠浪湖碼頭需要增加1個泊位，這樣效益最佳，一年可降低費用6 347.42萬元。根據(jù)未來鐵礦石船型大型化趨勢，鼠浪湖碼頭近期可增加1個(10～30)萬t級的泊位。

3.2 遠期優(yōu)化

遠期優(yōu)化需要對鼠浪湖碼頭的遠期吞吐量進行預(yù)測。由于鼠浪湖碼頭營運的時間不久，無法得到足夠的鐵礦石吞吐量數(shù)據(jù)樣本來進行預(yù)測，這里擬根據(jù)整個舟山港域鐵礦石吞吐量預(yù)測及近幾年鼠浪湖碼頭吞吐量占比分析推算鼠浪湖碼頭鐵礦石的吞吐量預(yù)測值。

根據(jù)舟山港域鐵礦石碼頭2006-2018年鐵礦石吞吐量數(shù)據(jù)，利用二次指數(shù)平滑預(yù)測和灰度預(yù)測所構(gòu)建的組合模型，在假定未來幾年內(nèi)經(jīng)濟結(jié)構(gòu)與需求趨勢不發(fā)生較大改變的前提下，預(yù)測舟山港域鐵礦石2025年吞吐量約為22 743萬t(篇幅限制，詳細預(yù)測過程略)。

鼠浪湖鐵礦石碼頭自2016年建成營運以來，由于其定位與硬件條件，其吞吐量在整個舟山港域鐵礦石吞吐量的占比快速增加(表6)。而這幾年中，港域內(nèi)內(nèi)其他鐵礦石碼頭占比則主要在下降。

表6 2016-2018年鼠浪湖碼頭鐵礦石吞吐量占比Tab.6 Proportion of iron ore throughput of Shulanghu wharf in 2016-2018

根據(jù)鼠浪湖碼頭從2016年投運至今的歷年吞吐量與其它礦石碼頭歷年吞吐能力的平衡分析，考慮碼頭現(xiàn)有岸線與航道資源容量上限，結(jié)合專家評估意見和政府規(guī)劃，推算2025年鼠浪湖碼頭的鐵礦石吞吐量將占舟山港域總吞吐量35%左右，約為7 960萬t。

3.2.2 遠期優(yōu)化計算

已知當前泊位的年設(shè)計通過能力為5 200萬t·a-1，預(yù)測2025年的鐵礦石吞吐量為7 960萬t，產(chǎn)生2 800萬t左右的能力缺口，因此需要增加泊位。設(shè)每臺裝卸機械的裝卸效率依然是3000 t·h-1，船舶在本港的平均裝卸量參照2018年的40 866 t·艘-1，營運期N=365 d，此時的船舶到達率λ=5.34艘·d-1，以S=7，8，9，10，11來進行測算，計算結(jié)果如表7、表8所示。

表7 M/M/S模型最優(yōu)泊位數(shù)計算結(jié)果Tab.7 Optimal berth number calculation results by M/M/S model

表8 不同泊位數(shù)產(chǎn)生的泊位閑置損失費和船舶待泊損失費用計算結(jié)果Tab.8 Calculation results of berth idle loss and berth waiting loss incurred by different berths

結(jié)合表7、表8中的計算結(jié)果可以看出，當泊位數(shù)S=9時，發(fā)生在船、港雙方的綜合費用T和船、港方綜合損失費用Cb、Cs最小，鑒于現(xiàn)有泊位為5個，因此需要增加4個泊位。

根據(jù)未來鐵礦石船型大型化趨勢，遠期可考慮增加4個10～30萬t級的泊位，泊位年設(shè)計吞吐能力達到 8 000 萬 t。

4 結(jié)論

(1)鼠浪湖鐵礦石碼頭排隊系統(tǒng)中，船舶到港和占用泊位服務(wù)時間服從泊松分布和負指數(shù)分布，適用M/M/S模型。

(2)以綜合費用最小為優(yōu)化目標，經(jīng)優(yōu)化計算，當前鼠浪湖碼頭泊位數(shù)量布局不能滿足要求，需要增加1個泊位，這樣每年可節(jié)省6 347.42萬元，使本港的效益最大化。綜合考慮未來鐵礦石船型大型化趨勢，鼠浪湖碼頭近期可考慮增加1個(10～30)萬t級的泊位。

(3)在假定未來幾年舟山港域腹地經(jīng)濟結(jié)構(gòu)和需求趨勢不發(fā)生較大變化的情況下，預(yù)測2025年舟山港域鐵礦石吞吐量為約為22 743萬t，鼠浪湖碼頭的鐵礦石吞吐量將占舟山港域總吞吐量35%左右，約為7 960 萬 t。

(4)按上述預(yù)測值，鼠浪湖碼頭未來(2025年前)泊位數(shù)量布局需要進一步優(yōu)化，在綜合費用最小原則下其泊位數(shù)量為9個。根據(jù)鼠浪湖碼頭吞吐量預(yù)測以及鐵礦石運輸船型分析，2025年之前需增加4個(10～30)萬t級的泊位，碼頭年設(shè)計吞吐能力需達到8 000萬t，以滿足對鐵礦石運輸?shù)囊蟆?/p>