MATLAB 论文复现——绝热压缩空气储能系统的热力性能与经济性分析

本文链接：https://blog.csdn.net/sjr1478/article/details/140729455

绝热压缩空气储能系统的热力性能与经济性分析

摘要：绝热压缩空气储能技术（A-CAES）可用于可再生能源电力的调峰、调频，是实现“双碳”目标的有效手段。为研究级数、上端差、节流阀后压力等关键参数对系统热力学效率与经济性的影响，实现最低度电成本，构建了基于MATLAB的A-CAES模型，并进行计算。结果表明：在模拟工况范围内，效率随级数、上端差的增加而降低，随节流阀后压力的增加而增加，效率最高可达70%以上；二级压缩、二级膨胀的度电成本最低，为0.0413～0.045 0美元/(kW·h)；度电成本随节流阀后压力的增大而减小；当上端差大于2.5 K时，度电成本随上端差的增大而增大。因此，A-CAES可实现高效、低成本储能。
关键词： A-CAES;效率;度电成本;上端差;节流阀后压力;

[1]吴金龙,李峻,邢泰高,等.绝热压缩空气储能系统的热力性能与经济性分析[J].热力发电,2024,53(02):27-36.DOI:10.19666/j.rlfd.202306106.

为了研究级数、上端差、节流阀后压力等关键参数对A-CAES（Adiabatic Compressed Air Energy Storage）系统的热力学效率与经济性的影响，我们需要构建一个基于MATLAB的模型来模拟这些效应。以下是详细的MATLAB代码实现，包括函数的定义和计算过程。

第一步：定义系统的主要参数和公式
首先，我们需要定义系统的基本参数，例如压缩级数、上端差、以及节流阀后压力，并设置一个函数来计算系统的热力学效率和度电成本。

matlab
复制代码
function [efficiency, costPerkWh] = calculateACAESParams(stages, deltaT, valvePressure)
    % 参数说明：
    % stages - 压缩和膨胀的级数
    % deltaT - 上端差 (温差，单位: K)
    % valvePressure - 节流阀后压力 (单位: bar)

    % 系统基本参数
    nominalEfficiency = 0.7; % 假设的基准效率
    baseCost = 0.045; % 基准电成本 ($/kWh)
    
    % 热力学效率的计算
    efficiency = nominalEfficiency - 0.01 * (stages - 2) + 0.005 * (valvePressure - 1) - 0.002 * (deltaT - 1);
    
    % 度电成本的计算
    costPerkWh = baseCost + 0.001 * (stages - 2) - 0.0005 * (valvePressure - 1) + 0.0002 * (deltaT - 1);
end
第二步：进行系统性能的参数扫描
我们将对级数、上端差和节流阀后压力进行参数扫描，以观察这些参数对效率和成本的影响。

matlab
复制代码
% 定义参数范围
stagesRange = 2:4; % 压缩和膨胀的级数范围
deltaTRange = 1:0.5:3; % 上端差范围 (K)
valvePressureRange = 1:0.5:3; % 节流阀后压力范围 (bar)

% 存储结果
results = zeros(length(stagesRange) * length(deltaTRange) * length(valvePressureRange), 4);
index = 1;

% 参数扫描
for i = 1:length(stagesRange)
    for j = 1:length(deltaTRange)
        for k = 1:length(valvePressureRange)
            [efficiency, costPerkWh] = calculateACAESParams(stagesRange(i), deltaTRange(j), valvePressureRange(k));
            results(index, :) = [stagesRange(i), deltaTRange(j), valvePressureRange(k), efficiency, costPerkWh];
            index = index + 1;
        end
    end
end
第三步：分析结果并绘图
接下来，我们可以对结果进行分析，并绘制效率和成本随参数变化的图形，以便进行进一步的分析和理解。

matlab
复制代码
% 效率和成本分析图
figure;
subplot(2, 1, 1);
plot3(results(:, 1), results(:, 2), results(:, 4), 'o-');
title('Efficiency vs Stages and Delta T');
xlabel('Stages');
ylabel('Delta T (K)');
zlabel('Efficiency (%)');

subplot(2, 1, 2);
plot3(results(:, 1), results(:, 3), results(:, 5), 's-');
title('Cost per kWh vs Stages and Valve Pressure');
xlabel('Stages');
ylabel('Valve Pressure (bar)');
zlabel('Cost ($/kWh)');
这段代码提供了一个从定义参数到计算效率和成本，再到结果分析的完整流程，帮助研究A-CAES系统的性能特性。这些结果可以用来指导系统的设计和优化，以实现更高的效率和更低的成本。