文／韓 鵬 馬海洋 馬 勇

針對(duì)中國航道船舶交通現(xiàn)狀，提出T形航路船舶交通通過能力的研究課題，建立T形航路船舶交通通過能力模型，最后借助計(jì)算機(jī)語言實(shí)現(xiàn)模型的求解。

In view of vessel traff i c situation in China, the article proposes a research subject about vessel transit capacity in T-type channel and sets up a transit capacity model. And results can be given by using computer language.

航道通過能力的研究主要分為以下三種：以單位時(shí)間內(nèi)通過的客運(yùn)人數(shù)/貨物噸數(shù)來表示基于吞吐量的航道通過能力；以航道所能通行的最大尺度的船型來表示基于單船的航道通過能力；以單位時(shí)間內(nèi)通過的船舶艘數(shù)來表示基于船舶交通的交通通過能力。其中航路船舶交通通過能力又分為單航路和多航路交通通過能力研究，單航路指單向或雙向航路，多航路包括T形及十字形航路等。

國外對(duì)航路船舶交通通過能力的研究起步較早，且大多以單、雙向航路船舶交通通過能力研究為主。對(duì)于多航路船舶交通，主要以日本稻石正明[1]等學(xué)者建立的船舶行為集群模型為代表，國內(nèi)主要以余勁[2]等船舶交通流的到船分布和船頭間距分布特性研究為代表。本文嘗試借鑒國內(nèi)外學(xué)者的先進(jìn)研究手段，展開T形航路船舶交通通過能力的研究。

一、T形航路船舶交通規(guī)律分析

1.排隊(duì)論的引入

A.k.Erlang于1909年創(chuàng)立了排隊(duì)論[3]。排隊(duì)是人們?nèi)粘Ｉ钪谐Ｒ姷默F(xiàn)象，多航路船舶交通亦是一種較為復(fù)雜的排隊(duì)現(xiàn)象。產(chǎn)生排隊(duì)現(xiàn)象的兩個(gè)主要因素是資源的有限性和人們對(duì)資源的需求。此外，服務(wù)機(jī)構(gòu)的管理水平、運(yùn)行效率等也會(huì)明顯影響排隊(duì)現(xiàn)象的產(chǎn)生。排隊(duì)論的研究目的就是尋找各種排隊(duì)系統(tǒng)的規(guī)律性，以調(diào)整和控制排隊(duì)系統(tǒng)，使之最大限度地發(fā)揮服務(wù)系統(tǒng)的效率。

2.T形航路船舶交通流概況

T形航路內(nèi)船舶交通情況復(fù)雜，一般分為以下幾種：圖1（1）為最簡單的航路船舶交通流分布情況，即簡化為兩支船舶交通流交匯情況，方向可向東或向西；圖1（2）為T形航路內(nèi)三支船舶交通流分布情況，有一支船舶交通流不受其他船舶交通流的干擾，航路內(nèi)有兩支船舶交通流交匯；圖2為本文所研究的T形航路基本交通流分布示意圖，在T形航路內(nèi)船舶為三支交通流交匯。

圖1 T形航路船舶交通流部分流向

圖2 T形航路內(nèi)船舶交通流分布

3.T形航路內(nèi)船舶到達(dá)規(guī)律

對(duì)于一組實(shí)測(cè)交通數(shù)據(jù)x1，x2，…，xn，一般滿足某種特定的分布。交通工程中常用χ2檢驗(yàn)加以驗(yàn)證，這種方法不論是連續(xù)型還是離散型分布都可應(yīng)用。本文采用χ2檢驗(yàn)擬合船舶到達(dá)規(guī)律，T形航路內(nèi)船舶到達(dá)為相互獨(dú)立事件，到達(dá)規(guī)律需在大量實(shí)測(cè)基礎(chǔ)上求得，一般服從相互獨(dú)立、同分布的Poisson分布，其分布的概率母函數(shù)、均值和方差分別是

4.T形航路內(nèi)船舶服務(wù)規(guī)律

T形航路交叉區(qū)域內(nèi)的船舶服務(wù)一般為單船服務(wù)，且為三支交通流循環(huán)服務(wù)，某一船舶在交叉區(qū)域駛過的距離除以其行駛速度即為該船舶在交叉區(qū)域的服務(wù)時(shí)間。服務(wù)的時(shí)間為一隨機(jī)變量，其分布的概率母函數(shù)、均值和方差分別為當(dāng)T形航路為某一交通流內(nèi)單船完成服務(wù)、轉(zhuǎn)向?yàn)槠渌?duì)列的船舶服務(wù)時(shí)，會(huì)出現(xiàn)一定時(shí)間的服務(wù)間歇，即交叉區(qū)域服務(wù)轉(zhuǎn)換時(shí)間，該隨機(jī)變量分布的概率母函數(shù)、均值和方差分別是本文結(jié)合T形航路內(nèi)船舶交通的實(shí)際特點(diǎn)，求得服務(wù)率β為

式中：k為船舶列數(shù)(也代表T形航路服務(wù)窗口數(shù)，取k=1)；m為船舶類型數(shù)； Vi為第i類船舶的船速；pi為第i類船舶出現(xiàn)的概率；Si為第i類船舶駛過交叉區(qū)域的通行距離。

二、T形航路船舶交通通過能力模型的建立

1.服務(wù)的時(shí)間Vi, m

假定在某一交通流服務(wù)期間（τm＝tm-tm-1）有ζi(τm)艘船舶進(jìn)入第i列交通流，記其中的第l艘船舶（l＝1,2,…,ζi(τm)）接受服務(wù)需要時(shí)間為Si,l，而在這Si,l期間內(nèi)又有船舶進(jìn)入第i列交通流。經(jīng)計(jì)算得，船舶進(jìn)入第i列交通流接受服務(wù)的時(shí)間Vi,m為

2.第 j 列交通流中平均船舶數(shù)

記μi(n)為交叉區(qū)域從第i列向第i+1列交通流服務(wù)轉(zhuǎn)換的時(shí)間，μj(μi)，ηj(vi)分別為在μi，vi時(shí)間內(nèi)進(jìn)入第j列交通流的船舶數(shù)量，則在tn+1時(shí)刻有

設(shè)在tn時(shí)刻第i列交通流開始接受服務(wù)，經(jīng)計(jì)算得第j列交通流中平均船舶數(shù)為

3.船舶接受服務(wù)的過程

若第i列交通流在tn時(shí)刻駛?cè)氩⒃趖n1時(shí)刻駛離交叉區(qū)域，由于在tn1時(shí)刻仍然有船舶到達(dá)第i列交通流，因此，在tn1時(shí)刻新進(jìn)入第i列交通流的船舶艘數(shù)為εi(n1)，而其駛離交叉區(qū)域的時(shí)刻為tn2，船舶接受服務(wù)的過程如圖3所示：

圖3 第i列船舶交通流接受服務(wù)的過程

4.T形航路內(nèi)船舶交通通過能力N

若第i列交通流中第l艘船舶在第τm時(shí)間段（τm＝τmτm-1）進(jìn)入該交通流，并在 時(shí)刻接受服務(wù)，則第l艘船舶等待接受服務(wù)的時(shí)間為

在此期間，有關(guān)系式：

式中：d(n)為第i列交通流在時(shí)刻tn再次接受服務(wù)的時(shí)間；ai,m(n)為在τm時(shí)間內(nèi)進(jìn)入第i列交通流的船舶數(shù)量。

同理，在第i列交通流接受服務(wù)期間，新進(jìn)隊(duì)列船舶等待接受服務(wù)時(shí)，有關(guān)系式:

由式 (6)、(7)計(jì)算得wi,1(n)的平均值，同理，可計(jì)算得wi,2(n)的平均值，由此，船舶在交通流中接受服務(wù)的時(shí)間為

在系統(tǒng)中，船舶到達(dá)服從泊松分布，系統(tǒng)中所能容納的船舶數(shù)量為

三、模型的實(shí)現(xiàn)

對(duì)于T形航路交通通過能力的計(jì)算，本文通過VC++編程加以實(shí)現(xiàn)。利用本模型可計(jì)算在一定約束性條件及相關(guān)航行參數(shù)下的T形航路船舶的交通通過能力。T形航路船舶交通通過能力求解流程如圖4所示。

依據(jù)本模型可求解不同船型（標(biāo)準(zhǔn)、中等、小型）在不同通過參數(shù)情況下對(duì)應(yīng)的T形航路船舶日服務(wù)艘數(shù)。如對(duì)于1000 t級(jí)中型船舶，若橫流速度為0.33 m/s，可得出該T形航路日船舶交通通過能力為152艘；選取5000 t級(jí)標(biāo)準(zhǔn)船型，則在不同的設(shè)定參數(shù)（橫流分別為0.755 m/s、0.13 m/s）下對(duì)應(yīng)的T形航路日船舶交通通過能力分別為110艘和124艘。經(jīng)驗(yàn)證，運(yùn)算結(jié)果與目前我國主要T形航路理論日服務(wù)艘數(shù)基本一致。另外，通過模型也可驗(yàn)證：T形航路服務(wù)水平越高，橫流越小，船舶交通通過能力越大。

圖4 T形航路船舶交通通過能力求解流程

[1] 稻石正明, 今宿智仁. 船舶群と船舶群遭遇シミユレ—シヨン[C].日本航海學(xué)會(huì)論文集, 2003:21-26.

[2] 余勁. 西江航道船舶流的概率分布特性[J]. 交通運(yùn)輸工程學(xué)報(bào),2006(2):19-20.

[3] 陸鳳山. 排隊(duì)論及其應(yīng)用[M]. 長沙:湖南科學(xué)技術(shù)出版社,1984:303-321.