In this section, we firstly introduce the proposed multi-timescale reactive power control framework for distribution system, then formulate the optimization model of day-ahead and intraday timescale, and lastly we propose the event-triggered real-time voltage zoning control strategy based on voltage sensitivity to solve the short-term fluctuation and sudden overrun of voltage.
3.1. Multi-Timescale Reactive Power Control Scheme for Distribution Networks
MPC is a fuzzy predictive control theory used in the optimization and regulation of dynamic systems [
31,
32]. The basic idea is to predict the future behavior of the system by combining the measured data at the current moment in order to optimize the control action in a future period of time. Considering the uncertainties that may exist in the input data, this control method optimizes its own prediction results by rolling iterations and feedback of the prediction results over multiple single time intervals. At each time interval the MPC builds a new optimization model based on the current operating parameters, so the set iteration length affects the MPC output results, and if the iteration length is shortened, the data sampling rate needs to be increased at the same time as the MPC accuracy is improved. MPC acts on the stochastic system by means of unfixed feedback control values. Through intra-day rolling optimization, the system is able to automatically generate an optimization model that gradually approaches the optimal control strategy by continuously reducing the error.
In the operation and control of distribution networks, the uncertainty of renewable energy power outputs means that the day-ahead optimal model may not achieve system optimization during intra-day operation. Therefore, the MPC method can be adopted to develop more effective strategies for multi-source interactive distribution networks. Considering the uncertainty of distributed new energy resources, this approach calculates the active and reactive regulation strategies of these resources in the distribution network and compares the roles of different distributed resources in the voltage regulation process. By utilizing the principles of model predictive control to coordinate the varying voltage regulation capabilities of each distributed resource, the method effectively addresses voltage overrun issues in the distribution network. With its multi-step optimization and prediction capabilities, MPC can determine a multi-scale voltage optimization model for distribution networks in a multi-source interactive environment.
The multi-timescale reactive power control framework based on MPC method for distribution system is shown in
Figure 4, which includes day-ahead prediction optimization, intra-day rolling optimization, and real-time voltage control. Intra-day rolling optimization uses the active and reactive power data of distributed resources obtained from day-ahead predictions to continuously perform feedback iterations. The day-ahead optimization layer aims to reduce network losses and maintain voltage stability, with a time scale of 1 hour, coordinated by the distributed resource voltage regulation model in
Section 2.2. The intra-day rolling optimization layer operates on a 15-minute time scale, focusing more on voltage stability control and the uncertainty of distributed resource output. Real-time voltage control, building on the previous layers, operates on a 5-minute time scale to keep voltage limits within the optimal operating range.
(1) Objective function
In order to maintain the stable operation of the distribution network, the day-ahead optimization model takes the optimal operation of the system as the objective, therefore the optimization objective function is modeled as minimizing the branch loss and voltage deviation of the distribution network, which is formulated as follows
where
denotes the line loss,
denotes the voltage of node
i at time slot
t,
denotes the reference voltage, which is generally set to 1,
and
denote the node set and branch set, respectively, and
and
are the weight coefficients.
(2) Constraints
There are the following constraints that need to be satisfied during the optimization:
1) Power flow constraints
For each node
i,
and
denote the active and reactive power injected into the node at the time slot
t ; for each line
,
and
denote the conductance and conductance of the line respectively,
denotes the active power flowing through the line at time slot
t, and
d enotes the difference in the phase between the end and the beginning of the line.
is the set of all nodes connected with node
i. In this paper, the polar coordinate form of the power flow model is used:
where
and
denote the PV and WT output of node
i at time slot
t,
and
denote the charging and discharging power of node
i at time slot
t,
denotes the user load of node
i at time slot
t,
and
denote the upper and lower voltage limits of the distribution node
i,
and
denote the current on the branch
and the maximum value allowed.
3) PV operation constraints
where
is the predicted active power of PV node
m,
is the capacity of the grid-connected PV converter, and
is the reactive power regulation strategy, with upper and lower bounds
and
, respectively:
4) WT operation constraints
where
is the active power generation at the WT node
n;
is the capacity of the grid-connected converter for WT, and
is the reactive power regulation strategy, and according to models given by Eqs. (3)-(8), the reactive power of WT is constrained by:
5) ESS operation constraints
where
and
are binary variables representing the charging and discharging state of the
k-th ES unit respectively,
and
are the charging and discharging power of the ES unit
k at the time slot
t,
and
are the upper limits of the charging and discharging power,
and
are the active and reactive power issued or absorbed by the
k-th ES at time slot
t,
is the inverter capacity, and
denotes the maximum capacity of the ES battery,
and
denote the minimum and maximum state of charge, respectively.
3.3. Intra-Day Rolling Optimization Model
Based on the day-ahead optimization results, the intraday rolling optimization model shortens the forecast time scale to cope with the fluctuation of renewable energy output. Therefore, the renewable energy outputs are forecasted and updated with a time resolution of 15 minutes, and the day-ahead optimization dispatch value is used as a reference value to calculate the regulation strategy of DERs in the control time horizon (4 hours) in a rolling manner.
When performing the intraday rolling optimization, the day-ahead operation plan of distributed resources must be followed, so the intraday objective function is to minimize the voltage deviation value and the penalty of distributed resource output regulation, which is formulated by:
where
denotes the voltage of node
i at time slot
t of intraday stage,
and
denote the active and reactive output correction compared with day-ahead optimization results of DERs,
and
are the corresponding penalty coefficients, and
and
are the weight coefficients.
(2) Constraints
In addition to power flow constraints, the following operation constraints of DERs also need to be satisfied:
where
and
are the active and reactive output correction compared with day-ahead optimization results of PV, similarly,
and
are the active and reactive output correction of WT, and
,
and
are the active and reactive output correction of ES.
3.4. Real Time Control Model
For voltage fluctuations and over-limits in a short period of time in the real-time stage, the optimization method is not used to solve the regulation instructions of each distributed resource, in this paper the ‘event-triggered’ voltage partition control strategy is directly used to obtain the control signal. And in order to further improve the rapidity and effectiveness of the distributed resource response to the voltage over-limit node, the distributed resource regulation priority is determined based on the active power-voltage sensitivity and reactive power-voltage sensitivity.
For each node in the system, the power-voltage sensitivity is used to measure the effect that a power change at one node has on the voltage at another node, as shown in equation (27):
where
,
are the active-voltage sensitivity and reactive-voltage sensitivity between nodes
i and
j respectively, representing the change in voltage of node
i when node
j changes unit active power and reactive power,
is the voltage change amount of the node
i,
and
are the active and reactive power change amount of the node
at time slot
t.
and
can be calculated from the modified equation of power flow calculation.
The distribution system can allocate the amount of active and reactive power to be regulated based on the power-voltage sensitivity between different distributed resource nodes to nodes of voltage violations. Nodes with higher sensitivity are better at voltage regulation. During the regulation process, nodes with voltage violations are adjusted according to the descending order of power-voltage sensitivity, which means DER nodes with high sensitivity being regulated first. Once the capacity of high sensitivity nodes is exhausted, regulation proceeds to lower sensitivity nodes. In addition, the regulation process also follows the voltage regulation strategy of “reactive power followed by active power, renewable power followed by energy storage system”, the specific strategy is shown in the figure below:
In the figure, the horizontal coordinate represents the current node voltage value,
and
represents the upper and lower limits of system voltage,
and
represent the upper and lower limits of the optimal voltage interval. The vertical coordinate represents the amount of active and reactive power that needs to be regulated by the system under the current voltage level. When the node voltage is not in the optimal voltage range, a “node with voltage violations” event is triggered, and the vertical coordinate corresponding to the current voltage level is the amount of active and reactive power to be regulated by multiple DERs. From
Figure 5, we can observe that the control area is divided into several different zones according to voltage intervals, and each zone corresponds to a corresponding distributed resource adjustment strategy, which is shown in
Figure 6:
For each voltage interval in the voltage partitioning control strategy of
Figure 5 and
Figure 6, the order of regulation of distributed resources is determined based on the active-voltage sensitivity or reactive-voltage sensitivity calculated in Eq. (39), and when the regulation capacity of all distributed resources is utilized, then it is transferred to the next interval, and so on.