Re: IEEE 6-Bus System

[Jose Luis Marin]

Hello Zahid,

Your questions below are on *State Estimation*, which is probably off-topic
for the this mailing list.  May I suggest that you direct this type of
questions to a forum such as Eng-Tips, for instance (
http://www.eng-tips.com/threadminder.cfm?pid=238).  Maybe others in the
mailing list can suggest other better forums.

Two good references on State Estimation are:
   * http://www.springer.com/engineering/electronics/book/978-0-7923-8519-6
   * http://www.crcpress.com/product/isbn/9780824755706

See also http://www.eng-tips.com/faqs.cfm?fid=1287

Hope it helps,

-- 
Jose L. Marin
Gridquant España SL
Grupo AIA



On Sat, Jan 3, 2015 at 5:34 AM, zahid khan  wrote:

> Dear Jose Luis Marin
>
> Thank you very much for your explanation.One thing is not clear to me that
> is regarding IEEE DATA set e.g IEEE 14 bus data or 1EEE 30 bus data
> etc.In state estimation of power system we need two type  of data sets one
> is Branch data and other is measurements like real and reactive powers
> voltages at nodes and we are interested to find unkown state vector which
> consist of voltage magnitudes and phase angles.In measurements data 
> signma(standard
> deviation ) is must given for weighting the observations.My question are
> 1:  How we can know about sigma in real data set?
> 2: Branch data or line data are same in all books written on power system
> analysis but measurements data is different in IEEE sources and other
> books why?
> 3:In topic of state estimation of Wood & Wollenberg book , 100  base is
> used for real and reactive power but what base is used for voltage
> measurements?
> i am requested to all MATPOWER discussion forum to help me in
> understanding of these issues.This is also beneficial for all using these
> data set because reliable data ensure reliable results.
>
> Your answers will be highly appreciated.
>
> THANKS
>
>
> On Tue, Dec 30, 2014 at 5:29 PM, Jose Luis Marin 
> wrote:
>
>>
>> As far as I know, there is no 6-bus test case from the IEEE.  But there
>> are actually two different 6-bus test cases out there:
>>
>>1. The one used by Wood & Wollenberg in their book (Section 6.15 of
>>the 3rd Ed.).  This one is "case6ww.m" in the MATPOWER distribution.
>>2. The one used by Garver, which appears in many articles on
>>Transmission Expansion Planning (http://lmgtfy.com/?q=6+bus+Garver)
>>
>>
>> Hope it helps,
>>
>> --
>> Jose L. Marin
>> Gridquant España SL
>> Grupo AIA
>>
>>
>> On Tue, Dec 30, 2014 at 2:29 AM, zahid khan 
>> wrote:
>>
>>> Thank you Murugan .Could you please tell me ,is it standard data like
>>> other available IEEE data set?I am using this book but little bit confuse
>>> about it.I am view author has taken hypothetical data for understanding of
>>> different topics.I am basically statistician and working in one of the
>>> project related to power system.
>>>
>>>
>>> Thanks
>>>
>>>
>>>
>>> On Mon, Dec 29, 2014 at 9:09 PM, murugan R 
>>> wrote:
>>>
 Yes you can get from generation ,transmission and control book written
 by WOOD AND WOOLENBERG

 From: zahid khan 
 Sent: Mon, 29 Dec 2014 13:38:04
 To: MATPOWER discussion forum 
 Subject: IEEE 6-Bus System

 Dear All

 Could anyone tell me where i can get line data for  IEEE 6 bus system?I
 am not from power system engineering background and hopefully users of this
 forum may help me in this matter.


 Thank you very much.

 zahid khan
 University Technology PETRONAS


 
 Get your own *FREE* website, *FREE* domain & *FREE* mobile app with
 Company email.
 *Know More >*
 
>>>
>>>
>>>
>>
>

Re: IEEE 6-Bus System

[zahid khan]

Dear Jose Luis Marin

Thank you very much for your help and suggestions.

Regards

zahid khan




On Tue, Jan 13, 2015 at 5:42 PM, Jose Luis Marin 
wrote:

> Hello Zahid,
>
> Your questions below are on *State Estimation*, which is probably
> off-topic for the this mailing list.  May I suggest that you direct this
> type of questions to a forum such as Eng-Tips, for instance (
> http://www.eng-tips.com/threadminder.cfm?pid=238).  Maybe others in the
> mailing list can suggest other better forums.
>
> Two good references on State Estimation are:
>*
> http://www.springer.com/engineering/electronics/book/978-0-7923-8519-6
>* http://www.crcpress.com/product/isbn/9780824755706
>
> See also http://www.eng-tips.com/faqs.cfm?fid=1287
>
> Hope it helps,
>
> --
> Jose L. Marin
> Gridquant España SL
> Grupo AIA
>
>
>
> On Sat, Jan 3, 2015 at 5:34 AM, zahid khan  wrote:
>
>> Dear Jose Luis Marin
>>
>> Thank you very much for your explanation.One thing is not clear to me
>> that is regarding IEEE DATA set e.g IEEE 14 bus data or 1EEE 30 bus data
>> etc.In state estimation of power system we need two type  of data sets one
>> is Branch data and other is measurements like real and reactive powers
>> voltages at nodes and we are interested to find unkown state vector
>> which consist of voltage magnitudes and phase angles.In measurements data
>> signma(standard deviation ) is must given for weighting the
>> observations.My question are
>> 1:  How we can know about sigma in real data set?
>> 2: Branch data or line data are same in all books written on power system
>> analysis but measurements data is different in IEEE sources and other
>> books why?
>> 3:In topic of state estimation of Wood & Wollenberg book , 100  base is
>> used for real and reactive power but what base is used for voltage
>> measurements?
>> i am requested to all MATPOWER discussion forum to help me in
>> understanding of these issues.This is also beneficial for all using these
>> data set because reliable data ensure reliable results.
>>
>> Your answers will be highly appreciated.
>>
>> THANKS
>>
>>
>> On Tue, Dec 30, 2014 at 5:29 PM, Jose Luis Marin 
>> wrote:
>>
>>>
>>> As far as I know, there is no 6-bus test case from the IEEE.  But there
>>> are actually two different 6-bus test cases out there:
>>>
>>>1. The one used by Wood & Wollenberg in their book (Section 6.15 of
>>>the 3rd Ed.).  This one is "case6ww.m" in the MATPOWER distribution.
>>>2. The one used by Garver, which appears in many articles on
>>>Transmission Expansion Planning (http://lmgtfy.com/?q=6+bus+Garver)
>>>
>>>
>>> Hope it helps,
>>>
>>> --
>>> Jose L. Marin
>>> Gridquant España SL
>>> Grupo AIA
>>>
>>>
>>> On Tue, Dec 30, 2014 at 2:29 AM, zahid khan 
>>> wrote:
>>>
 Thank you Murugan .Could you please tell me ,is it standard data like
 other available IEEE data set?I am using this book but little bit confuse
 about it.I am view author has taken hypothetical data for understanding of
 different topics.I am basically statistician and working in one of the
 project related to power system.


 Thanks



 On Mon, Dec 29, 2014 at 9:09 PM, murugan R >>> > wrote:

> Yes you can get from generation ,transmission and control book written
> by WOOD AND WOOLENBERG
>
> From: zahid khan 
> Sent: Mon, 29 Dec 2014 13:38:04
> To: MATPOWER discussion forum 
> Subject: IEEE 6-Bus System
>
> Dear All
>
> Could anyone tell me where i can get line data for  IEEE 6 bus
> system?I am not from power system engineering background and hopefully
> users of this forum may help me in this matter.
>
>
> Thank you very much.
>
> zahid khan
> University Technology PETRONAS
>
>
> 
> Get your own *FREE* website, *FREE* domain & *FREE* mobile app with
> Company email.
> *Know More >*
> 



>>>
>>
>

Re: Opf convergence problem

[Ray Zimmerman]

The constraints corresponding to columns 11-21 in the gen matrix are optional 
(i.e. the columns can contain all zeros), but those columns are not to be 
treated as optional columns. It turns out that some functions, such as runopf() 
will not cause an error and, yes, the results returned are correct. However, 
other MATPOWER functions will fail if there are missing columns (e.g. savecase).

I suppose I could just modify loadcase() to add any missing columns, so that we 
can explicitly treat them as optional columns … just added that to the to-do 
list.

Ray


> On Jan 12, 2015, at 1:00 PM, Simone Cochi  wrote:
> 
> Does the missing value in the gen struct affect the results? I thought they 
> were optional constraints
> 
>> Il giorno 12/gen/2015, alle ore 17:24, Ray Zimmerman > > ha scritto:
>> 
>> (First of all, this is not a valid MATPOWER case file … the gen matrix needs 
>> to have 21 columns).
>> 
>> I’m pretty sure these are simply numerical issues. Your network has a huge 
>> range for branch reactances (from 1e-4 up to 7.2415), and I’m sure the zero 
>> cost generation doesn’t help. It turns out that some solvers do converge 
>> (IPOPT, MINOPF, TRALM) using the default starting point, and others (MIPS 
>> and KNITRO) converge if you use the power flow solution to initialize the 
>> OPF.
>> 
>> Ray
>> 
>> E.g.
>> 
>> >> load mpc
>> >> mpopt = mpoption('verbose', 2, 'opf.ac.solver', 'MIPS', 
>> >> 'out.suppress_detail', 1, 'opf.init_from_mpc', 1);
>> >> r = runopf(mpc, mpopt);   
>> >> 
>> MATPOWER Version 5.0, 17-Dec-2014 -- AC Optimal Power Flow
>> MATLAB Interior Point Solver -- MIPS, Version 1.1, 17-Dec-2014
>>  itobjective   step size   feascond gradcond compcond 
>> costcond  
>>    -    
>> 
>>   0  2.71  0.003209520.005  394.161  
>>   0
>>   1 13.5613670.454670.0131881  49.7247  129.098   
>> 0.00108484
>>   2 2.5810166 0.2459  0.000347117  5.81263  30.6224   
>> 0.00109655
>>   3 3.0946997  0.0073205  1.34869e-07   0.5504  3.18123  
>> 5.13551e-05
>>   4 3.1139811 0.00033848  8.72888e-100.0214386 0.318187  
>> 1.92754e-06
>>   5 3.1489261 0.00061552  2.90301e-09 0.0183710.0318179  
>> 3.49342e-06
>>   6 3.0091397  0.00246214.645e-080.0117079   0.00598559  
>> 1.39742e-05
>>   7  2.710015  0.0068676  2.23731e-070.0199562  0.000603092  
>> 2.99035e-05
>>   8 2.7116241 2.7338e-05  5.10117e-12  8.49417e-06  6.03182e-05  
>> 1.60873e-07
>>   9 2.7101543 2.5902e-05  5.14787e-12  8.06346e-07  6.03182e-06  
>> 1.46939e-07
>>  10 2.7100156 2.4442e-06  1.27526e-13  7.54021e-09  6.03182e-07   
>> 1.3865e-08
>> Converged!
>> 
>> Converged in 0.06 seconds
>> Objective Function Value = 2.71 $/hr
>> 
>> | System Summary 
>>   |
>> 
>> 
>> How many?How much?  P (MW)Q (MVAr)
>> ----  -  
>> -
>> Buses192 Total Gen Capacity  25.0 -25.0 to 25.0
>> Generators 4 On-line Capacity25.0 -25.0 to 25.0
>> Committed Gens 4 Generation (actual)  3.5   1.7
>> Loads 73 Load 3.4   1.6
>>   Fixed   72   Fixed  3.3   1.6
>>   Dispatchable 1   Dispatchable   0.1 of 0.1   -0.0
>> Shunts 0 Shunt (inj) -0.0   0.0
>> Branches 191 Losses (I^2 * Z) 0.05  0.16
>> Transformers 191 Branch Charging (inj) -0.0
>> Inter-ties 0 Total Inter-tie Flow 0.0   0.0
>> Areas  1
>> 
>>   Minimum  Maximum
>>  -  
>> Voltage Magnitude   0.957 p.u. @ bus 1711.000 p.u. @ bus 1   
>> Voltage Angle  -2.53 deg   @ bus 1580.00 deg   @ bus 1   
>> P Losses (I^2*R) -  0.01 MW@ line 1-2
>> Q Losses (I^2*X) -  0.09 MVAr  @ line 1-2
>> Lambda P0.00 $/MWh @ bus 1  0.00 $/MWh @ bus 171 
>> Lambda Q0.00 $/MWh @ bus 1  0.00 $/MWh @ bus 171 
>> 
>> 
>> 
>>> On Jan 10, 2015, at 5:02 AM, Simone Cochi >> > wrote:
>>> 
>>> Dear all, 
>>> I can’t get the convergence for the OPF on a radial distribution network 
>>> described b

Re: Power system modelling with transformer ratio and shift.

[Mirish Thakur]

Dear Mr. Abhyankar .G.,

I'm not sure about it and I'm not getting correct way  how to extract the
elements from solution vector you are talking about. But what I've done is-
1) I didn't put tap ratio and shift values in branch matrix. I have created
another matrix data for transformer like mpc.tapdata in matpower case file
(case9.m) and assigned values with lower and upper bounds for both ratio
and shift. Also created idx_tapdata where I have defined constants for
named column.
2) Added these extra  variables and their constraints in opf_setup  by add
_vars and add_constraints. (Total variables were 24 and now they are 42 in
Jacobian matrix verified by spy function.)
3) Created new function where I get ratio and  shift values.

Problem what I'm facing is  how to extract  values of ratio and shift  from
solver to update the Y bus in opf_hessfcn  and opf_consfcn ? But facing the
problem in opf_args.

Thanks .
Mirish

On Fri, Jan 9, 2015 at 3:05 PM, Abhyankar, Shrirang G. 
wrote:

>  Are you setting the solution vector (x) size and the indexing
> information for extracting the elements from it correctly?
>
> On Jan 8, 2015, at 4:50 AM, "Mirish Thakur" 
> wrote:
>
>Dear Prof. R. Zimmerman,
>
>  I'm master student working on a  project dealing with power system
> optimization with additional variables  'tap ratio' and 'shift' on 9 bus
> system by  *mipsopf_solver*. I have made necessary changes in respective
> function files. Also I concern with only current flow constraints in
> *opf_hessfcn*. For that I used command like *>>* *mpopt =
> mpoption('OPF_FLOW_LIM', 2);* and* >> runopf(case9,mpopt);*. But I'm
> getting errors -
>
> Did not converge in 17 iterations.
> Error using opf_args (line 315)
> opf_args.m: Incorrect input parameter order, number or type
>
> Error in opf (line 199)
> [mpc, mpopt] = opf_args(varargin{:});
>
> Error in opf (line 262)
> opf(varargin)
>
> Error in runopf (line 96)
> [r, success] = opf(casedata, mpopt);
>
>  As I'm not much faliliar with MATPOWER so I could not understand what
> exactly I need to rectify in *opf_args.m* file. I need your help .Thanks
> in advance.
>
>
>  Sincerely
>  Mirish
>
>

Re: Power system modelling with transformer ratio and shift.

[Abhyankar, Shrirang G.]

Take a look at how Sbus is created in opf_consfcn by extracting Pg, and Qg from 
the solution vector x. You need to do something similar with Ybus.

Shri

From: Mirish Thakur mailto:mirishtha...@gmail.com>>
Reply-To: MATPOWER discussion forum 
mailto:matpowe...@list.cornell.edu>>
Date: Wed, 14 Jan 2015 00:12:16 +0100
To: MATPOWER discussion forum 
mailto:matpowe...@list.cornell.edu>>
Subject: Re: Power system modelling with transformer ratio and shift.

Dear Mr. Abhyankar .G.,

I'm not sure about it and I'm not getting correct way  how to extract the 
elements from solution vector you are talking about. But what I've done is-
1) I didn't put tap ratio and shift values in branch matrix. I have created 
another matrix data for transformer like mpc.tapdata in matpower case file 
(case9.m) and assigned values with lower and upper bounds for both ratio and 
shift. Also created idx_tapdata where I have defined constants for named column.
2) Added these extra  variables and their constraints in opf_setup  by add 
_vars and add_constraints. (Total variables were 24 and now they are 42 in 
Jacobian matrix verified by spy function.)
3) Created new function where I get ratio and  shift values.

Problem what I'm facing is  how to extract  values of ratio and shift  from 
solver to update the Y bus in opf_hessfcn  and opf_consfcn ? But facing the 
problem in opf_args.

Thanks .
Mirish

On Fri, Jan 9, 2015 at 3:05 PM, Abhyankar, Shrirang G. 
mailto:abhy...@mcs.anl.gov>> wrote:
Are you setting the solution vector (x) size and the indexing information for 
extracting the elements from it correctly?

On Jan 8, 2015, at 4:50 AM, "Mirish Thakur" 
mailto:mirishtha...@gmail.com>> wrote:

Dear Prof. R. Zimmerman,

I'm master student working on a  project dealing with power system optimization 
with additional variables  'tap ratio' and 'shift' on 9 bus system by  
mipsopf_solver. I have made necessary changes in respective function files. 
Also I concern with only current flow constraints in opf_hessfcn. For that I 
used command like >> mpopt = mpoption('OPF_FLOW_LIM', 2); and >> 
runopf(case9,mpopt);. But I'm getting errors -

Did not converge in 17 iterations.
Error using opf_args (line 315)
opf_args.m: Incorrect input parameter order, number or type

Error in opf (line 199)
[mpc, mpopt] = opf_args(varargin{:});

Error in opf (line 262)
opf(varargin)

Error in runopf (line 96)
[r, success] = opf(casedata, mpopt);

As I'm not much faliliar with MATPOWER so I could not understand what exactly I 
need to rectify in opf_args.m file. I need your help .Thanks in advance.


Sincerely
Mirish

error in installing matpower 5

[Zeya Aung [11547809d]]

Hi all

I can’t install matpower 5. says error below. Please help.

t_nested_struct_copyok
t_mpoption..ok
t_loadcase..Error using save
Unable to write file t_mat9_opf.mat: No such file or directory.

Error in t_loadcase (line 59)
eval(['save ' matfile '.mat baseMVA bus gen branch areas gencost']);

Error in t_run_tests (line 79)
feval( test_names{k}, ~verbose );

Error in test_matpower (line 146)
t_run_tests( tests, verbose );

Best
Zeya

   
-
Disclaimer:

This message (including any attachments) contains confidential information 
intended for a specific individual and purpose. If you are not the intended 
recipient, you should delete this message and notify the sender and The Hong 
Kong Polytechnic University (the University) immediately. Any disclosure, 
copying, or distribution of this message, or the taking of any action based on 
it, is strictly prohibited and may be unlawful.

The University specifically denies any responsibility for the accuracy or 
quality of information obtained through University E-mail Facilities. Any views 
and opinions expressed are only those of the author(s) and do not necessarily 
represent those of the University and the University accepts no liability 
whatsoever for any losses or damages incurred or caused to any party as a 
result of the use of such information.