Research on Optimal Access Point and Capacity of Multiscenario Wind Turbine Based on Voltage Sag

Aiming at reducing the impact on sensitive loads and the economic losses caused by voltage sags in the distribution network containing wind turbines, a multiobjective mathematical model of the optimal access point and access capacity is established from the three perspectives of the basic costs of wind turbines, the economic losses of voltage sags, and the economic losses of the power grid. A scene generation method based on the similarity of wind curve changes and a scene reduction algorithm based on DPC & k-means are proposed to solve the mathematical problem of wind turbine uncertainty. Based on the analysis of the voltage sag of the distribution network containing wind turbines, the evaluation indicator of the voltage sag depth of the bus and the evaluation indicator of the economic loss of the voltage sag of the sensitive load are proposed. -e simulation results of the IEEE 33-bus system show that the reasonable access of wind turbines can help reduce the power loss and voltage sag of the distribution network, thereby ensuring the safe and economical operation of the system.


Introduction
Wind turbine has the characteristics of low pollution and strong sustainability compared with traditional thermal power generation. e proportion of wind turbine in the distribution network is gradually increasing. Taking into account the intermittency and uncertainty of wind turbine output and the impact on the distribution network, the selection of the access point and capacity of wind turbines has become one of the research focuses. Most power quality problems are caused by voltage sags, which would make sensitive loads unable to operate normally and cause economic losses. e access of wind turbines has a mitigation effect on voltage sag, so, from the perspective of alleviating voltage sags, the selection of the wind turbine optimal access plan is of great significance.
Regarding the problem of access point and capacity of wind turbines, most researchers solve it from three aspects as follows.
(1) e model of wind turbines: most of the articles use wind turbines as a power source with constant output for modelling. For actual wind turbines, due to the influence of uncertain factors such as wind speed, there are uncertainties in the output of the wind turbine. In order to solve such problems more accurately, the authors in [1] proposed a new voltage stability model based on voltage, wind energy permeability, short-circuit capacity, and short-circuit impedance ratio (X/R), so that engineers can predict the voltage changes at the public access point of wind turbines. e authors in [2] use multiscenario probability models to analyze the time sequence of DG output and load demand. e authors in [3] proposed a new measurement method of wind power change rate and studied the influence of wind power change rate on the access capacity of independent wind power system. (2) Optimization goal: in [4], the loss sensitivity factor (LSF) is used to determine the installation position of DG with the objective of reducing the bus loss of radial distribution network. Intelligent water drop (IWD) algorithm is used to determine the installed capacity of DG. In [5], the optimal capacity of DG is determined by meeting the maximum allowable line capacity limit (MAlCl). e authors in [6] use weighted voltage support capability index (WVSAI) based on sensitivity analysis to quantify the impact of PV-STATCOM on voltage recovery efficiency under voltage sag. e authors in [7] proposed a multienergy power supply multipoint layout planning model based on complex adaptive system theory. e authors in [8] proposed a three-objective optimization planning model for wind turbines and capacitors considering the uncertainty of plug-in vehicle load demand, wind speed, and load. e specific optimization targets are voltage stability indicators, total costs, and greenhouse gas emissions.
(3) Optimization algorithm: e authors in [9] proposed a computer intelligent algorithm based on evolutionary programming (EP) to solve the problem of DG access point and capacity. Among them, the goal is to minimize the power loss of the system and to ensure the voltage status of the entire network bus. e authors in [10] use particle swarm optimization algorithm (PSO) to solve the problem. For the purpose of power system economy, the authors in [11] proposed an optical optimization algorithm (OIO) to improve the current problem. In [12], an improved long-horned beetle search algorithm is used to solve the problem of distribution network planning with DG. Among them, the improvement of the algorithm is based on the Monte Carlo rule in the simulated annealing algorithm. In [13], a fuzzy genetic algorithm (FGA) is proposed, which can maximize the economic benefits of the project.
It can be known that the consideration of wind speed uncertainty in the abovementioned references is only limited to fixed scenarios; the uncertainty of wind speed is not described in terms of actual conditions. In addition, in the actual problems of selection of access point and capacity of wind turbines, most of them determine the optimal access point and capacity from an economic point of view. e influence of wind turbine access on the reliability of the distribution network is not considered.
Wind power is uncertain and difficult to be expressed with actual mathematical models. Multiscenario analysis method is used to simulate the uncertainty of wind power in this paper. Scenes are generated based on the characteristics of wind curve change similarity, and a scene reduction algorithm based on DPC and k-means is used to delete the scene reasonably. In order to evaluate the mitigation degree of wind turbines access to voltage sag caused by fault, an evaluation indicator for the voltage sag depth is proposed. A DG multiobjective optimization model is established using the indicator of voltage sag depth. is model considers the impact of wind turbine access on the reliability of the distribution network while considering economy. An elite sampling algorithm suitable for engineering is used to solve the model, and example simulation was carried out in the IEEE 33-bus system to verify the feasibility and effectiveness of the method. e main contributions of this paper are as follows: (1) is paper proposes a scene reduction method based on the combination of DPC and K-means algorithm. is method reduces the number of scenes according to the demand, thereby simplifying the complexity of calculations.
(2) A multiobjective mathematical model for the optimization of wind turbine access point and capacity selection is established. is model considers the economic loss of voltage sag while considering the operation and maintenance costs of wind turbines and the economic loss of power grid loss. Combining the research methods in this paper can solve the problems related to voltage sag and economy in the distribution network containing wind turbine.
(3) is paper proposes the evaluation indicator of the voltage sag depth of the bus. is indicator is used to calculate the economic loss of sensitive loads caused by voltage sag and to determine the access point and capacity of the wind turbine.
is paper is organized as follows. Section 2 presents a wind power multiscenario analysis method based on wind speed curve change similarity. Section 3 describes the optimization of multitarget capacity and location selection of wind turbine. In Section 4, elite sampling algorithm is used to solve the mathematical model. In Section 5, a simulation based on IEEE 33-bus system is presented and an analysis of the results is also described in Section 5. At last, conclusions are discussed in Section 6.
Nomenclature list is shown in Table 1.

Wind Power Multiscenario Analysis Method Based on Wind Speed Curve Change Similarity
Due to the influence of the random change of wind speed, the output of wind turbine is uncertain. erefore, in the mathematical model of optimal access point and capacity of wind turbine, if the output of generator is regarded as a constant value, there may be a large error between the final calculation result and the actual result, which cannot simulate the real output of wind turbine. erefore, it is necessary to take the time sequence characteristics of wind power output into consideration when selecting optimal access point and capacity of wind turbines. e core content of multiscenario analysis method is to change the uncertain factors into many easily obtained and calculated deterministic factors. e multiscenario analysis method has high adaptability to random problem. erefore, this method is used to simulate the uncertainty of wind power output.

Scene Generation Based on Similarity of Wind Speed Curve
Change.
e traditional scene generation method classifies the wind speed or wind turbine output according to four seasons, but the wind speed data of different regions in the world extracted by HOMER shows that the wind speed will also change with the location and the wind speed changes in different regions in four seasons are not the same. e fixed mode of wind speed division cannot guarantee the correctness of the subsequent mathematical model; the traditional scheme of mathematical model has low adaptive ability. In order to better adapt to the characteristics of wind speed in different regions, this paper proposes a scene generation process based on the similarity of wind speed curve change trend on the basis of [14]. Comparing the daily wind speed changes in 365 days, those with similar wind speed change curves are regarded as one kind. When the number of wind speed curves of a certain type is small, the probability of this occurrence is considered small, and it is deleted from the original database.
rough the above classification process, the process of scene generation is further optimized, and the accuracy of subsequent calculation is improved.
In this paper, the absolute error of enclosing area between different wind speed curves is used to describe the trend similarity of curves [15]. When the trend similarity of two curves is smaller, the change trend of the two curves is closer. When the two curves coincide, the curve trend similarity is 0.
If there are two wind speed curves w1(t) and w2(t), the trend similarity of the two curves is CS, and the specific calculation formula is as follows.
(1) e optimized scene generation process is as follows.
Step 1: combined with the geographical and meteorological data, the wind speed changes of a certain area in each period of the whole day in a year are obtained through the HOMER.
Step 2: calculate the trend similarity between curves according to formula (1) and classify the wind speed according to the principle of curve trend similarity.
Step 3: assume that the error ΔD between the actual data and the typical data of each period obeys the corresponding probability distribution function (PDF). e wind speed data that have been classified according to the trend similarity principle are extracted. e probability density function PDF (t � 1,2, . . ., 23,24) is calculated by the error between the typical output value and the actual value of each type of wind speed curve. e random variable expression corresponding to each scene can be expressed by the following formulas: where wv x,t,s is the x-th wind speed actual value in scene s at moment t. wv o x,t is the x-th wind speed typical value in scene s at moment t. Δwv x,t,s is the difference between the x-th actual value and the x-th typical value in scene s at moment t. L x,t,s is the x-th load actual value in scene s at moment t. L o x,t is the x-th load typical value in scene s at moment t. ΔL x,t,s is the difference between the x-th actual value and the x-th typical value in scene s at moment t.
Step 4: divide the PDF generated by Step 3 into 7 equal parts, and the schematic diagram is shown as follows. In Figure 1, α i,t represents the probability of occurrence in the error interval of the wind speed i at moment t. Trace of the between-class deviation matrix Tr(S w ) Trace of the within-class deviation matrix Voltage matrix when x-type short-circuit fault occurs Voltage major axis matrix when x-type short-circuit fault occurs P (x) Voltage minor axis matrix when x-type short-circuit fault occurs M (x) Voltage sag matrix when x-type short-circuit fault occurs R (x) Voltage sag depth matrix when x-type short-circuit fault occurs λ x Failure probability of x-type short-circuit fault n e total number of buses F (x) j e depth indicator when x-type short-circuit fault occurs F j e average depth of voltage sag f sen e total economic loss of sensitive load caused by voltage sag f 1 e total cost of wind turbine accessing f i,rm e annual operation and maintenance cost of the wind turbine accessed to the node i f i,fault e failure cost of the wind turbine accessed to the node i P loss e network loss f 2 e total economic cost caused by network loss after considering multiscene Step 5: in order to facilitate the selection of scene, the interval in Step 4 is probability standardized.
Step 6: use Monte Carlo sampling/roulette selection method [16] to obtain the scene s.
(4) e scene s is a parameter consisting of 0 and 1. W x i,t,s , respectively, represents the selection of the error interval of the wind speed i in the scene s at the moment t.
Step 7: calculate the probability of scene s and standardize it.
where w(s) is the occurrence probability of each scene s.
where P r (s) is the normalized probability for each scene.

Scene Reduction Based on DPC and K-Means.
Generally, the number of the divided scenes should be large enough to meet the requirements of calculation accuracy, but the large number of scenes will inevitably increase the calculation complexity of the model. In order to solve the contradiction between calculation accuracy and calculation complexity, it is necessary to limit the number of scenes. Clustering by fast search and finding density peaks (DPC) [17] is a new clustering algorithm proposed by Alex and Alessandro in 2014. is algorithm can find the density peaks of arbitrary shape data sets, select the number of clusters according to the decision graph, and has an efficient sample point allocation strategy [18]. e advantage of this method is that it can select the cluster center of discrete points better, but the disadvantage is that it takes a long time to complete the whole process of clustering. K-means clustering algorithm is a traditional clustering algorithm based on iteration. It is widely used because of its simple principle, easiness to write algorithm, and good clustering effect. e main parameters of K-means algorithm are the number of initial clustering centers and the location of the centers. e performance of K-means clustering depends on the selection of the above parameters. By comparing the pros and cons of the two methods, this paper proposes a method based on the combination of DPC algorithm and K-means algorithm to optimize the clustering process.
Step 1: import the N s original scene data into the DPC algorithm to obtain the initial cluster center point set, namely, the preliminary cluster center scene set.
Step 2: import the initial cluster center scene set obtained in the first step into the k-means algorithm as the parameters for the first iteration. Finally, through the clustering algorithm based on the combination of DPC & K-means, we can get the central scene set and the number of scenes after clustering.
Step 3: assign the original scene data to the corresponding cluster center field according to the distance between the original scene data and the cluster center field and label the original scene by category, so as to complete the first clustering. e probability of the scene represented by the cluster center is calculated by the cumulative probability of all such scenes.
Step 4: considering the characteristics of clustering center in k-means algorithm, it is necessary to iterate the clustering center for many times. erefore, Step 2 and Step 3 are repeated until the final result meets the accuracy requirements or reaches the maximum number of iterations. e corresponding cluster centers are obtained through the above algorithm, but the quality of the results cannot be intuitively evaluated. In order to reflect the pros and cons of the selected clustering centers, this paper uses clustering validity indicator CH (+) [19] to determine the results of the scene reduction algorithm, so as to further improve the accuracy of the scene reduction algorithm. e cluster validity indicator CH (+) is the ratio of within cluster deviation to intercluster deviation, and its expression is shown in the following formula: In formula (7), N represents the number of original data, K represents the number of clusters, Tr(S B ) represents the trace of interclass deviation matrix, and Tr(S w ) represents the trace of intraclass deviation matrix. e calculation of Tr(S B ) and Tr(S w ) is shown in the following formulas: In formula (8), n i represents the number of data i and v j represents the center point of data j. v represents the central point of the entire data set; d(v j , v) is the distance between v j and v.  Mathematical Problems in Engineering In formula (9), x j represents the data j of the class I and v i represents the center point of data i. d(x j , v i ) is the distance between x j and v i . rough the analysis of the above indicators, it can be seen that the denominator of this indicator measures the tightness within the class and the numerator measures the separation between classes. erefore, it can be known that when the indicator CH (+) is larger, the data of the same classes are closer and the data of different classes are more scattered; that is, the clustering effect of the algorithm is better. When the number of clusters is x, its indicator CH (+) is the largest, which can prove that x is the best number of clusters.

Analysis and Evaluation of Voltage Sag Based on Voltage
Ellipse Parameters. e Clarke transformation was firstly proposed by Edith Clarke in 1951 [20]. e three-phase voltage can form space rotation vector E → (t) through Clarke transformation. e space rotation vector can be expressed as shown in Among them, E + � |E + |e jψ + E − � |E − |e jψ − . If the parameter c is defined as major axis and the parameter p is minor axis, the formula can be defined as follows.
According to the above analysis, the corresponding voltage ellipse parameters can be obtained by plotting the space rotation vector E → (t) in a complex plane. e voltage sag caused by different fault types can be expressed by corresponding voltage ellipse parameters [21,22]. erefore, voltage ellipse parameters can be used to classify and characterize voltage sag.
Compared with the voltage sag caused by line fault, the voltage sag caused by bus fault will have a greater adverse impact on the distribution network. erefore, this paper mainly discusses the impact of bus fault on the distribution network. Assuming that each bus in the distribution network has four types of short-circuit faults (all metallic shortcircuit faults), the voltage value of each bus is calculated by the method of fault point [23] and the voltage matrix is established under the given grid structure (where x � 1, 2, 3, and 4 represent single-phase grounding, two-phase short-circuit, two-phase grounding, and three-phase short-circuit, respectively).
U (x) in formula (12) is the voltage matrix under shortcircuit type x, which u (x) ij represents the three-phase voltage on bus j when fault type x occurs on bus i. Clarke transform is applied to the three-phase voltage of each bus to extract the major axis parameters and minor axis parameters of the corresponding bus.
e P (x) and C (x) in formulas (13) and (14) represent the voltage major axis matrix and the voltage minor axis matrix under the short-circuit type x, respectively. e p (x) ij and c (x) ij represent the major axis parameters and minor axis parameters of the bus j when the bus i has short-circuit fault type x.
According to the voltage ellipse parameters, the specific rules for judging whether there is voltage sag on the bus are as follows [21]: (1) When the major voltage axis parameters p (x) ij and minor voltage axis parameters c (x) ij of the bus are in the range [0.9333,1.06667], the bus is normal; that is, there is no voltage sag.
(2) When the major axis parameters p (x) ij and minor axis parameters c (x) ij of the bus are not in the range of [0.9333,1.06667], and the minor axis parameters c (x) ij are less than 0.9333, the voltage sag occurs on the bus. rough the above analysis, the voltage sag matrix with short-circuit type x is obtained.

Mathematical Problems in Engineering
e M (x) in formula (15) is the voltage sag matrix under short-circuit type x. m (x) ij represents the situation of bus j when short-circuit faults type x occurs on bus i. If voltage sag occurs on the bus j, let m (x) ij � 1. On the contrary, if there is no voltage sag on the bus j, let m (x) ij � 0. e voltage sag depth matrix R (x) can be obtained through the voltage matrix U (x) and the voltage sag matrix M (x) , the formula can be defined as follows.
In formula (16), u (x) ij,min represents the phase voltage with the smallest voltage value on bus j when an x-type fault occurs on bus i.

Analysis and Evaluation of Economic Loss Caused by
Voltage Sag of Sensitive Load. From the column j of the voltage sag depth matrix R (x) , the occurrence of voltage sags on bus j when other bus failures occur can be obtained.
erefore, from the column information, the depth indicator F (x) j of the voltage sag caused by the fault type m on the bus j by other buses can be obtained; see the following formula for details: Consider the average depth F j of voltage sag when four short-circuit faults occur on the bus j in turn. See the following formula for details: In order to measure the economic loss of sensitive load caused by voltage sag, this paper proposes the economic loss indicator of sensitive load based on depth indicators F (x) j and F j . (20) f sen is the total economic loss of sensitive load caused by voltage sag. τ j is the average economic loss of sensitive load on node j caused by a single voltage sag event. χ is the total number of nodes accessed to sensitive loads.
In order to further reflect the impact of wind turbine access on the economic loss of sensitive load voltage sag, the influence of wind power multiscenario is considered in formula (20), and then the calculation formula of sensitive load voltage sag economic loss is changed to the following formula: In formula (21), f sen,m,s,t is the economic loss of sensitive load voltage sag in scene s at moment t corresponding to the wind speed type m considering the time sequence of wind turbine output.

Objective Function.
e objective function in this paper includes the basic cost of wind turbine, the economic loss of voltage sag, and the economic loss of power grid.

Basic Cost of Wind Turbine.
e basic cost of wind turbine is mainly composed of three parts. One is the onetime investment and installation cost of wind turbine. e second is the operation and maintenance cost of wind turbine in the later stage. e third is the cost of failure cost. e establishment of the third part is to consider the cost incurring when the wind turbine is not put into the distribution network and the actual wind speed is lower than the cut-in wind speed of the wind turbine.
In formula (22), f 1 is the total cost of wind turbine accessing. f i,rm is the annual operation and maintenance cost of the wind turbine accessed to the node i, and the specific calculation formula is shown in formula (23). f i,inv is the one-time investment and installation cost of the wind turbine accessed to the node i. f i,fault is the failure cost of the wind turbine accessed to the node i, and the specific calculation formula is shown in (24). N DG is the total number of wind turbines accessed to the distribution network. NM is the total number of wind speed categories; NS m is the total number of scenes under wind speed type m; d w,m is the number of days under wind speed type m. f m,s,t,rm is the operation and maintenance cost of the wind turbine under the scene s of the wind speed type m at moment t. Considering that when the wind speed is lower than the cut-in wind speed, the wind turbine will exit the distribution network. At this time, the wind turbine cannot be used as a power source to supply power to the distribution network.
erefore, let f m,s,t,fault be the failure cost of the wind turbine under the scene s of the wind speed type m at moment t. P rw (s) is the normalized probability for wind speed scene s.

Economic Loss of Distribution Network.
Under a certain access scheme, the economic loss of power grid loss can be calculated by the following formulas: where P loss is the network loss, U i , U j are the voltage of node i and j, Y ij is the admittance of branch i-j, θ ij is the impedance angle of branch i-j, and δ ji is the phase angle difference of node voltage.
Among them, f 2 is the economic cost caused by network loss and δ is the economic cost caused by unit network loss.
After considering the influence of scene on network loss, the corresponding economic loss formula of network loss is calculated by the following formulas: (27) In formula (27), f 2 is the total economic cost caused by network loss after considering multiscene. Among them, LN is the total number of load curve categories and LC i is the total number of scenes under the load type m; d m is the number of days under the wind speed type m. f 2,m,s,t is the economic loss of network loss at the moment t in the scene s corresponding to the wind speed type m. In formula (28), f 2,m,s,i,c,t is the economic loss of network loss at moment t in scene s corresponding to wind speed type m and scene c corresponding to load type i. P L (c) is the standardized probability for load scene c. d l,m is the number of days under the load m.

Economic Loss of Sensitive Load Caused by Voltage Sag.
Under a certain access scheme, after considering the influence of multiscenes, the calculation formula of voltage sag economic loss of sensitive load is as follows: In formula (29), f sen,m,s,t is the economic loss of sensitive load voltage sag at the moment t corresponding to the scene s of the wind speed type m.

Constraint Condition
(1) Voltage variable inequality constraints: All node voltages must meet the In formula (30), U i min is the lower limit of allowable voltage value of node i. U i max is the upper limit of allowable voltage value of node i.
(2) Branch current constraint: In formula (31), P ij and P ij, max are, respectively, transmission power and the maximum allowed transmission power between node i and node j. (3) e total capacity constraint of the wind turbine installed in the system is as follows: Generally, the total installed capacity of wind turbines in distribution network should not exceed 10%∼25% of the maximum capacity of power grid. erefore, in formula (32), S i is the rated active power of the wind turbine accessed to the node i. η is the capacity penetration of distributed generation. S is the total rated active power of the load accessed to the distribution network. (4) e rated capacity constraints of wind turbine accessed to single node are as follows: In formula (33), P i is the rated active power of the wind turbine accessed to the node i. P i,max is the maximum rated active power of the wind turbine that the node i is allowed to access. e mathematical model of capacity and location selection of wind turbine optimization problem is to minimize the basic cost, network loss, and voltage sag economic loss of wind turbine under the above constraints.
Among them, [f 1 , f 2 , f sen ] is the objective function of basic cost, economic cost of network loss, and economic loss of voltage sag. g(x) and h(x) are inequality constrained function and equality constrained function, respectively.

Elite Sampling Algorithm
In this paper, elite sampling algorithm is used to solve the problem of optimal access point and capacity of wind turbine. Its main features are simplicity and easiness to understand, small amount of calculation, and high applicability [24]. So, it can be better used in optimization problems. In the elite sampling algorithm, when the number of solutions is different, the specific steps of the algorithm and the relevant parameters in the sampling process will also change. In this paper, the number of solutions is limited, so we only introduce the elite sampling algorithm when the number of solutions is finite.

e Basic Principle of Elite Sampling Algorithm.
Assuming that the total number of individuals in the whole solution space is M, all the solutions in the whole solution space are sorted by performance and then assume that the first n solution in the performance ranking is the elite solution. erefore, after generating N solutions randomly, the probability of not being in the first n solution is Suppose that the probability that none of the randomly generated N candidate solutions are in the first n is not greater than q%; then there is On the basis of formula (35), the minimum number N min ′ of samples needed for the problem can be obtained by iterative calculation. As long as the number of randomly generated candidate solutions is not less than N min ′ , there is a possibility of (1− q%) among the N min ′ candidate solutions, at least one of which is the first n satisfactory solution. e specific steps are as follows: Step 1: initialize parameter M and n, N � 1, t � 1 Step 2: calculate t � t(1 − n/M + 1 − N) Step 3: judge whether t ≤ q% is true? If yes, continue with Step 4; otherwise, N � N + 1; return to Step 2.
Step 4: N min ′ � N In order to improve its efficiency, parallel algorithm can be used in the process of analyzing and evaluating each candidate solution. In actual engineering practice, it is usually enough to obtain elite solutions with better effects and economic efficiency for solving practical problems.

Performance Ranking Selection Method.
e solution of multitarget optimization problem is not unique; it is generally a solution set, so the performance evaluation of candidate solutions is particularly critical.
ere are currently three typical traditional evaluation methods for the solution of multiobjective optimization problems: linear weighting method, reference target method, and weighted composite method. e disadvantage of these methods is that the weight setting is highly subjective and has a great influence on the results. In this paper, the performance ranking selection method is used to avoid the subjectivity of the weight setting of the traditional evaluation method, so as to make the problem solving more reasonable.
Before selecting candidate solutions, the performance ranking selection method needs to sort all candidate solutions under a single performance indicator to obtain the optimal solution under a single indicator. On the basis of sorting all performance indicators individually, the performance ranking selection is carried out. e solution with the worst ranking is selected. If there are multiple candidate solutions with the same worst ranking, the second worst ranking is considered. If the second worst ranking is also the same, the third worst ranking is considered, and so on, until the elite solution is selected.

Solution Strategy of Mathematical Model for Optimal
Access Point and Capacity of Wind Turbine. In this paper, the optimization of wind turbine capacity and location is an optimization problem with discrete control variables.
e specific solving steps are as follows.
Step 1: taking the access node and capacity of the wind turbine as variables, the solution space of the access scheme of the wind turbine to be optimized is determined. Considering that the total installed capacity of wind turbines is limited by the maximum capacity of the grid, this paper limits the total number of wind turbines in the early stage. Using p% elitist sampling method, the minimum number of samples N min ′ is obtained for the mathematical model proposed in this paper Step 2: randomly generate the wind turbine access scheme. e sampling strategy of each scheme is to randomly select the control variables within their variation range. e generated wind turbine access scheme needs to judge whether the scheme is the same as the existing scheme, and if it is the same, reextract it.
Step 3: calculate the power flow of the system corresponding to different randomly generated wind turbine access schemes, check the constraints such as electrical limit, and call the wind turbine access scheme that meets the constraints as the feasible wind turbine access scheme. e nonconforming schemes are directly deleted and resampled. Finally, based on the above N min ′ feasible wind turbine access schemes, the network loss under the corresponding scheme is calculated.
Step 4: combined with different wind turbine scenes, use formulas (21) and (22) to calculate the investment cost and operation and maintenance cost of each feasible wind turbine access scheme.
Step 5: combined with different wind turbine scenes and load scenes, the short-circuit calculation is carried out for the systems corresponding to the N min ′ feasible wind turbine access schemes. and then the economic loss of sensitive load voltage sag corresponding to each wind turbine access scheme is calculated by using formula (29). 8 Mathematical Problems in Engineering Step 6: calculate the corresponding objective function value of each feasible wind turbine access scheme, select the final elite solution according to the performance ranking selection method, and output the access point and access capacity of the optimal access scheme.
In summary, the flow chart of the research on the multitarget optimal access point and capacity of wind turbines is shown in Figure 2.

Example Analysis
e wind speed data used in this paper are from the simulation of HOMER. e wind speed data includes the wind speed of 24 hours in 365 days, and the wind speed of each hour is regarded as a constant value. After classifying the trend similarity of the wind speed curve, the wind speed curve is divided into 9 categories and a total of 328 wind speed data. e trend change of part of the first type wind speed curve obtained is shown in Figure 3.
Follow the scene generation steps described in Section 2.1 to get 500 sets of wind speed scenes. Because there are too many scenes, only the normalized probabilities of some scenes under the first wind speed curve are listed. See Table 2 for details.
If the wind speed curves are divided into m type, for each type of wind speed curve, it will generate 500 wind speed scenes randomly, and there will be 500 * n wind speed scenes. e increase of the number of scenes is helpful to the real response to the original situation and the accuracy of subsequent calculation. But, at the same time, it will also increase the amount of calculation in the subsequent calculation process, which is not enough to meet the needs of engineering. erefore, it is necessary to delete 500 sets of wind speed scenes under each type of wind speed in order to ensure a small amount of calculation and enough wind speed scenes. e final scene deletion results should be able to correctly and reasonably reflect the actual wind speed curve changes, and the number of scenes can be guaranteed in a reasonable range.
In order to ensure the validity of clustering analysis in the process of scene deletion, the actual number of clustering based on the combination of DPC & K-means is set to be 2∼20, and 19 clustering schemes are obtained. Select the number of clustering schemes according to Figure 4. e specific results are shown in the figure as follows.
It can be seen from Figure 4 that when the number of clusters is 2, the indicator CH (+) of DPC & k-means algorithm and K-means algorithm is the largest, which indicates that 2 is the optimal scene deletion result under this kind of wind speed. Moreover, when the number of clusters is 2, the indicator CH (+) obtained by DPC & k-means algorithm is larger than that obtained by K-means algorithm. It shows that DPC & k-means algorithm can cluster data better and get better results.
After the scene reduction process, we need to calculate the corresponding scene probability. e scene probability after scene reduction under the first wind speed curve is listed in Table 3. e optimization of optimal access point and capacity of wind turbine based on IEEE 33-bus system consists of 33 nodes, 32 branches, and 5 tie switch branches and operates in radial network. e five dashed lines in Figure 5 are the branches of network interconnection switch, and the Mathematical Problems in Engineering 9 switches on them are normally open.
e distribution system structure and line parameters in this example are from [25]. e inverter wind turbine is used for simulation verification. e inverter wind turbine adopts PV control, and the power factor is about 0.9. ere are 10,11,12,13,14,15,16,17,18,19, and 20 buses that can be accessed to the wind turbine, but only 4 buses are allowed to access the wind turbine.
e rated capacity of a single wind turbine is 100 kW. e installation cost of a single wind turbine is 100,000 yuan (CNY), and the hourly operation and maintenance cost of the equipment are 400 yuan. e access point of general load is 16, and the general load has time sequence, and the data is automatically generated by the HOMER. Sensitive load nodes are 11 and 25, and the capacity of sensitive load does not change with time.
e average economic loss caused by a single voltage sag event of sensitive load on bus 11 and bus 25 is 100,000 yuan (CNY) and 300,000 yuan (CNY), respectively. e economic cost of unit network loss is 1,000 yuan. e fault probabilities of single-     According to the number of wind turbines and the change of access points, different access schemes of wind turbines are constructed. Considering that the access capacity of wind turbine is limited by the total load capacity, the number of single bus accessed to wind turbine is 8. erefore, the number of wind turbines that can be accessed to each bus is between 0 and 8, and the value is randomly selected by Step 1. Finally, the total number of wind turbine access schemes (the total number of individuals in the solution space) can be calculated as M � (10 × 9) 4 , and the minimum number of samples N min ′ � 459 .
ere are many data about the comparison and ranking of the solutions obtained under different indicators, and only 20 groups of solutions are extracted. See Table 3 for details.

Conclusions
e access of wind turbines to the distribution network has a mitigation effect on voltage sag. is paper considers the time sequence characteristics of wind power output, and, from the perspective of alleviating voltage sags, an optimization algorithm for the selection of the access point and capacity of wind turbines based on voltage sags is proposed. Several detailed conclusions are obtained: (1) e multiscene analysis method is used to simulate the uncertainty of wind turbine output; the scenes are reasonably generated and deleted. It converts uncertain factors into deterministic factors that are easy to be obtained and calculated and solves practical problems on the basis of deterministic factors. is method effectively avoids mathematical modelling of high-complexity random problems and reduces the difficulty of solving problems. (2) e evaluation indicator of the voltage sag depth of the bus and economic loss evaluation indicator for sensitive load voltage sag are proposed. ese indicators can evaluate the mitigation degree of wind turbine access to voltage sag caused by fault. ese indicators provide help and basis for the reasonable selection of the access point and capacity of wind turbines.
(3) A multiobjective mathematical model for the selection of the access point and capacity of wind turbines is proposed in this paper. is model considers the traditional basic cost of wind turbine and the economic loss of network losses and also considers the voltage sag problem that has the most complaints in power quality. It provides a new idea for the reasonable selection of the access point and capacity of wind turbines. (4) e elite sampling algorithm is used to solve the multiobjective mathematical model proposed in this paper. It can obtain a more high-precision solution based on the reasonable solution space and avoid the situation of nonconvergence. is algorithm has the characteristics of simplicity and easiness to understand, small amount of calculation, and high applicability, and it is suitable for solving practical engineering problems. (5) e IEEE 33-bus system example is used to verify the optimization algorithm proposed in this paper. e results show that the wind turbine accessing scheme determined by the algorithm can not only reduce the network loss but also alleviate the voltage sag.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.