固体氧化物燃料电池与微型燃气轮机联合发电建模仿真研究

朱润凯，梁前超，闫 东，詹海洋

(海军工程大学 动力工程学院，湖北 武汉 430033)

固体氧化物燃料电池与微型燃气轮机联合发电建模仿真研究

朱润凯，梁前超，闫 东，詹海洋

(海军工程大学 动力工程学院，湖北 武汉 430033)

设计一种固体氧化物燃料电池与微型燃气轮机的联合发电系统，并对其进行建模仿真及性能研究。基于 Matlab/Simulink 仿真软件，采用模块化建模方法，再以拓扑结构连接各个子系统模型，搭建 SOFC-MGT 联合发电系统模型。仿真结果表明：模型满足仿真要求，系统具有较高的发电效率。在变工况，尤其是高背压条件下，该系统不能稳定工作，必须附加水处理系统才能保证该系统良好运行。

固体氧化物燃料电池；微型燃气轮机；高背压；联合发电；建模仿真

0 引 言

燃料电池是将储存在氧化剂与燃料内部的化学能直接转化为电能的发电装置，其中固体氧化物燃料电池（Solid oxide fuel cell，SOFC）是一种中高温燃料电池，具有诸多一般燃料电池所不备的优越性[1]，实验条件下已取得 60% 以上的发电效率，排气温度达 600 ℃以上[2-3]。将高品位的废热与燃气轮机组成联合动力系统，可以进一步提高设备发电效率。

计算机建模和仿真模拟技术，能够节约实验研究的成本，缩短研究周期，寻找规律，发现明显设计缺陷。本文重点考虑在高背压（1.7 kg/cm2）各工况工作时，对联合发电系统进行动态建模和性能分析，并使其稳定高效运行，为后期的实际应用开发提供理论指导。根据文献[4] 采用容阻特性建模，完成 SOFC 一维模型的快速动态仿真，在高背压 1.7 kg/cm2情况下，并有效结合微型燃气轮机（Micro gas turbine，MGT），对联合装置进行了性能仿真。本文主要抓住研究对象的稳态特性，采用 0 D 建模方法，得到仿真条件下装置的运行特性。

1 仿真数学模型

1.1 系统结构描述

本文综合考虑了建模对象的工作特点[5]，重点应用在高背压（1.7 kg/cm2）环境下，研究 SOFC-MGT 发电系统的运行工况及变工况性能。按照燃料气流经SOFC 内部各模块的顺序进行建模，并对部分复杂过程进行了简化，做出如下假设：

1）所有气体均为理想气体；

2）燃料选择为 100% 甲烷 CH4，流量为 0.007 2 kg/s；

3）系统与外界无传热传质过程；

4）重整反应及水气置换反应均处于平衡状态；

5）采用集总参数模型，独立模块内部各状态参数保持一致。

SOFC 模型，采用文献[6]中的外部预重整循环结构，MGT 主要用于对燃料电池进口气体加压。由于系统运行直接为工作工况，未考虑启动条件，因而省去了外加预热设备。图 1 为联合发电装置的拓扑结构[7]。

1.2 计算模型

采用模块化建模思路，按照图 1 所描述的拓扑结构构建计算模型。建模过程主要依据质量守恒和能量守恒定律将各模块内部反应前后变化表达出来[8]。

1.2.1 SOFC 模型

SOFC 工作原理同所有燃料电池一样，都是通过电解质层中的载流子传递电荷，实现 CH4和 O2电化学反应的过程。燃料首先经过压气机加压进入混合器，与阳极出口的部分尾气混合，形成含有 CH4和 H2O 的混合气体。随后进入预重整器进行蒸汽重整反应和水汽置换反应，生成 H2，重整后的燃料气进入 SOFC 阳极，此过程只能够实现 CH4的部分重整。另一方面，经压气机加压后的空气，在换热器中预热后进入SOFC 阴极。在电池内部高温环境下，CH4进一步重整生成 H2，同时，在电解质内部发生电化学反应，消耗两极送来的 H2和 O2，产生的电子经过外电路流出系统，即实现发电功能。

1）预重整器

经过混合器充分混合的气体中，主要成分是CH4、H2O 和 CO2，接着进入预重整器内，发生重整反应与水气置换反应，进一步增加氢气的含量。

依据质量守恒定律，反应过后，各组分摩尔分数变化情况为：

式中：i为各气体组分；P8,i为气体各组分在预重整器出口的分压，Pa；VrsVrs为预重整器体积，m3；T8为出口总温，K；R为气体普适常数，8.317 J/（mol·K）；G7,i，G8,i为各组分气体在预重整器进出口的流量，mol/s；X8,i为各组分气体在出口的摩尔分数，即X8,i=G8,i/G8；Rrs,i为各组分的反应速度，mol/s。

2）电堆模型

假设电堆内部各参数保持一致性，根据质量守恒定律，可以求出电堆内部的物质变化，以氧气 O2在阴极板中的反应为例。

式中各量与式（1）类似。

3）电化学模型

燃料电池实际运行中的电压Vfc为理想的开路可逆电压E减去各种极化损失引起的电压损失，可表示为[9]：

式中：VOM为欧姆极化；VCONC为浓度差极化，VACT,a、VACT,ca分别为阴极阳极的活化极化，V[10]。

根据 Nernst 方程可得到理想可逆电压[11]：

式中：E0为标准状况下燃料电池的理想电势，与温度TS有关；F为法拉第常数，96 485.338 C/mol；pref为标准压力，Pa。

得到单电池的工作电压后，便可求出 SOFC 电堆的输出电压VST和输出功率PST：

式中：m为电堆的单电池个数；i为电堆内部的电流密度，A。

4）温度模型

式中：k= 1，2，3；i分别代表 CH4、CO、CO2、H2、H2O，j分别代表 O2、N2；Cs为电堆气体的整体比热容，kJ/（kg·K）；h为焓，kJ/mol[12]；Qk分别为重整反应、水气置换反应、电化学反应的反应热，kJ/mol。

5）催化燃烧室模块

采用的燃烧室为一个气相反应室，尾气在里面不发生燃烧反应，而是通过催化剂来促进氧化反应进行，进一步将未充分利用的燃料气消耗，提升尾气温度，从而强化尾气做功能力。

根据质量守恒定律可得：

式中：Vb为燃烧室体积，m3；Rb,i为燃烧过程中消耗的气体流量，mol/s。

对于燃烧室出口气体的温度，假设与燃烧室内部温度一致，考虑到燃烧室的燃烧效率，根据能量守恒方程知：

式中，i分别代表 CH4，CO，CO2、H2，H2O；j分别代表 O2，N2；Cb为燃烧室内气体的比热容，kJ/（kg·K）；Qt分别为 H2，CO 的反应热，kJ/mol。

1.2.2 燃气轮机模型

SOFC 的阴极和阳极反应后的气体中仍含有部分O2和可燃气体成分（CO，H2），尾气送入催化燃烧室内，进一步将燃料燃烧利用完全。从催化燃烧室出来的尾气温度很高，在换热器中将压缩空气预热后，推动涡轮做功，带动发电机发电。

燃气轮机以已有 MGT 为依据，经过适当的数学变换得到压气机以及涡轮的数学模型[13]。

1）压气机模型

压气机是一个具有强非线性的部件，这里采用 C30微型燃机轮机作为建模基准[14]。其工作特性可由压比π、折合流量折合转速以及效率 ηC来表示[7]。

式中：qm为压气机的空气质量流量，kg/s；p为进气总压，Pa；T为进气总温，K；n为转速，r/min。

式中：T1，T2为压气机进出口温度，K；k为绝热系数；为进出口焓值，kJ/kg；PC为压气机耗功，kW。

2）涡轮模块

涡轮模型同样可以由相似理论进行特性分析。气体膨胀做功推动涡轮转动，压力下降，膨胀比为 ε。

3）转子模块

在仿真过程中转子转速为状态量，其工作特性可以由转子能量平衡得到[16]：

式中：r为转子的转动惯量，kg/m3；PMGT为转子不平衡功率，kW；ηm为涡轮做功效率；NT，Ng，Nf分别为涡轮总功率、MGT 输出功率和辅机负载功率，kW。

2 性能仿真分析

实现高背压下系统变工况运行，可以通过控制燃料流量、电流、转子转速、涡轮前温度等参数实现，并以联合发电系统输出总功率为工况选择参考点[16-17]。由稳态仿真计算结果可看出，燃料电池输出功率占系统总功率 80% 以上，因而本文将以满足燃料电池高效率工作为前提，采用燃气轮机转速恒定下，控制燃料流量实现系统输出功率变化的仿真策略。在保证系统正常运行基础上，燃气轮机将配合燃料电池工作，其工作效率将有所降低。

本文计算了 0.7～1.0 工况的系统运行情况。图 2 为SOFC-MGT 变工况仿真中，SOFC 电堆工作温度的变化曲线。可以看出，工况由 1.0 下降到 0.7，电堆内部反应强度减弱，工作温度下降约 30 K。燃机转速保持恒定，系统内部空气流量将基本不变，燃料流量随功率下降使燃空比降低，影响到电堆内部反应强度。

图3 显示了 SOFC-MGT 变工况运行中，SOFC 内部反应气体成分的变化特点。可以看出工况下降的过程中，H2成分显著下降，H2O，CO2明显增加，CH4和 CO 含量略有下降，基本保持不变。这是因为低工况蒸汽重整反应平衡点左移，水汽置换反应强度也随之降低。出口气体中可燃气体成分甲烷含量变多，进而使燃料利用率降低。

图4 为 SOFC 与 MGT 各系统的功率特性曲线。在转速恒定，空气流量基本不变的情况下，压气机压缩功也将基本保持不变，而涡轮入口气体温度下降，将使涡轮的膨胀功减少，因而 MGT 输出功率将会减小。另一方面，图 4 还显示在 SOFC-MGT 联合发电系统中，燃料电池提供主要的功率输出，所以整个系统功率特性与 SOFC 功率特性变化趋势接近。

图5 显示了变工况运行中，联合发电系统效率特性变化曲线。可以看出系统的工作效率会随着总功率下降而下降。在 0.72 工况处，系统效率下降到48.5%，非设计点系统工作性能比设计点差距较大，主要是由于低工况燃气轮机工作效率的迅速下滑。

3 结 语

在 Matlab/Simulink 仿真环境下，建立了 SOFCMGT 联合发电系统的仿真模型，2个子系统能够匹配运行。结果显示，在高背压（1.7 kg/cm2）条件下，其设计点工况输出功率为 223.6 kW，发电效率较低，仅为 57.9%，并且变工况效率急剧恶化。比较文献[4] 和[7] 常温常压下的仿真结果，表明：该背压系统必须加载专门的水处理系统，降低背压，才能满足 SOFC-MGT联合发电系统的高效运行。

[ 1 ]冯登满, 刘畅, 杨苗苗, 等. 固体氧化物燃料电池的发展与研究[J]. 科技展望, 2016(5): 60-60.

[ 2 ]耿孝儒, 吕小静, 翁一武. 基于生物质气的固体氧化物燃料电池-燃气轮机混合动力系统的性能分析[J]. 动力工程学报, 2015, 35(2): 166-172. GENG Xiao-ru, LV Xiao-jing, WENG Yi-wu. Performance analysis of an SOFC-GT hybrid system fueled with biomass[J]. Journal of Chinese Society of Power Engineering, 2015, 35(2): 166-172.

[ 3 ]李杨, 翁一武. 固体氧化物燃料电池-燃气轮机混合动力系统及控制策略分析[J]. 中国电机工程学报, 2010, 30(35): 94-100. LI Yang, WENG Yi-wu. Performance study and control strategies of temperature solid oxide fuel cell-gas turbine hybrid system[J]. Proceesings of the CSEE, 2010, 30(35): 94-100.

[ 4 ]黄刘松. 固体氧化物燃料电池-燃气轮机混合发电装置建模与仿真[D]. 上海: 上海交通大学, 2007.

[ 5 ]TIAN Yun. Modeling of solid oxide fuel cell-gas turbine hybrid power plant[D]. Singapore: Nanyang Technological University, 2002.

[ 6 ]GEORGE R A. Status of tubular SOFC field unit demonstrations[J]. Journal of Power Sources, 2000, 86(1-2): 134-139.

[ 7 ]吴小娟. 固体氧化物燃料电池/微型燃气轮机混合发电系统的建模与控制[D]. 上海: 上海交通大学, 2009.

[ 8 ]胡傲生, 吴锦超, 万昕. 固体氧化物燃料电池阳极的微观建模[J]. 电源技术, 2016, 40(2): 302-323. HU Ao-sheng, WU Jin-chao, WAN Xin. Anode micro modeling of solid oxide fuel cell[J]. Chinese Journal of Power Sources, 2016, 40(2): 302-323.

[ 9 ]徐泽亚, 杨晨. 管式 SOFC 数学模型及系统性能分析[J]. 计算机仿真, 2006, 23(6): 224-228. XU Ze-ya, YANG Chen. Mathematical model of tubular solid oxide fuel cell and system performance analysis[J]. Computer Simulation, 2006, 23(6): 224-228.

[10]AKKAYA A V. Electrochemical model for Performance analysis of a tubular SOFC[J]. International Journal of Energy Research, 2006(1), 31: 79-98.

[11]毛晓峰, 由宏新, 阿布里提·阿布都拉, 等. 电解质支撑固体氧化物燃料电池发电性能模拟[J]. 天然气工业, 2004, 24(6): 107-110. MAO Xiao-feng, YOU Hong-xin, Abuliti·Abudula, et al. Electrical performance simulation of electrolyte-supported solid oxide fuel cells[J]. Natural Gas Industry, 2004, 24(6): 107-110.

[12]刘彬, 卢荣, 物理化学[M]. 武汉: 华中科技大学出版社, 2008.

[13]刘尚明, 李忠义. 基于 SIMULINK 的单轴重型燃气轮机建模与仿真研究[J]. 燃气轮机技术, 2009, 22(3): 33-45. LIU Shang-ming, LI Zhong-yi. Modeling and simulation of single shafe heavy-duty gas turbine using SIMULINK[J]. Gas Turbine Technology, 2009, 22(3): 33-45.

[14]马丽娜. 微型燃气轮机 C30 气体排放实验研究[D]. 徐州: 中国矿业大学, 2014. MA Li-na. Experimental study on the exhaust of micro gas turbine C30[D]. Xuzhou: China University of Mining and Technology, 2014.

[15]THORUD. B. Dynamic modelling and characterisation of a solid oxide fuel cell integrated in a gas turbine cycle[D]. Trondheim: Norwegian University of Science and Technology, 2005.

[16]CALISEA F, ACCADIA M D D’, PALOMBO A, et al. Simulation and exergy analysis of a hybrid solid oxide fuel cell(SOFC)-gas turbine system[J]. Energy, 2006, 31(15): 3278-3299.

[17]张健. 固体氧化物燃料电池与燃气轮机联合循环特性分析[D]. 保定: 华北电力大学, 2007.

Modeling and simulation for solid oxide fuel cell and micro gas turbine combined power generation system

ZHU Run-kai, LIANG Qian-chao, YAN Dong, ZHAN Hai-yang
(College of Naval Architecture and Power, Naval University of Engineering, Wuhan 430033, China)

A solid oxide fuel cell and micro gas turbine combined power generation system was designed and considered the working performance through modeling and simulation. The model is established by modularization modelling method based on Matlab/Simulink. Then all submodulars are connected by topological structure and make up the SOFCMGT combined power generation system. After simulation, the result shows that the model can meet the requirements of the simulation precision and the system can operate with high efficiency. At off-design condition, especially when the high back pressure exists, the system can’t be stably operating. In order to ensure it running well, water-processing must be loaded.

solid oxide fuel cell；micro gas turbine；high back pressure；combined power generation；modeling and simulation

TM911；TK472

A

1672 - 7619(2017)04 - 0095 - 05

10.3404/j.issn.1672 - 7619.2017.04.019

2016 - 06 - 13；

2016 - 07 - 08

“十二五”国防预研基金资助项目(40103030403)；湖北省协同创新基金资助项目(HX2015B1003)

朱润凯(1993 - )，男，硕士研究生，研究方向为动力机械及热力系统的设计、仿真与优化。