Re: Zonal reserve with stochastic Wind generation in MOST

[Ray Zimmerman]

The zonal reserves in MOST are designed for deterministic problems. See 
footnote 19 at the bottom of page 35 in the MOST User’s Manual 
.

So, it is designed to handle either zonal reserve requirements or stochastic 
wind, but not both simultaneously. A clear understanding of MOST's problem 
formulation as described in the manual should help here.

   Ray


> On Sep 21, 2017, at 8:53 PM, Quarm JNR, Edward A 
>  wrote:
> 
> 
> Dear All,
> 
> I am running a multi-period simulation with :
> 
> Zonal reserve requirement
> Stochastic Wind generator
> Load profile 
> 
> My output for reserve generation always includes an extra row for the wind 
> generator but with zero output. For example, I have 5 generators in my system 
> and I add a Wind generator making it 6 generators. The output for reserve 
> generation will have 6 rows with the last row having zero output.
> 
> I would like to know, how does MOST consider the wind generator in running 
> unit commitment simulations with regards to zonal reserve.
> Is the wind generator constrained to have zero output or is it totally 
> ignored when optimizing the multi period reserve generation?
> What if the Wind generator belongs to a reserve zone?
> If the wind generators don't play any role in zonal reserve why are they 
> included in the output?
> 
> Thanks in advance for your clarification. 
> 
> Best regards,
> Edward

Re: Adding Wind Generation to Profiles

[Ray Zimmerman]

Josh,

Looks like you got it working by splitting out the individual wind profiles. It 
should also work to set the rows field to [1:7] and then define a matrix of 
profiles for each of the 7 generators by assigning an nt x 7 matrix to 
windprofile(:, 1, :).

Hope this helps,

Ray




> On Sep 23, 2017, at 9:20 PM, Joshua Sebben <joshkeep...@gmail.com> wrote:
> 
> Stephen,
>  
> Thank you so much! That works a treat!
>  
> Regards,
>  
> Josh
>  
> From: Stephen Suffian <mailto:ssuff...@villanova.edu>
> Sent: Saturday, 23 September 2017 3:03 AM
> To: MATPOWER discussion forum <mailto:matpowe...@list.cornell.edu>
> Subject: Re: Adding Wind Generation to Profiles
>  
> Check out the 'help get_profile'. Much like you pass the previous profiles 
> object on the load profile line, you can pass the previous profiles object on 
> the other lines:
> 
> profiles = getprofiles('wind_profile1',iwind);
> profiles = getprofiles('wind_profile2',profiles,iwind);
> profiles = getprofiles('load_profile_30',profiles)
> 
> See if that works!
>  
> On Thu, Sep 21, 2017 at 7:24 PM, Joshua Sebben <joshkeep...@gmail.com 
> <mailto:joshkeep...@gmail.com>> wrote:
> Stephen,
>  
> Thank you for your reply.
>  
> casefile = 'case30_modified_MedPen';
> mpc = loadcase(casefile);
> xgd = loadxgendata('xgd_uc_30_MedPen', mpc);
> [iwind, mpc, xgd] = addwind('wind_uc_30_MedPen', mpc, xgd);
> profiles = getprofiles('wind_profile1', iwind);
> profiles = getprofiles('wind_profile2', iwind);
> profiles = getprofiles('wind_profile3', iwind);
> profiles = getprofiles('wind_profile4', iwind);
> profiles = getprofiles('wind_profile5', iwind);
> profiles = getprofiles('wind_profile6', iwind);
> profiles = getprofiles('wind_profile7', iwind);
> profiles = getprofiles('load_profile_30', profiles);
> nt = size(profiles(1).values, 1);   % number of periods
>  
> I created the 7 wind profiles with each of the row counts changed to 1 to 7 
> for their respective generators.
>  
> This would then change my main code assuming to the snippet above? 
> This again also has the problem of only changing the last generator in the 
> list as each of the previous wind profile changes are overwritten by the new 
> getprofile?
>  
> Regards,
> Josh
>  
> On 22 September 2017 at 04:02, Stephen Suffian <ssuff...@villanova.edu 
> <mailto:ssuff...@villanova.edu>> wrote:
> 
> I also couldn't figure it out, but I got around it by adding a separate 
> profile for each generator and changing the rows value. So if a wind 
> generator is in row 5, you would have the struct look like this below.
> 
>  
> windprofile = struct( ...
> 'type', 'mpcData', ...
> 'table', CT_TGEN, ...
> 'rows', 5, ...
> 'col', PMAX, ...
> 'chgtype', CT_REL, ...
> 'values', [] );
>  
> And then for row 6 you would add another profile with the struct looking 
> tlike this:
>  
>  
> windprofile = struct( ...
> 'type', 'mpcData', ...
> 'table', CT_TGEN, ...
> 'rows', 6, ...
> 'col', PMAX, ...
> 'chgtype', CT_REL, ...
> 'values', [] );
>  
> It is a bit of a work around, but I believe it should work (I did something 
> similar to set the Pmax for each conventional generator in a profile, so I 
> imagine it will work the same).
>  
> On Thu, Sep 21, 2017 at 1:52 PM, Stephen Suffian <stephen.suff...@gmail.com 
> <mailto:stephen.suff...@gmail.com>> wrote:
> I also couldn't figure it out, but I got around it by adding a separate 
> profile for each generator and changing the rows value. So if a wind 
> generator is in row 5, you would have the struct look like this below.
> 
>  
> windprofile = struct( ...
> 'type', 'mpcData', ...
> 'table', CT_TGEN, ...
> 'rows', 5, ...
> 'col', PMAX, ...
> 'chgtype', CT_REL, ...
> 'values', [] );
>  
> And then for row 6 you would add another profile with the struct looking 
> tlike this:
>  
>  
> windprofile = struct( ...
> 'type', 'mpcData', ...
> 'table', CT_TGEN, ...
> 'rows', 6, ...
> 'col', PMAX, ...
> 'chgtype', CT_REL, ...
> 'values', [] );
>  
> It is a bit of a work around, but I believe it should work (I did something 
> similar to set the Pmax for each conventional generator in a profile, so I 
> imagine it will work the same).
>  
> On Thu, Sep 21, 2017 at 5:39 AM, Joshua Sebben <joshkeep...@gmail.com 
> <mailto:joshkeep...@gmail.com>> wrote:
> By the way,
>  
> I have also tried setting the row count to [1 2 3 4 5 6 7] for my 7 extra 
> generators that I want to add, however when I run my code I get an error:
>

RE: Adding Wind Generation to Profiles

[Joshua Sebben]

Stephen,

Thank you so much! That works a treat! 

Regards,

Josh

From: Stephen Suffian
Sent: Saturday, 23 September 2017 3:03 AM
To: MATPOWER discussion forum
Subject: Re: Adding Wind Generation to Profiles

Check out the 'help get_profile'. Much like you pass the previous profiles 
object on the load profile line, you can pass the previous profiles object on 
the other lines:
profiles = getprofiles('wind_profile1',iwind);
profiles = getprofiles('wind_profile2',profiles,iwind);
profiles = getprofiles('load_profile_30',profiles)
See if that works!

On Thu, Sep 21, 2017 at 7:24 PM, Joshua Sebben <joshkeep...@gmail.com> wrote:
Stephen, 

Thank you for your reply.

casefile = 'case30_modified_MedPen';
mpc = loadcase(casefile);
xgd = loadxgendata('xgd_uc_30_MedPen', mpc);
[iwind, mpc, xgd] = addwind('wind_uc_30_MedPen', mpc, xgd);
profiles = getprofiles('wind_profile1', iwind);
profiles = getprofiles('wind_profile2', iwind);
profiles = getprofiles('wind_profile3', iwind);
profiles = getprofiles('wind_profile4', iwind);
profiles = getprofiles('wind_profile5', iwind);
profiles = getprofiles('wind_profile6', iwind);
profiles = getprofiles('wind_profile7', iwind);
profiles = getprofiles('load_profile_30', profiles);
nt = size(profiles(1).values, 1);       % number of periods

I created the 7 wind profiles with each of the row counts changed to 1 to 7 for 
their respective generators.

This would then change my main code assuming to the snippet above? 
This again also has the problem of only changing the last generator in the list 
as each of the previous wind profile changes are overwritten by the new 
getprofile?

Regards,
Josh

On 22 September 2017 at 04:02, Stephen Suffian <ssuff...@villanova.edu> wrote:

I also couldn't figure it out, but I got around it by adding a separate profile 
for each generator and changing the rows value. So if a wind generator is in 
row 5, you would have the struct look like this below.

windprofile = struct( ...
    'type', 'mpcData', ...
    'table', CT_TGEN, ...
    'rows', 5, ...
    'col', PMAX, ...
    'chgtype', CT_REL, ...
    'values', [] );

And then for row 6 you would add another profile with the struct looking tlike 
this:


windprofile = struct( ...
    'type', 'mpcData', ...
    'table', CT_TGEN, ...
    'rows', 6, ...
    'col', PMAX, ...
    'chgtype', CT_REL, ...
    'values', [] );

It is a bit of a work around, but I believe it should work (I did something 
similar to set the Pmax for each conventional generator in a profile, so I 
imagine it will work the same).

On Thu, Sep 21, 2017 at 1:52 PM, Stephen Suffian <stephen.suff...@gmail.com> 
wrote:
I also couldn't figure it out, but I got around it by adding a separate profile 
for each generator and changing the rows value. So if a wind generator is in 
row 5, you would have the struct look like this below.

windprofile = struct( ...
    'type', 'mpcData', ...
    'table', CT_TGEN, ...
    'rows', 5, ...
    'col', PMAX, ...
    'chgtype', CT_REL, ...
    'values', [] );

And then for row 6 you would add another profile with the struct looking tlike 
this:


windprofile = struct( ...
    'type', 'mpcData', ...
    'table', CT_TGEN, ...
    'rows', 6, ...
    'col', PMAX, ...
    'chgtype', CT_REL, ...
    'values', [] );

It is a bit of a work around, but I believe it should work (I did something 
similar to set the Pmax for each conventional generator in a profile, so I 
imagine it will work the same).

On Thu, Sep 21, 2017 at 5:39 AM, Joshua Sebben <joshkeep...@gmail.com> wrote:
By the way,
 
I have also tried setting the row count to [1 2 3 4 5 6 7] for my 7 extra 
generators that I want to add, however when I run my code I get an error:
 
Error using apply_profile (line 148)
apply_profile: third dimension of profile.values should match length of pro=
file.rows
 
Error in loadmd (line 508)
    optab =3D apply_profile(profiles(p), optab);
 
Error in Test (line 33)
mdi =3D loadmd(mpc, transmat, xgd, [], [], profiles);

Regards,
Josh

On 20 September 2017 at 22:17, Joshua Sebben <joshkeep...@gmail.com> wrote:
Currently working on adding wind generator units to my model. I am working off 
the 30 bus example case in MOST.  I am trying to run the following code 


casefile = 'case30';
mpc = loadcase(casefile);
xgd = loadxgendata('xgd_uc', mpc);
[iwind, mpc, xgd] = addwind('wind_uc_30', mpc, xgd);
profiles = getprofiles('wind_profile', iwind);
profiles = getprofiles('load_profile', profiles);
nt = size(profiles(1).values, 1);       % number of periods

%%-  Full Transition Probabilities  -
transmat = ex_transmat(nt);
mdi = loadmd(mpc, transmat, xgd, [], [], profiles);
mdo = most(mdi, mpopt);
if verbose
    ms = most_summary(mdo);
end

However the added generator units of which there are 7 of them in wind_uc_30 
don't seem to get the profile set to them in the output. Instead only the first 
Generator unit gets the profile while the rest of them are set at PMAX.

windp

Re: Adding Wind Generation to Profiles

[Stephen Suffian]

Check out the 'help get_profile'. Much like you pass the previous profiles
object on the load profile line, you can pass the previous profiles object
on the other lines:

profiles = getprofiles('wind_profile1',iwind);
profiles = getprofiles('wind_profile2',profiles,iwind);
profiles = getprofiles('load_profile_30',profiles)

See if that works!

On Thu, Sep 21, 2017 at 7:24 PM, Joshua Sebben 
wrote:

> Stephen,
>
> Thank you for your reply.
>
> casefile = 'case30_modified_MedPen';
> mpc = loadcase(casefile);
> xgd = loadxgendata('xgd_uc_30_MedPen', mpc);
> [iwind, mpc, xgd] = addwind('wind_uc_30_MedPen', mpc, xgd);
> profiles = getprofiles('wind_profile1', iwind);
> profiles = getprofiles('wind_profile2', iwind);
> profiles = getprofiles('wind_profile3', iwind);
> profiles = getprofiles('wind_profile4', iwind);
> profiles = getprofiles('wind_profile5', iwind);
> profiles = getprofiles('wind_profile6', iwind);
> profiles = getprofiles('wind_profile7', iwind);
> profiles = getprofiles('load_profile_30', profiles);
> nt = size(profiles(1).values, 1);   % number of periods
>
> I created the 7 wind profiles with each of the row counts changed to 1 to
> 7 for their respective generators.
>
> This would then change my main code assuming to the snippet above?
> This again also has the problem of only changing the last generator in the
> list as each of the previous wind profile changes are overwritten by the
> new getprofile?
>
> Regards,
> Josh
>
> On 22 September 2017 at 04:02, Stephen Suffian 
> wrote:
>
>> I also couldn't figure it out, but I got around it by adding a separate
>> profile for each generator and changing the rows value. So if a wind
>> generator is in row 5, you would have the struct look like this below.
>>
>>
>> windprofile = struct( ...
>> 'type', 'mpcData', ...
>> 'table', CT_TGEN, ...
>> 'rows', 5, ...
>> 'col', PMAX, ...
>> 'chgtype', CT_REL, ...
>> 'values', [] );
>>
>> And then for row 6 you would add another profile with the struct looking
>> tlike this:
>>
>>
>> windprofile = struct( ...
>> 'type', 'mpcData', ...
>> 'table', CT_TGEN, ...
>> 'rows', 6, ...
>> 'col', PMAX, ...
>> 'chgtype', CT_REL, ...
>> 'values', [] );
>>
>> It is a bit of a work around, but I believe it should work (I did
>> something similar to set the Pmax for each conventional generator in a
>> profile, so I imagine it will work the same).
>>
>> On Thu, Sep 21, 2017 at 1:52 PM, Stephen Suffian <
>> stephen.suff...@gmail.com> wrote:
>>
>>> I also couldn't figure it out, but I got around it by adding a separate
>>> profile for each generator and changing the rows value. So if a wind
>>> generator is in row 5, you would have the struct look like this below.
>>>
>>>
>>> windprofile = struct( ...
>>> 'type', 'mpcData', ...
>>> 'table', CT_TGEN, ...
>>> 'rows', 5, ...
>>> 'col', PMAX, ...
>>> 'chgtype', CT_REL, ...
>>> 'values', [] );
>>>
>>> And then for row 6 you would add another profile with the struct looking
>>> tlike this:
>>>
>>>
>>> windprofile = struct( ...
>>> 'type', 'mpcData', ...
>>> 'table', CT_TGEN, ...
>>> 'rows', 6, ...
>>> 'col', PMAX, ...
>>> 'chgtype', CT_REL, ...
>>> 'values', [] );
>>>
>>> It is a bit of a work around, but I believe it should work (I did
>>> something similar to set the Pmax for each conventional generator in a
>>> profile, so I imagine it will work the same).
>>>
>>> On Thu, Sep 21, 2017 at 5:39 AM, Joshua Sebben 
>>> wrote:
>>>
 By the way,



 I have also tried setting the row count to [1 2 3 4 5 6 7] for my 7
 extra generators that I want to add, however when I run my code I get an
 error:



 Error using apply_profile (line 148)

 apply_profile: third dimension of profile.values should match length of
 pro=

 file.rows



 Error in loadmd (line 508)

 optab =3D apply_profile(profiles(p), optab);



 Error in Test (line 33)

 mdi =3D loadmd(mpc, transmat, xgd, [], [], profiles);


 Regards,

 Josh

 On 20 September 2017 at 22:17, Joshua Sebben 
 wrote:

> Currently working on adding wind generator units to my model. I am
> working off the 30 bus example case in MOST.  I am trying to run the
> following code
>
> casefile = 'case30';
> mpc = loadcase(casefile);
> xgd = loadxgendata('xgd_uc', mpc);
> [iwind, mpc, xgd] = addwind('wind_uc_30', mpc, xgd);
> profiles = getprofiles('wind_profile', iwind);
> profiles = getprofiles('load_profile', profiles);
> nt = size(profiles(1).values, 1);   % number of periods
>
> %%-  Full Transition Probabilities  -
> transmat = ex_transmat(nt);
> mdi = loadmd(mpc, transmat, xgd, [], [], profiles);
> mdo = most(mdi, mpopt);
> 

Zonal reserve with stochastic Wind generation in MOST

[Quarm JNR, Edward A]

Dear All,


I am running a multi-period simulation with :


  1.  Zonal reserve requirement
  2.  Stochastic Wind generator
  3.  Load profile


My output for reserve generation always includes an extra row for the wind 
generator but with zero output. For example, I have 5 generators in my system 
and I add a Wind generator making it 6 generators. The output for reserve 
generation will have 6 rows with the last row having zero output.


  1.  I would like to know, how does MOST consider the wind generator in 
running unit commitment simulations with regards to zonal reserve.
  2.  Is the wind generator constrained to have zero output or is it totally 
ignored when optimizing the multi period reserve generation?
  3.  What if the Wind generator belongs to a reserve zone?
  4.  If the wind generators don't play any role in zonal reserve why are they 
included in the output?

Thanks in advance for your clarification.

Best regards,
Edward

Re: Adding Wind Generation to Profiles

[Stephen Suffian]

I also couldn't figure it out, but I got around it by adding a separate
profile for each generator and changing the rows value. So if a wind
generator is in row 5, you would have the struct look like this below.


windprofile = struct( ...
'type', 'mpcData', ...
'table', CT_TGEN, ...
'rows', 5, ...
'col', PMAX, ...
'chgtype', CT_REL, ...
'values', [] );

And then for row 6 you would add another profile with the struct looking
tlike this:


windprofile = struct( ...
'type', 'mpcData', ...
'table', CT_TGEN, ...
'rows', 6, ...
'col', PMAX, ...
'chgtype', CT_REL, ...
'values', [] );

It is a bit of a work around, but I believe it should work (I did something
similar to set the Pmax for each conventional generator in a profile, so I
imagine it will work the same).

On Thu, Sep 21, 2017 at 1:52 PM, Stephen Suffian 
wrote:

> I also couldn't figure it out, but I got around it by adding a separate
> profile for each generator and changing the rows value. So if a wind
> generator is in row 5, you would have the struct look like this below.
>
>
> windprofile = struct( ...
> 'type', 'mpcData', ...
> 'table', CT_TGEN, ...
> 'rows', 5, ...
> 'col', PMAX, ...
> 'chgtype', CT_REL, ...
> 'values', [] );
>
> And then for row 6 you would add another profile with the struct looking
> tlike this:
>
>
> windprofile = struct( ...
> 'type', 'mpcData', ...
> 'table', CT_TGEN, ...
> 'rows', 6, ...
> 'col', PMAX, ...
> 'chgtype', CT_REL, ...
> 'values', [] );
>
> It is a bit of a work around, but I believe it should work (I did
> something similar to set the Pmax for each conventional generator in a
> profile, so I imagine it will work the same).
>
> On Thu, Sep 21, 2017 at 5:39 AM, Joshua Sebben 
> wrote:
>
>> By the way,
>>
>>
>>
>> I have also tried setting the row count to [1 2 3 4 5 6 7] for my 7 extra
>> generators that I want to add, however when I run my code I get an error:
>>
>>
>>
>> Error using apply_profile (line 148)
>>
>> apply_profile: third dimension of profile.values should match length of
>> pro=
>>
>> file.rows
>>
>>
>>
>> Error in loadmd (line 508)
>>
>> optab =3D apply_profile(profiles(p), optab);
>>
>>
>>
>> Error in Test (line 33)
>>
>> mdi =3D loadmd(mpc, transmat, xgd, [], [], profiles);
>>
>>
>> Regards,
>>
>> Josh
>>
>> On 20 September 2017 at 22:17, Joshua Sebben 
>> wrote:
>>
>>> Currently working on adding wind generator units to my model. I am
>>> working off the 30 bus example case in MOST.  I am trying to run the
>>> following code
>>>
>>> casefile = 'case30';
>>> mpc = loadcase(casefile);
>>> xgd = loadxgendata('xgd_uc', mpc);
>>> [iwind, mpc, xgd] = addwind('wind_uc_30', mpc, xgd);
>>> profiles = getprofiles('wind_profile', iwind);
>>> profiles = getprofiles('load_profile', profiles);
>>> nt = size(profiles(1).values, 1);   % number of periods
>>>
>>> %%-  Full Transition Probabilities  -
>>> transmat = ex_transmat(nt);
>>> mdi = loadmd(mpc, transmat, xgd, [], [], profiles);
>>> mdo = most(mdi, mpopt);
>>> if verbose
>>> ms = most_summary(mdo);
>>> end
>>>
>>> However the added generator units of which there are 7 of them in
>>> wind_uc_30 don't seem to get the profile set to them in the output. Instead
>>> only the first Generator unit gets the profile while the rest of them are
>>> set at PMAX.
>>>
>>> windprofile = struct( ...
>>> 'type', 'mpcData', ...
>>> 'table', CT_TGEN, ...
>>> 'rows', 1, ...
>>> 'col', PMAX, ...
>>> 'chgtype', CT_REL, ...
>>> 'values', [] );
>>>
>>>  windprofile.values(:, :, 1) = [
>>>   0.80;
>>>   0.65;
>>>   0.60;
>>>   0.82;
>>>   1.00;
>>>   0.70;
>>>   0.50;
>>>   0.85;
>>>   1.00;
>>>   1.10;
>>>   1.06;
>>>   0.95;
>>>
>>> Above is a snippet from the wind_profile.  I am assuming it has
>>> something to do with the row count, however, I haven't been able to find a
>>> solution.
>>>
>>> Could I please get some help. to ensure all the generator units receive
>>> the load profile.
>>>
>>> Thankyou,
>>>
>>> --
>>> J.Sebben
>>>
>>
>>
>>
>> --
>> J.Sebben
>>
>
>

Adding Wind Generation to Profiles

[Joshua Sebben]

Currently working on adding wind generator units to my model. I am working
off the 30 bus example case in MOST.  I am trying to run the following code

casefile = 'case30';
mpc = loadcase(casefile);
xgd = loadxgendata('xgd_uc', mpc);
[iwind, mpc, xgd] = addwind('wind_uc_30', mpc, xgd);
profiles = getprofiles('wind_profile', iwind);
profiles = getprofiles('load_profile', profiles);
nt = size(profiles(1).values, 1);   % number of periods

%%-  Full Transition Probabilities  -
transmat = ex_transmat(nt);
mdi = loadmd(mpc, transmat, xgd, [], [], profiles);
mdo = most(mdi, mpopt);
if verbose
ms = most_summary(mdo);
end

However the added generator units of which there are 7 of them in
wind_uc_30 don't seem to get the profile set to them in the output. Instead
only the first Generator unit gets the profile while the rest of them are
set at PMAX.

windprofile = struct( ...
'type', 'mpcData', ...
'table', CT_TGEN, ...
'rows', 1, ...
'col', PMAX, ...
'chgtype', CT_REL, ...
'values', [] );

 windprofile.values(:, :, 1) = [
  0.80;
  0.65;
  0.60;
  0.82;
  1.00;
  0.70;
  0.50;
  0.85;
  1.00;
  1.10;
  1.06;
  0.95;

Above is a snippet from the wind_profile.  I am assuming it has something
to do with the row count, however, I haven't been able to find a solution.

Could I please get some help. to ensure all the generator units receive the
load profile.

Thankyou,

-- 
J.Sebben

Re: Wind generation

[Abhyankar, Shrirang G.]

No.

> On May 22, 2016, at 5:25 PM, Mohamed Taleb  wrote:
> 
> Dears,
> 
> Greetings,
> I there any modeling for the wind turbine that simulates its behavior?
> 
> Regards,
> Mohamed Taleb

Wind generation

[Mohamed Taleb]

Dears,

Greetings,
I there any modeling for the wind turbine that simulates its behavior?

Regards,
Mohamed Taleb