基于航次營(yíng)運(yùn)效益的不定期船航速優(yōu)化研究

茆沐嘉,王 琪,馬來(lái)好,陸 萍

(1.江蘇航運(yùn)職業(yè)技術(shù)學(xué)院 輪機(jī)工程學(xué)院,江蘇 南通 226010;2.大連海事大學(xué) 輪機(jī)工程學(xué)院,遼寧 大連116026)

0 引言

在經(jīng)濟(jì)全球化的背景下,海運(yùn)已成為當(dāng)前國(guó)際商業(yè)貿(mào)易中重要的運(yùn)輸方式。站在全球航運(yùn)公司的角度,追求更高的船舶運(yùn)輸效益與盈利是其永恒的目標(biāo)。當(dāng)簽訂了航次租船合同并確定了相應(yīng)運(yùn)營(yíng)的船舶后,能夠直接影響船舶單航次的營(yíng)運(yùn)效益、并且可供航運(yùn)公司決策的只有航速這一個(gè)變量。一方面,一旦降低航速,航次的總時(shí)間會(huì)增加,可能導(dǎo)致船舶錯(cuò)過(guò)市場(chǎng)某些運(yùn)輸任務(wù),從而失去相應(yīng)的盈利機(jī)會(huì);另一方面,增加航速,則船舶主機(jī)的燃油消耗會(huì)急劇增長(zhǎng),致使船東難以承受大幅加劇的燃油成本[1]。經(jīng)過(guò)調(diào)研,目前國(guó)內(nèi)航運(yùn)公司調(diào)度人員主要借助傳統(tǒng)的excel表格計(jì)算并結(jié)合個(gè)人經(jīng)驗(yàn)對(duì)航速進(jìn)行決策,但此方法并沒(méi)有對(duì)航次整體經(jīng)濟(jì)效益進(jìn)行系統(tǒng)性的綜合考慮,缺乏科學(xué)性和準(zhǔn)確性。

國(guó)際上：FAGERHOLT等[2]以航次燃油消耗最小為優(yōu)化目標(biāo),建立了船舶航速優(yōu)化模型,并利用啟發(fā)式算法進(jìn)行計(jì)算;STALHANE等[3]借助分支定界法對(duì)港口裝卸時(shí)間約束下的不定期船運(yùn)輸問(wèn)題進(jìn)行求解;CHANG等[4]通過(guò)研究船舶降速與航運(yùn)成本間的關(guān)系,指出租船費(fèi)率和燃油價(jià)格對(duì)最優(yōu)航速的決策起到?jīng)Q定性的作用。國(guó)內(nèi)：中國(guó)船級(jí)社(CCS)修訂發(fā)布最新版的《智能船舶規(guī)范》(2020)把航速優(yōu)化作為智能船舶管理的重要組成部分[5];俞超等[6]建立了考慮滯期費(fèi)和速遣費(fèi)影響的不定期船航速優(yōu)化模型,并通過(guò)兩階段粒子群算法進(jìn)行求解;馬來(lái)好等[7]建立了船舶航次營(yíng)運(yùn)效益模型,通過(guò)實(shí)例對(duì)不同油價(jià)與運(yùn)價(jià)情形進(jìn)行計(jì)算。以上學(xué)者對(duì)于海船最優(yōu)航速問(wèn)題進(jìn)行了研究,建立了基于不同優(yōu)化目標(biāo)的數(shù)學(xué)模型,并提出了相關(guān)的優(yōu)化方法,但模型缺乏考慮船舶航速對(duì)航次凈收益、航次日收益、航次時(shí)間的影響和貨主的運(yùn)輸服務(wù)時(shí)間要求等。

本文針對(duì)不定期單船、單航次、單貨物的運(yùn)營(yíng)特征,提出科學(xué)的航速優(yōu)化方法,利用自主開發(fā)的優(yōu)化軟件基于優(yōu)化算法對(duì)模型進(jìn)行求解,實(shí)現(xiàn)不定期船的航次凈收益最大化或航次日收益最大化,以此提高船舶航次經(jīng)濟(jì)效益。

1 問(wèn)題描述

(1)不定期船運(yùn)輸?shù)臓I(yíng)運(yùn)核算是以航次為基礎(chǔ)。典型航次指船舶從始發(fā)港開航,航次在離港出發(fā)的那一刻即開始,接著壓載(空載狀態(tài))航行至合同指定的裝貨港裝貨,?？垦b貨港期間將一票貨物一次裝滿,然后裝載(滿載狀態(tài))航行至卸貨港卸貨,并在?？啃敦浉燮陂g一次卸完,本航次在卸完貨的那一刻結(jié)束。另外,航次中可能途徑加油港加油。

(2)考慮船舶航行安全性要求和主機(jī)的極限功率等因素,航速的選擇要在可行區(qū)間內(nèi)。

(3)在單航次運(yùn)輸過(guò)程中,船舶僅以壓載(空載)航速和裝載(滿載)航速定速航行。

(4)在海上航行期間,忽略外部環(huán)境因素及船舶自身污底等不確定性因素對(duì)船舶主機(jī)油耗和航速的影響[8]。

(5)船舶主機(jī)在低硫油限制區(qū)內(nèi)使用低硫油,區(qū)域外使用重油。

2 基于航次營(yíng)運(yùn)效益的數(shù)學(xué)模型

2.1 航次收入

Iv=(frq+dem-des)(1-cmr-r)

(1)

式中：Iv為航次凈收入,美元;fr為運(yùn)費(fèi)率(運(yùn)價(jià)),美元/t;q為載貨量,t;dem為滯期費(fèi),美元;des為速遣費(fèi),美元;cmr為傭金比例;r為營(yíng)業(yè)稅率。

2.2 航次時(shí)間

航次時(shí)間包括船舶壓載航行時(shí)間、裝載航行時(shí)間、裝貨港的停靠時(shí)間、卸貨港的停靠時(shí)間、加油港停靠時(shí)間、錨泊時(shí)間。

(2)

式中：tv為航次總時(shí)間,d;lb為壓載(空載)航程,n mile;ll為裝載航程,n mile;vb為壓載航速,kn;vl為裝載航速,kn;tp,l為船舶在裝貨港?？繒r(shí)間,即裝貨時(shí)間,d;tp,d為船舶在卸貨港?？繒r(shí)間,即卸貨時(shí)間,d;tp,o為船舶在加油港?？繒r(shí)間,d;tb為船舶航次中的錨泊時(shí)間,d。

2.3 航次成本

航次成本分為固定成本和變動(dòng)成本兩部分：固定成本指船員費(fèi)、潤(rùn)料費(fèi)、備件費(fèi)、折舊費(fèi)、修理費(fèi)、保險(xiǎn)費(fèi)、管理費(fèi)等不受運(yùn)輸貨量變化影響而相對(duì)保持固定不變的費(fèi)用;變動(dòng)成本包括主機(jī)、發(fā)電機(jī)、鍋爐燃油成本及港口使費(fèi)。

Cv=Cm+Cg+Cb+Cp+Co

(3)

式中：Cv為航次變動(dòng)成本,美元;Cm為主機(jī)燃油成本,美元;Cg為發(fā)電機(jī)燃油成本,美元;Cb為鍋爐燃油成本,美元;Cp為港口使費(fèi),美元;Co為其他費(fèi)用,美元。

主機(jī)燃油成本包括壓載航行和裝載航行主機(jī)燃油成本。若航線涉及低硫油限制區(qū),該成本也包括低硫油燃油費(fèi)用。假設(shè)主機(jī)重油油耗和低硫油油耗水平相同。

(4)

式中：lb,ls為壓載航程中涉及低硫油區(qū)航程,n mile;ll,ls為裝載航程中涉及低硫油區(qū)航程,n mile;pf為重油價(jià)格,美元/t;pls為低硫油價(jià)格,美元/t;qm,b、qm,l分別為壓載、裝載航行中的主機(jī)每天的燃油消耗量,t/d。

船舶主機(jī)每日油耗近似與船舶航速的三次方成正比[9]。本研究通過(guò)對(duì)實(shí)測(cè)的船舶主機(jī)油耗數(shù)據(jù)進(jìn)行曲線擬合的方式,得到船舶在壓載狀態(tài)、裝載狀態(tài)下的每日油耗和航速的關(guān)系。

(5)

(6)

發(fā)電機(jī)燃油成本Cg包括航行天、靠港期間、錨泊期間燃油成本,其中發(fā)電機(jī)靠港期間油耗與錨泊期間油耗近似相等。

(7)

式中：qg,v、qg,m分別為航行天和錨泊天發(fā)電機(jī)的每日油耗,t/d。

由于海上航行船舶僅使用廢氣鍋爐,因此鍋爐燃油消耗Cb主要發(fā)生在船舶停航階段。

Cb=pfqb(tp,l+tp,d+tp,o+tb)

(8)

式中：qb為鍋爐每日油耗,t/d。

港口使費(fèi)Cp包括航線涉及的裝、卸港口和掛靠加油港所產(chǎn)生的費(fèi)用。

(9)

式中：Cp,i為航次掛靠的裝、卸港口及可能掛靠加油港的港口使費(fèi),美元。

為保證模型的準(zhǔn)確性,在航次成本中增加一項(xiàng)其他費(fèi)用Co,用以調(diào)整和修正航次成本的估算。

2.4 目標(biāo)函數(shù)

綜上,對(duì)于不定期船,建立2個(gè)關(guān)于航次營(yíng)運(yùn)效益不同優(yōu)化目標(biāo)的數(shù)學(xué)模型,分別以航次凈收益最大和航次每日凈收益(日收益)最大為優(yōu)化目標(biāo)。

Cf

(11)

約束條件：

(1)抵港時(shí)間窗(受載期)約束,要求船舶必須在合同規(guī)定的時(shí)間區(qū)間內(nèi)抵達(dá)裝貨港：

Tl,min≤t1≤Tl,max

(12)

式中：t1為船舶實(shí)際抵達(dá)裝貨港的時(shí)間,d;Tl,min、Tl,max分別為合同規(guī)定的船舶抵達(dá)裝貨港的最早時(shí)間、最晚時(shí)間,d。

(2)定義船舶實(shí)際抵達(dá)裝貨港的時(shí)間：

(13)

式中：tm,b為壓載航程錨泊時(shí)間,d。

(3)限定船舶分別在壓載、裝載航行狀態(tài)下的航速要在其可行區(qū)間內(nèi)：

vb,min≤vb≤vb,max

(14)

vl,min≤vl≤vl,max

(15)

式中：vb,min、vb,max分別為船舶在壓載航行狀態(tài)下限定的最低航速和最高航速,kn;vl,min、vl,max分別為船舶在裝載航行狀態(tài)下限定的最低航速和最高航速,kn。

3 模型求解

2個(gè)數(shù)學(xué)模型的目標(biāo)函數(shù)、約束條件中均含有關(guān)于航速的非線性函數(shù)關(guān)系。因此,上述模型實(shí)質(zhì)是以壓載航速、裝載航速為自變量的具有非線性約束條件的非線性數(shù)學(xué)規(guī)劃問(wèn)題。根據(jù)2個(gè)決策變量之間的函數(shù)關(guān)系及目標(biāo)函數(shù)、約束條件的特征,同時(shí)考慮易于計(jì)算機(jī)編程,本研究選擇模式搜索法結(jié)合外罰函數(shù)法設(shè)計(jì)求解算法,利用VS作為開發(fā)平臺(tái),C#為編程語(yǔ)言,自主設(shè)計(jì)開發(fā)船舶營(yíng)運(yùn)效益優(yōu)化軟件對(duì)模型進(jìn)行求解。

4 算例分析

4.1 實(shí)船算例

選取某型散貨船作為目標(biāo)船,以其主流航線作為計(jì)算實(shí)例,從青島港始發(fā),途徑新加坡港進(jìn)行加油,壓載航行至巴西圖巴朗裝載鐵礦石,裝載航行回青島港卸貨,見(jiàn)圖1。

利用MATLAB對(duì)實(shí)測(cè)的船舶主機(jī)油耗數(shù)據(jù)進(jìn)行數(shù)學(xué)擬合得到目標(biāo)船的主機(jī)油耗,見(jiàn)圖2。

(16)

(17)

圖1 目標(biāo)船主流航線

圖2 目標(biāo)船主機(jī)油耗擬合曲線

算例涉及的船舶、港口、航線、貨物等具體參數(shù),見(jiàn)表1。

表1 算例具體參數(shù)

4.2 營(yíng)運(yùn)效益分布

將表1中數(shù)據(jù)代入所建立的2個(gè)基于航次營(yíng)運(yùn)效益的數(shù)學(xué)模型(不考慮受載期),利用MATLAB軟件得到目標(biāo)船航次凈收益和航次日收益在上述可行壓載、裝載航速范圍內(nèi)的分布情況,見(jiàn)圖3、圖4。從圖中可以看出,目標(biāo)船的航次凈收益較高值主要集中在其壓載、裝載航速都較低的區(qū)域,但并不是在航速最低值取得最大航次凈收益。整體上,航次凈收益隨著壓載航速、裝載航速的降低呈現(xiàn)迅速增加后緩慢降低的特點(diǎn);而目標(biāo)船的航次日收益較高值主要集中在壓載、裝載航速較高的區(qū)域,但不是在航速最高時(shí)取得最大日收益,并隨著壓載航速、裝載航速的降低呈現(xiàn)緩慢增加后迅速降低的特點(diǎn)。

圖3 目標(biāo)船航次凈收益隨航速變化分布圖

圖4 目標(biāo)船航次日收益隨航速變化分布圖

4.3 營(yíng)運(yùn)效益優(yōu)化對(duì)比

基于上述參數(shù),選擇受載期為32～37 d,利用開發(fā)的船舶營(yíng)運(yùn)效益優(yōu)化軟件進(jìn)行優(yōu)化計(jì)算。將優(yōu)化結(jié)果與正常營(yíng)運(yùn)航速下的相應(yīng)效益進(jìn)行對(duì)比,結(jié)果見(jiàn)表2。從表中可以看出,以航次凈收益最大為優(yōu)化目標(biāo)的壓載航速、裝載航速均小于正常營(yíng)運(yùn)航速,相應(yīng)航次時(shí)間高于正常營(yíng)運(yùn)時(shí)的航次時(shí)間,但航次凈收益比其提高了3.6%;而以航次日收益最大為優(yōu)化目標(biāo)的壓載航速、裝載航速卻均略高于船舶的正常營(yíng)運(yùn)航速,相應(yīng)航次時(shí)間也略低于正常營(yíng)運(yùn)時(shí)的航次時(shí)間,但日收益比其提高了0.53%。

表2 目標(biāo)船算例優(yōu)化對(duì)比結(jié)果

4.4 主要變動(dòng)因素分析

在上述算例基礎(chǔ)上,通過(guò)改變?nèi)加蛢r(jià)格進(jìn)一步比較分析優(yōu)化航速及相應(yīng)的營(yíng)運(yùn)效益,見(jiàn)表3。

橫向?qū)Ρ葋?lái)看,無(wú)論其以航次凈收益最大或航次日收益最大為優(yōu)化目標(biāo),當(dāng)燃油價(jià)格上升,對(duì)應(yīng)的船舶燃油成本也隨之增加,經(jīng)過(guò)航速優(yōu)化后的航次最大凈收益和航次最大日收益均下降。據(jù)此,建議船公司在攬貨談判,特別是長(zhǎng)期合同,應(yīng)將燃油價(jià)格的變動(dòng)機(jī)制在運(yùn)價(jià)的制定上予以體現(xiàn),以保證營(yíng)運(yùn)效益。在上述燃油價(jià)格增加幅度內(nèi),目標(biāo)船以航次凈收益最大為優(yōu)化目標(biāo),其優(yōu)化壓載航速保持不變,優(yōu)化裝載航速隨著燃油價(jià)格的增加而降低;以日收益最大為優(yōu)化目標(biāo),其優(yōu)化壓載航速和裝載航速均隨著燃油價(jià)格的增加而降低。2種優(yōu)化目標(biāo)分別對(duì)應(yīng)的優(yōu)化裝載航速均略小于優(yōu)化壓載航速,這與目標(biāo)船的營(yíng)運(yùn)收益的航速分布特點(diǎn)一致。

表3 不同燃油價(jià)格下的航次凈收益和

縱向?qū)Ρ葋?lái)看,在相同燃油價(jià)格下,以航次凈收益最大為優(yōu)化目標(biāo)取得的營(yíng)運(yùn)效益大于以航次日收益最大為優(yōu)化目標(biāo)取得的等效航次效益,但前者相對(duì)應(yīng)的航次時(shí)間始終大于后者。因此究竟選擇以何種目標(biāo)優(yōu)化后的航速進(jìn)行營(yíng)運(yùn),需充分兼顧航運(yùn)市場(chǎng)的供需：當(dāng)市場(chǎng)相對(duì)低迷,船多貨少,船期較為寬裕時(shí),建議航運(yùn)公司以航次凈收益最大為目標(biāo)的最優(yōu)航速進(jìn)行營(yíng)運(yùn);當(dāng)市場(chǎng)相對(duì)火爆,供不應(yīng)求,船期較為緊張,建議公司以航次日收益最大為目標(biāo)的最優(yōu)航速進(jìn)行營(yíng)運(yùn)。

5 結(jié)論

(1)本研究建立了以航次凈收益最大、航次日收益最大為優(yōu)化目標(biāo)的不定期船航速優(yōu)化數(shù)學(xué)模型,全面考慮了航速對(duì)航次凈收益、航次日收益、航次時(shí)間的影響及貨主的運(yùn)輸服務(wù)時(shí)間要求等。

(2)利用自主開發(fā)的優(yōu)化軟件基于優(yōu)化算法對(duì)模型進(jìn)行求解,能夠?qū)崿F(xiàn)不定期船的航次凈收益最大化或航次日收益最大化。

(3)通過(guò)實(shí)船算例分析,驗(yàn)證了本研究提供的不定期船最優(yōu)航速比航運(yùn)公司的指導(dǎo)航速取得更大的航次凈收益或航次日收益,能夠?yàn)楹竭\(yùn)公司在最優(yōu)航速?zèng)Q策領(lǐng)域提供輔助支持。