53_PowerFlowN-1/SecurityCertification.m at master · myncepu/53_PowerFlowN-1 · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
function [branchOff, flagVoltageViolation, voltageViolation, flagPowerViolation, powerViolation] = ...
    SecurityCertification(filename, iBranch)
% 蜿よ焚蜷ｫ荵
% 霎灘�蜿よ焚
% filename		逕ｵ鄂第焚謐ｮ譁�ｻｶ
% iBranch		譁ｭ蠑謾ｯ霍ｯ郛門捷
% 霎灘�蜿よ焚
% branchOff 		譁ｭ蠑謾ｯ霍ｯ襍ｷ蟋玖鰍轤ｹ郛門捷
% flagVoltageViolation	蛻､譁ｭ譁ｭ蠑謾ｯ霍ｯiBranch蜷守鳩蜴区弍蜷ｦ雜企剞�悟ｦよ棡雜企剞�悟�荳?�悟凄蛻吩ｸｺ0
% voltageViolation	隶ｰ蠖戊ｶ企剞闃らせ逕ｵ蜴
% flagPowerViolation	蛻､譁ｭ譁ｭ蠑謾ｯ霍ｯiBranch蜷守ｺｿ霍ｯ蜉溽紫譏ｯ蜷ｦ雜企剞�悟ｦよ棡雜企剞�悟�荳?�悟凄蛻吩ｸｺ0
% powerViolation	隶ｰ蠖戊ｶ企剞郤ｿ霍ｯ蜉溽紫諠���梧焚謐ｮ譬ｼ蠑丈ｸｺ
% PowerViolation =
%	蜷��蜷ｫ荵我ｸｺ襍ｷ蟋玖鰍轤ｹ?譛ｫ遶ｯ闃らせ縲∫ｺｿ霍ｯ鬥也ｫｯ豬∝髄譛ｫ遶ｯ隗�惠蜉溽紫?郤ｿ霍ｯ譛ｫ遶ｯ豬∝髄鬥也ｫｯ隗�惠蜉溽紫縲?
%	郤ｿ霍ｯ鬥也ｫｯ豬∝髄譛ｫ遶ｯ譛牙粥蜉溽紫縲∫ｺｿ霍ｯ鬥也ｫｯ豬∝髄譛ｫ遶ｯ譌蜉溷粥邇�?
%	郤ｿ霍ｯ譛ｫ遶ｯ豬∝髄鬥也ｫｯ譛牙粥蜉溽紫縲∫ｺｿ霍ｯ譛ｫ遶ｯ豬∝髄鬥也ｫｯ譌蜉溷粥邇?
%	[4.0000	9.0000	16.5119	16.4341	16.4326	1.6161	-16.4326	-0.2195
%	4.0000	9.0000	31.3803	30.6433	30.5782	7.0496	-30.5782	-1.9961
%	逶ｸ驍ｻ荳､蛻玲髪霍ｯ郛門捷�域ｯ剰｡悟燕荳､蛻暦ｼ臥嶌蜷鯉ｼ悟�蛻ｫ荳ｺ譁ｭ蠑?髪霍ｯiBranch蜑崎ｶ企剞謾ｯ霍ｯ蜉溽紫迥ｶ蜀ｵ�悟?
%	譁ｭ蠑iBranch蜷取髪霍ｯ蜉溽紫迥ｶ蜀?
%	6.0000	11.0000	7.7175	7.6241	6.9120	3.4327	-6.8607	-3.3253
%	6.0000	11.0000	16.4393	16.0894	15.8367	4.4102	-15.6039	-3.9227]

[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;
% IEEE14 蛹�性20荳ｪ謾ｯ霍ｯ�悟｢槫刈荳?ｸｪ謾ｯ霍ｯ
mpopt = mpoption('verbose', 0, 'out.all', 0); % 貎ｮ豬∬ｮ｡邂嶺ｸ崎ｾ灘�扈捺?
% mpopt = mpoption;
rOrigin = runpf(filename, mpopt); % 蜴滓擂逕ｵ鄂第ｽｮ豬∬ｮ｡邂
rOrigin.branch(iBranch, :) = [];

MPC = loadcase(filename);
branchOff = MPC.branch(iBranch, [F_BUS, T_BUS]);
MPC.branch(iBranch, :) = []; % 譁ｭ蠑謾ｯ霍ｯiBranch

% savecase('my_case14_SC.m', MPC);
rNow = runpf(MPC, mpopt); % 譁ｰ逕ｵ鄂第ｽｮ豬∬ｮ｡邂?
flagVoltageViolation = any(abs(rNow.bus(:, VM) - 1) > 0.5); % 蛻､譁ｭ譁ｭ蠑謾ｯ霍ｯ蜷守鳩鄂醍鳩蜴区弍蜷ｦ雜企?
voltageViolation = [];
for k = 1:size(MPC.bus)
    if abs(rNow.bus(k, VM) - 1) > 0.5
        voltageViolation = [voltageViolation; rNow.bus(k, [BUS_I, VM])];
    end
end
branchPowerOrigin = rOrigin.branch(:, [PF, QF, PT, QT]);
% branchPowerOrigin(iBranch, :) = [];
branchApparentPowerOrigin = [sqrt(branchPowerOrigin(:, 1) .^ 2 + branchPowerOrigin(:, 2) .^ 2) ...
    sqrt(branchPowerOrigin(:, 3) .^ 2 + branchPowerOrigin(:, 4) .^ 2)];
branchPowerNow = rNow.branch(:, [PF, QF, PT, QT]);
branchApparentPowerNow = [sqrt(branchPowerNow(:, 1) .^ 2 + branchPowerNow(:, 2) .^ 2) ...
    sqrt(branchPowerNow(:, 3) .^ 2 + branchPowerNow(:, 4) .^ 2)];
powerViolationFlag = branchApparentPowerNow > 1.5 * branchApparentPowerOrigin;
flagPowerViolation = any(any(powerViolationFlag));
powerViolationBranch = [rNow.branch(:, [F_BUS, T_BUS]), powerViolationFlag];
powerViolation = [];
for j = 1:size(MPC.branch, 1)
    if any(powerViolationBranch(j, 3:4) ~= 0)
%         powerViolation = [powerViolation; powerViolationBranch(j, :), ...
%             branchApparentPowerOrigin(j, :), branchApparentPowerNow(j, :)];
        powerViolation = [powerViolation; rOrigin.branch(j, [F_BUS, T_BUS]), ...
            branchApparentPowerOrigin(j, :), rOrigin.branch(j, [PF, QF, PT, QT]); ...
            rNow.branch(j, [F_BUS, T_BUS]), branchApparentPowerNow(j, :), ...
            rNow.branch(j, [PF, QF, PT, QT])];
    end
end