2.文章复现《热电联产系统在区域综合能源系统中的定容选址研究》（附matlab程序）

news2025/1/11 6:14:37

0.代码链接

1.简述

光热发电是大规模利用太阳能的新兴方式，其储热系统能够调节光热电站的出力特性，进而缓解光热电站并网带来的火电机组调峰问题。合理配置光热电站储热容量，能够有效降低火电机组调峰成本。该文提出一种光热电站储热容量配置方法，在对光热电站运行特性及其对电网调峰影响分析的基础上，综合考虑了火电机组向下调峰成本与储热成本对储热系统容量配置的影响。最后分别基于 IEEE 30 节点系统以及西北某省电网系统的仿真算例，验证了该文方法的可行性和有效性。

关键词：光热电站；储热系统；运行特性；调峰成本；最优储热容量配置

2.部分代码

clc
clear
tic
%% 参数定义
T=24; %调度时刻
yita_c = 0.9; %热储能充热效率
yita_f = 0.9; %热储能放热效率
yita_d = 0.4; %储热系统热电转换效率
N_i=6; %火电机组台数
a=[0.2625 1.224 4.375 0.5838 1.75 1.75];%机组燃料成本系数
b=[140 122.5 70 227.5 210 210];
S=[150 141 121 133 111 141]; %启停成本
P_max=[100 80 50 35 30 40]; %火电机组出力上限
P_min=[53 41 26 18 16 21]; %火电机组出力下限
k_G=230; %环境效益系数230元/MW
alpha=157; %正旋转备用系数
beta=134; %负旋转备用系数
gama=132; %事故旋转备用系数
k_s=120; %集热装置供热发电的成本系数
k_Ts=100; %储热装置供热发电的成本系数
c=420000; %单位储热容量配置成本
mpc=case30;
l=sum(mpc.bus(:,3));
A=[0.95 0.965 0.925 0.945 0.965 0.976 0.953 0.942 0.934 0.943 0.934 0.98 0.96 0.93 0.94 0.967 0.935 0.95 0.98 0.99 0.97 0.94 0.95 0.97];
load = A*l;
% load=[220 230 230 230 198 220 208 250 280 303 300 259 255 230 235 253 290 308 307 280 274 260 252 238];
%% 变量定义
P_cspr=[10.05 10.05 10.05 10.05 10.89 10.05 10.89 10.05 84.59 386.10 427.97 383.58 427.14 505.03 355.11 11.73 87.94 10.89 10.05 10.05 10.05 10.05 10.05 10.05]; %CSP集热装置吸收热功率
P_SFr=sdpvar(1,T,'full'); %集热装置直接发电的热功率
P_HT=sdpvar(1,T,'full'); %集热装置供给储热系统的热功率
P_SFd=sdpvar(1,T,'full'); %集热装置输出电功率
P_TH=sdpvar(1,T,'full'); %储热装置传递至传热流体的热功率
P_TSc=sdpvar(1,T,'full'); %储热装置的储热功率
P_TSf=sdpvar(1,T,'full'); %储热装置的放热功率
u_TSc=binvar(1,T,'full'); %储热装置的储热状态
u_TSf=binvar(1,T,'full'); %储热装置的放热状态
P_TSdf=sdpvar(1,T,'full'); %储热装置输出电功率
C_TSr=sdpvar(1,T,'full'); %储热量
P_G=sdpvar(1,T,'full'); %光热电站的输出功率
P=sdpvar(N_i,T,'full'); %火电机组功率
u=binvar(N_i,T,'full'); %火电机组运行状态
U=sdpvar(N_i,T,'full'); %火电机组正旋转备用容量
D=sdpvar(N_i,T,'full'); %火电机组负旋转备用容量
R=sdpvar(N_i,T,'full'); %火电机组事故旋转备用容量