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

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'); %火电机组事故旋转备用容量
 

3.运行结果