Pb1 y13 =-j10 Pb5. Pb4. y34 =-j10

Similar documents
Pb1 y13 =-j10 Pb5. Pb4. y34 =-j10

ECEN 460 Exam 1 Fall 2018

ECE 476 Power System Analysis Fall 2014 Exam #1, Thursday, October 2, :30AM - 10:50AM

SHORT QUESTIONS AND ANSWERS. Year/ Semester/ Class : III/ V/ EEE Academic Year: Subject Code/ Name: EE6501/ Power System Analysis

A Model for a Zonal Operating Reserve Demand Curve

Calculations of Capacitance for Transposed Bundled Conductor Transmission Lines

ECE 476. Exam #2. Tuesday, November 15, Minutes

THE UNIVERSITY OF NEW SOUTH WALES. School of Electrical Engineering & Telecommunications FINALEXAMINATION. Session

16PESGM2316 Characterizing Transmission System Harmonic Impedances with R-X Loci Plots. David Mueller

EE 6501 POWER SYSTEMS UNIT I INTRODUCTION

Software Tools: Congestion Management

EE 451 Power System Stability

Lecture 11 Transmission Line Parameters

The L-Shaped Method. Operations Research. Anthony Papavasiliou 1 / 38

Proper Security Criteria Determination in a Power System with High Penetration of Renewable Resources

Javier Contreras Sanz- Universidad de Castilla-La Mancha Jesús María López Lezama- Universidad de Antioquia Antonio Padilha-Feltrin- Universidade

Transmission and Distribution of Electrical Power

International Studies about the Grid Integration of Wind Generation

EVALUATION OF WIND ENERGY SOURCES INFLUENCE ON COMPOSITE GENERATION AND TRANSMISSION SYSTEMS RELIABILITY

ELG4125: Power Transmission Lines Steady State Operation

EE5250 TERM PROJECT. Report by: Akarsh Sheilendranath

A technique for simultaneous parameter identification and measurement calibration for overhead transmission lines

SSC-JE EE POWER SYSTEMS: GENERATION, TRANSMISSION & DISTRIBUTION SSC-JE STAFF SELECTION COMMISSION ELECTRICAL ENGINEERING STUDY MATERIAL

Generation Expansion Planning

Lecture (5) Power Factor,threephase circuits, and Per Unit Calculations

EE 581 Power Systems. Admittance Matrix: Development, Direct and Iterative solutions

KINGS COLLEGE OF ENGINEERING Punalkulam

EECE421 Power System Analysis. Chapter 4: Transmission Line Capacitance

B.E. / B.Tech. Degree Examination, April / May 2010 Sixth Semester. Electrical and Electronics Engineering. EE 1352 Power System Analysis

APPENDIX 7.4 Capacity Value of Wind Resources

Some Observations on the Locational Marginal Prices and Loads

A Nodal Capacity Market to Assure Resource Adequacy

ECE 325 Electric Energy System Components 5 Transmission Lines. Instructor: Kai Sun Fall 2015

System Planning Lecture 7, F7: Optimization

Capacity Planning with uncertainty in Industrial Gas Markets

Power System Analysis Prof. A. K. Sinha Department of Electrical Engineering Indian Institute of Technology, Kharagpur

A Benders Decomposition Approach to Corrective Security Constrained OPF with Power Flow Control Devices

Congestion Alleviation using Reactive Power Compensation in Radial Distribution Systems

Fault Analysis Power System Representation

The L-Shaped Method. Operations Research. Anthony Papavasiliou 1 / 44

The DC Optimal Power Flow

EE 742 Chapter 3: Power System in the Steady State. Y. Baghzouz

RTO Winter Resource Adequacy Assessment Status Report

DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING QUESTION BANK

EE2351 POWER SYSTEM ANALYSIS UNIT I: INTRODUCTION

Chapter 5. Transmission networks and electricity markets

3W Transformer Test Report Impedance Calculator Equations. Use these Equations for: Primary Secondary Primary Tertiary Secondary - Tertiary

A Decomposition Based Approach for Solving a General Bilevel Linear Programming

Increasing Transmission Capacities with Dynamic Monitoring Systems

CHAPTER 2 LOAD FLOW ANALYSIS FOR RADIAL DISTRIBUTION SYSTEM

Decision Models Lecture 5 1. Lecture 5. Foreign-Currency Trading Integer Programming Plant-location example Summary and Preparation for next class

Chapter 2. Planning Criteria. Turaj Amraee. Fall 2012 K.N.Toosi University of Technology

Analyzing the Effect of Loadability in the

Optimal Placement & sizing of Distributed Generator (DG)

Tutorial 2: Modelling Transmission

IGEE 402 Power System Analysis. FINAL EXAMINATION - SAMPLE Fall 2004

Impact Study on Power Factor of Electrical Load in Power Distribution System

SINGLE OBJECTIVE RISK- BASED TRANSMISSION EXPANSION

Performance Improvement of the Radial Distribution System by using Switched Capacitor Banks

In the previous chapter, attention was confined

Power Grid Partitioning: Static and Dynamic Approaches

J. Electrical Systems x-x (2010): x-xx. Regular paper

INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous)

Then the other two numbers may be represented as x + 1 and x + 2. Now we can represent 3 less than the second as (x + 1) 3.

Mixed Integer Linear Programming Formulation for Chance Constrained Mathematical Programs with Equilibrium Constraints

A three-level MILP model for generation and transmission expansion planning

Distributed vs Bulk Power in Distribution Systems Considering Distributed Generation

4. Duality and Sensitivity

Module 3 : Sequence Components and Fault Analysis

An Introduction to the Smith Chart for Amateur Radio. Jesse Sheinwald, N2CA

QUESTION BANK ENGINEERS ACADEMY. Power Systems Power System Stability 1

Research on Calculation Method of Power Supply Radius Based on Voltage Deviation in Rural Power Network. Ning Wang

Improvements to Benders' decomposition: systematic classification and performance comparison in a Transmission Expansion Planning problem

Stochastic Generation Expansion Planning

Branch Outage Simulation for Contingency Studies

A Unified Framework for Defining and Measuring Flexibility in Power System

The Power Flow & Optimal Power Flow Problems

Decomposition Techniques in Mathematical Programming

Evaluation of GA Performance in TNEP Considering Voltage Level, Network Losses and Number of Bundle Lines

ECE 325 Electric Energy System Components 7- Synchronous Machines. Instructor: Kai Sun Fall 2015

Power Engineering II. Fundamental terms and definitions

elgian energ imports are managed using forecasting software to increase overall network e 칁 cienc.

03-Economic Dispatch 1. EE570 Energy Utilization & Conservation Professor Henry Louie

A Merchant Mechanism for Electricity Transmission Expansion

Multi-Area Stochastic Unit Commitment for High Wind Penetration in a Transmission Constrained Network

A possible notion of short-term value-based reliability

Elevated Neutral to Earth Voltages Due to Harmonics A T&D Update

0-2 Operations with Complex Numbers

Fault Calculation Methods

Selecting the current rating for equipment

Outage Coordination and Business Practices

PROBLEM SOLUTIONS: Chapter 2

0-2 Operations with Complex Numbers

On-line supplement to: SMART: A Stochastic Multiscale Model for the Analysis of Energy Resources, Technology

Harmonic Modeling of Networks

TRANSMISSION LINES. All aluminum alloy conductor (AAAC) Aluminum conductor alloy reinforced (ACAR)

EE2351 POWER SYSTEM ANALYSIS

PROJECT ECONOMIC ANALYSIS

A PROBABILISTIC MODEL FOR POWER GENERATION ADEQUACY EVALUATION

Australian Journal of Basic and Applied Sciences. General Fault Admittance Method Solution of a Balanced Line-to-Line-to-Line Fault

Transcription:

EE 55, Exam, ake-home. Due Monday, April, 06, 5:00pm. You may use class notes or any reference materials (e.g., books, etc.) that you like; however, you must work alone, i.e., you should not be communicating with anyone else about how to work the problems on this exam.. Generation planning (5 pts): Develop the complete set of equations for generation expansion plan formulation -b (GEP-b), based on the following data and assumptions. a. Plan for 30 years, at 5 year intervals (thus you should have 6 investment periods (at years 5, 0, 5, 0, 5, and 30). b. Use annual discount rate of 9%. c. Maintain a reserve of 5% above peak load. d. Assume a CO cost of $5/M (you will need to add a term to the objective function to account for this). e. ransmission topology and data: wo comments follow Note the definitions of positive branch flow for each P bk, k=,,5. he admittances y ij are given on a 00 MVA power base. Pg Pg Pb3 y =-j0 Pb3 y4 =-j0 Pb y3 =-j0 Pb5 Pd y3 =-j0 4 Pb4 y34 =-j0 3 Pg4 Pd3 f. Line flow constraints: -500MW P bk 5000MW, k=,5. he per-unit expressions of these constraints, on a 00 MVA power base, are -500 P bk 500. Given the level of loads (see item (j) below), this implies the transmission has infinite capacity. g. Existing generation data: Bus Bus Bus 4 echnology PC (j=) NGCC (j=3) C (j=4) Capacity (MW) 00 00 50 h. echnology characteristics PC (j=) PC (j=) NGCC(j=3) C (j=4) Wind (j=5) Capacity credit.0.0.0.0 0.5 Investment cost (000$/MW) 3000 97 9 968 980 Fixed O&M costs (000$/MW/yr) 50 3 3 7 40 Var O&M ($/MWhr) 0 4.5 3.6 5.5 0 Full-load heat rate (MBU/MWhr) 0.0 8.8 7.05 -- Fuel cost ($/MBU) 3 3 4 4 0 CO emissions (lbs/mwhr) 0 87 85 64 0 Capacity factor 0.35 Salvage value (000$/MW) 50 50 40 0 0

i. Load duration curve and corresponding blocks. Assume a load growth rate of %/year, applied to each value d, d, d 3, for each of the two loads in the system, with the durations h, h, and h 3 remaining the same. herefore, for example, for the load at bus, d would be.0,.084,.95,.395,.4568,.6084,.7758 pu for years, 5, 0, 5, 0, 5, and 30, respectively. You will need to compute the values of d and d 3 for these years as well. hen you will need to do the same for the load at bus 3. d d L (MW) d3 h h+h h+h+h3 Number of hours that Load > L Initial load at bus (Pg) Block Block Block 3 Load (pu) d = d = d 3 =0.6 Duration (hrs) h =000 h =5000 h 3 =760 Initial load at bus 3 (Pg3) Block Block Block 3 Load (pu) d =. d =0.9 d 3 =0.7 Duration (hrs) h =900 h =500 h 3 =660. Production costing ( pts): a. (4 pts) A load duration curve characteries the next year for a power system is given below. he power system is supplied by one 3 MW unit. Identify the minimum load, the loss of load probability, the loss of load expectation, and the expected unserved energy, assuming the 3 MW unit is perfectly reliable. 0.6 F Dr ( d e ) 0. d e (MW) 3 4 5 6 7 8 Solution: Minimum Load= MW, LOLP=P(d>3)=0.333, LOLE=0.333*8760=97 hrs EUE=Area*=0.5()(.333)(8760)=0.666(8760)=460MWhrs

Solution: b. (5 pts) Now consider that the power system is supplied with one 4 MW unit having a forced outage rate of 0.5. he load is still characteried by the load duration curve above. Compute loss of load probability, loss of load expectation, and expected unserved energy. f Dj (d j ) 0.6 F Dr ( d e ) * 0. d e (MW) 0,6 0. 0.75 0.5 (b) C =4 3 4 5 6 7 8 3 4 5 6 7 8. 0.6 0. 0.75 (c) +. 0.6 0. 0.5 (d) 3 4 5 6 7 8 d e 3 4 5 6 7 8 d e =.0 0.6 F r ( () d D e ) (e) 0. 3 4 5 6 7 8 d e So LOLP=P(d>4)=0.5, LOLE=0.5*8760=90hrs EUE=8760{()(0.5)+0.5(3)(0.5)}=8760{.5+.375}=5475MWhrs c. (3 pts) Consider the two situations you have assessed in parts (a) and (b) of this problem. Complete the following sentence: Although situation has the larger LOLP; situation has the larger EUE because. Solution: Although situation (a) has a larger LOLP, situation (b) has a larger EUE because when load is interrupted in situation (a), it will only be for 0 to MW, whereas when load is interrupted in situation (b), it will be for 0 to 4 MW.

d. (4 pts) Compute the energy expected to be produced by the 4 MW unit of the year. Solution: E=A**Area where Area is computed from the previous load duration curve, i.e., the one before the unit was convolved in, according to Area x x j j F ( j) D e ( ) d his is the load duration curve given in the problem statement. In this case, x j- =0 and x j =4, and so the area is A=()()+0.5(3)()=.5. he availability is 0.75 so that E=0.75*8760*.5=645 MWhrs. e. (3 pts) Assume the 4 MW unit is a thermal unit. What is the minimum additional information you need to know to compute the cost incurred by unit 4 in producing the amount of energy developed in part (d)? Solution: We need the heat rate and the fuel price. f. (3 pts) After reviewing a production cost program which computes annual production costs in a way that is similar to the procedure suggested in parts (a)-(c) of this problem, an engineer decides to write a production cost simulation program that will allow representation of the power system with much greater operational fidelity. What is the single most important change that the engineer needs to make, relative to that suggested in parts (a)-(e)? Solution: he engineer needs to represent the power system chronologically, hour by hour. 3. Benders (7 pts): he relation between primal and dual optimiations can be stated as: Primal max s. t. x 0 c A x b subproblem A w * Consider the following problem: max subject to : x 4x 5w w 0, : x w 0 x 3w 8 w, x 0, w Dual subproblem : min s. t. A integer 0 c b A w* a. (4 pts) Express this problem in the form of Benders, i.e., express an appropriate master problem and an appropriate (primal) subproblem.

Solution: Master problem: max s.t 5w Q : w 0, Subproblem: x 0 * * max 4x subject to x 0 w x 8 3w M, b. (3 pts) Identify an initial solution to the master. Solution: * M, w 0 w 0, w integer c. (4 pts) Formulate the subproblem dual in terms of w*, and based on the dual formulation, and the solution to the master problem found in part (b), explain why you know that the primal subproblem is infeasible. Solution: We first express the subproblem dual. c 4, c 5, A, A, 3 b min s. t. 0 A c ( b A w*), 0 Substituting w*=0, the subproblem dual becomes: min s. t. 4, 0 0 0 8 3 min s. t. 0 5 4 0 8 0 * 8 3 w

Because the all coefficients in the objective are negative, and the variables are constrained to be non-negative, the subproblem dual is unbounded; an unbounded dual necessarily implies that the corresponding subproblem primal must be infeasible. d. (3 pts) Given the sub-problem primal is infeasible, explain why it is possible that we may still find a solution to this problem (and thus, that this solution will be feasible in the subproblem primal). Solution: he reason that we may still find a solution to this problem, in spite of the fact that the subproblem primal is infeasible in this step, is that the subproblem primal is expressed in terms of w*. If we further constrain w*, it will affect the subproblem primal s feasibility region. e. (3 pts) Find a constraint on w which will make the subproblem primal feasible. Solution: he feasibility constraint can be found by expressing (b-a w) λ 0, or 0 w 8 3 which results in from which we observe that 0 (0 w) (8 3 w) 0 0 w 0 83w 0 or 8 0 w w 3 And we conclude that w 8/3. 4. Planning principles (0 pts): A good paper on planning (you should read in its entirety - it can be found in IEEE XPlore) is: M. Awad, S. Broad, K. Casey, J. Chen, A. Geevarghese, J. Miller, A. Pere, A. Sheffrin, M. Zhang, E. oolson, G. Drayton, A. Rahimi, and B. Hobbs, he California ISO transmission economic assessment methodology (EAM): principles and application to Path 6, Power Engineering Society General Meeting, 006. his paper lists five principles to economic evaluation of proposed transmission upgrades. I have chosen from the paper an excerpt for each principle. For each one, identify how you would implement this principle. Where appropriate, base your answer on information, concepts, and/or methods we have covered in EE 55. I am expecting from you at minimum a -3 sentence paragraph for each principle. First key principle: Benefit Framework: Second key principle: Network Representation:

hird key principle: Market Prices: Fourth key principle: Uncertainty:

Fifth key principle: Resource alternatives to transmission expansion: Fifth key principle (continued):

5. ransmission line parameters (6 pts): A 765 kv transmission line has phase separation of 36 feet. Each phase is a 4-conductor bundle with 4-inch diameter. he conductor type is ern (795 kcmil). he following data is obtained from tables: Ind reactance ft spacing for 4-conductor bundle at 4 : X a =0.08Ω/mile Ind reactance spacing factor for 36 phase separation: X d =348Ω/mile Cap reactance at ft spacing for 4 conductor bundle at 4 : X a =0.005MΩ-mile Cap reactance spacing factor for 36 phase separation: X d =0.06MΩ-mile a. (6 pts) Compute the SIL for this line configuration with 4 conductors per bundle. Solution: Get per-unit length inductive reactance: So X L =X a +X d =0.08+348=68 ohms/mile. Now get per-unit length capacitive reactance. X C =X a +X d =0.005+0.06=0.3E6ohms-mile So =jx L =j376 Ohms/mile, and y=/-jx C =/-j(0.3 0 6 )=j8.9847 0-6 mhos/mile j.468 Z C 6.96ohms y -6 j8.9847 0 3 VLL 7650 PSIL.5785e + 009 Z 6.96 C he SIL with 4 conductors is 578 MW. b. (4 pts) A designer considers reducing the number of conductors per bundle from 4 to 3, while maintaining conductor spacing and phase spacing at 4 inches and 36 feet, respectively. Would you expect this change to increase or decrease the surge impedance loading? Support your answer. Solution: he 3 conductor bundle is an equilateral triangle with a GMR of (r *4*4) /3, whereas the 4 conductor bundle is a square with a GMR of (r *4*4*33.94) //4 where 33.94=sqrt(4 +4 ). herefore, the GMR R b of the 3 conductor bundle is slightly smaller than the GMR of the 4 conductor bundle. his causes the 3 conductor bundle to have larger inductive and capacitive reactances, per the below equations, X X L C 3 3.00 f ln.0 0 f ln Dm R b X X a d 6 6.7790 ln.7790 ln f c m R f b X X ' a d /mile D - mile and therefore a larger surge impedance and therefore a smaller SIL. c. (6 pts) he designer desires to maximie the power transfer capability of the line. Relative to the design of part (a) If the line is 30 miles long, what changes would you consider?

Solution: Use a conductor with a larger diameter and therefore higher ampacity. If the line is 300 miles long, what changes would you consider? Solution: Increase SIL by decreasing phase spacing, increasing the number of subconductors in a bundle, or increasing the bundle diameter.