储能策略的微网可靠性程序.rar_Microgrid_储能_微电网 可靠性_微电网储能_评估指标

load('data3'); MCTime=1; %系统模拟时间 Pbat=zeros(1,8760); %储能充放电量 lshed=zeros(8,8760);%初始削减负荷为空集,纵向为负荷编号，横向为时间顺序 lshedtdcs=zeros(1,8);%累计削减负荷的停电次数 lshedtdsj=zeros(1,8);%累计削减负荷的停电时间 %剩余负荷为所有负荷lremain矩阵% Tttf=zeros(1,15); Tttr=zeros(1,1); TIME=8760; N=100; Ebatmax=1.5; system=zeros(1,8760); Ebat=zeros(1,8760); Pl=zeros(1,8760);%总负荷 Alshedtdsj=zeros(1,8); Alshedtdcs=zeros(1,8); for i=1:8760 Pl(1,i)=LP1(1,i)+LP2(1,i)+LP3(1,i)+LP4(1,i)+LP5(1,i)+LP6(1,i)+LP7(1,i)+LP8(1,i); end for k=1:N %------------------------时序模拟start-------------------------% while MCTime<TIME %-----对元件正常工作时间抽样------% for i=1:15 x=rand(1,15);%(0,1)之间服从均匀分布的随机数 Tttf(1:i)=-log(x(i))/a(i); end Tw=min(Tttf);%选取故障元件 for i=1:15 if Tw==Tttf(1:i) n=i;%故障元件编号 break; end end y=rand();%(0,1)之间服从均匀分布的随机数 Tttr(1,1)=-log(y)/b(n);%抽样一个故障元件的修复时间 %-------对负荷进行分块，a类负荷记为0，b类负荷记为1----------% if n==1 for i=1:3 for j=1:4 lremain(i,MCTime+j-1)=0; end end elseif n==2 for j=1:4 lremain(3,MCTime+j-1)=0; end elseif n==3 for i=4:8 for j=1:4 lremain(i,MCTime+j-1)=0; end end elseif n==4 for j=1:4 lremain(7,MCTime+j-1)=0; lremain(8,MCTime+j-1)=0; end elseif n==5||n==6||n==7||n==8||n==9||n==10||n==11||n==12 for j=1:200 lremain(n-4,MCTime+j-1)=0; end elseif n==13 for j=1:12 Pd(1,MCTime+j-1)=0; end elseif n==14 for j=1:60 wind(1,MCTime+j-1)=0; end elseif n==15 %直接削减负荷 for t=1:50 c=t+MCTime; if c>8760 break; end if wind(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime) system(1,t+MCTime)=0;%系统可靠,削减了0次负荷 else lshed(7,t+MCTime)=1;%削减负荷78 lshed(8,t+MCTime)=1; lremain(7,t+MCTime)=0; lremain(8,t+MCTime)=0; end if wind(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime)-LP7(1,t+MCTime)-LP8(1,t+MCTime) system(1,t+MCTime)=1;%系统可靠，削减了1次负荷 else lshed(3,t+MCTime)=1;%削减负荷3 lremain(3,t+MCTime)=0; end if wind(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime)-LP7(1,t+MCTime)-LP8(1,t+MCTime)-LP3(1,t+MCTime) system(1,t+MCTime)=2;%系统可靠,削减了2次负荷 else lshed(1,t+MCTime)=1;%削减负荷12 lshed(2,t+MCTime)=1; lremain(1,t+MCTime)=0; lremain(2,t+MCTime)=0; end if wind(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime)-LP7(1,t+MCTime)-LP8(1,t+MCTime)-LP3(1,t+MCTime)-LP1(1,t+MCTime)-LP2(1,t+MCTime) system(1,t+MCTime)=3;%系统可靠,削减了3次负荷 else lshed(4,t+MCTime)=1;%削减负荷456 lshed(5,t+MCTime)=1; lshed(6,t+MCTime)=1; lremain(4,t+MCTime)=0; lremain(5,t+MCTime)=0; lremain(6,t+MCTime)=0; system(1,t+MCTime)=4;%系统可靠,削减了4次负荷 end %--------累计削减负荷的停电时间-----------% for i=1:8 if lshed(i,t+MCTime)==1 lshedtdsj(1,i)=lshedtdsj(1,i)+1; end end end %--------累计削减负荷的停电次数=停电总时间-停电连续时间-------% lshedtdcs=lshedtdsj; for i=1:8 for j=1:8759 if lshed(i,j)==1 if lshed(i,j+1)==1 lshedtdcs(1,i)=lshedtdcs(1,i)-1; end end end end end%-------负荷分块结束 t=-1; %------储能运行策略2--------% while t<Tw+Tttr t=t+1; if t+MCTime>8760 break; end if wind(1,t+MCTime)>=Pl(1,t+MCTime)*0.3 Pbat(1,t+MCTime)=wind(1,t+MCTime)-Pl(1,t+MCTime)*0.3; Ebat(1,t+MCTime+1)=Ebat(1,t+MCTime)+Pbat(1,t+MCTime); if Pd(1,t+MCTime)>=Pl(1,t+MCTime)*0.7 system(1,t+MCTime)=0;%系统可靠,削减了0次负荷 else %-------------------削减负荷---------------------% if Pwbat(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime) system(1,t+MCTime)=0;%系统可靠,削减了0次负荷 else lshed(7,t+MCTime)=1;%削减负荷78 lshed(8,t+MCTime)=1; lremain(7,t+MCTime)=0; lremain(8,t+MCTime)=0; end if Pwbat(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime)-LP7(1,t+MCTime)-LP8(1,t+MCTime) system(1,t+MCTime)=1;%系统可靠，削减了1次负荷 else lshed(3,t+MCTime)=1;%削减负荷3 lremain(3,t+MCTime)=0; end if Pwbat(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime)-LP7(1,t+MCTime)-LP8(1,t+MCTime)-LP3(1,t+MCTime) system(1,t+MCTime)=2;%系统可靠,削减了2次负荷 else lshed(1,t+MCTime)=1;%削减负荷12 lshed(2,t+MCTime)=1; lremain(1,t+MCTime)=0; lremain(2,t+MCTime)=0; end if Pwbat(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime)-LP7(1,t+MCTime)-LP8(1,t+MCTime)-LP3(1,t+MCTime)-LP1(1,t+MCTime)-LP2(1,t+MCTime) system(1,t+MCTime)=3;%系统可靠,削减了3次负荷 else lshed(4,t+MCTime)=1;%削减负荷456 lshed(5,t+MCTime)=1; lshed(6,t+MCTime)=1; lremain(4,t+MCTime)=0; lremain(5,t+MCTime)=0; lremain(6,t+MCTime)=0; system(1,t+MCTime)=4;%系统可靠,削减了4次负荷 end %--------累计削减负荷的停电时间-----------% for i=1:8 if lshed(i,t+MCTime)==1 lshedtdsj(1,i)=lshedtdsj(1,i)+1; end end %--------累计削减负荷的停电次数=停电总时间-停电连续时间-------% lshedtdcs=lshedtdsj; for i=1:8 for j=1:8759 if lshed(i,j)==1 if lshed(i,j+1)==1 lshedtdcs(1,i)=lshedtdcs(1,i)-1; end end end end end else if wind(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime) system(1,t+MCTime)=0;%系统可靠,削减了0次负荷 else Pbat(1,t+MCTime)=min(0.25*Ebatmax,Pl(1,t+MCTime)*0.3-wind(1,t+MCTime)); Ebat(1,t+MCTime+1)=Ebat(1,t+MCTime)-Pbat(1,t+MCTime); if wind(1,t+MCTime)+Pd(1,t+MCTime)+Pbat(1,t+MCTime)>=Pl(1,t+MCTime) system(1,t+MCTime)=0;%系统可靠,削减了0次负荷 else %----------------------------削减负荷-----------------------------% if Pwbat(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime) system(1,t+MCTime)=0;%系统可靠,削减了0次负荷 else lshed(7,t+MCTime)=1;%削减负荷78 lshed(8,t+MCTime)=1; lremain(7,t+MCTime)=0; lremain(8,t+MCTime)=0; end if Pwbat(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime)-LP7(1,t+MCTime)-LP8(1,t+MCTime) system(1,t+MCTime)=1;%系统可靠，削减了1次负荷 else lshed(3,t+MCTime)=1;%削减负荷3 lremain(3,t+MCTime)=0; end if Pwbat(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime)-LP7(1,t+MCTime)-LP8(1,t+MCTime)-LP3(1,t+MCTime) system(1,t+MCTime)=2;%系统可靠,削减了2次负荷 else lshed(1,t+MCTime)=1;%削减负荷12 lshed(2,t+MCTime)=1; lremain(1,t+MCTime)=0; lremain(2,t+MCTime)=0; end if Pwbat(1,t+MCTime)+Pd(1,t+MCTime)>=Pl(1,t+MCTime)-LP7(1,t+MCTime)-LP8(1,t+MCTime)-LP3(1,t+MCTime)-LP1(1,t+MCTime)-LP2(1,t+MCTime) system(1,t+MCTime)=3;%系统可靠,削减了3次负荷 else lshed(4,t+MCTime)=1;%削减负荷456 lshed(5,t+MCTime)=1; lshed(6,t+MCTime)=1; lremain(4,t+MCTime)=0; lremain(5,t+MCTime)=0; lremain(6,t+MCTime)=0; system(1,t+MCTime)=4;%系统可靠,削减了4次负荷 end %--------累计削减负荷的停电时间-----------% for i=1:8 if lshed(i,t+MCTime)==1 lshedtdsj(1,i)=lshedtdsj(1,i)+1; end end %--------累计削减负荷的停电次数=停电总时间-停电连续时间-------% lshedtdcs=lshedtdsj; for i=1:8 for j=1:8759 i