Two-Stage Active and Reactive Power Coordinated Optimal Dispatch for Active Distribution Network Considering Load Flexibility

: A high proportion of renewable energy connected to the power grid has caused power quality problems. Voltage-sensitive loads are extremely susceptible to voltage ﬂuctuations, causing power system safety issues and economic losses. Considering the uncertainty factor and the time-varying characteristic, a linearized random ZIP model (constant impedance (Z), constant current (I), and constant power (P)) with time-varying characteristics was proposed. In order to improve the voltage quality of the voltage-sensitive loads in the day-here stage in an active distribution network (ADN), a linearized two-stage active and reactive power coordinated stochastic optimization model was established. The day-ahead active and reactive power coordination optimization was to smooth the large voltage ﬂuctuation and develop a reserve plan to eliminate the unbalanced power caused by the prediction error in the day-here optimization. In the day-here real-time redispatch, the voltage was further improved by the continuous reactive power compensation device. Finally, the simulation results on the IEEE-33 bus system showed that the control strategy could better eliminate the unbalanced power caused by the prediction error and obviously improve the voltage of sensitive loads in the real-time stage on the premise of maintaining economic optimality.


Introduction
The rapid development worldwide has increased the demand for power construction. With the gradual expansion of the energy crisis and severe environmental problems, sustainable development has become the main goal of future power grid construction. The application of renewable energy source (RES) and distributed energy has developed rapidly. However, with the large amount of RESs connected to the grid, its volatilities and uncertainties have brought huge challenges to the operation and regulation to the active distribution network (ADN) [1]. Optimization and dispatch of the traditional ADN is carried out according to the method of generation load tracking. When a large proportion of RESs are connected, the ADN needs to reserve a large amount of power to eliminate the effect of these uncertainties, which not only reduces the efficiency of generators, but also increases the dispatching cost of the ADN [2,3]. Therefore, it is necessary to utilize the initiative and flexibility of the flexible load on the demand side. The economic cost of the distribution network for an ADN is reduced by optimizing the dynamic active and reactive power.
RESs will cause a certain degree of voltage fluctuations after being connected to the grid. The penetration rate of RESs is greater, and the voltage fluctuations within a short period are greater. Most of the equipment and electrical loads in the semiconductor manufacturing, information, computer, or electronic communications industries are highly sensitive loads [4]. The operating state of many industries like electronic communications, adjustable speed motors, direct current (DC) motor drives, alternating current (AC) contactors, programmable logic, and other equipment are affected by voltage fluctuations, which cause huge economic losses [5,6]. These highly sensitive loads have higher requirements for voltage quality and reactive power compensation [7,8]. Different types of highly sensitive loads do not have the same reference to sensitivity in power quality problems, which have great uncertain characteristics. Through the above analysis, it is necessary to comprehensively consider the flexibility of flexible loads and the loss cost of highly sensitive loads, and set up a reasonable active and reactive power optimal dispatch method to ensure the stability and the reasonable economic cost of an ADN.
There have been some studies on the optimal dispatch and voltage control of a distribution network. In [9], the division of an autonomous region active power control strategy was proposed to improve the voltage control ability of an ADN and make it work normally under large power change. In [10], a chance constrained programming (CCP) model based on the k-means clustering method was proposed considering the uncertainty of a multi-stage planning load and wind power in the ADN. In [11], a hierarchical coordination optimization model considering prediction error to eliminate risks based on chance-constrained programming was proposed to minimize the cost of exchange companies and achieve the best trend. In [12], the stochastic optimal Benders decomposition method considering the uncertainties of flexible load and intermittent energy in the distribution network was established for achieving optimal operation and efficient utilization of intermittent energy in the worst case. In [13], an optimization planning model considering network loss and system reliability was established based on the Pareto optimization method, which realizes the multi-objective optimization of ADN. In [14], a new distributed voltage control strategy based on a multi-agent system platform was proposed that could realize the overall network voltage management and maximize the distributed generation. In [15], multi-objective optimization algorithm models were applied to realize the economic and optimization of the distribution system to reduce the load loss and the cost of network investment. The above literature did not consider the impact of highly sensitive loads in an ADN, nor did it perform real-time regulation and dispatch based on the flexibility and randomness of the load. Therefore, on the basis of ensuring the stability of the highly sensitive load, two-stage optimization and dispatching of an ADN was carried out to optimize the economic cost of the distribution network.
In this paper, a linearized random voltage-sensitive load model was proposed to contain the time characteristics and the uncertainties of the load. A linearized two-stage active and reactive power stochastic coordinated optimization model was established to reduce the influence of voltage fluctuations and uncertainties on these sensitive loads. The day-ahead active and reactive power coordination optimization was to smooth large voltage fluctuations and develop a reserve plan to eliminate the unbalanced power caused by the prediction error in the day-here optimization. In the day-here real-time redispatch, the voltage was further improved by the continuous reactive power compensation devices. Finally, the simulation on the IEEE-33 bus system showed that the proposed coordinated optimal dispatch could effectively improve the voltage level of the sensitive load in the closed-loop method and minimize the operation cost of an ADN.
The rest of this paper is organized as follows. Section 2 presents the linearized random ZIP model (constant impedance (Z), constant current (I), and constant power (P)) considering the time-varying characteristics. In Section 3, the two-stage active and reactive power coordinated optimal dispatch for ADN are given. Simulation results in a modified IEEE 33-bus system are presented in Section 4. Finally, Section 5 presents the conclusion.

Linearized Random ZIP Model Considering Time-Varying Characteristics
The load in the power system is not all constant, and most of the load will be affected by the operating state of the system, especially the voltage of the bus [16]. In this paper, we mainly considered the static characteristics of the load, which can be described by the ZIP model [17]. The ZIP model divides the entire load into three parts including the constant impedance, constant power, and constant current, the traditional ZIP model cannot describe the random characteristics and the time-varying characteristics [18,19]. Therefore, an improved ZIP model was proposed, which is shown as follows: where a P,t , b P,t , c P,t , a Q,t , b Q,t , c Q,t represent the ratio of constant impedance a (CI), constant power b (CP), and constant current c (CC) in this load at time t, respectively, and a P,t + b P,t + c P,t = a Q,t + b Q,t + c Q,t = 1. This paper set a P,t = a Q,t , b P,t = b Q,t , c P,t = c Q,t ; P * m,i,t and Q * m,i,t are the predicted active value and reactive value of the load at bus i in time t, respectively; v i,t is the voltage of bus i in time t; v o is the unit value of power system voltage, v o = 1.0 p.u.; P m,i,t and Q m,i,t are the actual active/reactive power of load at bus i in time t, respectively; and δ P,m,i,t and δ Q,m,i,t are the uncertainty of the active/reactive power demand at bus i in time t.
Consider that there is a square term of voltage in Equations (1) and (2), which builds a nonlinear non-convex model. In the actual operation of the power system, the voltage is very close to the standard unit value 1.00 p.u. Therefore, this paper carried out Taylor's first-order expansion on the square term of voltage at v o = 1.00 p.u. and ignored the higher-order terms. Obtain v i 2 = 2v i − 1 and substitute it in Equations (1) and (2) as: Compared with the traditional ZIP model, the proposed improved ZIP model can well express the static voltage, time-varying, and random characteristics of sensitive loads.

Two-Stage Active and Reactive Power Coordinated Optimal Dispatch for an Active Distribution Network (ADN)
Considering that day-ahead scheduling formulates instructions based on predicted values, these predicted values have a certain degree error. Therefore, it is necessary to introduce the active and reactive power redispatch in the day-here optimization for eliminate the predicted error. Due to the predicted error, pursuing voltage quality in the day-ahead scheduling stage will increase the control cost and reduce the robustness of the power system. Therefore, this paper proposed a two-stage active and reactive power coordinated control strategy for ADN, which takes the efficient and reliable power supply of highly sensitive loads as the goal. The objective function of the proposed model was the sum of the day-ahead scheduling cost and the expected value in the real-time scheduling cost. The OLTC (on-load tap changer) was used to smooth out large voltage fluctuations for sensitive loads in the day-ahead scheduling stage [20]. In the real-time optimization stage, the continuous reactive power control devices distributed in the system were used to further control the bus voltage [21].

Objective Function
From the perspective of power grid operation, the dispatching of an ADN should pursue the maximization of social value and minimize the operating costs. Load utility represents the benefits brought to users by demand-side electricity consumption. Non-flexible loads do not have load utility because they do not have elastic utility. Only the flexible load is considered to provide backup increased load utility costs. The two-stage active and reactive power coordination optimization model proposed in this paper was based on the minimum operating cost of an ADN, which includes day-ahead scheduling costs and the expected value of the real-time scheduling costs. The day-ahead total costs include power purchase costs from the main system, fuel costs, and backup costs of the DG (distribution generator) power generation, backup costs of the flexible load, and day-ahead power network loss costs. The spinning reserve of flexible load and the DG determined by day-ahead are mainly used to balance the power imbalance caused by the prediction error in the real-time scheduling. Real-time scheduling costs mainly include the expected cost of the sum power network loss costs and highly sensitive load loss costs [22,23].
The objective function is shown in Equation (5): P + m,i,t,r − P − m,i,t,r = P m,max,r − P m,i,t,r In the above equations, t, i, g, e, l, m, r represent the indices of time, bus, generator, main grid, flexible load, sensitive load, and the random scene, respectively; N t represents the set of time period; Ω i represents the set of bus; Ω g,i represents the set of generator at bus i; Ω l,i represents the set of the flexible load at bus i; Ω e,i represents the set of the main grid at bus i; Ω line represents the set of line; and Ω m,i represents the set of the sensitive load at bus i.
Day-ahead dispatch cost F DA includes four parts: power generation and reserve cost, flexible load reserve cost, power purchase cost, and day-ahead power network loss cost. C g , C U g , and C D g , respectively, represent the marginal cost, up/down spinning reserve cost of power generator g; C B g and C S g represent the startup and shutdown costs of power generator g, respectively; α gn g,i,t and β gn g,i,t are the startup and shutdown binary integer commands of power generator g at bus i in time period t, respectively; P g,i,t is the output power of power generator g at bus i in time period t; P U g,i,t and P D g,i,t are the up-and down-reserved capacity of power generator g at bus i in time period t, respectively; C U l,t and C D l,t are the up/down spinning reserve cost of the flexible load l in time period t, respectively; P U l,i,t and P D l,i,t are the up/down spinning reserve of the flexible load l at bus i in time period t, respectively; P e,i,t and Q e,i,t are the active and reactive power purchased from the main grid at bus i in time period t, respectively; C e,p and C e,q are the unit purchased price of active power and reactive power, respectively; and C A loss and P A loss,ij,t are the unit electricity price of day-ahead power network loss and the day-ahead power network loss of line i j in time period t.
The expected real-time cost F DI is the product of the real-time cost of each scene and the scene probability π r ; C I loss and P I loss,ij,t,r are the unit electricity price of real-time power network loss and the real-time power network loss of line i j in time period t under the scene r, respectively; R is a set of real-time scenes. The set in this paper was 5. C m,t is the unit price of the highly sensitive load loss; and ∆P m,i,t,r is the real-time loss of a highly sensitive load in the time period t under the scene r.
In Equations (8) and (9), g ij is the conductance between bus i and bus j in the time period t; m ij,y is the linearized parameter that will be explained in detail later; Y variables are those needed to model the domain restriction of the function to the positive orthant. θ ij,y,t and θ ij,y,r,t are the segmented value of the voltage phase angle difference between the two ends of the bus i and bus j in the time period t in the day-head optimization or in the scenes r of the day-here optimization, respectively. In Equations (10)-(13), two slack variables P + m,i,t,r and P − m,i,t,r are introduced to obtain the change of highly sensitive load loss. P m,max,r and P m,i,t,r are the active power values of the high-sensitive load when the voltage is close to 1.0 p.u and the actual value of the high-sensitive load in time period t under the scene r.

1.
Constraints for the generator Equations (15) and (16) are the minimum startup time constraint of the generators; Equations (17) and (18) are the minimum shutdown time constraint of the generators; Equations (19) and (20) are the up-and down-ramp constraints of the generators; Equations (21) and (22) are the output active capacity constraints of the generators; Equation (23) is the reactive power constraint of the generators; and Equations (24) and (25) are the range of the upper and lower reserve of the generators. 2.
Constraints for the OLTC Equation (26) shows the voltage at both ends of the OLTC; Equation (27) shows the constraints of the tap; and Equation (28) shows the constraints of the allowable number of OLTC operations.
where T ij is the transformation ratio of OLTC, assuming that OLTC has an H tap, the transformation ratio is T ij ∈ { T ij,1 , T ij,2 , · · · , T ij,H }. In order to express the transformation ratio of the hth tap, a binary variable H is introduced at time t, which is β ij,1,t , β ij,2,t , · · · , β ij,H,t [24]. N T is the maximum allowable number of OLTC operations, which was set to 4 in this paper.

Constraints for the flexible load
The flexible load participates in the day-ahead dispatch of the power system by providing a reserve. Constraints (33) and (34) show the amount of up/down spinning reserve capability of the flexible load [25].
where P U l,i,t,max and P D l,i,t,max are the maximum schedulable reserve up and reserve down for the flexible load of the bus i in the time period t, respectively. P U l,i,t and P D l,i,t represent the up/down spinning reserve at the bus i in time, respectively.

4.
Constraints for the SVC Static var compensator (SVC) can not only output reactive power to the power grid, but also absorb reactive power from the power grid and realize continuous reactive power regulation. In order to simplify the model, the reactive power output of SVC connected to bus i in the time period t is expressed as Q s,i,t [24].
where Q s,i,max and Q s,i,min are respectively the upper and lower limits of output reactive power of the SVC connected to bus i, which represent the out reactive power of SVC at bus i in the time period t.

5.
Constraints for the output active/reactive power of RES The renewable energy considered in this paper was mainly wind power, which participates in the active and reactive power coordination optimization [26,27]. The output active power of the wind turbine is expressed as: where P w,i,t and P w o ,i are the output power of the wind turbine at the bus i in the time period t and rated power of the wind turbine at the bus i, respectively; v w,ip , v w,op , and v w o are the cut-in, cut-out, and rated wind speed of the wind turbine, respectively. a and b are the parameters related to the cut-in rated wind speed, which are fixed constants: a = v w,ip Wind turbines are subject to the maximum apparent capacity of the inverter during operation, which is shown as follows: where S w,i,max is the maximum apparent capacity and Q w,i,t is the output reactive power of the wind turbines at bus i in time t.
In order to simplify the calculation, linearize Equation (37) by using the hexagonal relaxation method: − 6.
The day-ahead power flow balance constraints For the radial ADN, Equations (37) and (38) are the power flow constraints of the system: In Equations (37) and (38), N is the total number of buses in the system; i and j are the bus numbers; Ω line,i is the set of lines connected with bus i; P ij,t and Q ij,t are the total active/reactive power of the line i j at time t; P w,i,t , P l,i,t , P m,i,t are the active power of the wind turbine injected into bus i in time t, the active power of the flexible load at bus i in time t, and the active power of the sensitive load at bus i in time t, respectively; Q g,i,t , Q w,i,t , Q s,i,t , Q l,i,t , and Q m,i,t are the reactive power of the generator injected into bus i in time t, the reactive power of the wind turbine injected into bus i in time t, the eactive power of the SVC injected into bus i in time t, the reactive power of the flexible load at bus i in time t, and the reactive power of the sensitive load at bus i in time t, respectively.
In Equations (39) and (40), P ij,t and Q ij,t are the active and reactive power of the line between bus i and bus j at time t. θ ij,t is the voltage phase angle difference between bus i and bus j at time t [28,29]. In Equations (39) and (40), θ 2 ij,t is a non-linear variable, but it can be linearized by piecewise fitting linearity. As the piecewise fitting linearity was not the focus of this article, please refer to [30,31] for details. 0 ≤ P U g,i,t,r ≤ P U g,i,t , ∀i, t, r 0 ≤ P D g,i,t,r ≤ P D g,i,t , ∀i, t, r (42) P g,i,t,r = P g,i,t + P U g,i,t,r − P D g,i,t,r , ∀i, t, r where P U g,i,t,r and P D g,i,t,r are respectively the upregulation and downregulation of the power generator under the scene r; P g,i,t,r is the active output power of the power generator under the scene r; and Q gn g,i,t,r is the reactive output power of the power generator g under scene r in the time period t.

2.
Constraints for the redispatch of flexible load Constraints (46) and (47) are the upregulation and downregulation of the reserve constraint in the real-time stage of the flexible load. In Equation (48), P l,i,t,r is the active power output of the flexible load in the real-time stage under the scene r.
0 ≤ P D l,i,t,r ≤ P D l,i,t , ∀i, t, r (47) P l,i,t,r = P l,i,t + P U l,i,t,r − P D l,i,t,r , ∀i, t, r where P U l,i,t,r and P D l,i,t,r represent the up/down spinning reserve deployed by the flexible load of the bus i in the time period t under the scene r.
In addition to the above Constraints (41)−(48), the day-here real-time optimization model also met the operational constraints of the SVC, RES, and power balance, which had a similar form to (31) and (34)-(38).

Calculation Examples and Analysis
The proposed two-stage active and reactive power coordination control strategy considering a random ZIP linearized model was tested in a modified IEEE 33-bus system. The specific modifications are as follows: (1) bus 33 was connected to a 0.5 MW micro gas turbine; (2) bus 16 and bus 29 were connected to 0.9 MW and 0.8 MW wind turbines, respectively; (3) bus 6 and bus 18 were connected to 0.3 Mvar SVC controller devices, respectively; (4) the load at bus 14 and bus 26 were regarded as sensitive load buses; (5) the load at bus 7 and bus 24 are treated as the flexible loads with up to 30% of the load being able to be shed; and (6) the range of allowable voltage fluctuation was 0.930−1.070. The proposed model was coded in the GAMS environment using CPLEX 12.5. All experiments were implemented on a desktop computer with a processor clocked at 3.10 GHz and 8 GB of RAM.
This paper mainly focused on the loss cost of highly sensitive loads. By comparing four cases that contained different reactive power regulation means, the effectiveness of the strategy proposed in this paper was illustrated: (1) Case 1: NO adjustment measure; (2) Case 2: RES adjustment measure; (3) Case 3: RES + SVC adjustment measure; and (4) Case 4: OLTC + RES + SVC adjustment measure.
The purpose of the day-ahead optimization was to provide sufficient reactive power reserve to reduce the effects of the uncertainty factors to the sensitive loads. Figure 1 shows the day-ahead voltage profile at the 14 bus and 26 bus under four cases. Case 4 was the best voltage regulation method compared to the other three cases. It can well regulate the voltage of different buses at different times. In Figure 1a, the overall voltage fluctuation ratio of 14 bus was 1.57% in Case 4, 2.05% in Case 3, 2.71% in Case 2, and 3.26% in Case 1, respectively. Case 4 can significantly reduce the overall voltage fluctuation compared with the other three cases. The minimum voltage of Case 1 at 16 h was 0.949, which was very close to the lower limit of the voltage regulation range of 0.930. When the reactive power compensation devices of RES, SVC, and OLTC were added, the voltage of the load was significantly improved.
Especially after adding OLTC, the voltage increased most obviously, rising to 0.986 at 16 h. In Figure 1b, the lowest voltage of the 26 bus was 0.965 in Case 1 at 16 h. After the regulation of Case 4, the voltage increased to 0.996 at 16 h. The overall voltage fluctuation ratio of bus 26 was 1.42% in Case 4; 2.01% in Case 3; 2.30% in Case 2, and 2.51% in Case 1, respectively. The results showed that Case 4 had a good control effect on voltage at different times and overall voltage fluctuation of bus 26. In the case of day-ahead scheduling, Case 4 could stabilize the overall voltage stability of different highly sensitive loads and also effectively adjust the voltage at different times. In the day-here real-time redispatch, the voltage was further improved by the continuous reactive power compensation device. Figure 2 shows the real-time voltage changes under five different scenarios. In Figure 2a, the maximum voltage of Case 4 was 1.012 at 8 h and the average value of voltage at 8 h in the typical scene was close to the reference value, 1.00 p.u, which means that the day-ahead optimization produced by Case 4 had relatively high robustness. Furthermore, the value of bus 14 at 7 h in the day-ahead stage was 1.046. While the related reactive control participated in the day-here stage, the voltage was further improved to 1.01. Furthermore, the overall voltage deviation ratio of bus 14 was 3.29% in Case 1, 2.78% in Case 2, 1.63% in Case 3, and 0.80% in Case 4, respectively. From Figure 2a at 10.00-20.00, the effect of these reactive power regulating devices was not obvious when the system load was heavy. For this situation, the voltage can be improved by making the day-ahead OLTC scheduling plan.  In the day-here real-time redispatch, the voltage was further improved by the continuous reactive power compensation device. Figure 2 shows the real-time voltage changes under five different scenarios. In Figure 2a, the maximum voltage of Case 4 was 1.012 at 8 h and the average value of voltage at 8 h in the typical scene was close to the reference value, 1.00 p.u, which means that the day-ahead optimization produced by Case 4 had relatively high robustness. Furthermore, the value of bus 14 at 7 h in the day-ahead stage was 1.046. While the related reactive control participated in the day-here stage, the voltage was further improved to 1.01. Furthermore, the overall voltage deviation ratio of bus 14 was 3.29% in Case 1, 2.78% in Case 2, 1.63% in Case 3, and 0.80% in Case 4, respectively. From Figure 2a at 10.00-20.00, the effect of these reactive power regulating devices was not obvious when the system load was heavy. For this situation, the voltage can be improved by making the day-ahead OLTC scheduling plan. In the day-here real-time redispatch, the voltage was further improved by the continuous reactive power compensation device. Figure 2 shows the real-time voltage changes under five different scenarios. In Figure 2a, the maximum voltage of Case 4 was 1.012 at 8 h and the average value of voltage at 8 h in the typical scene was close to the reference value, 1.00 p.u, which means that the day-ahead optimization produced by Case 4 had relatively high robustness. Furthermore, the value of bus 14 at 7 h in the day-ahead stage was 1.046. While the related reactive control participated in the day-here stage, the voltage was further improved to 1.01. Furthermore, the overall voltage deviation ratio of bus 14 was 3.29% in Case 1, 2.78% in Case 2, 1.63% in Case 3, and 0.80% in Case 4, respectively. From Figure 2a at 10.00-20.00, the effect of these reactive power regulating devices was not obvious when the system load was heavy. For this situation, the voltage can be improved by making the day-ahead OLTC scheduling plan. Figure 3. shows the outlet voltage of the OLTC. Furthermore, the voltage profile of the ADN was obviously improved from Figure 4 after executing the day-ahead OLTC scheduling plan at 16h. A similar conclusion can be drawn from Figure 2b. In summary, Case 4 had the best voltage regulation effect in both the day-ahead and realtime stages.      (7), we can conclude that the load fluctuation caused by the voltage is greater than the reference value if the voltage is greater than 1.00 p.u and less than the reference otherwise. The symbol N and P represent the load change ratios greater than or less than the reference value, respectively. Within 24 h, the average load fluctuation rate of bus 14 was 1.07% in Case 1 and 0.54% in Case 4; the average load fluctuation rate of bus 26 was 2.45% in Case 1 and 0.31% in Case 4. It was obvious that the rate of load fluctuation in Case 1 was much higher compared with Case 4. The maximum load fluctuation ratio of bus 14 was 8.02% at 16 h in Case 1 and 2.80% at 19 h in Case 4 and the biggest voltage deviation of bus 14 was 0.051 p.u. at 16 h in Case 1 and 0.019 in Case 4, respectively. From the constraint, there was a direct relationship between the load fluctuation with the voltage deviation. If the voltage value is closer to 1.0 p.u., the loss ratio of the bus will be lower. For the voltage sensitive loads of bus 26, the maximum loss ratio was 4.08% at 14 h in Case 1 and 2.36% at 20 h in Case 4. The corresponding voltage deviation of bus 26 was 0.035 at 14 h in Case 1 and 0.023 at 20 h in Case 4. Compared with Case 1, the ratio of the maximum value of voltage deviation in Case 4 was greatly reduced. This demonstrates that Case 4 can not only significantly reduce the loss ratio of highly sensitive loads at different times, but also the overall load loss ratio.     (7), we can conclude that the load fluctuation caused by the voltage is greater than the reference value if the voltage is greater than 1.00 p.u and less than the reference otherwise. The symbol N and P represent the load change ratios greater than or less than the reference value, respectively. Within 24 h, the average load fluctuation rate of bus 14 was 1.07% in Case 1 and 0.54% in Case 4; the average load fluctuation rate of bus 26 was 2.45% in Case 1 and 0.31% in Case 4. It was obvious that the rate of load fluctuation in Case 1 was much higher compared with Case 4. The maximum load fluctuation ratio of bus 14 was 8.02% at 16 h in Case 1 and 2.80% at 19 h in Case 4 and the biggest voltage deviation of bus 14 was 0.051 p.u. at 16 h in Case 1 and 0.019 in Case 4, respectively. From the constraint, there was a direct relationship between the load fluctuation with the voltage deviation. If the voltage value is closer to 1.0 p.u., the loss ratio of the bus will be lower.    (7), we can conclude that the load fluctuation caused by the voltage is greater than the reference value if the voltage is greater than 1.00 p.u and less than the reference otherwise. The symbol N and P represent the load change ratios greater than or less than the reference value, respectively. Within 24 h, the average load fluctuation rate of bus 14 was 1.07% in Case 1 and 0.54% in Case 4; the average load fluctuation rate of bus 26 was 2.45% in Case 1 and 0.31% in Case 4. It was obvious that the rate of load fluctuation in Case 1 was much higher compared with Case 4. The maximum load fluctuation ratio of bus 14 was 8.02% at 16 h in Case 1 and 2.80% at 19 h in Case 4 and the biggest voltage deviation of bus 14 was 0.051 p.u. at 16 h in Case 1 and 0.019 in Case 4, respectively. From the constraint, there was a direct relationship between the load fluctuation with the voltage deviation. If the voltage value is closer to 1.0 p.u., the loss ratio of the bus will be lower. Fo the voltage sensitive loads of bus 26, the maximum loss ratio was 4.08% at 14 h in Case 1 and 2.36% at 20 h in Case 4. The corresponding voltage deviation of bus 26 was 0.035 at 14 h in Case 1 and 0.023 at 20 h in Case 4. Compared with Case 1, the ratio of the maximum value of voltage deviation in Case 4 was greatly reduced. This demonstrates that Case 4 can not only significantly reduce the loss ratio of highly sensitive loads at different times, but also the overall load loss ratio. Figure 6 shows the daily power network loss changes of the IEEE-33 bus system. The power network loss in Case 2 was the smallest, and the total power network loss in Case 1 was the largest. The power network loss in Case 4 was slightly higher than that in Case 2. The average power network loss in the four cases were 7.916 in Case 4, 5.396 in Case 3, 13.363 in Case 2, and 19.045 in Case 1. Table 1 shows the five operating costs in four cases including total cost (TC), day-ahead cost (DAC), real-time cost (RTC), highly sensitive load loss cost (HSLLC), and real-time power network loss cost (RTPNLC). In Case 4, the total cost, day-ahead, and real-time cost were the smallest. In Case 3, the power network loss cost was minimal. However, compared with Case 1, the loss cost of the highly sensitive load of Case 4 was significantly reduced. This paper mainly focused on the loss cost of highly sensitive loads. Although the network loss cost produced by Case 4 was not the smallest, compared with other cases, the loss cost of highly sensitive loads was the smallest, which demonstrates that the load fluctuations caused by voltage fluctuations were the smallest. From the point of view that the total cost of the system and the highly sensitive loss cost were the main issues, the adjustment method of Case 4 was the best. Therefore, the coordinated control strategy proposed in this paper can not only significantly reduce the cost of high-sensitivity load loss caused by voltage fluctuations, but also reduce the total power cost and the daily dispatch cost of an ADN.  Figure 6 shows the daily power network loss changes of the IEEE-33 bus system. The power network loss in Case 2 was the smallest, and the total power network loss in Case 1 was the largest. The power network loss in Case 4 was slightly higher than that in Case 2. The average power network loss in the four cases were 7.916 in Case 4, 5.396 in Case 3, 13.363 in Case 2, and 19.045 in Case 1.  Table 1 shows the five operating costs in four cases including total cost (TC), day-ahead cost (DAC), real-time cost (RTC), highly sensitive load loss cost (HSLLC), and real-time power network loss cost (RTPNLC). In Case 4, the total cost, day-ahead, and real-time cost were the smallest. In Case 3, the power network loss cost was minimal. However, compared with Case 1, the loss cost of the highly sensitive load of Case 4 was significantly reduced. This paper mainly focused on the loss cost of highly sensitive loads. Although the network loss cost produced by Case 4 was not the smallest, compared with other cases, the loss cost of highly sensitive loads was the smallest, which demonstrates that the load fluctuations caused by voltage fluctuations were the smallest. From the point of view that the total cost of the system and the highly sensitive loss cost were the main issues, the adjustment method of Case 4 was the best. Therefore, the coordinated control strategy proposed in this paper can not only significantly reduce the cost of high-sensitivity load loss caused by voltage fluctuations, but also reduce the total power cost and the daily dispatch cost of an ADN.   Figure 6 shows the daily power network loss changes of the IEEE-33 bus system. The power network loss in Case 2 was the smallest, and the total power network loss in Case 1 was the largest. The power network loss in Case 4 was slightly higher than that in Case 2. The average power network loss in the four cases were 7.916 in Case 4, 5.396 in Case 3, 13.363 in Case 2, and 19.045 in Case 1.  Table 1 shows the five operating costs in four cases including total cost (TC), day-ahead cost (DAC), real-time cost (RTC), highly sensitive load loss cost (HSLLC), and real-time power network loss cost (RTPNLC). In Case 4, the total cost, day-ahead, and real-time cost were the smallest. In Case 3, the power network loss cost was minimal. However, compared with Case 1, the loss cost of the highly sensitive load of Case 4 was significantly reduced. This paper mainly focused on the loss cost of highly sensitive loads. Although the network loss cost produced by Case 4 was not the smallest, compared with other cases, the loss cost of highly sensitive loads was the smallest, which demonstrates that the load fluctuations caused by voltage fluctuations were the smallest. From the point of view that the total cost of the system and the highly sensitive loss cost were the main issues, the adjustment method of Case 4 was the best. Therefore, the coordinated control strategy proposed in this paper can not only significantly reduce the cost of high-sensitivity load loss caused by voltage fluctuations, but also reduce the total power cost and the daily dispatch cost of an ADN.

Conclusions
In order to consider the static voltage characteristics, time characteristics, and uncertainties of highly sensitive loads in scheduling, a linearized voltage sensitive load model with time characteristics was proposed. Furthermore, this paper proposed a linearized two-stage active and reactive power coordination stochastic optimization model to reduce the impact of the uncertainty of load and wind power generation to the sensitive load of these voltages. In the day-ahead optimization, the best OLTC tap and active backup plan was obtained by the active and reactive coordinate strategy to reduce the large voltage fluctuation and providing the sufficient active power reserve to eliminate the unbalance power caused by the predict error. In the day-here real-time redispatch, the voltage was further improved by the continuous reactive power compensation device. Finally, our results on the modified IEEE 33-bus system demonstrated the proposed closed-loop coordinate strategy can effectively reduce the loss cost of highly sensitive loads and minimize the economic cost of an ADN by rationally utilizing the reactive power regulation capability of the equipment and the potential of flexible load up-and downregulation.
In this paper, the model of flexible load participating in grid dispatching did not consider the influence of environmental and user willingness factors. However, these factors have an impact on the maximum potential and the response uncertainty of flexible load participation in dispatching. Therefore, how to build a coordinated scheduling model with multiple time scales considering the response uncertainty of a flexible load and the loss cost of highly sensitive loads is our future research direction. Furthermore, the scheduling cost of making a scheduling plan to meet all the uncertainties is very expensive. Therefore, how to introduce the chance constraints in a multiple time scales scheduling model is also our future research.