2. 武汉船用电力推进装置研究所,湖北 武汉 430064
2. Wuhan Institute of Marine Electric Propulsion, Wuhan 430064, China
容量配置会影响船舶动力系统的性能,其可以分为单目标和多目标2种,而单目标优化存在一定局限性:所选权重缺乏明确的物理依据,不同权重设定可能导致最优解结果出现较大偏差,所以进行多目标优化具有重要意义。在以往的研究中,Hatata等[1]提出了一种基于人工免疫系统的克隆进化算法;Xiao等[2]则针对风能、光伏、储能及柴油发电的容量配置问题,构建了一个多能源系统模型,通过设定权重系数将经济性与环保性目标合并,并采用水循环单目标优化算法进行求解;Eldeeb等[3]利用NSGA-II算法在最小化混合储能系统(HESS)成本、重量和体积的同时,最大化超级电容剩余循环寿命;冯龙祥等[4]将电池组与超级电容结合,以全生命周期行驶里程和动力系统成本为优化目标,通过NSGA-II算法获得了Pareto前沿解,即帕累托最优解构成的集合。
本文针对复合储能系统的容量配置难题,将锂电池与超级电容的串并联数量作为决策变量,并设定范围,构建以系统能量消耗与质量特性为优化目标、瞬时功率平衡、额定功率等为约束的模型,并采用改进型麻雀搜索算法求解,再通过优劣解距离法(TOPSIS)法进行多属性决策筛选最优配置参数,即距离正理想解最近,距离负理想解最远。该方法在维持动力性能的同时,能显著降低系统质量指标与能量损耗。
1 混合动力船复合储能系统 1.1 研究对象混合动力系统在功率波动大、航程短、排放受限的中小型船舶中优势明显。以“Alsterwasser”号渡轮为例,其主要动力来源于2台48 kW燃料电池,氢气储存于12个35 MPa高压储氢罐中,配有额定功率100 kW的推进电机和7组540 V、360 Ah的铅酸电池作为辅助动力源。
图1为该渡轮在典型工况下的负载曲线。船舶典型工况分为三阶段:经济巡航(0~90 s),功率稳定在45 kW;靠港(90~135 s),负载波动剧烈,峰值达100 kW;离港(160~190 s),功率升至110 kW。全程功率需求范围为0~110 kW,平均功率约42 kW。
|
图 1 母型船典型工况 Fig. 1 Typical operating conditions of the ship |
在复合储能系统的设计中,拓扑结构的选型不仅决定了能量传输效率和系统成本,还直接影响控制策略的实现难度。为此,本研究在原船动力系统模型中集成了复合储能模块,其拓扑形式须依照具体应用场景进行针对性设计。主流拓扑可按所配置的 DC/DC 变换器数量分为三类:被动式、半主动式和主动式,其能量流通路径及控制自由度依次增加。
综合权衡系统配置成本、能效表现及船舶运行需求,本文最终选用配备单级 DC/DC 变换器的半主动式拓扑,并在此基础上引入超级电容单元。该方案既可利用超级电容的高功率密度特点以应对起停和快速负载波动时的瞬态需求,又可通过精简 DC/DC 变换器数量来降低硬件投入及维护复杂性,兼顾了经济性与性能,实现结构与控制的优化平衡。
改进后的混合动力船动力系统结构如图2 所示。该系统集成了燃料电池、锂离子电池和超级电容3种能源单元,并分别通过2个独立的 DC/DC 升压转换器对燃料电池与超级电容进行功率调节,以期实现在多动力源并行供电下对负载需求的动态响应与高效协同。
|
图 2 混合动力船动力系统结构 Fig. 2 Hybrid propulsion system architecture for ships |
本文能量管理优化方法为等效最小氢耗策略,其核心在于通过优化算法在满足多种航行任务与约束条件的前提下,实现全时段氢气消耗最小化。该策略通过将燃料电池、超级电容与电池的能量输出转换为等效氢耗指标,计算系统在每个时刻的等效氢气消耗,并以此作为优化目标,此策略不仅有效降低了氢气消耗,也提升了复合储能系统在复杂工况下的运行效率。
2.1 ECMS原理将储能系统消耗和被补偿的能量转换成等效氢气消耗,与燃料电池的实际氢耗相加,将此作为优化问题所对应的目标值,每一瞬间等效氢气的消耗量为:
| $ \min C\left(t\right)={C}_{fc}+{k}_{1}\cdot {C}_{bat}+{k}_{2}\cdot {C}_{sc}。$ | (1) |
式中:
为了实现等效氢耗的实时最小化,本策略依据船舶当前航行工况及各动力源的运行状态动态计算等效氢气消耗,引入PI控制器对超级电容器的SOC进行实时监测,并将其SOC变化所对应的功率需求等效映射至锂电池能耗中,优化问题中的氢气消耗模型被修正为:
| $ \min C\left(t\right)={C}_{fc}+k\cdot {C}_{bat}。$ | (2) |
式中:
| $ {C}_{fc}\left(t\right)={\int }_{0}^{t}\left(\frac{{V}_{lpm}\cdot {P}_{a}}{R\cdot T}\cdot {M}_{{\mathrm{H}}_2}\right){\mathrm{d}}t。$ | (3) |
式中:
锂电池的等效氢气消耗
| $ {C}_{bat}\left(t\right)={P}_{bat}\cdot \sigma \cdot \frac{{C}_{fc.ave}}{{P}_{fc.ave}} 。$ | (4) |
其中,
| $ \sigma =\left\{\begin{aligned}&\frac{1}{{\eta }_{dis}\cdot {\eta }_{ch,ave}},{P}_{bat}\geqslant 0,\\& {\eta }_{ch}\cdot {\eta }_{dis,ave},{P}_{bat} \lt 0。\end{aligned}\right. $ | (5) |
式中:
本文基于模糊逻辑理论构建了等效因子计算模型,以应对复合储能系统在不同运行状态下的模糊性与不确定性。在控制器设计中,通过锂电池SOC和负载功率波动幅度与等效因子之间的模糊隶属函数及规则库,实现对系统能量分配的自适应调整。本文选择将负载
将锂电池SOC作为第一个输入变量,对应模糊子集为
|
图 3 输入及输出的隶属度函数 Fig. 3 Membership functions for inputs and outputs |
图中,横坐标表示变量的取值,也就是输入值,纵坐标表示对应横坐标值的隶属度。当锂电池SOC较低时,模糊控制器输出
输入为锂电池SOC和归一化后的负载功率
储能容量的合理设定在系统设计中仍面临显著挑战:储能系统的容量配置若过大,可能降低运行效率并增加系统成本;若过小,则可能无法满足动力需求,影响整体性能[5]。因此提出一种多目标容量优化模型,采用麻雀搜索算法(Sparrow Search Algorithm,SSA),结合多属性决策分析评估优化效果。最后通过分析动力性和经济性指标验证优化的有效性。
3.1 决策变量在本研究的多目标优化模型中,将超级电容并联个数
| $ {E}_{{bat},\mathrm{m}\mathrm{i}\mathrm{n}}=\mathrm{\Delta }E/\left[{\eta }_{{{bat}}}{\eta }_{dc}\right(SO{C}_\text{max}-SO{C}_\text{min}\left)\right],$ | (6) |
| $ {E}_{sc,\text{min}}=t\cdot {P}_{f.\text{max}}/{\eta }_{sc}。$ | (7) |
式中:
在混合动力船舶的性能评估中,动力性与经济性构成了2个核心维度,其表现不仅取决于航行时长、环境条件和载重等基本因素,还会受到具体工况与气温变化的影响。本文假设全工况时间、负载功率波动区间以及气象条件等外部参数在研究过程中保持不变,因而将降低系统整体能耗视为提升经济性的主要手段;与此同时,通过对典型工况下储能容量配置的优化,可在既定容量下延长连续运行时间,从而在确保动力性能的前提下,实现经济效益的最大化。
1)复合储能系统能量消耗
复合储能系统的能量消耗指标为:式中:
| $ \text{min}{E}_{c}={\sum }_{k=1}^{\mathrm{{T}}}\left[\Delta {E}_{bat}\right(k)+\Delta {E}_{sc}(k\left)\right] ,$ | (8) |
| $ \Delta {E}_{sc}\left(k\right)=0.5\cdot {C}_{sc}\cdot \left[{U}_{sc}^{2}\right(k)-{U}_{sc}^{2}(k-1\left)\right],$ | (9) |
| $ \Delta {E}_{bat}\left(k\right)=\left|{P}_{bat}\left(k\right)\right|\cdot \Delta t。$ | (10) |
式中:
2)复合储能系统质量特性优化
复合储能系统全体质量表达式为:
| $ {G}_\text{total}={G}_{fc}+{G}_{bat}+{G}_{sc}。$ | (11) |
式中:
1)瞬时功率平衡约束,动力系统需实时满足船舶动力需求,保证供需平衡,瞬时功率平衡约束表达式为:
| $ {P}_{bat}+{P}_{sc}+{P}_{fc}={P}_\text{load}。$ | (12) |
2)SOC约束,锂电池和超级电容过充过放会加速老化,影响使用寿命,应根据实际情况对荷电状态上下限进行约束,如下式:
| $\left\{\begin{split}&SO{C}_{bat.{\mathrm{min}}}\leqslant SO{C}_{bat}\leqslant SO{C}_{bat.{\mathrm{max}}} ,\\ &SO{C}_{sc.\text{min}}\leqslant SO{C}_{sc}\leqslant SO{C}_{sc.\text{max}} 。\end{split}\right.$ | (13) |
式 中:
3)储能容量约束
| $ {{E}_{HESS\_ref}\geqslant \dfrac{\text{max}\left\{\left|\dfrac{\text{max}\left(E\right(t\left)\right)}{{\eta }_\text{dis}{\eta }_\text{inv}}\right|,\left|\text{min}\left(E\right(t){\eta }_{ch}{\eta }_\text{inv}\right|\right\}}{SO{C}_\text{max}-SO{C}_\text{min}}}。$ | (14) |
式中:
4)额定功率约束
额定功率约为:
| $ { {P}_{HESS\_ref}\geqslant \dfrac{\text{max}\left\{\left|\dfrac{\text{max}\left(P\right(t\left)\right)}{{\eta }_\text{dis}{\eta }_{inv}}\right|,\left|\text{min}\left(P\right(t){\eta }_{ch}{\eta }_\text{inv}\right|\right\}}{SO{C}_\text{max}-SO{C}_\text{min}}}。$ | (15) |
式中:
麻雀搜索算法是一种模拟麻雀群体觅食与反捕食行为的群体智能优化方法[7]。算法中,种群个体被划分为“发现者”、“追随者”和“警戒者”三类角色:发现者负责在解空间中进行广域搜索,并在感知到潜在威胁时迅速引导整个种群朝向相对安全的区域迁移;跟随者则紧随发现者,动态更新自身位置以分享食源信息,当能量水平下降至阈值以下时可自主切换位置以获取更多“食物”;警戒者则占总群体的 15%~25%,其主要功能是在外围监测风险信号,向群内发出警报,并根据威胁程度调整自身位置或促使附近个体改变搜索方向。
在具体应用中,SSA 将每只“麻雀”的位置向量视为问题的一个候选解,其对应的“能量储备”与适应度函数值相映射。随着迭代进行,算法按照一定的更新公式不断调整发现者与跟随者的位置,并借助警戒者的风险预警机制在局部陷阱处跳出,从而逐步逼近全局最优解[8]。
发现者的位置更新式为:
| $ {X}_{i,j}^{t+1}=\left\{\begin{aligned}&{X}_{i,j}^{t}\cdot{\mathrm{ {exp}}}\left(\frac{-i}{\alpha \cdot i{t}_\text{max}}\right),{R}_{2} \lt ST,\\& {X}_{i,j}^{t}+Q\cdot L,{R}_{2}\geqslant ST。\end{aligned}\right. $ | (16) |
式中:
追随者的位置更新式为:
| $ {X}_{i,j}^{t+1}=\left\{\begin{aligned}&Q\cdot {\mathrm{{exp}}}\left(\frac{{X}_{w}^{t}-{X}_{i,j}^{t}}{{i}^{2}}\right),i \gt n/2,\\ &{X}_{p}^{t+1}+\left|{X}_{i,j}^{t}-{X}_{p}^{t+1}\right|\cdot {A}^+\cdot L,\text{其他}。\end{aligned}\right. $ | (17) |
式中:
而在模拟实验中,警戒者通过随机方式初始化其初始位置:
| $ {X}_{i,j}^{t+1}=\left\{\begin{aligned}&X_{i,j}^{t+1}+\beta \cdot \left|{X}_{i,j}^{t}-{X}_{b}^{t}\right|,{f}_{i} \gt {f}_{g},\\& {X}_{i,j}^{t+1}+m\cdot \left(\frac{\left|{X}_{i,j}^{t}-{X}_{w}^{t}\right|}{({f}_{i}-{f}_{w})+\delta }\right),{f}_{i}={f}_{g}。\end{aligned}\right. $ | (18) |
式中:
若麻雀检测到捕食者(适应度变化剧烈),进行逃逸,更新公式为:
| $ {X}_{i}^{t+1}={X}_{i}^{t}+\beta \cdot {\mathrm{sign}}\left({X}_{i}^{t}-{X}_{w}^{t}\right)。$ | (19) |
式中:因子
然而SSA算法的不足主要在于种群初始化随机性导致多样性不足、迭代后期易陷入局部最优以及位置更新参数随机性较大,影响搜索效果。针对不足,本文对麻雀搜索算法做出以下改进:
1)在种群初始化及后续迭代阶段引入折射反向学习机制,折射率
| $ n=\frac{\text{sin}\alpha }{\text{sin}\beta }\cdot \frac{{D}^{*}\left(\right(a+b)/2)-x}{D({x}^{*}-(a+b)/2}。$ | (20) |
令
| $ x'=\frac{a+b}{2}+\frac{{D}^{*}\left(\left(a+b\right)/2\right)-x}{2nj} 。$ | (21) |
式中:
2)引入差分变异、交叉与选择机制。生成个体的变异解,采用自适应交叉操作生成新的候选解,再利用贪婪选择机制保留更优的个体,实现算法性能的有效提升[9]。
3)引入动态步长因子。
| $ \beta ={f}_{b}-\left({\frac{i{t}_\text{max}}{i{t}_\text{max}}}^{1.5}\right)\cdot \left({f}_{b}-{f}_{w}\right),$ | (22) |
| $ n={\mathrm{{exp}}}\left(-20\cdot \text{tan}\left(\frac{t}{i{t}_\text{max}}\right)\right)\left(2{\mathrm{rand}}-1\right)\left({f}_{b}-{f}_{w}\right)。$ | (23) |
式中:
基于上述改进策略,提出一种改进型多目标麻雀搜索算法(MSSA)用于复合储能系统容量优化配置,整个优化流程如图4所示。
|
图 4 改进后的麻雀优化算法流程 Fig. 4 Flowchart of the improved sparrow search algorithm |
考虑到超级电容与直流母线直接相连,当其电压过低时可能会造成母线电压波动或电流激增,因而将其SOC初值设置为85%,正常工作范围设置为0.65~0.85。本文设置的仿真实验初始条件如表1和表2所示。
|
|
表 1 优化算法参数设置 Tab.1 Parameter settings for the optimization algorithm |
|
|
表 2 复合储能系统参数设置 Tab.2 Parameter settings for the composite energy storage system |
以设计的基于等效最小氢耗目标的模糊控制功率分配策略为基础,融合限值保护机制后获得的锂电池与超级电容输出功率如图5所示,这些数据作为优化模型的输入条件。
|
图 5 优化问题输入数据 Fig. 5 Input data for the optimization problem |
在Pareto前沿上,所有最优点对应的优化目标中,2个优化目标存在矛盾,无法同时满足,此处采用TOPSIS法,即最优方案应当距离在所有评价指标中最优的值所构成的向量最近,距离在所有评价指标中最差的值所构成的向量最远。本研究在考虑这2个目标的基础上,制定了5种不同的权重分配方案,如表3所示。
|
|
表 3 不同权重分配方案 Tab.3 Different weight allocation schemes |
在给定的目标权重条件下,能耗指标与质量特性指标对应的优化目标值如表4所示。
|
|
表 4 不同权重下子函数的目标值 Tab.4 Objective values of sub-functions under different weightings |
为兼顾能耗与系统寿命两方面性能,本文选取能量消耗与储能寿命质量目标函数权重均为0.5的均衡方案。最终配置方案如表5所示。
|
|
表 5 目标权重下的最优解 Tab.5 Optimal solutions under different objective weightings |
表6对比了优化前后复合储能系统的配置与性能指标。分析可知,复合储能系统的能量消耗降低了1.68%,质量特性降低了5.99%,系统的经济性优化了16.5%,对复合储能系统整体的经济性和动力性做出了优化。
|
|
表 6 优化前后储能系统容量对比 Tab.6 Capacity comparison of the energy storage system before and after optimization |
1)优化前后锂电池SOC对比
图6为优化前后锂电池 SOC 的实时变化曲线,优化后SOC曲线在整个运行过程中更加平缓。在 0~35 s、95~115 s 及 165~175 s 的进出港阶段,2种配置下的 SOC 波动幅度相当,但在经济巡航阶段,优化后下降速率明显减缓。容量优化不仅降低了最大放电深度,还能有效延长电池的循环寿命。
|
图 6 容量优化前后锂电池SOC Fig. 6 SOC of the lithium-ion battery before and after capacity optimization |
2)优化前后母线电压对比
图7为优化前后直流母线电压的实时曲线,在完整运行周期内,优化后的复合储能系统依然具备优秀的充放电性能,结合能量管理策略,可有效维持母线电压稳定并保障电能质量。
|
图 7 容量优化前后直流母线电压 Fig. 7 DC bus voltage before and after capacity optimization |
3)优化前后储能元件输出功率对比
图8和图9分别给出了优化前后锂电池与超级电容的实时输出功率曲线。可以看出,优化后的储能容量方案在经济巡航和高速巡航阶段与优化前几乎相同;在进出港阶段,锂电池略微承担更多功率,而超级电容为维持母线电压稳定,输出的高频功率略有减少,两者合计功率轻微下降。由此可见,适当削减超级电容容量并提升锂电池容量后,不仅保持了电能质量与系统响应性能,也实现了系统质量特性的显著优化。
|
图 8 容量优化前后锂电池输出功率 Fig. 8 Output power of the lithium-ion battery before and after capacity optimization |
|
图 9 容量优化前后超级电容输出功率 Fig. 9 Output power of the supercapacitor before and after capacity optimization |
4)等效氢气消耗对比
如图10所示,优化后典型工况下的等效氢耗约21.35 g,仅较优化前增加0.47%,几乎持平。其原因在于锂电池容量提升后略多承担功率,而系统整体功率消耗有所下降,这部分功率差由燃料电池补充,导致等效氢耗微增。
|
图 10 容量优化前后等效氢气消耗量 Fig. 10 Equivalent hydrogen consumption before and after capacity optimization |
综上所述,所提出的容量优化配置方案具有良好的有效性与准确性。在船舶动力性基本不受影响的情况下,降低复合储能系统的能量消耗,同时提升船舶的动力性。
5 结 语本文针对复合储能系统的容量配置问题进行了优化设计。选取锂电池和超级电容的串并联数量作为优化决策变量,并明确其可行范围;以复合储能系统的能量消耗和质量特性作为多目标优化函数,同时设置瞬时功率平衡、储能荷电状态(SOC)、最小储能容量以及额定功率为约束条件。为实现优化求解,采用改进的麻雀搜索算法(SSA),并结合TOPSIS方法进行多属性决策,最终确定最优容量配置参数。通过仿真验证表明,优化后的配置方案不仅有效保障了复合储能系统的动力性能,还显著降低了系统的质量特性和能量损耗。
| [1] |
HATATA A Y, OSMAN G, ALADL M M. An optimization method for sizing a solar/wind/battery hybrid power system based on the artificial immune system[J]. Sustainable Energy Technologies & Assessments, 2018, 27: 83−93.
|
| [2] |
XIAO Y, JIA L, PEI L, et al. Optimal scheduling of microgrid with distributed power based on water cycle algorithm[J]. Energies, 2018, 11(9): 2381. DOI:10.3390/en11092381 |
| [3] |
ELDEEB H H, ELSAYED A T, LASHWAY C R, et al. Hybrid energy storage sizing and power splitting optimization for plug-In electric vehicles[J]. IEEE transactions on industry applications, 2019, 55(3): 2252-2262. DOI:10.1109/TIA.2019.2898839 |
| [4] |
冯龙祥, 汤旭晶, 李昕, 等. 计及电池全寿命周期的船舶复合储能多目标容量配置研究[J]. 中国航海, 2023, 46(1): 30-38. FENG L X, TANG X J, LI X, et al. Configuring ship duplex energy storage when capacity life cycle of cell is taken into account[J]. Navigation of China, 2023, 46(1): 30-38. DOI:10.3969/j.issn.1000-4653.2023.01.005 |
| [5] |
宋传学, 周放, 肖峰, 等. 基于凸优化的车载复合电源参数匹配[J]. 机械工程学报, 2017, 53(16): 44-51. SONG CH X, ZHOU F, XIAO F, et al. Parameter matching of on-board hybrid energy storage system based on convex optimization method[J]. Journal of Mechanical Engineering, 2017, 53(16): 44-51. DOI:10.3901/JME.2017.16.044 |
| [6] |
肖能齐, 徐翔, 周瑞平. 船舶柴电混合动力系统能量管理控制策略[J]. 哈尔滨工程大学学报, 2020, 41(1): 153-160. XIAO N Q, XU X, ZHOU R P. Energy management and control strategy of ship diesel-electric hybrid power system[J]. Journal of Harbin Engineering University, 2020, 41(1): 153-160. DOI:10.11990/jheu.201808007 |
| [7] |
马纪梅, 张欣彤, 张政林, 等. 基于改进麻雀搜索算法的微网容量优化配置[J]. 电子测量技术, 2022, 45(8): 76-82. MA J M, ZHANG X T, ZHANG ZH L, et al. Microgrid capacity optimization based on improved sparrow search algorithm[J]. Electronic Measurement Technology, 2022, 45(8): 76-82. DOI:10.19912/j.0254-0096.tynxb.2024-0968 |
| [8] |
曾律淳, 尚前明, 孙钰波. 基于灰狼算法的混合动力船舶容量优化配置[J]. 舰船科学技术, 2025, 47(7): 116-124. ZENG L CH, SHANG Q M, SUN Y B. Hybrid propulsion ship capacity configuration optimization based on grey wolf algorithm[J]. Ship Science and Technology, 2025, 47(7): 116-124. DOI:10.3404/j.issn.1672-7649.2025.07.022 |
| [9] |
陈辉, 李维波, 孙万峰, 等. 基于优化功率分配的光伏混合储能系统能量管理策略[J]. 中国舰船研究, 2021, 16(2): 168-175. CHEN H, LI W B, SUN W F, et al. Energy management strategy of photovoltaic hybrid energy storage system based on optimal power distribution[J]. Chinese Journal of Ship Research, 2021, 16(2): 168-175. DOI:10.19693/j.issn.1673-3185.01869 |
2026, Vol. 48
