含风电的发电系统可靠性评估(MC法)matlab程序

%% 3.计算含风电场的发电系统可靠性指标（非序贯MC）

clc

clear

loadresult_WindFarmOutput

%文件“result_WindFarmOutput.mat”构成了风电场出力的状态模型【风力状态状态概率】相关状态计算查看百度文库“风电场出力模型matlab程序”

% 3.1 求出常规机组的出力模型，按类构成多状态模型

% RBTS发电系统中共有6类常规机组，%11台常规机组数据

% %2台5MW水电机组%% %1台10MW热电机组%% %4台20MW水电机组%% %1台20MW 热电机组%% %1台40MW水电机组%% %2台40MW热电机组%

Generator.Norm=[5 0.01

5 0.01

10 0.02

20 0.015

20 0.015

20 0.015

20 0.015

20 0.025

40 0.02

40 0.03

40 0.03];

save('process.mat');

% 3.2MC抽样机组确定机组状态

% 3.2.1计算含风电场的RBTS可靠性

% 共有7类机组，常规机组状态在StateNorm【出力概率】元胞数组中，风电状态在StateFORWeibull6【出力概率】

I=0 %I用来记录发生却负荷的次数

sumDNS=0;

DNS=zeros(200000,1);

K=rand(200000,12);%1-11常规12风电

pwind=zeros(200000,1);

for k=1:200000

Pout=zeros(12,1);

%得到一次抽样常规机组状态

fori=1:11

if K(k,i)>Generator.Norm(i,2)

Pout(i)=Generator.Norm(i,1);

else

Pout(i)=0;

end

end

%得到一次抽样风电场的出力

t=0;

fori=1:M

t=t+StateFORWeibull6(i,2);

if K(k,12)-t>=0

continue

else

Pout(12)=StateFORWeibull6(i,1);

pwind(k)= Pout(12);

break

end

end

sumPout=sum(Pout);%发电系统出力

% 统计系统状态参数

% 持续峰荷185MW

RBTSload=185;

DNS(k)=max([0;RBTSload-sumPout]);

if DNS(k)>0

I=I+1; %I用来记录发生却负荷的次数sumDNS=sumDNS+DNS(k);

end

end

LOLP1=I/200000;

EENS1=sumDNS/200000*8760;

LOLE1=I/200000*8760;

disp('LOLP1=');disp(LOLP1);

disp('EENS1=');disp(EENS1);

disp('LOLE1=');disp(LOLE1);

save('resultwind.mat','LOLP1','EENS1','LOLE1'); %% 3.2.2计算常规RBTS系统可靠性

clc

clear

load process

I=0;

sumDNS=0;

DNS=zeros(100000,1);

K=rand(100000,11);%1-11常规7风电

for k=1:100000

Pout=zeros(11,1);

%得到一次抽样常规机组状态

fori=1:11

if K(k,i)>Generator.Norm(i,2)

Pout(i)=Generator.Norm(i,1);

else

Pout(i)=0;

end

end

sumPout=sum(Pout);%发电系统出力

% 统计系统状态参数

% 持续峰荷185MW

RBTSload=185;

DNS(k)=max([0;RBTSload-sumPout]);

if DNS(k)>0

I=I+1; %I用来记录发生却负荷的次数sumDNS=sumDNS+DNS(k);

end

end

% 计算LOLP EENS LOLE

LOLP2=I/100000;

EENS2=sumDNS/100000*8760;

LOLE2=I/100000*8760;

disp('LOLP2=');disp(LOLP2);

disp('EENS2=');disp(EENS2);

disp('LOLE2=');disp(LOLE2);

save('result0.mat','LOLP2','EENS2','LOLE2');