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

: With the increasing penetration of distributed photovoltaics (PVs) in active distribution networks (ADNs), the risk of voltage violations caused by PV uncertainties is signiﬁcantly 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 di ﬀ erent devices. The mixed-integer voltage optimization model is ﬁrst 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 simpliﬁed 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.


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.  [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.
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.

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, The modified Newton Raphson power flow equation can be simplified as: where V is the diagonal matrix of voltage amplitude, V = diag(v 1 , v 2 , . . . , v n-1 ), ∆θ and ∆V are the corrections of the voltage phase angle and voltage amplitude respectively. The Jacobian matrix J is as follows: where Bcosθ can be expressed as Equation (3) and the structures of other terms B sin θ, Gcosθ, and Gsinθ are similar.
B n,1 cosθ n,1 · · · B n,j cosθ n,j · · · B n,n cosθ n,n In the distribution network, the value of θ ij is extremely small, and hence Equation (1) can be further simplified: where P and Q are diagonal matrices, and their diagonal elements are P i /V 2 i and Q i /V 2 i respectively. Using the Gaussian elimination method for Equation (4), the sensitivity of nodal voltages and injected power can be calculated as: Therefore, the active-voltage sensitivity and reactive-voltage sensitivity can be obtained:

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 tap 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: where ∆V d is voltage variation on the low-voltage side after a OLTC tap change.

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: P l = S m,PV cosϕ max (11) where, P t m,PV is the PV-m active power output. Q t,max m,PV is the maximum reactive power. S m,PV is the capacity of the PV-m, and ϕ max is the maximum power factor angle.

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.

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.

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: where P PV and P PV are the actual output and predicted output of PV active power. ∆N tap and ∆Q PV are the variations of OLTC tap and PVs reactive power respectively. c b and c PV 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 re 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: 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: In general, the objective function voltage optimization is established based on the chronological scenario: where t l and t s 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.

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.

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 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.

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

OLTC Tap Constraints
The OLTC tap range as a key parameter should be considered, as well as the amplitude of one tap change: where N max tap and N min tap are the lower limit and upper limit of the OLTC tap position.

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: (22,23) Equations(10), (11), (24) (27) where λ k is a dual variable generated by constraint (∆N t l tap The first line of the constraints is the form of Benders cut, where α is an auxiliary variable. α down 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.

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 k B 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 k B 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:

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.

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.  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.

Performance Comparison
Adjustment Costs PV reactive power outputs (for example, PV3) and OLTC tap actions of the three mentioned optimization strategies are shown in Figures 11 and 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 13. The nodal voltages of 11:00 and 24:00 from three strategies. 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. 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.

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 Figures 19-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.

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.