A. Mathematical Formulation

3 downloads 0 Views 415KB Size Report
Generation capacity for offshore wind plant g under scenario node (MW). .... Offshore generation investment co-optimized with transmissionΒ ...
A.

Mathematical Formulation

In this section the mathematical formulations for the perfect information and min-max regret planning problem are given. The mathematical nomenclature is shown in Section A.1. We then introduce the perfect information problem formulation, where a single scenario is considered and proceed with the min-max regret formulation expressed in terms of the perfect information problem solutions. Note that all problems presented in the paper were programmed using the commercial software (FICO Xpress, 2009) and solved using the in-built branch-and-bound algorithm with a stopping criterion of finding an integer solution within 0.1% of the best available bound.

A.1 Nomenclature Sets and indices Note that cardinality (number of elements) of a set Ξ© is denoted |Ξ©|. e(m) Epoch e associated with scenario tree node m (an epoch is a set of years) Set of countries in EU, indexed c. |Ω𝐢 | = 8. Ω𝐢 Set of epochs, indexed 𝑒. |Ω𝐢 | = 5. Ω𝐸 Set of onshore generators, indexed 𝑔. |Ω𝑔 | = 256. Ω𝐺 Set of all transmission corridors, indexed 𝑙. |Ω𝐿 | = 5. Ω𝐿 Ω𝐿𝐴𝐢

𝐴𝐢

Set of onshore AC transmission corridors, indexed 𝑙. |Ω𝐿 | = 28. 𝐷𝐢

Ω𝐷𝐢 Set of offshore DC transmission corridors, indexed 𝑙. |Ω𝐿 | = 106. 𝐿 Set of scenario tree nodes, indexed π‘š. For a single scenario, |Ω𝑀 | = 5. Ω𝑀 Set of system buses, indexed 𝑛. |Ω𝑁 | = 16. Ω𝑁 Set of scenarios, indexed 𝑠. Ω𝑆 Set of demand periods, indexed 𝑑. |Ξ© 𝑇 | = 200. Ω𝑇 Set of offshore wind clusters, indexed 𝑀. |Ξ©π‘Š | = 32. Ω𝑀 Set that contains all parent nodes of node π‘š, including π‘š Ξ¦π‘š Input parameters Susceptance of AC transmission corridor 𝑙 (p.u.) 𝑏𝑙 𝐷𝑑,𝑛 Demand at bus 𝑛 in period 𝑑 (MW) 𝐼𝑛,𝑙 Bus-to-corridor incidence matrix 0 Initial capacity of transmission corridor 𝑙 (MW) 𝐹𝑙 Maximum capacity per DC cable that can be added to corridor 𝑙 (MW) 𝐹̅𝑙 𝐡 Annuitized fixed cost of reinforcing corridor 𝑙 (€/year) πœ…π‘™ Μ… Annuitized fixed cost of a DC undersea cable used in corridor 𝑙 (€/year) πœ…π‘™π΅ 𝐹 Annuitized variable cost of reinforcing corridor 𝑙 (€/MW-year) πœ…π‘™ πœ…π‘”πΊ Operating cost of generating unit 𝑔 (€/MWh) Cost of offshore wind curtailment (€/MWh) πœ…π‘€ 𝑝̅𝑒,𝑔 Generation capacity for onshore plant g at epoch 𝑒 (MW) 𝐼 Cumulative discount factor for investment cost in epoch e (p.u.) π‘Ÿπ‘’ Cumulative discount factor for operation cost in epoch e (p.u.) π‘Ÿπ‘’π‘‚ Start node for corridor l 𝑒𝑙 Destination node for corridor l 𝑣𝑙 Wind availability factor at period 𝑑 (p.u.) π‘Šπ‘‘ Μ… πœ‰π‘š,𝑀 Generation capacity for offshore wind plant g under scenario node π‘š (MW) Time duration of period 𝑑 (hours) πœπ‘‘ Decision variables Binary variable signifying the choice of reinforcement scenario node π‘š for AC π΅π‘š,𝑙 corridor 𝑙 π΅Μ…π‘š,𝑙 Integer variable signifying the number of DC undersea cables used in reinforcement 1

scenario node π‘š for DC corridor 𝑙 Investment cost under scenario node π‘š (€) Operational cost under scenario node π‘š (€) Capacity added to corridor l in scenario node π‘š (MW) State variable of aggregate extra capacity available to corridor l in scenario node π‘š (MW) Power flow on corridor 𝑙 at operating point (π‘š, 𝑑) (MW) Output of onshore generation unit 𝑔 at operating point (π‘š, 𝑑) (MW) Bus angle of 𝑛 at operating point (π‘š, 𝑑) (rad) Output of offshore cluster 𝑀 at operating point (π‘š, 𝑑) (MW) Curtailed output of offshore cluster 𝑀 at operating point (π‘š, 𝑑) (MW)

𝐼 πΆπ‘š 𝑂 πΆπ‘š πΉπ‘š,𝑙

πΉΜƒπ‘š,𝑙 π‘“π‘š,𝑑,𝑙 π‘π‘š,𝑑,𝑔 πœƒπ‘š,𝑑,𝑛 πœ‰π‘š,𝑑,𝑀 𝐢 πœ‰π‘š,𝑑,𝑀

A.2 Mathematical Formulation – Perfect Information The perfect information problem for scenario 𝑠, comprised of all scenario tree nodes π‘š ∈ 𝑠, is a Mixed Integer-Linear Program (MILP) defined as follows. 𝑂 𝐼 𝐼 𝑂 Ο‰s = 𝐦𝐒𝐧 { βˆ‘ (π‘Ÿπ‘’(π‘š) π‘π‘š + π‘Ÿπ‘’(π‘š) π‘π‘š )}

(1)

βˆ€π‘šβˆˆπ‘ 

where 𝐼 π‘π‘š = βˆ‘ π΅π‘š,𝑙 πœ…π‘™π΅ + βˆ‘ π΅Μ…π‘š,𝑙 πœ…π‘™π΅Μ… + βˆ‘ πΉπ‘š,𝑙 πœ…π‘™πΉ π‘™βˆˆΞ©π΄πΆ 𝐿

π‘™βˆˆΞ©π·πΆ 𝐿

𝐢 𝑂 π‘π‘š = βˆ‘ πœπ‘‘ (βˆ‘ π‘π‘š,𝑑,𝑔 πœ…π‘”πΊ + βˆ‘ πœ‰π‘š,𝑑,𝑀 πœ…π‘€ ) βˆ€π‘‘

βˆ€g

,βˆ€π‘š ∈ 𝑠

(2)

,βˆ€π‘š ∈ 𝑠

(3)

π‘™βˆˆΞ©πΏ

βˆ€π‘€

Investment constraints ∈ Ω𝐿𝐴𝐢 ∈ Ω𝐷𝐢 𝐿 ∈ Ω𝐿𝐴𝐢 ∈ Ω𝐷𝐢 𝐿

π΅π‘š,𝑙 ∈ {0,1} π΅Μ…π‘š,𝑙 ∈ β„€ πΉπ‘š,𝑙 ≀ π΅π‘š,𝑙 𝑀 πΉπ‘š,𝑙 ≀ π΅Μ…π‘š,𝑙 𝐹𝑙

,βˆ€π‘š ∈ 𝑠, 𝑙 ,βˆ€π‘š ∈ 𝑠, 𝑙 ,βˆ€π‘š ∈ 𝑠, 𝑙 ,βˆ€π‘š ∈ 𝑠, 𝑙

πΉΜƒπ‘š,𝑙 = βˆ‘ 𝐹𝑗,𝑙

βˆ€π‘š ∈ 𝑠, 𝑙 ∈ Ω𝐿

(8)

0 ≀ π‘π‘š,𝑑,𝑔 ≀ π‘Μ…πœ€(π‘š),𝑔 Μ… π‘Šπ‘‘ 0 ≀ πœ‰π‘š,𝑑,𝑀 ≀ πœ‰π‘š,𝑀 𝐢 Μ… π‘Šπ‘‘ πœ‰π‘š,𝑑,𝑔 + πœ‰π‘š,𝑑,𝑀 = πœ‰π‘š,𝑀 π‘“π‘š,𝑑,𝑙 = 𝑏𝑙 (πœƒπ‘š,𝑑,𝑒𝑙 βˆ’ πœƒπ‘š,𝑑,𝑣𝑙 ) βˆ’(𝐹𝑙0 + πΉΜƒπ‘š,𝑙 ) ≀ π‘“π‘š,𝑑,𝑙 ≀ 𝐹𝑙0 + πΉΜƒπ‘š,𝑙

,βˆ€π‘š ∈ 𝑠, 𝑑 ,βˆ€π‘š ∈ 𝑠, 𝑑 ,βˆ€π‘š ∈ 𝑠, 𝑑 ,βˆ€π‘š ∈ 𝑠, 𝑑 ,βˆ€π‘š ∈ 𝑠, 𝑑

(9) (10) (11) (12) (13)

βˆ‘ π‘π‘š,𝑑,𝑔 + βˆ‘ πœ‰π‘š,𝑑,𝑀 + βˆ‘ 𝐼𝑛,𝑙 π‘“π‘š,𝑑,𝑙 = 𝐷𝑑,𝑛

,βˆ€π‘š ∈ 𝑠, 𝑑 ∈ Ω𝑇 , 𝑛 ∈ Ω𝑁

(14)

,βˆ€π‘š ∈ 𝑠, βˆ€π‘ ∈ Ω𝐢

(15)

π‘—βˆˆΞ¦π‘š

(4) (5) (6) (7)

Operation constraints

βˆ€π‘”βˆˆπ‘›

βˆ€π‘€βˆˆπ‘›

βˆ€π‘™

βˆ‘ πœπ‘‘ ( βˆ‘ π‘π‘š,𝑑,𝑔 + βˆ‘ πœ‰π‘š,𝑑,𝑀 βˆ’ βˆ‘ 𝐷𝑑,𝑛 ) = 0 βˆ€π‘‘

βˆ€π‘”βˆˆπ‘

βˆ€π‘€βˆˆπ‘

∈ Ω𝑇 , 𝑔 ∈ Ω𝐺 ∈ Ω𝑇 , 𝑀 ∈ Ξ©π‘Š ∈ Ω𝑇 , 𝑀 ∈ Ξ©π‘Š ∈ Ω𝑇 , 𝑙 ∈ Ξ©AC 𝐿 ∈ Ω𝑇 , 𝑙 ∈ Ω𝐿

βˆ€π‘›βˆˆπ‘

The objective function (1) minimizes the discounted system cost over the study horizon, defined as the sum of investment and operation costs accrued over all scenario tree nodes. In all presented studies, a discount rate of 5% has been assumed, which corresponds to the Transmission Owner companies’ weighted average cost of capital (WACC) and this is aligned with (National Grid, 2017). 2

Equation (2) defines the investment cost corresponding to scenario tree node π‘š as the summation of capital and fixed costs across all line reinforcements. Equation (3) defines the operation cost for scenario tree node π‘š as the summation of generation costs and curtailed wind production. Constraint (4) defines decision variables π΅π‘š,𝑙 , signifying the reinforcement of an onshore-to-onshore corridor, as binary. Constraint (5) defines decision variables π΅Μ…π‘š,𝑙 , signifying the number of times fixed cost has been paid for the reinforcement of a DC offshore corridor, as integer. Constraint (6) applies to AC corridors, where M is a large positive number enabling planners to reinforce by a large amount while subject to only one fixed cost payment. Constraint (7) applies to all subsea corridors and states that the amount of transmission capacity that can be built at node π‘š is bounded by the times fixed cost has been paid. Particularly in the case of offshore-offshore corridors this reflects the fact that individual offshore plants proceed with their connection-to-shore plans while coordinating with neighboring projects of the same cluster in order to reduce fixed cost elements such as permitting and use of cablelaying vessels (𝐹̅𝑙 being very large). However, under the β€˜Radial’ policy choice, 𝐹̅𝑙 for offshoreonshore corridors is equal to 500 MW (rather than a very large number, which is used in the other policies so as to pay the fixed cost only once per epoch, if there is investment in that epoch) meaning that corridors are reinforced in lumps of 500 MW, essentially capturing the effects of offshore noncoordination in terms of increased fixed costs. Constraint (8) defines the investment state variables πΉΜƒπ‘š,𝑙 as the summation of the corresponding control variables 𝐹𝑗,𝑙 . Note that commissioning delays are not considered in the present model and reinforcements are assumed to materialize instantaneously. Operation constant (9) provides limits on dispatch of conventional generation. In a similar vein, operation constraint (10) limits output of offshore wind clusters according to the period-variable wind factor π‘Šπ‘‘ . Equation (11) computes the amount of curtailed offshore wind as the difference between available and dispatched power. Equation (12) describes how power is distributed over all onshore links according to the dc, linear power flow formulation. Note that power flows over cross-border, onshore-offshore and offshore-offshore links are assumed to be fully controllable HVDC lines and not depend on bus angles and line reactance. Due to the complexity of the presented model, the susceptance value 𝑏𝑙 in (12) does not change with network reinforcements and this is suggested to be further developed. For this purpose, the disjunctive approach explained in (Binato et al., 2001) arises as a particularly interesting method to be explored and potentially applied within the presented model. Constraint (13) bounds line flows according to the existing and commissioned transmission capacity of the corridor. The system balance equation (14) states that at each bus, the local generation from conventional technologies and offshore wind and net incoming/outgoing power flows satisfy local demand. Equation (15) is activated solely for the energy neutral policy choices. It enforces energy neutrality at the member-state level, ensuring that year-round energy produced by domestic onshore and offshore generators is equal the country’s energy demand. Table 1 shows the changes in the above MILP model in each policy choice.

3

Policy choice

Model features

Radial

(1)-(15) Offshore-to-offshore investment prevented 𝐹𝑙 = 500 MW (1)-(15) Offshore-to-offshore investment prevented 𝐹𝑙 β†’large number (1)-(15) Offshore-to-offshore investment allowed 𝐹𝑙 β†’large number (1)-(14) [(15) relaxed] Offshore-to-offshore investment allowed 𝐹𝑙 β†’large number (1)-(14) [(15) relaxed] Offshore generation investment co-optimized with transmission investment Extra constraints to ensure that net installed capacity of offshore wind per epoch corresponds to those in S1-S4 Offshore-to-offshore investment allowed 𝐹𝑙 β†’large number

Hub

Energy-neutral integrated

Fully integrated

Proactive

Table 1: Features to model each policy choice

We indicatively mention that the MILP planning problem when applied to a single scenario spanning five epochs has 670 integer variables, over 470,000 continuous variables and over 600,000 constraints.

A.3 Mathematical Formulation – Minimax Regret Although perfect information problem formulations constitute a useful tool for exploring optimal system development under different scenarios, they disregard some aspects of real-world decisionmaking that become critical when the decisions to be taken are capital-intensive and irreversible. The first aspect relates to the lack of perfect information and the need to endogenously consider the possibility for alternative future developments. In view of this, a flexible decision framework enables a planner to move beyond the concept of a static investment plan and instead identify the optimal investment strategy which encapsulates a range of contingent courses of actions to be taken according to all possible paths of uncertainty evolution. Furthermore, a flexible decision framework must also be accompanied by a specific decision criterion that describes the planner’s attitude towards risk and ultimately determines which commitments optimally balance cost efficiency and risk exposure. In general, there are three main classes of decision criteria when facing uncertainty; stochastic (also known as probabilistic), risk-constrained and robust. Stochastic planning is the case where each scenario node is attributed a probability of occurrence; the planner’s objective is the minimization of expected system cost over all realisations as in (Konstantelos and Strbac, 2015) and (Munoz et al., 2014). In a similar vein, suitable constraints with respect to spectral risk measures such as the Conditional Value-at-Risk, as in (Rockafellar and Uryasev, 2000), can be included to render the planner risk-averse. Robust decision methods in the context of system planning can refer to two variants; uncertainty intervals and utilisation of the regret concept. The former guarantees optimal performance given a deterministic description of the uncertainty state space (Bertsimas et al., 2011). Naturally, such a formulation lends itself mostly to static descriptions of uncertainty and cannot take advantage of its inter-temporal resolution structure which is an important characteristic of dynamic system planning (Epstein et al., 1980). The latter identifies the optimal planning strategy so as to 4

minimize a planner’s worst-case regret under all considered realizations (Gorenstin et al., 1993); the regret related to a particular scenario is the difference between the realized system cost and the system cost that would have been incurred under the assumption of perfect foresight. The concept of regret has been successfully used in the past in the context of transmission planning problems (van der Weijde et al., 2012; Chen et al., 2014). For example, authors in (van der Weijde et al., 2012) develop a minimax regret transmission planning model to examine the cost of disregarding uncertainty associated with future economic, technology and regulatory developments in the UK. The attractiveness of adopting a β€˜minimax regret’ decision criterion lies in its ex-post decision-making nature; decisions are based on consequences of scenario occurrences. Probabilistic models, on the other hand, identify optimal decisions without considering the consequences after uncertainty has been resolved (Miranda and Proenca, 1998). The former philosophy is very much in line with existing regulatory frameworks which are based on ex-post performance evaluation against ex-ante-defined targets. Note that by definition, regret is defined in terms of the optimal scenario-specific solutions obtained under perfect foresight πœ”π‘  . It follows that the perfect information problem for each scenario must be solved before proceeding with the MILP min-max regret formulation shown below. 𝐦𝐒𝐧{π‘Ÿ} π‘Ÿ β‰₯ π‘Ÿπ‘  π‘Ÿπ‘  = βˆ‘

𝐼 (π‘Ÿπ‘’πΌ π‘π‘š

+

𝑂) π‘Ÿπ‘’π‘‚ π‘π‘š βˆ’

πœ”π‘ 

,βˆ€π‘  ∈ Ω𝑆

(16) (17)

,βˆ€π‘  ∈ Ω𝑆

(18)

,{βˆ€π‘š: 𝑒(π‘š) = 1, π‘š β‰  1}, 𝑙 ∈ Ω𝐿

(19)

βˆ€π‘šβˆˆπ‘ 

𝐹1,𝑙 = πΉπ‘š,𝑙

Subject to constraints (2)-(15) applied to all scenario trees nodes π‘š ∈ Ω𝑀 . The objective function (16) minimizes the maximum regret π‘Ÿ. Constraint (17) ensures that π‘Ÿ is equal or greater to the regret of each scenario π‘Ÿπ‘  . The combination of constraints (16) and (17) applies the min-max operator across all scenarios; the worst case realization does not have to be defined a priori but is a result of the optimization. Equation (18) computes the regret for each scenario by aggregating the investment and operation costs across all scenario tree nodes that belong to scenario 𝑠 and subtracting the optimal cost πœ”π‘  obtained from the deterministic problem associated to scenario 𝑠. Equation (19) imposes non-anticipativity; decisions taken in the first epoch are common across all scenarios. Note that although for simplicity solely first-stage branching is considered, the presented optimization framework employs a node-variable formulation that could be extended in a straightforward manner to scenario trees of arbitrary size and shape.

Complementary references Bertsimas, D., Brown, D. B., & Caramanis, C. (2011). Theory and applications of robust optimization. SIAM review, 53(3), 464-501. Binato, S., Pereira, M. V. F., & Granville, S. (2001). A new Benders decomposition approach to solve power transmission network design problems. IEEE Transactions on Power Systems, 16(2), 235-240. Chen, B., Wang, J., Wang, L., He, Y., & Wang, Z. (2014). Robust optimization for transmission expansion planning: Minimax cost vs. minimax regret. IEEE Transactions on Power Systems, 29(6), 3069-3077. Epstein, L. G. (1980). Decision making and the temporal resolution of uncertainty. International economic review, 269-283. 5

FICO Xpress (2009). Xpress-optimizer http://www.fico.com/en/node/8140?file=5097.

reference

manual.

[Online],

Available:

Gorenstin, B. G., Campodonico, N. M., Costa, J. P., & Pereira, M. V. F. (1993). Power system expansion planning under uncertainty. IEEE Transactions on Power Systems, 8(1), 129-136. National Grid (2017). [Online], Available: http://investors.nationalgrid.com/about-us/our-markets/ukprofile.aspx. Rockafellar, R. T., & Uryasev, S. (2000). Optimization of conditional value-at-risk. Journal of risk, 2, 21-42.

6