Next Article in Journal
Assessment of Boil-Off Losses and Their Cost Implication during Liquid Hydrogen Tank Filling with and without Precooling
Previous Article in Journal
Thermoecological Cost Analysis of Hydrothermal Carbonization for Valorization of Under-Sieve Fraction from Municipal Solid Wastes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Voltage-Triggered Flexibility Provision in a Distribution Network with Limited Observability

1
Croatian Distribution System Operator HEP ODS, Ltd., 48000 Koprivnica, Croatia
2
Faculty of Electrical Engineering and Computing, University of Zagreb, 10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Energies 2024, 17(16), 4091; https://doi.org/10.3390/en17164091
Submission received: 1 July 2024 / Revised: 8 August 2024 / Accepted: 13 August 2024 / Published: 17 August 2024

Abstract

:
Variable renewable energy sources (VRESs) are essential for decarbonizing the energy sector, but they introduce significant uncertainty into power grids. This uncertainty necessitates increased flexibility to ensure reliable and efficient grid operations, impacting both short-term strategies and long-term planning. Effective management of VRESs is particularly challenging for distribution system operators (DSOs) due to limited grid measurement and other data, complicating state estimation. This paper proposes a tractable framework that operates under low-observability conditions. The framework uses conservative linear approximations (CLAs) to manage grid constraints efficiently, requiring only the ranges of power injections typically available to operators. The objectives are twofold: first, to determine the amount and location of flexibility needed to prevent voltage violations and, second, to establish voltage measurement thresholds to trigger flexibility activation. Key contributions include the formulation of optimization problems to mitigate voltage issues, the introduction of flexibility provision triggered by voltage thresholds, and novel algorithms for determining flexibility and trigger points. The framework’s efficacy is demonstrated on IEEE 33-bus and UK 124-bus test systems, showing it can effectively mitigate grid voltage problems despite limited data.

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.

3. Conservative Linear Approximations

The solution to the aforementioned optimization problems leverages the recently proposed conservative linear approximations (CLAs) of the power flow equations. The detailed discussion on calculating the CLAs and the associated parameters can be found in Appendix A. Note here that we build upon a previous work and extend the CLAs to account for changes in the reference voltage as these can not be ignored if there is no on-load tap changer at the substation.

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 [ P m i n P m a x ] , [ Q m i n Q m a x ] , and the voltage range at the reference node [ V m a x r e f V m i n r e f ] , 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:
max P l l , Q l l , v , s 1 , s 2 , V r e f s 1 = 1 N b s 1
s . t . s 2 v m i n 2 v v m i n 2 + ( 1 s 1 ) or s 1 v m a x 2 v v m a x 2 + ( 1 s 2 )
s 1 + s 2 = 1
Q m i n Q l l , i Q m a x i N b { f l e x }
P m i n P l l , i P m a x i N b { f l e x }
V m i n r e f V l l r e f V m a x r e f
P l l , i = P u p , i i { f l e x }
Q l l , i = Q u p , i i { f l e x }
where
v = a ̲ 0 + a ̲ 1 V l l r e f P l l Q l l or v = a ¯ 0 + a ¯ 1 V l l r e f P l l Q l l
is an N b -sized vector of squared voltage variables, while parameters a ̲ 0 , a ̲ 1 T , a ¯ 0 , and a ¯ 1 T are CLA linear coefficients, as described in Appendix A. N b denotes the number of nodes in the network, while f l e x denotes a set of nodes where flexibility is to be activated. The u p index denotes the results from the upper-level problem. s 1 , s 2 { 0 , 1 } are vectors of size N b , while v m i n is a vector of nodal lower voltage limits.
If one or several variables in vector s 1 are set to one, the corresponding variables in vector s 2 are set to zero, causing v to violate the lower V m i n or upper V m a x voltage bound. Set of nodes with activated flexibility— f l e x —and the corresponding amount of flexibility— P u p , Q u p —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 l l (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 [ P m i n f l e x P m a x f l e x ] and [ Q m i n f l e x Q m a x f l e x ] , as well as the power injections determined by solving the lower-level problem, P l l , Q l l , and V l l r e f , we formulate the upper-level optimization problem as follows:
min P , Q , v s = 1 N b s
s . t . v m i n 2 v v m a x 2
Q l l + s ( Q m i n f l e x Q l l ) Q u p Q l l + s ( Q m a x f l e x Q l l )
P l l + s ( P m i n f l e x P l l ) P u p P l l + s ( P m a x f l e x P l l )
where
v = a ̲ 0 + a ̲ 1 V l l r e f P u p Q u p or v = a ¯ 0 + a ¯ 1 V l l r e f P u p Q u p ,
s { 0 , 1 } , while index l l 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:  N b , f l e x Ø , a ¯ i , 0 , a ¯ i , 1 , a ̲ i , 0 , a ̲ i , 1 , P m i n f l e x , P m a x f l e x , Q m i n f l e x , Q m a x f l e x , P m i n , P m a x , Q m i n , Q m a x , V m a x r e f , V m i n r e f
  2:
while true do
  3:
      the lower-level problems (1)–(8)
  4:
      if  s > = 1  then save the values of P l l , Q l l and V l l r e f .
  5:
      else
  6:
            break
  7:
      end if
  8:
      Solve the upper-level problems (10)–(13).
  9:
      Update { f l e x } with indices of s, where s 0 . Save P u p , Q u p
10:
end while
The algorithm is initialized with an empty set of f l e x , 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.

5. Critical Voltage Threshold Problem Formulation

5.1. Lower-Bound Voltage Threshold Calculation

Here, we formulate an optimization problem that calculates the threshold for the lower bound of critical voltage. The problem is structured so as to ensure that no violation of the lower voltage limit can occur in the network if critical voltage remains above this threshold level.
max P , Q , v , s 1 , s 2 , V r e f , v t h r s 1 = 1 N b s 1
s . t . s 2 v m i n 2 v v m i n 2 + ( 1 s 1 )
s 1 + s 2 = 1
Q m i n Q Q m a x
P m i n P P m a x
V m i n r e f V r e f V m a x r e f
v t h r = ( V t h r i n i t ( 1 + η ) ) 2
where
v = a ̲ 0 + a ̲ 1 V r e f P Q
t h r is an index of the node for which the threshold is calculated, V i n i t denotes an initial voltage, and η is a parameter which increases the initial voltage value with each iteration. Given the arbitrary location of the node ( t h r ) and the initial voltage V t h r i n i t , the program increases the critical voltage threshold value by η in each iteration while still trying to violate lower voltage limits ( v m i n ).

5.2. Upper-Bound Voltage Threshold Calculation

The threshold for an upper voltage bound of critical voltage is calculated in a similar fashion, as follows:
max P , Q , v , s 1 , s 2 , V r e f , v t h r s 1 = 1 N b s 1
s . t . s 1 v m a x 2 v v m a x 2 + ( 1 s 2 )
s 1 + s 2 = 1
Q m i n Q Q m a x
P m i n P P m a x
V m i n r e f V r e f V m a x r e f
v t h r = ( V t h r i n i t ( 1 η ) ) 2
where
v = a ¯ 0 + a ¯ 1 V r e f P Q .
Here, we only change Equations (25), (29) and (30) to account for the upper instead of the lower threshold value.

5.3. Initial Voltage Calculation

The initial voltage ( V i n i t ) is calculated by solving one instance of the lower-level problem formulated in Section 4.1, but with the following modifications. The objective function is minimized instead of maximized, and an additional constraint (33) is added to the original problem. The optimization problem now calculates a set of network voltages where only one is below (above) the lower (upper) voltage limit.
min P l l , Q l l , v , s 1 , s 2 , V r e f s 1 = 1 N b s 1
s . t . s 2 v m i n 2 v v m i n 2 + ( 1 s 1 ) or s 1 v m a x 2 v v m a x 2 + ( 1 s 2 )
s 1 = 1 N b s 1 1
s 1 + s 2 = 1
Q m i n Q Q m a x
P m i n P P m a x
V m i n r e f V r e f V m a x r e f

5.4. Critical Voltage Threshold Algorithm

We assume that the CLA parameters have already been determined, as outlined in Appendix A. The key steps of the Algorithm 2 are as follows:
Algorithm 2 Critical Voltage Threshold Algorithm
  1:
Input:  N b , η , t h r { N b } a ¯ i , 0 , a ¯ i , 1 , a ̲ i , 0 , a ̲ i , 1 , P m i n , P m a x , Q m i n , Q m a x , V m a x r e f , V m i n r e f
  2:
Initialize: Solve opt. problem from (31)–(37) for v i , i N b
  3:
Arbitrarily choose t h r and set: V t h r i n i t v t h r
  4:
while true do
  5:
      Find the lower-bound voltage threshold by solving (15)–(21)
  6:
      if  s 1 > = 1 then increase η
  7:
      else
  8:
            break
  9:
      end if
10:
end while
11:
Output:  v t h r l o w v t h r
12:
while true do
13:
      Find the upper-bound voltage threshold by solving (23)–(29)
14:
      if  s 1 > = 1 then increase η
15:
      else
16:
            break
17:
      end if
18:
end while
19:
Output:  v t h r u p p v t h r
Note here that the choice of voltage measurement location for which the threshold is being calculated is arbitrary. This allows the operator to calculate the voltage threshold for a location in the network where the measurement can be placed and monitored for a threshold value violation. Although there is no guarantee that the solution to the problem, given the predefined location, can be found, we did not encounter any problems in our test cases. However, further insights into whether this is broadly the case remains an open and interesting direction for further research.

6. Numerical Results

6.1. Test System and Software

The optimization problems and algorithms are coded in MATLAB version 9.5.0 (R2018b) using theYALMIP toolbox [17] version 20230622 and SCIP [18] version 4.0.1, as a (mixed integer) linear programming solver. Test system is a modest PC with 16 GB of RAM. Power flows are calculated using MATPOWER [19] version 7.1.

6.2. IEEE 33 Bus Test Case Analysis

We illustrate the proposed method on a well-known IEEE 33 bus test case from [16]. It is a 12.66 kV medium-voltage-level radial distribution network, which encompasses 33 buses and 32 lines. Node 1 is a reference node with voltage ranging from 0.97 to 1.03 p.u. We assume that every other node in the network (2−33) has independent active/reactive power injections within the range of −100% to 100% of nominal power. The range of possible power injections is shown in Figure 2. The range assumes some form of distributed generation connected in parallel to the load at every node in the network. Lower and upper voltage limits are set to 0.96 and 1.04 p.u., respectively, for all network nodes. Given the aforementioned assumptions, the maximum and minimum achievable node voltages are shown in the first graph in Figure 2, calculated by running ten thousand power flow simulations. We see that, given a certain combination of reference voltage and power injections, both lower and upper voltage limits can be violated.

6.3. Flexibility Provision

Here, we want to find the amount of flexibility per node that can ensure no voltage violations can occur. Firstly, we calculate CLA parameters as described in Appendix A. We generate one thousand power flow samples ( M = 1000 ), based on which we derive the parameters. The computational time needed for calculating a thousand power flow samples is approx. 17 s. The time needed for solving the regression problem for all CLA parameter calculations is around two minutes. By far, the most time consuming is the process of creating the optimization model for the regression problem; it takes approximately 10 min, due to YALMIP overhead.
We proceed by calculating nodal flexibility by using Algorithm 1. We assume active/reactive power flexibility availability ( P f l e x , Q f l e x ) of −50% to 50% of nominal power at all nodes. The algorithm is executed twice: once for lower and once for upper voltage issues mitigation. Total execution time for lower and upper voltage cases is 3 and 6 s, respectively. The results of the calculations can be seen in Figure 3 for the lower voltage case and in Figure 4 for the upper voltage case. From the figures, we can see that in order to elevate the voltages to above a 0.96 p.u. limit, flexibility needs to be used in nodes 7, 8, 10−15, 17–18, 24–25, and 29−32. In order to lower the voltages to below a 1.04 p.u. limit, flexibility needs to be used in nodes 7, 8, 13, 14, 16, 18, 24–25, and 29−33.

6.4. Critical Voltage Threshold

To determine the voltage threshold, we arbitrarily choose node 7 and set t h r = 7 . Network parameters and bounds remain the same. We a set step increase of η to 0.001. We then proceed by employing the proposed Algorithm 2 to find the upper and lower voltage thresholds for voltage at a specified node.
The algorithm required only one iteration for lower and upper voltage thresholds, each taking less than a second of computational time. The calculated lower and upper voltage thresholds are 0.9816 and 1.018 p.u., respectively.

6.5. Validating Results of IEEE 33 Bus Test Case

We validate the results using MATPOWER simulations again. We create samples by randomly setting reference voltage and power injections within the predefined ranges. Our validation code executes the following logic for every sample: Run power flow and check the voltage magnitude at node 7—if it is lower than the lower threshold of 0.9816 p.u., ‘activate’ pre-calculated flexibility, by setting the value of power injection at node 7. Then, run the power flow again for that sample and check again for the voltage at node 7. If it is still below lower threshold, we ‘activate’ flexibility at node 8 in the same manner as described previously. The process repeats until either the voltage at node 7 is above the threshold value or if there is no more flexibility that can be provided. If at any point the voltage at node 7 is above the threshold, move on to the next power flow sample. The logic is the same for the upper voltage bound. This procedure closely resembles the real-life handling of flexibility, where flexibility needs are determined in advance, reserved (as in our Flexibility Provision Problem), and then activated as needed.
We run five thousand samples in MATPOWER and apply the aforementioned logic to each one. The results of this simulation are shown in Figure 5. As can be seen, there are no instances of voltage magnitude crossing the 0.96. p.u. or 1.04 p.u. voltage limits.

6.6. UK 124 Bus Test Case Analysis

We now proceed by analyzing a more realistic and demanding test case from [20], as shown in Figure 6. It is a medium-voltage 6.6 kV weakly meshed synthetic UK distribution network with 124 buses and 124 lines. Node 7 is a reference node with voltage ranging from 0.98 to 1.01 p.u. We assume active/reactive power injections within the range of 80% to 120% of nominal power. Lower and upper voltage limits are set to 0.93 and 1.04 p.u., respectively, for all network nodes. Given the aforementioned assumptions, maximum and minimum achievable node voltages are shown in the first graph in Figure 7. We see that, given a certain combination of reference voltage and power injections, a lower voltage limit can be violated by a large margin, with some voltage magnitudes decreasing as low as 0.88 p.u.

6.7. Flexibility Provision

Again, we want to find the amount of flexibility per node, which can ensure no lower voltage limit violations can occur. We assume flexibility availability at nodes marked ‘FU’ in Figure 6 (nodes: 14, 20, 79, 83, 86, 99, 100, 106, and 107). Active power flexibility ( P f l e x ) ranges from 50% to 100% of nominal active power provided in the dataset. Reactive power flexibility ranges from −250 kVAr to 250 kVAr at FU nodes. CLA parameter calculations is carried out in the same manner as for the IEEE 33 bus test case. Given a much larger network, we resort to commercial solver CPLEX [21] version 12.9. to solve the regression problem for CLA parameter calculations. YALMIP takes approximately 120 s to construct the model, while CPLEX takes around 1000 s to solve the regression problem. The flexibility calculation requires only one iteration with an execution time of 4 s. The results of the calculations are presented in Figure 8. In order to increase all node voltage magnitudes to above the 0.93 p.u. limit, flexibility needs to be used in nodes 79, 83, 86, 99, 100, and 105.
Figure 7. Initial state of the UK 124 bus test network; (a) minimum (orange) and maximum (blue) achievable nodal voltage magnitudes. (b) Nodal minimum (grey) and maximum (orange) active power injection ranges. (c) Nodal minimum (grey) and maximum (orange) reactive power injection ranges.
Figure 7. Initial state of the UK 124 bus test network; (a) minimum (orange) and maximum (blue) achievable nodal voltage magnitudes. (b) Nodal minimum (grey) and maximum (orange) active power injection ranges. (c) Nodal minimum (grey) and maximum (orange) reactive power injection ranges.
Energies 17 04091 g007

6.8. Critical Voltage Threshold

We arbitrarily choose node 28, and set t h r = 28 and η to 0.001. Network parameters and bounds remain the same. Algorithm 2 finds that the lower voltage threshold for the voltage at node 28 is 0.9799 p.u. The algorithm required only one iteration, which took less than a second of computational time.

6.9. Validating Results of UK 124 Bus Test Case

The methodology used for validating results remains the same as for the IEEE 33 bus test case. We create random power flow samples with reference voltage and power injections within predefined ranges of 0.98 to 1.01 p.u. for the reference voltage and 80% to 120% of nominal active and reactive power injections, as defined in the dataset. We run power flow simulations in MATPOWER where, for each sample, we check if the calculated voltage magnitude at node 28 is below threshold of 0.9799 p.u. If so, we activate the flexibility sequentially, until the voltage is restored to above the threshold value. The results of this simulation are shown in Figure 9. As can be seen, there are no instances where the voltage magnitude is below the 0.93 p.u. limit.

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.

Author Contributions

Conceptualization, M.B.; methodology, M.B.; software, M.B.; validation, M.B. and T.C.; formal analysis, M.B.; investigation, M.B.; resources, T.C.; data curation, M.B.; writing—original draft preparation, M.B.; writing—review and editing, T.C.; visualization, M.B.; supervision, T.C.; project administration, T.C.; funding acquisition, T.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the European Union’s Horizon Europe Framework Programme: HORIZON-WIDERA-2023-ACCESS-04, Pathways to Synergies Coordination and Support Actions—under the project name SynGRID—Creating synergies in Widening countries on the topic of low-voltage grid management (grant number 101160145). The funding bodies had no involvement in the preparation of this manuscript.

Data Availability Statement

Data are contained within this article.

Conflicts of Interest

Author Martin Bolfek was employed by the company Croatian Distribution System Operator HEP ODS, Ltd. The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The funders had no role in the design of this study; in the collection, analysis, or interpretation of data; in the writing of this manuscript; or in the decision to publish the results.

Appendix A. Conservative Linear Approximations (CLA)

At the core of the optimization problems proposed in this paper is the novel approach known as conservative linear approximations, first introduced in [24] and subsequently exploited in [25,26]. The authors propose a form of data-driven power flow approximations that establish linear relationships between the voltage magnitudes at a particular node and the power injections at all nodes. Instead of using the classical approach, whereby a relationships between the quantities of interest is derived by solving a set of nonlinear power flow equations, the aforementioned paper assumes a linear relationship between the voltage and a set of power injections at each node in the network:
V i g i ( P , Q ) 0 ( 0 ) ; for over ( under ) estimation
where
g i ( P , Q ) = a i , 0 + a i , 1 T P Q
In order to form the relationship between voltage and power injections, a number of nonlinear power flow samples are run on a given network. A sample is derived by randomly generating power injections for each node and voltage magnitude at a reference node, which are within the predetermined bounds. For a more sophisticated approach to scenario generation and sample selection, see [24]. The results of the calculations, in these case node voltage magnitudes, are saved for the next step where a constrained linear regression problem in the form of (A3)–(A5) is solved, as follows:
min g ( · ) m = 1 M L ( V i , m g i ( P m , Q m ) )
s . t . V i , m g i ( P m , Q m ) 0 ; for overestimation
V i , m g i ( P m , Q m ) 0 ; for underestimation
where m denotes a nonlinear power flow sample, M is a total number of samples, i is a node index, and L is a chosen loss function, in this case a 1-norm, which represents the approximation error. Vectors are expressed in bold font. P and Q represent nodal active and reactive power injection vectors, respectively.
These approximations can be further improved by using CLA approximations of the squared voltage magnitudes ( v 2 ), which are better in terms of accuracy than approximations of the voltage magnitudes themselves. Hence, the formulation in (A1) becomes
v ¯ i 2 = a i , 0 + a i , 1 P Q v m a x 2
v ̲ i 2 = a i , 0 + a i , 1 P Q v m i n 2
where v ¯ i 2 and v ̲ i 2 denote over- and underestimating CLAs, while v m i n and v m a x represent the minimum and maximum voltage limits which can be different for each node. The superscript ⊤ denotes the transpose.
The approximations are called conservative because with regard to the construction (inequalities in the formulation), they overestimate or underestimate a certain quantity of interest, in our case the node voltage magnitude. This property of CLAs could be extremely useful to the SOs who are accustomed to a more conservative approach to day-to-day network operation.

Appendix A.1. Extending the CLAs

The CLAs can be easily extended to any quantity of interest (i.e., line current instead of node voltage). Here, we extend the CLAs to account for the change in voltage at the reference node, which can vary substantially and have a great impact on voltage magnitudes on the rest of the network. Therefore, we expand on Equations (A6) and (A7) by adding V r e f to the vector of node power injections:
v ¯ i 2 = a i , 0 ¯ + a ¯ i , 1 T V r e f P Q v m a x 2
v ̲ i 2 = a ̲ i , 0 + a ̲ i , 1 V r e f P Q v m i n 2
The regression problem from (A3) is modified accordingly, and it becomes
min a 0 , a 1 m = 1 M L v m 2 a 0 + a 1 V r e f P Q
s . t . v i , m 2 a ¯ i , 0 + a ¯ i , 1 V r e f P Q 0 ; for overestimation
v i , m 2 a ̲ i , 0 + a ̲ i , 1 V r e f P Q 0 ; for underestimation
This extension necessitates an increase in the number of samples of power flows needed to adequately capture all of the relationships. However, this is not a major drawback since the calculation of the power flows and linear coefficients is performed offline.

Appendix A.2. CLA Parameter Calculation Algorithm

Algorithm A1 CLA parameter calculation
1:
Input: Network model parameters, bounds on power injections and reference voltage
2:
Generate M random samples of P [ P m i n P m a x ] , Q [ Q m i n Q m a x ] and V r e f [ V m i n r e f V m a x r e f ]
3:
Solve the power flow calculation for each sample m M . Save nodal voltage magnitudes ( v i , m ).
4:
Solve the regression problem of (A10)–(A12) to obtain parameters a ¯ i , 0 , a ¯ i , 1 , a ̲ i , 0 and a ̲ i , 1

References

  1. Akrami, A.; Doostizadeh, M.; Aminifar, F. Power system flexibility: An overview of emergence to evolution. J. Mod. Power Syst. Clean Energy 2019, 7, 987–1007. [Google Scholar] [CrossRef]
  2. Cochran, J.; Miller, M.; Zinaman, O.; Milligan, M.; Arent, D.; Palmintier, B.; O’Malley, M.; Mueller, S.; Lannoye, E.; Tuohy, A.; et al. Flexibility in 21st Century Power Systems; Technical Report; National Renewable Energy Laboratory (NREL): Golden, CO, USA, 2014. [Google Scholar] [CrossRef]
  3. Prettico, G.; Vitiello, A. Distribution System Operator Observatory 2020; Technical Report; Publications Office of the European Union: Luxembourg, 2020. [Google Scholar] [CrossRef]
  4. Vijaychandra, J.; Prasad, B.R.V.; Darapureddi, V.K.; Rao, B.V.; Knypiński, Ł. A Review of Distribution System State Estimation Methods and Their Applications in Power Systems. Electronics 2023, 12, 603. [Google Scholar] [CrossRef]
  5. Pokhrel, B.R.; Bak-Jensen, B.; R. Pillai, J. Integrated Approach for Network Observability and State Estimation in Active Distribution Grid. Energies 2019, 12, 2230. [Google Scholar] [CrossRef]
  6. Millar, B.; Jiang, D.; Haque, M.E. Constrained coordinated distributed control of smart grid with asynchronous information exchange. J. Mod. Power Syst. Clean Energy 2015, 3, 512–525. [Google Scholar] [CrossRef]
  7. Hussain, S.; Lai, C.; Eicker, U. Flexibility: Literature review on concepts, modeling, and provision method in smart grid. Sustain. Energy Grids Netw. 2023, 35, 101113. [Google Scholar] [CrossRef]
  8. Nazir, N.; Almassalkhi, M. Voltage Positioning Using Co-Optimization of Controllable Grid Assets in Radial Networks. IEEE Trans. Power Syst. 2021, 36, 2761–2770. [Google Scholar] [CrossRef]
  9. Lee, D.; Turitsyn, K.; Molzahn, D.K.; Roald, L. Robust AC Optimal Power Flow with Robust Convex Restriction. IEEE Trans. Power Syst. 2021, 36, 4953–4966. [Google Scholar] [CrossRef]
  10. Zhou, S.; Zhao, Y.; Gu, W.; Wu, Z.; Li, Y.; Qian, Z.; Ji, Y. Robust energy management in active distribution systems considering temporal and spatial correlation. IEEE Access 2019, 7, 153635–153649. [Google Scholar] [CrossRef]
  11. Evangelopoulos, V.A.; Avramidis, I.I.; Georgilakis, P.S. Flexibility Services Management under Uncertainties for Power Distribution Systems: Stochastic Scheduling and Predictive Real-Time Dispatch. IEEE Access 2020, 8, 38855–38871. [Google Scholar] [CrossRef]
  12. Molzahn, D.K.; Roald, L.A. Grid-aware versus grid-agnostic distribution system control: A method for certifying engineering constraint satisfaction. In Proceedings of the Annual Hawaii International Conference on System Sciences, Maui, HI, USA, 8–11 January 2019; Volume 2019, pp. 3445–3454. [Google Scholar] [CrossRef]
  13. Ross, S.C.; Mathieu, J.L. A Method for Ensuring a Load Aggregator’s Power Deviations Are Safe for Distribution Networks. Electr. Power Syst. Res. 2020, 189, 106781. [Google Scholar] [CrossRef]
  14. Girigoudar, K.; Roald, L.A. Identifying Secure Operating Ranges for der Control Using Bilevel Optimization. IEEE Trans. Smart Grid 2024, 15, 2921–2933. [Google Scholar] [CrossRef]
  15. Lesage-Landry, A.; Callaway, D.S. Batch reinforcement learning for network-safe demand response in unknown electric grids. Electr. Power Syst. Res. 2022, 212, 108375. [Google Scholar] [CrossRef]
  16. Baran, M.; Wu, F. Network reconfiguration in distribution systems for loss reduction and load balancing. IEEE Trans. Power Deliv. 1989, 4, 1401–1407. [Google Scholar] [CrossRef]
  17. Löfberg, J. YALMIP: A toolbox for modeling and optimization in MATLAB. In Proceedings of the IEEE International Symposium on Computer-Aided Control System Design, Taipei, Taiwan, 2–4 September 2004; pp. 284–289. [Google Scholar] [CrossRef]
  18. Maher, S.J.; Fischer, T.; Gally, T.; Gamrath, G.; Gleixner, A.; Gottwald, R.L.; Hendel, G.; Koch, T.; Lübbecke, M.E.; Miltenberger, M.; et al. The SCIP Optimization Suite 4.0. Technical Report, Optimization Online. 2017. Available online: http://www.optimization-online.org/DB_HTML/2017/03/5895.html (accessed on 15 May 2018).
  19. Zimmerman, R.D.; Murillo-Sánchez, C.E.; Thomas, R.J. MATPOWER: Steady-state operations, planning, and analysis tools for power systems research and education. IEEE Trans. Power Syst. 2011, 26, 12–19. [Google Scholar] [CrossRef]
  20. Martinez-cesena, E.; Churkin, A. Synthetic Electricity Distribution Network from UK for Flexibility Analysis—ATTEST Project. 2022. Available online: https://zenodo.org/records/7502542 (accessed on 2 July 2024).
  21. IBM. IBM ILOG CPLEX Optimization Studio, Version 12.9. 2019. Available online: https://www.ibm.com/products/ilog-cplex-optimization-studio (accessed on 13 June 2024).
  22. Lin, Z.; Wu, Q.; Chen, H.; Ji, T.; Xu, Y.; Sun, H. Scenarios-Oriented Distributionally Robust Optimization for Energy and Reserve Scheduling. IEEE Trans. Power Syst. 2023, 38, 2943–2946. [Google Scholar] [CrossRef]
  23. Ding, B.; Li, Z.; Li, Z.; Xue, Y.; Chang, X.; Su, J.; Jin, X.; Sun, H. A CCP-based distributed cooperative operation strategy for multi-agent energy systems integrated with wind, solar, and buildings. Appl. Energy 2024, 365, 123275. [Google Scholar] [CrossRef]
  24. Buason, P.; Misra, S.; Molzahn, D.K. A sample-based approach for computing conservative linear power flow approximations. Electr. Power Syst. Res. 2022, 212, 108579. [Google Scholar] [CrossRef]
  25. Gupta, R.K.; Buason, P.; Molzahn, D.K. Fairness-Aware Photovoltaic Generation Limits for Voltage Regulation in Power Distribution Networks Using Conservative Linear Approximations. In Proceedings of the 2024 IEEE Texas Power and Energy Conference (TPEC), College Station, TX, USA, 13–14 February 2024; pp. 1–6. [Google Scholar] [CrossRef]
  26. Buason, P.; Misra, S.; Talkington, S.; Molzahn, D.K. A data-driven sensor placement approach for detecting voltage violations in distribution systems. Electr. Power Syst. Res. 2024, 232, 110387. [Google Scholar] [CrossRef]
Figure 1. Minimum achievable voltages across the IEEE 33 bus network with 50% of nominal power injection variability and reference voltage variability between [0.97 and 1.03] p.u.
Figure 1. Minimum achievable voltages across the IEEE 33 bus network with 50% of nominal power injection variability and reference voltage variability between [0.97 and 1.03] p.u.
Energies 17 04091 g001
Figure 2. Initial state of the IEEE 33 bus test network; (a) minimum (orange) and maximum (blue) achievable nodal voltage magnitudes. (b) Nodal active power injection ranges. (c) Nodal reactive power injection ranges.
Figure 2. Initial state of the IEEE 33 bus test network; (a) minimum (orange) and maximum (blue) achievable nodal voltage magnitudes. (b) Nodal active power injection ranges. (c) Nodal reactive power injection ranges.
Energies 17 04091 g002
Figure 3. Output of Algorithm 1. (a) Minimum achievable nodal voltage magnitudes, represented by orange circles. (b) Corresponding node power injections. Blue bars represent real, orange bars represent reactive power, and yellow/purple combination represents calculated active/reactive flexibility provision.
Figure 3. Output of Algorithm 1. (a) Minimum achievable nodal voltage magnitudes, represented by orange circles. (b) Corresponding node power injections. Blue bars represent real, orange bars represent reactive power, and yellow/purple combination represents calculated active/reactive flexibility provision.
Energies 17 04091 g003
Figure 4. Output of Algorithm 1 for IEEE 33 bus test case. (a) Maximum achievable nodal voltage magnitudes, represented by orange circles. (b) Corresponding node power injections. Blue bars represent real, orange bars represent reactive power, and yellow/purple combination represents calculated active/reactive flexibility provision.
Figure 4. Output of Algorithm 1 for IEEE 33 bus test case. (a) Maximum achievable nodal voltage magnitudes, represented by orange circles. (b) Corresponding node power injections. Blue bars represent real, orange bars represent reactive power, and yellow/purple combination represents calculated active/reactive flexibility provision.
Energies 17 04091 g004
Figure 5. Nodal voltage magnitudes of IEEE 33 bus test case as a result of five thousand MATPOWER power flow samples.
Figure 5. Nodal voltage magnitudes of IEEE 33 bus test case as a result of five thousand MATPOWER power flow samples.
Energies 17 04091 g005
Figure 6. Topology of UK 124 bus test case.
Figure 6. Topology of UK 124 bus test case.
Energies 17 04091 g006
Figure 8. Output of Algorithm 1 for UK 124 bus test case. (a) Minimum achievable nodal voltage magnitudes, represented by orange circle. (b) Corresponding node power injections. Blue bars represent real, orange bars represent reactive power, and yellow/purple combination represents calculated active/reactive flexibility provision.
Figure 8. Output of Algorithm 1 for UK 124 bus test case. (a) Minimum achievable nodal voltage magnitudes, represented by orange circle. (b) Corresponding node power injections. Blue bars represent real, orange bars represent reactive power, and yellow/purple combination represents calculated active/reactive flexibility provision.
Energies 17 04091 g008
Figure 9. Nodal voltage magnitudes of UK 124 bus test case as a result of five thousand MATPOWER power flow samples.
Figure 9. Nodal voltage magnitudes of UK 124 bus test case as a result of five thousand MATPOWER power flow samples.
Energies 17 04091 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bolfek, M.; Capuder, T. Voltage-Triggered Flexibility Provision in a Distribution Network with Limited Observability. Energies 2024, 17, 4091. https://doi.org/10.3390/en17164091

AMA Style

Bolfek M, Capuder T. Voltage-Triggered Flexibility Provision in a Distribution Network with Limited Observability. Energies. 2024; 17(16):4091. https://doi.org/10.3390/en17164091

Chicago/Turabian Style

Bolfek, Martin, and Tomislav Capuder. 2024. "Voltage-Triggered Flexibility Provision in a Distribution Network with Limited Observability" Energies 17, no. 16: 4091. https://doi.org/10.3390/en17164091

APA Style

Bolfek, M., & Capuder, T. (2024). Voltage-Triggered Flexibility Provision in a Distribution Network with Limited Observability. Energies, 17(16), 4091. https://doi.org/10.3390/en17164091

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop