-
为了应对全球气候变化，满足能源与环境的可持续发展，煤电逐渐被风能、太阳能等可再生能源发电技术取代[1]。我国已提出了“碳达峰、碳中和”的发展战略，要求2060年非化石能源比例上升到80%，实现以新能源为主体、多种形式能源共同构成的多元化能源系统[2]。由于风能和太阳能的不稳定性和间歇性，大规模可再生能源纳入电网会影响现有电网的稳定性，增加了电网的调控压力。面对这些问题，储能技术应运而生，解决了可再生能源电站输出功率不稳定、连续性差的问题，提高了电力系统的效率和安全性[3-5]。
储能技术按照运行原理可分为机械储能、电磁储能和化学储能，压缩空气储能属于机械储能的一种[6]。早期的压缩空气储能系统（Compressed Air Energy Storage System）大多是补燃式，通过与燃气轮机相配合，展现出了较高的可靠性、稳定性和灵活性，但是化石燃料燃烧会产生温室气体，对环境造成污染。为此，更多的CAES（Compressed-Air Energy Storage, CAES）技术，如先进绝热CAES（Advanced Adiabatic-CAES，AA-CAES）、超临界CAES （Supercritical-CAES，SCCAES）、水下CAES（Under Water-CAES，UWCAES）、等温CAES（Isothermal CAES，I-CAES）等被提出。其中，AA-CAES通过回收利用压缩过程中产生的热量，不仅实现了储能、释能全过程无污染，还提高了系统的能效[7–9]。
目前，AA-CAES成为全球广泛关注和研究的热点之一，已建立多座AA-CAES示范系统。在我国，国家电网公司、中国科学院工程热物理研究所、清华大学等也取得了一系列重要成果，如兆瓦级超临界CAES综合实验平台、500 kW AA-CAES实验系统等，为我国AA-CAES技术的发展积累了宝贵的经验[10]。国内外学者对AA-CAES进行了大量的理论和实验研究，Guo等[11]构建了AA-CAES热力学模型，并分析了影响系统效率和储能密度的关键因素。张新敬[12]研究了AA-CAES系统在运行时，压比、温度等参数变化对系统性能的影响。Chen等[13]分析了等压AA-CAES系统的热力学特性，探究了关键参数变化对系统㶲效率的影响。韩中合等[14]分析了系统储能效率和储能密度随着压缩级数、压缩结构以及不同储热介质变化的规律。刘畅等[15]构建了压缩空气储能系统经济性模型，并对某蓄热式压缩空气储能电站进行了盈亏平衡分析以及敏感性分析。
相较于大型压缩空气储能系统需要地下岩层、废弃矿井等较大的储气室，微小型压缩空气储能系统常采用储气罐来存储压缩空气，从而摆脱了对特定地理环境的依赖，使用区域可以更广泛[16]，并且微小型压缩空气储能系统规模较小，建设和运营成本相对较低，更适合于分布式能源系统和微电网的应用。基于上述优势，微小型压缩空气储能系统有更广阔的应用前景，比如在一些偏远地区、岛屿或临时设施中，电力供应可能不稳定或难以通过传统方式获得[17]，微小型压缩空气储能系统就可以作为一种可靠的电力保障方案，为这些地区提供持续、稳定的电力供应。此外，它还可以用于应急电源、备用电源等场景，确保关键设施和设备的连续运行[18-19]。
近年来，诸多国内外学者对微小型压缩空气储能系统的运行特性进行了研究。王成山等[20]建立了微型压缩空气储能系统的热力学模型，研究了压缩空气的压力、温度等因素对系统效率的影响。Khamis等[21]搭建了一台微型压缩空气储能试验系统，研究了透平入口压力对透平出口温度、转速、输出轴功和系统效率的影响。Martínez等[22]建立了微型压缩空气储能系统的动态数学模型，并利用Matlab软件进行了仿真分析。褚晓广等[23]建立了微型压缩空气储能系统动态数学模型，并通过仿真与实验的方法对涡旋压缩机效率优化控制进行了研究。Lemofouet等[24]构建了一个由微型压缩空气储能与超级电容组成的混合储能系统，提出了一种系统最大效率跟踪的控制方法，并通过实验验证了该方法的可行性。
鉴于微小型压缩空气储能系统在电力系统中的重要作用以及对其变工况特性研究的不足，开展系统的变工况特性研究具有重要的理论价值和实践意义。文章建立了小型先进绝热压缩空气储能系统的静态模型和动态模型，以10 kW级储能系统为对象，分析不同工况对系统性能和储气罐内空气动态变化特性的影响。
-
AA-CAES系统与传统压缩空气储能系统相比，取消了燃烧室，引入了冷/热储能罐，用于储存压缩过程产生的热能。如图1所示，该系统主要由压缩机、膨胀机、换热器、储气罐、节流阀、储热罐和储冷罐组成。储能时，周边大气经多级压缩升至高温高压状态，来自储冷罐中的换热介质依次通过间冷器吸收各级压缩热，并集热于储热罐，压缩空气降温进入储气罐；释能时，储气罐中的高压空气从储热罐中依次吸热、经多级膨胀做功，带动发电机发电，来自储热罐中的换热介质被冷却后进入储冷罐，末级膨胀结束排废气于大气中，系统完成一次工作或循环。
Figure 1. AA-CAES system flowchart
-
为了分析不同工况对AA-CAES系统效率的影响，需对系统中关键设备建立静态模型和动态模型，数值分析其热力学过程，获得压缩机、膨胀机、换热器、储气罐中空气及介质参数变化对系统性能的影响趋势。为简化模型和分析，做以下假设：
1）把空气视为理想气体，其性质满足理想气体状态方程。
2）忽略所有连接部件管道中压力和热量的损失。
3）不计储热罐中的热量损失。
4）压缩与膨胀过程采用等压比和等膨胀比方式。
-
压缩机是压缩空气储能系统的重要部件之一，对于微型压缩空气储能系统，压缩机常采用容积式结构。认为压缩过程为绝热过程，但级间进行气体冷却。压缩空气进、出口比焓与等熵效率的关系见式（1）。
$$ {h_{{\text{c,out}}}} = {h_{{\text{c,in}}}} + \frac{{{h_{{\text{c,out},i}}} - {h_{{\text{c,in}}}}}}{{{\eta _{\text{c}}}}} $$ (1) 式中：
$ {h}_{{\mathrm{c,in}}} $、$ {h}_{{\mathrm{c,out}}} $ −实际等熵压缩过程中压缩机进、出口处的比焓（J/kg）；
$ {h}_{{\mathrm{c,out},i}} $ −理想等熵过程中压缩机出口的比焓（J/kg）；
ηc −压缩机的等熵效率。
第i级压缩机的功率可表达为：
$$ {P_{{\text{c},i}}} = ({h_{{\text{c,out}}}} - {h_{{\text{c,in}}}}){{{q}}_{\text{m}}} $$ (2) 多级压缩机的总功耗可按下式计算：
$$ {P_{\text{c}}} = \sum\limits_{{i = 1}}^{{N_{\text{c}}}} {{P_{{\text{c},i}}}} $$ (3) 式中：
$ {q}_{{\mathrm{m}}} $ −进入压缩机的空气质量流量（kg/s）。
NC −压缩机的级数。
-
在换热器中，换热器效能是影响换热器性能的关键指标，定义如下：
$$ \varepsilon =\frac{{\left|{t}^{{'}}-{t}^{''}\right|}_{\text{max}}}{\left|{t}_{1}^{{'}}-{t}_{2}^{{'}}\right|} $$ (4) 式中：
$ {\left|{t}^{{'}}-{t}^{''}\right|}_{\mathrm{max}} $−换热器中高温介质或低温介质温差较大的值（K）；
$ \left| {t_1^{'} - t_2^{'}} \right| $ −换热器理想状态下实现的最大温差（K）。
储能时，将第i级压缩机的出口温度视为第i级换热器的进口温度，第i级换热器的出口温度视为第i+1级压缩机的进口温度。所以在压缩阶段，空气在换热器出口温度为：
$$ T_{{\text{c},i + 1}}^{{\text{in}}} = {\varepsilon _{\text{c}}}T_{{\text{ch},i}}^{{\text{in}}} + (1 - {\varepsilon _{\text{c}}})T_{{\text{c},i}}^{{\text{out}}} $$ (5) 式中：
$ {{T}}_{{\mathrm{c},i}+1}^{{\mathrm{in}}} $ −第i级换热器出口的空气温度（K），也是i+1级压缩机的进口温度（K）；
$ {\varepsilon }_{{\mathrm{c}}} $ −压缩过程中换热器的效能；
$ {{T}}_{{\mathrm{ch},i}}^{{\mathrm{in}}} $ −第i级换热器换热介质的进口温度（K）；
$ {{T}}_{{\mathrm{c},i}}^{{\mathrm{out}}} $ −第i级压缩机的出口空气温度（K）。
根据文献[14]，当换热介质的热容率（比热容与质量流量的乘积）与空气热容率相等时，系统储能效率最高，所以设换热介质热容率等于空气热容率，则换热器出口的换热介质温度为：
$$ T_{{\text{ch},i}}^{{\text{out}}} = {\varepsilon _{\text{c}}}T_{{\text{c},i}}^{{\text{out}}} + (1 - {\varepsilon _{\text{c}}})T_{{\text{ch},i}}^{{\text{in}}} $$ (6) 根据能量守恒定律，压缩空气释放的能量等于换热介质吸收的能量，有：
$$ Q = (h_{{\text{air},i}}^{{\text{in}}} - h_{{\text{air},i}}^{{\text{out}}}){q_{\text{m}}} = (h_{{\text{ch},i}}^{{\text{out}}} - h_{{\text{ch},i}}^{{\text{in}}}){q_{{\text{ch},i}}} $$ (7) 忽略换热过程中压力损失和流量损失，若给定压缩机的质量流量为$ {q}_{{\mathrm{m}}} $，则换热器的质量流量有：
$$ {{{q}}_{{\text{ch},i}}} = \frac{{(h_{{\text{air},i}}^{{\text{in}}} - h_{{\text{air},i}}^{{\text{out}}}){q_{\text{m}}}}}{{(h_{{\text{ch},i}}^{{\text{out}}} - h_{{\text{ch},i}}^{{\text{in}}})}} $$ (8) $$ {q_{{\text{ch}}}} = \sum\limits_{{i = 1}}^{{N_{\text{c}}}} {{q_{{\text{ch},i}}}} $$ (9) 式中：
$ {h}_{{\mathrm{ch},i}}^{{\mathrm{in}}} $、$ {h}_{{\mathrm{ch},i}}^{{\mathrm{out}}} $−换热介质在第i级换热器的进、出口比焓（J/kg）；
$ {q}_{{\mathrm{ch},i}} $ −第i级换热器中换热介质的质量流量（kg/s）。
膨胀级间换热关系与压缩过程相类似，依照上式可以求得，这里不再赘述。
-
储气罐与压缩机、膨胀机之间不仅有物质交换，还有能量交换，储气罐内空气的状态参数不断发生变化，属非稳定热力学的开口系统问题。在系统储能和释能的过程中，储气罐内的空气状态满足质量守恒方程和能量守恒方程：
$$ \frac{{{\mathrm{d}}m}}{{{\mathrm{d}}t}} = {q_{\text{m}}} - {q_{\text{e}}} $$ (10) $$ \frac{\mathrm{d}(mu)}{\mathrm{d}t}=q_{\text{m}}h_{\mathrm{in}}-q_{\text{e}}h_{\text{out}}-HA(T-T_{\text{env}}) $$ (11) 式中：
t −时间（s）；
$ {q}_{{\mathrm{m}}} $、$ {q}_{{\mathrm{e}}} $ −储能和释能时的空气质量流量（kg/s）；
u −单位质量空气的内能（J/kg）；
$ {{h}}_{{\mathrm{in}}} $、$ {{h}}_{{\mathrm{out}}} $−储能和释能过程空气的比焓（J/kg）；
H −储气罐与外界的对流换热系数（W/m2/K）；
A −储气罐与外界的换热面积（m2）；
T −储气罐内空气温度（K）；
$ {{T}}_{{\mathrm{env}}} $ −外界环境温度（K）。
结合理想气体状态方程：
$$ \frac{\mathrm{d}p}{p}+\frac{\mathrm{d}V}{V}=\frac{\mathrm{d}T}{T}+\frac{\mathrm{d}m}{m} $$ (12) 代入计算可以得出储气罐中空气的温度和压力随时间变化的微分方程：
$$ \begin{split} \frac{\mathrm{d}T}{\mathrm{d}t}= & \frac{1}{m}(q_{\text{m}}\frac{c_{\text{p}}}{c_{\text{v}}}T_{\text{in}}+\left[(1-\frac{c_{\text{p}}}{c_{\text{v}}})q_{\text{e}}-q_{\text{m}}\right]T)- \\ &\frac{HA}{mc_{\text{v}}}(T-T_{\text{env}}) \end{split} $$ (13) $$ \frac{\mathrm{d}p}{\mathrm{d}t}=\frac{R}{V}\frac{c_{\text{p}}}{c_{\text{v}}}(q_{\text{m}}T_{\text{in}}-q_{\text{e}}T)-\frac{HA}{V}\frac{R}{c_{\text{v}}}(T-T_{\text{env}}) $$ (14) 式中：
$ {c}_{{\mathrm{v}}} $、$ {c}_{{\mathrm{p}}} $−空气的定容比热容和定压比热容（J/kg/K）；
$ {T}_{{\mathrm{in}}} $ −进入储气罐的空气温度（K）；
R −气体常数，8.314 J/mol/K；
V −储气罐的容积（m3）。
-
将空气在节流阀处的节流过程视为绝热过程，节流前后空气的焓值相等：
$$ {h_{{\text{before}}}} = {h_{{\text{after}}}} $$ (15) -
膨胀机是CAES系统中对外做功的部件，本文的每一级膨胀对应一台膨胀机。膨胀机的入口压力、温度以及等熵效率对系统效率有重要影响。设定末级的膨胀机出口压力为1.5 bar，各级膨胀比相同，则膨胀机的出口压力为：
$$ {p_{{\text{e},i + 1}}} = {p_{{\mathrm{e},i}}}\,{\beta _{i}} $$ (16) 利用等熵效率的定义，可以求出膨胀机的功率：
$$ {P_{{\text{e},i}}} = ({h_{{\text{e,in}}}} - {h_{{\text{e,out}}}}){q_{\text{e}}} = ({h_{{\text{e,in}}}} - {h_{{\text{e,out,s}}}}){q_{\text{e}}}{\eta _{\text{e}}} $$ (17) 膨胀机产生的总功率为：
$$ {P_{\text{e}}} = \sum\limits_{i = 1}^{{N_{\text{e}}}} {{P_{{\text{e},i}}}} $$ (18) 式中：
$ {\beta }_{{\mathrm{i}}} $ −第i级膨胀机的膨胀比；
$ {{q}}_{{\mathrm{e}}} $ −进入膨胀机的空气质量流量（kg/s）；
$ {h}_{{\mathrm{e,in}}} $、$ {h}_{{\mathrm{e,out}}} $−实际膨胀过程中膨胀机进出口处的比焓（J/kg）；
$ {h}_{{\mathrm{e,out,s}}} $ −等熵膨胀过程中出口的比焓（J/kg）；
$ {\eta }_{{\mathrm{e}}} $ −膨胀机等熵效率；
Ne −膨胀机级数。
-
储能效率定义为释能过程中膨胀机输出的总功与储能过程中压缩机的总耗功之比:
$$ \eta_{\text{elc}}=\dfrac{\displaystyle\int_0^{t_{\mathrm{e}}}P_{\text{e}}\mathrm{d}t}{\displaystyle\int_0^{t_{\mathrm{c}}}P_{\text{c}}\mathrm{d}t} $$ (19) 储能密度定义为释能过程中膨胀机输出的总功与储气罐容积之比：
$$ {\text{ESD}}=\frac{\displaystyle\int_0^{t_{\text{e}}}P_{\text{e}}\mathrm{d}t}{V} $$ (20) -
为了深入研究参数变化对系统效率的影响，本文以10 kW功率等级的AA-CAES系统为对象，利用Mtalab软件构建压缩空气储能仿真系统。如图1所示，本系统压缩机级数为4级，膨胀机级数为3级；采用固定容积的储气罐储气，换热器中的换热介质是1 MPa的高压热水（为防止换热过程中热水汽化，所以选用高压热水），系统的设计工况参数如表1所示。
部件 参数 数值 环境 温度/K 298 压力/MPa 0.1 压缩机 功率/kW 10 等熵效率/% 85 空气质量流量/（kg·s−1） 0.02 储气罐 体积/m3 13 进口温度（K） 318 最大压力/MPa 5 最小压力/MPa 3 储能时间/h 4 释能时间/h 2 换热器 效能 0.85 膨胀机 功率/kW 10 等熵效率/% 85 空气质量流量/（kg·s−1） 0.04 储冷罐 温度/K 298 储热罐 温度/K 411 Table 1. Design working condition parameters of AA-CAES
-
构建AA-CAES系统静态模型操作简单，且能很直观的反应不同因素对系统效率的影响，因此常用于分析系统变工况运行特性以及系统效率优化问题。储气罐采用等温模型，在储能和释能过程中储气罐中空气与环境换热充分，储气罐内空气温度等于环境温度。为了探究不同工况对系统性能的影响，以压缩机进口温度、总压比、膨胀段入口压力、膨胀级数为变量，分析其对系统储能效率和储能密度的影响。
-
压缩机入口处空气的温度即为环境温度，随着当地环境的变化，压缩机入口处的温度也会变化。在分析入口温度的影响时，不考虑压缩机等熵效率的变化，改变入口空气温度，分析系统储能效率和储能密度的变化。
图2为各级压缩机和换热器出口空气温度随压缩机入口空气温度的变化曲线，可以看出：当入口空气温度由273 K到增加到308 K，对第一、二级压缩机的出口温度影响较大，分别增加了47.9 K和9.8 K，经过换热器降温后，一、二级换热器出口的空气温度增加了7.2 K和1.5 K，三、四级压缩机和换热器出口空气温度基本不变。这是因为随着压缩进气温度的增加，压缩机出口温度逐渐增大，进入换热器的空气温度也随之升高，但是经过4级换热器充分换热逐渐缩小了气体的温差，致使压缩机最后出口温度几乎不变。
Figure 2. The air temperature of the compressed section varies with the inlet air temperature
图3表示储热罐温度随压缩机进气温度的变化情况。随进气温度的升高，换热介质吸收的热量增加，储热罐中水的温度相应增加，由403.2 K增加到415.9 K。
Figure 3. The temperature of the hot water tank changes with the temperature of the inlet air
系统的储能效率和储能密度随压缩机入口温度变化如图4所示，系统储能效率略微下降，储能密度逐渐增大，这是因为储热罐中介质温度的提升增加了膨胀机的输出功，提高了储能密度，但压缩机入口温度的增加提高了压缩机耗功，且提高比例高于膨胀功，故降低了系统的储能效率。
Figure 4. Changes in energy storage efficiency and energy storage density of the system
-
在储能过程中，压缩机末级的出口压力要大于储气罐中的空气压力，但是过大的压力会造成压缩机的耗功增大。在储气罐最大储气压力不变的情况下，改变压缩机总压比，可计算出系统的储能效率和储能密度，结果如图5所示。
Figure 5. Variation of system performance indicators with the total pressure ratio of the compressor
随着压缩机总压比增大，系统的储能效率由52.9%下降到49.5%，储能密度由5.52 MJ/m3增加到5.62 MJ/m3，这是因为压缩机总压比增大，单位质量压缩空气的耗功增加，虽压缩机出口温度上升增加了换热介质的温度，增强了释能过程膨胀机入口空气的加热效果，单位质量的空气做功能力增强，但膨胀机输出功的增加幅度小于压缩功耗上升的幅度，所以系统的储能效率下降。系统储能密度的增加，是由于空气密度随压力上升而增加所致。
-
为了使系统输出功率恒定，不随储气罐压力降低而波动，可以在膨胀段前加节流阀（本文暂不考虑膨胀机滑压运行），高压空气经节流降至所需压力，再由换热器加热后进入膨胀机做功。在储气罐储气压力为5 MPa，忽略节流损失，分析不同节流后压力对系统性能影响，结果如图6、图7所示。
Figure 6. Variation of system performance with post-throttling pressure
Figure 7. Variation of energy storage efficiency with energy storage density
如图6所示，随着节流后压力逐渐增大，储能密度先升高后降低，这是由于节流后压力升高，单位质量空气输出的功就越多，但是从储气罐中释放的空气质量越来越少，膨胀机输出的总功先增加到最大值后逐渐减小，当节流后压力为1.35 MPa时，储能密度达到最大值8.15 MJ/m3。节流后压力越高，循环过程中储气罐释放的空气质量越少，系统的储能效率逐渐增大。图7为系统储能密度随着储能效率变化的曲线，可以看出随着储能效率增加，储能密度会先增大后减小，最大值对应的储能效率为40.17%。在优化系统效率时，也要权衡储能密度的变化，结合实际条件和需求选择最优设计点。
储气罐储气压力改变时，不同节流后压力对储能密度影响见图8所示，储能密度变化都是先增大到最大值后减小，当储能压力由3 MPa增加到6 MPa 时，系统最大储能密度由3.95 MJ/m3增加到10.50 MJ/m3，增加了2.66倍。
Figure 8. Variation of system energy storage density with post-throttling pressure
-
在CAES系统中，进行多级膨胀和级间再热有助于减少能量在膨胀过程中的损失，从而提高系统的效率。多级膨胀有诸多益处，但非级数越多越好，还要考虑级数增加的成本上升和连接管道增加带来的流动损失。在压缩机级数不变的情况下，改变膨胀级数，可分析在不同储能压力下膨胀级数对系统性能的影响。
图9所示，在相同储能压力下，随着膨胀级数的增加，系统储能效率逐渐增大，且增长速率逐渐变小。这是因为越多级的膨胀和再热，系统膨胀段越接近于等温膨胀，单位质量空气输出的功越多，但受限于储能时吸收的热量有限，换热器用于加热空气的能量就有限，所以系统储能效率的增幅减小。膨胀级数增加，系统输出功增加，储能密度也会增加。在设计工况下，随着储能压力由3 MPa增加到6 MPa，系统的储能效率由44.37%增长到53.39%，增长了9.02%，系统的储能密度增长了1.72倍。
Figure 9. Effect of the number of expansion stages on system performance
-
建立AA-CAES系统的动态模型，可更好反映系统中各关键物理量的动态变化在不同工况对系统性能的影响，其中储气罐的动态特性是系统动态分析的关键，通过微分方程构建储气罐动态模型，分析压力和温度的动态变化。系统储能阶段完成的标志是罐中压力达到设定的最大压力，此时系统给出反馈信号，压缩机停止工作，系统储能过程结束。在系统释能过程中，当罐中压力低于设定的最小压力，此时膨胀机停止工作，释能过程结束。分别改变储气罐的总传热系数、储气压力、储气温度，探究对系统性能的影响。
-
对于微小型压缩空气储能系统，储气装置通常采用钢制的储气罐，可以通过改变储气罐表面的保温层来改变储气罐与环境的换热系数。文献[25]中通过实验模拟出储气罐总传热系数为定值34.4 W/m2/K，本章节在此基础上探究不同换热情况对系统性能的影响，在设计工况下，由理想的绝热情况逐渐增加储气罐与环境的换热系数，选取换热系数为0 W/m2/K、5 W/m2/K、10 W/m2/K、20 W/m2/K、50 W/m2/K，分析储气罐换热对储气罐温度、压力以及系统性能的影响，见图10。储能过程中，绝热情况下储气罐内空气的温度随时间持续增长，由298 K增加到340.7 K。当储气罐与环境有换热时，储气罐温度随时间先增加后趋于平缓，不同换热系数下储气罐的最终温度不同，分别为314.5 K、307.1 K、302.7 K、299.9 K。储气罐压力随着时间近似线性增加，换热系数越小，压力增长的趋势越缓慢，储气时间越长。这是因为随着压缩空气不断进入储气罐，储气罐内空气被压缩，空气的温度和压力均增加，与外界环境的换热温差增大，换热量增大，最终温度趋于稳定值。与外界换热系数越大，交换的热量就越多，储气罐内压缩空气的终温就越接近环境温度。
Figure 10. Variation of temperature and pressure of gas storage tanks with different heat transfer coefficients during energy storage process
图11表示释能过程储气罐内压缩空气温度、压力随时间的变化曲线，储能过程结束时储气罐温度为释能过程的初始温度。可以看出，在储气罐绝热的情况下，随着释能过程的进行，储气罐温度持续下降，到释能过程结束时，温度下降为268.8 K。随着换热系数增大，储气罐内温度先缓慢下降后趋近于一个稳定值，且越大的换热系数对应的储气罐终温越接近于环境温度。储气罐压力随时间近似线性减少，当换热系数为5 W/m2/K时，压力曲线下降的斜率最大。当系统开始释能时，空气不断流出，储气罐内压力和温度都下降，随着释能过程的进行，储气罐内外温差逐渐加大，压缩空气与环境换热量增加，所以储气罐内空气温度最终趋于稳定。
Figure 11. Variation of temperature and pressure of gas storage tanks with different heat transfer coefficients during energy release
图12表示换热系数对系统储能效率和储能密度的影响。如图所示，储能密度随着换热系数的增大而增大，这是因为压缩空气与外界环境交换的能量越多，储能和释能过程储气罐终了温度越接近环境温度，则系统中循环的空气质量越多，所以储能密度会持续增加。而且越高的换热系数会降低储气罐储能终了温度，有利于储气罐的安全运行。系统的储能效率随着换热系数先降低后增加，在换热系数为5 W/m2/K有最小值，因为此时储气罐压力下降速率最快，释能时间最短，膨胀机输出功最小，随着换热系数继续增加，储气罐终温升高，释能时间延长，膨胀段输出更多的功，储能效率逐渐增加。
Figure 12. Variation of energy storage efficiency and energy storage density with heat transfer coefficient
-
在分析储气压力影响时，不考虑压力波动的影响。储气罐最高储气压力由4 MPa增加到8 MPa，储气罐换热系数取10 W/m2/K时，系统性能的变化如图13所示。
Figure 13. Impact of the system on the maximum pressure
由图13可见，随着最高储气压力增加，储能和释能时间增加，膨胀机总膨胀比增大，做功增加，储能密度上升。储能效率先上升后下降，压力为6 MPa时对应最大储能效率为50.46%。在储能效率最高点之前，随着储气罐最大压力的增加，膨胀机做功量的增长速率超过了压缩机耗功的增长速率，从而使得储能效率逐步提升，当储气罐的最大压力超过6 MPa后，压缩机的耗功增长显著，这导致了储能效率的下降。
当储气罐最高储气压力为5 MPa时，改变储气罐的最低储气压力，系统性能变化如图14所示。随着储气罐最低储气压力由0.5 MPa增加到4 MPa，系统的储能、释能时间近似线性下降，膨胀机做功先增大后减小，压缩机耗功持续下降。这是因为随着储气罐最低压力的增加，膨胀段的膨胀比增加，单位质量空气做功能力增大，但是从储气罐中释放的空气质量越来越少，所以膨胀功先增加后缓慢减小，系统的储能密度先增加后减小，在储气罐最低压力为1 MPa时，系统获得最大储能密度为7.48 MJ/m3。
Figure 14. Impact of system with minimum tank pressure
对于系统储能效率来说，储气罐最低压力增加，膨胀段压比增大，膨胀功增加，由于储气罐最高压力不变，储能过程压缩机的功率不变，所以系统的储能效率逐渐增加。
-
在末级压缩机与储气罐之间安装后冷器，进入储气罐前降低高压空气的温度，储气罐的初始温度为室温，取换热系数为10 W/m2/K，储气罐进气温度由298K增至318K，储气罐内空气的温度、压力变化见图15，储能效率和密度变化见图16。
Figure 15. Variation of temperature and pressure of gas storage tank with different inlet temperatures
Figure 16. Effect of storage tank inlet temperature on system performance
由图15、图16可知，进气温度越高，储气罐内空气压力增加越快，系统储、释能时间更短，储气罐的终温越高。当进气温度为298 K时，储气罐终温为307.1 K，进气温度增加到318 K时，储气罐终温为309.2 K，增加了2.1 K。随着储气罐进气温度增加，系统的储能效率基本不变，储能密度逐渐减小。这是因为越高的储气罐进气温度，系统储、释能时间越短，储气室终温越高，虽然较高的进气温度可以提高单位质量空气的做功能力，但由于储、释能时间的缩短，循环中用于做功的空气质量减少，这对储能密度的影响更大，从而使系统储能密度下降。
-
本文构建了压缩空气储能系统的静态模型和动态模型，以10 kW级别压缩空气储能系统为对象，仿真分析了系统性能变化趋势，主要结论如下：
1） 压缩机进口空气温度变化对前两级压缩和换热器出口温度影响较大，经过4级换热器充分换热后，对末级影响较小。随着压缩机进口温度升高，系统储能效率略微下降，储能密度呈上升趋势。
2） 压缩段总压比增加，系统的储能效率下降，但是储能密度上升；膨胀级数增加，系统的储能效率和储能密度增加，但受限于储能时吸收的热量有限，系统储能效率和储能密度增加受限。
3） 随着节流后压力的增加，系统储能密度先增加后逐渐减小，在设计工况下，当节流后压力为1.35 MPa时，储能密度达到最大值8.15 MJ/m3，储能效率持续增长但增长速率逐渐变缓。
4） 储气罐换热系数（即保温效果）对储气罐内气体温度、压力有直接影响，换热系数越大，储气罐终温越接近于环境温度；与绝热情况对比，当储气罐换热系数增加到50 W/m2/K时，系统储能密度增长56.17%，储能效率增长7.17%。
5） 增加储气室最大储气压力，系统储能效率先增加后减小，在6 MPa时有最大效率50.46%，系统储能密度逐渐增加；储气罐最低储气压力增加，系统的储能效率逐渐增加，储能密度呈先增加后减小的趋势。
Simulation and Dynamic Analysis of Small Advanced Insulated Compressed Air Energy Storage System
doi: 10.16516/j.ceec.2024-173
- Received Date: 2024-05-24
- Rev Recd Date: 2024-06-29
- Available Online: 2025-03-26
- Publish Date: 2025-03-30
-
Key words:
- compressed air energy storage system /
- dynamic variation characteristics /
- compressor /
- tank /
- energy storage efficiency /
- energy storage density
Abstract:
| Citation: | WANG Yuxuan, ZHANG Yufeng, LI Liansheng. Simulation and dynamic analysis of small advanced insulated compressed air energy storage system [J]. Southern energy construction, 2025, 12(2): 145-157 doi: 10.16516/j.ceec.2024-173
|
DownLoad:
