👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》
论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html
论文参考:
计及调峰影响系数的费用分摊及随机经济调度
计及调峰主动性的风光水火储多能系统 互补协调优化调度
考虑广义储能与火电联合调峰的日前–日内 两阶段滚动优化调度
针对风电、光伏等大规模可再生能源并网后电力系统调峰能力不足的问题,在分析调峰补偿与调峰分摊的基础上,考虑火电机组调峰主动性约束,提出一种风光水火储多能系统互补协调优化调度略。考虑多能系统电源结构复杂,涉及变量及约束条件较多,因此采用分层优化调度方案。上层模型以净负荷波动最小和储能系统运行收益最大为优化目标,旨在充分利用储能装置削峰填谷特性,降低负荷峰谷差,提高可再生能源的消纳空间;下层模型以火电机组运行成本最小和可再生能源弃电量最小为优化目标,考虑调峰主动性约束,旨在充分发挥火电机组深度调峰能力,优化可再生能源消纳能力和火电机组经济运行。基于分解协调思想,实现上、下层问题的协调和下层问题的交替迭代求解。最后,以改进的IEEE30节点系统为例,进行多种场景的仿真计算分析,结果表明所提策略能有效提升系统可再生能源的消纳能力和系统运行的经济性,验证了该模型的有效性。
部分代码展示:
% %% 导入ieee30节点网络
% caseName = case30;
% fsc=5.6196;
%% 参数
iv = 10; %光伏接入节点
iw = 3; %风电接入节点
ie = 15; %EES接入节点
idg = [1;2;5;8;11]; %DG接入节点
is=13;
igd = [1]; %大电网接入节点
Horizon = 96;
nbus = size(caseName.bus, 1);
ngen=5;%发电机数量
nv = 1;%光伏数量
nw = 1;%风电数量
nees=1;%储能数量
%电价
%global prg prw prgd price
prg=375;prw=570;qw=536;qv=536;prv=921;%发电机上网价格 风电上网 弃风 弃光 光伏上网价格
prgd=500;%深度调峰
Pdeepdn=100;%调峰下限
Pdeepup=2000;%调峰上限
alfaqf=500;
yitac=0.9;
yitad=0.9;
price=[0.3177.*ones(1,24),0.6062.*ones(1,10),0.8948.*ones(1,10),0.6062.*ones(1,4),0.3177.*ones(1,4),0.6062.*ones(1,6),0.8948.*ones(1,26),0.6062.*ones(1,8),0.3177.*ones(1,4)];
%L_Horizon = r_base.*repmat(L_base,1,Horizon) + r_fluc.*repmat(L_base,1,Horizon).*repmat(L_profile,nbus,1); %IEEE30节点的负荷
%整体负荷
plt=[615 611 605 603 603 612 613 615 616 617 618 619 632 645 658 669 675 689 691 701 712 713 756 769 789 812 845 863 852 861 872 879 882 891 905 912 924 935 945 951 962 974 982 990 1052 1152 1061 1032 1021 1016 996 985 963 942 931 921 910 892 885 872 865 851 840 831 820 812 802 792 782 764 753 744 738 725 716 706 715 726 735 749 755 768 771 785 792 806 816 823 829 830 821 811 854 831 765 700];
ru=plt.*0.1;%备用
%% 决策变量
x_theta = sdpvar(nbus, Horizon,'full');%直角坐标角度
x_P_grid = sdpvar(1, Horizon,'full');
x_P_g = sdpvar(ngen, Horizon,'full');%发电机
x_P_ch = sdpvar(nees, Horizon);%蓄电池
x_P_dis = sdpvar(nees, Horizon);%蓄电池
x_P_w = sdpvar(nw, Horizon,'full');%风电
x_P_v = sdpvar(nv, Horizon,'full');%光伏
x_u_g = binvar(ngen, Horizon,'full');%发电机状态
x_u_ch = binvar(nees, Horizon);%蓄电池状态,下同
x_u_dis = binvar(nees, Horizon);
x_P_sg=sdpvar(1, Horizon,'full');%水电
gn=5;%分段函数线性化,下同
x_pf=sdpvar(ngen, Horizon,'full');
gw1=sdpvar(gn+1,Horizon,'full');
gw2=sdpvar(gn+1,Horizon,'full');
gw3=sdpvar(gn+1,Horizon,'full');
gw4=sdpvar(gn+1,Horizon,'full');
gw5=sdpvar(gn+1,Horizon,'full');
gw6=sdpvar(gn+1,Horizon,'full');
gz1=binvar(gn, Horizon,'full');gz2=binvar(gn, Horizon,'full');gz3=binvar(gn, Horizon,'full');gz4=binvar(gn, Horizon,'full');gz5=binvar(gn, Horizon,'full');
%% 约束条件生成
cons = [];
% DG
Pmax = [460;300;243;120;130];%最大发电功率
Pmin = 0.5.*Pmax;%最小发电功率
rg = 0.5.*Pmax;%爬坡
Pa=0.45.*Pmax;
Pb=0.3.*Pmax;
Psmax=100;%水电机组
效果展示:
123号资源-源程序:论文可在知网下载《计及调峰主动性的风光水火储多能系统互补协调优化调度》本人博客资源-CSDN文库https://download.csdn.net/download/LIANG674027206/89401591👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》
论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html