1. Introduction
As global concerns about energy, the environment, and climate issues escalate, the advancement of clean and renewable energy sources, including wind power and photovoltaic (PV) energy, has emerged as a crucial endeavor across numerous countries [
1]. In this context, distributed PV power generation and other related technologies have experienced rapid development. However, compared with traditional thermal power and hydropower, etc., PV power generation is characterized by intermittency, stochasticity, and uncertainty [
2]. Because of the unpredictable nature of load demand and renewable energy sources like PV power generation, it may lead to intermittent voltage violations and other problems in the distribution system, which significantly affects its safe operation [
3,
4]. Therefore, studying the voltage optimization problem within distribution systems with high penetration of PV energy is of utmost importance.
In view of the above problems, the traditional voltage regulation method has certain defects. For example, transformer tap switch regulation and capacitor regulation are discrete voltage regulation methods [
5], which have high operating costs and cannot be operated frequently. Compared with the traditional voltage regulation methods, the inverter of the PV power generation system has a preset redundant capacity, which can be used for reactive power regulation more quickly and flexibly [
6]. At the same time, the energy storage system (ESS) is connected to the distribution grid through an inverter device, which has the ability of fast and flexible regulation in four quadrants [
7]. Therefore, more flexible and fast voltage regulation can be realized by utilizing the PV power generation system and ESS, which has prompted extensive research. Ref. [
8] mitigates the violation of voltage by adjusting the output of the PV inverter, realizing voltage regulation through operation modes such as preventive control, emergency control, and restoration control. Ref. [
9] proposes a short-term ESS for improving the voltage quality problem in distributed wind power systems. Ref. [
10] proposes a distributed control strategy for coordinating multiple ESSs for voltage regulation in active distribution networks with high penetration of fluctuating renewable energy sources. Ref. [
11] investigates a method to improve voltage quality in low-voltage distribution networks by utilizing the reactive power capability of single-phase PV inverters.
With the development of active distribution network technology, new technical means have emerged for voltage control in the distribution network, and flexible distribution network reconfiguration can effectively improve the voltage distribution. For example, Ref. [
12] establishes a two-stage optimization model to reduce the impact of PV power generation on voltage by reducing the equivalent impedance of the Davignan model at the common connection point through distribution network reconfiguration. Ref. [
13] proposes a mixed-integer programming model for the reconfiguration of active distribution systems that considers voltage dependency, the type of loads, and renewable sources. Ref. [
14] achieves preventive voltage control in distribution networks by reconfiguring the network. This approach addresses voltage security issues and ensures the safe operation of distribution networks. Ref. [
15] improves voltage distribution within the distribution network by coordinating the regulation of both active and reactive power of distributed generations (DGs). Ref. [
16] proposes coordinating the operation of multiple voltage regulators and DGs in the distribution system to achieve effective voltage control under changes in network structure and DG availability.
Utilizing the flexibility of multiple source-network-load-storage resources of the distribution network, and considering the coordination of the above resources, the voltage optimization of the distribution network can be realized more effectively [
17,
18]. For example, Ref. [
19] proposes a real-time voltage regulation strategy for reconfigurable low-voltage distribution networks. By coordinating the optimal feedforward control of DGs with the operation of line switches, pre-compensation of impending voltage deviations during network reconfiguration-assisted load restoration is achieved. Ref. [
20] explores techniques for managing voltage fluctuations in low-voltage distribution systems with a significant presence of PV penetration. Where, a wide range of applications of three types of flexible resources, namely source, network, and storage, are considered for voltage regulation of low-voltage distribution networks. Ref. [
21] considers the coordination of electric vehicles and PVs to achieve voltage optimization in distribution networks with active network losses and voltage fluctuations as optimization objectives. Ref. [
22] achieves distribution network voltage quality optimization through dynamic reconfiguration of the distribution network in concert with energy storage and other synergies. Ref. [
23] presents a control framework based on source-grid-load collaboration to enhance power quality in distribution networks with high penetration of low-carbon energy sources. Ref. [
24] proposes optimal locations for distributed PV generation to identify and mitigate voltage limit violations.
Although many relevant studies have been carried out on the voltage optimization problem of distribution networks under high PV penetration, existing studies have not simultaneously considered the coordination of PV generation systems, network reconfiguration, load demand, and ESS in different time scales. So, achieving efficient utilization of renewable energy and ensuring the economic operation of the grid becomes challenging. Therefore, this paper proposes an active distribution network voltage optimization method based on the coordination and interaction of source-network-load-storage resources. Where, in the day-ahead prediction stage, based on the forecasted scenarios of PV generation and load demand, the reconfiguration scheme for the distribution network is formulated, and the voltage distribution is improved by optimizing the topology of the network. In the real-time operational phase, taking into account the variability of load demand and PV generation across different scenarios, the flexibility of PV inverters and ESSs is harnessed to make instantaneous adjustments. The active and reactive power outputs of the PV generation systems, as well as the charging and discharging power of the ESSs, are fine-tuned in real-time. This dynamic adaptation is implemented according to the differing scenarios of load demands and PV outputs. Consequently, the voltage quality of the distribution network is further enhanced, and the capacity to utilize renewable energy is effectively increased.
The main contributions of this paper are summarized as:
- (1)
An active distribution network voltage optimization method based on source-network-load-storage coordination and interaction is proposed.
The method leverages the flexibility of PV generation systems, ESSs, load demand, and network reconfiguration to optimize the voltage quality of distribution networks. By formulating a distribution network reconfiguration plan in the day-ahead prediction stage and adjusting the active and reactive power outputs of PV systems and charging/discharging powers of ESSs in real-time, the method effectively improves voltage distribution and ensures safe and economic operation.
- (2)
A 0–1 mixed-integer second-order cone programming model based on the scenario-based stochastic optimization for voltage optimization is established.
Taking into account the uncertainties associated with PV output and load demand, a scenario-based stochastic optimization model is devised to optimize voltage levels. This model encompasses network topology constraints, power flow and security constraints, load demand constraints, PV generation system operational constraints, and ESS operational constraints, among others. Structured as a 0–1 mixed-integer second-order cone programming problem, this model offers enhanced extensibility, making it easier to incorporate additional voltage regulation measures in the future.
- (3)
The scenario decomposition algorithm is introduced as an efficient solution.
To tackle the complexity of the large-scale stochastic optimization problem, an efficient solution algorithm based on the scenario decomposition algorithm is employed. This algorithm decomposes the original large-scale problem into multiple smaller, independent subproblems, each considering a single scenario. This provides potential for parallel computing on multiple computers.
The rest of the paper is organized as follows:
Section 2 describes the model framework and the formulation of the proposed model;
Section 3 elaborates on the model solution method;
Section 4 verifies the effectiveness of the model using the IEEE 33-node distribution system; and
Section 5 summarizes this paper.
3. Solution Method
Constraints (14), (21), (32), and (39) contain quadratic terms, while the remaining constraints are linear. Furthermore, constraints (14), (21), (32), and (39) can be re-expressed in the form of the second-order cone, which are shown as constraints (45)–(48). Consequently, the original optimization problem is transformed into a 0–1 mixed-integer second-order cone programming (MISOCP) problem.
Further, considering the two-stage form characteristic of the model, the scenario decomposition algorithm (SDA) is adopted to decompose the large-scale problem into multiple small-scale problems [
29]. The detailed procedure for finding the solution is outlined below:
Step 1: Let the number of scenarios be
N and
k denote a particular scenario number (
). Then the original model can be decomposed into the following
N subproblems.
Step 2: Since subproblem SPk only considers scenario k to form the corresponding constraints instead of all scenarios, which is essentially equivalent to relaxing the constraints of the initial problem, the objective obtained from solving the N subproblems individually will be the lower bound of the objective, which is denoted as LB.
Step 3: In order to obtain the upper bound UB of the objective , substitute the optimal solution of SPk into all the other subproblems SPm(m ≠ k). Since the optimal solution of the subproblem SPk is not necessarily optimal for the other subproblems SPk, the final objective function value obtained by substituting it into the other subproblems will be the upper bound. For each subproblem SPk, a corresponding upper bound will be obtained. For each solution of the subproblem SPk, substituting it into other subproblems will obtain a corresponding upper bound, denoted as UBk. Then .
Step 4: When
UB and
LB do not converge, it means that the current solution (denoted as
) is not optimal for the initial problem. Therefore, it is necessary to eliminate the current solution. Given that
is a 0–1 variable, the elimination of the non-optimal solution constraint is shown in constraint (49). Then, the constraint (49) is added to the subproblem, and Step 3 and Step 4 are repeated until the
UB and
LB converge, and the corresponding subproblem solution is optimal for the initial problem.
The use of relaxed scenario subproblems reduces the problem size compared to the original model. It is noteworthy that the benefit of utilizing the SDA lies in that the subproblems of each scene are independent, and these subproblems can be taken to be computed in parallel, which opens up the possibility of solving large-scale optimization problems.
The flowchart describing the steps proposed is shown in
Figure 2.
4. Case Study
4.1. System Parameters
In this section, the IEEE 33-node distribution system is used to verify the effectiveness of the proposed method [
30], of which the topology is shown in
Figure 3. The system contains a total number of 33 nodes and 37 distribution lines, with Node 1 being the substation bus and industrial loads at other nodes. The reference voltage is 12.66 kV and the voltage range is [0.95, 1.05] p.u. The ESS is located at Nodes 11 and 31, each having a capacity of 500 kWh. The maximum charge/discharge power is 100 kW, and the efficiency of charging/discharging is 94%. The PV generation system is located at Nodes 15 and 24, with a capacity of 300 kVA.
,
,
, and
are all 0.25.
and
are set to be 0.4 and 0.6, respectively.
,
,
,
, and
cost coefficients are set to be 100, 1, 100, 100, 100, and 1, respectively.
The PV and load demand prediction curves in the day-ahead prediction stage are shown in
Figure 4, respectively. The load demand and solar power output data are based on actual measurements and reliable estimates derived from historical data. The load demand is higher during peak hours (e.g., mornings and evenings) and lower during off-peak hours (e.g., overnight). Conversely, solar power output is typically higher during daylight hours, peaking around midday, and lower or non-existent at night. This alignment with general patterns ensures that the scenarios are representative of real-world conditions.
A total of four possible PV output and load demand scenarios are considered. The output of PVs and demand of loads for each scenario in the real-time operation stage are, respectively, shown in
Figure 5,
Figure 6 and
Figure 7. The probability of each scenario is 0.4, 0.1, 0.2, and 0.3, respectively. Other parameters are detailed in Ref. [
31].
The algorithms are compiled in MATLAB (R2014a), modeled using YALMIP (R20230622), and solved using GUROBI (9.5.2). The computer is configured with an Intel Core i5-10400F CPU (Intel, St. Clara, CA, USA) operating at 2.90 GHz and features 16.0 GB of RAM.
4.2. Result Analysis
The result of network reconfiguration is shown in
Figure 8. The calculation time is 25,459 s, which is considered acceptable given that the ultimate outcome of this solution is a network reconfiguration plan for the day-ahead prediction stage. It should be noted that the voltage at Node 1, representing the substation node, is not assumed to be a constant but rather dynamically optimized through the coordination and optimization of source-network-load-storage resources.
In
Figure 8, the reconfigured network topology shows the optimized connection between nodes and lines to optimize the voltage across the distribution system.
It is worth mentioning that, the objective function minimizes costs related to wasted PV energy, load loss, ESS operation, and voltage deviation. Constraints ensure stability and safety by regulating network topology, power flows, loads, PV, ESS operations, and so on. During optimization, the objective function adapts to each scenario. For example, in Scenario 4, with high PV output and large load fluctuations, the objective prioritizes minimizing PV energy waste (Equation (41)) by adjusting ESS charging/discharging (Equations (33)–(39)) and PV reactive power (Equations (31) and(32)) to balance supply and demand, while reducing load loss (Equation (43)) and voltage deviation (Equation (44)).
Taking Node 2 as an example (
Figure 9), despite fluctuations in PV output and load demand across scenarios, the voltage remains stable within the allowed range (0.95–1.05 p.u.), indicating that the optimization method successfully adjusts PV output, ESS charging/discharging, and network topology to satisfy voltage constraints and minimize costs. Similarly, for Node 17 (
Figure 10), the optimization ensures voltage stability across diverse scenarios, demonstrating the effectiveness of the proposed method in balancing costs and ensuring system stability under high PV penetration.
To elaborate on the results, taking Scenario 4 as an example, the active output of each resource is shown in
Figure 11; the total active power of the ESSs and the total actual active PV power are shown in
Figure 12; and the total reactive PV power is shown in
Figure 13.
Combining
Figure 6,
Figure 11 and
Figure 12, it is observable that during the period from 0 h to 8 h, the ESSs are primarily in a charging mode because the active load demand is relatively low. When moving into the 8 h–17 h period, the active load demand starts to increase notably. However, concurrently, the active PV output also experiences a surge. To ensure efficient utilization of PV energy and minimize energy wastage, the ESSs refrain from engaging in significant charging and discharging activities during this timeframe. In the subsequent period, from 17 to 24 h, the active load demand continues to escalate, whereas the PV system’s output starts to decline. To maintain equilibrium between power supply and demand, the ESSs transition into a discharging state. Notably, in Scenario 4, due to the substantial fluctuations in load during this period, the discharging power of the ESSs also exhibits corresponding variations.
Combining
Figure 11,
Figure 12 and
Figure 13, it is observable that, to achieve the effective utilization of PV, the reactive power output is reduced accordingly during the higher active power output period (i.e., time 8 h–17 h) to avoid the waste of PV energy.
Combining
Figure 7 and
Figure 13, it is obvious that the reactive power output of the ESSs is maintained at a low level during the time period when the load reactive power demand is relatively low (i.e., time 0 h–7 h). However, as the load level increases, the reactive output of the ESS increases accordingly to ensure the voltage quality of the power system.
In order to further verify the effectiveness of the proposed method, the following comparison cases are designed to show its superiority over traditional methods:
- (1)
Case 1: the proposed model is adopted.
- (2)
Case 2: the proposed model is adopted without considering network reconfiguration, i.e., the initial topology is fixed.
- (3)
Case 3: the proposed model is adopted without considering the charging/discharging process of ESSs.
- (4)
Case 4: the proposed model is adopted without considering the PV output regulation.
The percentages of average voltage deviation across all nodes in the system display notable variations across different scenarios, as shown in
Figure 14,
Figure 15,
Figure 16 and
Figure 17. Where, “voltage deviation percent” refers to the percentage deviation of the actual voltage from the nominal or reference voltage at each node in the distribution network. Case 1 exhibits lower voltage deviations in all scenarios, particularly during periods of significant load fluctuations (such as time 17 h–24 h in Scenario 4), where the voltage deviation control effect is particularly notable. In Case 2, the voltage deviation is higher than in Case 1, especially when the network topology is unfavorable for voltage regulation, and the problem of voltage violations is more prominent. In Case 3, although it shows some improvement compared to Case 2, due to the lack of flexible regulation capability of the ESS, the voltage deviation is still relatively large. In Case 4, the absence of photovoltaic output leads to a reduction in system adjustable resources, limited voltage regulation capability, and further increased voltage deviation.
The comparative analysis reveals that the crucial elements for achieving voltage optimization in active distribution networks encompass the holistic consideration of coordinated interactions among diverse source-network-load-storage resources, implementing network reconfiguration during the day-ahead phase, and dynamically modulating PV output, ESS operations, and load demand in real-time. The presented approach demonstrates substantial benefits in enhancing voltage quality.
It should be noted that, although the voltage remains within the acceptable limits of 0.95 to 1.05 p.u. in all cases, this is merely the minimum requirement. Even within the acceptable range, deviations from the nominal value can still negatively impact grid stability and equipment operation. Therefore, the optimization ensures that the voltage remains as close as possible to the nominal value (i.e., 1 p.u.), thereby improving voltage quality.
Moreover, the flexibility of the PV system is primarily realized through the reactive power output adjustment of its inverters. The reactive power capability of the inverters enables rapid response to grid voltage fluctuations, optimizing voltage profiles. This adjustment is highly dynamic, responding to real-time grid demands and operational conditions. Due to the inherent uncertainty in the active power output of solar power systems, such as fluctuations caused by weather changes and day-night variations, the integration of ESSs plays a pivotal role. ESSs can absorb excess energy during peak solar generation periods and release it during low-generation times, ensuring a stable power supply and addressing the intermittency of solar power.