- 1College of Energy and Electrical Engineering, Hohai University, Nanjing, China
- 2Smart Grid Research Institute, Nanjing Institute of Technology, Nanjing, China
With a high proportion of photovoltaic (PV) connected to the active distribution network (ADN), the correlation and uncertainty of the PV output will significantly affect the grid dispatching operation. Therefore, this paper proposes a novel robust ADN dispatching model, which considers the dynamic spatial correlation and power uncertainty of PV. First, the dynamic spatial correlation of PV output is innovatively modeled by dynamic conditional correlation (DCC) generalized autoregressive conditional heteroskedasticity (DCC-GARCH) model. DCC can accurately represent and forecast the spatial correlation of the PV output and reflect its time-varying characteristics. Second, a time-varying ellipsoidal uncertainty set constructed using the DCC, is introduced to bound the uncertainty of the PV outputs. Subsequently, the original mixed integer linear programming (MILP) model is transformed into the mixed integer robust programming (MIRP) model to realize robust optimal ADN dispatching. Finally, a numerical example is provided to demonstrate the effectiveness of the proposed method.
1 Introduction
As a clean and sustainable renewable energy source, photovoltaic (PV) generation has become one of the world’s fastest-growing energy sources, progressively becoming the primary source of electricity in power system (Calcabrini et al., 2019). However, owing to the correlation and uncertainty of PV output, the high proportion of PV in the active distribution network (ADN) significantly influences ADN operation and increases the complexity of dispatching (Yu et al., 2015; Haque and Wolfs, 2016). Thus, a reasonable consideration of the correlation and uncertainty of PV outputs in the ADN dispatching model can help enhance PV consumption and promote the balanced development of the ADN.
First, the modeling of PV-output spatial correlation is considered. According to earlier research, there is a spatial correlation between the PV output attributed to locational considerations, micrometeorological circumstances (Ding and Mather, 2017). Currently, common methods to establish a PV-output spatial-correlation model mainly include correlation coefficient matrices, covariance matrix, Copula function, and deep neural network. In (Wu et al., 2021), the Pearson correlation coefficient matrix was used to calculate the spatial correlation of the two PVs, from which empirical distributions of spatial correlation coefficients and distances were obtained. In (Luo et al., 2020), the Kendall rank correlation coefficient, which can measure the correlation of nonlinear variables, was used as the parameter of the Frank-Copula function. Various correlation coefficient matrices were used to evaluate and compare the numerical values of PV-output spatial correlation. In (Wu et al., 2022), a covariance-based spatiotemporal correlation model was proposed to quantify and exploit the PV output spatial correlation. In (Pan et al., 2019), a variety of Copula functions were selected to obtain an appropriate PV spatial correlation expression. In (Zamee and Won, 2020), Spearman rank-order correlation and an artificial neural network (ANN) were combined to characterize PV output. Each method has its own advantages and disadvantages, as listed in Table 1. Overall, neither the numerical expression nor the model establishment of PV-output spatial correlation can measure or represent the state of spatial correlation at each moment; that is, they cannot reflect the time-varying characteristics of PV-output spatial correlation. Furthermore, there is a lack of prediction models for multidimensional PV output spatial correlations in existing studies.
Second, the modeling of PV-output uncertainty in the optimal dispatching problem of an ADN is investigated. Stochastic optimization (SO) and robust optimization (RO) are the two primary forms of modeling optimization methodologies for PV-output uncertainty (Aharon and Laurent El, 2009). For instance, in (Liu et al., 2019), a day-ahead economic scheduling method based on chance-constrained programming was proposed considering the uncertainty of PV output. In (Vilaça Gomes et al., 2019), a new approach was formulated to model the uncertainty of the wind-sun-hydrothermal system by generating several representative scenarios. However, SO requires estimating the probability distribution of variables from historical data, which is prone to large errors in actual situations. In addition, considering the accurate representation of uncertainty, many scenarios may need to be considered, which will increase the computational complexity. Conversely, RO does not need to know the specific probability distribution, but given the values ranges of variables, that is, the uncertainty set (El-Meligy et al., 2022). In (Xu et al., 2020; Choi et al., 2022), the box uncertainty sets were used to restrict the upper and lower bounds of variables, which tend to be over-conservative. (Ji et al., 2019). and (Aghamohamadi et al., 2021) used the polyhedral uncertainty sets, which have a linear structure and can easily control the uncertain budget, to describe the uncertainties of PV outputs. However, it is difficult to depict the correlation between these parameters using box or polyhedral uncertainty sets. Significantly, the ellipsoidal uncertainty sets can effectively describe multi-type sets to facilitate data input and indicate the correlation between uncertain multi-variables (Chassein and Goerigk, 2016; Golestaneh et al., 2018). Considering the dynamic spatial correlation between multiple PV outputs, it is better to use an ellipsoidal uncertainty set to characterize the PV-output uncertainty in the robust optimization dispatching model of ADN.
Based on the above discussions, this study establishes a robust dispatching model for an ADN that fully considers the spatial correlation and uncertainty of PV output. In this model, the dynamic conditional correlation (DCC) generalized autoregressive conditional heteroskedasticity (DCC-GARCH) model is introduced to construct and predict multidimensional dynamic correlation coefficient models of PV outputs, representing temporal changes in spatial correlation. After modeling the dynamic spatial correlation, a time-varying ellipsoidal uncertainty set is introduced to bound the uncertainty of PV outputs. In addition, to build a highly reasonable and widely applicable ADN dispatching model, the coordination optimization of active and reactive power is considered, and a variety of methods are adopted to linearize the nonlinear constraints for an easy solution. Finally, the effectiveness and rationality of the proposed methodologies are demonstrated considering a rural ADN in China as an example. The major contributions of this study are as follows:
1) To characterize the time-varying characteristics of PV-output spatial correlation, the DCC-GARCH model is proposed in this study to establish a dynamic model of the correlation. The model also allows for intraday short-term forecasting of DCC, which can measure the spatial correlation.
2) To address the uncertainty problem of PV output in the optimization dispatching of an ADN, a time-varying ellipsoidal uncertain set based on the DCC-GARCH model is proposed to improve the robustness of photovoltaic output modeling. By introducing the uncertainty set, the MILP model for ADN optimal dispatching can be transformed into mixed integer robust optimization (MIRP).
The remainder of this paper is organized as follows. Section 2 describes the dynamic spatial correlation and proposes a time-varying spatial correlation model for the PV outputs. Section 3 explains the time-varying ellipsoidal uncertainty set of the PV outputs and constructs a robust dispatch model of ADN. In Section 4, the numerical results for a rural ADN in China are presented to verify the effectiveness of the proposed model. Finally, the conclusions are presented in Section 5.
2 Multidimensional dynamic spatial correlation model of PV output based on DCC-GARCH
2.1 Dynamic changes of PV output spatial correlation
Owing to the influence of micro-meteorology, PV power stations in the same area show different output variations. Relevant research has shown that there is a high correlation between the outputs of these PV stations, which is known as spatial correlation. The five PV stations in Zhangjiagang, Suzhou is taken as an example, which are shown in Figure 1. The PV output data in 2018 is used to draw Figure 2, which has been standardized by the installed capacity. And the sampling interval is 5 min.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1012581/fenrg-10-1012581-HTML-r1/image_m/fenrg-10-1012581-g002.jpg)
FIGURE 2. Power distribution and confidence interval of PV1-PV2 and PV1-PV2-PV3 at partial time. (A) PV1-PV2-January, (B) PV1-PV2-March, (C) PV1-PV2-May, (D) PV1-PV2-PV3 April 6th, (E) PV1-PV2-PV3 April 7th, and (F) PV1-PV2-PV3 April 8th.
Figure 2 displays the scatter plots of the power distribution of PVs at various time scales and dimensions, as well as the accompanying confidence intervals. As shown in Figure 2, the confidence ellipses of the same PVs at different times are considerably different, which implies that the spatial correlation of the same PVs varies over the time. In Figures 2A–C, the time scale is a month. The confidence ellipse for January is longer and thinner than those for March and May, resulting in a stronger spatial correlation. Figures 2D–F reflects the changes of three-dimensional (3D) confidence ellipses of PV1, PV2, and PV3 outputs across three consecutive days from April 6th to April 8th. On April 8th, the ellipsoid is the most elongated and has the shape of a prolate ellipsoid. Thus, the spatial correlation among the three PV stations is the strongest for the 3 days. In addition, on comparing Figures 2D–F with Figures 2A–C, it can be observed that the smaller the time scale is, the more noticeable is the fluctuation of spatial correlation of PV output.
Figure 2 illustrates that the spatial correlation of PV output has a significant time-varying feature, with the smaller time scale indicating a larger variation. However, most existing studies have built static models of spatial correlation and do not predict the spatial correlation itself. Based on this, a dynamic spatial-correlation model based on the DCC-GARCH model is proposed in this paper to realize the numerical characterization and prediction of spatial correlation with time-varying properties.
2.2 Dynamic spatial-correlation model and prediction based on DCC-GARCH
Suppose that the output time-series of k PV stations is
Equations 1, 2 establish the ARMA-GARCH model for the output of each PV station. Equation 1 is the mean equation using the ARMA model, and Equation 2 refers to the variance equation of the GARCH model, where the residual term {
Equations 3–8 build the DCC-GARCH model. In Eq. 3,
The DCC (1,1)-GARCH (1,1) model is more frequently. In the simplified model, Eqs 2, 7 are simplified to Eqs 9, 10, respectively. The DCC can then be written as Eq. 11.
where
Consider the covariance matrix of r-step ahead prediction is
where
So far, a dynamic spatial correlation model of PV output can be obtained. In addition to accurately describing the PV output, the model can also assess, compute, and predict the DCC of different PV outputs. In this study, the relative value
3 Robust dispatching model for ADN considering the dynamic PV spatial correlation
In contrast to the transmission network, the resistance and reactance values of the ADN lines are close to each other, and the coupling between active and reactive power is strong (Sun et al., 2022). It is not sufficient to establish a unilateral active- or reactive-power dispatching model based on the traditional active- and reactive-power decoupling theory. The reactive-power resources about the high proportion of PV will affect the network loss and voltage quality (Antoniadou-Plytaria et al., 2017; Hu et al., 2022). Active power optimization can reduce generation costs, and reactive power regulation can reduce network losses, together achieving the goal of minimizing operating costs. Therefore, first, this section establishes the active- and reactive-power coordination dispatching model of the ADN. Thereafter, to facilitate the solution, the nonlinear constraints are linearized to form the MILP model. Finally, considering the uncertainty of PV output, a time-varying ellipsoid uncertainty set is proposed to construct a robust dispatching model of ADN.
3.1 Objective function
Because economy is a significant evaluation indicator for ADN, the lowest operating cost of ADN is chosen as the objective function in this paper. The operational expense consists of the cost of purchasing electricity from the grid, the dispatching cost of curtailable loads, and the lifespan-loss cost of the energy storage system (ESS), which is expressed as:
where Fbuy(t) and FCL(t) denote the cost of electricity purchase and dispatching cost of the curtailable load at time t, respectively, and FEss denotes the lifespan-loss cost of energy storage. In Eq. 15, mbuy(t) denotes the unit cost of electricity purchased from the upper grid. In Eq. 16,
3.2 Constraints
The normal operation of the ADN must satisfy power flow constraints. To achieve reasonably coordinated dispatching within the system, it is also necessary to consider the operating constraints of each device. In this study, we consider the demand-side response constraints and the operating constraints of an ESS, on-load tap changer (OLTC), capacitor bank (CB), static var generator (SVG), and distributed generation (DG), as well as the impact of ESS lifespan losses on operational costs.
3.2.1 Power flow constraints
The injected and outflow powers at each node must be equal for the ADN. Therefore, the branch flow model can be expressed as
Eq. 18 gives the branch power flow equations, where θij,t is the difference in the voltage phase angle between nodes i and j. gij and bij are the electric conductance and susceptance of branch ij, respectively. Eq. 19 shows the node power-balance equation, where Gij and Bij are the conductance and susceptance of the node to the ground, respectively. ϕ is the power-factor angle of the load being reduced. Eq. 20 constraints the capacity of branch ij.
3.2.2 OLTC operation constraints
OLTC can regulate the voltage by adjusting the position of the tap, which is an important component for maintaining voltage stability. A virtual node m can be added to the branches that contain OLTC to separate the transformer branch into an ideal transformer section and a lossy section.
The voltage of the ideal transformer on the secondary side can be expressed as
where δij,t represents the turn ratio of the transformer branch at time t, which can be defined as a linear combination of the following constraints.
where
3.2.3 ESS operation constraints
Constraints (25) and (26) limits the active power of the ESS at node i at time t. Constraint (27) limits the state of charge (SOC) at node i at time t. In constraint (27),
3.2.4 ESS cycle life loss constraints
The ESS loss cost accounts for a significant portion of the ADN operational cost. Because the ESS lifespan-loss might influence economic efficiency during operation, cyclic-life loss constraints for the ESS are constructed in this paper using the method proposed in the literature (Wang et al., 2016).
Assuming that the charging and discharging procedures have the same impact on the cycle lifespan degradation, a full charging and discharging cycle is divided into two distinct processes. The daily degradation is the sum of the degradations during each time interval, as shown in (28). For each period, the cycle lifespan degradation can be calculated by deducting the two regular degradations.
where degt denotes the cycle loss percentage corresponding to the SOC at time t, which can be obtained from the degradation-SOC curve, as shown in Wang et al. (2016).
Given that all the absolute values in Eq. 29 are less than 1, the nonlinear function can be converted into a linear inequality constraint by adding two binary variables d1 and d2. The functions are as follows:
3.2.5 CB and SVG constraints
Reactive-power-compensation components in an ADN typically fall into one of two types: the discrete component CB and the continuous component SVG. The constraints of the CB are as follows:
Constraint (32) constrains the number of operations, where where
SVG is a continuous reactive power compensation device that can effectively respond to the sudden changes in voltage or overvoltage caused by DG fluctuations in an ADN. The constraints of an SVG are as shown in Constraint (33):
where
3.2.6 PV-output constraints
The PV output constraints are presented in Constraints (34) and (35). Constraint (34) denotes the maximum PV active power at time t, whereas
3.2.7 Demand-response constraints
Demand response (DR) is crucial for ADN dispatch and optimization. It can reduce the uneven tide distribution caused by DG uncertainty. By tracking DG generation, DR can provide some regulation capability for the ADN. In this study, the main DR we considered is the curtailable load with the following constraints:
where
3.3 Model linearization
3.3.1 Successive linear approximation of power flow
Despite the high computational accuracy of AC power flow, its non-convex nonlinear properties render it unsuitable for integration into complex distribution network optimization issues. Consequently, this study linearizes the power flow equation to improve the solution efficiency while maintaining appropriate computational accuracy. It then creates a linearized optimization model for an ADN.
Constraint (18) is nonlinear because of the product of a trigonometric function with the voltage. Suppose that an initial point (vi,j,k, θi,j,k) is provided. The first-order Taylor series expansion of the sine and cosine functions is expressed as:
where,
The power flow equation is formulated by incorporating (42) and (43) into Equation 18.
With
The vi,tvj,t can be uncoupled as (46). Subsequently, the first-order Taylor expansion can be used to further approximate the linearization of
The linearized power flow constraints can be obtained by substituting (45)–(47) into (44), as follows:
where
In Yang et al. (2016) and Yang et al. (2017), it has been demonstrated by the examination of numerous instances that successive linear approximation of power flow has high accuracy. Its efficiency has a sizable advantage over heuristic methods owing to the rapid development of commercial linear-programming tools. Furthermore, compared with the second-order cone relaxation, the successive linear approximation, which is based on the Taylor series, has unrestricted objects and superior scalability.
3.3.2 Linearization of line capacity and DG capacity
Constraints (20), (35) are elliptical nonlinear constraints on capacity. Linearization can be realized using a linear approximation with multiple rectangular constraints.
The line capacity constraint is transformed into (50). The capacity constrained of DG inverter is transformed into (51).
3.3.3 Linearization of ESS cycle lifespan loss constraints
The constraints of the ESS cycle lifespan loss are non-convex and nonlinear and have integer variables. Through the second type of special-order sets (SOS2), the degradation-SOC curve can be linearized and is expressed as:
where m is the number of SOC-curve segments of the ESS. Δt,m and δt,m are the SOS2 variables at time t.
3.4.4 Linearization of OLTC operating constraints
There are integer variables in the OLTC constraints. To avoid the dimensional disaster caused by the non-deterministic polynomial (NP) problem, the OLTC constraints are linearized using following procedure:
Suppose λij,t,n is a binary variable,
where
While defining mij,t = δij,tUj,t, hij,t = λij,t,nUj,t, and gij,t,n = λij,t,nmij,t, Eqs 59, 60 can be respectively obtained by multiplying both sides of Equation 22 by Uj,t and mij,t.
Using the Big-M method, two equivalences can be achieved by introducing a large number M:
Thus far, OLTC operation constraints have been linearized, converting the model to a mixed integer model. Additionally, daily modifications of the OLTC must be restricted, as shown in (63).
where δij,t,up and δij,t,do are binary variables that define the increasing and decreasing of OLTC ratio at time t.
3.4 Time-varying ellipsoidal uncertainty set of PVs
It is vital to describe the uncertainty set for many optimizations dispatching problems based on robust optimization. Among the different ways to contribute to uncertainty sets, the box uncertainty set is over conservative, and the polyhedral uncertainty set cannot express the correlation between uncertain parameters. Therefore, a more flexible ellipsoidal uncertainty set is adopted to describe the uncertainty of PV output. The specific ellipsoidal uncertainty set is as follows
where
In this study, the ARMA model is used to obtain the PV output predicted value
3.5 Solution method
Through the modeling process of Sections 3.1–Sections 3.2, we obtained the active and reactive power coordination dispatching model of the ADN. After the linearization process of the constraints mentioned in Section 3.3, the model proposed in this study is transformed into a mixed integer linear programmed (MILP) optimization model.
Considering the uncertainty of PV output, the time-varying ellipsoidal uncertainty set was developed in Section 3.4. Introducing this set, a robust constraint, into the MILP model, we can obtain a mixed integer robust programmed (MIRP) optimization model that can be computed using solvers. The detailed flowchart of modeling process in this study is presented as Figure 3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1012581/fenrg-10-1012581-HTML-r1/image_m/fenrg-10-1012581-g003.jpg)
FIGURE 3. Flowchart of the robust dispatching model of ADN considering PV time-varying spatial correlation
4 Case studies
4.1 Dynamic spatial correlation model and prediction based on DCC-GARCH
First, the dynamic spatial correlation model of PV output is analyzed, and the data are from the output of PV system in Figure 1 in Suzhou, China. The data of March 1–24 in 2018, during which the weather was sunny, were selected for modeling and prediction. The interval was set to 20 min from 8:00 to 16:00, which can display the changes in the spatial correlation between PVs. After the Augmented Dickey-Fuller (ADF) test, normality, and Lagrange multiplier (LM) test, the results indicate that the historical PV data can be modeled using ARMA and DCC-GARCH, with the ARMA model of PV output being shown in Supplementary Table S1. The parameters of the five-dimensional DCC-GARCH model being shown in Table 2.
The standardized residual series were re-evaluated after the model was constructed. The results indicated no correlation between variables, proving that the model had successfully eliminated the autocorrelation and heteroscedasticity of the data itself. As shown in Table 2, α and β are positive. The sum of α and β (0.974) is less than 1, which indicates that the model is stable, and there is a valid dynamic correlation between the five PVs. Parameter α indicates the influence of the current residual on the time series volatility at the next moment, and a higher α indicates that the time series is more sensitive to the recent residual. Parameter α +β indicates the disappearance speed of current fluctuations in the future, namely, the duration of the spatial correlation of PV output. The higher the value α +β is, the longer is the correlation duration. From Table 2, the result of fitting β is larger than α, indicating that the current dynamic heteroscedasticity of each series arises mainly from the residual of the previous period. The sum of α+β is close to 1, which implies that the spatial correlation has strong continuity.
The DCC matrix Rt of the PV output was calculated using the DCC-GARCH model, with the DCC curve between two pairs under 20 min time resolution being drawn as shown in Figure 4 below.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1012581/fenrg-10-1012581-HTML-r1/image_m/fenrg-10-1012581-g004.jpg)
FIGURE 4. DCC diagrams of PV output. (A) DCC12, DCC13, (B) DCC14, DCC15, (C) DCC23, DCC24, (D) DCC25, DCC34, and (E) DCC35, DCC45.
According to the comparison between Figure 1 and Table 3, it can be observed that distance is one of the factors that affect DCC. The closer the PVs are, the higher is the DCC value. Otherwise, owing to the dynamic change of DCC, the spatial correlation between the two PVs shows different magnitude relationships at different times. This phenomenon is mainly related to temperature and weather.
Comparing Figure 4 with the weather conditions shown in Supplementary Table S2, the maximum and minimum temperatures varied significantly from March 9 to 13. However, the temperature difference remained largely the same owing to five consecutive days of partly cloudy weather. Therefore, the DCC was mainly above 0.9, indicating a high spatial correlation among the PVs during this period. Additionally, the DCC varied dramatically owing to the more pronounced temperature difference fluctuations on March 1–7 and March 5–21. On these days, the weather in the location was cloudy and rainy.
Therefore, the following conclusions can be drawn: When the temperature difference fluctuates more steadily, the DCC changes smoothly, and the DCC is more prominent when the weather is fine. When the temperature difference fluctuates sharply, usually accompanied by inclement weather, the change in solar radiation causes a considerable change in the PV output, which in turn causes a variation in the DCC.
Data from March 1 to 24 were used as the test set, and the forecast step was set to 24. The DCC-GARCH model was used for the rolling forecast to obtain the DCC from 8:00 to 16:00 on March 25 (with an interval of 20 min). The predicted values of the DCC are compared with the actual values in Figure 5, and the mean square error (MSE) is shown in Supplementary Table S3. It can be observed that the errors are within the acceptable range, and the prediction results are accurate. However, in most instances, the prediction values are more prominent than the actual values because of the rapid variation in cloud cover.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1012581/fenrg-10-1012581-HTML-r1/image_m/fenrg-10-1012581-g005.jpg)
FIGURE 5. Comparison of DCC predicted value and actual value. (A) PV1-PV2, PV1-PV3, (B) PV1-PV4, PV1-PV5, (C) PV2-PV3, PV2-PV4, (D) PV2-PV5, PV3-PV4, and (E) PV3-PV5, PV4-PV5.
4.2 Case description
Robust ADN dispatching is verified in the next section. The grid structure is abstracted from a rural distribution network in Anhui, China, as shown in Figure 6. There are five PV stations (different from Figure 1), one ESS, one SVG, and three CBs. The case uses 24 h as the dispatch period and 60 min as the dispatch interval. The capacities of the PV stations were set to 3MW, 2.4MW, 2MW, 2.4MW, and 4 MW. The parameters of the ESS are as follows: The ESS capacity is 1MW; the initial charge state SOCi,0 is 0.5;
To verify the significance of considering the time-varying characteristic of the PV-output spatial correlation, the numerical tests were mainly conducted in the following two cases:
Case A. Robust dispatching model for the ADN without considering the time-varying characteristics of PV-output spatial correlation.
Case B. Robust dispatching model for the ADN with time-varying characteristic consideration of PV-output spatial correlation.
4.3 Dispatching results for case A
When the time-varying characteristic of PV-output spatial correlation is not considered, the correlation coefficients among the five PVs can easily be derived from the Pearson coefficient. In this case, the ellipsoidal uncertainty set is static, as shown in Eq. 64. The objective function and related parameters can be determined as follow:
Figure 7A depicts the PV output curves, and Figure 7B illustrates the active power changes of PCC and the SOC curves of ESS for Case A. It can be observed that when the time-varying spatial correlation is ignored, the PV-output curves are jagged. Some PV-power valleys occur around midday, which leads to the apparent fluctuations of the PCC active power in Figure 7B. To secure the balance of power supply for the entire system, the SOC state of ESS begins charging immediately at 12:00 when the PV output is low and reaches full power at 13:00. The full power condition lasts for nearly 7 h to guarantee that the load is supplied throughout the night so that the objective function is minimized to the greatest extent.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1012581/fenrg-10-1012581-HTML-r1/image_m/fenrg-10-1012581-g007.jpg)
FIGURE 7. Partial-component-change curves for Case A. (A) PV output, (B) Active power changes of PCC and SOC changes of ESS.
The maximum number of OLTC changes within the dispatching period was set to no more than two. Three CBs were connected to nodes 9, 27, and 44. The maximum number of access CB groups was set to five, and the maximum number of changes was two. As shown in Supplementary Figure S1, the variations in OLTC and CB both satisfied the limits. The maximum number of access CB groups during the dispatching period was five. This means that a high proportion of PV in the ADN will cause reactive power shortages and other problems, requiring timely action of reactive-power compensation equipment to ensure stable system operation. In the proposed coordinated dispatching model of ADN, OLTC and CB operate together with the ESS to achieve active management of the ADN and always ensure the voltage quality of the entire system supply.
4.4 Dispatching results for case B
Considering the dynamic trait of PV-output spatial correlation, the covariance matrix Qt, interconnected with the DCC by Eq. 11, can be used in the uncertainty set, where the time-varying ellipsoidal uncertainty set is obtained, as shown in Eq. 65. In addition, because the dispatching interval is 1 h, the PV output time was set as 8:00–16:00, according to the actual situation. The covariances between PVs are shown in Supplementary Figure S2, which indicates the spatial correlation variation. The covariance curves for different PVs share a similar trend: a strong correlation at noon and a weak correlation in the morning and evening, which is caused by the illumination intensity and temperature.
Subsequently, a model for the robust dispatching of ADN that considers time-varying spatial correlation was developed. Considering the influence of different uncertainty budgets on the ellipsoidal uncertainty set, the objective function and related parameters were obtained under different uncertainty budgets, as listed in Table 5, where the PV consumption rate is the ratio of the actual output to the ideal output.
First, the effect of time-varying spatial correlation is considered. Comparing Table 5 with Table 4, without time-varying spatial correlation, it can be observed that the total cost of operation increases owing to the reduction in PV output, which increases the cost of power purchase from the upper grid. On the other hand, the decrease in PV output affects the charging and discharging status of the ESS, thus reducing the ESS benefits.
Figure 8 illustrates the PV output, and Figure 9 reveals the active power changes at the point of common coupling (PCC) at node 1 and the SOC tendencies of the ESS. Compared with Figure 7, under any uncertainty budget, the curves of the PV output and PCC output are smoother, and the value of the PV output at each node is higher when the time-varying spatial correlation is considered. In Figure 8, the PCC output is in the trend of the morning and evening peaks and noon trough. Because most of the loads are residential loads, the peak period of electricity consumption is from 20:00 to 22:00. During this period, the PCC output reaches its peak, and the ESS is in a discharged state to ensure the load supply. The SOC curves indicate that the ESS is releasing energy during the morning peak of electricity consumption after 7 a.m. and charges during the peak of the PV generation. In conclusion, considering the time-varying spatial correlation of PVs is helpful for comprehensively evaluating the output characteristics of PV and optimizing the dispatching results.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1012581/fenrg-10-1012581-HTML-r1/image_m/fenrg-10-1012581-g008.jpg)
FIGURE 8. PV output under different uncertainty budget for Case B. (A) uncertainty budget 99%, (B) uncertainty budget 80%.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1012581/fenrg-10-1012581-HTML-r1/image_m/fenrg-10-1012581-g009.jpg)
FIGURE 9. Active power changes of PCC and SOC changes of ESS for Case B. (A) uncertainty budget 99%, (B) uncertainty budget 80%.
Secondly, we consider the influence of the uncertainty budget. As can be observed from Table 5, as the uncertainty budget increases, the overall operating cost decreases, and the PV consumption rate increases progressively. This is caused by the selection of the uncertainty budget influencing the robustness of the system operation: the greater the uncertainty budget is, the more points are included in the ellipsoidal uncertainty set, and similarly, the more conservative are the results (Wu et al., 2022). Because the total cost of operation only includes the cost of power purchase and energy storage operation, not the cost of PV operation, to minimize the objective function, a high PV output caused by the high uncertainty budget will become the main power supply of the grid, which will reduce the cost of power purchase from the upper grid and the total cost of operation. Under different circumstances, the network loss varies based on the dispatching status of the PV output and the upper-grid supply.
Under uncertainty budgets of 99% and 80%, Figures 8, 9 shows quite different curves. When the uncertain budget is 99%, the output value of the PV at each moment is slightly larger than that under the 80% budget, and the output curve of the PCC is smoother. Because the total operating cost considers the lifespan loss cost of the ESS, to minimize the objective function as much as possible, the charging and discharging times of the ESS are reduced under an uncertainty budget of 99% with higher conservatism.
Because the case is a rural distribution network, a low load demand is accompanied by high PV penetration. In the high PV output period, combining the peak–valley spread revenue of the ESS, surplus power will be sold to the upper grid at the purchase price. This approach can improve the PV consumption rate and result in PCC power for negative values from 8:00 to 16:00. In the entire dispatching process with an uncertainty budget of 99%, the ESS benefit is 649.108 yuan, which appropriately reduces the total operation cost.
In the daytime, the voltage of each node connected to the PV system increases significantly, as shown in Supplementary Figure S3. Owing to the PV power cut-in at 8:00 and cut-out at 16:00, all five node voltages exhibited considerable fluctuations. Additionally, more PVs with larger capacities were connected to Line I, resulting in a higher overall node voltage than that of Line III.
5 Conclusion
In this study, a robust dispatching model for an ADN was proposed, which considers PV-output spatial correlation and uncertainty. A dynamic spatial correlation model was proposed to characterize the time-varying characteristics of PV output spatial correlation. Based on this, a time-varying ellipsoidal uncertain set was constructed and applied to the ADN dispatching model, which improved the robustness of PV-output model. Case studies on a rural distribution network in China demonstrated the validity and rationality of the proposed methods. The results indicated that:
(1) The DCC derived using the DCC-GARCH model correctly characterized the spatial correlation of PV output and reflected the time-varying properties of the spatial correlation.
(2) The value of DCC is related to the weather conditions and temperature difference in the area where the PV is located: When the temperature difference swings more consistently, the DCC changes more smoothly; When the weather is favorable, the DCC is greater. When the temperature difference varies dramatically, which is generally accompanied by severe weather, the DCC also changes.
(3) The time-varying ellipsoidal uncertainty set constructed using the DCC can be well applied to the optimal dispatching model of ADN. Considering the time-varying characteristics of the spatial correlation, it is advantageous to build a complete PV output model and optimize the dispatching results.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
XM: Methodology, analysis, and writing. HW: Conceptualization, software, review, and editing. YY: supervision, review, and editing.
Funding
This work is supported by the Major Basic Research Project of the Natural Science Foundation of the Jiangsu Higher Education Institutions of China (22KJD470003).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fenrg.2022.1012581/full#supplementary-material
References
Aghamohamadi, M., Mahmoudi, A., and Haque, M. H. (2021). Two-stage robust sizing and operation co-optimization for residential PV–battery systems considering the uncertainty of PV generation and load. IEEE Trans. Ind. Inf. 17, 1005–1017. doi:10.1109/TII.2020.2990682
Aharon, B.-T., and Laurent El, G. (2009). Robust optimization. Princeton: Princeton University Press.
Antoniadou-Plytaria, K. E., Kouveliotis-Lysikatos, I. N., Georgilakis, P. S., and Hatziargyriou, N. D. (2017). Distributed and decentralized voltage control of smart distribution networks: Models, methods, and future research. IEEE Trans. Smart Grid 8, 2999–3008. doi:10.1109/TSG.2017.2679238
Calcabrini, A., Ziar, H., Isabella, O., and Zeman, M. (2019). A simplified skyline-based method for estimating the annual solar energy potential in urban environments. Nat. Energy 4, 206–215. doi:10.1038/s41560-018-0318-6
Chassein, A., and Goerigk, M. (2016). Min-max regret problems with ellipsoidal uncertainty sets. arXiv: Optimization and Control. doi:10.48550/arXiv.1606.01180
Choi, J., Lee, J.-I., Lee, I.-W., and Cha, S.-W. (2022). Robust PV-BESS scheduling for a grid with incentive for forecast accuracy. IEEE Trans. Sustain. Energy 13, 567–578. doi:10.1109/TSTE.2021.3120451
Ding, F., and Mather, B. (2017). On distributed PV hosting capacity estimation, sensitivity study, and improvement. IEEE Trans. Sustain. Energy 8, 1010–1020. doi:10.1109/TSTE.2016.2640239
El-Meligy, M. A., El-Sherbeeny, A. M., and Anvari-Moghaddam, A. (2022). Transmission expansion planning considering resistance variations of overhead lines using minimum-volume covering ellipsoid. IEEE Trans. Power Syst. 37, 1916–1926. doi:10.1109/TPWRS.2021.3110738
Engle, R., and Sheppard, K. (2001). Theoretical and empirical properties of dynamic conditional correlation multivariate GARCH. Cambridge, MA: National Bureau of Economic Research.
Golestaneh, F., Pinson, P., Azizipanah-Abarghooee, R., and Gooi, H. B. (2018). Ellipsoidal prediction regions for multivariate uncertainty characterization. IEEE Trans. Power Syst. 33, 4519–4530. doi:10.1109/TPWRS.2018.2791975
Haque, M. M., and Wolfs, P. (2016). A review of high PV penetrations in LV distribution networks: present status, impacts and mitigation measures. Renew. Sustain. Energy Rev. 62, 1195–1208. doi:10.1016/j.rser.2016.04.025
Hu, R., Wang, W., Wu, X., Chen, Z., Jing, L., Ma, W., et al. (2022). Coordinated active and reactive power control for distribution networks with high penetrations of photovoltaic systems. Sol. Energy 231, 809–827. doi:10.1016/j.solener.2021.12.025
Ji, H., Wang, C., Li, P., Ding, F., and Wu, J. (2019). Robust operation of soft open points in active distribution networks with high penetration of photovoltaic integration. IEEE Trans. Sustain. Energy 10, 280–289. doi:10.1109/TSTE.2018.2833545
Liu, H., Zhang, Y., Ge, S., Gu, C., and Li, F. (2019). Day-ahead scheduling for an electric vehicle PV-based battery swapping station considering the dual uncertainties. IEEE Access 7, 115625–115636. doi:10.1109/ACCESS.2019.2935774
Luo, Y., Yang, D., Yin, Z., Zhou, B., and Sun, Q. (2020). Optimal configuration of hybrid-energy microgrid considering the correlation and randomness of the wind power and photovoltaic power. IET Renew. Power Gener. 14, 616–627. doi:10.1049/iet-rpg.2019.0752
Pan, C., Wang, C., Zhao, Z., Wang, J., and Bie, Z. (2019). A copula function based monte Carlo simulation method of multivariate wind speed and PV power spatio-temporal series. Energy Procedia 159, 213–218. doi:10.1016/j.egypro.2018.12.053
Sun, X., Qiu, J., Tao, Y., Ma, Y., and Zhao, J. (2022). Coordinated real-time voltage control in active distribution networks: An incentive-based fairness approach. IEEE Trans. Smart Grid 13, 2650–2663. doi:10.1109/TSG.2022.3162909
Vilaça Gomes, P., Saraiva, J. T., Carvalho, L., Dias, B., and Oliveira, L. W. (2019). Impact of decision-making models in transmission expansion planning considering large shares of renewable energy sources. Electr. Power Syst. Res. 174, 105852. doi:10.1016/j.epsr.2019.04.030
Wang, Y., Zhou, Z., Botterud, A., Zhang, K., and Ding, Q. (2016). Stochastic coordinated operation of wind and battery energy storage system considering battery degradation. J. Mod. Power Syst. Clean. Energy 4, 581–592. doi:10.1007/s40565-016-0238-z
Wu, H., Yuan, Y., Zhang, X., Miao, A., and Zhu, J. (2022). Robust comprehensive PV hosting capacity assessment model for active distribution networks with spatiotemporal correlation. Appl. Energy 323, 119558. doi:10.1016/j.apenergy.2022.119558
Wu, H., Yuan, Y., Zhu, J., Qian, K., and Xu, Y. (2021). Potential assessment of spatial correlation to improve maximum distributed PV hosting capacity of distribution networks. J. Mod. Power Syst. Clean Energy 9, 800–810. doi:10.35833/MPCE.2020.000886
Xu, S., Sun, H., Zhao, B., Yi, J., Hayat, T., Alsaedi, A., et al. (2020). The integrated design of a novel secondary control and robust optimal energy management for photovoltaic-storage system considering generation uncertainty. Electronics 9, 69. doi:10.3390/electronics9010069
Yang, Z., Bose, A., Zhong, H., Zhang, N., Xia, Q., and Kang, C. (2017). Optimal reactive power dispatch with accurately modeled discrete control devices: A successive linear approximation approach. IEEE Trans. Power Syst. 32, 2435–2444. doi:10.1109/TPWRS.2016.2608178
Yang, Z., Zhong, H., Xia, Q., Bose, A., and Kang, C. (2016). Optimal power flow based on successive linear approximation of power flow equations. IET Gener. Transm. Distrib. 10, 3654–3662. doi:10.1049/iet-gtd.2016.0547
Yu, Y., Konstantinou, G., Hredzak, B., and Agelidis, V. G. (2015). Operation of cascaded H-Bridge multilevel converters for large-scale photovoltaic power plants under bridge failures. IEEE Trans. Ind. Electron. 62, 7228–7236. doi:10.1109/TIE.2015.2434995
Zamee, M. A., and Won, D. (2020). Novel mode adaptive artificial neural network for dynamic learning: Application in renewable energy sources power generation prediction. Energies 13, 6405. doi:10.3390/en13236405
Nomenclature
A. Sets
BNode set of all nodes
T set of dispatching time
BOLTC set of OLTC nodes
BESS set of ESS nodes
BCB set of CB nodes
BSVG set of SVG nodes
BPV set of PV nodes
B. Parameters
Sij,max maximum capacity of branch ij
SOCi,0, SOCi,set, SOCi,end the initial value, the set initial value, and the end value of SOC during the dispatch cycle
C. Variables
Keywords: time-varying spatial correlation, DCC-GARCH, correlation prediction, robust dispatching, time-varying ellipsoidal uncertainty set
Citation: Ma X, Wu H and Yuan Y (2023) Robust dispatching model of active distribution network considering PV time-varying spatial correlation. Front. Energy Res. 10:1012581. doi: 10.3389/fenrg.2022.1012581
Received: 05 August 2022; Accepted: 02 September 2022;
Published: 09 January 2023.
Edited by:
Yuchen Zhang, University of New South Wales, AustraliaCopyright © 2023 Ma, Wu and Yuan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Han Wu, wuhanichina@vip.qq.com