MATLAB 论文复现——集成恒压压缩空气储能的燃气轮机CCHP系统特性研究

集成恒压压缩空气储能的燃气轮机CCHP系统特性研究

摘要：燃气轮机冷热电三联供(CCHP)技术发展前景广阔，但存在部分负荷运行效率不高、热-电输出耦合性强、热电比调节不灵活等问题。构建了集成恒压压缩空气储能(CAES)的新型燃气轮机CCHP系统，建立了压气机、燃烧室、透平、余热锅炉、换热器、喷射器、水泵/水轮机等部件的热力学模型，分析了采用抽气-释气流量调节策略的系统运行特性。结果表明：抽气、释气系数在0～0.2内时，新型系统热电比调节范围为0.58～2.27,一次能源利用率稳定在58.1%～59.5%,■效率不低于32.2%,可实现宽工况高效运行。

关键词：    压缩空气储能;燃气轮机;CCHP;调节策略;

[1]张斌,王光磊,徐震,等.集成恒压压缩空气储能的燃气轮机CCHP系统特性研究[J].动力工程学报,2024,44(03):369-375.DOI:10.19805/j.cnki.jcspe.2024.230592.
 

为了模拟和分析集成恒压压缩空气储能（CAES）的新型燃气轮机CCHP系统，并研究抽气-释气流量调节策略对系统性能的影响，下面将详细介绍MATLAB代码的编写过程。

第一步：定义系统组件的热力学模型
首先，我们需要为燃气轮机CCHP系统中的各个组件建立热力学模型。我们将定义压气机、燃烧室、透平、余热锅炉、换热器、喷射器以及水泵/水轮机的基本功能。

matlab
复制代码
function [workCompressor, heatCombustor, workTurbine, heatRecBoiler, heatExchanger, workInjector, workPump] = componentModels(airFlow, bleedRatio, heatInput)
    % 压气机模型
    workCompressor = 0.8 * airFlow;  % 假设压气机工作效率为80%
    
    % 燃烧室模型
    heatCombustor = 0.9 * heatInput;  % 假设燃烧效率为90%
    
    % 透平模型
    workTurbine = 1.2 * airFlow;  % 假设膨胀做功系数为1.2
    
    % 余热锅炉模型
    heatRecBoiler = 0.5 * heatCombustor;  % 假设余热回收效率为50%
    
    % 换热器模型
    heatExchanger = 0.4 * heatRecBoiler;  % 假设换热效率为40%
    
    % 喷射器模型
    workInjector = bleedRatio * workTurbine;  % 假设喷射功率与透平功率成比例
    
    % 水泵/水轮机模型
    workPump = 0.7 * (workCompressor - workInjector);  % 假设水泵效率为70%
end
第二步：系统性能分析函数
接下来，根据定义的模型计算整个系统在不同抽气和释气系数下的性能。

matlab
复制代码
function [thermalElectricRatio, totalEfficiency, exergyEfficiency] = systemPerformance(airFlow, heatInput, bleedRatio)
    % 调用组件模型
    [workCompressor, heatCombustor, workTurbine, heatRecBoiler, heatExchanger, workInjector, workPump] = componentModels(airFlow, bleedRatio, heatInput);
    
    % 总能量输出
    totalEnergyOutput = workTurbine + heatRecBoiler + heatExchanger;
    
    % 热电比
    thermalElectricRatio = (heatRecBoiler + heatExchanger) / workTurbine;
    
    % 一次能源利用率
    totalEfficiency = totalEnergyOutput / heatInput;
    
    % ■效率（Exergy效率）
    exergyEfficiency = (workTurbine - workCompressor + workPump) / heatInput;
end
第三步：系统性能评估
创建一个高级函数来评估系统在不同操作条件下的性能，包括热电比调节、一次能源利用率和■效率。

matlab
复制代码
function evaluateCCHPSystem()
    % 假设流量和热输入参数
    airFlow = 100; % 假设空气流量（单位：任意）
    heatInput = 1000; % 热输入（单位：kJ）
    bleedRatios = linspace(0, 0.2, 20); % 抽气系数从0到0.2
    
    % 存储结果
    results = zeros(length(bleedRatios), 3);
    
    for i = 1:length(bleedRatios)
        [thermalElectricRatio, totalEfficiency, exergyEfficiency] = systemPerformance(airFlow, heatInput, bleedRatios(i));
        results(i, :) = [thermalElectricRatio, totalEfficiency, exergyEfficiency];
    end
    
    % 绘制结果
    figure;
    plot(bleedRatios, results(:, 1), 'b-o', 'DisplayName', '热电比');
    hold on;
    plot(bleedRatios, results(:, 2) * 100, 'r-s', 'DisplayName', '一次能源利用率 (%)');
    plot(bleedRatios, results(:, 3) * 100, 'g-^', 'DisplayName', '■效率 (%)');
    title('CCHP系统性能分析');
    xlabel('抽气系数');
    ylabel('性能指标');
    legend;
    grid on;
end
第四步：运行评估函数
最后，在MATLAB中运行evaluateCCHPSystem函数来查看系统性能随抽气系数变化的效果。

matlab
复制代码
evaluateCCHPSystem();
这段代码提供了一个全面的框架来模拟和分析集成CAES的新型燃气轮机CCHP系统，帮助研究其在不同工况下的性能表现。