1. Introduction
Variable renewable energy sources (VRESs), characterized by uncertainty, are expected to play a crucial role in the decarbonization of the energy sector. The integration of these sources into the power grid requires increased flexibility to ensure reliable, efficient, and economical operation. The flexibility of the power system is critical in managing the variability and uncertainty inherent to renewable energy generation, affecting both short-term operational strategies and long-term planning. Several studies highlight the importance of developing flexible systems to accommodate the rapid fluctuations in power supply and demand, mitigate the challenges posed by increased renewable penetration, and maintain system stability and reliability [
1,
2]. Note here that in this paper, although flexibility can be achieved through various means, including advanced forecasting, improved ramping capabilities, and the strategic deployment of energy storage and responsive load management techniques, we are primarily focused on load and/or generation flexibility facilitated by the technical capability to adjust grid users’ power consumption or generation upon request. While rather abstract, this definition is sufficiently specific, remaining agnostic to technical details of
how flexibility is provided or if there are certain temporal or technical constraints regarding the available flexibility. In this paper, we take the perspective of a distribution system operator, while assuming a flexibility provider can always provide an amount of flexibility that has been agreed on previously with the system operator.
The introduction of VRESs into the distribution grid represents a significant operational paradigm shift for the distribution system operators who face significant challenges due to the lack of measurements from the power grid [
3]. The lack of measurement data complicates accurate state estimation, making it difficult to verify real-time network constraints and leading to overly cautious operational strategies [
4]. With limited network observability, the operators struggle to monitor and control VRESs effectively, which can result in inefficient grid management, increased risk of instability, and sub-optimal decision making. The scarcity of precise, timely data also forces system operators to rely heavily on pseudo-measurements and forecasts, which can introduce significant errors and uncertainties into the state estimation process [
5,
6].
Using flexibility of load and/or generation to mitigate voltage or current constraint violations has been extensively researched [
7]. However, most of these approaches make the assumption of full grid observability, often using exact data for demand and generation in their numerical examples. More recently, the lack of measurements and how this lack can reflect on system operations is being recognized and addressed. In [
8], the authors address the challenge of maintaining acceptable voltage levels in distribution grids with increasing solar PV penetration, focusing on optimizing the schedule for transformers and capacitor banks and minimizing the use of responsive inverter-interfaced resources for real-time conditions. It presents a convex optimal power flow (OPF) formulation to minimize the voltage deviations and integrates discrete mechanical devices into the optimization problem. However, the convex inner approximation developed in the paper is unlikely to be sufficiently robust to varying network topology and is not applicable to more realistic grids. In [
9], the authors consider the robust AC OPF problem, which minimizes the generation cost while requiring a certain level of system security in the presence of uncertainty. They extend a previously developed convex restriction to a robust convex restriction: a convex inner approximation of the non-convex feasible region of the AC OPF problem that accounts for uncertainty in power injections. However, the proposed method is better suited for transmission networks and, given the nature of convex restriction, could yield overly conservative results in distribution networks. In [
10], a two-stage min–max–min robust energy management model considering the correlation constraints and the uncertainty of solar energy is proposed. However, the applicability of the paper to grid management is rather questionable, as certain parameters regarding the treatment of uncertainties are not readily available nor intuitive for system operators (e.g., distribution of uncertainties). Ref. [
11] introduces a holistic approach to optimal operational planning and operational management of active distribution systems under uncertainties to avoid line congestion and voltage limit violations using a two-stage stochastic programming model based on weighted scenarios. While comprehensive, the method is computationally demanding.
Furthermore, the issue of
when to use the available flexibility, given the lack of measurements, remains unexplored. In [
12], while not considering flexibility, the authors use convex relaxation to certify grid voltage constrain satisfaction. However, their approach results in overly conservative voltage estimation, precisely because of the use of relaxations. In [
13], the authors propose a method for constraining the control actions of third-party aggregators to ensure safe operation of the distribution network. While somewhat conceptually similar, our paper differs by both the nature of the problems it addresses and the class of underlying optimization problems. A bi-level optimization-based framework to determine the aggregate power flexibility that can be obtained from an unbalanced distribution grid while ensuring that there is no solution that leads to grid constraint violations is proposed in [
14]. While somewhat conceptually similar to our paper, this work is primarily focused on determining the aggregate flexibility of a distribution network, rather than using it to mitigate problems in the distribution grid itself. In a similar vein, authors in [
15] formulate a batch reinforcement learning-based demand response approach to prevent distribution network constraint violations in unknown grids. While accounting for uncertainties in the network, the paper’s concern is that of a demand response for frequency regulation.
While developing our optimization framework, we selected a single-objective mixed-integer linear programming (MILP) model due to its computational efficiency, robustness, and practicality. MILP models allow for efficient problem solving using commercial solvers, which is essential given the significant computational time required for large network regression problems. The linear nature of MILP constraints and objective functions simplifies integration with existing grid management systems and ensures scalability. This reliability and ease of implementation are critical for real-world applications, where more complex models may pose significant deployment challenges. Additionally, while learning-based optimization frameworks and multi-objective mixed-integer nonlinear programming (MINLP) frameworks offer potential benefits, they require extensive and high-quality datasets for training and validation—resources that are currently limited in power network flexibility and injection scenarios. A single-objective MILP provides clear, interpretable solutions crucial for decision making in power grid management, facilitating better understanding and trust among operators and stakeholders. Our future research will focus on integrating advanced optimization techniques and data-driven methods as more comprehensive datasets become available and computational capabilities improve. This approach ensures our method remains efficient, robust, and practical for real-world power network optimization.
In order to solve the problem of mitigating grid constraints in an unobservant grid, authors often rely on overly conservative or computationally expensive methods that require certain techniques to circumvent the issue of non-convexity of the AC power flow equations or design methods that accommodate a small subset of distribution grids. This paper tries to address this gap by designing a tractable framework that is robust with regard to demand and generation uncertainties; conservative, as it relies on conservative linear approximations, therefore providing the distribution system operator with sufficiently tight lower and upper bounds on decision variables by construction; data-efficient, in a sense that it only requires knowledge of ranges of power injections, usually readily available to the grid operator; temporally agnostic, in that it is able to encompass an arbitrarily large operational time span, given its nature of dealing with ranges of power injection; and is practical, from an engineering perspective, by effectively avoiding complex mathematical transformations.
The objective of this paper is twofold. Firstly, we want to determine the amount and the location of flexibility that is needed to ensure no voltage constraints are violated in a distribution network with limited observability. Here, we assume that only the range of possible power injections and the range of voltage at the reference node are known. Secondly, we want to determine, given a specific location, a voltage measurement threshold that can serve as a proxy for the voltage condition of the overall network which, if crossed, would mean at least one voltage constraint violation is occurring somewhere in the network. This value then serves as the trigger for activating flexibility.
In summary, the main contributions of this paper are as follows:
- 1.
This paper leverages recently proposed conservative linear approximations (CLAs) to formulate tractable, optimization-based problems that are able to mitigate voltage problems in a distribution network with low observability. This paper extends the proposed CLAs to include changes in voltage at the reference node.
- 2.
To the best of the authors’ knowledge, this is the first paper that establishes a concept of flexibility provision triggered by a pre-calculated voltage threshold in a network with limited observability and establishes a tractable framework to derive and validate the solutions within this scheme.
- 3.
In order to find the amount of flexibility and voltage that triggers it, two novel, optimization-based, iterative algorithms are proposed.
2. Description of Problems
Let us consider a simplified model of a medium-voltage, radial, distribution network from [
16]. We assume no measurements are available in the network; however, based on the billing data, typical customer load profiles, and DG historical data, we can assume bounds on node power injections to vary within 50% to 150% of real and reactive power nominal value, as defined by the test case. Voltage at a reference node varies between 0.97 and 1.03 p.u. A reference node is located on the secondary side of a substation that feeds the network. Reference voltage magnitude is independent of power injections and is usually measured and set by the tap changer in a substation.
With these assumptions, we run power flow simulations by varying power injections and reference voltage values within the aforementioned limits. As a result, we obtain the magnitudes of the minimum achievable voltages for each node in the network, as shown in
Figure 1. We see that voltages in some nodes in the network can be as low as 0.87 p.u.
Traditionally, the distribution system operator must comply with a 0.90 p.u. at the low-voltage customers’ grid connection points. Since we are dealing with a medium-voltage network, maintaining a 0.92 p.u. at this voltage level allows for an additional voltage drop of 0.02 p.u. at a lower voltage level. In order to maintain the voltage, and in the absence of other voltage control methods (such as on-load tap changing equipment or capacitor banks), the distribution system operator can turn to flexibility providers. These are the customers that have the technical ability to provide ancillary services to the system operator by adjusting their power usage at the point of common coupling. When procuring flexibility, the distribution system operator needs to know in advance where, how much, and when the flexibility should be procured in order not to violate the grid’s voltage limits. To answer the aforementioned question of where and how much, we define the Flexibility Provision Problem as the optimization problem of determining the optimal location and the amount of flexibility that would guarantee no voltage limit violation, irrespective of the power injections in the rest of the distribution network.
Since the operator has to deal with low distribution network observability, it is important to have at least one voltage measurement that can serve as a proxy for overall voltage awareness in the network; this is the key to resolving the issue of ‘when’ to use flexibility. The placement of this voltage measurement, which we refer to as critical voltage, ensures that when the measured value of the critical voltage is below the predefined threshold value, the flexibility, determined previously by solving the Flexibility Provision Problem, will be activated by the distribution system operator. The Critical Voltage Threshold Problem is an optimization problem for determining the threshold value of critical voltage that guarantees no lower (or upper) voltage violation can occur in the entire network if the measured voltage magnitude is above (below) the pre-calculated threshold value. The flexibility provision, which is calculated by solving the Flexibility Provision Problem, should be seen as an upper limit of flexibility that is needed to provide robustness against a realization of the worst-case scenario. Knowing and setting the critical voltage threshold could enable autonomous flexibility control, where critical voltage provides a feedback for power injection control without the need for operator intervention.
4. Flexibility Provision Problem Formulation
In order to find the location and the amount of flexibility provision, we propose a two-stage, optimization-based method that loops between two levels: an upper level, which finds locations and the appropriate amount of flexibility, and a lower level, which finds a combination of power injections that maximizes the number of voltage bound violations. A detailed formulation of each level is described in the following subsections.
4.1. Lower-Level Problem
Given the range of nodal power injections
, and the voltage range at the reference node
, this optimization problem maximizes the number of lower or upper voltage limit violations. Depending on which voltage bound we want to violate (upper or lower), appropriate constraints in Equations (
3) and (
7) are selected. The full optimization problem is structured as follows:
where
is an
-sized vector of squared voltage variables, while parameters
,
,
, and
are CLA linear coefficients, as described in
Appendix A.
denotes the number of nodes in the network, while
denotes a set of nodes where flexibility is to be activated. The
index denotes the results from the upper-level problem.
are vectors of size
, while
is a vector of nodal lower voltage limits.
If one or several variables in vector are set to one, the corresponding variables in vector are set to zero, causing v to violate the lower or upper voltage bound. Set of nodes with activated flexibility——and the corresponding amount of flexibility——are the results of the upper-level problem. The decision variables of the optimization problem that are propagated to the upper-level problem are the nodal power injections and reference voltage, both denoted by the index (lower level). The optimization procedure is designed to find the set of nodal power injections that represent the worst-case scenario with regard to the voltage bound violation.
Note here that the optimization problem is structured to find the maximum number of voltage violations. A viable alternative to this would be to minimize (or maximize) each node voltage and then to verify if any voltage violates the limits. In contrast to our method, this would require substantially more computational time on larger networks. Another viable option is to minimize (or maximize) the sum of all nodal voltages. While computationally more efficient, this method is generally less conservative, which makes it less preferred by the distribution system operator.
4.2. Upper-Level Problem
Given the bounds on active and reactive power flexibility
and
, as well as the power injections determined by solving the lower-level problem,
, and
, we formulate the upper-level optimization problem as follows:
where
, while index
denotes parameters from the lower-level problem.
The reasoning behind the previous problem is now inverted; instead of trying to violate voltage constraints, we enforce the satisfaction of voltage constraints in Equation (
11) by ‘procuring’ the flexibility from the providers. This is accomplished by setting the variables in vector
s to one, meaning that the flexibility at a specific node should be used. If a variable in the vector is set to zero, the flexibility is not used, and power injections from lower-level problem are used instead. The decision variables of the optimization problem are nodal power injections for nodes where the flexibility is procured.
Note here that the optimization problem is not structured to find the minimum amount of flexibility; rather, it is trying to limit the number of flexibility providers, a feature which could provide sub-optimal results with regard to the total amount of flexibility needed to mitigate potential voltage issues. While sub-optimal, this design choice is intentional, as the distribution system operators prefer to minimize their interactions to as few entities as possible. Furthermore, once an operator knows when to activate flexibility, they do so in a sequential manner, until the critical voltage is restored to a predefined value. This means that the flexibility calculated here will always be the upper bound to what is actually used by the distribution system operator.
4.3. Flexibility Provision Algorithm
We assume that the CLA parameters have already been determined as outlined in
Appendix A. The key steps of the Algorithm 1 that connects the lower and upper optimization problems are as follows:
Algorithm 1 Flexibility Provision Algorithm |
- 1:
Input: , , ,,,, , ,, - 2:
while true do - 3:
the lower-level problems ( 1)–( 8) - 4:
if then save the values of and . - 5:
else - 6:
break - 7:
end if - 8:
Solve the upper-level problems ( 10)–( 13). - 9:
Update with indices of s, where . Save , - 10:
end while
|
The algorithm is initialized with an empty set of
, meaning no flexibility is provided. The algorithm then iterates between the lower- and upper-level problems, until the objective of a lower-level problem is zero, meaning that no further voltage bound violations can occur in the network. The algorithm is run once for each (lower and upper) voltage bound violation, with the only difference being the choice of constraints in Equation (
2) and CLA parameters used in Equations (
9) and (
14).
The convergence of the algorithm is not guaranteed. Given the topology, power, and voltage limits, as well as current operating point of the network, it is not unrealistic that there is not enough flexibility in the network that would make the network robust for the entire range of possible power injections and reference voltage variations. However, the algorithm can, even in those cases, provide insight into how this robustness can be achieved. For example, flexibility bounds can be relaxed to ensure convergence.
7. Discussion
If we disregard YALMIP overhead, the computational time needed for solving the regression problem of (
A3)–(
A5) for larger networks can be substantial; hence, using a commercial MILP solver is recommended in this case. This calculation can be completed offline, so we do not consider it a major drawback. However, further research is needed to reduce the number of samples used in the regression problem (M), while maintaining the accuracy of the CLAs.
Realistic data regarding flexibility availability and ranges of power injections are needed to conduct further research into the potential downsides of the proposed method.
In order to avoid overly conservative results, a more complex construction of an uncertainty set, which is specially tailored to operators intuition and grid data availability, such as in [
22,
23], is needed and is an interesting avenue for future research.
We emphasize the need for further research into selecting the appropriate location for voltage measurement placements that can serve as a proxy for voltage conditions of the entire network. The problem with arbitrarily choosing a location is that a voltage measurement is not guaranteed to be sensitive to changes in power injections stemming from the use of flexibility. This could lead to a situation where the activation of flexibility mitigates voltage problems in a part of the network, but this is not reflected in voltage measurements that are monitored, leading the operator to unnecessarily activate additional flexibility.
8. Conclusions
The provision of flexibility will likely play a key role in future power grids. From a technical perspective, operators are currently unable to exploit the flexibility potential of the consumer, mostly because of the lack of tools that incorporate system uncertainties, some of which stem from the lack of medium- and low-voltage measurements. To address the issue, this paper proposes a method to determine the location and the amount of flexibility needed to mitigate potential voltage problems, given limited network observability. The method leverages recently proposed conservative linear approximation to establish tractable, mixed-integer linear optimization problems that are used to determine the required flexibility. Furthermore, this paper proposes a framework that enables the operator to know when to use the required flexibility by determining the critical voltage threshold value. As shown on 33 and 124 bus distribution system test cases, the iterative algorithms proposed are able to determine the amount of flexibility and the voltage threshold at which this flexibility is to be activated efficiently and accurately, requiring only the knowledge of the ranges of nodal power injections.
To enhance the robustness of the proposed method, future research should focus on optimizing the sample size used in the regression problem while ensuring the accuracy of the CLAs, as well as on developing advanced techniques for constructing uncertainty sets tailored to grid operators’ intuition and data availability. Additionally, further research is required to identify optimal voltage measurement locations that accurately reflect network-wide voltage conditions, thereby preventing the unnecessary activation of flexibility and ensuring efficient grid management.