杜文龍，郭鳳祥，陳俐*

1 上海交通大學(xué) 動力裝置及其自動化研究所，上海 200240

2 上海交通大學(xué) 海洋工程國家重點實驗室，上海 200240

0 引 言

北極所蘊藏的豐富油氣儲量和礦產(chǎn)資源加速了各國開展極地運輸和科考活動的步伐[1-2]。作為關(guān)鍵設(shè)備，極地破冰船承擔(dān)著破冰護航、緊急救援的重要角色，但也面臨著能耗高、排放污染高和噪聲嚴(yán)重等問題[3]。目前，現(xiàn)役極地破冰船的動力系統(tǒng)以柴電推進系統(tǒng)為主[4-5]，但在破冰船航行或作業(yè)過程中，柴油發(fā)電機組須持續(xù)提供動力，所以無法在港口和極區(qū)水域?qū)崿F(xiàn)局域零排放或超低噪聲運行，這將對生態(tài)環(huán)境造成嚴(yán)重影響。

近年來，隨著儲能電池和柴油機混合動力系統(tǒng)的迅速發(fā)展，因其巨大的節(jié)能減排潛力，已在陸基微電網(wǎng)、新能源以及中小型船舶等領(lǐng)域開展了廣泛的研究工作[6-8]。然而，由于破冰船的極地航行周期長、環(huán)境惡劣，冰載荷變化頻繁且幅度較大[9]，目前尚無柴電混合動力系統(tǒng)優(yōu)化設(shè)計方面的文獻報告。

鑒于鋰電池能量密度和循環(huán)壽命的提高，及其在低溫工況下的研究成果問世，增加了其應(yīng)用于極地船舶的可行性[10-11]。針對低溫循環(huán)工況，劉中孝等[12]和馬瑞駿等[13]發(fā)現(xiàn)磷酸鐵鋰電池在–10 ℃低溫下的放電容量可達常溫電池容量的90.14%，在理論上可以滿足極地船舶的夏季能耗需求。如果將鋰電池應(yīng)用于極地破冰船的柴電混合動力系統(tǒng)，并利用儲能裝置在負載變化時的“削峰填谷”作用來提高柴油機的燃油效率，且支持零排放的純電推進航行模式，必將進一步提升極地破冰船的節(jié)能減排潛力。雖然大容量蓄電池可以支持長航程純電航行以及大功率電力輸出[14]，但其存在成本高、重量大等問題[15]。為了兼顧環(huán)境效益和極地破冰船的能耗與成本需求，須合理設(shè)計柴電混合動力系統(tǒng)的參數(shù)。

通過參考“雪龍”號科考船的航行路線，本文擬基于冰載荷變化等級建立阻力模型，并采用反向建模法建立柴電混合動力系統(tǒng)的能量流模型。鑒于破冰船混合動力系統(tǒng)的經(jīng)濟性和環(huán)保性，本文擬以極地破冰船的年油耗量和生命周期總成本為優(yōu)化目標(biāo)，以減速箱減速比、電機轉(zhuǎn)子體積、主機排量、螺旋槳直徑、電池組數(shù)量、最大荷電狀態(tài)以及最小荷電狀態(tài)等7 個參數(shù)為設(shè)計變量，采用快速非支配排序遺傳優(yōu)化算法（nondominated sorting genetic algorithm-II，NSGA-II）得 到Pareto優(yōu)化解集，并基于優(yōu)劣解距離（technique for order preference by similarity to an ideal solution，TOPSIS）綜合評價法從中得到最優(yōu)設(shè)計方案。最后，通過對比該最優(yōu)解與傳統(tǒng)柴電推進系統(tǒng)的年油耗量和生命周期總成本，用以驗證混合動力系統(tǒng)雙目標(biāo)優(yōu)化模型的合理性，同時分析優(yōu)化目標(biāo)對7 個設(shè)計參數(shù)的敏感性，相關(guān)研究成果可為混合動力系統(tǒng)的參數(shù)設(shè)計提供參考。

1 柴電混合動力系統(tǒng)模型

柴電混合動力系統(tǒng)和傳統(tǒng)柴電推進系統(tǒng)的架構(gòu)分別如圖1（a）和圖1（b）所示。傳統(tǒng)柴電推進系統(tǒng)包括4 臺柴油機、4 臺發(fā)電機、2 臺驅(qū)動電機、2 臺減速箱以及2 個螺旋槳，其中柴油機旋轉(zhuǎn)帶動發(fā)電機工作，將電力輸出至船舶電網(wǎng)，而驅(qū)動電機從船舶電網(wǎng)取電，并通過齒輪箱驅(qū)動螺旋槳。柴電混合動力系統(tǒng)則增加了1 套可充放儲能電池，用于存儲岸電或柴油發(fā)電機的多余電量，同時供給動力或用于滿足日用負載的功率需求。

圖1 動力系統(tǒng)架構(gòu)示意圖Fig. 1 Schematic diagram of propulsion system

1.1 破冰船負載及阻力模型

不同于常規(guī)船舶較小的航行負載變化量，極地破冰船既要滿足敞水區(qū)的長時間航行，還需滿足破冰區(qū)多種冰況下的負載功率需求。敞水區(qū)的負載特點為：變化范圍較小，周期長；而破冰區(qū)的負載特點為：阻力變化范圍大，受不同冰況碎冰尺寸的影響，其負載遠大于敞水區(qū)。參考“雪龍”號科考船的航行路線[16]，以及極區(qū)航行作業(yè)的復(fù)雜性，本文建立了破冰船單次航行周期的工況表，包括不同階段的冰況等級和對應(yīng)目標(biāo)航速，如表1 所示。根據(jù)冰厚及碎冰尺寸，可以將冰況分為0～5 級，用于計算不同海域的航行阻力[17-18]。其中：0 級表示無冰；1 級表示碎冰1 區(qū)（密集度為50%，尺寸為1.8 m）；2 級表示冰厚0.5 m；3 級表示冰厚1 m；4 級表示冰厚1.5 m；5 級表示碎冰2 區(qū)（密集度為80%，尺寸為2.5 m）。根據(jù)冰況變化曲線，本文還設(shè)定了單個周期的破冰船目標(biāo)航速曲線，如圖2 所示，整個航行周期為74 天，分為30 天敞水區(qū)I、14 天破冰區(qū)和30 天敞水區(qū)II 這3 個階段。

表1 破冰船單次航行周期工況Table 1 Operating profile of the icebreaker in single sailing cycle

圖2 單次航行周期的航速與冰級Fig. 2 Speed and ice scale in single sailing cycle

根據(jù)表1 的冰況等級，破冰船所受阻力可分為敞水區(qū)阻力和破冰區(qū)阻力。根據(jù)Holtrop 法，敞水區(qū)阻力的數(shù)學(xué)模型為

1） 當(dāng)冰級不大于3 級時，將采用Jeong 阻力公式：

1.2 柴油發(fā)電機組

本文將采用Willans Line 模型建立主機伸縮模型，僅關(guān)注柴油機能量流層面的功率輸入/輸出變化，而不涉及具體的參數(shù)尺寸和架構(gòu)。通過將參數(shù)化可伸縮模型與基于試驗的經(jīng)驗?zāi)Ｐ拖嘟Y(jié)合，可以大為簡化發(fā)動機建模的復(fù)雜度，其具體驗證過程參見文獻[19-20]。以某船用四沖程柴油機為基準(zhǔn)模型，其比油耗（specific fuel oil consumption，SFOC）圖譜如圖3 所示。圖3 中，紅線和藍線分別表示基準(zhǔn)主機的最大扭矩曲線和SFOC 等高線。Willans Line 模型的核心為引入平均制動有效壓力pme和 平均燃油有效壓力pmf，用以分別表征柴油機的輸出機械能和輸入化學(xué)能?；诎l(fā)動機的結(jié)構(gòu)參數(shù)和運行參數(shù)，pme和pmf的計算公式分別為

圖3 基準(zhǔn)主機的SFOC 圖譜Fig. 3 SFOC map of the baseline diesel engine

1.3 儲能鋰電池

鑒于鋰離子電池的高能量密度及其電壓曲線的平滑特性，本文將采用Rint 模型對鋰電池進行建模[21-22]，具體如下：

1.4 電機

1.5 螺旋槳

驅(qū)動電機和定距螺旋槳通過減速箱相連，其中減速箱的數(shù)學(xué)模型為

圖4 基準(zhǔn)電機的效率圖譜Fig. 4 Efficiency map of the baseline motor

1.6 船舶的縱向動力學(xué)模型

破冰船在航行過程中，其縱向受力模型為

2 雙目標(biāo)優(yōu)化

雖然配置儲能電池可以有效降低混合動力系統(tǒng)的油耗，減少系統(tǒng)排放，但電池的購置和維護將顯著增加系統(tǒng)成本，因此，本文的研究目標(biāo)是在滿足破冰船目標(biāo)航速和冰況變化的前提下，使破冰船混合動力系統(tǒng)的生命周期總成本和年油耗量實現(xiàn)最低，從而獲得良好的經(jīng)濟性和環(huán)境效益。本文將破冰船動力系統(tǒng)的參數(shù)設(shè)計轉(zhuǎn)化為雙目標(biāo)優(yōu)化問題，其優(yōu)化目標(biāo)為破冰船混合動力系統(tǒng)的生命周期總成本C和年油耗量mlife，具體船舶尺寸和模型參數(shù)如表2 所示。設(shè)計變量包括減速箱減速比i、電機轉(zhuǎn)子體積Vr、主機排量Vd、螺旋槳直徑D、電池組數(shù)量Numb、最大荷電狀態(tài)SOCmax以 及最小荷電狀態(tài)SOCmin等7 個參數(shù)，各優(yōu)化變量的取值范圍如表3 所示。

表2 船舶尺寸及模型參數(shù)Table 2 ship size and model parameters

表3 設(shè)計變量的變化范圍Table 3 variation range of design variables

優(yōu)化目標(biāo)函數(shù)f1，f2為

由此可見，破冰船的參數(shù)優(yōu)化設(shè)計被轉(zhuǎn)化為了一個具有約束條件的雙目標(biāo)、多約束耦合的非線性優(yōu)化問題。由于本文的目標(biāo)函數(shù)包含了動力學(xué)模型、能量管理策略等非線性環(huán)節(jié)，其數(shù)學(xué)模型較為復(fù)雜，故難以采用強理論保證的經(jīng)典數(shù)學(xué)優(yōu)化方法予以求解[23-24]。而遺傳算法、粒子群等基于進化的智能算法，可以通過特定的搜索機制來逼近最優(yōu)解，僅需利用目標(biāo)函數(shù)的信息且無需導(dǎo)數(shù)，其通用性較好，目前已在船舶、航天、汽車等領(lǐng)域得以廣泛應(yīng)用[25-27]。其中，非支配排序遺傳算法（NSGA-II）具有計算復(fù)雜度低、解集收斂性好等優(yōu)點[28]，適用于求解高度非線性問題，已被大量應(yīng)用于混合動力發(fā)電系統(tǒng)以及相關(guān)電力優(yōu)化模型[25,29]，所以本文將采用NSGA-II 算法進行柴電混合動力系統(tǒng)的雙目標(biāo)優(yōu)化設(shè)計。整個優(yōu)化流程如圖5 所示：首先，NSGA-II 算法通過交叉、變異產(chǎn)生父種群，即可得到更新后各設(shè)計變量的值；然后，將其代入混合動力系統(tǒng)模型中，根據(jù)目標(biāo)函數(shù)計算擁擠度；最后，結(jié)合擁擠度比較算子和精英策略并更新父種群，以此迭代得到最優(yōu)解集。圖5 中，Gen為進化代數(shù)。

圖5 優(yōu)化流程圖Fig. 5 Optimization flow chart

3 仿真及結(jié)果分析

3.1 最優(yōu)解集

在Matlab 中分別對傳統(tǒng)柴電動力系統(tǒng)和柴電混合動力系統(tǒng)進行建模，通過分析2 套系統(tǒng)方案的生命周期總成本和年油耗量的變化規(guī)律，從而選取最優(yōu)方案，如圖6 所示。

圖6 優(yōu)化結(jié)果Fig. 6 Optimization results

圖6（a）中的藍色點集為柴電混合動力系統(tǒng)的仿真結(jié)果，其年油耗量和總成本為負相關(guān)關(guān)系：總成本最大值為6.05×107元，此時年油耗量為3 098.83 t，為所有優(yōu)化方案的最低值，記為方案R；所有方案的最低成本為4.003 6×107元，低于方案R的66.18%，而此時年油耗量為3 172.14 t，比方案R增加了2.37%，記為方案Q。由此可見，柴電混合動力系統(tǒng)的儲能電池有助于減少破冰船的年油耗量，從而降低碳氮硫排放量；然而，電池的購買、更換、維修等將導(dǎo)致柴電混合動力系統(tǒng)方案的總成本有所升高，因此，柴電混合動力系統(tǒng)通過犧牲系統(tǒng)經(jīng)濟性而實現(xiàn)了年油耗量的降低。

本文將采用TOPSIS 法[30]來選取柴電混合動力系統(tǒng)的最優(yōu)方案，其基本思想為：當(dāng)無法獲得多個目標(biāo)屬性均為最優(yōu)的方案時，則通過評估所需方案與理想解、非理想解（圖6 中所示）的距離進行排序。當(dāng)評價方案點與理想解的距離最近，且同時與負理想解的距離最遠時，此方案即為最優(yōu)。通過對2 個優(yōu)化目標(biāo)進行歸一化處理，即可采用TOPSIS 法選擇混合動力系統(tǒng)Pareto 解集中的最優(yōu)方案，即圖6（a）中的方案B。

與之相比，傳統(tǒng)柴電推進系統(tǒng)的最優(yōu)解集為圖6（b），其年油耗量和總成本整體為正相關(guān)關(guān)系。成本變化來源于各方案的年油耗量差異，不同方案的設(shè)計變量值使得動力系統(tǒng)子部件無法同時處于最佳工作狀態(tài)，最終導(dǎo)致了油耗量增加。柴電推進系統(tǒng)的最優(yōu)方案即方案P（圖6（b）），其年油耗量為3 171.80 t，總成本為4.003 3×107元，為所有方案最小值。

將圖6（b）左側(cè)的總成本、年油耗量相對偏低的3 個方案點作于圖6（a），即圖6（a）右下方的紅色點集。對比柴電混合動力系統(tǒng)的優(yōu)化方案集和方案P，可知方案P 的年油耗量基本高于柴電混合動力系統(tǒng)的所有方案，且總成本趨近于圖6（a）的最低值，即方案Q。表4 所示為圖6（a）中4 種方案的仿真結(jié)果，與傳統(tǒng)柴電推進系統(tǒng)的最優(yōu)方案P 相比，柴電混合動力系統(tǒng)最優(yōu)方案B 的年油耗量為3 111.89 t，節(jié)約了1.89%的油耗。由于儲能電池的加入，方案B 支持純電航行模式，其中554.53 h的純電航行總里程占整個航程的31.22%，從而有利于提高船員舒適度。然而，儲能電池的購買和置換也增加了系統(tǒng)總成本，增幅達17.30%，所以降低了系統(tǒng)的經(jīng)濟性。

表4 設(shè)計方案的仿真結(jié)果對比Table 4 Comparison of simulation results of design schemes

3.2 參數(shù)敏感性分析

為了進一步定性判斷各設(shè)計變量對2 個目標(biāo)性能的影響，本節(jié)將開展參數(shù)敏感性分析，圖7所示為優(yōu)化目標(biāo)年油耗量和生命周期總成本隨著7 個參數(shù)的變化規(guī)律。

圖7（a）和圖7（b）分別為減速箱減速比和電機轉(zhuǎn)子體積對目標(biāo)函數(shù)的影響，未呈現(xiàn)明顯規(guī)律，整體波動較小，其中減速器減速比集中分布于31.5，電機轉(zhuǎn)子體積集中分布于50.1 L。

圖7（c）記錄了系統(tǒng)年油耗量和生命周期總成本隨柴油機排量的變化規(guī)律，其中系統(tǒng)年油耗量隨著發(fā)動機排量的增加而減小，而系統(tǒng)總成本則呈相反規(guī)律。發(fā)動機排量的下限值為59.03 L，這是因為柴油機需滿足破冰船航行過程的最大運行功率要求。

圖7（d）顯示破冰船年油耗量隨著電池組數(shù)量的增加呈下降趨勢，而系統(tǒng)總成本與電池組數(shù)量則近似為線性關(guān)系，這表明增加電池組數(shù)量來代替柴油機工作將有助于降低系統(tǒng)油耗量，但電池的購買和置換也導(dǎo)致了系統(tǒng)總成本升高。需注意的是，當(dāng)柴電混合動力系統(tǒng)方案的電池組數(shù)量趨于0 時，其年油耗量和總成本與圖6（b）中的方案P 近似相等，這也進一步驗證了本文仿真模型的合理性。

圖7（e）為2 個目標(biāo)函數(shù)關(guān)于電池荷電量最大值SOCmax和最小值SOCmin的變化規(guī)律，圖中點集顯示：對于年油耗量較低的方案集，SOCmax和SOCmin的區(qū)間范圍更大，更有利于電池容量的充分發(fā)揮，從而減少油耗量；對于總成本較高的方案集，SOCmax和SOCmin的區(qū)間范圍也偏大，與年油耗量呈現(xiàn)相反規(guī)律，這與3.1 節(jié)的分析結(jié)果一致。

由圖7（f）可知，螺旋槳直徑的優(yōu)化值穩(wěn)定于5 m，而優(yōu)化約束設(shè)定螺旋槳的直徑上限也是5 m（表2），所以推斷螺旋槳直徑對目標(biāo)函數(shù)基本無影響。

綜上所述，2 個優(yōu)化目標(biāo)（年油耗量和生命周期總成本）對發(fā)動機排量、電池組數(shù)量、最大荷電量以及最小荷電量4 個變量的敏感程度較高；而減速箱減速比、電機轉(zhuǎn)子體積和螺旋槳直徑這3 個變量對其影響則很小，可以簡化模型以降低后續(xù)優(yōu)化工作的計算成本。

4 結(jié) 論

針對復(fù)雜的極地航行工況，提出由柴油發(fā)電機組和儲能電池組成的破冰船柴電混合動力系統(tǒng)。本文采用反向建模法建立了柴電混合動力系統(tǒng)的能量流模型，然后以破冰船年油耗量及生命周期總成本為優(yōu)化目標(biāo)，基于快速非支配排序遺傳算法對動力系統(tǒng)部件進行了優(yōu)化分析，得到如下結(jié)論：

1） 柴電混合動力系統(tǒng)的最優(yōu)設(shè)計方案比傳統(tǒng)柴電推進系統(tǒng)節(jié)約了1.89%油耗；同時，支持純電航行模式，其中純電航行總里程占比31.22%，有利于實現(xiàn)局域低噪聲、零污染航行，但系統(tǒng)經(jīng)濟性則有所降低。

2） 參數(shù)敏感性的分析結(jié)果表明，混合動力系統(tǒng)的年油耗量和生命周期總成本對主機容量、電池組數(shù)量、電池荷電狀態(tài)邊界較敏感，而對減速器減速比、電機轉(zhuǎn)子體積和螺旋槳尺寸則相對不敏感，在后續(xù)設(shè)計研究中可簡化參數(shù)以降低模型的復(fù)雜度和計算量。