Developing a Mathematical Model for Wind Power Plant Siting and Sizing in Distribution Networks

Wind Power Plants (WPPs) are generally located in remote areas with weak distribution connections. Hence, the value of Short Circuit Capacity (SCC), WPP size and the short circuit impedance angle ratio (X/R) are all very critical in the voltage stability of a distribution system connected WPP. This paper presents a new voltage stability model based on the mathematical relations between voltage, the level of wind power penetration, SCC and X/R at a given Point of Common Coupling (PCC) of a distribution network connected WPP. The proposed model introduces six equations based on the SCC and X/R values seen from a particular PCC point. The equations were developed for two common types of Wind Turbine Generators (WTGs), including: the Induction Generator (IG) and the Double Fed Induction Generator (DFIG). Taking advantage of the proposed equations, design engineers can predict how the steady-state PCC voltage will behave in response to different penetrations of IGand DFIG-based WPPs. In addition, the proposed equations enable computing the maximum size of the WPP, ensuring grid code requirements at the given PCC without the need to carry out complex and time-consuming computational tasks or modelling of the system, which is a significant advantage over existing WPP sizing approaches.


Introduction
Wind power is one of the most abundant, sustainable, cost effective and clean fuel energy sources [1,2]. The majority of small wind farms are being connected to distribution networks [3]. Increasing the penetration of Wind Power Plants (WPPs) in a distribution grid is subject to providing the voltage stability requirements defined by the grid codes. For example, the Australian and UK grid code requires continuous control of steady-state voltage at the Point of Common Coupling (V PCC ) of a distribution network connected WPP, with a set point voltage ranging from 95% to 105% of the grid rated voltage [4]. Furthermore, the step-V PCC variation in response to the changes in wind power penetration should generally be maintained at less than 3% [5]. Nevertheless, two main reasons have incurred challenges in meeting the requirements for grid code: the provisions of WPP design regarding the placement of the Point of Common Coupling (PCC), and the limitations in WPP capability to control the terminal voltage [6].
In a distribution grid connected WPP, the overall system impedance angle (X/R) ratio seen at the PCC (X/R PCC ) and the Short Circuit Capacity (SCC) are important parameters that affect V PCC stability [7,8]. Short Circuit Ratio (SCR), defined as the grid's SCC divided by the power generated by WPP (P wind ) [9], is another determining factor. A SCR greater than 20 signifies that the stability requirements of the grid codes have been satisfied [10]. At a given point of a distribution system, SCC is proportional to the square of the system nominal voltage and inverse of the magnitude of system short circuit impedance (Z sc ) seen at this point [11]. Generally, the wind velocity is high at sites located far from the distribution substation. Therefore, WPPs are usually connected to the distribution systems through long lines, which makes the Z sc value measured at the PCC large and the values of SCC and SCR small. The SCR value in many grid-connected WPPs is less than 10. The small values of the SCR leads to substantial difficulties to keep the V PCC profile and step-V PCC variations in the limits of steady-state standard defined by the grid codes [12]. Hence, it is critical to determine the optimal size of WPP and select the best site for the interconnection of the WPP to the grid in a way that the grid code requirements, in regards to steady-state voltage stability, are satisfied.
The optimal siting and sizing of Distributed Generators (DGs) in distribution networks has been addressed in literature by the use of analytical methods through improving the voltage stability and minimising the power loss [13][14][15]. The analytical approaches mainly rely on calculating the bus impedance matrix, the inverse of the bus admittance matrix, and the Jacobean matrix. The large size and complicated structure of distribution systems, however, detrimentally affect the efficiency of the algorithms developed for calculating these matrices. Moreover, the assumptions used for the simplifications of calculating Jacobean, impedance and the inverse of bus admittance matrices are not valid in distribution systems [16]. In particular, the inverse of admittance matrix is not applicable to distribution systems consisting of overhead lines as the shunt admittance of the overhead lines is ignorable resulting in singular admittance matrix [17].
A. Onlam et al. proposed an optimization algorithm for power loss reduction in a distribution system through network reconfiguration and DG installation [18]. Their algorithm is based on the simulation model of a 16-bus distribution network. However, the efficiency of the algorithm proposed in [18] is adversely impacted when it is applied into complex networks with a large number of buses.
Artificial intelligence has been widely used in the optimal sizing and interconnection siting of DG units. The efficiency of a method based on the Genetic Algorithm (GA) was compared with that of an analytical method by Pisica et al. [19] in the optimal siting and sizing of DG units in an IEEE 69-bus distribution network. The analytical method was based on nonlinear optimization. The results show that GA can provide a higher accuracy compared with the non-linear optimisation method when the number of DGs connected to the test distribution system is large. On the other hand, the results achieved by both approaches had similar accuracy when the number of DGs is small. Therefore, the results demonstrated that the analytical technique cannot effectively deal with the system complexity as the number of DGs increases, while GA does not require any computational derivatives. However, the main drawback of the GA-based method proposed in [19] is that it requires modelling the entire power system. Other Artificial intelligence-Based approaches for optimal sizing and siting of DG units presented in literature includes Particle Swarm Optimisation (PSO) [20], simulated annealing (SA) [21], the hybrid genetic algorithm (HGA) [21], an Artificial Bee Colony (ABC) [22] and the decision-making algorithm [23]. These approaches require modelling and simulating complicated test distribution systems.
The authors of [24] presented a strategy for the optimal siting of wind farms considering the local wind potential of a given interconnection site, wind speed correlations to other interconnection points and their installed capacities. The main advantage of their work is to simplify WPP site allocation by removing the need to model the system. However, their studies lack the investigation of the impact of WPP interconnection on the voltage stability. Golieva proposed an appropriate approach for the siting of wind farms into distribution networks by developing a numerical relation between V PCC and the system SCR. The proposed relation enables to predict the voltage value at a given distribution feeder and find the optimal size and interconnection site of WPP while ensuring that the voltage regulation requirements are provided without carrying out complex calculations and the computational modelling of test systems. However, the author did not take into account all the steps required for validating the developed relations. For instance, the analysis was performed on a test system having SCR in the range of 0 and 2.5. According to the guidelines provided by Australia Energy Market Operator (AEMO) [25], the wind turbine generator may trip when SCR < 2, indicating that the SCR value should not be less than 2 in actual systems. This demonstrates that the numerical relations developed in [26] cannot be applied in practice. In addition, the main disadvantage of the aforementioned work is that it does not consider the relation between the V PCC and X/R PCC ratio. As discussed earlier, in distribution networks, the V PCC stability has a high dependency on X/R PCC . Therefore, the validation of the equation proposed in [26] is adversely affected due to the lack of consideration of the relation between V PCC and X/R PCC .
Referring to the previous works reviewed in this section, the majority of existing techniques applied for finding optimal size allocation of DG systems require complicated computational tasks and/or designing the simulation model of the entire distribution systems. Furthermore, the limited number of previous works conducted to simplify the WPP siting and sizing did not consider the effect of all critical parameters on the voltage stability at the PCC point. These issues adversely impact the usability of the existing approaches for the sizing and siting of WPPs in distribution networks. Sifting the literature, it was observed that there is a knowledge gaps regarding: • A simple WPP sizing and allocation approach: As discussed, the existing DG sizing and allocation methods suffer from complex calculations, modelling, and simulation issues. Therefore, a simple approach that does not require the analytical or simulation model of the distribution system is still a noticeable gap in the literature • A holistic mathematical relation for WPP sizing and siting: The problems of existing DG placement and sizing methods in regard to calculating large and complex dimensional matrices and modelling the entire test system can be removed using a holistic mathematical model between steady-state V PCC and the key parameters of distribution systems. However, only one reference [26], addressed this issue by developing a mathematical relation between V PCC and SCR. The work did not consider the necessary factors needed for validating the results, such as consideration of X/R as one of the parameters of the proposed equation or the realistic range of SCR ratio.
The novelty of the work is to develop a new mathematical approach in response to the existing challenges of analysing the optimal siting and sizing of WPP as mentioned above. This work stems from our previous research [27,28]. The voltage stability analysis carried out in [27] demonstrated that there is a relation between V PCC , SCC and X/R PCC and P wind . However, in [27], the relation between the aforementioned parameters were presented using power-voltage and reactive-power-voltage curve characteristics and the work has not proposed a mathematical relation between the parameters. The analysis studies carried out in [27] was further developed in our work in [28] by proposing an initial voltage stability mathematical model to show the relation between V PCC and key distribution network connected to Induction Generator (IG)-based WPP. However, the model presented in [28] lacks the investigation of WPPs that are based on the Double Fed Induction Generator (DFIG), which is one of the most common type of WTGs. In addition, the mathematical model developed in [28] has not considered all the steps required for validating the developed relations such as verifying the model using test systems that are different from those used for developing the model.
The hypothesis in this research was to create a holistic mathematical voltage stability model to simplify the sizing and siting of IG-and DFIG-based WPP in distribution network. The proposed model relies on the mathematical relation between PCC bus voltage and the PCC characteristics of a distribution network penetrated by IG-and DFIG-based WPP. The objectives of this research were to: • Propose a novel voltage stability mathematical model demonstrating the mathematical relations between V PCC , P wind , SCC, and X/R PCC .

•
Validate the accuracy of the proposed analytical model in predicting three important voltage stability criteria at a given connection point of a distribution network penetrated by wind power, including: V PCC profile, step variation of V PCC due to the change of P wind (∆V PCC ), and the WPP maximum permissible size ensuring the grid code requirements (P max ), using standard test systems.
To achieve the aforementioned aims, a sensitivity analysis was first carried out to gather datasets, which were used to identify a numerical relation between the voltage and X/R ratio at the PCC of various IEEE test systems with different SCR values. Analysis was carried out based on both IG and DFIG. For each generator type, the obtained V-X/R data points were then used to develop general forms of equations capable of modelling the relations between the voltage, P wind , and the PCC parameters. A Genetic Algorithm (GA)-based approach was then used to identify the coefficients of the developed equations. The accuracy of the proposed equations was then evaluated using different scenarios involving a wide range of operating conditions. The methodology proposed is based on all critical parameters that impact the voltage stability at potential distribution system WPP interconnection point. Therefore, the model enables to predict the important voltage stability criteria at a given PCC site. Using the predicted criteria, the design engineers can select the optimal site and size of WPPs in distribution grids, ensuring that the grid code requirements in regard to the voltage stability are provided.
Achieving the aforementioned aims made this work a leader in the research deals with modelling the mathematical relations between the key parameters of distribution network connected WPPs and main voltage stability criteria at the PCC. The mathematical relations enable prompt computation of the size of a WPP. It also enables the identification of the appropriate PCC buses in which it is assured that V PCC profile and ∆V PCC are in the standard range defined by the grid codes given the SCC and X/R PCC values. In the proposed mathematical model, the key PCC parameters, i.e., SCC, X/R PCC , are the only unknown of the developed equations. These parameters are usually available as the baseline characteristics of the distribution grid connected WPPs or can be easily calculated at any point looking back to the distribution substation. Hence, the proposed equations enable WPP siting and sizing without the need to solve complex and time-consuming computational tasks, which is the main advantage of the presented method over the existing approaches.
The paper is organized as follows: Section 2 presents the methodology and the specific steps carried out in this study to propose the mathematical model, including: the modelling test systems used for providing data needed for developing a series of alternative equations, identifying the coefficient of the developed equations using a GA-based approach, determining the most accurate equations among the developed alternative equations and proposing the final form of the mathematical model. Section 3 demonstrate the significance of the methodology proposed by comparison with the main previous works available in literature. The accuracy of the proposed model in predicting the critical voltage stability criteria was investigated in Section 4, and, then, Section 5 summarizes the major conclusions of this study.

Developing Mathematical Formulations
As discussed in Section 1, this paper aims to develop the mathematical relationships between V PCC , WPP size and inherent PCC characteristics of a distribution system, including: the SCC and X/R PCC ratio through voltage stability analysis studies. This was primarily conducted by investigating how V PCC behaves in response to changes in X/R PCC for networks with different SCC and SCR ratios. Analysis was preformed based on two types of generators most commonly used in wind turbines: the Induction Generator (IG) and the Double Fed Induction Generator (DFIG). Different steps followed for developing the proposed mathematical model have been presented in the sub-sequent sections.

Methodology
The overall methodology used in developing the mathematical formulations is presented below. Modelling the test systems and measuring the SCR ratio for each system: Four distribution system models defined by the Institute of Electrical and Electronics Engineers (IEEE) were simulated and considered to provide data required for developing the mathematical relations. After modelling the test systems using MATLAB/Simulink (version 2019 developed by MathWorks, Natick, MA, USA) the short circuit current value was obtained at the PCC of each system. The I sc values were obtained and then the grid's SCC was calculated using (1). Finally, the SCR value was given by (2) regarding the amount of power injected by the wind farm in each test system.
where V rated is the nominal voltage of the distribution feeder, I sc is the short circuit current at the PCC and P wind is the amount of power generated by WPP. Plotting V PCC -X/R PCC characteristics: For each test feeder modeled in the previous step, the X/R PCC ratio was measured and varied under a fixed SCC and SCR network. The X/R PCC ratio was changed throughout the course of the analysis by changing the X/R ratios of the distribution lines. Then, V PCC value was recorded taking note of the variations in response to the change of the overall X/R PCC ratio. The V PCC versus X/R PCC characteristic has been plotted for each test feeder using the obtained data.
Developing alternative mathematical functions modelling the relations between V PCC , X/R PCC , SCC and P wind : The voltage behaviour, with respect to changes in the X/R PCC ratio, was then analysed for all test systems to mathematically correlate V PCC , X/R PCC , P wind and SCC. At the end of this step, the general forms of alternative equations, which can be used to describe the relation between the aforementioned parameters were developed. Then, the coefficients of the equations were identified using an Artificial Intelligence (AI)-based approach.
Determining the best fit function: The accuracy of each alternative equation was investigated in order to identify the best fit equations.
Developing the final form of the functions: Finally, it was investigated how the functions presented in the previous step can be further developed for different operating conditions. Details about the aforementioned steps are provided in Sections 2.2-2.6.

Modelling Test Systems
In this paper, four test distribution networks were modelled and simulated to develop the mathematical formulations. The adoption of four systems was to take into account various values of SCC and X/R PCC and hence develop a formulation with higher validity. The modelling of two test systems, referred to as Test 1 and Test 2, relied on the 37-bus IEEE distribution network (Figure 1), while the modelling of other two test systems, referred to as Test 3 and Test 4, was based on the 9-bus IEEE distribution network ( Figure 2). The PCC was considered to be at Bus 9 in the 9-bus system and at Bus 6 in the 37-bus system. The size of the WPP and the length of the lines are different amongst the four test systems, resulting in different values of SCC and SCR.
In each model, the distribution system is connected to a 120 kV, 50 Hz external source of infinite short circuit current through a Yg/∆ configuration of a three-phase distribution transformer. Furthermore, in each model, a WPP with the total rated capacity of 9 MVA, provided by three 3 MVA WTGs, is connected to the distribution system at the PCC. A three-phase transformer with Yg/Yn configuration is connected to each WTG to boost the generator terminal voltage to the rated voltage of distribution system, i.e., 22 kV. The parameter values of different components of the simulated models are presented in Tables A1-A4 in Appendix A.  The feeder steady-state voltage before the WPP connections, called V initial , has a significant impact on the V PCC stability after the WPP connection. Developing the mathematical relations for a constant V initial value hinders the application of the relations for any distribution network, as loading conditions and various voltage regulator set-point values impact the V initial value. In this study, it was firstly assumed that V initial is around 0.98 pu for each test distribution feeder. Therefore, the mathematical relations were firstly developed regarding V initial = 0.98 pu. Then, in Section 2.6, the proposed relations were further developed such that it satisfies a wide range of V initial values. The adoption of various V initial values in the proposed methodology was to reduce the uncertainty due to the load deviations from the values used in the test systems. Table 1 presents the values of PCC parameters, including: I sc , SCC, P wind , and SCR, for each test system. Among the cases presented in Table 1, Test 1 with the least SCR value, is the weakest system. Moreover, the lowest value of SCC and P wind is seen for Test 4 with the greatest SCR value, which is due to the small P wind . This results in the PCC point in Test 4 being the stiffest among other test systems.

Characteristics of Voltage-X/R Ratio at the PCC
In this paper, the proposed methodology is based on the V PCC -X/R PCC data obtained from simulation models. The reason for using simulation models for obtaining V PCC -X/R PCC data is the fact that the development of mathematical relations with an appropriate accuracy requires a large number of data points. The lack of such a data base adversely impacts the accuracy of the mathematical relations. The use of simulation models enables to provide a large number of data by carrying out sensitive analysis through changing the X/R PCC ratios and monitoring V PCC taking note of the variations against the change of the X/R PCC ratio. The X/R PCC is an inherent characteristic of a power system and cannot be dynamically changed in an actual network. Therefore, having access to a large quantity of V PCC -X/R PCC data points obtained from actual systems is a complicated process.
From a practical perspective, the application of the proposed analytical model to the real-world distribution systems may impose additional complexity and challenges due to the uncertainty of the simulation models. To enhance the accuracy of the proposed model against the issues concerned with the uncertainty of the simulated test systems, the range of the X/R PCC ratio, analysed in this study, was derived from an analysis of actual distribution networks performed in [29]. Reginato et. all [29] carried out a research study concerned with an analysis of the effect of X/R PCC and the inverse of SCR, the so-called integration level (p), on V PCC using data from actual systems. The analysis was performed to determine a potentially viable X/R PCC range, ensuring the grid code requirements with regard to the steady-state V PCC stability for a given integration level. The analysis results were presented using p vs X/R PCC graphs, while 0.95 pu < VPCC < 1.05 pu. Although the authors in [29] did not present the value of V PCC for a given X/R PCC and p, the realistic range of X/R PCC has been identified for different integration levels using data obtained by actual systems. In the V PCC -X/R PCC characteristics presented in this section, the range of X/R PCC and the corresponding SCR value are based on the p-X/R graphs presented in [29] to decrease the issues related to the uncertainty of the data obtained from simulation models, i.e., V PCC -X/R PCC data points, in case of real application.
To draw the V PCC /X/R PCC curve characteristics, the X/R PCC ratio of the test systems in Section 2.2 was changed to observe the V PCC values for each X/R PCC with SCC and SCR values of the test system being constant. The data points collected were utilised to plot the characteristic of V PCC against the X/R PCC ratio. The results are presented in Figure 3. For the IG-based WPP, Figure 3, in X/R PCC around 2, the V PCC values for all four systems are nearly equal. This observation is in accordance with [29], which states that the lowest variation in voltage happens at X/R PCC = 2 irrespective of the SCR value in an IG-based WPP, which is another piece of evidence to confirm that the methodology results are based on realistic cases.
The PCC voltage for both IG-and DFIG-based WPPs is high at the points with a low ratio of X/R PCC (X/R PCC < 2). This is particularly apparent in small SCR ratios. For example, from Figure 3, it can be seen that the PCC voltage of Test 1 with an SCR ratio equal to 3 is greater than the maximum acceptable voltage level defined by the grid codes (V PCC > 1.05 pu) when the X/R PCC ratio is smaller than 0.5. This is in contrast with Test 4 with the maximum SCR ratio equal to 7, where the low values of the X/R PCC ratio do not cause a noticeable increase in the voltage. When dealing with large X/R PCC ratios, Figure 3 demonstrates that the voltage variation due to changes in wind power penetration is negligible at the connection point of a DFIG-based WPP with X/R PCC > 2. Furthermore, Figure 3 shows that voltage at the PCC of IG-based WPP is small when the X/R PCC ratio is large (X/P PCC > 2). The variations in voltage according to a rise in X/R PCC ratio is more remarkable at the weaker test feeders. In Figure 3, Test 1 shows a V PCC value less than the minimum limit of the allowable range, i.e., V PCC < 0.95 pu, when X/R PCC is around 2.5. The voltage for Test 4 is shown to be within the standard range, even in X/R PCC ratios greater than 5.

Developing Alternative Functions
In this section, the main aim is to present the mathematical relations which can model the V PCC -X/R PCC characteristics shown in the previous section. The relations were developed for the values of X/R PCC that imposes Power Quality (PQ) concerns and voltage regulation problems at the PCC. Therefore, both small and large X/R PCC ratio ranges were considered for developing the formulation for the IG-based WPP. Furthermore, only small X/R PCC ratios were considered for expressing the mathematical relations for the DFIG-based WPP.

General Form of Alternative Functions for IG-Based WPPs
Mathematical formulas were developed to introduce the best fit curve for the V PCC -X/R PCC data points presented in Figure 3. Two key features in Figure 3 were considered in developing the equations, which are presented in Table 2.

Feature Description
Feature 1 Voltage at feeders with a low X/R PCC ratio is increased as SCR is decreased. Feature 2 Voltage at feeders with a high X/R PCC ratio is decreased as SCR is decreased.
Two possible functions, including an exponential function and a polynomial function with an order of 2, were proposed and investigated to cover the characteristics listed in Table 2. The general forms of the functions are as follows: Function 2-Polynomial with an order of 2 where Z 0 , Z 1 , Y 0 , Y 1 and Y 2 , are positive coefficients. γ and δ represent the exponential and polynomial decay coefficients, respectively.
Mathematically, it can be tested and confirmed that the V PCC -X/R PCC characteristic curves given by (3) and (4) have similar patterns for all X/R PCC ranges and different decay coefficients. Hence, both functions can be applied to model the relation between V PCC and X/R PCC . Figure 4 shows the general V PCC -X/R PCC curve of the considered alternative functions for different values of decay coefficients.   (4) shows that smaller decay coefficients (γ and δ) increase the PCC voltage. Looking at the first characteristic presented in Table 2, small SCR values increase voltage at the PCC points with a small X/R PCC ratio. Therefore, it can be concluded that SCR and the decay coefficients are directly related when X/R PCC is small (X/R PCC < 2), expressed by (5) and (6).
where K 1 and K 2 are considered positive. Figure 4 demonstrates that higher values of decay coefficients results in a decrease in V PCC . Referring to the second characteristic explained in Table 1, small SCR values decrease voltage at the feeder with a large X/R PCC ratio. Hence, SCR and decay coefficients have an inverse relationship when X/R PCC is large (X/R PCC > 2). This relationship is given by (7) and (8): δ= where K 3 and K 4 are considered positive. Noting the relation between decay coefficients and SCR given in (5) to (8), the alternative functions (3) and (4) can be rewritten as (9) to (12).
where Z s0 , Z s1 , Y s0 , Y s1 and Y s2 are positive coefficients in the equations representing the small X/R PCC range (X/R PCC < 2) and Z l0 , Z l1 , Y l0 , Y l1 and Y l2 are positive coefficients in the equations representing the large X/R PCC range (X/R PCC > 2).

General Form of Alternative Functions for DFIG-Based WPPs
Comparing the results obtained for IG-and DFIG-based WPPs shown in Figure 3, the V PCC variation against changes in X/R PCC follows a similar pattern in both IG-and DFIG-based WPPs when X/R PCC < 2. Hence, the V PCC -X/R PCC curves plotted for DFIG and IG are similar to each other in the small X/RPCC ratio, where DFIG cannot efficiently regulate the voltage. This signifies that the alternative functions applied for the DFIG-based WPP are similar to those developed for the IG-based WPP with small X/R PCC . The general forms of the functions are expressed in (13) and (14).
where W s0 , W s1 , W s0 , U s0 , U s1 , U s2 and K 5 are positive coefficients. Due to the crucial effect of the coefficients on the accuracy of Equations (9)- (14), the next section is devoted to determining these coefficients by the use of a GA-based approach.

Determining the Coefficients Values Using GA
The Genetic Algorithm uses an exploratory approach to optimise a given objective function utilising the ideas of evolutionary processes [30]. Although GA is available in the MATLAB Optimization Toolbox, coding is required for defining the objective function and GA parameters.
The value/type of the most important GA parameters considered in developing the mathematical model presented in this paper is tabulated in Table 3. The default number of generations and population size in the MATLAB GA toolbox are 100 and 20, respectively. As shown in Table 3, the generation and population sizes considered in this paper are noticeably larger than the default values increasing the accuracy of the GA outputs. From Table 3, the selection of individuals in current generation as the parents for the next generation is based on stochastic uniform, which provides the highest accuracy compared with other search strategies defined in GA [31]. Furthermore, the type of crossover parameter, which classifies two parents of the intermediate generation to build a new individual for the next generation, is 'scattered' as it provides higher flexibility compared with the other types of the crossover functions [32]. It can also be seen from Table 3 that the type of mutation function is constrained dependent. Another important GA parameter is mutation, which makes small randomly changes in the genes of parents to create mutation individuals for the next generation [33]. Mutation function is categorised into two types: unconstrained and constrained dependent. As shown in Table 3, in this study, a constrained dependent mutation function has been used and the constraints are related to the upper and lower boundaries of the input parameters of the objective functions.
A separate objective function has been coded for each Equation in (9)- (14), which resulted in 6 objective functions, as shown in Table 4. Table 4. Characteristics of the objective functions.

Objective
Function Name

Input of Objective
Function Output of Objective Function Coefficients of (9), i.e., Z s0 , Z s1 and K 1 Error between the V PCC -X/R PCC graphs predicted by (9) and the reference curves presented in Figure 3 for IG-based WPP and X/R PCC < 2 (10) Func2 Coefficients of (10), i.e., Y s0 , Y s1 and Y s2 Error between the V PCC -X/R PCC graphs predicted by (10) and the reference curves presented in Figure 3 for IG-based WPP and X/R PCC < 2 (11) Func3 Coefficients of (11), i.e., Z l0 , Z l1 and K 3 Error between the V PCC -X/R PCC graphs predicted by (11) and the reference curves presented in Figure 3 for IG-based WPP and X/R PCC > 2 Coefficients of (12), i.e., Y l0 , Y l1 and Y l2 Error between the V PCC -X/R PCC graphs predicted by (12) and the reference curves presented in Figure 3 for IG-based WPP and X/R PCC > 2 Coefficients of (13), i.e., W s0 , W s1 and K 5 Error between the V PCC -X/R PCC graphs predicted by (13) and the reference curves presented in Figure 3 for DFIG-based WPP and X/R PCC < 2 Coefficients of (14), i.e., U s0 , U s1 and U s2 Error between the V PCC -X/R PCC graphs predicted by (14) and the reference curves presented in Figure 3 for DFIG-based WPP and X/R PCC < 2 From Table 4, the values of the input parameters of Objective Functions 1 to 6 (Func1 to Func6) represent the values of the coefficients of (9) to (14), respectively. In each objective function, the coefficient of the corresponding equation was considered as the input parameters, while the error between predicted graph and the reference V PCC -X/R PCC curve was considered as the output parameter. The main idea was to find the coefficient of each alternative equation in a way that the curve predicted and plotted by the equation has the lowest error with respect to the corresponding reference V PCC -X/R PCC curves shown in Figure 3. The error was calculated using Standard Deviation (SD) formulas.
For clarifying the operation of the objective functions, Figure 5 presents the flowchart diagram for one of the objective functions used for finding the coefficients of (9), named Func1. Once the objective function has been coded and GA parameters have been defined, the optimal values of the function input parameters was determined by running a separate GA for each objective function. GA was responsible to find values of the input variables for which the output of the objective function becomes minimum. The values are as shown in Tables 5-7.  (9)-(14) were rewritten to give Equations (15)- (20).

Finding the Fit Functions
The best fit function with the least amounts of errors compared with the reference V PCC -X/R PCC data points was obtained for Equations (15)- (20). As stated in Section 2.4.3, the output parameter for GA objective functions was considered to be the error between reference and predicted values and was calculated by SD formulas. In addition to SD, the Mean of Relative Error (MRE) and the Mean of Absolute Error (MAE), as two frequently used criteria, were considered for evaluation of the error in this section.
For each alternative function and test system, the error between the reference and predicted V PCC -X/R PCC curve characteristics were calculated using SD, MRE and MAE formulas, as shown in (21) to (23).
where the reference values of V PCC is expressed as V k , and its predicted values are expressed asV k . Moreover, the total number of reference values is given as m.
The values of the three factors used for evaluating the error between the reference and predicted V PCC -X/RPCC values have been presented in Table 8 for each alternative function and test system. Furthermore, the difference between the reference V PCC -X/R PCC curve characteristics obtained by the simulation results and the V PCC -X/R PCC curve characteristics plotted by the developed alternative equations can be compared for each test system using graphical representation, as shown in Figure 6a-l.    (15) and (16) for Test 1; (b) V PCC -X/R PCC curves obtained by (15) and (16) for Test 2; (c) V PCC -X/R PCC curves obtained by (15) and (16) for Test 3; (d) V PCC -X/R PCC curves obtained by (15) and (16) for Test 4; (e) V PCC -X/R PCC curves obtained by (17) and (18) for Test 1; (f) V PCC -X/R PCC curves obtained by (17) and (18) for Test 2; (g) V PCC -X/R PCC curves obtained by (17) and (18) for Test 3; (h) V PCC -X/R PCC curves obtained by (17) and (18) for Test 4; (i) V PCC -X/R PCC curves obtained by (19) and (20) for Test 1; (j) V PCC -X/R PCC curves obtained by (19) and (20) for Test 2; (k) V PCC -X/R PCC curves obtained by (19) and (20) for Test 3; (l) V PCC -X/R PCC curves obtained by (19) and (20) for Test 4.
Based on Table 8, in IG-based WPP with X/R PCC < 2 for all test cases, the evaluation criteria calculated for (15) is less than the criteria obtained for (16). This is confirmed in Figure 6a-d, where the curve for (15) shows a smaller error than the curve for (16) compared with the reference curve. Analysis of IG-based WPP with X/R PCC > 2 for all test cases based on Figure 6e-h and Table 8 reveals a slightly higher accuracy for (18) compared to (17).
In the case of DFIG-based WPP, Figure 6i-l and Table 8 show that the exponential function provides less error compared to the polynomial function. Therefore, according to the statistical and graphical results presented, the best fit functions are as listed in Table 9. Table 9. Best fit functions considering WPP type and X/R PCC range.

WPP Type X/RPCC Range Proposed Equation
The mathematical relations developed in Section 2.2 is based on an initial voltage of 0.98 pu for the test systems, hence Equations (15), (18) and (19) would be only used for the feeders with V initial of 0.98 pu. Nevertheless, the initial voltage varies from the default value because of the loading conditions and values of voltage regulator set-point. Consequently, another set of equations is required to calculate the voltage stability criteria at a certain test system with any initial voltage value.

Developing the Final Form of the Functions
To establish the final form of the developed mathematical equations, the best fit functions presented in Table 9 were further developed for various initial voltage levels. The impact of variations in the initial voltage on the P wind versus V PCC curve was evaluated using simulation systems characterized in Table 1. The initial voltage was varied in the simulations from the default value (V inital = 0.98 pu) to plot the variations of P wind against V PCC for each new voltage value. The shape of the curves remained the same with both the default and varied initial voltages. The only difference was a downward shift for initial voltages less than 0.98 pu and an upward shift for the initial voltages greater than 0.98 pu. Defining ∆V initial = V initial -0.98, it can be approximated that for a given initial voltage, the addition of ∆V initial to the PCC voltage value calculated by the equations in Table 9 gives the voltage value at PCC point. Hence, V initial parameter was incorporated in (15), (18) and (19) to give (24)-(26) in Table 10, respectively.
As the SCR is equal to the ratio of the grid's SCC to P wind , Equations (21)-(23) can be rewritten to calculate optimal value of wind power generation (P max ) such that the voltage value at PCC point is kept in the standard range. To do so, in Equation (25) related to WPP with large X/R PCC, the voltage value at the PCC point was assumed to be 0.95 pu, which is the minimum allowable voltage level defined by the Australian grid code. In Equations (24) and (26), developed for WPP with a small X/R PCC ratio, the voltage at PCC point was assumed to be equal to the maximum limit of the admissible steady-state voltage.
After determining the V PCC values, the equations were rewritten in terms of P wind as a function of PCC parameters and V initial as can be found in (27) to (29) in Table 10.

Equation Equation Number Application
Plotting PV curve for IG WPP with X/R PCC < 2.

The Significance of Proposed Model
In this section, the proposed mathematical model is compared with the main previous works available in literature to demonstrate the significance of the research performed in this paper and its contribution to the existing knowledge.
The majority of the previous works dealt with the effect of SCC, SCR and X/R ratios on V PCC in general terms and using graphical representation and different curve characteristics such as V PCC versus X/R PCC , V PCC versus SCR etc. For example, Kothari and et al. investigated the general relation between V PCC , X/R PCC and SCR in distribution systems connected WPP [34]. Their work presented the analysis results using the curve characteristic of V PCC versus SCR plotted for different X/R PCC ratios. For the X/R PCC = 0.5 case, the results in [34] demonstrated that the difference between the initial value of V PCC (V initial ) when P wind = 0 and the V PCC value obtained for high wind power penetration is around 7%. Referring to (2), the high penetration of wind power makes the SCR value small. Hence, in [34], it was demonstrated that the difference between V initial and V PCC after WPP connection is around 7% (∆V PCC = 7%) when X/R PCC is 0.5 and SCR is noticeably small.
The results presented in [34] were compared with the results obtained by the equations proposed in this paper. For this purpose, the difference between the V initial and V PCC values after WPP connection was calculated using proposed equations considering that X/R PCC = 0.5 and SCR is small. The calculations were performed assuming that SCR is slightly greater than 2, i.e., SCR = 2.5. This is because, as mentioned in Section 1, AEMO documentation indicates that an SCR < 2 should be avoided in WPPs, as it adversely impacts the voltage stability in the steady-state operation and may lead to generator tripping [25]. The comparison results have been shown in Table 11.
As shown in Table 11, the results gained by the proposed equations confirm the discussion presented in [34] that the difference between V initial and V PCC under high wind power penetration is around 7% when X/R PCC = 0.5. Although the authors in [34] investigated the general relation between V PCC , SCR and X/R PCC using graphical representation, they did not propose any mathematical relation between these parameters. This shortcoming hinders to investigate the voltage variation for any given X/R PCC and SCR values. Table 11. ∆V PCC obtained by (24) and (26) for X/R PCC = 0.5 and SCR = 2.5. In [12,35], it was demonstrated that the steady-state V PCC tends to increase for small X/R PCC ratios, while it decreases for large X/R PCC ratios in weak distribution networks penetrated by DGs. Furthermore, the authors in [36] showed that the voltage drop, due to an increase in X/R PCC , is more serious in distribution systems with low SCR values (weak systems) compared with that in the distribution feeders with larger SCR values (stiff systems). The results presented in [12,35,36] confirm the general characteristics of the functions given by the equations presented in Table 10. However, the authors in these works, only discussed the relation between V PCC , X/R PCC and SCR in a general manner and using graphical representation and their works do not propose a mathematical relation between these parameters.

WPP
The author of [26] offer one of the most valuable preliminary works, discussing the mathematical formulation of any formula to express the relationships between V PCC and SCR for the steady-state operation of the WPPs. It was demonstrated that the relation between V PCC and SCR can be expressed through a polynomial function with an order two as shown in (30) [26]. V PCC = V initial × 0.0067 × SCR 2 − 0.063 × SCR + 1.1142 (30) As claimed by the author in [26], taking advantage of this relation, the voltage profile can easily be predicted using the SCR value at each distribution feeder. This is a great achievement as it would allow WPP planning engineers to select the best site for the interconnection of a WPP without the need to carry out complex and time-consuming computational tasks and modelling test systems. However, the accuracy of the equation proposed in [26] is adversely impacted due to the lack of considering the relation between X/R PCC , which is one of the most important characteristics of distribution networks, and V PCC . In addition, the equation was developed using data provided by an invented test system with unrealistic SCR range (0 < SCR < 2.5), which adversely impact the accuracy of (30) in real distribution systems.
Sifting the literature, it is clear that a holistic relation between V PCC and the key parameters of the distribution systems, i.e., the SCR and X/R PCC ratio, is still a noticeable gap. This shortcoming has been eliminated in this paper using the proposed mathematical model presented in Table 10. Furthermore, the proposed equations have been developed using data obtained by IEEE standard distribution networks.
Referring to Table 10, the proposed model enables design engineers to plot the curve characteristic of P wind versus voltage at a potential distribution network WPP interconnection point using (24) and (26), considering the WTG type and the values of SCC, V initial and X/R at that point. Taking the advantage of the P wind versus voltage curve characteristic achieved by the proposed equations, the engineers can predict the V PCC value for different wind power penetration at a potential interconnection point. The predicted V PCC values enable the grid codes' compliance check, i.e., to verify if V PCC is between 0.95 and 1.05 pu (compliance with the grid codes) or if V PCC is not maintained within 0.95 and 1.05 pu (grid codes violation), at the interconnection of WPP to a potential interconnection site. The potential site can then be selected for connecting WPP if the grid code requirements are met or identified as an inappropriate PCC point if the grid code requirements are not provided.
As another advantage of the proposed mathematical model, the engineers can estimate the step-V PCC variations in response to changes in wind power penetration at a given PCC point using (24) to (26) depending on the WTG type and the X/R PCC ratio. The predicted results obtained for the step-V PCC variations enables the engineers to decide how many WTGs can be simultaneously switched on, ensuring that the step-V PCC variation satisfies the grid code requirements (∆V PCC ≤ 3%).
Furthermore, the proposed model enables the network and grid interconnection engineers to find optimal WPP size and estimate the maximum power injection by a WPP to ensure voltage stays within acceptable steady-state limits at a given PCC point using (27) to (29). Therefore, the proposed mathematical model provides design engineers with a series of informative equations to find the optimal size and site of WPP in a distribution network by predicting three critical voltage stability criteria, i.e., V PCC profile, step-V PCC variation and maximum permissible wind power generation, ensuring the grid code requirements.
As discussed in Section 1, the majority of existing approaches used for WPP siting and sizing are based on the modelling and simulation of the whole system and/or calculating the bus impedance, inverse of the bus admittance, and the Jacobean matrices. However, the simulation of distribution networks is a demanding process due to the size and complexity of these networks. Moreover, the assumptions used for the simplifications of calculating Jacobean, impedance and inverse of admittance matrices are not valid in distribution systems. These issues are removed using the proposed mathematical model presented in Table 10. The proposed equations require only three input parameters (V initial , X/R PCC and SCC) which are usually available as the baseline characteristics of the distribution systems or can easily be calculated at any point looking back to the distribution substation. Hence, the proposed equations enable WPP siting and sizing by promptly conducting three important voltage stability criteria at a potential distribution network interconnection point. Consequently, it removes the need to solve complex and time-consuming computational tasks or modelling the entire test distribution system, which is the main advantage of the presented over the existing WPP sizing and siting approaches. The proposed analytical model removes the need to investigate the effect of distribution network configuration and its component specifications on PCC bus voltage stability. This is due to the fact that the effect of these factors has been considered and modelled in the proposed analytical approach using SCC and X/R PCC parameters.
The validation of the proposed mathematical model in predicting the aforementioned three voltage stability criteria at a potential distribution network WPP interconnection point has been conducted in the subsequent sections of this paper.

Validation Results
This section presents the analyses performed to validate the proposed functions, listed in Table 10, which can be used to project three important V PCC stability criteria, including: V PCC , ∆V PCC , and P max for a potential connection site in a distribution network, given the V initial value. In this respect, the accuracy of the proposed functions was investigated using two case studies. Case study 1 provides initial validation studies by assessing the accuracy of the proposed mathematical model using the four test systems presented in Section 2 with different X/R ratios and V initial values. Then, in Case study 2, the functions developed were applied for a new test system based on a 30-bus IEEE distribution network.

Case Study 1-Validation for Different Vinitial Values
As discussed in Section 2.6, loading conditions and various voltage regulator set-point values may result in different V initial values. Therefore, the effect of load variations and any other factors that change operating voltage at a given distribution feeder before the WPP connection has been modelled using V initial parameter. V initial has a profound effect on the P wind versus V PCC characteristic and on the steady-state voltage stability after the connection of WPP. In this regard, the proposed mathematical model should ideally work for any V initial value to ensure that the equations developed can be applied for a wide range of distribution networks with different loading conditions. Hence, as an initial assessment of the accuracy of the proposed mathematical model, the equations presented in Table 10 are verified for different V initial values. The validation study has been performed using the four test systems applied for developing the mathematical model with details presented in Section 2.2. However, the validation studies were performed considering that the loading conditions and V initial values of the four test systems are different from the default value (V initial_default = 0.98 pu) considered in Section 2.2.
Two test systems were simulated based on a 37-bus IEEE distribution network and two test systems were simulated based on a 9-bus IEEE distribution network using MATLAB/Simulink as shown in Figures 1 and 2, respectively. For the 37-bus model, the potential distribution network WPP interconnection point was considered to be at Bus 6. Furthermore, the potential site for connecting WPP was considered to be at Bus 9 for the 9-bus model. The topology and/or length of distribution lines are different among the test systems considered making their SCC values to be different. The SCC value and topology of each test system are as listed in Table 1 in Section 2.2.
It is noted that the analysis performed in Case study 1 is the initial stage of validating the proposed mathematical model. In this stage, the SCC values and topology of the test systems used for validation studies are similar to the SCC values and topology of the test systems used for developing the mathematical model. However, the V initial values are different from that considered for developing the proposed equations. Therefore, the mathematical model is validated using new operating conditions. In addition, the proposed model is validated for various X/R PCC ratios.
In this section, six scenarios were analysed: two scenarios for the IG-based WPP with X/R PCC < 2, two scenarios for the IG-based WPP with X/R PCC > 2 and two scenarios for the DFIG-based WPP. Table 12 shows the test system, X/R PCC ratio and V initial value considered in each scenario. Furthermore, the SCC value and topology of each test system have been re-written in Table 12. For each scenario, the P wind versus V PCC curve characteristic obtained by the simulation results is compared with the P wind versus V PCC curve characteristic given by the corresponding proposed equation presented in Table 10   When predicting V PCC value, the largest error between the reference and predicted results occurs in Scenarios 1 and 5 when wind power is not connected (P wind = 0). Referring to Table 12, both the SCC and X/R PCC ratio are small in Test 1 and 5. Therefore, the accuracy of the proposed mathematical model is slightly impacted when both the SCC and X/R PCC ratio are small. However, as shown in Figure 7a,e, the highest error is 1%. In the other scenarios considered, the error between the simulation and predicted results is less than 0.5%.
In the case of the prediction of ∆V PCC , the simulation results in Figure 7a-f demonstrate that the proposed equations accurately verify the compliance or violation of ∆V PCC with respect to the grid code requirements. As defined by the grid codes, the ∆V PCC value must not exceed 3%. Comparing the curve characteristics shown in Figure 7a-f, it can be observed that both simulated and predicted results demonstrate that the highest ∆V PCC occurred in Scenario A, where the SCC and X/R PCC are both small. According to Figure 7a, both simulation and predicted results show that the V PCC variation in response to an increase in P wind from 0 to 3 MVA is greater than 3%. Therefore, the results confirm that ∆V PCC violates the grid code requirement when only one 3 MVA generator is connected to the grid. However, for the same increase in P wind , ∆V PCC does not exceed the standard range in other scenarios. In the case of scenarios with large X/R PCC (Scenarios 3 and 4), both simulation and predicted results presented in Figure 7c,d demonstrate that ∆V PCC is more serious in Scenario 3, where the SCC value of the test system considered is smaller.
Generally, comparing the simulation and predicted results in Figure 7a-f, it is clear that, in all scenarios, the proposed equations were accurate in estimating ∆V PCC , projecting the correct grid code compliance or violation outcome. Therefore, two voltage stability criteria (V PCC value and step-V PCC variation) could be predicted using the proposed mathematical model with an appropriate approximation for all scenarios considered for the 9-bus and 37-bus test systems.
Another key criterion for voltage stability is P max , which is possible to be predicted from the curve characteristics presented in Figure 7a-f. However, because there is no straight way to estimate the value of P max using PV characteristics, Equations (27) to (29) were utilized for the estimation of the P max value in Scenarios 1 to 6. For each scenario, the maximum permissible wind power generation obtained from both simulation and equations is illustrated in Figure 8. As shown in Figure 8, for the IG-based scenarios (Scenarios 1 to 4), the greatest discrepancy between the simulation and predicted results is less than 1 MW (Percent error = 13%) in Scenario 1, where the values of SCC and X/R PCC are small. The error is much less, around 0.5 MW, in scenarios with greater SCC or X/R PCC > 2. Therefore, the IG-based WPP with a small X/R PCC ratio slightly affects the accuracy of the predictions by the mathematical model. For DFIG-based scenarios (Scenarios 5 and 6), the highest error is related to Scenario 5 where the SCC value is smaller than another scenario. Therefore, small SCC values slightly impact the accuracy of the equation proposed for DFIG-based WPP. However, as presented in Scenario 5 in Figure 8, the highest error is less than 1 MW (percentage error = 13%). Therefore, in both IG-and DFIG-based Scenarios, the worst-case accuracy is 87%. It can be seen from the results that the values of P max calculated from the presented mathematical model well agree with the values obtained from simulation with only a slight error margin.

Case Study 2-Validation for a New Test System
In this section, the accuracy of the final proposed functions presented in Table 10 is validated using a new test system with different characteristics and specifications from the test systems used for developing the mathematical model to increase the validity of the developed model by considering various system structures and PCC parameters. In Case study 2, the new test system has different topology and PCC characteristics and analysis was performed for different operating conditions, i.e., different load conditions and V initial values, various SCC and X/RPCC ratios, compared to the models used for developing the proposed equations. The new system is based on the 30-bus IEEE distribution network simulated in MATLAB, as shown in Figure 9. The values of the 30-bus test system parameters, including rated voltage, frequency, the characteristics of substation transformer and specifications of WPPs, are the same as those used for the 9-bus and 37-bus test systems presented in Section 2.2. The validation studies were carried out based on three scenarios shown in Table 13. The length of lines and load demand are different among Scenarios A to B resulting in different values of X/R PCC and V initial . For each scenario, the type of WPP, PCC site location, X/R PCC ratio, grid's SCC and V initial value are as shown in Table 13. For each scenario presented in Table 13, the simulation curve characteristic of P wind versus V PCC was obtained using simulation results. In addition, the equations proposed for predicting V PCC profile and ∆V PCC , i.e., (24) to (26), were used to plot the characteristic of P wind versus V PCC . Comparing the simulation and predicted results enable one to investigate the accuracy of the developed functions in projecting the V PCC value and step-V PCC variation for different wind power penetration. The results for IG-based scenarios (Scenarios A and B) and DFIG-based scenario (Scenario C) are as shown in Figure 10a,b, respectively. From Figure 10a,b, it can be observed that the highest error between the simulation and predicted results is less than 1% in all three scenarios. Furthermore, as shown in Figure 10b, the simulation and predicted curves intersect each other at some points in Scenario C. Therefore, similar to the validation results presented in the previous section, the functions proposed for predicting the voltage value at a given connection point, i.e., (24), (25) and (26), provide significant accuracy.
As discussed in the previous section, the P wind versus V PCC characteristics enable an investigation on validating whether the proposed functions can accurately foresee ∆V PCC and project compliance or violation of the grid codes. For the scenario with the lowest SCC and small X/R PCC ratio (Scenario A with SCC = 18 and X/R PCC = 0.4), both simulation and predicted results presented in Figure 10a demonstrate that ∆V PCC violates the grid code requirements (∆V PCC > 3%) when only one 3 MVA generator is connected to the grid. Furthermore, both the simulation and predicted curves plotted for another IG-based scenario (Scenario B with SCC = 46 and X/R PCC = 3.5) demonstrate that ∆V PCC is less than 3% for different wind power penetration. In the case of DFIG-based scenario, both simulation and predicted results in Figure 10b demonstrate that ∆V PCC satisfies the standard range if the maximum increase in P wind is 6 MW, for example, when P wind increases from 0 to 3 MVA (one 3 MVA generator is connected to the grid) or from 3 to 9 MVA (two 3 MVA generators are simultaneously connected to the grid). However, ∆V PCC violates the grid code requirements when three or more 3 MVA generators are simultaneously connected to the grid. Hence, the results demonstrate that the proposed equations are accurate in predicting ∆V PCC for different wind power penetration levels.
To validate the proposed mathematical model in estimating P max , the simulation results shown in Figure 10a,b were used to determine the maximum value of wind power generation ensuring 0.95 < V PCC < 1.05. The simulation results were then compared with the values given by (27) to (29). Figure 11 presents the simulation and predicted values of P max for Scenarios A to C. As shown in Figure 11, the difference between the simulation and predicted results is negligible in all three scenarios. The highest difference occurred in Scenario A, where the predicted value is around 0.2 MW greater than the reference simulation value, which is corresponding with a 3.8% error. Therefore, the results demonstrate that the proposed equations can provide a very high accuracy in estimating the P max value for the IEEE 30-bus test system.

Conclusions
A mathematical model has been developed for the initial analysis of siting and sizing of potential distribution feeders for the interconnection of WPPs. The formulae consider a weak distribution network penetrated by WPP to determine the relation between V PCC , SCC, the overall system X/R PCC ratio and P wind . The equations were developed using the V PCC -X/R PCC data points given by four test systems with different SCC and SCR values. The test systems were modelled and simulated based on IEEE standard networks. For each test system, a voltage stability hypothesis was developed based on the SCC and X/R PCC ratio measured at a potential PCC bus. In this respect, for each test system, the V PCC -X/R PCC characteristic was obtained and plotted under a specific SCR value considering the realistic range of X/R PCC and V PCC . The analysis studies were performed based on two types of WTGs commonly used in the WPPs, including: IG and DFIG. Alternative mathematical approximations were formulated for each WTG type and the corresponding V PCC -X/R PCC characteristic. The coefficients of the developed formula were determined using the GA optimization method. The accuracy of the alternative equations was then tested using the error evaluation criteria to determine the most accurate equations for each WPP type. The equations with the lowest error were then developed to calculate V PCC for any V initial values. Finally, six equations were proposed to calculate the three critical voltage stability criteria at a given PCC point, including the V PCC profile, step-V PCC variation in response to changes in wind power penetration and the maximum permissible wind power, ensuring that the steady-state V PCC requirements defined by the grid codes would be satisfied. The validation studies conducted for the V PCC profile demonstrated that the highest error amongst the all scenarios considered was less than 1%. Furthermore, validation studies confirmed the accuracy of the proposed method in predicting the step-voltage variation grid codes compliance check, i.e., to verify ∆V PCC ≤ 3% (compliance with the grid codes) or ∆V PCC ≥ 3% (grid codes violation). For maximum permissible wind power, the verification results showed that the accuracy of the proposed relations was slightly impacted in IG-and DFIG-based WPPs with small SCC and X/R PCC ratio. However, the worst-case accuracy amongst the all scenarios investigated was around 87%. Hence, the validity of the equations developed was confirmed through comparison with simulation results, which showed a highly accurate result for the voltage stability criteria at PCC points with various SCC and X/R PCC values.
The formulae could be used to assess the voltage stability in various distribution systems connected WPP as the basis of developing equations was the key parameters that affect the voltage stability at PCC points. The developed mathematical formulations would alleviate the challenges of the siting and sizing of WPPs through eliminating complex calculations and simulations of the existing methods. The research presented can be regarded as a first step towards the WPP sizing and siting using mathematical relations. However, the work can be further extended by removing some of the scope limitations assumed in this research to make the proposed mathematical model more comprehensive and reliable.
One possible research idea for future work is concerned with predicting other voltage stability criteria that were not considered in this research, specifically voltage stability limit. For this purpose, it is required to develop the proposed equations for modeling the relationships between voltage and the WPP reactive power at a given interconnection site. In this research, the mathematical formulation was developed as a polynomial function with the order of 2 and an exponential function for two X/R PCC regions, (X/R PCC < 2 and X/R PCC > 2). As another extension to this work, the proposed formulation can be further developed as a single function, such as a polynomial function with a high rank, which satisfies the whole X/R region and removes the need for dividing the X/R region into two parts. This may enhance the applicability of the mathematical model for estimating V PCC at feeders with X/R PCC around 2. However, the accuracy of such a function must be compared with the equations developed in this work to ensure that the error is not high. Moreover, the proposed model was developed and validated using IEEE distribution system models. Although, the IEEE standard models have been commonly used in electrical engineering studies, the validation of the proposed mathematical relations using actual cases is important. From a practical perspective, the application of the proposed model to the actual cases requires data obtained from real world distribution systems, such as the values of the X/R ratio and SCC measured at the buses of test system, which is not currently available to the authors. Hence, as part of future work, the authors are considering applying the obtained functions to actual cases to further complement this research.