Abstract
The penetration rate of intermittent renewable energy in power system is gradually increasing, which brings challenges to the optimal dispatch of power system. To solve these problems, a multi-regional interconnected transmission network optimization method based on analytical target cascading is proposed. Firstly, the reactive power regulation characteristics of renewable energy is investigated, and the models of wind turbine and photovoltaic generation are established. Secondly, the power system is decomposed into multiple sub-systems using bus tearing method, and an optimal transmission switching (OTS) model with renewable energies is established. Then, the analytical target cascading (ATC) approach is employed to decompose the model into the main problem and several sub-problems for parallel computation to achieve coordinated optimization of the complex transmission network. Finally, the IEEE 14-bus system and the IEEE 118-bus system are used to verify the proposed model. The results show that the proposed method can deal with the coupling nonlinear problem well and promote the consumption of renewable energy.
1 Introduction
In the context of the Energy Internet, the interconnections of the power grid facilitate the long-distance transmission of electrical energy and the improvement of energy utilization efficiency, which in turn supports the optimal management of high-capacity and high-efficiency units. However, the increasing contradiction of source–load imbalance and various intermittent renewable energy sources account for an increasing proportion of power generation, which makes it difficult to predict the uncertainties at the grid side and the contradictions of the source network become prominent (). In addition, information transmission between dispatching centers of different transmission networks is often inaccurate and untimely, which brings some difficulties to the development of the power transmission plan of boundary nodes, network blockage, and other problems that often occur during transmission grid operation (). Therefore, it is of great significance to enhance interconnection and non-homology effects between regions of the power system.
The non-synchronous phenomenon in the power transmission network can be alleviated to a certain extent by optimal transmission switching (OTS) through line cut-off and open-loop operation of the electromagnetic loop network on the strength of the practical working conditions of the power network, without additional investment in power transmission network construction. This can eliminate network congestion, alleviate conflicts between sources, networks, and loads, and improve the flexibility and economy of the system. Based on this, OTS has become an effective method to change the tide distribution by changing the opening and closing states of some lines in the power system (). Some scholars have studied the OTS, and in ), the OTS based on the generalized short-circuit ratio sensitivity analysis was proposed, which helps reconcile the conflict between system strength and short-circuit current levels. In ), an SSDS system framework of the UHV power grid considering the risk of the communication system was proposed, which is profitable in improving the stability of the UHV system. The economic scheduling model considering OTS involves mixed-integer programing, which introduces a significant number of discrete variables into the model, resulting in increased complexity when solving the model. While previous studies have demonstrated the effectiveness of OTS in improving system economy, there is a lack of research on optimizing the solving speed of the model.
In order to improve the efficiency of solving, many scholars have carried out thorough research around the solving algorithm. In ), the distributed optimal power flow (OPF) calculation for alternating current systems was achieved by using the alternating direction method of multipliers (ADMM). The ADMM is cannot eliminate the drawbacks of the first-order algorithm because it is based on the augmented Lagrange multiplier and the proximal point algorithm (). The distributed interior point algorithm was applied in ) to decentralize the solution of the OPF of the grid and to achieve simultaneous iterations in each region. In ), the synchronous ADMM was used to solve the coordinated optimization of numerous sub-systems, and then the overall optimization was determined. In ), based on the research of vulnerability analysis of critical nodes in the complex network, a cascading fault model was established. In ), an energy trading model based on stochastic programing (SP) was established, and it used distributed alternating search procedures to accelerate the calculation process of the Nash equilibrium. The aforementioned research optimization algorithms can effectively improve the speed of the model solution but do not consider the complexity of the power system network, which is affected by zonal management and non-sharing of information, thus making it hard to set up a general solution of multi-regional interconnected transmission.
Analytical target cascading (ATC) is a distributed algorithm suitable for solving interconnection and coupling problems. It does not depend on additional information input, can be solved efficiently and accurately only by decoupling variables of the initial problem, and can be used to handle the large-scale renewable energy grid of entropy increase, accelerating the use of centralized ways to work out the majorization of the large-scale optimization process. In ), ATC is applied to seek the optimal plan for day-ahead scheduling. During iteration, only the voltage and phase angle of the boundary are shared between the main problem and sub-problems. As a result, the amount of data exchanged per iteration is reduced. In ), a double-deck scheduling solution was established by ATC, and then, a linearized optimization model for active distribution networks is established due to existence of interaction variables between the upper and lower layers. The aforementioned research provides a good reference for this work to conduct OTS.
In addition, the large-scale grid integration of renewable energy gradually changes the balance pattern of load generation in the traditional grid, and its related research discusses methods of optimizing the consumption rate of renewable energy in the power system (; ; ). A large number of studies on OTS that consider renewable energy sources have been conducted. In ), an OTS model with renewable energy considering structural optimization was proposed, which reduces the system operation cost and cost of loss due to network blockage after failure by OTS. In ), a two-stage stochastic optimization model that considers OTS and energy storage allocation was constructed to improve the economy of the power grid by incorporating renewable energy sources. In ), to reduce the fluctuation influence of large-scale wind power in the grid, an OTS model for integration of large-scale renewable energy based on probabilistic power flow was proposed, and the superiority of which has been verified by numerical examples. However, the aforementioned studies do not consider the reactive voltage characteristics of wind and photovoltaic generators, cannot realize the deep excavation of renewable energy consumption capacity of the power system, and have certain limitations. As fossil energy generation gradually gives way to renewable energy generation, on the one hand, the active power balance gets worse regularly and absorption capacity of the renewable energy becomes insufficient. On the other hand, the proportion of synchronous generators with active excitation regulation in the system decreases, while the asynchronous power generation forms such as wind power and photovoltaic that consist of voltage support increase. Therefore, reactive power voltage regulation resources in the system should be fully tapped to support the consumption of renewable energy. In ), an optimization model based on the fruit fly optimization technique (FFOT) is established, which solves the coupling problem of output adjustment of multi-area and multi-renewable energy power. The advantages of ATC over the aforementioned research are that it allows for a systematic step-by-step approach to address the interdependencies between different components or sub-systems. By breaking down the problem into smaller and more manageable goals, it becomes easier to identify and address the specific coupling issues. This method also promotes modularity and flexibility, as changes or updates in an area can be made without affecting the entire system. Ultimately, the goal-oriented cascade method helps to reduce complexity and improve overall system performance.
To tap the renewable energy consumption potential of multi-regional interconnected power systems and design a model that conforms to the dispatching mode, an OTS model with renewable energy based on ATC is established. The main contributions are as follows:
1) The ATC-based model can effectively enhance the economic efficiency of operating interconnected grids, fully leverage the system’s flexible operational capabilities, and optimize resource utilization across the entire network.
2) An OTS model considering wind and photovoltaic energy is proposed, and the reactive voltage regulation characteristics of wind and photovoltaic are considered to measure uncertainty caused by the change of entropy, and then the renewable energy generation component model is established. By incorporating binary variables that indicate the status of the lines into the decision variables, the operation economy of power system with wind and photovoltaic generations can be improved by optimizing the power transmission network structure and the power output of the generators.
3) For the model of multi-regional interconnected transmission, the bus-tearing method is used to divide the electrical grids into separate zones, and independent economic dispatch is performed within each zone. ATC is used to solve the proposed model, which can partly reduce the difficulty of calculation. The applicability of the proposed model is verified based on numerical studies.
2 Renewable energy generation model
To deal with the uncertainty caused by the entropy increase of the interconnected grid such as renewable energy, it is crucial to take the development status and output characteristics of renewable energy generator components into account, mainly wind generators and photovoltaic, and realize their fine modeling, so as to provide effective guarantee for the resource allocation of OTS. Section 2 considers the characteristics of the current power system application technology development to build renewable energy generator models.
2.1 Doubly fed induction wind turbine
One of the key advantages of the doubly fed wind turbine is that its rotor winding has an additional power supply that can adjust the amplitude and phase angle, which can carry out the regulation of reactive power and emit reactive power. At the same time, the control of active and reactive power can be separated by regulating the excitation current of the rotor. The operating limit of the doubly fed induction wind generator is shown in Figure 1, where represents the maximum active power output by the doubly fed induction wind generator under a certain wind speed. The area surrounded by ABCD in the figure is the safe operation area of the unit. Similar to the synchronous generator, the regulation radius of the doubly fed induction wind turbine gradually increases with the active power output and decreases with the reactive power regulation range.
FIGURE 1
The doubly fed induction wind turbine is a popular model, as it uses a double PWM fully controlled converter to regulate the excitation current, allowing for variable speed and constant frequency operation. This type of wind turbine can adjust its active and reactive power output by modifying the amplitude and phase of the excitation current. Its power characteristics can be summarized as follows:
2.2 Photovoltaic generation
Photovoltaic generation output limits are shown in Figure 2. OA and OB correspond to the power factor constraint, and the photovoltaic generation can operate safely in the area surrounded by OABC within the condition of satisfying the maximum active output constraint. Under a certain active output, the reactive power regulation capability depends on the system maximum capacity. The grid-connected photovoltaic generator system has a poor reactive voltage support capacity, and the system has a higher possibility of voltage crossing limits and voltage instability.
FIGURE 2
The output limits of photovoltaic generation can be described as
3 Analytical target cascading
In ATC, the coupled information in the upper system is first passed to the sub-systems, and then each sub-system is solved separately to acquire the optimal solution (). Similar to the ADMM and auxiliary problem principle (APP) (), APP and ADMM methods utilize duality and penalty functions to decompose the primary optimization problem into multiple sub-problems. While the entire system is decomposed into multiple sub-systems by ATC, the constraint relaxation is utilized to enable information interconnection between sub-systems, and each problem is solved in parallel ().
To ensure the effectiveness of the amount of information shared between interconnected regions, in this work, the inter-regional tie line exchange power and voltage phase angle are selected as the shared information variables, and the shared information between adjacent layers is modeled by setting the target and response variables to form the objective function and constraints associated with each independent region. The objective variable is the shared information transmitted from the upper layer (main problem layer) to the lower layer (sub-problem layer), and response variables are the shared information transmitted from the lower layer to upper layer.
4 OPF model with renewable energy based on ATC
Considering that the economic dispatch of renewable energy and OPF is a non-linear optimization of complex networks, for this reason, the OTS model is established based on the bus-tearing method and ATC. Without loss of generality, this work divides the system into three parts by the bus-tearing method and builds a corresponding OPF model. This section takes region A as an example to construct a multi-regional economic dispatch model of the transmission grid based on ATC, and similarly, the models of region B and region C can be established.
4.1 Objective function
The objective function is to minimize the operation cost of the generators.
4.2 Constraints
(1) Power flow constraints
(2) Transmission line current constraints
(3) Generator output constraints
(4) Branch power constraints
(5) Voltage phase angle constraints
(6) Node power balance constraints
(7) Region-coupled constraints
In addition, constraints (1)–(8) related to wind and photovoltaic generations are also included in the proposed OPF model.
4.3 Model simplification
The proposed model is decomposed by taking the upper and lower systems shown in Figure 3. As seen in Figure 3A, to pursue the minimization of the total system operation cost, the upper–lower system information transfer variables are optimized iteratively. As shown in Figure 3, to illustrate the solution process of the proposed model, the relationship of each area in ATC is depicted with area A as the upper system.
FIGURE 3
To facilitate the description of ATC, the matrix function of the model is established as follows:
5 Model transformation and solution
5.1 Objective function
To facilitate the formation of objective functions and constraints about each region, two sets of variables were set. According to the structural characteristics of ATC, constraint (29) is solved in the upper and lower system models.
The optimization model of the upper–lower system can be described as follows:
By relaxing the coupling constraints through (30)–(32), to ensure the convergence of the upper and lower layer problems, the penalty term which includes the augmented Lagrangian function is added to the relevant objective function, by which, only the regional constraints and regional decision variables that have to be satisfied in the regional variables remain in the whole system, thus realizing the decoupling process.
5.2 Decoupled model
The objective function in optimal scheduling of the interconnected grid based on ATC can be described as follows:
When region A of the upper layer solves its own model, it transfers the optimized value of virtual load to the lower layer system in the form of parameters. When updating the upper layer system, it is necessary to optimize the linkage of virtual load and optimize the virtual generators in each area.
It can be seen from (30) that similar to the upper system, when the lower system is independently optimized, it is necessary to combine virtual generator and virtual load for optimization, relax the coupling constraint by introducing penalty function, and add it into the objective function of the lower system. The objective function of the lower system in region B can be described as follows:
By updating multipliers, the objective function meets the convergence condition and tends to be optimal. At this point, the upper and lower systems can be updated and solved independently.
5.3 Convergence criterion
The convergence criterion of the proposed model based on ATC is as follows:
If both (35) and (36) cannot be satisfied, the multipliers are updated according to (37).
5.4 Optimization process of ATC
The optimization processes of the transmission network structure based on ATC are shown in Figure 4. The approach steps are as follows:
FIGURE 4
Step 1Set the iteration limit of the ATC model as , set current iteration , input the thermal generator parameters, and initialize the values of each optimization variables and penalty function multipliers.
Step 2Solve the lower system variables. According to (28) and (34), parallel optimization is conducted for each region, and the power obtained from the solution is transferred to the upper system.
Step 3Solve the upper layer system variables. Optimization is conducted according to (27) and (33), and the power obtained is transferred to the lower system.
Step 4Determine whether (35) and (36) are satisfied. If not, evolve according to the update policy with regard to (37), set the number of iterations , and finally, return to Step 2; or else, terminate the iteration and output the results. If is greater than the maximum number of iterations of the ATC model, record that this iteration fails to converge and terminate the calculation.
6 Numerical studies
The improved IEEE 14-bus and IEEE 118-bus systems are used to simulate the model, and optimization of the results are analyzed. The computer is configured with Intel i7-6300 processor and 8 G memory. The GAMS software is used for programing, and the CPLEX tool is used to solve the model. The detailed line parameters are shown in Tables 1–5.
TABLE 1
| Unit | Upper limit of output (MW) | Lower limit of output (MW) | Generator cost factor ($/MWh) |
|---|---|---|---|
| GA1 | 285 | 0 | 1.06 |
| GA2 | 90 | 0 | 5.25 |
| GA3 | 85 | 0 | 3.12 |
| GB1 | 150 | 0 | 1.724 |
| GB2 | 285 | 0 | 2.011 |
| GC1 | 200 | 0 | 3.561 |
IEEE 14-bus generator parameters.
TABLE 2
| Sub-region | Transmission line | Reactance (p.u.) | Transmission capacity (MW) | Sub-region | Transmission line | Reactance (p.u.) | Transmission capacity (MW) |
|---|---|---|---|---|---|---|---|
| B | 1–2 | 0.1739 | 150 | A | 1–2 | 0.0592 | 80 |
| 1–3 | 0.171 | 200 | 1–101 | 0.223 | 70 | ||
| 2–4 | 0.0421 | 150 | 2–101 | 0.198 | 80 | ||
| 3–5 | 0.2091 | 70 | 101–201 | 0.1763 | 150 | ||
| 3–6 | 0.5562 | 150 | C | 1–3 | 0.1558 | 150 | |
| 301–1 | 0.252 | 150 | 2–3 | 0.1303 | 60 | ||
| 102–301 | 0.1989 | 150 | 202–1 | 0.1762 | 60 | ||
| - | - | - | 202–2 | 0.11 | 150 |
IEEE 14-bus line parameters.
TABLE 3
| Sub-region | Node number | Load (MW) | Sub-region | Node number | Load (MW) |
|---|---|---|---|---|---|
| B | 1 | 80 | A | 1 | 150 |
| 2 | 24 | 2 | 60 | ||
| 3 | 24 | 101 | 70 | ||
| 4 | 60 | 201 | 60 | ||
| 5 | 30 | C | 1 | 30 | |
| 6 | 30 | 2 | 50 | ||
| 102 | 50 | 3 | 40 | ||
| 301 | 50 | 202 | 60 |
IEEE 14-bus load parameters.
TABLE 4
| Node number | Region A (Node 2) | Region B (Node 2) | Region C (Node 2) |
|---|---|---|---|
| Maximum output (MW) | 60 | 50 | 60 |
| Minimum output (MW) | 0 | 0 | 0 |
Wind generator parameters.
TABLE 5
| Node number | Region A (Node 1) | Region B (Node 4) | Region C (Node 3) |
|---|---|---|---|
| Maximum output (MW) | 60 | 50 | 60 |
| Minimum output (MW) | 0 | 0 | 0 |
Photovoltaic power generator parameters.
6.1 IEEE 14-bus system
There are five generators and 20 lines in the IEEE 14-bus system, and the other parameters are given in Appendix A. Considering the randomness of photovoltaic and wind power generation, the random variables in the upper and lower limits of their output are selected for testing. The system is decomposed into a three-zone system by the busbar-tearing method, as shown in Figure 5. In addition, the number of open power lines in region A and region B is 1.
FIGURE 5
6.1.1 Algorithm effectiveness
The results obtained aim to verify the solving advantages of the adopted ATC algorithm in comparison with the APP algorithm and augmented Lagrange relaxation (ALR) algorithm, as shown in Table 2 and Figure 6.
FIGURE 6
As shown in Table 6 and Figure 6, the operating costs of the solution using the APP algorithm and ALR algorithm are $18,212.23 and $18,103.25, respectively, while the operating cost obtained by the ATC algorithm is $18,394.46. The reason for the operating cost of the ATC algorithm being a little higher than that of the APP algorithm and ALR algorithm is because of the setting of the initial parameters and the change in the network topology. The number of iterations and convergence time of the ATC algorithm are less than that of the APP algorithm and ALR algorithm, which are 34 times and 7.71 s, respectively, indicating that the ATC algorithm has a better convergence performance. In addition, the complexity of the OTS problem itself makes it difficult for the APP algorithm and ALR algorithm to obtain a high-precision solution. It is evident that the ATC algorithm is highly scalable and can be well adapted to the OTS optimization problem of the system.
TABLE 6
| Algorithm | Operating cost ($) | Number of iterations | Iteration time (s) |
|---|---|---|---|
| ATC | 18,394.46 | 34 | 7.71 |
| APP | 18,212.23 | 43 | 8.97 |
| ALR | 18,103.25 | 35 | 10.31 |
Results of different algorithms.
6.1.2 OTS and renewable energy access impact analysis
In this study, the following four cases are set up and compared for analysis in order to demonstrate the effect of OTS and the impact of wind power and photovoltaic power generation on system operation.
Case 1No consideration of OTS and renewable energy generation.
Case 2Without considering OTS and considering renewable energy generation.
Case 3Considering OTS without considering renewable energy generation.
Case 4
Considering both OTS and renewable energy generation.
1) Analysis of the impact of not considering OTS
From
Table 7, due to the consideration of renewable energy generation access, case 2 greatly releases the generation capacity of GA1, which increases the output of GA1 with a better economy and reduces the output of GA2 with the worst economy. The cost of each region in case 1 is $8,476.62, $7,363.38, and $4,363.69. The operating cost of each region in case 2 is reduced by 19.97%, 1.39%, and −0.13% when compared with case 1, where the cost of region A is reduced more, while the cost of region B and region C remain less changed. This indicates that the lower region B and region C can achieve distribution autonomy by virtue of their own generation resources, thus greatly reducing the power relief of the upper region A and making the economy of region A significantly improved.
2) Impact analysis of OTS
TABLE 7
| Region | Units | Case 1 | Case 2 | ||
|---|---|---|---|---|---|
| Cost ($) | Output power (MW) | Cost ($) | Output power (MW) | ||
| Region A | GA1 | 8,476.62 | 86.01 | 6,783.88 | 118.65 |
| GA2 | 79.34 | 33.44 | |||
| GA3 | 85.00 | 85.00 | |||
| Region B | GB1 | 7,363.38 | 63.00 | 7,261.27 | 63.00 |
| GB2 | 185.00 | 185.00 | |||
| Region C | GC1 | 4,363.69 | 126.41 | 4,369.31 | 129.96 |
Operating results of case 1 and case 2.
TABLE 8
| Regions | Units | Case 1 | Case 2 | ||
|---|---|---|---|---|---|
| Cost ($) | Output power (MW) | Cost ($) | Output power (MW) | ||
| Region A | GA1 | 6,783.88 | 118.65 | 5,814.34 | 169.99 |
| GA2 | 33.44 | 21.64 | |||
| GA3 | 85.00 | 85.00 | |||
| Region B | GB1 | 7,261.27 | 63.00 | 6,086.17 | 69.00 |
| GB2 | 185.00 | 166.00 | |||
| Region C | GC1 | 4,369.31 | 129.96 | 4,261.58 | 121.00 |
Operating results of case 2 and case 4.
TABLE 9
| Region | Transmission line | Transmission power (MW) | ||
|---|---|---|---|---|
| Case 2 | Case 3 | Case 4 | ||
| Region A | 1-2 | 67.34 | 58.88 | 89.49 |
| 1-101 | 41.65 | 34.90 | 0.00 | |
| 2-101 | 3.35 | 0.00 | 28.59 | |
| 101-201 | 56.72 | 61.65 | 56.38 | |
| Region B | 1-2 | 84.00 | 77.00 | 77.00 |
| 1-3 | −150.00 | −110.00 | −110.00 | |
| 2-4 | 60.00 | 54.00 | 54.00 | |
| 3-5 | 70.00 | 70.00 | 70.00 | |
| 3-6 | 60.00 | 54.00 | 54.00 | |
| 301-1 | −86.00 | −91.00 | −91.00 | |
| 102-301 | −40.66 | −47.33 | −43.63 | |
| Region C | 1-3 | 90.00 | 68.77 | 93.89 |
| 2-3 | −50.00 | −57.58 | −59.15 | |
| 202-1 | −26.41 | −31.14 | −32.40 | |
| 202-2 | 0.00 | 0.00 | 0.00 | |
Transmission line operation considering OTS.
6.2 IEEE 118-bus system
The IEEE 118-bus system structure diagram and partitions are shown in Figure 7 in the Appendix, and the system generator set characteristics data, transmission line data, and system load data are shown in ).
FIGURE 7
6.2.1 Algorithm performance comparison
In the IEEE 118-bus system, ATC, APP, and ALR are compared in order to testify ATC in the large-scale system. The results are shown in Figure 8.
FIGURE 8
As shown in Figure 8, in the wake of improvement in convergence accuracy, ATC tends to the optimum at the fastest speed, showing good convergence performance. In addition, the number of iterations increases with increase of convergence accuracy. When compared with APP and ALR, ATC always maintains a good convergence performance and optimization results under any accuracy requirement. This shows that the ATC still has good applicability to large-scale systems.
6.2.2 Comparison of optimization results
For the IEEE 118-bus system, the aforementioned four schemes are used to test the proposed model, and the results are shown as follows:
Table 10 shows that the total cost of case 4 decreases by 18.99%, 17.65%, and 2.50% when compared with cases 1–3, respectively, and the operating costs of the three regions are $47,311.34, $63,384.13, and $52,117.87, all of which have different degrees of reduction when compared to cases 1–3. This shows that OTS can not only ensure the safe operation of the system but also improve the acceptance level of renewable energy and economy of system operation. The operating cost of region A and region B in case 4 decreases more than that of region C when compared with cases 1–3, which is due to the more serious blockage in region A and region B, so the cost saving accounts for the more obvious. In summary, the simultaneous application of renewable energy and open transmission lines in the grid to realize the cooperation between system regions can reduce the cost of power generation and heighten the utilization efficiency of electric energy, thus enhancing flexibility and economy of overall system operation.
TABLE 10
| Case | Cost ($) | |||
|---|---|---|---|---|
| Region A | Region B | Region C | In total | |
| Case 1 | 64,348.42 | 77,384.13 | 59,244.67 | 200,977.22 |
| Case 2 | 61,694.56 | 76,478.54 | 58,536.23 | 197,709.33 |
| Case 3 | 48,324.48 | 66,536.32 | 52,135.30 | 166,996.10 |
| Case 4 | 47,311.34 | 63,384.13 | 52,117.87 | 162,813.34 |
Results of different cases.
7 Conclusion
In this study, research was conducted on the optimization problem of multi-regional complex network system dispatching, and an OTS model based on the ATC that was proposed to account for renewable energy. The conclusions are as follows:
1) The ATC adopted in this study has a fast convergence speed, which has notable efficiency advantages when applied to the distributed solution of multi-regional interconnected transmission. It can also realize the decoupling of interconnected information, which not only ensures the privacy requirement of multi-agent operation mode but also has good adaptability to large-scale systems. The ATC-based model can effectively improve the economy of the operation of the interconnected grid, bring the flexible operation capability of the system into full play, and efficiently utilize the resources of the whole network.
2) The proposed model adopts OTS in the system containing renewable energy generation, improves the distribution of power flow in transmission networks to optimize the entropies of complex networks, and actively cooperates with renewable energy generators that do not have controllability to realize source–grid co-optimization, which is a proven solution to improve renewable energy consumption capacity and system voltage level.
3) The proposed model adds discrete variables, indicating whether the transmission lines are operating or not, which can be transferred to the dispatch center as decidable information, making it possible to reasonably decide the required transmission grid architecture of the system. The combination of grid structure optimization, renewable energy, and ATC greatly frees up the transmission capacity of the grid and improves the generality of the model.
Future research can focus on enhancing the computational efficiency and accuracy of the model.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, editors, and reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Nomenclature
| Sets | |
| Inequality constraint vectors of the system | |
| Equation constraint vectors of the system | |
| Inequality constraint vectors of the sub-system | |
| Equation constraint vectors of the sub-system | |
| Thermal generators in region A | |
| Thermal generators in region B | |
| Thermal generators in region C | |
| Transmission lines | |
| Time periods | |
| Branches of region A | |
| Thermal generators of the grid in region A | |
| Wind generators of the grid in region A | |
| Photovoltaic generators of the grid in region A | |
| Branch nodes of the grid in region A | |
| Loads of the grid in region A | |
| Adjacent regions of region A | |
| Tie lines of region A and B | |
| First and last nodes of tie lines of regions A and B | |
| Inequality constraints to be satisfied by region A in the upper system | |
| Equation constraints to be satisfied by region A in the upper system | |
| Inequality constraints of the lower system n | |
| Equality constraints of the lower system n | |
| Coupling constraints between upper system region variables and lower system n region variables | |
| Local variables in region A | |
| Local variables in region B | |
| Local variables in region C | |
| Maximum active power output by the doubly fed generator w | |
| Stator reactance | |
| Excitation reactance | |
| Upper limits of the slew rate | |
| Lower limits of the slew rate | |
| Maximum values of the stator currents | |
| Maximum values of the rotor currents | |
| Maximum active power output by photovoltaic generator | |
| Maximum current output by photovoltaic generator | |
| Minimum power factor allowed | |
| Positive weight coefficients for the response variables | |
| Positive weight coefficients for the coupling variables | |
| g | Index of generators |
| i | Index of first nodes of line l |
| j | Index of last nodes of line l |
| l | Index of line |
| t | Index of time |
| Conductance parameters of line l | |
| Susceptance parameters of line l | |
| n | Minimum technical output power of generator g |
| Maximum allowable thermal current of transmission line l | |
| Upper limits of the output power of generator g in region A | |
| Lower limits of the output power of generator g in region A | |
| Ramp-up rates of generator g | |
| Ramp-down rates of generator g | |
| Allowable ramp time of the generator | |
| Susceptance of branch l in region A | |
| A very large constant and | |
| Maximum transmission capacity of line l in region A | |
| Maximum voltage phase angle of node i | |
| Minimum voltage phase angle of node i | |
| Maximum allowed number of line openings in region A | |
| Hadamard product | |
| Active power output by wind generator w | |
| Reactive power output by wind generator w | |
| Voltage at the machine end of wind generator w | |
| Slew rate | |
| Active power of the photovoltaic generator | |
| Reactive power of the photovoltaic generator | |
| Output voltage of the photovoltaic system | |
| Upper system design variables | |
| Upper system responses | |
| Deviation between the objective and response of system | |
| Euclidean parametrization, which can be chosen to calculate the difference between the objective and response | |
| Deviation calculated using 2-parametrization | |
| Sub-system design variables | |
| Cost functions of generator g in region A | |
| Cost functions of generator g in region B | |
| Cost functions of generator g in region C | |
| Generator g output power in region A | |
| Generator g output power in region B | |
| Generator g output power in region C | |
| Binary variable about the operating status of line l at time t | |
| Voltage amplitude of node d in time t | |
| Voltage phase angle of node d in time t | |
| Active power of line l in time t | |
| Reactive power of line l in time t | |
| Amplitude of the current flowing on the transmission line l in time t | |
| Voltage phase angle of node i in region A | |
| Voltage phase angle of node j in region A | |
| Binary variable corresponding to the state of branch l in region A | |
| Active power of branch l in region A | |
| Values of wind generator output in region A | |
| Values of photovoltaic generator output in region A | |
| Active load of node d in region A | |
| Active power in the positive direction of branches | |
| Active power in the negative direction of branches | |
| Power differences of tie lines | |
| Phase angle differences of tie lines | |
| Decision variables of the upper layer except tie line | |
| Regional variables of the upper layer | |
| Regional constraints to be satisfied by the variables of the lower system n | |
| Regional variables of the lower system n | |
| Multipliers of first terms of penalty function | |
| Multipliers of second terms of penalty function | |
| Penalty function about relaxing coupling constraints | |
| Multipliers of the Lagrangian primary terms | |
| Multipliers of the Lagrangian second terms | |
| Virtual generator of the lower system | |
| Virtual load of the upper system |
References
1
Ahmed NawaZ.Shah ValareR. A. Z. M. A. O. (2021). Optimal transmission switching based on probabilistic load flow in power system with large-scale renewable energy integration. Electr. Eng.104, 883–898. 10.1007/s00202-021-01344-z
2
BaiY.ZhongH. W.XiaQ.KangC.XieL. (2015). A decomposition method for network-constrained unit commitment with ac power flow con-straints. Energy, Aug.88, 595–603. 10.1016/j.energy.2015.05.082
3
DashtdarM.BelkhierY.BaiaiM. “Protection of DC microgrids based on frequency domain analysis using fourier transform,” in Proceedings of the 2022 IEEE 3rd KhPI Week on Advanced Technology (KhPIWeek), Kharkiv, Ukraine, October 2022, 1–6.
4
EresghsT. (2014). Distributed optimal power flow using ADMM. IEEE Trans. Power Syst.29 (5), 2370–2380. 10.1109/tpwrs.2014.2306495
5
GaoX.PengM. F.TseC. K. (2023). Cascading Failure analysis of cyber–physical power systems considering routing strategy. IEEE Trans. Circuits Syst. II Express Briefs70 (1), 136–140. 10.1109/tcsii.2021.3071920
6
HuangL. B.XinH. H.LiZ. Y.LiZ.JuP.YuanH.et al (2020). Identification of generalized short-circuit ratio for on-line stability monitoring of wind farms. IEEE Trans. Power Syst.35 (4), 3282–3285. 10.1109/tpwrs.2020.2975413
7
JiX. Q.ZhangY. M.HanX. S.YeP.XuB.YuY. (2021). Multi-level interactive unit commitment of regional power system. Int. J. Electr. Power and Energy Syst.125, 106464. 10.1016/j.ijepes.2020.106464
8
KargarianA.MehrtashFalahatiM. B. (2017). Decentralized implementation of unit commitment with analytical target cascading: A parallel approach. IEEE Trans. Power Syst.33 (4), 3981–3993. 10.1109/tpwrs.2017.2787645
9
KhanabadiM.GhasemiH.DoostizadehM. (2013). Optimal transmission switching considering voltage security and n-1 contingency analysis. IEEE Trans. Power Syst.28 (1), 542–550. 10.1109/tpwrs.2012.2207464
10
LiH. X.HaoG. F.MaoM. X.ZhouL. (2022). Decentralized communication based two-tier voltvar control strategy for largescale cen-tralized photovoltaic power plant. IEEE Trans. Sustain. Energy13 (1), 592–606. 10.1109/tste.2021.3121536
11
LiZ. M.WuL.XuY. (2021). Risk-averse coordinated operation of a multi-energy microgrid considering voltage/var control and thermal flow: An adaptive stochastic approach. IEEE Trans. Smart Grid. Early access12 (5), 3914–3927. 10.1109/tsg.2021.3080312
12
LiZ.WuL.XuY. (2023). Distributed tri-layer risk-averse stochastic game approach for energy trading among multi-energy microgrids. Appl. Energy331. 10.1016/j.apenergy.2022.120282
13
LiZ.XuY.FangS.WangY.ZhengX. (2019). Multiobjective coordinated energy dispatch and voyage scheduling for a multi-energy ship microgrid. IEEE Trans. Industrial Appl.56 (2), 989–999. 10.1109/tia.2019.2956720
14
LiuB.LiZ.ChenX.HuangY.LiuX.ChenX.et al (2017). Recognition and vulnerability analysis of key nodes in power grid based on complex network centrality. IEEE Trans. Circuits Syst. II Express Briefs65 (3), 346–350. 10.1109/tcsii.2017.2705482
15
LiuY. M.WangGuT. X. P. (2022). Robust optimal transmission switching for wind farm‐integrated power systems from a perspective of steady-state security region. IET Renew. Power Gener.16 (5), 945–955. 10.1049/rpg2.12404
16
Naga Sai KalyanC. H.Srikanth GoudB.NagendraM. “Performance enhancement of combined LFC and AVR system with the integration of HVDC line,” in Proceedings of the 2023 IEEE IAS Global Conference on Renewable Energy and Hydrogen Technologies (GlobConHT), Male, Maldives, March 2023, 1–6.
17
OgundairoO.KamalasadanS.NairA. R.SmithM. (2022). Oscillation damping of integrated transmission and distribution power grid with re-newables based on novel measurement-based optimal controller. IEEE Trans. Industry Appl.58 (3), 4181–4191. 10.1109/tia.2022.3162565
18
OstrowskiJ.WangJ.LiuC.LiuC. (2014). Transmission switching with connectivity-ensuring constraints. IEEE Trans. Power Syst.29 (6), 2621–2627. 10.1109/tpwrs.2014.2315434
19
SahriY.BelkhierY.TamalouztS.UllahN.ShawR. N.ChowdhuryM. S. (2021). Energy management system for hybrid PV/wind/battery/fuel cell in microgrid-based hydrogen and economical hybrid battery/super capacitor energy storage. Energies14 (18), 5722. 10.3390/en14185722
20
ShayestehE.HobbsB.SöderL.AmelinM. (2015). ATC-based system reduction for planning power systems with correlated wind and loads. IEEE Trans. Power Syst.30 (1), 429–438. 10.1109/tpwrs.2014.2326615
21
SrivastavaP.BajajM.BelkhierY. “Power quality issues with penetration of wind energy and the role of ldc: A mini review,” in Proceedings of the International Conference on Computational Intelligence and Sustainable Technologies: ICoCIST 2021, Singapore, June 2022Springer Nature Singapore, 625–635.
22
WangP.WuQ.HuangS.ZhouB. (2022). ADMM-based distributed active and reactive power control for regional ac power grid with wind farms. J. Mod. Power Syst. Clean Energy10 (3), 588–596. 10.35833/mpce.2020.000918
23
WeiY.WenL. L.LiW.ChungC.WongK. (2011). Decomposition–coordination interior point method and its application to multi-area optimal re-active power flow. Int. J. Electr. Power and Energy Syst.33 (1), 55–60. 10.1016/j.ijepes.2010.08.004
24
ZhaiJ. Y.ZhouM.LiJ. N.ZhangX. P.ZhangW. (2021). Decentralised and distributed day-ahead robust scheduling frameworks for bulk AC/DC hybrid interconnected systems with a high share of wind power. Electr. Power Syst. Res.201. 10.1016/j.epsr.2021.107492107492
Summary
Keywords
multi-regional interconnected transmission, analytical target cascading, renewable energy, optimal transmission switching, optimal power flow
Citation
Li R, Ye P, Jiang D, Ma L and Sun C (2023) Optimal transmission switching for power system integrating renewable energy based on analytical target cascading. Front. Energy Res. 11:1239232. doi: 10.3389/fenrg.2023.1239232
Received
13 June 2023
Accepted
17 July 2023
Published
27 July 2023
Volume
11 - 2023
Edited by
Yuqing Dong, The University of Tennessee, United States
Reviewed by
Lili Hao, Nanjing Tech University, China
Xiangjun Zeng, China Three Gorges University, China
Xingxu Zhu, Northeast Electric Power University, China
Youcef Belkhier, Maynooth University, Ireland
Updates
Copyright
© 2023 Li, Ye, Jiang, Ma and Sun.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pingfeng Ye, ypf@sdust.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.