Next Article in Journal
Investigation of the Internal Displacement Chamber Pressure of a Rotary Vane Pump
Next Article in Special Issue
PCA Forecast Averaging—Predicting Day-Ahead and Intraday Electricity Prices
Previous Article in Journal
Electromagnetic and Calorimetric Validation of a Direct Oil Cooled Tooth Coil Winding PM Machine for Traction Application
Previous Article in Special Issue
Impact of Long-Term Water Inflow Uncertainty on Wholesale Electricity Prices in Markets with High Shares of Renewable Energies and Storages
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solving the Stochastic Generation and Transmission Capacity Planning Problem Applied to Large-Scale Power Systems Using Generalized Shift-Factors

Department of Electrical Engineering, Universidad Técnica Federico Santa María, Valparaíso 2390123, Chile
*
Author to whom correspondence should be addressed.
Submission received: 5 June 2020 / Revised: 22 June 2020 / Accepted: 27 June 2020 / Published: 30 June 2020
(This article belongs to the Special Issue Uncertainties and Risk Management in Competitive Energy Markets)

Abstract

:
In this study, we successfully develop the transmission planning problem of large-scale power systems based on generalized shift-factors. These distribution factors produce a reduced solution space which does not need the voltage bus angles to model new transmission investments. The introduced formulation copes with the stochastic generation and transmission capacity expansion planning problem modeling the operational problem using a 24-hourly load behaviour. Results show that this formulation achieves an important reduction of decision variables and constraints in comparison with the classical disjunctive transmission planning methodology known as the Big M formulation without sacrificing optimality. We test both the introduced and the Big M formulations to find out convergence and time performance using a commercial solver. Finally, several test power systems and extensive computational experiments are conducted to assess the capacity planning methodology. Solving deterministic and stochastic problems, we demonstrate a prominent reduction in the solver simulation time especially with large-scale power systems.

1. Introduction

The capacity expansion planning (CEP) problem solves the following issues—(1) how much generation and transmission capacity requires the electrical power system, (2) when and what type of new power units and new transmission elements are needed and (3) where locate these new elements (generation units and transmission transformer and lines) to supply adequately the long-term energy of the customers.

1.1. Technical Literature Review

Assuming that the electricity market is centrally operated, for instance the Chilean situation, the CEP problem is formulated as a cost minimization problem considering both investment and operational components. Notice that the operational problem has been traditionally formulated using a linear optimal power flow.
Even though the expansion planning problem was carried out in 1957 [1], this optimization problem is presently figured out by innovative computational and/or decomposition techniques [2]. In the technical literature, different CEP methodologies have been applied to solve short- and medium-scale power systems adopting different objective functions, sets of decision variables and different constraints among other technical considerations [3].
The CEP formulation [4,5] solves the generation capacity expansion planning (GCEP) problem and the transmission capacity expansion planning (TCEP) problem simultaneously. In the CEP problem, the number of decision variables and constraints increases exponentially solving medium-scale power systems as well as modeling different energy customers estimations (load duration curves or an hourly data).
Traditionally, the CEP problem is divided using two optimization frameworks: (1) the GCEP problem and (2) the TCEP problem. This is accomplished not only because the investment infrastructure involved in each type of decision is different but also because the way of the new transmission decisions modify the transmission network approach. On the one side, the GCEP problem produces local modifications in the operational problem. On the other side, the TCEP problem makes important changes in the bus admittance matrix.
To solve the TCEP problem, the most popular mixed-integer linear problem (MILP) is the classical disjunctive Big M model [6]. Using the voltage bus angles as decision variables, the Big M formulation incorporates two groups of binary inequality constraints for modeling new transmission investments. The main disadvantage is the higher number of decision variables and constraints which depends on the electrical power system, the investment transmission candidates and the generation/load time-steps. Although the disjunctive model displays an improvement in linearity [7], this framework enables coping with stochastic problems solving medium- and large-scale power systems [8].
To solve the GCEP problem, Hinojosa et al. [9,10,11] accomplish several studies using traditional shift-factors without include the voltage bus angles to model the bulk transmission network [12,13]. Nevertheless, these linear factors cannot be applied directly to solve the TCEP problem because nonlinear terms appear in candidate transmission circuits.
To keep the power systems topology and generalized shift-factors (GSF) unchanged, the line outage is modeled with a virtual generation located in the sending bus m and a virtual load located in the receiving bus n. Using this GSF-based methodology, the power injection to bus m and withdrawal from bus n is canceled out each other with no impact on other transmission elements. Review reference [14] for a detailed information about the GSF applied to the security N-1 power flow problem.
In [15,16], authors formulate the TCEP problem taking into account GSF. These factors are calculated considering that existing and future transmission elements are included in the transmission network. The transmission investment decisions are added in the planning problem using the GSF-based methodology based on virtual power injections and a bilinear formulation. Applying a disjunctive technique, it is possible to transform one nonlinear power flow equation to two linear binary inequality constraints. Notice that this methodology eliminates the voltage bus angles as decision variables. Therefore, there is a very important reduction in decision variables and constraints.
In [15], the load of the customers is modeled in the multi-stage problem by only one yearly value. Besides, authors solve up to IEEE 118-bus power system using a deterministic approach. In our opinion, there is necessary to conduct several simulations in order to find out the performance of the GSF-based methodology. We propose to solve the multi-stage CEP problem applied to large scale problem as well as to model an hourly behaviour. Accordingly, we implement an hourly load modeling to obtain a resilient electrical power system that will be flexible and robust to withstand, for example, renewable power variability.
In the state-of-the-art, a few studies adequately model the long-term uncertainty in the planning problem because of the high computational complexity [17]. Although the idea of implementing stochastic-mixed integer programming is not new (over 15 years) [18], recently MILP improvements as well as high computational capacity have done possible to solve medium-scale power systems. For more information about uncertainty implementations, review references [2,10,11,16,17].
In this study, we carry out the stochastic CEP problem based on a two-stage MILP problem [19]. The mathematical formulation uses scenarios to represent uncertainty in the input variables. To obtain a flexible electrical power system, the load uncertainty is implemented by an hourly formulation using representative days for expansion decisions [20].
To accomplish effectively the stochastic CEP problem, it is necessary to include simultaneously both GCEP and TCEP formulations in order to evaluate performance and simulation time. We want to emphasize that there are no studies solving both planning problems [15,16]. Additionally, we conduct several analyses solving large-scale power systems.

1.2. Contributions

Based on the literature review, we could mention that development of new planning models is still an open issue in the power system field. The proposed CEP methodology has not been developed in the technical literature. Therefore, it would be very interesting to determine the performance of the methodology solving large-scale power systems. The following contributions are expected to accomplish in this study: (1) the stochastic CEP problem is solved using the GSF and (2) a very large-scale power system is successfully solved in a reasonable simulation time. In our opinion, the achieved stochastic formulation could bring better performance and practical advantages solving very large-scale problems.
This study is organized as follows. In Section 2, the CEP formulation is shown and applied to an example 4-bus power system. Section 3 presents results and comparisons using several power systems. Last, Section 4 concludes the study.

2. Introducing the Stochastic Capacity Expansion Planning Formulation

From a central planner point of view, the stochastic CEP problem is formulated as the minimization (1a) of the following costs: (1) generation C I G (1b) and transmission C I T (1c) investment costs; and (2) operational and unserved energy costs C O P (1d). Notice that the operational problem is modeled using a DC-network through a linear optimal power flow and the transmission DC power losses are not included in the optimization problem.
The mathematical formulation uses scenarios to represent uncertainty. Therefore, the stochastic problem is carried out from a deterministic MILP formulation as the minimization of the expected value of the deterministic operational costs (1a) for E possible realizations of the uncertainty parameter with their respective probabilities of occurrence ( P r e ) [19].
In the two-stage CEP problem, the first stage variables, which are generation and transmission investment candidates, are decisions variables with uncertainty. These variables are the same for all different scenarios using non-anticipativity constraints [11]. The second stage variables, which are the power generation of each unit and virtual candidate power flows [16]), must supply each load scenario. Mathematically, the stochastic CEP problem is formulated as follows:
min t T R t · [ C I G + C I T + e E P r e · s S H s · C O P ]
where:
C I G = g G B ( C g · ( n g , t , e n g , t 1 , e ) + O M g · n g , t , e ) · P g M
C I T = k m L B C k m y x k m M ( x k m , t , e , y x k m , t 1 , e , y )
C O P = g ( G E G B ) F C g · p g , s , t , e + V o L L · g G R p g , s , t , e
The optimization problem is subject to the following technical constrains:
g G E P g M + g G B P g M · n g , t , e D t , e M · S R
g ( G E G B G R ) h s · p g , s , t , e = D s , t , e t o t a l
| b B G S F k m , b · ( p b , s , t , e D b , s , t , e ) + k m L B y X k m M ( G S F k m , k G S F k m , m ) · f ˜ k m , s , t , e , y | F k m M k m L E
| f ˜ k m , s , t , e , y b B G S F k m , b · ( p b , s , t , e D b , s , t , e ) k m L B y X k m M ( G S F k m , k G S F k m , m ) · f ˜ k m , s , t , e , y | x k m , t , e , y · F k m M k m L B
| f ˜ k m , s , t , e , y | ( 1 x k m , t , e , y ) · M k m L B
x k m , t , e , y x k m , t 1 , e , y t T , t > 1
n g , t , e n g , t 1 , e t T , t > 1
n g , t , e n g , t 1 , e N g M t T , t > 1
y X k m M ( x k m , t , e , y x k m , t 1 , e , y ) X k m M t T , t > 1
0 p g , s , t , e P g M g ( G E G R )
0 p g , s , t , e P g M · n g , t , e g G B
n g , t , 1 = = n g , t , e g G B
x k m , t , 1 , y = = x k m , t , e , y k m L B
0 x k m , t , e , y 1
In (2), the planning reserve margin measures the amount of generation capacity available for unexpected generation/load events. Equation (3) represents the total load balance. Equations (4)–(6) accomplish a binary disjunctive formulation using GSF and virtual power flows [15,16]. Constraints (7) and (8) ensure coherence for both transmission and generation investment variables, respectively. Equations (9) and (10) limits the generation and transmission investment variables. Equation (11) shows the maximum power generation for existing power units. Equation (12) includes a disjunctive constraint to model investment generation capacity. Last, the non-anticipativity constraints for both generation and transmission variables are shown in (13) and (14), respectively. A binary variable models the candidate transmission elements (15).
For new transmission decisions, variable x k m represents k m line status (ON/OFF) as follows: x k m = 1 indicates that a new line is built and x k m = 0 implies that there is no investment decision.

The Proposed Formulation Applied to an Example System (4-Bus)

In Figure 1, we present a 4-bus power system used to apply the deterministic CEP optimization problem. Transmission reactance and thermal limits (MW) are also included.
Table 1 gives information for existing and candidate power units. We also include investment and operational costs.
Table 2 shows transmission investment data. Moreover, we can see investment cost for transmission future options.
These artificial data are employed to realize effects and insights solving the introduced methodology.
Table 3 shows the GSF matrix using bus 1 as the slack bus. Notice that shift-factor values consider both existing and future transmission elements.
For illustrative purposes, generation and transmission inequality constraints (4)–(6) and (12) are developed in this section. We use M = 2 π * S b a s e ; where S b a s e = 100.
Using (4) for the existing 1 2 , 1 3 and 2 3 transmission elements, the following constraints must be formulated (16), (17) and (18), respectively.
203.45 0.420 p B 0.420 p C 0.360 p D + 0.420 f ˜ 12 , 1 0.060 f ˜ 24 , 1 0.060 f ˜ 24 , 2 + 0.040 f ˜ 34 , 1 + 0.040 f ˜ 34 , 2 + 0.040 f ˜ 34 , 3 3.45
193.10 0.159 p B 0.159 p C 0.283 p D + 0.159 f ˜ 12 , 1 + 0.124 f ˜ 24 , 1 + 0.124 f ˜ 24 , 2 0.083 f ˜ 34 , 1 0.083 f ˜ 34 , 2 0.083 f ˜ 34 , 3 6.90
127.59 0.034 p B + 0.034 p C 0.069 p D 0.034 f ˜ 12 , 1 + 0.103 f ˜ 24 , 1 + 0.103 f ˜ 24 , 2 0.069 f ˜ 34 , 1 0.069 f ˜ 34 , 2 0.069 f ˜ 34 , 3 72.41
Using (5) and (6), we develop disjunctive constraints for the new transmission candidates as follows: (a) for line 1–2 (19) and (20); (b) for the first circuit of line 2–4 (21) and (22); (c) for the second circuit of line 2–4 (23) and (24); (d) for the first circuit of line 3–4 (25) and (26); (e) for the second circuit of line 3–4 (27) and (28); and (f) for the third circuit of line 3–4 (29) and (30).
| 0.420 p B + 0.420 p C + 0.360 p D + 0.580 f ˜ 12 , 1 + 0.060 f ˜ 24 , 1 + 0.060 f ˜ 24 , 2 0.040 f ˜ 34 , 1 0.040 f ˜ 34 , 2 0.040 f ˜ 34 , 3 + 100 x 12 , 1 | 103.45
| f ˜ 12 , 1 628.32 x 12 , 1 | 628.32
| 0.062 p B 0.062 p C + 0.324 p D + 0.062 f ˜ 12 , 1 + 0.614 f ˜ 24 , 1 0.386 f ˜ 24 , 2 0.076 f ˜ 34 , 1 0.076 f ˜ 34 , 2 0.076 f ˜ 34 , 3 + 150 x 24 , 1 | 89.66
| f ˜ 24 , 1 628.32 x 24 , 1 | 628.32
| 0.062 p B 0.062 p C + 0.324 p D + 0.062 f ˜ 12 , 1 0.386 f ˜ 24 , 1 + 0.614 f ˜ 24 , 2 0.076 f ˜ 34 , 1 0.076 f ˜ 34 , 2 0.076 f ˜ 34 , 3 + 150 x 24 , 2 | 89.66
| f ˜ 24 , 2 628.32 x 24 , 2 | 628.32
| 0.040 p B + 0.040 p C + 0.120 p D 0.040 f ˜ 12 , 1 0.080 f ˜ 24 , 1 0.080 f ˜ 24 , 2 + 0.720 f ˜ 34 , 1 0.280 f ˜ 34 , 2 0.280 f ˜ 34 , 3 + 150 x 34 , 1 | 6.90
| f ˜ 34 , 1 628.32 x 34 , 1 | 628.32
| 0.041 p B + 0.041 p C + 0.117 p D 0.041 f ˜ 12 , 1 0.076 f ˜ 24 , 1 0.076 f ˜ 24 , 2 0.283 f ˜ 34 , 1 + 0.717 f ˜ 34 , 2 0.283 f ˜ 34 , 3 + 150 x 34 , 2 | 6.90
| f ˜ 34 , 2 628.32 x 34 , 2 | 628.32
| 0.041 p B + 0.041 p C + 0.117 p D 0.041 f ˜ 12 , 1 0.076 f ˜ 24 , 1 0.076 f ˜ 24 , 2 0.283 f ˜ 34 , 1 0.283 f ˜ 34 , 2 + 0.717 f ˜ 34 , 3 + 150 x 34 , 3 | 6.90
| f ˜ 34 , 3 628.32 x 34 , 3 | 628.32
Using (12) for candidate power generation units, constraint (31) model new power unit G C and constraint (32) model new power unit G D .
p C 250 n 2 0
p D 100 n 4 0
These mathematical constraints are presented to clearly explain the accomplished CEP problem.
The static planning problem is solved using the maximum load. Reserve requirements are not considered and the MILP gap is 0.0001%.
Using G u r o b i (9.0.0) [21] to solve the CEP problem, the optimal cost is 33.7288MM$ (1a). The investment planning problem includes one power unit at bus 4 and one transmission line 2 4 .
For the transmission planning problem, new line 2 4 (investment decision x 24 , 1 = 1) means the virtual power flow is f ˜ 24 , 1 = 0 MW. The optimal solution does not require other transmission investments to supply the load of the customers. That is, x 12 , 1 = 0, x 24 , 2 = 0, x 34 , 1 = 0, x 34 , 2 = 0 and x 34 , 3 = 0. Additionally, G u r o b y displays the following virtual power flows: f ˜ 12 , 1 = 100 MW, f ˜ 24 , 2 = 90 MW, f ˜ 34 , 1 = 40 MW, f ˜ 34 , 2 = 40 MW and f ˜ 34 , 3 = 40 MW. Notice that these virtual power flows satisfy candidate transmission constraint (5) when decision variable is zero.
Figure 2 displays the power flow solution for existing and investment lines. We do not include all virtual investment power flows because these values are injected (generation) in the sending bus and withdrawal (load) in the receiving bus with no impact on the other transmission power flows.
It should be mentioned that the classical Big M formulation also obtains the same solution using a higher number of decision variables and constraints.

3. Simulation Results

In Figure 3, we present a diagram showing power system requirements as well as planning variables and parameters. Using this information, both generation and transmission decisions and hourly operational solution can be reported using the introduced formulation.
Two mathematical formulations are applied to solve the CEP problem—(1) the introduced formulation (GSF) and (2) the classical disjunctive Big M formulation (CD). The CD formulation validates optimality and performance solving the CEP problem. Several simulations are conducted using short- and large-scale electrical power systems to find out convergence and simulation time.
Matlab [22] (Version 2015b, UTFSM, Chile) develops the mathematical capacity planning formulation and G u r o b i solves the optimization problem on a computer with the following characteristics: Intel Xeon E5-2630 (2.20 GHz) with a RAM of 32 GB.

3.1. Garver Test Power System

This classical power system [23] is traditionally used to present new planning formulations as well as to establish performance and optimality. Review reference [24] to obtain more information.
Table 4 shows generation data including buses where the power units are located, the maximum power generation, the investment and operation & maintenance costs, and the fuel cost for existing and future power units.
The CEP problem considers a reserve margin of 20%; that is 1.2 times the yearly peak load (760 MW). The V o L L value is 10,000 [$/MWh]. For the multi-stage problem, the interest rate is 10%. For the MILP solver, the optimality gap is 0.0001%.

3.1.1. Static Deterministic Analysis

Applying the GSF-based formulation, the optimal cost is 475, 809, 470.91$; the transmission investment cost is 110 MM$ and the generation investment cost is 247.2 MM$. The yearly operational cost is 118.6095 MM$. The optimal solution does not display unserved energy.
The generation investment plan considers two units G 4 (bus 3) and two units G 6 (bus 6). Though the investment cost of power unit G 5 is the lowest, that power unit is not needed to supply the load of the customers because of the highest fuel cost.
The power generation solution is p 1 = 90.0 MW, p 2 = 60.0 MW, p 3 = 70.61 MW, p 4 = 240.0 MW and p 6 = 299.40 MW. In the original TCEP problem [23], the maximum capacity at bus 6 is 600 MW. However, a lower power generation produced by G 6 is enough to supply load demand.
The transmission investment plan incorporates three circuits in corridor x 2 6 and one circuit in corridor x 3 5 . This solution is well known in the technical literature [23,24].
Table 5 displays simulation results for both GSF and CD formulations. The authors include the number of decision variables (v), the number of equality and inequality constraints (Θ), the non-zeros elements ( n Z ), the objective function ( O F ) and the average solver time obtained by G u r o b i ( t s ) simulating 200–trials.
Results show that both formulations achieve the same optimal solution. Consequently, the introduced and the Big M formulations are equivalent using different solution spaces.
Figure 4 displays the power flow solution. We do not include virtual power flows when there is no investment decision.

3.1.2. Multi-Stage Deterministic and Stochastic Analyses Using an Hourly Load Modeling

Tenth investment years are incorporated to carry out the multi-stage CEP problem. Review reference [16] to have more information about the long-term load forecasting.
Figure 5 shows the hourly load behaviour. To formulate the multi-stage CEP problem, the yearly load demand is proportionally divided using these hourly data.
For the stochastic analysis, ten load scenarios E = 10 are created using the following factors: 1.0490, 0.9867, 0.9726, 0.9952, 0.9310, 0.9636, 1.0943, 0.8529, 1.0490, and 0.9404. To create five scenarios E = 5 , we use the first five values. We assume each scenario has the same probability of occurrence with P r = 0.2 for E = 5 and P r = 0.1 for E = 10 .
Table 6 shows simulation results performed by GSF-based formulation. We show the optimal objective function and the solver simulation time t s achieved by the deterministic ( E = 1 ) problem and the two stochastic ( E = 5 and E = 10 ) problems. Results obtained by the Big M formulation (CD) are also included to verify optimality.
For the deterministic problem, the generation planning solution includes two G 4 units one in the second year ( t 2 ) and one in the seventh year ( t 7 ); and two G 6 units one in the fourth year ( t 4 ) and one in the ninth year ( t 9 ).
The investment transmission solution considers one corridor of x 3 5 in the eighth year ( t 8 ); and three corridors of x 2 6 , two in the fourth year ( t 4 ) and one in the eighth year ( t 8 ). Regarding the operational problem, there is no unserved energy.
Although the same optimal investment plan is also obtained by the stochastic formulation, there are differences in the expected operational costs.
In the GSF-based formulation, decisions variables increase. However, there is a little reduction in the number of constraints. Also, the non-zero elements increase. Notice that simulation time grows exponentially when different load scenarios are added in the optimization problem.
In these simulations, both formulations obtain the same solution. Review reference [25] to obtain more simulation results.
Even though the lower simulation time is obtained by the traditional Big M formulation, we conduct more analyses not only to clarify advantages or disadvantages of the proposed formulation but also to determine the performance using large-scale power systems.

3.2. IEEE 300-Bus Test Power System

This study case (http://www.ee.washington.edu/research/pstca/) includes several voltages (345, 230, 115, 86, 66, 27, 20, 16.5, 13.8, 6.6, 2.3, and 0.6 kV). The power system contains sixty-nine generators and four hundred eleven transmission lines. Data for existing power units can be reviewed using MatPower [26].
To model the stochastic CEP problem, tenth investment long-term periods and five load scenarios are modeled using hourly data. For each load scenario, we select the hourly behaviour of bus 1 (Figure 5).
To determine generation and transmission candidates, we solve an optimal power flow using the maximum load of the customers. Future planning options are selected based on unserved energy and transmission congestion problems as follows:
  • For the power unit candidates, we include fifteen new power generation options with a maximum number of N g M = 1.
  • For the transmission network candidates, there are thirty investments options where each line has a maximum power flow of 390 MW and the maximum number of new elements is X l M = 2 .
The CEP problem considers a reserve margin of 20%. The V o L L value is 10,000 $/MWh, the interest rate is 10% and the MILP gap is 1%.
Table 7 shows results and time performance for both deterministic (one scenario) and stochastic (five scenarios) CEP formulations.
The optimal investment planning solution consist of thirty-one transmission elements and four new power generation units.
In the deterministic analysis (E = 1), the optimal cost is the same for both GSF and CD formulations. More importantly, we achieve the same CEP solution using a reduced solution space. With respect to the classical Big M solution space, the introduced formulation has 79.7% and 50.0% lower number of decision variables and constraints, respectively. The outstanding result is the lower simulation time obtained by the GSF-based formulation, which is 95.4%—21.8 times—faster in comparison with the CD formulation.
Considering 5 load scenarios and using the traditional Big M stochastic formulation, we are not able to obtain a reasonable solution after a simulation time of 3.13 weeks. At that time, the G u r o b i gap is approximately 6%. Therefore, we decide to stop the computer simulation because the solution achieved by the GSF-based formulation had a MILP gap lower than 1% in 272,240 s (3.2 days).
The hourly optimal solutions show there is no load shedding in the electrical power system.

3.3. The PEGASE 1354-Bus Test Power System

This case accurately represents the size and complexity of part of the European high voltage transmission network. The network contains 1354 buses, 260 generators, and 1991 branches, and these lines operates at 380 and 220 kV. Data is created from the Pan European Grid Advanced Simulation and State Estimation (PEGASE) project [27]. Power system information can be reviewed using M a t P o w e r .
Considering not only infinite simulation time but also memory M a t l a b problems, we are able to model only the deterministic planning problem taking into account two investment periods. We include the end-of-time operational component using five years as well.
To determine generation and transmission candidates, we solve an optimal power flow with the maximum load.
For the planning generation plan, we propose twenty-eight power generation options with N g M = 1 .
For the transmission network plan, there are one-hundred and seven candidate lines. These transmission elements have a thermal limit of 500 MW and the maximum number of transmission candidates is X l M = 2 .
Because M a t p o w e r does not include fuel variable costs, we use random numbers provided by a uniform discrete distribution function using as limits 40 and 100 $/MWh.
In the mathematical formulation, the reserve margin is 20%, the V o L L value is 10,000 $/MWh, the interest rate is 10% and the optimality gap is 3%. The hourly load behaviour is based on bus 1.
In this very large-scale problem, the G u r o b y solver time is approximately three weeks using the GSF-based formulation. We want to highlight that the traditional Big M formulation does not carry out the CEP optimization problem because of memory issues.
Using M a t l a b to formulate the planning problem, the accomplished CEP methodology is able to solve up to 1354-bus power systems using a 24-hourly load behaviour. Considering higher power system dimension and memory issues, stochastic generation and transmission planning problems should be modeled and simulated by decomposition techniques such as the Benders approach [28].

4. Conclusions

This study introduced an efficient planning methodology to formulate large-scale transmission networks based on generalized shift-factors. Several analyses were conducted to determine effects and insights solving stochastic planning problems with different uncertainty scenarios. Although the traditional disjunctive Big M formulation displays an improvement modeling the transmission planning problem, this methodology only had a better performance solving small-scale power systems. With reference to the stochastic CEP problem, the Big M formulation did not obtain a high quality solution applied to the IEEE 300-bus test power system in comparison with the solution achieved by the developed formulation. Furthermore, the Big M formulation was not able to model the deterministic CEP optimization problem applied to the PEGASE 1354-bus power system because of M a t l a b memory issues.
The outstanding result of this study is the best simulation time obtained by the generalized shift-factors formulation because the optimization problem was modeled using a lower number of decision variables and constraints without sacrificing optimality. Results lead to the conclusion that the proposed stochastic formulation could be applied to formulate real large-scale power systems. Even though this methodology reduced the computational burden of stochastic planning methodologies, we were limited by M a t l a b memory capacity to solve up to 1354-bus power systems. However, we are evaluating to implement the achieved formulation using P y t h o n .
Considering the best simulation time and practical advantages obtained by the proposed formulation, the research group will include operational generation Unit Commitment constraints to the stochastic CEP problem.

Author Contributions

In this study, all the authors were involved in the mathematical formulation, simulation, results analysis and conclusions as well as manuscript preparation. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded and supported by the Universidad Técnica Federico Santa María, Chile.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mass, P.; Gibrat, R. Application of linear programming to investment in the electric industry. Manag. Sci. 1957, 3, 149–166. [Google Scholar] [CrossRef]
  2. Akhavizadegan, F.; Wang, L.; McCalley, J. Scenario selection for iterative stochastic transmission expansion planning. Energies 2020, 13, 1203. [Google Scholar] [CrossRef] [Green Version]
  3. Quintero, J.; Zhang, H.; Chakhchoukh, Y.; Vittal, V.; Heydt, G. Next generation transmission expansion planning framework: Models, tools, and educational opportunities. IEEE Trans. Power Syst. 2014, 29, 1911–1918. [Google Scholar] [CrossRef]
  4. Khodaei, A.; Shahidehpour, M. Microgrid-based co-optimization of generation and transmission planning in power systems. IEEE Trans. Power Syst. 2013, 28, 1582–1590. [Google Scholar] [CrossRef]
  5. Aravena, I.; Cárdenas, R.; Gil, E.; Hinojosa, V.; Araneda, J.C.; Reyes, P. Co-optimization of generation and transmission investment decisions under hydro uncertainty using stochastic mixed-integer programming. In Proceedings of the 10th Latin-American Congress on Electric Power Generation, Transmission and Distribution (CLAGTEE 2013), Vina del Mar, Chile, 6–9 October 2013; pp. 1–8. [Google Scholar]
  6. Villanasa, R. Transmission Network Planning Using Linear and Mixed Linear Integer Programming. Ph.D. Thesis, Ressenlaer Polythechnic Institute, Troy, NY, USA, 1984. [Google Scholar]
  7. Moreira, A.; Pozo, D.; Street, A.; Sauma, E. Reliable renewable generation and transmission expansion planning: Co-optimizing system’s resources for meeting renewable targets. IEEE Trans. Power Syst. 2017, 32, 3246–3257. [Google Scholar] [CrossRef]
  8. Lumbreras, S.; Ramos, A.; Banez-Chicharro, F. Optimal transmission network expansion planning in real-sized power systems with high renewable penetration. Electr. Power Syst. Res. 2017, 149, 76–88. [Google Scholar] [CrossRef]
  9. Hinojosa, V.H.; Velásquez, J. Improving the mathematical formulation of security-constrained generation capacity expansion planning using power transmission distribution factors and line outage distribution factors. Electr. Power Syst. Res. 2016, 140, 391–400. [Google Scholar] [CrossRef]
  10. Hinojosa, V.H.; Velásquez, J. Stochastic security-constrained generation expansion planning based on linear distribution factors. Electr. Power Syst. Res. 2016, 140, 139–146. [Google Scholar] [CrossRef]
  11. Hinojosa, V.H. A generalized stochastic N-m security-constrained generation expansion planning methodology using partial transmission distribution factors. In Proceedings of the 2017 IEEE Power & Energy Society General Meeting, Chicago, IL, USA, 16–20 July 2017. [Google Scholar]
  12. Hinojosa, V.H. Comparing Corrective and Preventive Security-Constrained DCOPF Problems Using Linear Shift-Factors. Energies 2020, 13, 516. [Google Scholar] [CrossRef] [Green Version]
  13. Hinojosa, V.H.; Gonzales-Longatt, F. Preventive security-constrained DC-OPF formulation using power transmission distribution factors and line outage distribution factors. Energies 2018, 11, 1497. [Google Scholar] [CrossRef] [Green Version]
  14. Wood, A.J.; Wollenberg, B.F. Power Generation, Operation and Control; John Wiley & Sons: Hoboken, NJ, USA, 1996. [Google Scholar]
  15. Rahmani, M.; Kargarian, A.; Hug, G. Comprehensive power transfer distribution factor for large-scale transmission expansion planning. IET Gen. Trans. Dist. 2016, 10, 2981–2989. [Google Scholar] [CrossRef]
  16. Hinojosa, V.H.; Sepúlveda, J. Solving the stochastic transmission expansion capacity planning methodology based on a generalized partial transmission distribution factors. In Proceedings of the IEEE PowerTech Conference, Milan, Italy, 23–27 June 2019. [Google Scholar]
  17. Liang, Z.; Chen, H.; Wang, X.; Idris, I.I.; Tan, B.; Zhang, C. An extreme scenario method for robust transmission expansion planning with wind power uncertainty. Energies 2018, 11, 2116. [Google Scholar] [CrossRef] [Green Version]
  18. Ryan, S.; McCalley, J.; Woodruff, D. Long Term Resource Planning for Electric Power Systems under Uncertainty. 2011. Available online: http://works.bepress.com/sarah_m_ryan/27/ (accessed on 30 June 2020).
  19. Birge, J.; Louveuax, F. Introduction to Stochastic Programming; ser. Springer Series in Operations Reserach and Financial Engineering; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  20. García-Cerezo, A.; Baringo, L.; García-Bertrand, R. Representative days for expansion decisions in power systems. Energies 2020, 13, 335. [Google Scholar] [CrossRef] [Green Version]
  21. Gurobi Optimization. Available online: https://www.gurobi.com/ (accessed on 1 June 2020).
  22. Matlab. Available online: https://www.mathworks.com/ (accessed on 1 June 2020).
  23. Garver, L.L. Transmission network estimation using linear programming. IEEE Trans. Power Apparat. Syst. 1970, PAS-89, 1688–1697. [Google Scholar] [CrossRef]
  24. Hinojosa, V.H.; Galleguillos, N.; Nuques, B. A simulated rebounding algorithm applied to the multi-stage security-constrained transmission expansion planning in power systems. Int. J. Elect. Power Energy Syst. 2013, 47, 168–180. [Google Scholar] [CrossRef]
  25. Sepúlveda, J. Planificación de la Capacidad de Generación y Transmission Usando Factores de Sensibilidad Considerando Incertidumbre en las Variables de Entrada, Contingencias y Restricciones de Flexibilidad. Bachelor’s Thesis, Departamento of Electrical Engineering, Universidad Técnica Federico Santa María, Valparaíso, Chile, 2019. (In Spanish). [Google Scholar]
  26. Zimmerman, R.D.; Murillo-Sanchez, C.E.; Thomas, R.J. MATPOWER: Steady-state operations, planning, and analysus tools for power systems research and education. IEEE Trans. Power Syst. 2011, 26, 12–19. [Google Scholar] [CrossRef] [Green Version]
  27. Fliscounakis, S.; Panciatici, P.; Capitanescu, F.; Wehenkel, L. Contingency ranking with respect to overloads in very large power systems taking into account uncertainty, preventive and corrective actions. IEEE Trans. Power Syst. 2013, 28, 4909–4917. [Google Scholar] [CrossRef] [Green Version]
  28. Majidi-Qadikolai, M.; Baldick, R. A generalized decomposition framework for large-scale transmission expansion planning. IEEE Trans. Power Syst. 2018, 33, 1635–1649. [Google Scholar] [CrossRef]
Figure 1. Power system example.
Figure 1. Power system example.
Energies 13 03327 g001
Figure 2. Power flow solution.
Figure 2. Power flow solution.
Energies 13 03327 g002
Figure 3. Capacity expansion planning framework.
Figure 3. Capacity expansion planning framework.
Energies 13 03327 g003
Figure 4. Power flow solution for the Garver power system.
Figure 4. Power flow solution for the Garver power system.
Energies 13 03327 g004
Figure 5. Hourly load demand for a typical day.
Figure 5. Hourly load demand for a typical day.
Energies 13 03327 g005
Table 1. New and future power generation data.
Table 1. New and future power generation data.
# P g M C g FC g N M
(MW)($)($/MWh)
G A 250existing unit10-
G B 100existing unit14-
G C 1503,000,000162
G D 1102,500,00083
Table 2. Investment data for the transmission system.
Table 2. Investment data for the transmission system.
#( from ) to C km X M
($)
L I 1210,000,0001
L I I 246,000,0002
L I I I 345,000,0003
Table 3. Generalized shift-factors for the power system example.
Table 3. Generalized shift-factors for the power system example.
Line bus 1 bus 2 bus 3 bus 4
1–20−0.420−0.317−0.360
1–30−0.159−0.366−0.283
2–400.034−0.138−0.069
2–400.062−0.248−0.324
3–40−0.0410.166−0.117
Table 4. Generator data for the Garver system.
Table 4. Generator data for the Garver system.
#Bus P g M C g OM g FC g
(MW)($/MW)($/MW)($/MWh)
G 1 13 × 30existing unit-14.08
G 2 11 × 60existing unit-22.11
G 3 32 × 60existing unit-25.95
G 4 32 × 120300,000900020.41
G 5 61 × 120250,000750025.95
G 6 62 × 240350,00010,50014.08
G 7 G 9 2,4,5999virtual unit- V o L L
Table 5. Simulation results solving the Garver power system.
Table 5. Simulation results solving the Garver power system.
Modelv Θ nZ OF t s  (s)
GSF2361238213475,809,470.910.030
CD269132778475,809,470.910.021
Table 6. Simulation data applied to the Garver system using the multi-stage analysis.
Table 6. Simulation data applied to the Garver system using the multi-stage analysis.
ModelEv Θ nZ OF t s  (s)
GSF156,90816,8951,959,394621,509,628.5979.55
CD155,93818,830187,419621,509,628.592.39
GSF5282,35582,1099,790,325625,989,931.68793
CD5277,54591,790932,616625,989,931.61043
GSF10564,153163,62419,578,388629,964,367.335,692
CD10554,543182,9901,863,994629,964,367.39801
Table 7. Simulation results obtained for the IEEE 300-bus power system.
Table 7. Simulation results obtained for the IEEE 300-bus power system.
ModelEv Θ OF t s  (s)
GSF143,470312,90052,024,910,450.44032
CD1214,110627,69052,024,910,450.487,797
GSF5217,3501,567,500257,487,249,145.5272,240
CD51,070,5503,141,450--

Share and Cite

MDPI and ACS Style

Hinojosa, V.H.; Sepúlveda, J. Solving the Stochastic Generation and Transmission Capacity Planning Problem Applied to Large-Scale Power Systems Using Generalized Shift-Factors. Energies 2020, 13, 3327. https://0-doi-org.brum.beds.ac.uk/10.3390/en13133327

AMA Style

Hinojosa VH, Sepúlveda J. Solving the Stochastic Generation and Transmission Capacity Planning Problem Applied to Large-Scale Power Systems Using Generalized Shift-Factors. Energies. 2020; 13(13):3327. https://0-doi-org.brum.beds.ac.uk/10.3390/en13133327

Chicago/Turabian Style

Hinojosa, Victor H., and Joaquín Sepúlveda. 2020. "Solving the Stochastic Generation and Transmission Capacity Planning Problem Applied to Large-Scale Power Systems Using Generalized Shift-Factors" Energies 13, no. 13: 3327. https://0-doi-org.brum.beds.ac.uk/10.3390/en13133327

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop