Next Article in Journal
Design Analysis of Micro Gas Turbines in Closed Cycles
Previous Article in Journal
Simplified Assessment on the Wind Farm Noise Impact of the E2O Experimental Offshore Station in the Asturian Coast
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Decoupling Rolling Multi-Period Power and Voltage Optimization Strategy in Active Distribution Networks

1
Electric Power Research Institute of State Grid Zhejiang Electric Power Company, Hangzhou 310014, China
2
Department of Electrical Engineering, Southeast University, Nanjing 210096, China
3
State Grid Zhejiang Electric Power Corporation, Hangzhou 310007, China
*
Author to whom correspondence should be addressed.
Energies 2020, 13(21), 5789; https://doi.org/10.3390/en13215789
Submission received: 7 September 2020 / Revised: 28 October 2020 / Accepted: 1 November 2020 / Published: 5 November 2020
(This article belongs to the Section A1: Smart Grids and Microgrids)

Abstract

:
With the increasing penetration of distributed photovoltaics (PVs) in active distribution networks (ADNs), the risk of voltage violations caused by PV uncertainties is significantly exacerbated. Since the conventional voltage regulation strategy is limited by its discrete devices and delay, ADN operators allow PVs to participate in voltage optimization by controlling their power outputs and cooperating with traditional regulation devices. This paper proposes a decoupling rolling multi-period reactive power and voltage optimization strategy considering the strong time coupling between different devices. The mixed-integer voltage optimization model is first decomposed into a long-period master problem for on-load tap changer (OLTC) and multiple short-period subproblems for PV power by Benders decomposition algorithm. Then, based on the high-precision PV and load forecasts, the model predictive control (MPC) method is utilized to modify the independent subproblems into a series of subproblems that roll with the time window, achieving a smooth transition from the current state to the ideal state. The estimated voltage variation in the prediction horizon of MPC is calculated by a simplified discrete equation for OLTC tap and a linearized sensitivity matrix between power and voltage for fast computation. The feasibility of the proposed optimization strategy is demonstrated by performing simulations on a distribution test system.

1. Introduction

The rapid penetration of distributed photovoltaics (PVs) in active distribution networks (ADNs) has considerably challenged the operation mode of the traditional grid [1]. Some research indicates that one of the major constraints for promoting PV integration into ADNs is voltage stability [2,3]. The voltage at the end of feeder will rise if the intermittent production of PV generation causes power flow reversal [4,5]. Since undesirable voltage violations may cause serious negative effects [6], an appropriate effective reactive power and voltage strategy is required to solve these problems.
The on-load tap changer (OLTC) in distribution networks is considered to be a typical control device to prevent the violations of voltages, which keeps the substation secondary bus voltage constant by adjusting the tap position [7]. However, due to physical constraints and low cost benefit ratios, it cannot change flexibly to achieve real-time adjustment [8,9]. Recently, the continuous, quick and frequent reactive power outputs of distributed generation are appealing to mitigate voltage violations and fluctuations [10,11], which avoids adding extra voltage devices and reactive compensators in ADN. The authors of [12,13] presented an optimized decentralized control approach for each PV inverter based on voltage sensitivity analysis by implementing the control action at the point of common coupling (PCC) [14] as one of ancillary services to the ADN. A multi-mode adaptive local reactive power control method based on Q(P) characteristics is discussed in [15] considering control parameter optimization. Reference [16] proposed a combination of centralized-remote control and local-decentralized control to ensure that the nodal voltages remain within regulatory limits. Although these studies have fully demonstrated that distributed PVs have the ability to enhance power quality as well as reliability of the system [17], their reactive outputs are limited by capacity and power factor. Moreover, the line equivalent impedance directly affects the voltage regulation. the R/X ratio determines the influence of the active power and reactive power injected by distributed generation on the node voltage [18]. For the low-reactance cables, it is more difficult to try to keep the voltage within the allowable range by adjusting the reactive power output of distributed generations. The fluctuation of PV active power output greatly affects the node voltage. When serious violations occur, it is insufficient to correct serious violations by distributed PVs alone, which requires cooperation with traditional adjustment devices [19]. It follows that the reactive power and voltage optimization model of ADN is a complex mixed integer nonlinear programming (MINP) problem, not only dealing with continuous control variables such as of PV outputs, but also handling discrete variables such as OLTC tap changes. In [20,21,22], a control method for coordinating the operation of OLTC and distributed generations (DGs) based on the initial optimal power flow (OPF) calculation in order to maintain the load end voltage of the feeder within the allowed band. The authors of [23,24] incorporated a semi-definite power flow relaxation to tackle such complex problems, while [25,26] decomposed the problem into subproblems by Benders decomposition method and alternate direction method of multipliers respectively. In [27,28,29,30], the discrete variables in optimization models are relaxed into continuous variables and then rounded to close integers. However, these solutions are the static optimization results, regardless of the strong time coupling between the adjustments of different devices. The authors in [31] proposed a two-stage dynamic coordination dispatch method based on the Haar wavelet transform to adjust the OLTC tap position and the reactive power of distributed generation. [32] implemented another two-stage dynamic reactive power dispatch method, which first performed a day-ahead dispatch plan for switched capacitor banks and OLTCs by Heuristic search, and then recalibrated the intra-day dispatch results of PVs. These two-stage methods mentioned pre-optimize discrete devices in advance without considering the interaction between OLTC tap changes and PV outputs. Consequently, the pre-planned scheduling may have large errors and lose dynamic optimality.
Furthermore, previous research has drawn attention to the solutions of PV generation uncertainties. In fact, most optimal dispatch methods based on the predicted scenario can only work in shorter periods, since accurate weather and load forecasts are only valid for a limited period in advance. A multi-step optimization method based on the model predictive control (MPC) principle has been verified to be well suited for this issue by determining the appropriate prediction horizon and control horizon [33,34]. Solving the rolling optimization problems of specific horizons can achieve a dynamic transition from the measured to the target values, which is different from single step optimizations. Notably, varying the horizon of prediction and control has a significant impact on the optimal operation of the current time window. The authors in [35] developed a planning method that leverages MPC principle to decide the optimal location and size of storage devices in ADN. [36] proposed a control strategy based on MPC to balance the intermittency of PV outputs by little adjustments of conventional generation and load, as well as the storage devices. In [9], an MPC-based corrective control is utilized to correct voltages out of limits by solving the OLTC voltage set-point and PV outputs. However, due to the different operation characteristics of various devices, the requirement of subdividing the control horizon was ignored. On the one hand, limited forecasts cannot support the long-period scheduling. On the other hand, frequent short horizons will increase the modeling time and make the optimization problem computationally intractable. Table 1 gathers the main features of ADN optimization methods found in the existing literature.
In this paper, a decoupling rolling multi-period (DRM) reactive power and voltage optimization method of ADN is proposed. The main contribution is two-fold. First, for the mixed-integer reactive power and voltage optimization problem, we utilize Benders decomposition algorithm to decouple the coordination between discrete OLTC tap changers and continuous PV power outputs in the original problem into the master problem and the subproblem. In terms of time, our method designs the master problem of OLTC tap as a long-period optimization model, and the subproblem of PV power is a series of short-period optimization models. This allows us to address PV and load uncertainty by dividing the fluctuating PV and load condition into various short-period scenarios, avoiding frequent OLTC tap changes. Secondly, these short-period independent subproblems are incorporated into the MPC strategy and solved on a rolling basis, because predicting PV power generation of all short periods in the current time window may lead to large errors and additional adjustment costs. Based on the linearized voltage prediction model, the interrelation between rolling subproblems are considered to realize the gradual evolution of the ideal voltage in a long period, i.e., the prediction and control horizon of MPC. Besides, we present how to formulate Benders decomposition algorithm for this reactive power and voltage optimization problem, including the solution process, achieving the multi-period decoupling rolling optimization effect of different devices.
The rest of the paper is organized as follows: In Section 2, the model of voltage sensitivity is established, as well as the modeling of OLTC and PV. In Section 3, the proposed DRM architecture is illustrated. Section 4 formulates the decomposed optimization models and provides the solution process. In Section 5, the simulation results on the distribution test system are presented to validate the proposed strategy. The conclusions are drawn in Section 6.

2. Modeling of Active Distribution Network

2.1. Sensitivity between Node Voltage and Power

In order to discuss the influence of different voltage regulation devices on nodal voltages, the distribution network can be locally linearized to calculate the sensitivity model (the relationship between nodal voltage deviations and output changes of regulation devices) [37].
Assuming that a distribution network system contains n nodes, the voltage of Node i are expressed in polar coordinate form, V ˙ i = V i e j θ i = V i c o s θ i + j s i n θ i . The modified Newton Raphson power flow equation can be simplified as:
Δ P Δ Q = J Δ θ Δ V
where V is the diagonal matrix of voltage amplitude, V = diag(v1, v2,…, vn-1), Δ θ and Δ V are the corrections of the voltage phase angle and voltage amplitude respectively. The Jacobian matrix J is as follows:
J = H N K L = V V B c o s θ G c o s θ G c o s θ B c o s θ G c o s θ B s i n θ B c o s θ G s i n θ V V
where B c o s θ can be expressed as Equation (3) and the structures of other terms B sin θ , G c o s θ , and G s i n θ are similar.
B c o s θ = B 1 , 1 c o s θ 1 , 1 B 1 , j c o s θ 1 , j B 1 , j c o s θ 1 , n B i , 1 c o s θ i , 1 B i , j c o s θ i , j B i , n c o s θ i , n B n , 1 c o s θ n , 1 B n , j c o s θ n , j B n , n c o s θ n , n
In the distribution network, the value of θ i j is extremely small, and hence Equation (1) can be further simplified:
B + Q G P G P B Q V Δ θ Δ V = Δ P / V Δ Q / V
where P and Q are diagonal matrices, and their diagonal elements are P i / V i 2 and Q i / V i 2 respectively. Using the Gaussian elimination method for Equation (4), the sensitivity of nodal voltages and injected power can be calculated as:
Δ V = B + Q G P 1 B Q + G + P 1 Δ P
Δ V = G P B + Q 1 G + P + B + Q 1 Δ Q
Therefore, the active-voltage sensitivity and reactive-voltage sensitivity can be obtained:
S P = B + Q G P 1 B Q + G + P 1
S Q = G P B + Q 1 G + P + B + Q 1

2.2. Modeling of OLTC

The discussed OLTC acts on the voltage set-point of the secondary side by changing the tap position. The sensitivity V / V t a p between nodal voltages and OLTC set-point can be approximated from the solution of two power flow runs with a single tap position difference [9], as follows:
S t a p = V V t a p Δ V Δ V d
where ΔVd is voltage variation on the low-voltage side after a OLTC tap change.

2.3. Modeling of PV

In this paper, PVs are set to adjust the reactive power output in the maximum power point tracking (MPPT) mode [7]. The maximum reactive power is limited by active power, power factor and capacity:
Q m , P V t , m a x = P m , P V t t a n φ m a x , 0 < P m , P V t < P l Q m , P V t , m a x = S m , P V 2 P m , P V t 2 , P l < P m , P V t < P m a x
P l = S m , P V c o s φ m a x
where, P m , P V t is the PV-m active power output. Q m , P V t , m a x is the maximum reactive power. S m , P V is the capacity of the PV-m, and φ m a x is the maximum power factor angle.

3. Decoupling Rolling Multi-Period Architecture

In view of nodal voltage deviations caused by the high proportion of PV generation, the practical implementation is to adjust the OLTC tap positions and the active/reactive power outputs of distributed PVs. The OLTC tap changes are discrete, while the PV power outputs can continuously, quickly, efficiently and accurately respond to operation instructions. The DRM strategy focuses on the temporal interaction between PV and OLTC and explores their cooperation performance on voltage optimization.
The reactive power and voltage optimization model is a MINP problem. The Benders decomposition method is a distributed algorithm that decomposes the original MINP problem into the master problem and the subproblem to obtain the optimal operation results of different controllable devices in the same time step, as shown in Figure 1a. Due to the limitation of OLTC tap changes, the synchronized optimization framework cannot take full advantage of PV flexibility potential. The improved strategy in Figure 1b further decomposes the long-period subproblem of PV reactive power into multiple short-period problems, not only realizing rapid adjustment of PV reactive power, but also avoiding excessive and frequent OLTC tap changes. However, these parallel independent short-period subproblems of PV outputs are still based on the definite long-period PV power generation forecasts. Taking into account the uncertainty of PV power generation prediction, the prediction error becomes larger when the period is longer. The DRM method proposed in this paper rolls a series of short-period subproblems in the prediction horizon to smoothly adjust the PV power, as shown in Figure 1c.
The detailed DRM reactive power and voltage optimization strategy is shown in the Figure 2. The decomposed long-period master problem is to solve OLTC tap changes. Since the increasing number of tap actions will shorten the life of OLTC, the optimization function is to minimize the adjustment cost of OLTC and set at most one tap position change per long period to reduce frequent actions. The voltage is calculated by the sensitivity between the nodal voltage and the OLTC tap in Equation (9) to constrain the range of tap change. After the master problem is solved, the optimal solution of OLTC tap will be sent to subproblems for feasibility verification.
The subproblem is a series of PV active power and reactive power optimization problems of the prediction and control horizon considering the constraints of PV power factor and capacity. As the PV capacity limits the effectiveness of reactive power optimization, it is considered to reduce PV active power outputs in critical situations. Compared with the reactive power adjustment costs, the unit cost of PV active power regulation is significantly higher. Based on the MPC algorithm, a voltage prediction model is established from the sensitivity matrix of nodal power and voltage to ensure that the voltage is always in a safe range. Taking the current measured value as the initial state, the time window is rolled to complete the optimal solution in a long period, achieving a dynamic transition to the ideal result and suppressing large fluctuations of the PV reactive power. The first short-period optimal set of PV reactive power instructions is issued and operated to minimize the sum of deviations between the future nodal voltages and the reference voltage.
Based on the iterative solution step of the Benders decomposition method, the global optimal solution of the DRM optimization model can be found [38]. After the master problem is solved, the OLTC tap position result is sent to the corresponding subproblems in the current time window, based on the duality theory. The subproblems implements rolling solution, which feedback the modified Benders cut constraint to the master problem, according to the feasibility verification. Then the master problem is re-solved by adding the new constraint from subproblems. In this way, the master problem and subproblem are iterated repeatedly until the convergence condition is satisfied and the optimal solution is obtained.

4. Reactive Power and Voltage Optimization Model

This section first establishes a common model of the ADN reactive power and voltage optimization problem, which is later decomposed in to the master problem and rolling subproblems according to the DRM method. At last, the solution process is introduced.

4.1. Objective Function

The objective function of the proposed DRM method is composed of three parts, which are the adjustment cost, the voltage deviation and the active power loss. Both OLTC and PV adjustment costs are considered in the objective function:
C t a p = c b Δ N t a p 2
C P V = c P , P V P P V P ˜ P V 2 2 + c Q , P V Δ Q P V 2 2
where P P V and P ˜ P V are the actual output and predicted output of PV active power. Δ N t a p and Δ Q P V are the variations of OLTC tap and PVs reactive power respectively. c b and c P V are the unit cost coefficients of OLTC and PVs. Please note that the reactive power adjustment cost of PV is calculated by the PV reactive power output changes, while the PV active adjustment cost is the power loss caused by the PV active power curtailment when necessary.
The nodal voltage deviation is one of the most significant reliability factors and service quality indices [39,40], because it is harmful to the operation life and efficiency of electrical devices when the voltage deviation beyond a reasonable range. The fluctuations of load and PV generation increases the risk of node voltage exceeding the limit. However, considering the node voltage limit as a constraint may lead to the voltages rising to their maximum limits after optimization. Choosing the voltage deviations from the ideal voltage reference V r e f as an objective function can remove this problem, as shown in Equations (14). If the voltage deviation of the optimization result is large, the voltage is more likely to exceed the limit:
V s u m , d e = V V r e f 2 2
Furthermore, distributed PV generation will affect the power flow direction and change the node voltages. This paper considers minimizing active power losses to improve the economy of ADN operation. The network active power loss can be obtained by:
P l o s s t = i , j N G i j V i 2 + V i 2 2 V i V j cos θ i j Δ t
In general, the objective function voltage optimization is established based on the chronological scenario:
J * = t l = 1 N T l c b Δ N t a p t l 2 + t s = 1 N T s c Q , P V Δ Q P V t l , t s 2 2 + c P , P V P P V t l , t s P ˜ P V t l , t s 2 2 + V t l , t s V r e f 2 2 + P l o s s t l , t s
where tl and ts mean the long-period time window and rolling short-period time window, while N T l is the total number of long periods and N T s is the total short periods in a long period.

4.2. Constraints

Some following constraints are required to set in the ADN reactive power and voltage optimization model. First, based on the sensitivity model for voltage calculation in Section 1, voltage security of each node should be guaranteed. Besides, the PV reactive power output cannot exceed its configured maximum range. Similarly, the OLTC tap is limited by its maximum tap position, as well as the restriction on tap changes of one action.

4.2.1. Voltage Constraints

The nodal voltages of each short period V t l , t s can be estimated by the voltages in the previous period and the predicted variations in the current period with respect to OLTC tap changes Δ V t a p t l , t s and PV and load power Δ V P V , l o a d t l , t s , as shown in Equations (17) and (18). Then, the nodal voltages should be limited within a safe range:
Δ V t a p t l , t s = S t a p Δ V d Δ N t a p t l , t s
Δ V P V , l o a d t l , t s = S P Δ P P V t l , t s + Δ P l o a d t l , t s + S Q Δ Q P V t l , t s + Δ Q l o a d t l , t s
V t l , t s = V t l , t s 1 + Δ V t a p t l , t s + Δ V P V , l o a d t l , t s
V min V t l , t s V max
where V min and V max are the lower limit and upper limit of the voltage amplitude. Δ N t a p t l , t s is the OLTC tap position change. Δ P P V t l , t s , Δ Q P V t l , t s , Δ P l o a d t l , t s and Δ Q l o a d t l , t s are the power variations of PV and load, which can be obtained as follows:
Δ N t a p t l , t s = N t a p t l , t s N t a p t l 1 , t s i f t s = 1 N t a p t l , t s = N t a p t l , 1 i f t s 1
Δ P P V t l , t s = P P V t l , t s P P V t l , t s 1 Δ Q P V t l , t s = Q P V t l , t s Q P V t l , t s 1
Δ P l o a d t l , t s = P l o a d t l , t s P l o a d t l , t s 1 Δ Q l o a d t l , t s = Q l o a d t l , t s Q l o a d t l , t s 1
It is worth noting that the OLTC tap is only changed once in a long period, which is assumed to occur in the first short period, as in Equation (21). The PV active outputs and load conditions in Equations (22) and (23) are based on the current forecast.

4.2.2. PV Power Constraints

Due to the limitations of PV operation mode and capacity in Equations (10) and (11), the reactive adjustment constraint is:
P P V min P P V t l , t s P P V max Q P V min Q P V t l , t s Q P V max
where Q P V max and Q P V min are the reactive power up-limit matrix and down-limit matrix. P P V min and P P V max are the active power up-limit matrix and down-limit matrix.

4.2.3. OLTC Tap Constraints

The OLTC tap range as a key parameter should be considered, as well as the amplitude of one tap change:
N t a p min N t a p t l N t a p max
Δ N t a p t l 0 , ± 1
where N t a p max and N t a p min are the lower limit and upper limit of the OLTC tap position.

4.3. Model Decomposition

In this section, the Benders decomposition method is used to decompose the original optimization problem into a master problem and several MPC-based rolling subproblems to achieve the decoupling of different action periods, as well as the decoupling of discrete and continuous control variables.
(1)
MPC-Based subproblem
The rolling subproblems are based on the MPC algorithm to solve the optimization problem of short-period continuous control variables, PV active and reactive power outputs. We set a predictive and control horizon consisting of H s short periods. According to load and PV forecasts in the prediction and control horizon, the optimal operation sequence of PV power outputs can be obtained. The specific subproblem model solved by the k-th iteration is:
J s u b t l , k Δ N t a p t l , k = m i n t s = 1 N T s ( h s = 1 H s c Q , P V Δ Q P V t l , t s 2 2 + c P , P V P P V t l , t s P ˜ P V t l , t s 2 2 + V t l , t s , h s V r e f 2 2 + P l o s s t l , t s ) S . T . Δ N t a p t l = Δ N t a p t l , k : λ k E q u a t i o n s 17 20 , 22 , 23 E q u a t i o n s 10 , 11 , 24
where λ k is a dual variable generated by constraint ( Δ N t a p t l = Δ N t a p t l , k ), which means that the result Δ N t a p t l of the k-th master problem is used as a parameter when the subproblem is solved.
(2)
Master Problem
The long-period master problem considers the adjustment cost of OLTC and retains the discreteness of the control variable, the OLTC tap. The objective function is to minimize the adjustment cost caused by OLTC tap changes. The master problem is expressed as:
J m t l , k + 1 J s u b t l , k , Δ N t a p t l , k , λ k = min c b Δ N t a p t l 2 + α S . T . J s u b t l , k + λ k Δ N t a p t l Δ N t a p t l , k α α α d o w n E q u a t i o n s ( 19 ) , ( 21 )
The first line of the constraints is the form of Benders cut, where α is an auxiliary variable. α d o w n is a sufficiently small positive number used as a substitute for the objective function value of the subproblem when the master problem is solved for the first time, i.e., k = 1. Afterwards, the master problem is re-solved according to the new constraints in each iteration.

4.4. Solution Algorithm

In this paper, the global optimal solution is obtained through multiple iterations between the master problem and subproblems. The detailed solution process based on Benders decomposition algorithm is shown in Figure 3.
After collecting the parameters and forecasts in the optimization model, initialization is required to enter the Benders loop. The master problem in Equation (28) of OLTC tap position is first solved, and the lower bound L B k of the original problem is updated according to Equation (29). Then, the result of the master problem is substituted into subproblems for feasibility test. If the OLTC tap result satisfies the constraints, the subproblems in Equation (27) are continued to be solved, and the upper bound U B k is updated according to Equation (30). If not, the Benders cut constraint is modified to re-solve the master problem. When the difference between the upper bound and the lower bound is less than the threshold ε in Equation (31), the global optimal solution of the original problem has been found, otherwise the loop needs to be entered again:
L B k = c b Δ N t a p t l 2 + α k
U B k = J s u b t l , k + c b Δ N t a p t l 2
U B k L B k U B k ε

5. Results and Discussion

5.1. Simulation Results of Case 1

A modified PG&E 69-bus distribution network in Figure 4 with specific configurations in Table 2 is selected to test the optimization algorithms performance. More detailed topology information can be obtained by referring to [41]. Figure 5 shows the forecasted aggregate power profile of PV and total daily load profile. The example simulation models are implemented in the MATLAB R2018b environment, and the optimization problems are solved by a commercial solver Gurobi8.9.0 [42]. The upper and lower levels of voltage are 1.05 p.u. and 0.95 p.u. of its base value [43]. The minimum power factor of PV is set to 0.95. Moreover, the unit cost coefficient of OLTC tap is $6 [44] and the upper limit for the daily OLTC tap is set to 20 changes [45,46]. The estimated unit cost coefficient of PV active power and reactive power output are 0.8$/MW and 0.1$/MVar [47] for the economy estimation.
To verify the advantages of the proposed DRM strategy, two other strategies were investigated to compare and analyze the performance in different cases, namely: single-period decoupling strategy (S1) and multi-period decoupling and non-rolling strategy (S2). The strategies are illustrated in Figure 1a,b. Table 3 further presents the differences between them in detail.

5.1.1. Optimal Results

This paper takes 11:00 and 24:00 results as examples to represent voltage deviations of the proposed DRM strategy, as shown in Figure 6. Although some nodes near the end of feeders appear overvoltage (over 1.05 p.u.) at 11:00 due to the high PV outputs, the nodal voltages return to the safe state after optimization. For the result of 24:00 without PV generation, it not only shows that each feeder has a voltage drop due to line impedance, but also illustrates the whole ADN operates within the acceptable range. The optimized overall voltage level is closer to the reference value. The voltage deviation of the Node 16 and Node 54 for a day are shown in Figure 7. The proposed DRM strategy can greatly reduce the voltage amplitudes by adjusting OLTC tap changes and PV reactive power outputs. It is worth noting that Node16, the grid-connected point of PV3, severely exceeded the upper limit at around 11:00, and the voltage drop after optimization is more obvious than that of Node 54.
Figure 8, Figure 9 and Figure 10 show the optimal operation results of OLTC tap and PV power outputs after adopting the DRM strategy. When PV generation is excessive, OLTC raises the tap to reduce the substation secondary side voltage, whereas PVs absorb reactive power to regulate voltages of the surrounding nodes. Remarkably, at around 11:00, even though the voltage has exceeded the limit, the OLTC tap is not further changed, while the active power output of PV4 in Figure 10 is reduced at noon to decrease the voltage level, due to the limited reactive power capacity. The results indicate that the adjustment of PV output can be used to address local voltage problems in feeders. However, because of the limited solar radiation in the morning from 1:00 to 6:00, PV generation is nonexistent. The OLTC tap can only be used to stabilize the voltage level. The numerical results of DRM strategy have been presented in Appendix A for this paper.

5.1.2. Performance Comparison

Adjustment Costs

PV reactive power outputs (for example, PV3) and OLTC tap actions of the three mentioned optimization strategies are shown in Figure 11 and Figure 12, respectively. Since the master problems of OLTC in the three strategies are set within the same long period, the results and adjustment costs of OLTC tap is the same. Table 4 shows that the adjustment costs of S1 and S2 are much higher than that of DRM strategy, which is because of the differences in PV reactive power outputs. The long-period PV forecasts applied by S1 usually present uncertainties, which leads to large changes in PV reactive power outputs of neighbored long periods, and hence the adjustment costs increase significantly. Compared with DRM, S2 lacks the rolling optimization setting of MPC algorithm based on the prediction and control horizon. The optimization function of the rolling subproblems is to gradually reduce the difference between the current state and the ideal state in the long prediction and control horizon. A series of short-period PV reactive output instructions are obtained, although only the first short-period instructions are executed. As shown in Figure 11, the PV reactive power results of several independent short-period subproblems in S2 fluctuate obviously, raising the adjustment costs of PVs, while the based-MPC DRM strategy indicates a significant reduction in cost.

Voltage Deviations

In order to verify the superiority of the proposed DRM strategy in terms of voltage deviations, the voltage results at 11:00 and 24:00 from different strategies is shown in Figure 13. Three strategies can effectively reduce the voltage level, the voltage amplitude of S1 is higher than that of S2 and DRM at 11:00. The voltage results of some nodes by S2 and DRM are similar, while the voltages of the feeder from Node 13 to Node 27 by S2 are higher than the DRM strategy. The difference between the three strategies does not exist at night, because PVs do not generate electricity.
Figure 14 compares the daily voltages of Node 16 and Node 54 optimized by S1, S2 and DRM strategy from the perspective of time changes. The differences between S1 and DRM strategy are still existed, while the voltage results of S2 and DRM are similar. More accurate data are given in Table 5, which is also shows this phenomenon. For the all-day voltage deviations, compared with the original profile, the three optimization methods all achieve a voltage deviation reduction of more than 60%. The S1 voltage deviation is 0.88 p.u. more than that of the DRM due to the errors in the long-period PV forecasts, while the voltage deviation difference between S2 and DRM is 0.23 p.u.

Active Power Losses

Table 6 shows the simulation results of the active power loss from the three strategies. Before PVs are connected to the distribution, the active power loss is as high as 4.8305 MWh a day. The distributed PV power generation can greatly reduce the network losses, only 2.0546 MWh. Due to the PV power transmission, the network loss results of the three optimization strategies are slightly higher than the original data. However, since some original nodal voltages exceed the upper limit, the voltage optimization strategy has to be implemented. Among the three strategies, S1 has the largest active power losses, because the forecast of load and PV power generation is uncertain, and the results of long-period strategy cannot achieve the dynamic optimization, resulting in a large amount of power loss. The active power loss of S2 considering different optimization periods is a little more than that of DRM strategy. The reason is that S2 is the static optimization method of each short period based on forecasts with errors, while the DRM strategy using dynamic real-time forecasts achieves the gradual approximation of the ideal value by rolling the time window.

Computation Time

The computation time of the three strategies is shown in Table 7. The present case only focuses on the differences in computation time caused by the optimization strategy, including modeling and solution, without considering the time spent on PV and load forecasting. Based on the 69-node distribution system configuration of Case 1 and the long period setting of S1, the 1-h daytime optimization result including OLTC tap and PV power requires 21.46 s, and the result of one day requires 5.33 min. Both S2 and the DRM strategy set the short optimization period about 15 min for PV output, and hence there are 96 changes in PV output per day. S2 takes 9.07 s to solve a short-period subproblem of PV power during the day. Since the linear sensitivity between node voltage and power is used instead of the nonlinear power flow constraints, although the optimization problems and variables have been increased by almost 4 times compared with S1, the computation time has only increased by around 2 times. The 1-h optimization requires 48.94 s, and the total computation time of one day is 11.64 min. The DRM strategy applies MPC algorithm to improve the subproblems of S2, and a short-period PV results by rolling optimization only increases by 1.36 s. The optimal results of 96 times a day takes 2.75 min longer than S2 in total. Comparing the computation time of a 1-h optimization, it can be seen that the computation time mainly depends on the number of iterations between master problems and subproblems, while the additional time required for rolling optimization of continuous PV outputs is not much. When the scale of the distribution network expands, the significant increase in nodes and discrete devices will result in more computation time.

5.2. Simulation Results of Case 2

Compared with Case 1, Case 2 is a system configuration with uneven distribution of PVs, as shown in Figure 15. The feeder from Node 28 to Node 35 is completely passive and is not connected to any PV. The parameters of each PV source and OLTC remain unchanged, in Table 2. The forecasted aggregate power profile of PV and total daily load profile are shown in Figure 16.
The voltage results of all nodes at 11:00 and 24:00 are shown in Figure 17. The 11:00 voltage is significantly different from that of Case 1, and a slight overvoltage occurs in the system. In the original profile, the feeder voltage without PV connection from Node 28 to Node 35 gradually decrease due to the power loss caused by the line impedance. Compared with Case 1, the sum of original voltage deviations is smaller, and it is further reduced after optimization by DRM strategy.
Figure 18 reflects the differences in voltage changes at Node 27 and Node 35, the ends of two feeders. The voltage of Node 27 varies considerably, which is affected by the adjustment of the PV power on the feeder between Node 4 and Node 27. The Node 35 is only affected by load fluctuations and OLTC tap adjustments, and hence the changes are small.
The OLTC tap position and PV outputs of the DRM optimization strategy are shown in Figure 19, Figure 20 and Figure 21. Compared with Figure 8 of Case 1, the number of OLTC tap changes is significantly decreased, because the original node voltages of Case 2 do not excessively exceed the upper limit. The frequent OLTC tap changes are avoided due to the impact of OLTC on the substation secondary bus voltage and its higher adjustment cost. We can see that DRM strategy mainly optimizes the voltage by adjusting PV outputs. In Figure 20, PV1 and PV4 absorb a large amount of reactive power to decrease the voltage deviation. The reactive power absorption of PV9 at Node 50 is less because the end node voltage of the feeder is relatively low. From Figure 21, the PV active power outputs are still at the predicted level. It can be explained that PVs provide sufficient reactive power capacity at this time, avoiding the PV active power curtailment.

6. Conclusions

In this paper, a novel DRM method is applied to the power and voltage optimization of ADN cooperating distributed PV and OLTC. According to the discrete and continuous characteristics of control variables, the proposed method utilizes a Benders decomposition algorithm to design a long-period master problem of OLTC tap and short-period subproblems of PV power outputs, achieving the decoupling of different devices actions. The different optimization time step considering the operation principles of devices not only avoids the frequent OLTC tap changes, but also further reduces voltage deviations caused by the stochastic intermittence and fluctuations of PV. The independent subproblems are then further modified into a series of rolling optimization models within long predictive and control horizons based on MPC algorithm considering the uncertainty of PV power generation and load. The rolling optimization of PV power achieves the gradual evolution of the ideal voltage in a long period without causing violent fluctuations in PV reactive power outputs, which greatly decreases adjustment costs. In the case study of the modified PG&E 69-node distribution network, the other two methods are compared to demonstrate the effectiveness and economy of DRM method. The future work would move forward to the implementation of the proposed optimization method considering the high uncertainty of renewable energy generation in large scale ADN. Besides, the impact of energy storage system on the dynamic voltage optimization is introduced. These issues would be established as subsequent research objectives.

Author Contributions

Conceptualization, X.G. and L.S.; methodology, L.S. and C.Z.; software, L.S.; validation, L.S., X.G. and P.L.; investigation, X.G.; resources, X.G.; data curation, C.Z.; writing—original draft preparation, L.S.; writing—review and editing, L.S.; supervision, X.D.; project administration, X.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Technology Program of State Grid Corporation of Zhejiang Province under Grand 5211DS17001Z.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

α auxiliary variable
α d o w n initial lower bound of α
λ dual variable
θ , θ i , i N nodal voltage phase angle
φ max PV maximum power factor angle
B , G set of conductance and susceptance
c b unit cost coefficients of OLTC
c P , P V , c Q , P V unit cost coefficients of PV active power and reactive power
H , N , K , L Elements in the Jacobian matrix
J Jacobian matrix from Newton Raphson power flow equation
J s u b , J m objective function of subproblems and the master problem
L B k , U B k lower and upper bound of the original problem
k number of iterations
S P , S Q power-voltage sensitivity matrix
N t a p , Δ N t a p OLTC tap position and OLTC tap changes
N t a p max , N t a p min lower and upper bound of OLTC tap position
S m , P V PV-m capacity
t l long-period time window
t s rolling short-period time window
N T l total number of long periods
N T s total number of short periods in a long period
V , V i , i N nodal voltage amplitude
V t a p OLTC set-point
Δ V d voltage variation on the low-voltage side after a OLTC tap change
V min , V max lower and upper bound of nodal voltage amplitude
P , Q set of nodal power and reactive power
P l o s s active power loss
P l o a d , Q l o a d set of load power
P ˜ P V set of PV predicted active power outputs
P P V , Δ P P V set of PV actual active power and active power changes
Q P V , Δ Q P V set of PV reactive power and reactive power changes
Q P V min , Q P V max lower and upper bound PVs reactive power

Abbreviations

ADNActive distribution network
DGDistributed generation
DRMDecoupling rolling multi-period
MINPMixed integer nonlinear programming
MPCModel predictive control
MPPTMaximum power point tracking
OPFOptimal power flow
OLTCOn-load tap changer
PCCPoint of common coupling
PVPhotovoltaic

Appendix A

Table A1. Voltage results obtained by DRM strategy of Case 1 (p.u.).
Table A1. Voltage results obtained by DRM strategy of Case 1 (p.u.).
Node0:001:002:003:004:005:006:007:00
11.0500001.0500001.0500001.0500001.0500001.0500001.0500001.050000
21.0242971.0243491.0243311.0243161.0242881.0242581.0242931.036911
31.0242041.0243081.0242711.0242421.0241851.0241261.0241961.036784
41.0241481.0242621.0242211.0241901.0241271.0240581.0241251.036702
51.0236091.0238821.0237831.0237091.0235581.0233741.0234721.035926
61.0186921.0210931.0202271.0195721.0182491.0167981.0182171.029404
71.0138011.0183101.0166831.0154531.0129691.0103101.0132071.023211
81.0126511.0176571.0158511.0144861.0117281.0087891.0120571.021791
91.0120831.0173391.0154431.0140091.0111141.0080261.0114481.021032
101.0067021.0145761.0117351.0095881.0052491.0007261.0061851.014380
111.0056411.0140191.0109971.0087121.0040950.9992671.0050261.012910
121.0027281.0125521.0090091.0063301.0009150.9954191.0027331.010033
130.9997311.0111421.0070281.0039160.9976230.9912920.9999681.006460
140.9968401.0097741.0051111.0015840.9944500.9873610.9974931.003277
150.9940611.0084501.0032640.9993390.9914010.9836350.9953191.000496
160.9935441.0082041.0029200.9989230.9908340.9829430.9949160.999981
170.9927021.0078051.0023620.9982430.9899090.9817110.9938170.998542
180.9926921.0078001.0023550.9982350.9898980.9816970.9938050.998526
190.9921041.0075271.0019680.9977620.9892520.9808570.9931380.997652
200.9917261.0073511.0017200.9974590.9888360.9803170.9927090.997090
210.9911141.0070671.0013180.9969680.9881650.9794450.9920170.996183
220.9910981.0070611.0013080.9969550.9881470.9794220.9920050.996166
230.9909711.0070071.0012280.9968550.9880060.9792680.9919940.996153
240.9906941.0068911.0010540.9966370.9876990.9789320.9919720.996124
250.9902871.0067321.0008060.9963210.9872450.9785130.9923100.996566
260.9901181.0066661.0007030.9961900.9870580.9783400.9924480.996748
270.9900711.0066481.0006740.9961540.9870050.9783250.9926300.996990
281.0241831.0242941.0242541.0242241.0241631.0240971.0241641.036745
291.0239241.0241221.0240511.0239971.0238881.0237441.0237761.036263
301.0231081.0238081.0235541.0233631.0229811.0224831.0226131.034712
311.0229641.0237521.0234661.0232521.0228211.0222601.0224071.034438
321.0222451.0234751.0230281.0226931.0220201.0211461.0213811.033068
331.0207841.0228911.0221261.0215521.0204001.0189581.0195481.030635
341.0189141.0221641.0209841.0200981.0183201.0163681.0182081.028856
351.0180581.0218181.0204541.0194281.0173711.0149051.0163291.026352
361.0241401.0242561.0242141.0241821.0241191.0240491.0241151.036691
371.0239481.0240971.0240431.0240021.0239211.0238231.0238861.036430
381.0237561.0239681.0238911.0238331.0237181.0235771.0236651.036159
391.0237561.0239681.0238911.0238331.0237181.0235821.0236931.036194
401.0126111.0176321.0158201.0144511.0116841.0087621.0121251.021888
411.0125621.0176101.0157881.0144121.0116311.0088031.0125581.022472
421.0114891.0170351.0150341.0135211.0104661.0072091.0108281.020245
431.0108341.0166941.0145801.0129811.0097531.0063231.0101911.019443
441.0099681.0162521.0139851.0122711.0088091.0051531.0093841.018424
451.0091571.0158451.0134321.0116081.0079231.0039631.0082361.016940
461.0050101.0139821.0107451.0082981.0033540.9976981.0022741.009103
471.0029661.0130651.0094211.0066671.0011030.9946110.9993361.005240
481.0021761.0127121.0089101.0060371.0002320.9934150.9981981.003744
491.0013971.0124021.0084311.0054300.9993660.9922250.9971531.002350
501.0002811.0119041.0077111.0045410.9981370.9905620.9956551.000383
511.0001551.0118431.0076261.0044380.9979990.9903580.9954001.000047
521.0000101.0117711.0075281.0043200.9978410.9901250.9951090.999666
530.9993001.0114191.0070461.0037410.9970650.9889850.9936840.997795
540.9988001.0111491.0066931.0033260.9965220.9882090.9927240.996549
551.0055251.0139731.0109251.0086211.0039670.9990691.0047731.012572
561.0055241.0139721.0109241.0086201.0039650.9990671.0047701.012568
571.0024521.0124001.0088111.0060981.0006150.9952661.0034111.010974
581.0024511.0123991.0088111.0060981.0006140.9952671.0034181.010982
591.0241801.0242911.0242511.0242211.0241591.0240931.0241591.036739
601.0238781.0240811.0240081.0239521.0238411.0236921.0237201.036201
611.0236051.0239201.0238061.0237201.0235471.0233371.0234531.035876
621.0235261.0238731.0237471.0236521.0234631.0232351.0233761.035782
631.0235221.0238711.0237441.0236491.0234581.0232291.0233711.035775
641.0223181.0231641.0228571.0226261.0221641.0216141.0219731.034042
651.0218401.0228761.0225001.0222171.0216511.0209871.0214561.033407
661.0217771.0228381.0224531.0221631.0215831.0209041.0213881.033324
671.0217631.0228301.0224421.0221511.0215681.0208851.0213721.033305
681.0215971.0227291.0223181.0220091.0213911.0206671.0211891.033080
691.0215961.0227291.0223181.0220091.0213901.0206671.0211891.033080
Node8:009:0010:0011:0012:0013:0014:0015:00
11.0500001.0500001.0500001.0500001.0500001.0500001.0500001.050000
21.0368521.0368621.0243151.0243971.0243431.0369191.0369541.036933
31.0366661.0366871.0242401.0244041.0242971.0368021.0368711.036829
41.0365681.0365871.0241551.0243351.0242191.0367151.0367821.036732
51.0355811.0356041.0234441.0238621.0236021.0359221.0360241.035891
61.0265561.0269871.0188421.0225381.0201251.0296911.0310711.030044
71.0179361.0188361.0147421.0217111.0171191.0238891.0266801.024804
81.0159501.0169691.0138401.0215851.0164731.0225741.0257171.023645
91.0148951.0159621.0133081.0214381.0160711.0218481.0251431.022966
101.0053091.0070761.0094661.0215571.0134321.0157801.0210171.017874
111.0032371.0051091.0084361.0212741.0126351.0143651.0199031.016555
120.9988821.0013311.0075511.0225251.0122771.0120551.0190131.015247
130.9935500.9965361.0057961.0230121.0110541.0088891.0171441.012866
140.9887300.9922961.0045161.0238791.0102481.0061831.0158071.011078
150.9844450.9886331.0037271.0251371.0098721.0039551.0150231.009905
160.9836500.9879541.0035841.0253741.0098041.0035431.0148801.009690
170.9816220.9859971.0023951.0247741.0087621.0020661.0135801.008200
180.9815990.9859751.0023821.0247681.0087511.0020491.0135661.008184
190.9803370.9847991.0017641.0245541.0082081.0011921.0128981.007405
200.9795250.9840431.0013671.0244171.0078581.0006421.0124691.006904
210.9782140.9828231.0007291.0241991.0072970.9997541.0117791.006099
220.9781870.9828011.0007251.0242071.0072930.9997401.0117751.006093
230.9781220.9828031.0008771.0244351.0074220.9997901.0119531.006275
240.9779790.9828081.0012081.0249321.0077030.9998981.0123411.006673
250.9783460.9835281.0024801.0264161.0087951.0007071.0137991.008226
260.9784960.9838251.0030041.0270281.0092451.0010401.0144011.008867
270.9787700.9841901.0034531.0274821.0096311.0013871.0149111.009420
281.0366161.0366361.0241991.0243741.0242611.0367591.0368261.036780
291.0360071.0360111.0237041.0240091.0238291.0362371.0362851.036188
301.0338131.0338281.0224041.0234941.0228451.0346291.0348231.034476
311.0334261.0334431.0221741.0234041.0226711.0343461.0345641.034174
321.0314901.0315161.0210271.0229501.0218031.0329271.0332741.032664
331.0279951.0281091.0192321.0225491.0205431.0305061.0312551.030249
341.0251061.0256291.0191441.0243631.0210751.0292151.0311321.029787
351.0217601.0221001.0163371.0223101.0186441.0263441.0279891.026271
361.0365551.0365741.0241451.0243261.0242091.0367051.0367701.036720
371.0362491.0362631.0238781.0241151.0239661.0364361.0364951.036419
381.0359041.0359161.0236251.0239771.0237571.0361771.0362401.036113
391.0359471.0359641.0236771.0240361.0238081.0362281.0363031.036174
401.0160591.0171131.0140431.0218221.0166751.0227261.0259461.023885
411.0167461.0179421.0150531.0229111.0176601.0235401.0270781.025086
421.0137661.0148931.0128161.0214151.0157481.0211141.0246011.022305
431.0126061.0138091.0123721.0214831.0154841.0203911.0241161.021697
441.0111121.0124291.0118661.0216701.0152171.0194941.0235621.020982
451.0090651.0103791.0105481.0209851.0141611.0179411.0220951.019299
460.9980830.9993971.0037221.0177451.0088031.0097541.0144521.010444
470.9926700.9939841.0003581.0161511.0061651.0057191.0106861.006080
480.9905720.9918870.9990561.0155371.0051441.0041571.0092281.004390
490.9885620.9899080.9979601.0151831.0043401.0027411.0079911.002923
500.9857740.9871540.9963591.0145551.0031261.0007311.0061931.000819
510.9853270.9866830.9959871.0142811.0028031.0003491.0057781.000355
520.9848180.9861490.9955641.0139701.0024360.9999141.0053060.999828
530.9823290.9835320.9934911.0124441.0006350.9977861.0029960.997247
540.9806840.9818050.9921151.0114200.9994340.9963741.0014680.995544
551.0027831.0046311.0080591.0210001.0123091.0139781.0194811.016082
561.0027781.0046251.0080551.0209971.0123051.0139731.0194761.016076
570.9999541.0027021.0094051.0246381.0141191.0134691.0210991.017444
580.9999641.0027141.0094201.0246531.0141331.0134811.0211151.017462
591.0366101.0366291.0241921.0243681.0242541.0367531.0368191.036772
601.0359371.0359331.0236241.0239421.0237511.0361761.0362111.036095
611.0354911.0355121.0234161.0239221.0236081.0358811.0359951.035829
621.0353621.0353901.0233561.0239171.0235671.0357961.0359321.035753
631.0353541.0353811.0233511.0239141.0235621.0357901.0359261.035746
641.0330291.0331001.0219591.0233311.0224741.0340861.0344361.034005
651.0321771.0322751.0214831.0231701.0221131.0334801.0339321.033410
661.0320641.0321661.0214211.0231491.0220651.0334001.0338661.033331
671.0320391.0321411.0214061.0231441.0220541.0333821.0338501.033313
681.0317391.0318501.0212331.0230801.0219211.0331661.0336681.033098
691.0317381.0318491.0212331.0230801.0219211.0331661.0336681.033098
Node16:0017:0018:0019:0020:0021:0022:0023:00
11.0500001.0500001.0500001.0500001.0500001.0500001.0500001.050000
21.0369201.0369031.0368641.0368601.0368401.0368691.0242591.024297
31.0368041.0367691.0366901.0366831.0366441.0367021.0241281.024204
41.0367121.0366781.0366041.0366041.0365661.0366301.0240631.024148
51.0358871.0358291.0357201.0357681.0357061.0358611.0234061.023609
61.0296171.0289111.0273271.0273131.0264931.0278601.0169091.018692
71.0238401.0224311.0191981.0189941.0173441.0199111.0104511.013801
81.0225401.0209571.0173101.0170461.0151901.0180401.0089311.012651
91.0218101.0201501.0163311.0160611.0141171.0171101.0081781.012083
101.0158701.0132801.0071861.0065141.0034331.0079201.0008501.006702
111.0144451.0117081.0052801.0046021.0013481.0061240.9994131.005641
121.0123691.0089811.0008260.9995860.9954891.0010980.9954221.002728
131.0093341.0053540.9956840.9940930.9892370.9957670.9912360.999731
141.0068131.0022150.9909260.9888800.9832240.9906380.9872050.996840
151.0048240.9995820.9865660.9839590.9774610.9857190.9833360.994061
161.0044570.9990940.9857560.9830450.9763900.9848050.9826170.993544
171.0029260.9974820.9839940.9813800.9746380.9833110.9814420.992702
181.0029100.9974640.9839740.9813600.9746170.9832930.9814280.992692
191.0020450.9965180.9828360.9802280.9733840.9822420.9806050.992104
201.0014900.9959090.9821050.9794990.9725910.9815670.9800760.991726
211.0005940.9949280.9809230.9783220.9713080.9804750.9792190.991114
221.0005830.9949120.9808960.9782910.9712720.9804440.9791960.991098
231.0006700.9949490.9807750.9780880.9709940.9802080.9790140.990971
241.0008620.9950290.9805120.9776450.9703870.9796950.9786170.990694
251.0019020.9958140.9804730.9771130.9694670.9789200.9780220.990287
261.0023310.9961380.9804570.9768930.9690880.9786000.9777770.990118
271.0027460.9964900.9805950.9768860.9689810.9785100.9777080.990071
281.0367581.0367241.0366491.0366471.0366091.0366711.0241021.024183
291.0362071.0361691.0361371.0362011.0361791.0362881.0237791.023924
301.0345271.0343961.0342511.0344651.0343741.0347571.0225971.023108
311.0342311.0340821.0339181.0341591.0340561.0344871.0223881.022964
321.0327491.0325171.0322541.0326271.0324641.0331371.0213451.022245
331.0302711.0298101.0291491.0296451.0292741.0304251.0192431.020784
341.0292161.0281721.0261221.0261691.0251441.0269211.0165361.018914
351.0260681.0251131.0235191.0241211.0232781.0253351.0153061.018058
361.0367011.0366671.0365931.0365931.0365561.0366201.0240551.024140
371.0364161.0363851.0363171.0363361.0363011.0363831.0238391.023948
381.0361251.0360961.0360051.0360311.0359871.0361041.0236011.023756
391.0361771.0361451.0360311.0360401.0359871.0361041.0236011.023756
401.0227231.0211061.0173531.0170221.0151191.0179781.0088801.012611
411.0236771.0219241.0177371.0171141.0150181.0178921.0088121.012562
421.0210791.0193191.0152771.0149901.0129371.0160941.0073681.011489
431.0203721.0184911.0141691.0138351.0116461.0149831.0064791.010834
441.0195071.0174531.0127231.0123081.0099221.0135001.0052991.009968
451.0178731.0157571.0109741.0107161.0082911.0121001.0041861.009157
461.0092581.0067671.0016131.0021940.9995171.0046270.9983421.005010
471.0050121.0023370.9969990.9979940.9951931.0009440.9954621.002966
481.0033671.0006210.9952110.9963660.9935170.9995160.9943461.002176
491.0018950.9990430.9934650.9947280.9917870.9980530.9932191.001397
500.9998010.9968180.9910470.9924720.9894200.9960380.9916441.000281
510.9993830.9964100.9906960.9921920.9891630.9958180.9914701.000155
520.9989090.9959480.9902980.9918730.9888700.9955670.9912711.000010
530.9965840.9936810.9883440.9903130.9874360.9943360.9902950.999300
540.9950460.9921800.9870440.9892670.9864700.9935020.9896240.998800
551.0140211.0112951.0049271.0043241.0010941.0059090.9992461.005525
561.0140161.0112901.0049231.0043201.0010911.0059060.9992441.005524
571.0140621.0103811.0013050.9994740.9949611.0006410.9950521.002452
581.0140761.0103931.0013100.9994750.9949601.0006390.9950511.002451
591.0367511.0367171.0366431.0366421.0366051.0366661.0240981.024180
601.0361261.0361041.0360751.0361451.0361251.0362371.0237301.023878
611.0358311.0357721.0356451.0356971.0356331.0358051.0233741.023605
621.0357461.0356761.0355211.0355681.0354901.0356801.0232721.023526
631.0357391.0356691.0355141.0355611.0354831.0356741.0232671.023522
641.0339741.0337951.0334001.0335031.0333071.0337671.0217001.022318
651.0333531.0331241.0326071.0327111.0324571.0330211.0210831.021840
661.0332711.0330351.0325021.0326061.0323451.0329221.0210021.021777
671.0332521.0330151.0324791.0325831.0323201.0329011.0209841.021763
681.0330301.0327771.0322021.0323081.0320271.0326431.0207701.021597
691.0330301.0327761.0322011.0323071.0320261.0326421.0207691.021596
Table A2. Voltage results obtained by S1 of Case 1 (p.u.).
Table A2. Voltage results obtained by S1 of Case 1 (p.u.).
Node0:001:002:003:004:005:006:007:00
11.0500001.0500001.0500001.0500001.0500001.0500001.0500001.050000
21.0242971.0243491.0243311.0243161.0242881.0242581.0242931.036912
31.0242041.0243081.0242711.0242421.0241851.0241261.0241971.036786
41.0241481.0242621.0242211.0241901.0241271.0240581.0241261.036706
51.0236091.0238821.0237831.0237091.0235581.0233771.0234831.035945
61.0186921.0210931.0202271.0195721.0182491.0168221.0182891.029551
71.0138011.0183101.0166831.0154531.0129691.0103551.0133431.023488
81.0126511.0176571.0158511.0144861.0117281.0088401.0122081.022100
91.0120831.0173391.0154431.0140091.0111141.0080791.0116071.021357
101.0067021.0145761.0117351.0095881.0052491.0008001.0064161.014868
111.0056411.0140191.0109971.0087121.0040950.9993461.0052751.013435
121.0027281.0125521.0090091.0063301.0009150.9955141.0030441.010688
130.9997311.0111421.0070281.0039160.9976230.9913821.0002671.007153
140.9968401.0097741.0051111.0015840.9944500.9874470.9977801.004008
150.9940611.0084501.0032640.9993390.9914010.9837170.9955941.001267
160.9935441.0082041.0029200.9989230.9908340.9830240.9951881.000758
170.9927021.0078051.0023620.9982430.9899090.9817900.9940860.999330
180.9926921.0078001.0023550.9982350.9898980.9817750.9940740.999314
190.9921041.0075271.0019680.9977620.9892520.9809330.9934040.998449
200.9917261.0073511.0017200.9974590.9888360.9803920.9929730.997893
210.9911141.0070671.0013180.9969680.9881650.9795170.9922790.996996
220.9910981.0070611.0013080.9969550.9881470.9794950.9922660.996980
230.9909711.0070071.0012280.9968550.9880060.9793390.9922550.996970
240.9906941.0068911.0010540.9966370.9876990.9790000.9922290.996950
250.9902871.0067321.0008060.9963210.9872450.9785760.9925600.997412
260.9901181.0066661.0007030.9961900.9870580.9784010.9926960.997602
270.9900711.0066481.0006740.9961540.9870050.9783850.9928760.997848
281.0241831.0242941.0242541.0242241.0241631.0240971.0241631.036745
291.0239241.0241221.0240511.0239971.0238881.0237391.0237501.036234
301.0231081.0238081.0235541.0233631.0229811.0224731.0225661.034660
311.0229641.0237521.0234661.0232521.0228211.0222501.0223571.034382
321.0222451.0234751.0230281.0226931.0220201.0211331.0213131.032994
331.0207841.0228911.0221261.0215521.0204001.0189371.0194361.030514
341.0189141.0221641.0209841.0200981.0183201.0163311.0180081.028642
351.0180581.0218181.0204541.0194281.0173711.0148681.0161291.026138
361.0241401.0242561.0242141.0241821.0241191.0240491.0241171.036695
371.0239481.0240971.0240431.0240021.0239211.0238241.0238841.036434
381.0237561.0239681.0238911.0238331.0237181.0235791.0236551.036159
391.0237561.0239681.0238911.0238331.0237181.0235851.0236801.036193
401.0126111.0176321.0158201.0144511.0116841.0088131.0122761.022198
411.0125621.0176101.0157881.0144121.0116311.0088541.0127071.022784
421.0114891.0170351.0150341.0135211.0104661.0072651.0109931.020583
431.0108341.0166941.0145801.0129811.0097531.0063821.0103641.019796
441.0099681.0162521.0139851.0122711.0088091.0052161.0095661.018797
451.0091571.0158451.0134321.0116081.0079231.0040291.0084301.017331
461.0050101.0139821.0107451.0082981.0033540.9977781.0025101.009571
471.0029661.0130651.0094211.0066671.0011030.9946980.9995931.005746
481.0021761.0127121.0089101.0060371.0002320.9935050.9984641.004264
491.0013971.0124021.0084311.0054300.9993660.9923180.9974281.002887
501.0002811.0119041.0077111.0045410.9981370.9906600.9959501.000953
511.0001551.0118431.0076261.0044380.9979990.9904570.9956951.000617
521.0000101.0117711.0075281.0043200.9978410.9902240.9954041.000236
530.9993001.0114191.0070461.0037410.9970650.9890840.9939790.998366
540.9988001.0111491.0066931.0033260.9965220.9883080.9930200.997121
551.0055251.0139731.0109251.0086211.0039670.9991481.0050221.013097
561.0055241.0139721.0109241.0086201.0039650.9991451.0050191.013094
571.0024521.0124001.0088111.0060981.0006150.9953801.0037941.011736
581.0024511.0123991.0088111.0060981.0006140.9953811.0038011.011745
591.0241801.0242911.0242511.0242211.0241591.0240931.0241601.036741
601.0238781.0240811.0240081.0239521.0238411.0236931.0237161.036202
611.0236051.0239201.0238061.0237201.0235471.0233391.0234461.035877
621.0235261.0238731.0237471.0236521.0234631.0232361.0233681.035783
631.0235221.0238711.0237441.0236491.0234581.0232311.0233631.035777
641.0223181.0231641.0228571.0226261.0221641.0216181.0219561.034045
651.0218401.0228761.0225001.0222171.0216511.0209921.0214341.033411
661.0217771.0228381.0224531.0221631.0215831.0209091.0213651.033327
671.0217631.0228301.0224421.0221511.0215681.0208901.0213491.033308
681.0215971.0227291.0223181.0220091.0213911.0206731.0211651.033084
691.0215961.0227291.0223181.0220091.0213901.0206721.0211651.033084
Node8:009:0010:0011:0012:0013:0014:0015:00
11.0500001.0500001.0500001.0500001.0500001.0500001.0500001.050000
21.0368531.0368571.0243211.0244001.0243491.0369191.0369541.036930
31.0366701.0366781.0242521.0244091.0243091.0368011.0368701.036822
41.0365741.0365791.0241701.0243431.0242341.0367161.0367841.036729
51.0356071.0356031.0234951.0239001.0236541.0359371.0360501.035910
61.0267581.0268851.0192651.0228511.0205971.0297921.0312331.030090
71.0183201.0186261.0155511.0223071.0180261.0240771.0269801.024871
81.0163781.0167331.0147421.0222491.0174851.0227831.0260491.023716
91.0153461.0157161.0142561.0221381.0171371.0220701.0254931.023042
101.0059901.0066431.0109071.0227121.0151761.0161301.0215221.017915
111.0039701.0046521.0099791.0225361.0145271.0147511.0204501.016604
120.9997971.0007661.0094581.0241761.0147121.0125641.0197001.015300
130.9945370.9957041.0079211.0249991.0140031.0094101.0178031.012717
140.9897900.9911911.0068621.0262091.0137181.0067171.0164381.010723
150.9855800.9872521.0062991.0278171.0138741.0045031.0156261.009339
160.9847980.9865211.0061971.0281201.0139051.0040941.0154781.009085
170.9827910.9844971.0050621.0276481.0130361.0026301.0141771.007544
180.9827680.9844741.0050501.0276441.0130261.0026141.0141631.007528
190.9815230.9832401.0044781.0275421.0126321.0017691.0134941.006704
200.9807220.9824471.0041101.0274771.0123791.0012271.0130651.006175
210.9794290.9811661.0035201.0273751.0119731.0003511.0123741.005324
220.9794030.9811411.0035181.0273881.0119761.0003381.0123701.005316
230.9793450.9811161.0036911.0276691.0121761.0003941.0125471.005477
240.9792190.9810621.0040671.0282831.0126111.0005141.0129341.005829
250.9796210.9816551.0054371.0300171.0140361.0013501.0143901.007284
260.9797860.9818991.0060021.0307331.0146231.0016941.0149901.007884
270.9800680.9822351.0064731.0312441.0150841.0020461.0154991.008415
281.0366181.0366231.0242071.0243751.0242691.0367561.0368221.036769
291.0359681.0359471.0236591.0239551.0237831.0361901.0362221.036109
301.0337481.0336891.0223411.0233931.0227731.0345411.0347011.034311
311.0333571.0332911.0221091.0232951.0225941.0342501.0344331.033994
321.0313981.0312981.0209461.0228001.0217031.0327951.0330911.032408
331.0278471.0277331.0191121.0223001.0203871.0302841.0309481.029812
341.0248481.0249361.0189521.0239171.0208101.0288161.0305801.028987
351.0215011.0214041.0161451.0218641.0183791.0259441.0274361.025469
361.0365611.0365661.0241601.0243341.0242251.0367051.0367731.036717
371.0362531.0362521.0238971.0241201.0239841.0364311.0364951.036420
381.0359011.0358961.0236591.0239761.0237831.0361501.0362331.036124
391.0359411.0359421.0237141.0240331.0238361.0361951.0362941.036188
401.0164891.0168701.0149521.0224851.0176931.0229321.0262761.023951
411.0171811.0176721.0159861.0235741.0186951.0237391.0274001.025130
421.0142361.0146281.0138041.0221321.0168461.0213401.0249621.022374
431.0130971.0135241.0134061.0222201.0166201.0206201.0244891.021758
441.0116341.0121141.0129641.0224331.0164051.0197291.0239511.021032
451.0096131.0100561.0116901.0217761.0153861.0181901.0225071.019354
460.9987460.9989951.0050751.0186431.0102031.0100471.0149401.010472
470.9933900.9935431.0018151.0171011.0076501.0060341.0112111.006095
480.9913140.9914301.0005531.0165061.0066631.0044791.0097671.004399
490.9893300.9894310.9995061.0161761.0058991.0030721.0085451.002923
500.9865890.9866640.9979821.0155961.0047511.0010871.0067881.000828
510.9861420.9861930.9976111.0153231.0044291.0007051.0063731.000363
520.9856340.9856580.9971881.0150121.0040621.0002701.0059020.999837
530.9831470.9830400.9951181.0134871.0022630.9981421.0035930.997255
540.9815040.9813120.9937451.0124651.0010640.9967311.0020660.995552
551.0035161.0041741.0096021.0222631.0142021.0143641.0200281.016132
561.0035111.0041681.0095981.0222601.0141981.0143591.0200231.016126
571.0010021.0022151.0115301.0264491.0167481.0140931.0219471.017639
581.0010131.0022281.0115461.0264661.0167631.0141061.0219641.017657
591.0366131.0366191.0242041.0243731.0242661.0367521.0368181.036765
601.0359371.0359191.0236421.0239441.0237671.0361651.0362071.036094
611.0354921.0354851.0234481.0239221.0236331.0358621.0359851.035823
621.0353631.0353591.0233921.0239151.0235951.0357741.0359221.035744
631.0353551.0353511.0233861.0239121.0235901.0357671.0359161.035737
641.0330301.0330251.0220421.0233221.0225351.0340331.0344081.033979
651.0321791.0321801.0215871.0231581.0221871.0334131.0338971.033377
661.0320661.0320691.0215271.0231371.0221411.0333311.0338301.033297
671.0320411.0320441.0215131.0231311.0221311.0333131.0338141.033279
681.0317411.0317451.0213471.0230651.0220021.0330921.0336281.033061
691.0317401.0317451.0213471.0230661.0220031.0330921.0336291.033062
Node16:0017:0018:0019:0020:0021:0022:0023:00
11.0500001.0500001.0500001.0500001.0500001.0500001.0500001.050000
21.0369221.0369001.0368661.0368611.0368401.0368691.0242591.024297
31.0368081.0367641.0366941.0366861.0366441.0367021.0241281.024204
41.0367191.0366751.0366091.0366071.0365661.0366301.0240631.024148
51.0359201.0358371.0357381.0357781.0357061.0358611.0234061.023609
61.0298671.0289131.0274751.0273971.0264931.0278601.0169091.018692
71.0243141.0224251.0194801.0191561.0173441.0199111.0104511.013801
81.0230671.0209481.0176251.0172271.0151901.0180401.0089311.012651
91.0223651.0201431.0166621.0162511.0141171.0171101.0081781.012083
101.0167081.0132471.0076821.0068021.0034331.0079201.0008501.006702
111.0153471.0116831.0058101.0049101.0013481.0061240.9994131.005641
121.0134951.0089701.0014790.9999660.9954891.0010980.9954221.002728
131.0105451.0052340.9963830.9945180.9892370.9957670.9912360.999731
141.0081091.0019840.9916710.9893500.9832240.9906380.9872050.996840
151.0062090.9992380.9873580.9844750.9774610.9857190.9833360.994061
161.0058580.9987300.9865570.9835690.9763900.9848050.9826170.993544
171.0043530.9971000.9848030.9819140.9746380.9833110.9814420.992702
181.0043360.9970820.9847830.9818940.9746170.9832930.9814280.992692
191.0034940.9961200.9836510.9807700.9733840.9822420.9806050.992104
201.0029520.9955020.9829230.9800470.9725910.9815670.9800760.991726
211.0020800.9945050.9817480.9788780.9713080.9804750.9792190.991114
221.0020690.9944880.9817210.9788480.9712720.9804440.9791960.991098
231.0021670.9945180.9816020.9786490.9709940.9802080.9790140.990971
241.0023800.9945820.9813450.9782140.9703870.9796950.9786170.990694
251.0034660.9953330.9813180.9777000.9694670.9789200.9780220.990287
261.0039140.9956430.9813060.9774870.9690880.9786000.9777770.990118
271.0043400.9959870.9814470.9774840.9689810.9785100.9777080.990071
281.0367591.0367151.0366511.0366491.0366091.0366711.0241021.024183
291.0361531.0361181.0361161.0361951.0361791.0362881.0237791.023924
301.0344321.0342871.0342201.0344591.0343741.0347571.0225971.023108
311.0341281.0339641.0338861.0341521.0340561.0344871.0223881.022964
321.0326101.0323481.0322131.0326201.0324641.0331371.0213451.022245
331.0300451.0295201.0290881.0296351.0292741.0304251.0192431.020784
341.0288191.0276381.0260211.0261571.0251441.0269211.0165361.018914
351.0256691.0245781.0234181.0241081.0232781.0253351.0153061.018058
361.0367081.0366631.0365981.0365971.0365561.0366201.0240551.024140
371.0364231.0363771.0363211.0363401.0363011.0363831.0238391.023948
381.0361351.0360741.0360051.0360391.0359871.0361041.0236011.023756
391.0361881.0361191.0360301.0360481.0359871.0361041.0236011.023756
401.0232511.0210921.0176701.0172041.0151191.0179781.0088801.012611
411.0242111.0218911.0180621.0173021.0150181.0178921.0088121.012562
421.0216541.0193041.0156241.0151881.0129371.0160941.0073681.011489
431.0209701.0184671.0145331.0140441.0116461.0149831.0064791.010834
441.0201391.0174181.0131141.0125301.0099221.0135001.0052991.009968
451.0185331.0157251.0113831.0109471.0082911.0121001.0041861.009157
461.0100361.0067191.0021131.0024710.9995171.0046270.9983421.005010
471.0058481.0022810.9975440.9982930.9951931.0009440.9954621.002966
481.0042251.0005610.9957730.9966730.9935170.9995160.9943461.002176
491.0027790.9989770.9940480.9950460.9917870.9980530.9932191.001397
501.0007350.9967590.9916640.9928060.9894200.9960380.9916441.000281
511.0003180.9963520.9913130.9925260.9891630.9958180.9914701.000155
520.9998440.9958890.9909150.9922080.9888700.9955670.9912711.000010
530.9975210.9936220.9889620.9906480.9874360.9943360.9902950.999300
540.9959850.9921210.9876640.9896020.9864700.9935020.9896240.998800
551.0149231.0112701.0054571.0046321.0010941.0059090.9992461.005525
561.0149181.0112651.0054531.0046291.0010911.0059060.9992441.005524
571.0153561.0104611.0020510.9998960.9949611.0006410.9950521.002452
581.0153711.0104731.0020570.9998970.9949601.0006390.9950511.002451
591.0367551.0367121.0366471.0366451.0366051.0366661.0240981.024180
601.0361321.0360921.0360771.0361491.0361251.0362371.0237301.023878
611.0358391.0357481.0356491.0357041.0356331.0358051.0233741.023605
621.0357551.0356491.0355261.0355761.0354901.0356801.0232721.023526
631.0357481.0356421.0355191.0355691.0354831.0356741.0232671.023522
641.0339931.0337281.0334121.0335221.0333071.0337671.0217001.022318
651.0333771.0330401.0326221.0327341.0324571.0330211.0210831.021840
661.0332951.0329491.0325181.0326301.0323451.0329221.0210021.021777
671.0332771.0329291.0324951.0326071.0323201.0329011.0209841.021763
681.0330561.0326841.0322181.0323341.0320271.0326431.0207701.021597
691.0330561.0326841.0322171.0323331.0320261.0326421.0207691.021596
Table A3. Voltage results obtained by S2 of Case 1 (p.u.).
Table A3. Voltage results obtained by S2 of Case 1 (p.u.).
Node0:001:002:003:004:005:006:007:00
11.0500001.0500001.0500001.0500001.0500001.0500001.0500001.050000
21.0242971.0243491.0243311.0243161.0242881.0242581.0242931.036913
31.0242041.0243081.0242711.0242421.0241851.0241261.0241961.036788
41.0241481.0242621.0242211.0241901.0241271.0240571.0241251.036711
51.0236091.0238821.0237831.0237091.0235581.0233731.0234781.035970
61.0186921.0210931.0202271.0195721.0182491.0167951.0182611.029687
71.0138011.0183101.0166831.0154531.0129691.0103041.0132901.023740
81.0126511.0176571.0158511.0144861.0117281.0087831.0121491.022379
91.0120831.0173391.0154431.0140091.0111141.0080201.0115451.021649
101.0067021.0145761.0117351.0095881.0052491.0007161.0063481.015232
111.0056411.0140191.0109971.0087121.0040950.9992581.0052061.013813
121.0027281.0125521.0090091.0063301.0009150.9954101.0029711.011108
130.9997311.0111421.0070281.0039160.9976230.9912631.0001881.007631
140.9968401.0097741.0051111.0015840.9944500.9873110.9976951.004542
150.9940611.0084501.0032640.9993390.9914010.9835640.9955031.001856
160.9935441.0082041.0029200.9989230.9908340.9828680.9950961.001358
170.9927021.0078051.0023620.9982430.9899090.9816310.9939930.999940
180.9926921.0078001.0023550.9982350.9898980.9816160.9939810.999924
190.9921041.0075271.0019680.9977620.9892520.9807720.9933100.999068
200.9917261.0073511.0017200.9974590.9888360.9802290.9928780.998517
210.9911141.0070671.0013180.9969680.9881650.9793510.9921830.997629
220.9910981.0070611.0013080.9969550.9881470.9793290.9921700.997613
230.9909711.0070071.0012280.9968550.9880060.9791720.9921580.997608
240.9906941.0068911.0010540.9966370.9876990.9788300.9921320.997596
250.9902871.0067321.0008060.9963210.9872450.9784010.9924600.998076
260.9901181.0066661.0007030.9961900.9870580.9782240.9925950.998274
270.9900711.0066481.0006740.9961540.9870050.9782060.9927750.998524
281.0241831.0242941.0242541.0242241.0241631.0240961.0241631.036748
291.0239241.0241221.0240511.0239971.0238881.0237361.0237481.036253
301.0231081.0238081.0235541.0233631.0229811.0224671.0225641.034689
311.0229641.0237521.0234661.0232521.0228211.0222441.0223551.034413
321.0222451.0234751.0230281.0226931.0220201.0211241.0213091.033033
331.0207841.0228911.0221261.0215521.0204001.0189221.0194301.030574
341.0189141.0221641.0209841.0200981.0183201.0163031.0179981.028743
351.0180581.0218181.0204541.0194281.0173711.0148401.0161181.026239
361.0241401.0242561.0242141.0241821.0241191.0240481.0241161.036701
371.0239481.0240971.0240431.0240021.0239211.0238201.0238821.036461
381.0237561.0239681.0238911.0238331.0237181.0235641.0236501.036259
391.0237561.0239681.0238911.0238331.0237181.0235671.0236741.036314
401.0126111.0176321.0158201.0144511.0116841.0087551.0122161.022481
411.0125621.0176101.0157881.0144121.0116311.0087941.0126471.023075
421.0114891.0170351.0150341.0135211.0104661.0072031.0109221.020891
431.0108341.0166941.0145801.0129811.0097531.0063171.0102821.020123
441.0099681.0162521.0139851.0122711.0088091.0051481.0094701.019152
451.0091571.0158451.0134321.0116081.0079231.0039601.0083201.017700
461.0050101.0139821.0107451.0082981.0033540.9977031.0023491.009989
471.0029661.0130651.0094211.0066671.0011030.9946190.9994061.006188
481.0021761.0127121.0089101.0060371.0002320.9934250.9982671.004716
491.0013971.0124021.0084311.0054300.9993660.9922370.9972201.003349
501.0002811.0119041.0077111.0045410.9981370.9905770.9957181.001441
511.0001551.0118431.0076261.0044380.9979990.9903730.9954631.001105
521.0000101.0117711.0075281.0043200.9978410.9901410.9951721.000724
530.9993001.0114191.0070461.0037410.9970650.9890010.9937470.998855
540.9988001.0111491.0066931.0033260.9965220.9882250.9927870.997611
551.0055251.0139731.0109251.0086211.0039670.9990601.0049531.013475
561.0055241.0139721.0109241.0086201.0039650.9990571.0049501.013471
571.0024521.0124001.0088111.0060981.0006150.9952721.0037221.012157
581.0024511.0123991.0088111.0060981.0006140.9952731.0037291.012166
591.0241801.0242911.0242511.0242211.0241591.0240921.0241591.036745
601.0238781.0240811.0240081.0239521.0238411.0236871.0237141.036239
611.0236051.0239201.0238061.0237201.0235471.0233291.0234431.035937
621.0235261.0238731.0237471.0236521.0234631.0232261.0233651.035850
631.0235221.0238711.0237441.0236491.0234581.0232201.0233591.035844
641.0223181.0231641.0228571.0226261.0221641.0215941.0219471.034195
651.0218401.0228761.0225001.0222171.0216511.0209611.0214231.033596
661.0217771.0228381.0224531.0221631.0215831.0208781.0213541.033517
671.0217631.0228301.0224421.0221511.0215681.0208591.0213381.033500
681.0215971.0227291.0223181.0220091.0213911.0206401.0211531.033289
691.0215961.0227291.0223181.0220091.0213901.0206391.0211531.033289
Node8:009:0010:0011:0012:0013:0014:0015:00
11.0500001.0500001.0500001.0500001.0500001.0500001.0500001.050000
21.0368541.0368621.0243151.0244021.0243471.0369201.0369551.036935
31.0366701.0366861.0242401.0244141.0243041.0368031.0368731.036832
41.0365771.0365871.0241571.0243521.0242321.0367191.0367861.036741
51.0356281.0356131.0234601.0239381.0236601.0359441.0360531.035936
61.0268641.0270411.0189701.0230791.0205731.0298491.0312691.030329
71.0185121.0189361.0149851.0227271.0179681.0241861.0270501.025331
81.0165911.0170811.0141101.0227151.0174181.0229051.0261291.024231
91.0155681.0160821.0135941.0226251.0170661.0221971.0255761.023580
101.0062751.0072581.0098821.0233421.0150011.0163551.0216771.018765
111.0042701.0053121.0088881.0231921.0143341.0149951.0206171.017509
121.0001411.0015991.0081281.0249151.0144481.0128781.0199141.016405
130.9949070.9966631.0063771.0258481.0136801.0097641.0180631.014021
140.9901870.9922791.0051021.0271641.0133341.0071121.0167451.012228
150.9860010.9884721.0043231.0288771.0134261.0049401.0159801.011046
160.9852240.9877671.0041811.0291981.0134451.0045381.0158411.010830
170.9832210.9857681.0030051.0287461.0125641.0030821.0145491.009327
180.9831980.9857451.0029921.0287411.0125551.0030661.0145351.009311
190.9819570.9845321.0023861.0286561.0121501.0022281.0138741.008520
200.9811590.9837531.0019961.0286011.0118901.0016891.0134491.008011
210.9798700.9824941.0013701.0285171.0114741.0008211.0127661.007193
220.9798440.9824701.0013661.0285301.0114761.0008081.0127631.007187
230.9797880.9824551.0015231.0288191.0116711.0008661.0129431.007363
240.9796650.9824221.0018661.0294481.0120961.0009941.0133371.007747
250.9800740.9830601.0031631.0312171.0134981.0018431.0148091.009269
260.9802430.9833231.0036981.0319461.0140761.0021931.0154161.009897
270.9805260.9836691.0041521.0324661.0145321.0025481.0159291.010443
281.0366191.0366311.0241951.0243811.0242661.0367581.0368241.036780
291.0359841.0359461.0236331.0239851.0237891.0361951.0362261.036131
301.0337681.0337121.0222751.0234431.0227661.0345531.0347141.034370
311.0333771.0333181.0220351.0233481.0225861.0342641.0344471.034060
321.0314211.0313461.0208371.0228711.0216841.0328151.0331131.032507
331.0278801.0278291.0189191.0224121.0203421.0303201.0309881.029988
341.0248971.0251331.0185881.0241131.0207121.0288841.0306581.029322
351.0215501.0216021.0157791.0220601.0182811.0260111.0275151.025805
361.0365651.0365731.0241451.0243451.0242231.0367081.0367751.036729
371.0362761.0362471.0238641.0241621.0239951.0364371.0364991.036446
381.0359871.0358521.0235611.0241241.0238361.0361691.0362431.036200
391.0360461.0358881.0235991.0242111.0239021.0362181.0363051.036279
401.0167041.0172211.0143101.0229581.0176241.0230561.0263581.024475
411.0173991.0180431.0153111.0240621.0186161.0238681.0274891.025684
421.0144641.0150111.0131251.0226441.0167731.0214611.0250431.022926
431.0133311.0139251.0127091.0227621.0165441.0207351.0245671.022325
441.0118761.0125431.0122391.0230171.0163251.0198341.0240251.021622
451.0098561.0105001.0109691.0223801.0153071.0182801.0225721.019937
460.9989780.9995481.0043211.0193231.0100961.0100871.0149801.011079
470.9936170.9941501.0010451.0178181.0075281.0060491.0112391.006713
480.9915390.9920570.9997761.0172371.0065341.0044851.0097901.005023
490.9895520.9900850.9987181.0169231.0057621.0030681.0085641.003555
500.9868120.9873450.9972021.0163791.0046151.0010551.0067891.001448
510.9863650.9868740.9968301.0161061.0042931.0006731.0063741.000984
520.9858570.9863400.9964071.0157951.0039261.0002391.0059031.000458
530.9833700.9837230.9943361.0142721.0021270.9981111.0035940.997878
540.9817280.9819970.9929611.0132501.0009280.9967001.0020670.996176
551.0038171.0048341.0085111.0229181.0140091.0146081.0201961.017036
561.0038121.0048281.0085071.0229151.0140051.0146041.0201911.017031
571.0013691.0031361.0101081.0271941.0164491.0144491.0221761.018808
581.0013801.0031501.0101231.0272111.0164641.0144621.0221941.018827
591.0366161.0366261.0241911.0243801.0242631.0367541.0368211.036777
601.0359681.0359091.0236011.0239991.0237831.0361731.0362121.036127
611.0355421.0354711.0233751.0240131.0236581.0358761.0359951.035883
621.0354191.0353441.0233101.0240181.0236211.0357901.0359321.035812
631.0354111.0353361.0233041.0240151.0236171.0357831.0359261.035805
641.0331541.0329971.0218521.0235531.0225911.0340691.0344341.034137
651.0323311.0321461.0213511.0234431.0222551.0334581.0339301.033574
661.0322221.0320341.0212851.0234291.0222111.0333781.0338631.033499
671.0321981.0320091.0212691.0234251.0222011.0333601.0338481.033482
681.0319091.0317081.0210871.0233801.0220781.0331421.0336651.033279
691.0319081.0317071.0210871.0233811.0220781.0331421.0336651.033279
Node16:0017:0018:0019:0020:0021:0022:0023:00
11.0500001.0500001.0500001.0500001.0500001.0500001.0500001.050000
21.0369221.0369011.0368651.0368601.0368401.0368691.0242591.024297
31.0368071.0367651.0366931.0366841.0366441.0367021.0241281.024204
41.0367201.0366711.0366091.0366061.0365661.0366301.0240631.024148
51.0359291.0358081.0357461.0357771.0357061.0358611.0234061.023609
61.0299021.0287951.0274861.0273741.0264931.0278601.0169091.018692
71.0243761.0222191.0194941.0191101.0173441.0199111.0104511.013801
81.0231351.0207221.0176391.0171751.0151901.0180401.0089311.012651
91.0224361.0199081.0166761.0161971.0141171.0171101.0081781.012083
101.0167541.0129941.0076661.0067011.0034331.0079201.0008501.006702
111.0153891.0114231.0057891.0048011.0013481.0061240.9994131.005641
121.0135201.0087021.0014310.9998250.9954891.0010980.9954221.002728
131.0105541.0049540.9963340.9943500.9892370.9957670.9912360.999731
141.0081031.0016950.9916200.9891550.9832240.9906380.9872050.996840
151.0061860.9989440.9873040.9842530.9774610.9857190.9833360.994061
161.0058320.9984350.9865030.9833420.9763900.9848050.9826170.993544
171.0043240.9968040.9847470.9816810.9746380.9833110.9814420.992702
181.0043070.9967860.9847270.9816620.9746170.9832930.9814280.992692
191.0034620.9958230.9835950.9805330.9733840.9822420.9806050.992104
201.0029190.9952040.9828670.9798070.9725910.9815670.9800760.991726
211.0020440.9942070.9816910.9786340.9713080.9804750.9792190.991114
221.0020330.9941900.9816640.9786040.9712720.9804440.9791960.991098
231.0021290.9942190.9815450.9784020.9709940.9802080.9790140.990971
241.0023390.9942830.9812870.9779630.9703870.9796950.9786170.990694
251.0034200.9950350.9812580.9774390.9694670.9789200.9780220.990287
261.0038650.9953450.9812460.9772220.9690880.9786000.9777770.990118
271.0042890.9956890.9813860.9772170.9689810.9785100.9777080.990071
281.0367581.0367151.0366501.0366471.0366091.0366711.0241021.024183
291.0361561.0360941.0361271.0361951.0361791.0362881.0237791.023924
301.0344321.0342621.0342311.0344541.0343741.0347571.0225971.023108
311.0341281.0339391.0338961.0341461.0340561.0344871.0223881.022964
321.0326071.0323231.0322231.0326101.0324641.0331371.0213451.022245
331.0300361.0294931.0290971.0296151.0292741.0304251.0192431.020784
341.0287951.0276111.0260271.0261161.0251441.0269211.0165361.018914
351.0256461.0245511.0234241.0240681.0232781.0253351.0153061.018058
361.0367091.0366581.0365991.0365951.0365561.0366201.0240551.024140
371.0364291.0363401.0363381.0363411.0363011.0363831.0238391.023948
381.0361591.0359281.0360751.0360461.0359871.0361041.0236011.023756
391.0362171.0359421.0361141.0360581.0359871.0361041.0236011.023756
401.0233191.0208631.0176861.0171521.0151191.0179781.0088801.012611
411.0242761.0216611.0180771.0172451.0150181.0178921.0088121.012562
421.0217381.0190631.0156351.0151331.0129371.0160941.0073681.011489
431.0210701.0182201.0145401.0139871.0116461.0149831.0064791.010834
441.0202591.0171621.0131151.0124721.0099221.0135001.0052991.009968
451.0186741.0154701.0113761.0108911.0082911.0121001.0041861.009157
461.0102471.0064971.0020521.0024080.9995171.0046270.9983421.005010
471.0060941.0020750.9974570.9982280.9951931.0009440.9954621.002966
481.0044841.0003620.9956760.9966070.9935170.9995160.9943461.002176
491.0030530.9987880.9939390.9949780.9917870.9980530.9932191.001397
501.0010470.9965710.9915380.9927400.9894200.9960380.9916441.000281
511.0006300.9961640.9911870.9924600.9891630.9958180.9914701.000155
521.0001560.9957010.9907890.9921420.9888700.9955670.9912711.000010
530.9978340.9934330.9888360.9905820.9874360.9943360.9902950.999300
540.9962980.9919320.9875370.9895360.9864700.9935020.9896240.998800
551.0149651.0110101.0054361.0045231.0010941.0059090.9992461.005525
561.0149601.0110051.0054321.0045201.0010911.0059060.9992441.005524
571.0153721.0101971.0019760.9997420.9949611.0006410.9950521.002452
581.0153871.0102091.0019820.9997430.9949601.0006390.9950511.002451
591.0367551.0367101.0366471.0366431.0366051.0366661.0240981.024180
601.0361391.0360421.0361011.0361501.0361251.0362371.0237301.023878
611.0358511.0356661.0356881.0357061.0356331.0358051.0233741.023605
621.0357681.0355581.0355691.0355781.0354901.0356801.0232721.023526
631.0357611.0355511.0355621.0355711.0354831.0356741.0232671.023522
641.0340211.0335281.0335071.0335251.0333071.0337671.0217001.022318
651.0334111.0327931.0327401.0327391.0324571.0330211.0210831.021840
661.0333301.0326961.0326391.0326351.0323451.0329221.0210021.021777
671.0333121.0326741.0326161.0326121.0323201.0329011.0209841.021763
681.0330931.0324111.0323491.0323391.0320271.0326431.0207701.021597
691.0330931.0324111.0323481.0323381.0320261.0326421.0207691.021596

References

  1. Cossent, R.; Gómez, T.; Olmos, L.; Gómez, T. Large-scale integration of renewable and distributed generation of electricity in Spain: Current situation and future needs. Energy Policy 2011, 39, 8078–8087. [Google Scholar] [CrossRef]
  2. Schonardie, M.F.; Martins, D.C. Application of the dq0 transformation in the three-phase grid-connected PV systems with active and reactive power control. In Proceedings of the 2008 IEEE International Conference on Sustainable Energy Technologies, Singapore, 24–27 November 2008; pp. 18–23. [Google Scholar]
  3. Reference Technical Rules for the Connection of Active and Passive Consumers to the HV and MV Electrical Networks of Distribution Company Issue; Italian Technical Guidelines CEI 0-16 2012-12; Comitato Elettrotecnico Italiano: Milan, Italy, 2012.
  4. Thomson, M.; Infield, D. Impact of widespread photovoltaics generation on distribution systems. IET Renew. Power Gener. 2007, 1, 33. [Google Scholar] [CrossRef]
  5. Tonkoski, R.; Lopes, L.A.C.; El-Fouly, T.H.M. Coordinated Active Power Curtailment of Grid Connected PV Inverters for Overvoltage Prevention. IEEE Trans. Sustain. Energy 2010, 2, 139–147. [Google Scholar] [CrossRef]
  6. Zhao, C.; Gu, C.; Li, F.; Dale, M. Understanding LV network voltage distribution- UK smart grid demonstration experience. In Proceedings of the 2015 IEEE Power & Energy Society Innovative Smart Grid Technologies Conference (ISGT), Washington, DC, USA, 18–20 February 2015; pp. 1–5. [Google Scholar]
  7. Dou, X.; Duan, X.; Hu, Q.; Shen, L.; Wu, Z. A nonintrusive control strategy using voltage and reactive power for distribution systems based on PV and the nine-zone diagram. Int. J. Electr. Power Energy Syst. 2019, 105, 89–97. [Google Scholar] [CrossRef]
  8. Zhu, Y.; Tomsovic, K. Adaptive power flow method for distribution systems with dispersed generation. IEEE Trans. Power Deliv. 2002, 17, 822–827. [Google Scholar] [CrossRef]
  9. Valverde, G.; Van Cutsem, T.H. Model Predictive Control of Voltages in Active Distribution Networks. IEEE Trans. Smart Grid 2013, 4, 2152–2161. [Google Scholar] [CrossRef]
  10. Robbins, B.A.; Hadjicostis, C.N.; Domínguez-García, A.D. A Two-Stage Distributed Architecture for Voltage Control in Power Distribution Systems. IEEE Trans. Power Syst. 2012, 28, 1470–1482. [Google Scholar] [CrossRef]
  11. Sulc, P.; Backhaus, S.; Chertkov, M. Optimal Distributed Control of Reactive Power via the Alternating Direction Method of Multipliers. IEEE Trans. Energy Convers. 2014, 29, 968–977. [Google Scholar] [CrossRef] [Green Version]
  12. Calderaro, V.; Galdi, V.; Massa, G.; Piccolo, A. Distributed generation management: An optimal sensitivity approach for decentralized power control. In Proceedings of the 2012 3rd IEEE PES Innovative Smart Grid Technologies Europe (ISGT Europe), Berlin, Germany, 14–17 October 2012; pp. 1–8. [Google Scholar] [CrossRef]
  13. Calderaro, V.; Conio, G.; Galdi, V.; Massa, G.; Piccolo, A. Optimal Decentralized Voltage Control for Distribution Systems with Inverter-Based Distributed Generators. IEEE Trans. Power Syst. 2013, 29, 230–241. [Google Scholar] [CrossRef]
  14. Turitsyn, K.; Sulc, P.; Backhaus, S.; Chertkov, M. Options for Control of Reactive Power by Distributed Photovoltaic Generators. Proc. IEEE 2011, 99, 1063–1073. [Google Scholar] [CrossRef] [Green Version]
  15. Cai, Y.; Tang, W.; Li, L.; Zhang, B.; Zhang, L.; Wang, Y. Multi-mode adaptive local reactive power control method based on PV inverters in low voltage distribution networks. IET Gener. Transm. Distrib. 2020, 14, 542–551. [Google Scholar] [CrossRef]
  16. Molina-Garcia, A.; Mastromauro, R.A.; Garcia-Sanchez, T.; Pugliese, S.; Liserre, M.; Stasi, S. Reactive Power Flow Control for PV Inverters Voltage Support in LV Distribution Networks. IEEE Trans. Smart Grid 2016, 8, 447–456. [Google Scholar] [CrossRef] [Green Version]
  17. Dickert, J.; Domagk, M.; Schegner, P. Benchmark low voltage distribution networks based on cluster analysis of actual grid properties. In Proceedings of the 2013 IEEE Grenoble Conference, Grenoble, France, 16–20 June 2013; pp. 1–6. [Google Scholar]
  18. Kabir, M.N.; Mishra, Y.; Ledwich, G.; Dong, Z.Y.; Wong, K.P. Coordinated Control of Grid-Connected Photovoltaic Reactive Power and Battery Energy Storage Systems to Improve the Voltage Profile of a Residential Distribution Feeder. IEEE Trans. Ind. Inform. 2014, 10, 967–977. [Google Scholar] [CrossRef]
  19. Mumtaz, F.; Syed, M.H.; Al Hosani, M.; Zeineldin, H.H. A Novel Approach to Solve Power Flow for Islanded Microgrids Using Modified Newton Raphson with Droop Control of DG. IEEE Trans. Sustain. Energy 2015, 7, 493–503. [Google Scholar] [CrossRef] [Green Version]
  20. Biserica, M.; Berseneff, B.; Besanger, Y.; Kieny, C. Upgraded Coordinated Voltage Control for Distribution Systems. In Proceedings of the IEEE Trondheim PowerTech, Trondheim, Norway, 19–23 June 2011; pp. 1–6. [Google Scholar] [CrossRef]
  21. Zhang, B.; Lam, A.Y.S.; Dominguez-Garcia, A.D.; Tse, D. An Optimal and Distributed Method for Voltage Regulation in Power Distribution Systems. IEEE Trans. Power Syst. 2014, 30, 1714–1726. [Google Scholar] [CrossRef] [Green Version]
  22. Ahmed, M.; Bhattarai, R.; Hossain, S.J.; Abdelrazek, S.; Kamalasadan, S. Coordinated voltage control strategy for voltage regulators and voltage source converters integrated distribution system. In Proceedings of the 2017 IEEE Industry Applications Society Annual Meeting, Harbor, MD, USA, 14–17 January 2017; pp. 1–8. [Google Scholar]
  23. Castillo, A.; Gayme, D.F. Profit maximizing storage allocation in power grids. In Proceedings of the 52nd IEEE Conference on Decision and Control, Florence, Italy, 10–13 December 2013; pp. 429–435. [Google Scholar]
  24. Bose, S.; Gayme, D.F.; Topcu, U.; Chandy, K.M.; Bose, S. Optimal placement of energy storage in the grid. In Proceedings of the 2012 IEEE 51st IEEE Conference on Decision and Control (CDC), Maui, HI, USA, 10–13 December 2012; pp. 5605–5612. [Google Scholar]
  25. Nasrolahpour, E.; Kazempour, S.; Zareipour, H.; Rosehart, W.D. Strategic Sizing of Energy Storage Facilities in Electricity Markets. IEEE Trans. Sustain. Energy 2016, 7, 1462–1472. [Google Scholar] [CrossRef] [Green Version]
  26. Nick, M.; Cherkaoui, R.; Paolone, M. Optimal siting and sizing of distributed energy storage systems via alternating direction method of multipliers. Int. J. Electr. Power Energy Syst. 2015, 72, 33–39. [Google Scholar] [CrossRef]
  27. Yan, W.; Yu, J.; Yu, D.; Bhattarai, K. A new optimal reactive power flow model in rectangular form and its solution by predictor corrector primal dual interior point method. IEEE Trans. Power Syst. 2006, 21, 61–67. [Google Scholar] [CrossRef]
  28. Esmin, A.; Lambert-Torres, G.; De Souza, A.Z. A hybrid particle swarm optimization applied to loss power minimization. IEEE Trans. Power Syst. 2005, 20, 859–866. [Google Scholar] [CrossRef]
  29. Li, Y.; Li, M.; Wu, Q. Optimal reactive power dispatch with wind power integrated using group search optimizer with intraspecific competition and lévy walk. J. Mod. Power Syst. Clean Energy 2014, 2, 308–318. [Google Scholar] [CrossRef] [Green Version]
  30. Farivar, M.; Low, S.H. Branch Flow Model: Relaxations and Convexification—Part I. IEEE Trans. Power Syst. 2013, 28, 2554–2564. [Google Scholar] [CrossRef]
  31. Wang, C.; Xiang, L.; Deng, Z.; Liang, Z. Reactive Power Optimization Dispatch Strategy in Distribution Network with Distributed Generators. J. Eng. 2017, 2017, 1418–1422. [Google Scholar] [CrossRef]
  32. Chen, L.; Deng, Z.; Xu, X. Two-Stage Dynamic Reactive Power Dispatch Strategy in Distribution Network Considering the Reactive Power Regulation of Distributed Generations. IEEE Trans. Power Syst. 2018, 34, 1021–1032. [Google Scholar] [CrossRef]
  33. Maciejowski, J.M. Predictive Control with Constraints; Prentice-Hall: Englewood Cliffs, NJ, USA, 2002. [Google Scholar]
  34. Glavic, M.; Hajian, M.; Rosehart, W.; Van Cutsem, T.H. Receding-Horizon Multi-Step Optimization to Correct Nonviable or Unstable Transmission Voltages. IEEE Trans. Power Syst. 2011, 26, 1641–1650. [Google Scholar] [CrossRef] [Green Version]
  35. Fortenbacher, P.; Ulbig, A.; Andersson, G. Optimal Placement and Sizing of Distributed Battery Storage in Low Voltage Grids Using Receding Horizon Control Strategies. IEEE Trans. Power Syst. 2017, 33, 2383–2394. [Google Scholar] [CrossRef] [Green Version]
  36. Hug-Glanzmann, G. Coordination of intermittent generation with storage, demand control and conventional energy sources. In Proceedings of the 2010 IREP Symposium Bulk Power System Dynamics and Control—VIII (IREP), Rio de Janeiro, Brazil, 1–6 August 2010; pp. 1–7. [Google Scholar] [CrossRef]
  37. Bolognani, S.; Zampieri, S. On the Existence and Linear Approximation of the Power Flow Solution in Power Distribution Networks. IEEE Trans. Power Syst. 2015, 31, 163–172. [Google Scholar] [CrossRef] [Green Version]
  38. Yang, Z.; Wu, R.; Yang, J.; Long, K.; You, P. Economical Operation of Microgrid with Various Devices Via Distributed Optimization. IEEE Trans. Smart Grid 2015, 7, 1–11. [Google Scholar] [CrossRef]
  39. Nick, M.; Cherkaoui, R.; Paolone, M. Optimal Allocation of Dispersed Energy Storage Systems in Active Distribution Networks for Energy Balance and Grid Support. IEEE Trans. Power Syst. 2014, 29, 2300–2310. [Google Scholar] [CrossRef]
  40. Niknam, T.; Zare, M.; Aghaei, J. Scenario-Based Multiobjective Volt/Var Control in Distribution Networks Including Renewable Energy Sources. IEEE Trans. Power Deliv. 2012, 27, 2004–2019. [Google Scholar] [CrossRef]
  41. Baran, M.; Wu, F. Optimal capacitor placement on radial distribution systems. IEEE Trans. Power Deliv. 1989, 4, 725–734. [Google Scholar] [CrossRef]
  42. GUROBI. Available online: http://www.gurobi.com (accessed on 22 February 2018).
  43. National Bureau of Quality Supervision. Inspection and Quarantine and National Standardization Management Committee. National Standard of the People’s Republic of China GB/T12325-2008 Power Quality and Voltage Deviation; National Bureau of Quality Supervision: Beijing, China, 2008.
  44. Zhang, Y.-J.; Ren, Z. Optimal Reactive Power Dispatch Considering Costs of Adjusting the Control Devices. IEEE Trans. Power Syst. 2005, 20, 1349–1356. [Google Scholar] [CrossRef]
  45. Wong, T. Co-ordination of transformer tap and capacitor operation for reactive power voltage control in a distribution primary substation. In Proceedings of the APSCOM 2000—5th International Conference on Advances in Power System Control, Operation and Management, Hong Kong, China, 30 October–1 November 2000; Volume 2000, pp. 479–485. [Google Scholar]
  46. Sun, C.; Xu, J.; Cai, B.; Yang, Z.; Li, X. Research on efficiency increase of the automatic voltage control system in Shanghai urban power grid. In Proceedings of the 8th Renewable Power Generation Conference (RPG 2019), Shanghai, China, 24–25 October 2019; pp. 1–6. [Google Scholar]
  47. Khan, I.; Li, Z.; Xu, Y.; Gu, W. Distributed control algorithm for optimal reactive power control in power grids. Int. J. Electr. Power Energy Syst. 2016, 83, 505–513. [Google Scholar] [CrossRef]
Figure 1. Comparison of the reactive power and voltage optimization strategy framework.
Figure 1. Comparison of the reactive power and voltage optimization strategy framework.
Energies 13 05789 g001
Figure 2. The detailed decomposition diagram of DRM strategy.
Figure 2. The detailed decomposition diagram of DRM strategy.
Energies 13 05789 g002
Figure 3. Flowchart of the proposed DRM strategy to solve the reactive power and voltage problem.
Figure 3. Flowchart of the proposed DRM strategy to solve the reactive power and voltage problem.
Energies 13 05789 g003
Figure 4. Modified PG&E 69-bus system with PVs.
Figure 4. Modified PG&E 69-bus system with PVs.
Energies 13 05789 g004
Figure 5. The aggregate daily load profile and solar power profile of Case 1.
Figure 5. The aggregate daily load profile and solar power profile of Case 1.
Energies 13 05789 g005
Figure 6. The nodal voltages of 11:00 and 24:00.
Figure 6. The nodal voltages of 11:00 and 24:00.
Energies 13 05789 g006
Figure 7. The daily voltages of Node16 and Node54.
Figure 7. The daily voltages of Node16 and Node54.
Energies 13 05789 g007
Figure 8. OLTC tap changes with DRM strategy.
Figure 8. OLTC tap changes with DRM strategy.
Energies 13 05789 g008
Figure 9. PV reactive power outputs with DRM strategy.
Figure 9. PV reactive power outputs with DRM strategy.
Energies 13 05789 g009
Figure 10. PV active power outputs with DRM strategy.
Figure 10. PV active power outputs with DRM strategy.
Energies 13 05789 g010
Figure 11. PV3 reactive power outputs with three strategies.
Figure 11. PV3 reactive power outputs with three strategies.
Energies 13 05789 g011
Figure 12. OLTC tap changes with three strategies.
Figure 12. OLTC tap changes with three strategies.
Energies 13 05789 g012
Figure 13. The nodal voltages of 11:00 and 24:00 from three strategies.
Figure 13. The nodal voltages of 11:00 and 24:00 from three strategies.
Energies 13 05789 g013
Figure 14. The daily voltages of Node16 and Node54 from three strategies.
Figure 14. The daily voltages of Node16 and Node54 from three strategies.
Energies 13 05789 g014
Figure 15. Modified PG&E 69-bus system with uneven PVs.
Figure 15. Modified PG&E 69-bus system with uneven PVs.
Energies 13 05789 g015
Figure 16. The aggregate daily load profile and forecasted PV power profile of Case 2.
Figure 16. The aggregate daily load profile and forecasted PV power profile of Case 2.
Energies 13 05789 g016
Figure 17. The nodal voltages at 11:00 and 24:00 of Case 2.
Figure 17. The nodal voltages at 11:00 and 24:00 of Case 2.
Energies 13 05789 g017
Figure 18. The daily voltages at Node27 and Node35 of Case 2.
Figure 18. The daily voltages at Node27 and Node35 of Case 2.
Energies 13 05789 g018
Figure 19. OLTC tap changes with DRM strategy of Case 2.
Figure 19. OLTC tap changes with DRM strategy of Case 2.
Energies 13 05789 g019
Figure 20. PV reactive power outputs with DRM strategy of Case 2.
Figure 20. PV reactive power outputs with DRM strategy of Case 2.
Energies 13 05789 g020
Figure 21. PV active power outputs with DRM strategy of Case 2.
Figure 21. PV active power outputs with DRM strategy of Case 2.
Energies 13 05789 g021
Table 1. Summary of the literature related to ADN optimization method.
Table 1. Summary of the literature related to ADN optimization method.
DevicesReferenceMethod/AlgorithmRemarks
OLTC[7]Mixed integer nonlinear programming (MINP)
  • Disadvantages: due to physical constraints and low cost benefit ratios, OLTC cannot change flexibly to achieve real-time adjustment
DG[12,13]Decentralized control based on voltage sensitivity
  • Advantages: The continuous, quick and frequent reactive power outputs of DGs avoids adding extra voltage devices and reactive compensators.
  • Disadvantages: DG reactive outputs are limited by capacity and power factor. It is insufficient to correct serious violations by DGs alone.
[15]A multi-mode adaptive local control method based on Q(P) characteristics
[16]a combination of centralized remote control and local decentralized control
DG, OLTC[20,21,22]Optimal power flow (OPF) calculation
  • Disadvantages: These solutions are the static optimization results, regardless of the strong time coupling between the adjustments of different devices.
[23,24]Semi-definite relaxation
[25,26]Benders decomposition
[27,28,29,30]Relaxing discrete variables into continuous variables and rounding to close integers.
[31,32]Two-stage dynamic coordination dispatch
  • Advantages: the dynamic optimization results.
  • Disadvantages: The pre-planned scheduling may have large errors and lose dynamic optimality.
Storage device, OLTC, DG[9,35,36]Model predictive control (MPC)
  • Advantages: Rolling optimization of the specific horizon can achieve a dynamic transition from the measured to the target values.
  • Disadvantages: The requirement of subdividing the control horizon for different devices is ignored.
Table 2. The basic installation parameters of OLTC and PVs.
Table 2. The basic installation parameters of OLTC and PVs.
DevicesParametersValues
OLTCTap range−8–8
PVCapacity500 kWP
Table 3. Comparison of three voltage optimization strategies.
Table 3. Comparison of three voltage optimization strategies.
CharacteristicS1S2DRM
Decoupling of discrete/continuous device actions
Subdividing optimization periods×
MPC-based rolling optimization××
Table 4. Adjustment costs of three strategies.
Table 4. Adjustment costs of three strategies.
Adjustment CostS1S2DRM
OLTC$24$24$24
PVs$1.17$1.64$0.63
Total$25.17$25.64$24.63
Table 5. Voltage deviation results of three strategies.
Table 5. Voltage deviation results of three strategies.
OriginalS1S2DRM
All nodes at 11:003.4951 p.u.1.6899 p.u.1.6806 p.u.1.5839 p.u.
All nodes at 24:002.4852 p.u.0.9467 p.u.0.9467 p.u.0.9467 p.u.
Node 16 all day1.9247 p.u.1.0189 p.u.1.0010 p.u.0.9740 p.u.
Node 54 all day1.5351 p.u.0.6928 p.u.0.6899 p.u.0.6881 p.u.
All nodes all day299.9501 p.u.122.7256 p.u.122.0783 p.u.121.8496 p.u.
Table 6. Active power losses from three strategies.
Table 6. Active power losses from three strategies.
No PVOriginalS1S2DRM
Active power losses4.8305 MWh2.0546 MWh2.7143 MWh2.6610 MWh2.6472 MWh
Table 7. Computation time of three strategies.
Table 7. Computation time of three strategies.
S1S2DRM
Long period (1 h) optimization21.46 s48.94 s50.23 s
Short period (15 min) PV optimization-9.07 s10.43 s
One day optimization5.33 min11.64 min14.39 min
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ge, X.; Shen, L.; Zheng, C.; Li, P.; Dou, X. A Decoupling Rolling Multi-Period Power and Voltage Optimization Strategy in Active Distribution Networks. Energies 2020, 13, 5789. https://doi.org/10.3390/en13215789

AMA Style

Ge X, Shen L, Zheng C, Li P, Dou X. A Decoupling Rolling Multi-Period Power and Voltage Optimization Strategy in Active Distribution Networks. Energies. 2020; 13(21):5789. https://doi.org/10.3390/en13215789

Chicago/Turabian Style

Ge, Xiaohui, Lu Shen, Chaoming Zheng, Peng Li, and Xiaobo Dou. 2020. "A Decoupling Rolling Multi-Period Power and Voltage Optimization Strategy in Active Distribution Networks" Energies 13, no. 21: 5789. https://doi.org/10.3390/en13215789

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