ee/econ 458, syllabus section v...ga gc gb gd ly lz branch reactance [ p.u.] capacit y [ mw] 1-2 0.2...

70
© 2005 D. Kirschen 1 EE/Econ 458, Syllabus Section V.C Power Trading Subject to Transmission Constraints Leigh Tesfatsion Last Revised: 3 November 2011

Upload: others

Post on 26-Jan-2021

2 views

Category:

Documents


0 download

TRANSCRIPT

  • © 2005 D. Kirschen 1

    EE/Econ 458, Syllabus Section V.C

    Power Trading Subject to Transmission Constraints

    Leigh Tesfatsion

    Last Revised: 3 November 2011

  • © 2005 D. Kirschen 2

    Taking Transmission Network Effects Into Account ( Fig. Source: ABB c/o Jamie Austen, PacifiCorps, SSGwi Conference, 2008 )

  • © 2005 D. Kirschen 3

    Traditional Trade-Off Between Market & Network Models ( Fig. Source: ABB c/o Jamie Austen, PacifiCorps, SSGwi Conference, 2008 )

  • © 2005 D. Kirschen 4

    New agent-based modeling tools are increasingly being

    used to study dynamic power market/grid systems … (see www.econ.iastate.edu/tesfatsi/aelect.htm for tools & surveys )

  • © 2005 D. Kirschen 5

    Glossary of Basic Terms for K/S Chapter 6

    (Important for EE/Econ 458 Students who did not take EE 303)

    • Load: Electric energy user (e.g., company, appliance, household)

    • Demand: How much power a load requires (e.g., power to run a

    production process, to turn on a light, to heat a home, etc.)

    • Bulk Transmission Grid: The portion of the transmission grid operating at relatively high voltage levels (138, 161, 230, 345, 500,

    or 765 kV) that is used to transfer electric energy from generators to

    major load centers.

    • Line: A conductor for transmitting electric energy

    • Branch: K/S use “branch” interchangeably with “line”

  • © 2005 D. Kirschen 6

    Glossary of Basic Terms … Continued

    • Bus: A physical transmission element that connects two or more branches in a network representation of a transmission grid

    • Node: K/S use “node” and “bus” interchangeably and assume all buses/nodes are “pricing nodes” with LMPs. U.S. power markets are

    more complicated. (Example: MISO has over 33,000 buses but only

    about 1,540 pricing nodes)

    • Locational Marginal Price (LMP) at Bus k in

    Hour H: The least cost of servicing an additional increment of demand at bus k during hour H while respecting all system constraints.

    • Nodal Price: K/S use “nodal price” interchangeably with LMP.

  • © 2005 D. Kirschen 7

    Glossary of Basic Terms … Continued:

    • System Operator: Organization charged with the primary responsibility of maintaining bulk transmission grid reliability = [ adequacy (balance of power withdrawals and injections) + security (safe operating range) ]

    • Load-Serving Entity: An electric utility, transmitting utility, or Federal power marketing agency that has an obligation under Federal, State, or local law, or under long-term contracts, to provide electric energy to retail customers or to other LSEs with retail customers.

    • (Wholesale) Generator: A company that produces and sells electric energy at the wholesale level through the bulk transmission grid.

    • Electric Energy Demand: The planned purchase of electric energy by one or more buyers. Can be expressed in various forms, including:

    Fixed demand bid = A bid to buy a specified amount of electric energy during a particular time period at a particular bus, regardless of price.

    Price-sensitive demand bid = A bid to buy electric energy during a particular time period at a particular bus that expresses the buyer’s maximum willingness to pay for each successive block of energy.

  • © 2005 D. Kirschen 8

    Glossary of Basic Terms … Continued:

    • Impedance: Takes the complex form z = r + jx, where r (ohms) = resistance, x (ohms) = reactance, j = √-1, & ohm = 1 V/A.

    • Admittance: Inverse of impedance. Takes the complex form y = g + jb, where the conductance g is given by g = r/[r2 + x2] (mhos)

    and the susceptance b is given by b = -x/[r2 + x2] (mhos).

    • Base Voltage Vo and Base Apparent Power So: Used to convert voltage (kV) and power (MVA) magnitudes into dimensionless

    units. Example: Vo = 138 kV, So = 100 MVA.

    • Base Impedance Zo: Given a balanced 3-phase network with Vo in

    line-to-line kV and So in 3-phase MVA, then Zo = Vo2/So (ohms)

    • Per Unit (pu) Normalization: Divide all power system voltage, power, & derived quantities by their base values. Example:

    reactance xkm on branch k→m in p.u. = xkm p.u. = xkm/Zo (dimensionless)

  • © 2005 D. Kirschen 9

    Two Fundamental Laws Governing Current and Power

    Flows in Electrical Networks (cf. K/S Section 6.2.2.1)

    • Kirchhoff’s Current Law (KCL)

    The sum of all currents Imk entering a bus k along branches m→k must equal the sum of all currents Ikm exiting bus k along branches k→m, where

    m is any other bus and I denotes current in Amperes (A).

    KCL implies at any bus k that the inflow/outflow of active (or real) power

    (MW) and of reactive power (MVAR) must each separately be in balance.

    • Kirchhoff’s Voltage Law (KVL)

    The sum of the voltage drops along the branches comprising any closed loop in an electrical network must be zero.

    Implies that the voltage drops along parallel branches must be equal.

    • Complex power flowing on branch k→m: Skm = Pkm + j Qkm = VkI*km

    • Ohm’s Law (AC): Ikm = [Vk – Vm]/zkm (zkm = rkm+ j xkm = impedance)

    voltage drop from k to m (cf. K/S p. 143) Current flowing from k to m

    Text Box Omit

  • © 2005 D. Kirschen 10

    1 2

    zA

    zB

    Application of Kirchhoff’s Voltage Law (KVL): Power

    Transfer Distribution Factors (PTDFs) for Parallel Paths (cf. K/S Section 6.2.2.1)

    IB

    IA

    Impedance

    Impedance

    Bus 1 Bus 2 I12 I12

    Current flowing from Bus 1 into

    Bus 2

    Current flowing

    into Bus 2 from

    Bus 1

    Current I12 = [IA + IB] (Amperes) flows from Bus 1 into Bus 2 along two

    parallel paths with parallel flows IA and IB and impedences zA and zB, (Ohms).

    By Ohm’s Law and KVL, IA = [V1 – V2]/zA and IB = [V1 – V2]/zB .

    Hence zAIA = zBIB = zB[ I12 – IA] [zA + zB]IA = zB I12 IA = I12 • zB/[zA + zB]

    Similarly, IB = I12 • zA/[zA + zB] Text Box Omit

  • © 2005 D. Kirschen 11

    1 2

    xA

    xB

    Power Transfer Distribution Factors (PTDFs) for Parallel

    Paths Neglecting Reactive Power Flows: Corollary of KVL (cf. K/S Section 6.2.2.1)

    FB

    FA

    FA

    x B

    x A x BP F

    Bx A

    x A x BP

    Reactance

    Reactance

    Path Flows (Neglecting Reactive Power and Losses): K/S pp. 143-144

    Bus 1 Bus 2 P P

    Injection of active power

    at Bus 1

    Withdrawal of

    active power

    at Bus 2

    upper path flow lower path flow

    PTDFA PTDFB Text Box Omit

  • © 2011 D. Kirschen and L. Tesfatsion 12

    Effects of Transmission Network Constraints

    on Electric Power Market Operation

    (Kirschen/Strbac Chapter 6, Sections 6.1-6.3.2.9)

    Important Acknowledgement:

    The following slides are based on materials originally prepared by

    Daniel Kirschen (U of Manchester) with edits by Leigh Tesfatsion

    (Iowa State U).

    Last Revised: 3 November 2011

  • © 2011 D. Kirschen and L. Tesfatsion 13

    K/S Chapter 6: Introduction

    • Simple Economic Dispatch formulations examined to date assume

    a given demand that must be met

    all GenCos and LSEs are connected to the same bus

    • Not sufficient for real-world restructured wholesale power markets

    • Need to consider:

    Constraints on transmission lines that restrict power flows

    Line congestion and losses

    Possibility of strategic demand as well as strategic supply behaviors

    • Two forms of trading over a bulk transmission grid will be discussed:

    Bilateral (decentralized) trading between a GenCo and an LSE

    Pool (centralized) trading among GenCos and LSEs

  • © 2011 D. Kirschen and L. Tesfatsion 14

    Bilateral Trading between a GenCo and LSE

    • GenCo/LSE agree on price, quantity, & other conditions

    • GenCo/LSE must schedule their transactions with the system operator, and pay/receive LMP at their buses

    • GenCo/LSE agree to “make whole” any discrepancies between LMPs and their contractual price level

    • System operator must maintain reliability of the grid

    Buys or sells limited amounts of energy to keep power withdrawals

    and injections in balance on the transmission grid (adequacy)

    Limits amount of power that GenCos/LSEs can withdraw/inject at

    some buses if reliability cannot be ensured by other means

  • © 2011 D. Kirschen and L. Tesfatsion 15

    G1

    G2

    L1

    L2

    Bus 1 Bus 2

    G3

    Example of Bilateral Contract Trading:

    • Three GenCos {G1,G2,G3} and two LSEs {L1,L2)

    • G1 sells 300 MW to L1 (planned generation=demand)

    • G2 sells 200 MW to L2 (planned generation=demand)

    • Contractual prices are a private matter

    • Planned quantities (300 MW & 200 MW) must be reported to system operator to ensure grid reliability.

    Line corridor 1→2

  • © 2011 D. Kirschen and L. Tesfatsion 16

    G1

    G2

    L1

    L2

    Bus 1 Bus 2

    G3

    Example of Bilateral Contract Trading … Continued

    • G1 injects 300 MW at Bus 1 to meet 300 MW demand of L1 at Bus 2

    • G2 injects 200 MW at Bus 1 to meet 200 MW demand of L2 at Bus 2

    • If line capacity 1→2 ≥ 500 MW No problem

    • If line capacity 1→2 < 500 MW Some of these transactions

    might have to be curtailed

    Line Corridor 1→2

  • © 2011 D. Kirschen and L. Tesfatsion 17

    But curtail which one? Generation? Demand?

    • Could use administrative Transmission Load Relief

    (TLR) procedures

    • TLR procedures consider:

    Firm (i.e., guaranteed) vs. non-firm transactions

    Order in which transactions were registered

    Historical considerations

    • TLR procedures do not consider economic benefits

    • Economically inefficient

    • How about letting the participants themselves decide?

  • © 2011 D. Kirschen and L. Tesfatsion 18

    But curtail which one … Continued

    • What about permitting traders to purchase a “right” to use

    the network when arranging a bilateral energy trade?

    Transmission right = Right to use the transmission system for a particular transaction.

    Physical transmission right = Transmission right meant to support the actual transmission of a certain amount of power

    over a given transmission line.

    Financial transmission right (FTR) for 1 MWh from Bus k to Bus m in Hour H = Entitlement to receive a payment

    ($/h) equal to [LMPm - LMPk]•1 MWh in hour H if positive, or an

    obligation to pay this amount to system operator in hour H if

    negative.

  • © 2011 D. Kirschen and L. Tesfatsion 19

    G1

    G2

    L1

    L2

    Bus 1 Bus 2

    G3

    Physical Transmission Rights:

    • G1 sells 300 MW to L1 at 30.00 $/MWh

    • G2 sells 200 MW to L2 at 32.00 $/MWh

    • G3 sells electric power at 35.00 $/MWh

    • L1 not willing to pay more than 5 $/MWh for transmission rights

    • L2 not willing to pay more than 3 $/MWh for transmission rights

    • Otherwise, why not buy from G3 instead ?

    Line Corridor 1→2

  • © 2011 D. Kirschen and L. Tesfatsion

    20

    Problems with Physical Transmission Rights

    • Not practical to implement (parallel path problem)

    Path that electric power takes through a network from any point

    of injection is determined by physical laws and not by the wishes

    of market traders

    • Market power

    The assignment of physical transmission rights can potentially

    increase GenCo market power despite parallel path problems

    Can permit strategic GenCos to create “load pockets” by limiting

    the flow of power into their region

  • © 2011 D. Kirschen and L. Tesfatsion 21

    Parallel Path Problem: Three-Bus Example

    1 2

    3

    GC GA

    GB

    GD

    LY

    LZ

    Branch Reactance

    [p.u.]

    Capacity

    [MW]

    1-2 0.2 126

    1-3 0.2 250

    2-3 0.1 130

    Important Technical Remark:

    The reactance on a branch

    k-m does not depend on the

    direction of the power flow on

    the branch, i.e., xkm = xmk .

  • © 2011 D. Kirschen and L. Tesfatsion 22

    Parallel Path Problem … 3-Bus Example: Continued

    1 2

    3

    GC GA

    GB

    GD

    LY

    LZ

    I

    II

    A 400 MW sale between GenCo GB and LSE LY would

    require the purchase of transmission rights on all three

    branches: 1-3 (direct route) and 1-2 & 2-3 (indirect route)!

  • © 2011 D. Kirschen and L. Tesfatsion 23

    400 MW transfer between GB and LY physically possible?

    1 2

    3

    GC GA

    GB

    GD

    LY

    LZ

    Branch Reactance

    [p.u.]

    Capacity

    [MW]

    1-2 0.2 126

    1-3 0.2 250

    2-3 0.1 130

    I

    II

    sum of reactances for 1-2 and 2-3

    Not possible -- exceeds capacities of branches 1-2 and 2-3 !

    F I

    0 . 2

    0.5 400 160 MW

    F II

    0 . 3

    0.5 400 240 MW

    sum of all reactances reactance for 1-3

    F12

    F23 F13

  • © 2011 D. Kirschen and L. Tesfatsion 24

    Counter-Flows: 200MW Trade Between GD & LZ

    1 2

    3

    GC GA

    GB

    GD

    LY

    LZ

    Flow effects resulting from

    the 200 MW transaction

    between GD and LZ III

    IV

    F III

    0 . 2

    200 80 MW

    0.5

    F IV

    0 . 3 200 120 MW

    0.5

    F12

    F23 F13

  • © 2011 D. Kirschen and L. Tesfatsion 25

    Resultant Flows:

    1 2

    3

    GC GA

    GB

    GD

    LY

    LZ

    F12 F23 FIFIII

    160 80 80 MW

    F13 FII

    FIV

    240 120 120 MW

    The resultant flows are now within the branch capacity limits.

    Branch Reactance

    [p.u.]

    Capacity

    [MW]

    1-2 0.2 126

    1-3 0.2 250

    2-3 0.1 130

    1→2→3 3→2→1

    1→3 3→1

  • © 2011 D. Kirschen and L. Tesfatsion 26

    Physical Transmission Rights and Parallel Paths

    • Counter-flows create additional physical transmission rights

    • Economic efficiency requires that these rights be considered

    • Decentralized trading:

    System operator only checks overall feasibility

    Participants trade physical transmission rights bilaterally

    Theory:

    • Enough participants market discovers optimum

    Practice:

    • Complexity and amount of information involved are such that it is

    unlikely that this optimum can be found in time to be practical.

  • © 2011 D. Kirschen and L. Tesfatsion 27

    Physical Transmission Rights and Market Power

    • G3 – the only GenCo at Bus B -- is in a potential “load pocket”

    • G3 purchases physical transmission rights from Bus A to Bus B

    • G3 does not use or resell these rights

    • Effectively reduces transmission capacity from Bus A to Bus B

    • Allows G3 to increase price at Bus B charged to L1 and L2

    • “Use them or loose them” provision for physical transmission rights is difficult to enforce in a timely manner

    G1

    G2

    L1

    L2

    Bus A Bus B

    G3

    Max Capacity = 250 MW

  • © 2011 D. Kirschen and L. Tesfatsion 28

    Centralized (Pool) Trading: Day-Ahead Market

    (Cf. K/S Section 6.3)

    • At the beginning of each day D, LSEs & GenCos submit demand bids & supply offers to an ISO for 24 successive hours of trade in the Day-Ahead Market for day D+1

    • ISO’s objective:

    Optimally balance generation and load for each hour, conditional on GenCo reported cost and operating capacity attributes

    Respect network transmission constraints (e.g., capacity limits)

    • No branch congestion/losses in hour H Uniform market price ($/MWh) determined across the grid for hour H

    • Branch congestion/losses in hour H LMPs for hour H will typically separate across the grid (different market prices set for LSEs/GenCos located at different buses)

  • © 20011 D. Kirschen and L. Tesfatsion 29

    DS= 1500 MW

    Syldavia

    DB= 500MW

    Borduria

    Example: Borduria-Syldavia Interconnection (Cf. K/S Section 6.3.1)

    • Assume installed generation capacity exceeds demand in each

    country by a significant margin.

    • Assume no congestion or losses within each country.

    • Assume a uniform market price for electric energy is set within

    each country at the LMP, i.e., at the least cost of servicing the

    next MW of power for that country while respecting all system

    constraints.

  • © 2011 D. Kirschen and L. Tesfatsion 30

    MW

    $/MWh

    13

    1500

    43

    DS= 1500 MW

    Syldavia

    DB= 500MW

    Borduria

    Borduria/Syldavia Demand and Supply Functions in the

    Absence of Any Interconnection Between the Countries

    10

    15

    500 MW

    $/MWh

    B MCB 10 0.01PB [$ / MWh] S MC S 13 0.02PS [$ / MWh]

    S MC S 13 0.02 1500 43 $ /MWh

    Borduria Market Price Syldavia Market Price

    S

    S D

    D

    B

    MC B 10 0 . 01 500 15 $ / MWh

  • © 2011 D. Kirschen and L. Tesfatsion 31

    DB= 500MW

    Borduria

    DS= 1500 MW

    Syldavia

    Now suppose an interconnection is built …

    • Without any interconnection, countries operate where

    MCB = 15 $/MWh < MCS = 43 $/MWh.

    • Economic effects of introducing an interconnection?

    Max Cap = 1600 MW

  • © 2011 D. Kirschen and L. Tesfatsion 32

    DB= 500MW

    Borduria

    DS= 1500 MW

    Syldavia

    Can Borduria GenCos meet entire demand?

    • GenCos in Borduria CAN meet entire demand (2000 MW).

    • But then Syldavia MCs = 13 $/MWh < Borduria MCB = 30 $/MWh.

    • Having Borduria generate all power is not economically tenable.

    • At operating point the load would want to buy from Syldavia instead.

    PB 2000MW

    PS 0MW

    MCB 30$ / MWh

    MCS 13$ / MWh

    Max Cap = 1600 MW

  • © 2011 D. Kirschen and L. Tesfatsion 33

    DB= 500MW

    Borduria

    DS= 1500 MW

    Syldavia

    Two-Country Market Clearing Point

    B S

    PB PS DB DS 500 1500 2000MW

    PB 1433MW

    PS 567MW

    B MCB 10 0.01PB [$ / MWh] S MC S 13 0.02PS [$ / MWh]

    Max Cap = 1600 MW

    B

    S

    24 . 30 $ / MWh

  • © 2011 D. Kirschen and L. Tesfatsion 34

    More precisely, 5 equations and 5 unknowns:

    • 5 Unknowns

    Market price π

    Bordurian price πB and generation level PB

    Syldavian price πS and generation level PS

    • 5 Equations in the 5 Unknowns (Can Solve by Substitution)

    1. π = πB

    2. π = πS

    3. PB + PS = 2000 MW

    4. πB = 10 + 0.01 PB [ i.e., Bordurian price = MCB(PB) ]

    5. πS = 13 + 0.02 PS [ i.e., Syldavian price = MCS(PS) ]

  • © 2011 D. Kirschen and L. Tesfatsion 35

    DB= 500MW

    Borduria

    DS= 1500 MW

    Syldavia

    Interconnection Flow at the Market Clearing Point

    PB 1433MW PS 567MW

    Flow of power across the

    interconnection B→S

    resulting power

    flow 933 MW →

    F BS P B D B D S P S 933 MW

    Max Cap = 1600 MW

  • © 2011 D. Kirschen and L. Tesfatsion 36

    Graphical Representation of Market Clearing

    = 567 MW

    24.30 $/MWh

    = 1433 MW

    = 2000 MW

    = 500 MW = 1500 MW

    24.30 $/MWh

    = 933 MW

    Supply curve for

    Syldavia

    Supply curve for

    Borduria

    D B D S

    D B D S

    Borduria Syldavia

    F BS

    S

    MC S

    B

    MC B 30.00 $/MWh

    13.00 $/MWh

    P B

    P S

  • © 2011 D. Kirschen and L. Tesfatsion 37

    Constrained Transmission

    • What if the interconnection can only carry 400 MW?

    • Borduria must cut back its production so that at most 400 MW is

    transmitted to Syldavia over the interconnection:

    PB = 500 MW + 400 MW = 900 MW

    • Syldavia must then generate any additional power needed to

    meet it’s local load of 1500 MW:

    PS = 1500 MW - 400 MW = 1100 MW

    • Results in persistent market price differential of [35 – 19] $/MWh

    = 16 $/MWh between Borduria and Syldavia

    • The two bus prices 19 $/MWh and 35 $/MWh are called

    Locational marginal Prices (LMP) or nodal prices

    B MCB 10 0.01 900 19 $ / MWh

    S MC S 13 0.02 1100 35 $ /MWh

    Market prices in

    Borduria & Syldavia

    set by price = MC

  • © 2011 D. Kirschen and L. Tesfatsion 38

    = 1100 MW

    35 $/MWh

    = 900 MW

    = 2000 MW

    = 500 MW = 1500 MW

    = 400 MW

    16 $/MWh

    Graphical Representation: Constrained Transmission

    S MC S

    B MC B

    FBS

    D B D S

    D B D S

    Borduria Syldavia

    13 $/MWh

    19 $/MWh

    P B

    P S

  • © 2011 D. Kirschen and L. Tesfatsion 39

    Separate markets Single market Single market

    with conges tion

    PB [MW] 500 1,433 900

    B [$/MWh] 15 24.33 19

    RB [$/h] 7,500 34,865 17,100

    E B [$/h] 7,500 12,165 9,500

    PS [MW] 1500 567 1100

    S [$/MWh] 43 24.33 35

    RS [$/h] 64,500 13,795 38,500

    E S [$/h] 64,500 36,495 52,500

    FBS [MW] 0 933 400

    RTOTAL R B R S 72,000 48,660 55,600

    E TOTAL E B E S 72,000 48,660 62,000

    Summary of Results for Borduria and Syldavia

    500 1433 900

    15.00 24.33 19

    7,500 34,865 17,100 (buyer payments)

    (seller revenues)

    (seller revenues)

    (buyer payments)

    7,500 12,165 9,500

    1,500 567 1,100

    43 24.33 35

    64,500 13,795 38,500

    64,500 36,495 52,500

    0 933 400

    72,000 48,660 55,600

    72,000 48,600 62,000

    Why? Why?

  • © 2011 D. Kirschen and L. Tesfatsion 40

    Winners and Losers from Interconnection

    • Winners? (Knowing prices, payments, revenues not enough! )

    Economies of both countries? (Is total net surplus higher?)

    Bordurian power sellers ? (Is their total net seller surplus higher?)

    Syldavian power buyers ? (Is their total net buyer surplus higher?)

    • Losers ?

    Bordurian power buyers ? (Is their total net buyer surplus lower?)

    Syldavian power sellers ? (Is their total net seller surplus lower?)

    • Congestion in the interconnection reduces its benefits

  • © 2011 D. Kirschen and L. Tesfatsion 41

    ISO Net Surplus (or “Congestion Rent”) Defined: (K/S also refer to it as “merchandizing surplus”, Section 6.3.1.3)

    E TOTAL B DB S D S

    RTOTAL B PB S PS B (DB FBS ) S (DS FBS )

    E TOTAL RTOTAL S D S B DB S PS B PB

    S (D S PS ) B (DB PB )

    S FBS B ( FBS )

    ( S B ) FBS

    Buyers’ Payments (Collected by ISO):

    Sellers’ Revenues (Paid by ISO):

    ISO Net Surplus (“Congestion Rent”):

    Tends to increase as

    price separation increases

  • © 2011 D. Kirschen and L. Tesfatsion 42

    ISO Net Surplus as a Function of Power Flow FBS

    0

    10000

    20000

    30000

    40000

    50000

    60000

    70000

    80000

    0 100 200 300 400 500 600 700 800 900 1000

    Flow on the Interconnection [MW]

    Pa

    ym

    ents

    and

    Re

    ve

    nu

    es [

    $/h

    ]

    Buyers’ payments

    Sellers’ revenues

  • © 2011 D. Kirschen and L. Tesfatsion 43

    ISO Net Surplus Summary:

    • ISO Net Surplus = [buyer payments – seller revenues]

    collected in pool-based trading by an independent

    system operator (ISO), who also acts as market operator

    • ISO Net Surplus is typically positive in loss-less systems

    whenever one or more lines are congested.

    • Should not be kept by the ISO because it gives a

    perverse incentive to the ISO to maintain or encourage

    grid congestion.

    • Should not be returned directly to buyers and sellers

    users because that would blunt the economic incentive

    provided by nodal pricing.

  • © 2011 D. Kirschen and L. Tesfatsion 44

    Pool Trading in Three-Bus Examples: K/S 6.3.2

    1 2

    3

    C A

    B

    D

    50 MW 60 MW

    300 MW

    Assumptions:

    Constant branch capacity limits restricting branch power flows. Resistance of lines is negligible and can be ignored.

    GenCos report their true capacities and marginal costs to the ISO.

    The load’s 60MW demand is fixed (price insensitive).

    = GenCo

    = Load

  • © 2011 D. Kirschen and L. Tesfatsion 45

    Pool Trading in a Three-Bus Systems: K/S 6.3.2

    1 2

    3

    C A

    B

    D

    Branch Reactance

    [p.u.]

    Capacity

    [MW]

    1-2 0.2 126

    1-3 0.2 250

    2-3 0.1 13050 MW 60 MW

    300 MW

    Gene r a t o r C a pac it y

    [ MW ]

    M a rg i na l Co s t

    [ $ / MWh ]

    A 140 7 .50

    B 285 6.00

    C 90 14.00

    D 85 10.00

    Note: Each GenCo assumed

    to have a constant MC!

  • © 2011 D. Kirschen and L. Tesfatsion 46

    Problem Solution: First consider simple economic dispatch with no binding transmission constraints

    1 2

    3

    C A

    B

    D

    50 MW 60 MW

    300 MW

    125 MW

    285 MW

    0 MW

    0 MW

    F13 F23

    F12

    Total Demand (MW): 410 = [50+300+60]

    Cheapest GenCos B&A:

    Capacity B = 285 MW Capacity A = 140 MW

    Tot Cap A+B = 425 MW

    = enough to meet demand

    PB=285 (max), PA=125 Are branch limits exceeded?

  • © 2011 D. Kirschen and L. Tesfatsion 47

    Solving for Power Flows: What Do We Know?

    • Overall power balance equation (assuming no branch losses):

    (a) P1 + P2 + P3 = Demand1 + Demand2 + Demand3

    • Power balance equation for each bus, from KCL (assuming no branch losses) – compare K/S (6.30)-(6.32)

    (b) P1 – Demand1 = F12 + F13

    (c) P2 – Demand2 = F21 + F23 (Note: F21 = - F12)

    (d) P3 – Demand3 = F31 + F32 (Note: F31 = - F13 , and F32 = - F23)

    NOTE: Summing equations (b), (c), and (d) gives equation (a), so any one of these 4 equations can be omitted with any loss of information.

  • © 2011 D. Kirschen and L. Tesfatsion 48

    Solving for Power Flows for Given Power Injections & Demands:

    • Overall power balance equation (assuming no branch losses):

    (a) P1 + P2 + P3 = Demand1 + Demand2 + Demand3

    410 + 0 + 0 = 50 + 60 + 300 (MW)

    • Power balance equation for each bus, from KCL (assuming no branch losses) – compare K/S (6.30)-(6.32)

    (b) P1 – Demand1 = 360 = F12 + F13

    (c) P2 – Demand2 = - 60 = F21 + F23 (Note: F21 = - F12)

    (d) P3 – Demand3 = -300 = F31 + F32 (Note: F31 = - F13 , and F32 = - F23)

    NOTE: Since equations (a), (b), and (c) sum to zero, we have only two

    independent equations in 3 unknowns {F12,F13,F23} ! We need an

    additional equation.

  • © 2011 D. Kirschen and L. Tesfatsion 49

    1

    60 MW

    2

    3

    300 MW

    360 MW

    1

    60 MW

    2

    3

    60 MW

    1 2

    3

    300 MW

    300 MW

    Calculating Branch Flows Via Superposition:

    F12 = F1A + F2

    A = 120+36 = 156 MW

    F1

    A = PTDF1A • 300

    = ( x13/ [ ∑ xiJ ] ) • 300

    = (.2/.5) • 300 = 120MW F2

    A = PTDF2A• 60 = ([x13+x32]/ [ ∑ xiJ ]) • 60MW

    = ([.2+.1])/.5) • 60 = 36MW

    Additional equation from linear

    superposition (Fkm = PTDFkm• Pk):

    F12

    F1A

    F2A

  • © 2011 D. Kirschen and L. Tesfatsion 50

    Solving for Power Flows for Given Bus Power Injections & Demands:

    • Overall Power Balance Equation (Assuming No Losses):

    (a) P1 + P2 + P3 = Demand1 + Demand2 + Demand3

    410 + 0 + 0 = 50 + 60 + 300 (MW)

    • Power Balance Equation for Each Bus from KCL (Assuming No Losses) – compare K/S (6.30)-(6.32)

    (b) P1 – Demand1 = 360 = F12 + F13

    (c) P2 – Demand2 = - 60 = F21 + F23 (Note: F21 = - F12)

    (d) P3 – Demand3 = -300 = F31 + F32 (Note: F31 = - F13 , and F32 = - F23)

    • Implied Solution: Compare K/S, (6.45)-(6.47), p. 160

    F12 = 156 MW (superposition) Violates branch 1-2 Max Cap = 126!

    F13 = [ 360 – 156 ] = 204 MW from equation (b)

    F23 = [ F31 + 300 ] = [ - 204 + 300 ] = 96 MW from equation (d)

  • © 2011 D. Kirschen and L. Tesfatsion 51

    Power Flows Implied by Simple Economic Dispatch

    1 2

    3

    C A

    B

    D

    50 MW 60 MW

    300 MW

    125 MW

    285 MW

    0 MW

    0 MW

    156 MW

    96 MW 204 MW

  • © 2011 D. Kirschen and L. Tesfatsion 52

    Branch 1-2 is overloaded by 30 MW !

    1 2

    3

    C A

    B

    D

    50 MW 60 MW

    300 MW

    125 MW

    285 MW

    0 MW

    0 MW

    156 MW

    96 MW 204 MW

    Max Cap = 126 MW

  • © 2011 D. Kirschen and L. Tesfatsion 53

    1

    1 MW

    2

    3

    1 MW

    Correcting the Economic Dispatch: Find the Least-Cost Dispatch

    Modification that will Remove the Branch Overload

    Suppose 1 MW of power is injected at Bus 2

    and withdrawn at Bus 1:

    0.6 MW

    0.4 MW

    PTDF21 = [ X13 + x23]/[x12 + x13 + x23] = 0.3/0.5 = 0.6

    PTDF231 = [ X12]/[x12 + x13 + x23] = 0.2/0.5 = 0.4

    X12 = 0.2

    X13 = 0.2 X23 = 0.1

  • © 2011 D. Kirschen and L. Tesfatsion 54

    1

    60

    MW

    2

    3

    300 MW

    360 MW 156 MW

    204 MW 96 MW

    1

    50 MW

    2

    3

    50 MW 30 MW

    20MW

    1

    10

    MW

    2

    3

    300 MW

    310 MW 126 MW

    184 MW 116 MW

    Suppose 50 MW of power is injected at Bus 2 and

    withdrawn at Bus 1: Solve by Superposition

    NOTE: 50 x 0.6 = 30

    50 x 0.4 = 20

    This gives a feasible solution !

  • © 2011 D. Kirschen and L. Tesfatsion 55

    1

    1 MW

    2

    3

    1 MW

    Alternatively, suppose 1 MW of power is injected at

    Bus 3 and withdrawn at Bus 1:

    0.6 MW

    0.4 MW

    X12 = 0.2

    X13 = 0.2 X23 = 0.1

    PTDF321 = 0.2/0.5 = 0.4 PTDF31 = 0.3/0.5 = 0.6

  • © 2011 D. Kirschen and L. Tesfatsion 56

    1

    60 MW

    2

    3

    300 MW

    360 MW 156 MW

    204 MW 96 MW

    1

    60 MW

    2

    3

    225 MW

    285 MW 126 MW

    159 MW 66 MW

    1

    75 MW

    2

    3

    75 MW

    30 MW

    45 MW

    Suppose an additional 75 MW of power is injected at

    Bus 3 and withdrawn at Bus 1: Solve by Superposition

    NOTE: 75 x 0.4 = 30

    75 x 0.6 = 45

    This also is a feasible solution !

  • © 2011 D. Kirschen and L. Tesfatsion 57

    Compare the Total Avoidable Costs Resulting from

    Each of the Two Feasible Solutions: (See K/S Section 6.3.2.2) Non-Feasible Simple Economic Dispatch (NOTE: Here MC is constant):

    MCA•125 + MCB•285 + MCC•0 + MCD•0 = 937.50 + 710 = 2,647.50 $/h

    Cost with [-50] MW Redispatch of GenCo A at Bus 1 and [+50] MW Redispatch of GenCo C at Bus 2:

    MCA•75 + MCB•285 + MCC•50 = 562.50 + 1710 + 700 = 2,972.50 $/h

    Cost with [-75] MW Redispatch of GenCo A at Bus 1 and a [+75] MW Redispatch of GenCo D at Bus 3:

    MCA•50 + MCB•285 + MCD•75 = 375 + 1710 + 750 = 2,835.00 $/h

    Least-Cost Solution: Redispatch GenCos A & D

    Cost of Having a Max Cap of 126 MW on Line 1-2:

    [ 2,835.00 – 2,647.50 ] = 187.50 $/h

  • © 2011 D. Kirschen and L. Tesfatsion 58

    Transmission-Constrained Economic Dispatch Solution:

    1 2

    3

    C A

    B

    D

    50 MW 60 MW

    50 MW

    285 MW

    0 MW

    75 MW

    126 MW

    66 MW 159 MW

    300 MW

  • © 2011 D. Kirschen and L. Tesfatsion 59

    Locational Marginal Prices (Nodal Prices):

    • Least cost of servicing an additional MW of demand at a

    particular bus at a particular time subject to system

    feasibility constraints

    • These system feasibility constraints typically include

    power balance constraints, GenCo operating capacity

    constraints, and branch capacity constraints.

    • K/S refer to branch capacity constraints as security

    constraints.

  • © 2011 D. Kirschen and L. Tesfatsion 60

    LMP: Transmission-Constrained Economic Dispatch

    3-Bus Example

    • Bus 1 LMP?

    • GenCo A at Bus 1 is not at its max or min operating capacity limits (i.e., Genco A is marginal )

    • GenCo A cheaper than all other GenCos not capacity constrained

    • Therefore the LMP at Bus 1 is

    1 2

    3

    C A

    B

    D

    50 MW 60 MW

    50 MW

    285 MW

    0 MW

    75 MW

    126

    MW

    66 MW

    159 MW

    300 MW

    1

    MC A 7 . 50 $ / MWh

  • © 2011 D. Kirschen and L. Tesfatsion 61

    LMP Example … Continued

    • Bus 3 LMP ?

    • Genco A cheaper than GenCo D

    • But increasing A’s generation would overload line 1-2

    • D cheaper than C and is at Bus 3. Therefore LMP at Bus 3 is

    1 2

    3

    C A

    B

    D

    50 MW

    60 MW

    50 MW

    285 MW

    0 MW

    75 MW

    126

    MW

    66 MW

    159 MW

    300 MW

    3

    MC D 10 $ / MWh

  • © 2011 D. Kirschen and L. Tesfatsion 62

    LMP Example … Continued

    • Bus 2 LMP ?

    • GenCo C is very expensive

    • GenCo B is already at its max operating capacity

    • Need to consider redispatch of marginal GenCos A and/or D

    1 2

    3

    C A

    B

    D

    50 MW 60 MW

    50 MW

    285 MW

    0 MW

    75 MW

    126

    MW

    66 MW

    159 MW

    300 MW

  • © 2011 D. Kirschen and L. Tesfatsion 63

    1

    1 MW

    2

    3

    1 MW

    0.6 MW

    0.4 MW

    1

    1 MW

    2

    3

    0.2 MW

    0.8 MW

    1 MW

    Bus 2 LMP ... Continued

    Option 1: 1 MW for Bus 2

    generated at Bus 1

    by GenCo A

    Option 2: 1 MW for Bus 2

    generated at Bus 3

    by GenCo D

    Either option increases flow on branch 1-2, violating max cap!

    Branch Reactance

    [p.u.]

    Capacity

    [MW]

    1-2 0.2 126

    1-3 0.2 250

    2-3 0.1 130

  • © 2011 D. Kirschen and L. Tesfatsion 64

    LMP at Bus 2 … Continued

    Consider simultaneous redispatch ΔP1 for GenCo A at

    Bus 1 and ΔP3 for GenCo D at Bus 3 satisfying:

    (i) ΔP1 + ΔP3 = ΔP2 ; (ii) ΔP2 = 1 MW ; (iii) ΔF12 = 0 ;

    (iv) ΔF12 = [ 0.6 • ΔP1 + 0.2 • ΔP3 ]

    1

    1 MW

    2

    3

    F 12 = 0

    P 3

    P 1

    P 2

  • © 2011 D. Kirschen and L. Tesfatsion 65

    1

    1 MW

    2

    3

    0.2 MW

    0.8 MW

    1 MW

    LMP Solution Using Superposition

    P1 P3 P2 1 MW

    0.6 P1 0.2 P3 F12 0 MW

    P1 0.5 MW

    P3 1.5 MW

    2 1.5 MCD 0.5 MCA 11.25 $ / MWh

    1

    1 MW

    2

    3

    1 MW

    0.6 MW

    0.4 MW

    (10 $/MWh =

    LMP at Bus 1)

    (7.50 $/MWh =

    LMP at Bus 3)

    (LMP at

    Bus 2)

  • © 2011 D. Kirschen and L. Tesfatsion 66

    General LMP Observations

    A GenCo is marginal if it is operating at a point p that lies

    strictly between its min and max operating capacity limits:

    CapL < p < CapU.

    The LMP at a bus where there is a marginal GenCo is determined by this GenCo’s marginal cost.

    The LMP at a bus with no marginal GenCo is a linear

    combination of the marginal costs of the marginal GenCos.

    For the 3-bus example at hand:

    GenCos A and D are marginal.

    GenCos B & C are not marginal, but for different reasons.

    • GenCo B is an inframarginal unit at its max operating capacity

    • GenCo C is an extramarginal unit too expensive to dispatch, hence it

    is at its min operating capacity limit of 0

  • © 2011 D. Kirschen and L. Tesfatsion 67

    Outcome Summary for the Three-Bus Example:

    Bu s 1 Bu s 2 Bu s 3 S ys t e m

    Purchase [ MW ] 50 60 300 410

    P rodu ct ion [ MW ] 335 0 75 410

    LMP [ $ / M W h] 7 . 50 11 . 25 10 . 00 -

    [ $ / h] 375 . 00 675 . 00 3 , 000.00 4 , 050.00

    P rodu c er r evenue s [ $ / h] 2 , 512.50 0 . 00 750 . 00 3 , 262.50

    [ $ / h ] 787 . 50

    Buyer Payments

    ISO Net Surplus*

    * In ISO-managed wholesale power markets, ISO Net Surplus is also known as “congestion rent” or “merchandising surplus.” The latter usage is

    somewhat confusing for those familiar with standard economic terminology

    since it represents an amount of money collected by the market operator and

    not by “merchants” (private traders).

  • © 2011 D. Kirschen and L. Tesfatsion 68

    Power flows from

    high price to low

    price!

    Counter-Intuitive Flows: (cf. K/S 6.3.2.5)

    1 2

    3

    C A

    B

    D

    50 MW 60 MW

    50 MW

    285 MW

    0 MW

    75 MW

    126 MW

    66 MW 159 MW

    300 MW

    2=11.25 $/MWh

    3=10.00 $/MWh

    3=7.50 $/MWh

  • © 2011 D. Kirschen and L. Tesfatsion 69

    Counter-Intuitive Prices and Strategic Behavior:

    • LMPs at buses without a marginal GenCo can be higher or

    lower than LMPs at other buses.

    • LMPs at such buses can even be negative!

    • Determining LMPs in realistically sized and structured

    grids is complicated, requiring specialized software.

    • Strategically placed GenCos can exert control over LMPs

    through appropriate choice of reported supply offers.

    • Transmission branch congestion can facilitate the exercise

    of GenCo market power (i.e., it can help GenCos exert

    advantageous control over LMPs).

  • © 2011 D. Kirschen and L. Tesfatsion

    70

    General Method for Computing LMPs:

    Bid/Offer-Based Optimal Power Flow (OPF)

    • Optimisation Problem:

    Objective: Maximisation of Total Net Surplus

    Constraints: Power flow equations (bus balance constraints, generation operating capacity limits, branch capacity limits, …)

    The Lagrange multiplier associated with the bus power balance constraint for any bus k gives the LMP solution value for this bus.

    OPF solution often approximated by a “dc optimal power flow” formulation (discussed in K/S Section 6.3.4.4)

    • In North American Restructured Wholesale Power Markets:

    Hourly LMPs are typically calculated each operating day D for all 24 hours of a “day-ahead market” on day D+1

    LMPs are also calculated/updated every five minutes for a “real-time market” (managed spot market) that remains open throughout day D.