????????????????????????????????????????下載資源鏈接??????????????????????????????????????????文章來(lái)源地址http://www.zghlxwxcb.cn/news/detail-778728.html
摘要:分布式電源的接入使得配電系統(tǒng)從放射狀無(wú)源網(wǎng)絡(luò)變?yōu)榉植加兄行⌒碗娫吹挠性淳W(wǎng)絡(luò)。帶來(lái)了使單向流動(dòng)的電流方向具有了不確定性等等問(wèn)題,使得配電系統(tǒng)的控制和管理變得更加復(fù)雜。但同時(shí),分布式電源又具有提高電網(wǎng)可靠性,綠色節(jié)能,等等優(yōu)點(diǎn),所以為更好的利用分布式電源為人類(lèi)造福,我們必須對(duì)其進(jìn)行研究與分析。
本文采取通過(guò)利用仿真軟件Matlab編寫(xiě)計(jì)算潮流程序模擬分布式電源接入配電網(wǎng)的模型進(jìn)行潮流計(jì)算的方法對(duì)分布式電源的穩(wěn)態(tài)影響進(jìn)行探索與分析。選取了34節(jié)點(diǎn)的配電網(wǎng)網(wǎng)絡(luò)模型,通過(guò)對(duì)單個(gè)以及多個(gè)分布式電源的接入位置以及容量的不同情況對(duì)34節(jié)點(diǎn)配電網(wǎng)的網(wǎng)損以及節(jié)點(diǎn)電壓狀況進(jìn)行了分析。
部分代碼展示:
clc
clear
close all
%%運(yùn)行潮流
t0 = clock;
%% 列索引
[PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
[F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, RATE_C, TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST, ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idx_brch;
[GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, GEN_STATUS, PMAX, PMIN, MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN, PC1, PC2, QC1MIN, QC1MAX, QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = idx_gen;
casename = 'case9';
mpopt = mpoption; %% use default options
%% 打印選項(xiàng)
verbose = mpopt(31);
%% 1.初始潮流
mpc = loadcase(casename);
mpc.branch
[baseMVA, bus, gen, branch] = deal(mpc.baseMVA, mpc.bus, mpc.gen, mpc.branch);
[ref, pv, pq] = bustypes(bus, gen);
[Ybus, Yf, Yt] = makeYbus(baseMVA, bus, branch);
Sbus = makeSbus(baseMVA, bus, gen);
V0 = ones(size(bus, 1), 1);
[V, success, iterations] = newtonpf(Ybus, Sbus, V0, ref, pv, pq, mpopt);
V1=sqrt(abs(V).^2);
figure
plot(V1,'-*')
[bus, gen, branch] = pfsoln(baseMVA, bus, gen, branch, Ybus, Yf, Yt, V, ref, pv, pq);
mpc.et = etime(clock, t0);
mpc.success = success;
%% ----- 輸出結(jié)果 -----
[mpc.bus, mpc.gen, mpc.branch] = deal(bus, gen, branch);
printpf(mpc, 1, mpopt);
%% 2.加入分布式電源后潮流
% 讀數(shù)據(jù)
fprintf('請(qǐng)輸入分布式電源的節(jié)點(diǎn)位置,有功與無(wú)功大小 bus Pg Qg ')
DGs2=[0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0]
DGbus=input('請(qǐng)輸入DG節(jié)點(diǎn)位置的值:');%4
DGPg=input('請(qǐng)輸入DG有功大小的值:');%30
DGQg=input('請(qǐng)輸入DG無(wú)功大小的值:');%30
DGs1=[DGbus DGPg DGQg]
DGs=[DGs1,DGs2]
mpc.gen=[mpc.gen;DGs]
%% 轉(zhuǎn)換為內(nèi)部輸入
[baseMVA, bus, gen, branch] = deal(mpc.baseMVA, mpc.bus, mpc.gen, mpc.branch);
%% 轉(zhuǎn)換為pq,pv,參考節(jié)點(diǎn)
[ref, pv, pq] = bustypes(bus, gen);
%%電源信息
on = find(gen(:, GEN_STATUS) > 0); %% 哪個(gè)電源在工作
gbus = gen(on, GEN_BUS); %% 在那條母線(xiàn)
%% 初始化
V0 = ones(size(bus, 1), 1); %% flat start
% V0 = bus(:, VM) .* exp(sqrt(-1) * pi/180 * bus(:, VA));
V0(gbus) = gen(on, VG) ./ abs(V0(gbus)).* V0(gbus);
%% 建立導(dǎo)納矩陣
[Ybus, Yf, Yt] = makeYbus(baseMVA, bus, branch);
%% 計(jì)算母線(xiàn)功率注入(發(fā)電 - 負(fù)載)
Sbus = makeSbus(baseMVA, bus, gen);
%% 潮流計(jì)算
[V, success, iterations] = newtonpf(Ybus, Sbus, V0, ref, pv, pq, mpopt);
V2=sqrt(abs(V).^2);
效果展示:
資源鏈接地址https://download.csdn.net/download/LIANG674027206/88691937文章來(lái)源:http://www.zghlxwxcb.cn/news/detail-778728.html
????????????????????????????????????????下載資源鏈接??????????????????????????????????????????
到了這里,關(guān)于008文章解讀與程序——《含分布式電源接入的配電網(wǎng)可靠性評(píng)估》已提供下載資源的文章就介紹完了。如果您還想了解更多內(nèi)容,請(qǐng)?jiān)谟疑辖撬阉鱐OY模板網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章,希望大家以后多多支持TOY模板網(wǎng)!