This is the script I'm using to split it up into independent areas. Is it
the right approach?
% Case 2 area 1
close all
clear all
OPT = mpoption('OUT_ALL', 0);
mpc = loadcase('case39');
mpc1=mpc;
mpc1.bus=mpc.bus([1:5,18,25,30,37,39],:);
mpc1.bus(10,2)=3;
mpc1.gen=mpc.gen([1,8,10],:);
mpc1.branch=mpc.branch([1:4,6:8,41,5],:);
for i=1:length(mpc1.bus(:,1))
mpc1.bus(i:7,8)=1;
mpc1.bus(i,9)=0;
end
mpc1.bus(1,3)=0;mpc1.bus(1,4)=0;mpc1.bus(2,3)=0;mpc1.bus(2,4)=0;mpc1.bus(3,3)=322;mpc1.bus(3,4)=2.4;mpc1.bus(4,3)=229.027;mpc1.bus(4,4)=132.433;
mpc1.bus(5,3)=-135.472;mpc1.bus(5,4)=46.542;mpc1.bus(6,3)=-53.531;mpc1.bus(6,4)=46.542;mpc1.bus(7,3)=296.4;mpc1.bus(7,4)=-87.947;
mpc1.bus(8,3)=0;mpc1.bus(8,4)=0;mpc1.bus(9,3)=0;mpc1.bus(9,4)=0;mpc1.bus(10,3)=980.795;mpc1.bus(10,4)=15.692;
for i=1:length(mpc1.branch(:,1))
if mpc1.branch(i,9)==0
mpc1.branch(i,9)=1;
end
end
for i=1:length(mpc1.gen(:,1))
mpc1.gen(i,3)=0;
mpc1.gen(i,4)=300;
mpc1.gen(i,5)=-300;
end
Thanks
On Sun, Dec 20, 2015 at 2:34 PM, Abhyankar, Shrirang G. <abhy...@anl.gov>
wrote:
> Yes. However, note that each area needs to have enough generation capacity
> to supply the load.
>
> Shri
>
> > On Dec 19, 2015, at 12:42 PM, "Owusuboateng, Prince Charles" <
> powusuboateng8...@mytu.tuskegee.edu> wrote:
> >
> > Hello,
> >
> > I have the 39 bus system divided into 4 independent areas. Can I run
> individual optimal power flow analysis on each of these areas independently?
> >
> > thank you
>
>
>
