Integration Capability Evaluation of Wind and Photovoltaic Generation in Power Systems Based on Temporal and Spatial Correlations

With the rapid growth of renewable energy generation, it has become essential to give a comprehensive evaluation of renewable energy integration capability in power systems to reduce renewable generation curtailment. Existing research has not considered the correlations between wind power and photovoltaic (PV) power. In this paper, temporal and spatial correlations among different renewable generations are utilized to evaluate the integration capability of power systems based on the copula model. Firstly, the temporal and spatial correlation between wind and PV power generation is analyzed. Secondly, the temporal and spatial distribution model of both wind and PV power generation output is formulated based on the copula model. Thirdly, aggregated generation output scenarios of wind and PV power are generated. Fourthly, wind and PV power scenarios are utilized in an optimal power flow calculation model of power systems. Lastly, the integration capacity of wind power and PV power is shown to be able to be evaluated by satisfying the reliability of power system operation. Simulation results of a modified IEEE RTS-24 bus system indicate that the integration capability of renewable energy generation in power systems can be comprehensively evaluated based on the temporal and spatial correlations of renewable energy generation.


Introduction
To cope with the impending fossil energy crisis, environmental pollution, and the greenhouse effect, renewable energy is being developed rapidly, especially with regards to power generation.Given the superiorities of clean energy and its low marginal price, renewable energy generation, including wind power and photovoltaic (PV) power, has occupied a larger and larger proportion of power system generation in recent years.While the global cumulative installed capacity of wind power in 2008 amounted to only 120.8 GW, the following decade saw a rapid growth rate and the number was predicted to reach 592 GW in 2018 [1].There is also currently a distinct global development in PV generation.With a total installed capacity of 15.84 GW in 2008, it was forecast to total 501.6 GW in 2018 [2].
However, the generation outputs of wind turbines and PV panels depend heavily on the weather, which leads to the nondeterministic nature of wind power and PV power.For example, wind power generation relies largely on the wind speed and wind direction [3], while solar irradiance and ambient temperature are the key factors affecting PV power [4].Due to the inherent uncertainty of wind and PV power, it is difficult to make full use of these renewable generations.In addition, most existing power system networks were designed and built years ago without enough foresight to imagine such rapid growth in renewable generation [5,6].Moreover, most renewable generation farms are built at the end of power grids.Existing factors, such as transmission capacities and reserve capacities, constrain the capability of power systems from fully adopting renewable generation, thereby resulting in severe wind power and PV power curtailments [7].Therefore, it is essential to evaluate the integration capability of renewable energy generation in power systems, which can be utilized in power system planning for updating networks or building new renewable generation farms.
Much research has been done on evaluating the integration capacity of either wind power or PV power generation in power systems [8].The impact of increasing PV penetration on the static performance and transient stability of transmission systems was analyzed in [9].In [10], optimal levels of plug-in electric vehicles and renewable distributed generation penetration with respect to the location, size, and year of installation were optimized based on a non-dominated sorting genetic algorithm.Technologies for accommodating renewable resources penetration have been proposed in cooperation with energy storage systems and grid flexibility [11,12].Distribution system feeders have also been investigated to cope with the possible voltage problems related to the integration sizes and locations of PV generation and loads [13].The flexibility requirements for wind and solar power in Europe were quantitatively evaluated in terms of ramping capability and geography in [14].Voltage source inverters control strategies in renewable energy integration in power systems have also been widely studied [15].A comprehensive review for analysis of renewable energy integration into multi-energy systems was presented in [16].However, most existing research either focuses on penetration capability analysis with individual renewable resources or combinations of multiple renewable resources without considering their correlations in time and space.It is possible for outputs of wind generation and PV generation to have large temporal and spatial correlations because of their shared dependence on weather forecasts.This correlation may distort the aggregated outputs of wind power and PV power compared with their direct summation results.The correlation between wind and PV power has been explored in a number of studies.For example, the sample correlation coefficient between wind and PV power was calculated in [17].In [18], a Monte Carlo based method was adopted for the assessment of wind and PV power correlations.However, the joint probability distribution of wind and solar power has not been formulated in these studies.Without distribution estimations, it is difficult to take wind and PV power correlations into account in evaluating renewable generation penetration in power system planning.
In this paper, a novel approach for evaluating the integration capability of renewable energy generation in transmission systems is proposed based on temporal and spatial correlations.Firstly, the correlations between wind and PV power series in different time periods and locations are analyzed based on the cross-correlation function (CCF).Secondly, an aggregated generation output distribution model of wind and PV power is formulated from individual predictive distributions based on the copula function.In addition, aggregated generation outputs are segmented based on different correlation degrees of wind and PV power series.Thirdly, an optimization model with the objective of maximizing the loss of load probability (LOLP) index is established based on the aggregated renewable generation distribution corresponding to the varying installed capacities of wind farms and PV power farms.Lastly, the optimal installed capacities of wind power and PV power is shown to be able to be achieved by solving the optimization model.

Correlation Analysis
A correlation analysis is critical for estimating the probability of renewable energy power output and determining a reasonable strategy of renewable energy generation aggregation.The statistical properties of wind power and PV power series are explained in this section.
To intuitively present the correlations between wind power and PV power, the different time series of normalized wind power generation and PV power generation are shown in Figure 1.The wind farm and PV farm shown in this figure are from the same area.It is clearly shown in this figure that wind power and PV power are correlated in some degree when in the same area.To intuitively present the correlations between wind power and PV power, the different time series of normalized wind power generation and PV power generation are shown in Figure 1.The wind farm and PV farm shown in this figure are from the same area.It is clearly shown in this figure that wind power and PV power are correlated in some degree when in the same area.To give a mathematical description of the correlations between wind power and PV power, we adopt the Auto-Correlation Function (ACF) and the Cross-Correlation Function (CCF) [19] to explore the temporal and spatial dependence of wind power and PV power.Given a set of sampled wind power and PV power series, the ACF of wind power and PV power series can be defined as where a is the lag time between the two series; t g is the power generation of the PV or wind farm at time t; and μ and 2 σ are the mean value and variance of the time series, respectively.
The CCF between two times series is defined as where , , pv i t g and , , wind j t a g − are the power of sampled PV series i and sampled wind power series j at time t, respectively; a is the lag time between the two series; i μ and i σ are the mean value and variance of sampled PV series i , respectively; and j μ and j σ are the mean value and variance of sampled PV series j , respectively.
The ACF of the PV and wind power generation, as well as the CCF values between the series of the PV farm and wind farm are presented in Figure 2. Figure 2a,d show that PV power and wind power have strong temporal dependence while PV power also has strong auto-correlation at the same time for neighboring days.Figure 2b,c show relatively large CCF values, which verifies the strong spatial dependence relation between the generation output series of wind power and PV power in the same area.However, as shown in Figure 1, PV power output at night is zero, which may reduce the cross-correlation between PV and wind power.For this reason, we tested the CCF values between PV and wind power by removing the time segments when PV output was zero in both the PV and wind power series.The results are shown in Figure 3, which indicates that the spatial dependence between PV and wind farms is more obvious.Because of this, the aggregated distribution of the output of PV farms and wind farms has only been modeled for day time periods.To give a mathematical description of the correlations between wind power and PV power, we adopt the Auto-Correlation Function (ACF) and the Cross-Correlation Function (CCF) [19] to explore the temporal and spatial dependence of wind power and PV power.Given a set of sampled wind power and PV power series, the ACF of wind power and PV power series can be defined as where a is the lag time between the two series; g t is the power generation of the PV or wind farm at time t; and µ and σ 2 are the mean value and variance of the time series, respectively.The CCF between two times series is defined as where g pv,i,t and g wind,j,t−a are the power of sampled PV series i and sampled wind power series j at time t, respectively; a is the lag time between the two series; µ i and σ i are the mean value and variance of sampled PV series i, respectively; and µ j and σ j are the mean value and variance of sampled PV series j, respectively.The ACF of the PV and wind power generation, as well as the CCF values between the series of the PV farm and wind farm are presented in Figure 2. Figure 2a,d show that PV power and wind power have strong temporal dependence while PV power also has strong auto-correlation at the same time for neighboring days.Figure 2b,c show relatively large CCF values, which verifies the strong spatial dependence relation between the generation output series of wind power and PV power in the same area.However, as shown in Figure 1, PV power output at night is zero, which may reduce the cross-correlation between PV and wind power.For this reason, we tested the CCF values between PV and wind power by removing the time segments when PV output was zero in both the PV and wind power series.The results are shown in Figure 3, which indicates that the spatial dependence between PV and wind farms is more obvious.Because of this, the aggregated distribution of the output of PV farms and wind farms has only been modeled for day time periods.The spot forecasts of PV power and wind power were achieved based on the Sparse Bayesian Learning (SBL) model.SBL is a parametric forecasting method that assumes wind power and PV power output follow Gaussian distributions.Although the statistical distributions of power generation time series are non-Gaussian, some research has verified its effectiveness in modeling the uncertainty characteristics of renewable generation [20].In this paper, we only utilize the SBL model to forecast the expectations of wind power and PV power.It should be noted that any forecasting  The spot forecasts of PV power and wind power were achieved based on the Sparse Bayesian Learning (SBL) model.SBL is a parametric forecasting method that assumes wind power and PV power output follow Gaussian distributions.Although the statistical distributions of power generation time series are non-Gaussian, some research has verified its effectiveness in modeling the uncertainty characteristics of renewable generation [20].In this paper, we only utilize the SBL model to forecast the expectations of wind power and PV power.It should be noted that any forecasting The spot forecasts of PV power and wind power were achieved based on the Sparse Bayesian Learning (SBL) model.SBL is a parametric forecasting method that assumes wind power and PV power output follow Gaussian distributions.Although the statistical distributions of power generation time series are non-Gaussian, some research has verified its effectiveness in modeling the uncertainty characteristics of renewable generation [20].In this paper, we only utilize the SBL model to forecast the expectations of wind power and PV power.It should be noted that any forecasting method which has good forecasting accuracy can be utilized to obtain spot forecast results of wind and PV power in this model.Thus, the forecast error series of PV power and wind power can be calculated by a comparison with the real measurements of PV power and wind power.
In addition, cross-correlation between the PV forecast error and wind forecast error series was tested.The results show significant cross-correlation relationships.Therefore, the PV forecast error and wind forecast error series have been utilized in the power aggregation model.

Basics of the Copula Function
The copula function is a function widely utilized for calculating the joint distributions of random variables given all marginal distributions [21].A bivariate copula function can be described as F e pv , e wind = C F pv e pv , F wind (e wind ) (3) where C(•) represents the copula function; F pv (•) and F wind (•) are the marginal cumulative distribution functions (CDFs) of PV and wind power, respectively; e pv and e wind stand for the random variables of PV and wind power forecast errors, respectively; and F(•) is the joint CDF of the random variables, that is, in this case it is the joint CDF of PV and wind power.
Based on Sklar's theorem [22], a unique copula function C(•) can be calculated from its continuous marginal CDFs.Therefore, in this paper, the dependence structure of wind power and PV power forecast errors may be modeled based on an appropriate copula function.

Copula-Based Renewable Generation Aggregation
A Gaussian copula function has been utilized in this paper.As copula functions construct the joint CDF only for uniform random variables, it is necessary to transform the marginal CDFs of wind and PV power into uniform distributions.It is indicated in [23] that for uniform random variables with a symmetrical dependence structure, a Gaussian copula function can estimate their joint CDF with high accuracy.For this reason, a scatter plot was utilized to explore the dependence structure of wind and PV power after uniform transformation.Figure 4 presents the scatter plot results of the pairs of PV and wind power forecast errors, which were obtained by calculating the difference between the sampled data and the results of the SBL prediction model mentioned in Section 2.1.In Figure 4a, the wind power and PV power forecast error series are plotted.In Figure 4b, the scatter plot of the transformed uniform random variables of wind power and PV power forecast error series are displayed.The transformation was implemented by probability integral transform rules [24].With regards to the scatter plot results, the clearly symmetrical dependence structure demonstrates that the dependence relation between PV and wind power forecast error series can be modeled by a Gaussian copula function.
Because of this, the dependence relationship between the forecast error series of wind power and PV power was modeled by the two-dimensional Gaussian copula (as defined by [23]) where Φ represents a univariate standard Gaussian CDF with Φ −1 as its inverse function; and Φ R stands for a bivariate Gaussian joint CDF with zero means and unit marginal variances whose covariance matrix is R.
To obtain the joint CDF described in (4), the marginal CDFs of the forecast error series of wind power and PV power, as well as the corresponding correlation matrix R, must be calculated.In this paper, the marginal CDFs of the forecast error series were obtained based on the empirical CDF.R can be determined based on the maximum likelihood estimation [25] ln c(F pv e pv,i , F wind (e wind,i ), R) (5) where c(•) is the joint density function corresponding to C(•); V is the number of the sampled pairs of PV and wind power forecast errors and e pv,i , e wind,i is the i th pair.It can be seen that R calculated by ( 5) maximizes the probability of the occurrence of the sampled data, which is a reasonable estimation of the parameters of the Gaussian joint CDF.According to the spot forecast results, a CDF of the spot forecast error, and the correlation matrix of the errors, V groups of the formulated two-dimensional scenarios of the aggregated output of wind power and PV power can be generated.The process is as follows: • Step 1: Sample V groups of the two-dimensional uniform random variables PV and wind power forecast errors based on the joint CDF estimated by the Gaussian copula model in (4).The groups of these samples can be represented by Step 2: Transform the above V groups of sample u i = u pv,i , u wind,i , i = 1, • • • V to V groups of error samples e i = e pv,i , e wind,i , i = 1, • • • V through the inverse transformation e pv,i , e wind,i = F −1 pv u pv,i , F −1 wind (u wind,i ) where F −1 pv (•) and F −1 wind (•) are the inverse functions of F pv (•) and F wind (•), respectively.

•
Step 3: Generate V groups of the two-dimensional scenarios of PV and wind power generation.
The jth scenario s j = P PV j , P W j can be achieved by where p i is the forecasted PV and wind power generated by the SBL prediction model.
The V groups of the two-dimensional scenarios of PV and wind power generation are respectively normalized values based on the maximum PV and wind power generation output.Therefore, it is possible to use the scenarios to simulate PV and wind power generation scenarios by multiplying different integration capacities of PV farms and wind farms.

Optimization Model for Exploiting the Integration Capacity of Wind Power and PV Power
With increasing integration capacity of wind power and PV power, possible variations of renewable generation output will be enlarged due to deviations between predictive results and real observations.To cope with the increasing uncertainty of renewable generation, the integration capacity of wind and PV power must satisfy system reliability requirements.Thus, the objective of the proposed evaluation model was to obtain the maximum wind power and PV power integration capacities considering the required system operating constraints and reliability requirements.By evaluating the system reliability corresponding to the integrations of wind power and PV power, it is possible to achieve the maximum integration capacity of wind farms and PV farms.In this reliability evaluation model, the N-1 contingency scenarios were considered.For each contingency scenario, the simulated V groups of wind power and PV power scenarios were utilized to represent the uncertainty of renewable generation.Therefore, the number K of the system scenarios was where N g is the number of generators; N l is the number of branches; and N pd is the number of loads.Thus, the above optimization model was formulated as Obj.: min where ∆LC k,i is load curtailment at bus i in scenario k; and B is the set of system buses, such that Energies 2019, 12, 171 8 of 12 where P G k,g represents the scheduled generation of unit g in scenario k; G i , S i and W i are the set of generators, PV farms, and wind farms that are connected with bus i, respectively; P L k,i is the load at bus i in scenario k; P G,max g and P G,min g represent the maximum and minimum generation of unit g; P PV k,s and P W k,w indicate the simulated PV and wind power in scenario k; ∆LC k,i is the load curtailment at bus i in scenario k; P ij and P ij are, respectively, the lower and upper bounds of power flow in transmission line from bus i to bus j; P k,ij is the power flow in transmission line from bus i to bus j; V k,i and δ k,i are, respectively, the amplitude and angle at bus i; and G k,ij and B k,ij are the admittances.
The constraints were constructed based on simulated wind power and PV power scenarios.The system balance constraint is presented in (10).The transmission line capacity limit is presented in (11).Equations ( 12) and ( 13) describe the amplitude and phase angle constraints.The generation capacity constraints of generating units are given in (14).The installed capacities of wind power and PV power should be within the limits of pre-set maximum capacities as shown in ( 15) and ( 16), respectively.
The reliability represented by LOLP should be within the required limit, as shown in (17).
where 1(right) = 1 and 1(wrong) = 0; and p k is the probability of scenario k.Here, we assume p k = 1 K in simulations.It should be noted that p k can be achieved from experience values from other papers or actual practices in real system operation.
By solving the optimization model as given in ( 9)~( 17), the LOLP indexes corresponding to different installed capacities of wind farms and PV farms can be achieved.Based on estimated LOLP values, the optimal wind power and PV power integration capacity can be analyzed.

Overall Simulation Frame Work of the Proposed Model
The pseudocodes for the overall simulation process are shown in Algorithm 1. Firstly, the SBL model was adopted to generate wind and PV power spot forecasting results based on sampled historical data.Following this, the joint probability distribution of wind and PV power was formulated by a Gaussian copula function based on the spot forecasting results.Then, the aggregated wind power and PV power scenarios could be generated by the joint probability distribution.Lastly, the reliability of power systems was evaluated by solving the optimization model constructed in Section 3.

Algorithm 1 Pseudocodes for the overall simulation process
Function: spot forecast (PV power and wind power) Model: SBL model; Inputs: historical wind power/PV power measurements in the last 4 h at t − 1, t − 2, t − 3 and t − 4; Outputs: wind power/PV power spot forecasting results at time t.
Function: copula-based distribution of PV power and wind power Model: Gaussian copula; Inputs: wind power/PV power spot forecasting results; forecast error series of wind power and PV power; correlation matrix; Outputs: aggregated wind power and PV power scenarios.
Function: reliability evaluation model Model: shown in equations in ( 8)~( 17); Inputs: varying integration capacities of wind and PV farms; aggregated wind power and PV power scenarios; data of RTS-24 system; Outputs: LOLP values.

Case Study
The modified IEEE RTS-24 bus system was utilized to verify the effectiveness of the proposed model.The parameters of the system are described in detail in [26].In this paper, a wind farm and a PV farm were connected to bus 7 and bus 18 to replace two generators connected to the two buses, respectively.Real measurements of the outputs of a wind farm and a PV farm in the same area of Zhejiang Province were utilized to formulate the aggregated distribution model of wind power and PV power.

Reliability Evaluation
The system reliability represented by indices LOLP respect to different installed capacities of wind power and PV power are presented in Figure 5a.In addition, the LOLP variations corresponding to individual renewable generation resources are shown in Figures 5b and 5c, respectively.Figure 5b,c were obtained by projecting the curved surface in Figure 5a onto the "wind power-LOLP" space and "PV power-LOLP" space respectively.It should be noted that the total installed generation capacity was regarded as a constant.Consequently, as the installed capacity of the renewable farms increased, the capacity of traditional generators decreased accordingly.Based on the proportional allocation rule of capacity, the decreased generation amount was shared by all traditional generators.Clearly, the LOLP increased with the growth of the installed capacities of wind power and PV power.Via consideration of this figure, we can obtain the optimal installed renewable generation capacity respect to the acceptable reliability level in power system operation.
To indicate the effect of considering the correlation of wind and PV power, we also evaluated the system reliability without taking into account the correlation of wind and PV power for comparison's sake.The evaluation result is presented in Figure 6.It can be seen that with the same installed capacities of wind and PV power, LOLP is smaller compared to that considering correlation.This result indicates that the correlation of wind and PV power has an impact on the reliability result, the ignoring of which leads to a relatively optimistic estimate for the system reliability and therefore more load curtailment in the actual power system operation.To indicate the effect of considering the correlation of wind and PV power, we also evaluated the system reliability without taking into account the correlation of wind and PV power for comparison's sake.The evaluation result is presented in Figure 6.It can be seen that with the same installed capacities of wind and PV power, LOLP is smaller compared to that considering correlation.This result indicates that the correlation of wind and PV power has an impact on the reliability result, the ignoring of which leads to a relatively optimistic estimate for the system reliability and therefore more load curtailment in the actual power system operation.To indicate the effect of considering the correlation of wind and PV power, we also evaluated the system reliability without taking into account the correlation of wind and PV power for comparison's sake.The evaluation result is presented in Figure 6.It can be seen that with the same installed capacities of wind and PV power, LOLP is smaller compared to that considering correlation.This result indicates that the correlation of wind and PV power has an impact on the reliability result, the ignoring of which leads to a relatively optimistic estimate for the system reliability and therefore more load curtailment in the actual power system operation.

Voltage Variation
The variation of the system voltage amplitude corresponding to the installed capacities of wind power and PV power is shown in Figure 7.The voltage amplitude with respect to a specific pair of wind and PV power is the voltage of the specific bus whose voltage deviated most from the reference value of 1 in p.u. in all normal scenarios without component failure.In China, the normalized bus voltage amplitude is required to be within 0.95-1.05 in a normal operation state of a power system.As shown in Figure 7, the blue box in three-dimensional space is the acceptable voltage amplitude corresponding to the optimal installed capacities of wind power and PV power.The red circles with a deep gray color are the optimal simulation results.With the proposed model, we can identify the optimal installed capacities of wind power and PV power considering the possible bus voltage variations of the power system.

Voltage Variation
The variation of the system voltage amplitude corresponding to the installed capacities of wind power and PV power is shown in Figure 7.The voltage amplitude with respect to a specific pair of wind and PV power is the voltage of the specific bus whose voltage deviated most from the reference value of 1 in p.u. in all normal scenarios without component failure.In China, the normalized bus voltage amplitude is required to be within 0.95-1.05 in a normal operation state of a power system.As shown in Figure 7, the blue box in three-dimensional space is the acceptable voltage amplitude corresponding to the optimal installed capacities of wind power and PV power.The red circles with a deep gray color are the optimal simulation results.With the proposed model, we can identify the optimal installed capacities of wind power and PV power considering the possible bus voltage variations of the power system.

Conclusions
In this study, a novel integration capability evaluation model based on temporal and spatial correlations among different renewable generations have been proposed using a copula model.The aggregation distribution model for wind power generation and PV power generation was formulated based on their temporal and spatial correlations.The integration capacities of wind power and PV power were evaluated based on the reliability analysis of power systems.Two-dimensional evaluation results were produced based on the temporal and spatial model.The optimal integration capacities of wind power and PV power in a power grid can be determined according to the requirements of system reliability.The simulation results demonstrate that the proposed model can give a comprehensive evaluation of the renewable generation integration ability of power systems.
Author Contributions: H.Z., H.W., C.Y., S.X., J.Z., X.H., and B.W. proposed and developed the basic models and methods, analyzed the data and wrote the paper together.

Conclusions
In this study, a novel integration capability evaluation model based on temporal and spatial correlations among different renewable generations have been proposed using a copula model.The aggregation distribution model for wind power generation and PV power generation was formulated based on their temporal and spatial correlations.The integration capacities of wind power and PV power were evaluated based on the reliability analysis of power systems.Two-dimensional evaluation results were produced based on the temporal and spatial model.The optimal integration capacities of wind power and PV power in a power grid can be determined according to the requirements of system reliability.The simulation results demonstrate that the proposed model can give a comprehensive evaluation of the renewable generation integration ability of power systems.

Figure 1 .
Figure 1.Photovoltaic (PV) power and wind power time series.

Figure 1 .
Figure 1.Photovoltaic (PV) power and wind power time series.

Figure 2 .
Figure 2. The Auto-Correlation Function (ACF) of PV power generation and wind power generation and Cross-Correlation Function (CCF) values between the series of PV farms and wind farms.

Figure 3 .
Figure 3.The ACF of PV power generation and wind power generation and CCF values between the series of PV farms and wind farms without zero PV output time.

Figure 2 .
Figure 2. The Auto-Correlation Function (ACF) of PV power generation and wind power generation and Cross-Correlation Function (CCF) values between the series of PV farms and wind farms.

Figure 2 .Figure 3 .
Figure 2. The Auto-Correlation Function (ACF) of PV power generation and wind power generation and Cross-Correlation Function (CCF) values between the series of PV farms and wind farms.

Figure 3 .
Figure 3.The ACF of PV power generation and wind power generation and CCF values between the series of PV farms and wind farms without zero PV output time.

Figure 4 .
Figure 4. Scatter plot results of the dependence relationship between the random variables.(a) Scatter plot between the forecast errors of wind power and PV power.(b) Scatter plot between the transformed forecast errors of wind power and PV power.

Figure 4 .
Figure 4. Scatter plot results of the dependence relationship between the random variables.(a) Scatter plot between the forecast errors of wind power and PV power.(b) Scatter plot between the transformed forecast errors of wind power and PV power.

Figure 5 .
Figure 5. Loss of load probability (LOLP) variation with respect to installed capacities of wind power and PV power, considering correlations.(a) LOLP with respect to wind power and PV power; (b) LOLP with respect to wind power; (c) LOLP with respect to PV power.

Figure 5 .
Figure 5. Loss of load probability (LOLP) variation with respect to installed capacities of wind power and PV power, considering correlations.(a) LOLP with respect to wind power and PV power; (b) LOLP with respect to wind power; (c) LOLP with respect to PV power.

Figure 5 .
Figure 5. Loss of load probability (LOLP) variation with respect to installed capacities of wind power and PV power, considering correlations.(a) LOLP with respect to wind power and PV power; (b) LOLP with respect to wind power; (c) LOLP with respect to PV power.

Figure 6 .
Figure 6.LOLP variation with respect to installed capacities of wind power and PV power, without considering correlations.(a) LOLP with respect to wind power and PV power; (b) LOLP with respect to wind power; (c) LOLP with respect to PV power.

Figure 7 .
Figure 7. Voltage amplitude variations corresponding to the varying installed capacities of wind power and PV power.

Figure 7 .
Figure 7. Voltage amplitude variations corresponding to the varying installed capacities of wind power and PV power.