目录
主要内容
该模型采用序贯蒙特卡洛模拟法对电力系统网络的可靠性进行评估,在可靠性评估方向,序贯蒙特卡洛模拟法是常用方法,该方法是一种全值估计方法,可以处理非线性、大幅波动及“肥尾”问题,比历史模拟方法更精确和可靠,程序采用matlab编写,利用节点影响分析法判断受影响的负荷,通过序贯蒙特卡洛算法完成配电网可靠性评估,运行结果包括各节点年平均停电时间、次数等评估结果以及误差分析,程序注释清楚,方便学习!
程序代码
%统计故障次数 load_fault=load_fault+equipment_load(index,:)+operator_load(index,:); %统计故障时间 load_fault_time=load_fault_time+TTR1*equipment_load(index,:)+TTR2*operator_load(index,:); %统计失负荷量 load_power=load_power+power'.*(equipment_load(index,:)+operator_load(index,:)); %累加无故障工作时间 load_normal=load_normal+min_TTF; %累加仿真时间 time=time+min_TTF+max([TTR1,TTR2]); end %负荷停电次数 Lambda_load=load_fault./load_normal; %负荷停电时间 Gamma_load=load_fault_time./load_fault; U_load=Lambda_load.*Gamma_load; disp('**************系统可靠性指标**************') SAIFI=Lambda_load*load_users/(sum(load_users)); disp(['SAIFI=',num2str(SAIFI),'次/(年·户)']) SAIDI=U_load*load_users/(sum(load_users)); CAIDI=SAIDI/SAIFI; disp(['SAIDI=',num2str(SAIDI),'小时/(年·户)']) EENS=load_power*load_users/(sum(load_users))/N; disp(['EENS=',num2str(EENS),'MW/(年·户)']) ASAI=1-U_load*load_users/(sum(load_users*8760)); disp(['ASAI=',num2str(ASAI*100),'%']) figure(1) bar(Lambda_load,'stack') title('负荷可靠性评估结果') xlabel('负荷节点');ylabel('年平均停电次数/次'); figure(2) bar(U_load,'stack') title('负荷可靠性评估结果') xlabel('负荷节点');ylabel('年平均停电时间/h'); Z=[SAIFI SAIDI CAIDI ASAI]; figure(3) bar(Z,0.2) title('系统可靠性评估结果') set(gca,'XTickLabel',{'SAIFI','SAIDI','CAIDI','ASAI'}); Lambda_load0=[1.6725 1.6725 1.6725 1.6725 1.6725 1.7115 1.7213 2.5370 2.5890 2.5370 2.5370 2.5370 1.6725 1.7115 1.6725 2.5110 2.5598 2.5110 2.5110 2.5110 2.2250 2.2250 2.2250 ]; U_load0=[8.4015 8.4015 8.4015 8.4015 8.4015 8.5965 8.6453 12.7240 12.9840 12.7240 12.7240 12.7240 11.2875 11.4825 11.2875 15.4800 15.7238 15.4800 15.4800 15.4800 14.0500 14.0500 14.0500 ]; SAIFI0=1.97781318681319; SAIDI0=11.074659340659341; CAIDI0=5.599446608253094; ASAI0=0.998735769481660; figure(4) bar(Lambda_load-Lambda_load0,'stack') title('模拟法误差分析-负荷年平均停电次数') xlabel('负荷节点');ylabel('年平均停电次数的误差'); figure(5) bar(U_load-U_load0,'stack') title('模拟法误差分析-负荷年平均停电时间') xlabel('负荷节点');ylabel('年平均停电时间的误差'); figure(6) Z0=[SAIFI0 SAIDI0 CAIDI0 ASAI0]; bar(Z-Z0,0.4) title('系统可靠性评估误差') set(gca,'XTickLabel',{'SAIFI误差','SAIDI误差','CAIDI误差','ASAI误差'});