👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》
论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html
论文可参考:
考虑多能灵活性的综合能源系统多时间尺度优化调度
计及碳捕集电厂低碳特性的含风电电力系统 源–荷多时间尺度调度方法
含电热气多种能源的综合能源系统中,复杂的能量转换关系以及可再生能源和负荷的波动 性,给综合能源系统的灵活安全运行带来了挑战。为了减小新能源和负荷不确定性对系统的影响,同时提升系统的功率平抑能力,提出了考虑多能灵活性的综合能源系统多时间尺度优化调度策略。首先,在日前调度中建立多能灵活性状态方程,以日运行成本最小为目标,构建波动场景下考虑系统多能灵活性的调度模型;日内调度则以日前调度计划为基础,考虑电能和热气能在不同时间尺度上的响应能力,建立分时间尺度的滚动优化模型,对日前计划进行修正,平抑功率波动。算例结果表明:所提出的日前模型可以通过协调设备的出力来提升系统的运行灵活性,日内模型能够平抑电热气能在不同时间尺度上的功率波动。
部分代码展示:
clear all
%% 输入原始数据 光照强度:w/m^2 风速:m/s 电热负荷:kw等
Solar=[0 0 0 0 0 9 115 282 454 620 744 801 793 709 593 457 297 127 15 0 0 0 0 0];%光照强度
Wind=[7.14 7.20 7.22 6.89 7.18 7.07 6.43 5.97 6.33 6.56 6.63 6.70 6.65 6.92 6.86 6.85 6.89 7.01 7.02 6.97 7.05 7.06 7.03 7.10];%风速
Load_E =[30420.49 32008.64 32673.34 33327.69 29716.69 29540.05 27917.86 37536.115 44387.93 42669.255 44470.04 32343.06 34827.06 41421.16 41072.135 43739.79 38737.75 33883.05 36011.1 34328.19 34977.48 32621.59 32920.36 32306.72];%电负荷
Load_H =[27394.52 35116.88 36361.76 37704.68 38214.26 40975.34 42859.76 50462.42 55585.58 58076.48 59895.92 56950.16 57761.84 58141.46 57643.28 59874.26 55563.92 55509.2 48967.88 44603.96 43942.76 36773.13 34553.72 30860.12];%热负荷
Price_E=[0.35 0.35 0.35 0.35 0.35 0.35 0.35 0.61 0.97 0.97 0.97 0.61 0.35 0.61 0.97 0.97 0.97 0.97 0.97 0.97 0.97 0.61 0.61 0.35];%分时电价,元/kwh
Price_Gas=2.7;%天然气价格,元/m^3
E=0;
a_Gas=35600;%天然气热值,35800kJ/m^3
alpha=0.08;%碳捕集设备选择捕集产生的二氧化碳比例, 燃烧产生的二氧化碳不一定要全部捕集,这里0.08=8%,至少捕集8%产生的二氧化碳
%% 输入设备模型
%光伏
S_PV=40000;%光伏面积
Q_PV=0.15*S_PV;%一平米折合约0.15kw
Pe_PV=0.001*0.157*S_PV*Solar;%光伏功率,效率0.157,0.001将光照强度:w/m^2转换成kw/m2,光照强度已载入
PV_omfix=0.02*Q_PV;%固定运维成本,按天结算,如设备清灰这样的人工费
PV_omwork=0.039*sum(Pe_PV);%可变运维成本,与发电量有关
PV_om=PV_omfix+PV_omwork;%总运维成本
%风电WP
N_WP=4;%总台数
Q_sWP=2000;%单台额定容量2MW
for i=1:24 %风机,单台额定2MW,切入风速3,额定12,切出20
if 3<=Wind(i)&&Wind(i)<=12
Pe_swind(i)=222.2*(Wind(i)-3); %2000=2220*(12-3)因为上限容量为2MW,所以风速大于12时满额输出
elseif 12<Wind(i)&&Wind(i)<=20
Pe_swind(i)=Q_sWP;
else
Pe_swind(i)=0; %其余情况风力不发电
end
end
Pe_WP=N_WP*Pe_swind;%总出力
WP_om=0.07*sum(Pe_WP);%风电运维
%热电联产CHP
Q_CHP=25000;%CHP容量,kw
Qg_CHP= sdpvar(1, 24);%CHP耗气速率,m3/s sdpvar-创建实数型决策变量,等待求解
Prg_CHP=Qg_CHP*a_Gas;%CHP燃烧天然气产生功率
Pe_CHP=0.31*Prg_CHP;%CHP输出电功率,产电效率0.31
Ph_CHP=0.5*(1-0.31)*Prg_CHP;%CHP产热功率;余热回收效率0.5 将发电剩下的能量进行余热回收供热,能量阶梯化
CHP_omfix=0.02*Q_CHP;
CHP_omwork=0.03*sum(Pe_CHP);
CHP_om=CHP_omfix+CHP_omwork;
%燃气锅炉GB
Q_GB=18000;%燃气锅炉容量,kw
Qg_GB= sdpvar(1, 24);%燃气锅炉耗气速率,m3/s
Ph_GB=0.9*Qg_GB*a_Gas;%燃气锅炉产热,效率0.9
GB_omfix=0.03*Q_GB;
GB_omwork=0.02*sum(Ph_GB);
GB_om=GB_omfix+GB_omwork;
%电锅炉
Q_EB=16000;%电锅炉容量,kw
Pe_EB=sdpvar(1,24);%电锅炉电功率 例如:x = sdpvar(1,1); %x为1×1的实数变量
Ph_EB=0.95*Pe_EB;%电热锅炉产热,效率0.95
EB_omfix=0.03*Q_EB;
EB_omwork=0.04*sum(Ph_EB);
EB_om=EB_omfix+EB_omwork;
效果展示:
109号资源-源程序:论文可在知网下载《综合能源系统低碳优化调度》本人博客有解读资源-CSDN文库https://download.csdn.net/download/LIANG674027206/89315328👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》
论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html