Front. Energy Res. Frontiers in Energy Research Front. Energy Res. 2296-598X Frontiers Media S.A. 1097185 10.3389/fenrg.2022.1097185 Energy Research Original Research Equivalent modeling and multi-parameter coupling optimization for DFIG-based wind farms considering SSO mode Zhang et al. 10.3389/fenrg.2022.1097185 Zhang Qian 1 * Jin Xuyang 1 Zhang Fan 2 Yuan Hao 2 Zhou Baorong 2 1 State Key Laboratory of Power Transmission Equipment and System Security and New Technology, Chongqing University, Chongqing, China 2 State Key Laboratory of HVDC, Electric Power Research Institute CSG, Guangzhou, China

Edited by: Wei Gan, Cardiff University, United Kingdom

Reviewed by: Hui Hou, Wuhan University of Technology, China

Yue Xiang, Sichuan University, China

 *Correspondence: Qian Zhang, zhangqian@cqu.edu.cn

This article was submitted to Process and Energy Systems Engineering, a section of the journal Frontiers in Energy Research

 18 01 2023 2022 10 1097185 13 11 2022 17 11 2022 Copyright © 2023 Zhang, Jin, Zhang, Yuan and Zhou. 2023 Zhang, Jin, Zhang, Yuan and Zhou

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

As a low-carbon and environmentally friendly renewable energy source, wind power has been globally recognized as the best solution to achieve energy saving and emission reduction and promote low-carbon economic growth. With the increase of wind power penetration, wind power has a great impact on sub-synchronous state stability and dynamic characteristics of the grid-connected system. Aiming at the fact that the correlation between clustering indexes and sub-synchronous oscillation (SSO) mode and the difference of the contribution to the clustering results are seldom considered in the current equivalent modeling of doubly-fed induction generator (DFIG)-based wind farm, this paper proposes a clustering method based on the index dimension reduction and weighted fuzzy C-means (WFCM) clustering algorithm. Besides, for the SSO study of the grid-connected system without sufficiently considering the coupling effects between controller parameters, a multi-parameter coupling optimization design strategy combining orthogonal experiment method (OEM) and response surface method is proposed. Firstly, the dominant variables of SSO mode of the DFIG-based wind farm connected to weak grid by series compensation system are taken as the initial clustering indexes. After dimension reduction by principal component analysis, the WFCM algorithm is utilized to cluster the wind farm. Then, the proportional and integral coefficients of the grid-side controller, rotor-side controller and phase-locked loop are optimized to achieve the simultaneous optimization of the SSO characteristics and dynamic characteristics of the system. Finally, the interaction between control parameters and the influence degree and trend on the system performance are quantitatively evaluated, and the optimal parameter combination is obtained. The proposed strategy can mitigate SSO more effectively while improving anti-interference than the particle swarm optimization based on OEM.

 decarbonization sub-synchronous oscillation (SSO) equivalent modeling multi-parameter coupling optimization DFIG-based wind farm 1 Introduction

Building a new power system with new energy as the main body and non-carbon energy as the carrier, and constructing a low-carbon and high-efficiency energy system are important ways to accelerate the construction of ecological civilization and realize the harmonious coexistence between human and nature (Huang et al., 2022). As a low-carbon and environmentally friendly renewable energy source, wind power has become an effective way to alleviate energy scarcity, combat climate change, and promote low-carbon economic growth worldwide (Xia et al., 2020). Doubly-fed induction generator (DFIG)-based wind farms are generally located in remote areas far from load centers. The strength of the grid is usually classified according to the short-circuit ratio (SCR). Generally, the grid with the value greater than 20 is considered as the strong grid and the value less than six is weak grid (Miao, 2012). As the grid-connected capacity of DFIGs continues to increase, the interaction between grid impedance and the equivalent impedance of DFIG-based wind farm can cause oscillations with weak grid operating condition. In addition, the series compensation capacitor will lead to the risk of sub-synchronous oscillation (SSO) (Liu et al., 2020; Leon et al., 2021). With the increasing penetration of wind power, equivalent modeling of wind farms has attracted attention at home and abroad. If each wind turbine is modeled in detail, it will not only increase the complexity of the power system model, but also be highly prone to the problem of dimensional disaster, which cannot reproduce SSO correctly (Han et al., 2020). Therefore, it is necessary to study dynamic equivalent modeling methods for wind farms to simplify the complexity of the model and reduce the computational effort.

Wind farm dynamic equivalence methods are generally divided into the single-machine equivalence method and the multi-machine equivalence method. The single-machine equivalence method equates all wind turbines into one unit, and the equivalent parameters are usually obtained by the weighted summation method (Pan et al., 2015). However, the actual operating conditions of the units within a large wind farm vary due to factors such as topography, wake effect, and time lag. When studying the dynamic characteristics of wind power grid-connected systems, the accuracy of the single-machine equivalence model cannot meet the actual demand. The multi-machine equivalence method clusters wind farms into multiple groups based on the principle that the operating conditions of wind turbines are the same or similar, and equates the wind turbine parameters of the same cluster, and finally uses several units to characterize the wind farm (Lei et al., 2013). The core idea of the multi-machine equivalence is to find the attribute that characterizes the key features of DFIG and use this attribute for clustering. In (Shabanikia et al., 2021), a weighted dynamic aggregation method is proposed to obtain equivalent turbines for the problem of wind speed variations in different regions and different turbine parameters in wind farm aggregation modeling. A dynamic equivalence method for DFIG-based wind farms based on density-peak clustering algorithm is proposed in (Liu et al., 2019), which performs a two-stage clustering of wind farms based on the protection state of the crowbar and the transient characteristics of each unit under faults. In (Zhou et al., 2019), variables such as steady-state wind speed, slip rate, active power output, and rotor current are selected as indicators and Gaussian density distance clustering algorithm is used for clustering.

However, the current equivalence model does not consider the correlation between clustering indexes and SSO mode. It is inconclusive whether it is applicable to the analysis of SSO characteristics of DFIG-based wind farm connected to weak grid by series compensation systems. Moreover, the data set consisting of multiple DFIGs and clustering indexes may contain strongly correlated variables, making the effect of a particular index magnified and resulting in a bias in clustering results. Because of high dimensionality, the data is complicated and redundant, which leads to a long time spent in clustering. Most of the clustering algorithms currently used focus on the optimization of clustering centers and the determination of the number of clusters, ignoring the differences and correlations between clustering indexes.

The control system of DFIG consists of several PI controllers and the effects of different control parameters on the system performance vary greatly. The SSO can be effectively mitigated by reasonably adjusting control parameters to optimize the overall dynamic performance of the system. The main methods for optimizing DFIG control parameters are engineering tuning (Wanigasekara et al., 2022), small-signal analysis (Hu et al., 2022), and intelligent algorithm (Han et al., 2021). The engineering tuning relies too much on experience and there is “over-adjustment” or “blind adjustment” in the process of parameter adjustment; the small-signal analysis is not universal and is greatly influenced by the system operation state; the intelligent algorithm is difficult to obtain the optimization range and the global optimum cannot be satisfied. In (Bakir et al., 2020), the offline optimization of PI controllers at the rotor side of the DFIG wind energy system is achieved by bacterial foraging algorithm. In (Chen et al., 2019), based on participation factors and small signal model eigenvalues, the system damping is increased by exploring the optimal range of PI parameters. Most of the present studies on the optimization of DFIG parameters adopt intelligent algorithms. There is a lack of analysis of the interaction between parameters, and it is easy to fall into local optimum due to the lack of clear optimization direction.

To address the above problems, the dominant variables of SSO mode of the DFIG-based wind farm connected to weak grid by series compensation system are taken as the initial clustering indexes. Considering the adaptation to dynamic and transient process analysis under multi-disturbance scenarios, a clustering method for DFIG-based wind farm based on the principal component analysis (PCA) and weighted fuzzy C-means (WFCM) clustering algorithm is proposed. The main contributions are as follows:

1) The proposed clustering method solves the problem that the correlation between clustering indexes and SSO mode, the contribution degree of different clustering indexes to the clustering results, and the inaccurate clustering caused by data redundancy are seldom considered in the equivalent modeling of the wind farm in the current study of the DFIG-based wind farm connected to weak grid by series compensation system.

2) A multi-parameter coupling optimization strategy combining orthogonal experiment method (OEM) and response surface method (RSM) is proposed to quantitatively analyze the interaction between each control parameter and the sensitivity and influence trend on the system performance, which solves the problem that the intelligent algorithm is easy to fall into local optimum.

3) Aiming at the simultaneous optimization of the sub-synchronous state stability and dynamic characteristics of the system, the control parameters of the clustered equivalent units are optimized. The optimal parameter combination is obtained, and the effective mitigation of SSO and the improvement of system anti-interference are achieved.

The rest of this paper is organized as follows: Section 2 introduces the DFIG-based wind farm clustering method based on the PCA and WFCM algorithm. The multi-parameter coupling optimization strategy based on the OEM-RSM is explained in Section 3. Section 4 gives the results of equivalence modeling and parameter optimization of the wind farm, and analyzes the interaction between control parameters and the influence degree on the performance of the grid-connected system. Conclusions and outlooks are drawn in Section 5.

 2 Sub-synchronous oscillation equivalence of DFIG-based wind farm 2.1 Dimension reduction of clustering indexes based on PCA 2.1.1 Selection of clustering indexes for DFIG-based wind farm

In the DFIG-based wind farm connected to weak grid by series compensation system, the analysis of the SSO mechanism shows that the phase-locked loop (PLL) phase angle may affect the output voltages of grid-side controller (GSC) and rotor-side controller (RSC) and thus the output current of DFIG leading to system destabilization oscillation. And the fluctuation of DC-link voltage may affect the output current of DFIG and thus the grid voltage through GSC and RSC output voltages (Dong et al., 2016). Also, the grid strength, series compensation degree (SCD), DFIG d-axis and q-axis stator current and d-axis and q-axis rotor current are the main influencing factors of SSO (Liu et al., 2017). In order to better analyze SSO characteristics of the system, improve the effectiveness of clustering, and then design more targeted mitigation measures, this paper takes the dominant variables of SSO mode as clustering indexes: grid current angular frequency ω s , DC-link voltage U d c , d-axis stator current i s d , q-axis stator current i s q , d-axis rotor current i r d , q-axis rotor i r q .

The initial operating point of DFIG-based wind farm is an important basis for characterizing its steady-state properties. Different initial operating points affect the dynamic response of the wind turbine during SSO, resulting in differences in the state. Considering the adaptation to the analysis of dynamic and transient processes under multi-disturbance scenarios, this paper selects one time point of data before DFIG-based wind farm is connected to weak grid by series compensation ( ω s 1 , U d c 1 , i s d 1 , i s q 1 , i r d 1 , i r q 1 ), one time point of data during DFIG-based wind farm is connected to weak grid by series compensation ( ω s 2 , U d c 2 , i s d 2 , i s q 2 , i r d 2 , i r q 2 ), one time point of data during the fault ( ω s 3 , U d c 3 , i s d 3 , i s q 3 , i r d 3 , i r q 3 ), and one time point of data after fault removal ( ω s 4 , U d c 4 , i s d 4 , i s q 4 , i r d 4 , i r q 4 ). A total of 24 indexes are taken as initial clustering indexes.

 2.1.2 Dimension reduction by PCA

The 24 initial clustering indexes may contain strongly correlated variables, which makes the weight and influence of a certain index magnified, resulting in a tendency of clustering results and affecting the accuracy of clustering. Moreover, the dimension of the data set is high and the data are complicated and redundant, so it takes longer time for clustering. Therefore, before using the initial clustering index to cluster DFIGs, it is necessary to process the data to eliminate the strong correlation and redundancy among them and reduce the negative impact on the clustering results caused by the interference of the data itself.

As a dimension reduction method for unsupervised learning, PCA is mainly used to reduce the dimension of dataset by measuring the amount of information through variance to eliminate the interaction between the original data (Sun et al., 2021). The magnitude of carrying the original feature variance information can be characterized by the contribution rate of principal component (PC). A larger contribution rate indicates that the PC is more capable of explaining the original feature information. In this paper, PCA is used to reduce the dimension of 24 initial clustering indexes in five steps (Li et al., 2019):

1) Build the sample matrix of the original data. The initial clustering index data extracted from the wind farm are constructed into a sample matrix X w as shown in Eq. 1. Where n is the number of samples, i.e., the number of DFIGs in the wind farm; p is the number of indexes in each sample, and p is 24.

X w = x 11 x 12 ⋯ x 1 p x 21 x 22 ⋯ x 2 p ⋮ ⋮ ⋱ ⋮ x n 1 x n 2 ⋯ x n p = x 1 , x 2 , … , x p

2) Calculate the sample covariance matrix. The data in X w can be normalized by calculating the mean and standard deviation of the sample.

 x i j * = x i j − x j ¯ s j 2 i = 1,2 , … , n , j = 1,2 , … , p x j ¯ = 1 n ∑ i = 1 n x i j s j 2 = 1 n − 1 ∑ i = 1 n x i j − x j ¯ 2 Where x i j is the original sample data of the j-th index of the i-th DFIG; x i j * is the normalized data of x i j ; x j ¯ and s j 2 are the mean and standard deviation of the j-th index, respectively. X w is normalized to X w * and the covariance matrix R w can be obtained by Eq. 5. The correlation coefficient r e j between x e * and x j * in X w * can be calculated from Eq. 6, where x e * ¯ and x j * ¯ are the mean values of the e-th and j-th indexes; x i e * and x i j * are the normalized data of the e-th and j-th indexes of the i-th DFIG, respectively. R w = 1 n X w * T X w * = r 11 r 12 ⋯ r 1 p r 21 r 22 ⋯ r 2 p ⋮ ⋮ ⋱ ⋮ r p 1 r p 2 ⋯ r p p r e j = ∑ i = 1 n x i e * − x e * ¯ x i j * − x j * ¯ ∑ i = 1 n x i e * − x e * ¯ 2 ∑ i = 1 n x i j * − x j * ¯ 2 e , j = 1,2 , … , p

3) Calculate the eigenvalues of R w and their corresponding eigenvectors. Establish the eigen equation λ I − R w = 0 and apply the Jacobi method to solve it to obtain the eigenvalues. Arrange the eigenvalues in descending order to make λ 1 ≥ λ 2 ≥ ⋯ ≥ λ p ≥ 0 . The corresponding eigenvectors can be obtained from the arranged eigenvalues as in Eq. 7, and then the eigenvector matrix V w can be constructed.

 ν 1 = ν 11 , ν 21 , ⋯ , ν p 1 T ν 2 = ν 12 , ν 22 , ⋯ , ν p 2 T ⋮ ν p = ν 1 p , ν 2 p , ⋯ , ν p p T

4) Calculate the PC contribution rate and cumulative contribution rate.

 δ j = λ j ∑ e = 1 p λ e j = 1,2 , … , p μ j = ∑ j = 1 j λ j ∑ j = 1 p λ j

Where δ j is the variance contribution rate of the j-th index; μ j is the cumulative variance contribution rate of the first j indexes; λ j is the eigenvalue of the j-th index. The variance contribution rate represents the ability to describe the original sample matrix. The larger the value is, the more information the original sample matrix contains, and the closer the reduced matrix is to the original matrix. The number of PCs is determined by cumulative contribution rate. To ensure the validity of data after dimension reduction and subsequent clustering, the first l PCs can reflect the information of p indexes when μ l ≥ 95 % , so the first l PCs play the main role (Ma et al., 2021).

5) Construct the new sample matrix after dimension reduction. After determining the number of PC, take the first l rows of V w to form a new matrix V w l , and the new sample matrix after dimension reduction can be calculated as Y w = X w * V w l .

Compared with the initial data of DFIG-based wind farm, the reliability and validity of the data after dimension reduction by PCA are improved. Taking them as the clustering indexes of DFIG-based wind farm can improve the accuracy of clustering results.

 2.2 Clustering of DFIG-based wind farm based on WFCM algorithm 2.2.1 WFCM clustering algorithm

FCM clustering algorithm is a method to describe and partition things with ambiguity or uncertainty (Zhao et al., 2022). When the traditional FCM algorithm is used to cluster wind turbines, the contribution degree of different clustering indexes to the clustering results is ignored, resulting in unreasonable clustering results. According to PCA, the contribution of each index is greatly different. Therefore, this paper proposes an FCM clustering algorithm with fixed weights of indexes by assigning corresponding weights to each index. The weights can be calculated by Eq. 10. P C = λ 1 λ P C 1 + λ 2 λ P C 2 + ⋯ + λ l λ P C l λ = λ 1 + λ 2 + ⋯ λ l ω k = λ k λ k = 1,2 , … , l Where P C l is the l-th PC; λ l is the eigenvalue of the l-th PC; ω k is the weight of the k-th PC; and the weight vector of l PCs is W = ω 1 , ω 2 , … , ω l .

The reduced dimensional PC is taken as the clustering index of the DFIG-based wind farm. Y w is the sample data set, where y i = y i 1 , y i 2 , … , y i l is the i-th sample with l PCs. The WFCM algorithm can be calculated as follows:

1) Initialize the number of clusters c c ∈ 2 , n , the maximum number of iterations T m , and the fuzzy index M . M is a constant that portrays the degree of fuzzification and takes a range of values from 1,2.5 . Generally, M is 2 (Wang et al., 2020).

2) Determine the affiliation matrix U and the clustering center matrix C . U is a matrix of n × c , where u i h is the affiliation of the i-th sample y i in Y w to the h-th clustering center, as shown in Eq. 11. C is a matrix of c × l , where C h k is given in Eq. 12.

 u i h = 1 ∑ m = 1 c ∑ k = 1 l ω k y i k − C h k 2 ∑ k = 1 l ω k y i k − C m k 2 M − 1 i = 1,2 , … , n , h = 1,2 , … , c C h k = ∑ i = 1 n u i h M y i k ∑ i = 1 n u i h M h = 1,2 , … , c , k = 1,2 , … , l Where y i k is the i-th sample in Y w ; d 2 y i k , C h k = y i k − C h k 2 is the Euclidean distance from the i-th sample y i k to the h-th clustering center.

3) Construct the objective function F U , C , W as in Eq. 13 and the constraints as in Eq. 14.

 min ⁡ F U , C , W = ∑ i = 1 n ∑ h = 1 c ∑ k = 1 l u i h M ω k d 2 y i k , C h k 0 ≤ u i h ≤ 1 ∑ h = 1 c u i h = 1 ∑ k = 1 l ω k = 1

4) Update the clustering center C h k and the affiliation degree u i h so that the objective function is the minimum, and keep repeating the above process until the number of iterations reaches T m .

2.2.2 Clustering validity test

Clustering validity research is the process of evaluating the quality of clusters and determining the optimal number of clusters by establishing validity indexes (Chen et al., 2018). In order to ensure the accuracy and objectivity of the clustering results, the Xie-Beni(XB) index, which is widely used in fuzzy clustering, is chosen to determine the optimal number of clusters (Geng et al., 2019). V X B = ∑ i = 1 n ∑ h = 1 c ∑ k = 1 l u i h M ω k d 2 C h k , y i k n min h ≠ a ∑ k = 1 l ω k d 2 C h k , C a k a = 1,2 , … , c

The numerator of V X B reflects the degree of compactness within a cluster, the smaller the value the more the compact; the denominator reflects the degree of separation between clusters, the larger the value the better the separation. Therefore, the smaller the value of V X B , the better the clustering effect. The c corresponding to the smallest V X B is the optimal clustering number.

To verify the rationality of the wind farm clustering indexes and further illustrate the effectiveness of WFCM algorithm and the rationality of clustering results, the Silhouette coefficient (SC) index and Calinski-Harabasz index (CHI) are selected for evaluation. The SC index is defined as follows: S C i = β i − α i max α i , β i i = 1,2 , … , n S C m = 1 n ∑ i = 1 n S C i Where α i is the intra-class compactness, which is the average distance between the sample y i and all other samples in the same cluster, and the smaller the value, the more compact the cluster; β i is the inter-class dispersion, which is the average distance between the sample y i in any cluster and all samples in the nearest cluster, and the larger the value, the more disperse the cluster. S C i is the profile coefficient of the i-th unit and S C i ∈ − 1,1 . S C i < 0 indicates that the clustering is not reasonable, and the larger the value of S C i , the more reasonable the clustering. S C m is the average profile coefficient, and the larger the value, the more effective the clustering result.

The CHI takes into account the inter-class dispersion B and intra-class compactness W. Larger I C H I indicates better inter-class dispersion and intra-class compactness. Where y ¯ is the mean value of all samples; ε h , i is the affiliation of the i-th sample to the h-th cluster. I C H I = B / c − 1 W / n − c B = ∑ h = 1 c ∑ i = 1 n ε h , i C h − y ¯ 2 W = ∑ h = 1 c ∑ i = 1 n ε h , i y i − C h 2 ε h , i = 1 , y i ∈ C h 0 , y i ∈ C h

The flow chart of the proposed WFCM algorithm is demonstrated in Figure 1.

Flow chart of the WFCM algorithm.

 2.3 Calculation of equivalent parameters for DFIG-based wind farm

According to the principle of constant output characteristics before and after equivalence, all the units in the cluster can be equivalent using the model of one unit. Assuming that there are d DFIGs of the same type in the same cluster after clustering by WFCM algorithm, the equivalent parameters can be calculated as follows: S e q = ∑ i = 1 d S i , P e q = ∑ i = 1 d P i , Q e q = ∑ i = 1 d Q i H e q = ∑ i = 1 d H i , D e q = ∑ i = 1 d D i , K e q = ∑ i = 1 d K i R e q = R s d , X e q = X s d , v e q = 1 d ∑ i = 1 d v i 3 1 / 3 , Z T e q = Z T d where the subscript eq is the equivalent parameter; the subscript i is the i-th wind turbine in the same cluster; S, P, Q, H, D, K, R s , X s , v, and Z T are the rated capacity, active power, reactive power, inertia time constant, damping coefficient of shafting, stiffness coefficient of shafting, stator resistance, stator reactance, effective input wind speed, and transformer impedance, respectively.

According to the constant voltage difference before and after equivalent, the equivalent impedance of the trunk type collector line of the wind farm Z e q 1 can be calculated by Eq. 21, and the equivalent impedance of the radial collector line Z e q 2 can be calculated by Eq. 22. Where Z f is the impedance of the f-th line; Z i is the impedance of the i-th turbine line. Z e q 1 = ∑ i = 1 d ∑ f = 1 i Z f ∑ b = f n P b P i ∑ i = 1 d P i 2 Z e q 2 = ∑ i = 1 d P i 2 Z i ∑ i = 1 d P i 2

3 Multi-parameter coupling optimization design of DFIG-based wind farm based on OEM-RSM

In (Zhang et al., 2022), we have optimized the proportional and integral coefficients of GSC, RSC and PLL for a single DFIG connected to weak grid by series compensation system. Considering the SSO characteristics and dynamic characteristics of the system, the integrated average value of DC-link voltage overshoot Δ U d c , active power overshoot Δ P s , and grid current angular frequency overshoot Δ ω s when the disturbance occurs have been chosen as the performance indexes. After weighting by the analytic hierarchy process-coefficient of variation compound weighting method, the minimization of the weighted sum OD has been taken as the optimization objective. The control parameters have been optimized by OEM and particle swarm algorithm (PSO). As a result, the sensitivity and influence trend of control parameters on the system performance have been revealed, and the SSO has been effectively mitigated while ensuring strong anti-interference. However, the strategy is unable to analyze the interaction between parameters and is prone to fall into local optimum due to the use of intelligent algorithm. Therefore, this paper improves it and proposes a multi-parameter coupling optimization strategy for DFIG-based wind farms combining OEM and RSM.

The single-factor experiment method is to conduct experiments on only one factor while keeping all other factors constant, so as to explore the optimal range of values for each factor. The RSM uses a multiple quadratic regression equation to establish the functional relationship between the influencing factors and the response values based on the data obtained by the experimental design method, and the optimal parameter combination can be found by analyzing the equation (Si et al., 2018). It is generally utilized for experimental designs with few factors (generally 2–4) and few factor levels. In this paper, control parameters of each equivalent wind turbine are changed simultaneously and preliminary optimized by OEM to obtain the three factors that have the most significant influence on SSO and dynamic characteristics of the grid-connected system and seven insignificant factors. The seven insignificant factors can be valued according to the optimization results of OEM. Based on the values of three significant factors and the influence trend of them on OD, the single-factor experiment method is utilized to further determine the appropriate range of values. Set the values of each factor in RSM, a 3-factor, 3-level response surface experiment can be designed by Box-Behnken central combination. The experimental data of OD can be fitted according to the model with the best fit, and the model can be analyzed by analysis of variance (ANOVA) and constructed statistics. When all the indexes of the model meet the criteria, the regression equation between each factor and the evaluation index can be fitted, which can be expressed by Eq. 23. Then, the three-dimensional response surface can be plotted to obtain the interaction between the factors and the optimal OD of the system and the corresponding levels of each factor. y = c 0 + ∑ i = 1 N c i ξ i + ∑ i = 1 N c i i ξ i 2 + ∑ i < j ∑ c i j ξ i ξ j + ε , i , j = 1,2 , ⋅ ⋅ ⋅ , N

Where ε is the fitting error, c 0 is a constant; c i , c i i , and c i j are the coefficients of the first-order, second-order, and interaction terms, respectively; ξ i and ξ j are experiment factors.

 4 Case study 4.1 Parameter setting

A time-domain simulation model of DFIG-based wind farm connected to weak grid by series compensation is built in Matlab/Simulink as shown in Figure 2. The wind farm contains 25 DFIGs of the same type, and the wind farm parameters and DFIG control parameters are shown in Table 1 and Table 2. The distance between each row of DFIGs is 500 m. The DFIG-based wind farm is connected to the 35 kV medium voltage bus through terminal transformers and overhead lines, and then the main transformer is connected to the 220 kV grid through the transmission line. According to the wind speed, wind direction, and wake effect, the wind speed of each unit in the wind farm can be calculated as shown in Table 3.

Topology structure of DFIG-based wind farm.

Parameters of DFIG-based wind farm.

Parameters Value
Base capacity S B 0.8 MVA
Rated voltage U B 690 V
Rated frequency f B 50 Hz
Stator resistance R s 0.03 p.u
Rotor resistance R r 0.016 p.u
Stator leakage inductance L l s 0.18 p.u
Rotor leakage inductance L l r 0.16 p.u
Stator-rotor mutual inductance L m 2.9 p.u
Grid-side filter inductance L l 0.3 p.u
Grid-side filter inductance parasitic resistance R l 0.02 p.u
DC-link capacitor C d c 0.01 F
DC-link voltage U d c _ r e f 1150 V
Terminal transformer rated capacity 1.67 MVA
Terminal transformer voltage ratio 0.69/35 kV
Terminal transformer X T 6%
Main transformer rated capacity 50 MVA
Main transformer voltage ratio 35/220 kV
Main transformer X T 10.5%
Line resistance R L 0.115 Ω/km
Line inductance L L 1.05e-3H/km
Line length 30 km
Distance between unit rows 0.5 km

DFIG control parameters.

Parameters k p _ d c k i _ d c k p _ i l k i _ i l k p _ p
Value 6 400 0.8 5 1.6
Parameters k i _ p k p _ i r k i _ i r k p _ p l l k i _ p l l
Value 100 0.6 8 5 1.1

Wind speed of each unit in DFIG-based wind farm (m/s).

W T i j 1 2 3 4 5
1 13.0000 13.0000 13.0000 13.0000 13.0000
2 13.0000 11.1719 11.1719 11.1719 13.0000
3 11.1719 11.1719 9.3866 9.3866 13.0000
4 9.3866 11.1719 9.3866 7.9655 13.0000
5 7.9655 11.1719 9.3866 7.9655 7.1540
4.2 Equivalent modeling results for DFIG-based wind farm

The initial SCR is set to 3 and SCD is 20%. The SCR is changed to 2 and SCD is changed to 60% at 2.1 s. The SCR and SCD are changed back to initial values at 3.1 s. A three-phase ground short-circuit fault is occurred at 4 s and removed at 4.15 s. The data collected at 2s are ω s 1 , U d c 1 , i s d 1 , i s q 1 , i r d 1 , i r q 1 , the data collected at 3 s are ω s 2 , U d c 2 , i s d 2 , i s q 2 , i r d 2 , i r q 2 , the data collected at 4.1s are ω s 3 , U d c 3 , i s d 3 , i s q 3 , i r d 3 , i r q 3 , and the data collected at 4.6s are ω s 4 , U d c 4 , i s d 4 , i s q 4 , i r d 4 , i r q 4 . The data set consisting of 25 DFIGs and 24 initial clustering indexes is the sample matrix X w of 25 × 24 .

4.2.1 Results of dimension reduction by PCA

The dimension of X w is reduced by PCA. The eigenvalues of each component are arranged in descending order, and the corresponding variance contribution rate and cumulative contribution rate are calculated as shown in Table 4.

Contribution rates of PCs.

PC Eigenvalue Variance contribution rate Cumulative contribution rate PC Eigenvalue Variance contribution rate Cumulative contribution rate
1 12.34236 0.514265 0.514265 13 2.79E-06 1.16E-07 1
2 8.155651 0.339819 0.854084 14 1.11E-06 4.63E-08 1
3 2.717655 0.113236 0.967319 15 2.33E-07 9.71E-09 1
4 0.567313 0.023638 0.990957 16 6.1E-08 2.54E-09 1
5 0.157528 0.006564 0.997521 17 8.24E-09 3.43E-10 1
6 0.040513 0.001688 0.999209 18 6.79E-10 2.83E-11 1
7 0.015618 0.000651 0.99986 19 3.16E-11 1.32E-12 1
8 0.002584 0.000108 0.999967 20 8.03E-12 3.34E-13 1
9 0.000376 1.57E-05 0.999983 21 1.56E-13 6.52E-15 1
10 0.000266 1.11E-05 0.999994 22 1.57E-14 6.51E-16 1
11 9.14E-05 3.81E-06 0.999998 23 4.85E-15 2.05E-16 1
12 4.52E-05 1.88E-06 1 24 -1.6E-16 -6.8E-18 1

It can be found that the 24 initial clustering indexes contain strongly correlated indexes, which will cause bias in the clustering results, so the dimension of data need to be reduced. The cumulative contribution rate of the first three PCs has reached 96.7319%, greater than 95%, which meets the requirement for the selection of PCs in PCA, while the remaining indexes have little influence on contribution. The contribution rates of the first three PCs are 51.4265%, 33.9819%, and 11.3236%, respectively. Therefore, the first three PCs can be selected for dimension reduction, i.e., the initial clustering indexes are reduced to 3 dimensions. The component with the eigenvalue of 12.34236 is the first PC, the component with the eigenvalue of 8.155651 is the second PC, and the component with the eigenvalue of 2.717655 is the third PC, which basically retains the information of the original data. The contribution rates of these three PCs are given in Figure 3.

Contribution rates of PCs.

So far, the dimension of 24 initial clustering indexes are reduced to obtain 3 PC indexes by PCA, and the new sample matrix after dimension reduction Y w is of 25 × 3 . The weight vector of the 3 PCs can be calculated as W = 0.5316,0.3513,0.1171 according to (10).

 4.2.2 Clustering results for DFIG-based wind farm

Y w is clustered by WFCM clustering algorithm. c is set as 2, T m is set as 100, and M is set as 2. According to the iteration of (15), the optimal cluster for 25 DFIGs c o p = 4 . Thus the clustering results with 3 PCs as multiple clustering indexes can be obtained in Table 5; Figure 4.

The clustering results for DFIG-based wind farm by WFCM.

Cluster Clustering results
1 W T 14 , W T 24 , W T 15 , W T 25 , W T 35 , W T 45
2 W T 44 , W T 54 , W T 55
3 W T 41 , W T 51 , W T 33 , W T 43 , W T 53 , W T 34
4 W T 11 , W T 21 , W T 31 , W T 12 , W T 22 , W T 32 , W T 42 , W T 52 , W T 13 , W T 23

The clustering results for DFIG-based wind farm by WFCM.

In order to verify the effectiveness of WFCM algorithm and the reasonableness of clustering results, Figure 5A shows the profile coefficients of 25 DFIGs when the clusters is 4, Figure 5B shows the results of average profile coefficient, and Figure 5C shows the results of CHI.

Clustering evaluation indexes. (A) Profile coefficient; (B) Average profile coefficient; (C) CHI.

It can be observed from Figure 5A that the profile coefficients of 25 DFIGs are all greater than 0 when the cluster is 4, indicating the clustering is reasonable and effective. From Figure 5B,C, the average profile coefficient and CHI are both maximum at the cluster of 4, which is consistent with XB index. Therefore, the 25 DFIGs built in this paper can be divided into four equivalent clusters.

 4.3 Multi-parameter coupling optimization results for DFIG-based wind farm

The equivalent wind farm after clustering are preliminarily optimized by OEM, and the values of each control parameter can be obtained as shown in Table 6. From the results of OEM, the three factors that have the most significant effect on SSO characteristics and dynamic characteristics of DFIG-based wind farm connected to weak grid by series compensation system are GSC DC voltage loop proportional coefficient k p _ d c , RSC current loop proportional coefficient k p _ i r , and PLL proportional coefficient k p _ p l l . And the optimized values are 10, 0.5, and 3, respectively. The single-factor experiments are conducted for these three factors to narrow their range of values. The corresponding generalized Nyquist curves of the system scaled up at (−1, j0) are drawn in Figure 6.

Optimal parameter combination based on OEM.

Parameters k p _ d c k i _ d c k p _ i l k i _ i l k p _ p
Value 10 380 0.6 5 1.6
Parameters k i _ p k p _ i r k i _ i r k p _ p l l k i _ p l l
Value 100 0.5 8 3 0.8

Generalized Nyquist curves of the system with (A) different k p _ d c , (B) different k p _ i r , and (C) different k p _ p l l .

According to the generalized Nyquist stability criterion, if the Nyquist curves do not enclose (−1, j0), the system is stable and the further away from that point, the more stable the system is. It can be found that the optimal value exists between 45 and 55 for k p _ d c , 0.45 to 0.55 for k p _ i r , and 2 to 3 for k p _ p l l . And the Nyquist curves do not enclose (−1, j0) within the above values, indicating that the system is stable. So the value ranges of each parameter determined by the single-factor experiment method is reasonable and valid. The response surface experiment is designed in Design-Expert 11, and the level values of k p _ d c are set as 45, 50, and 55, the level values of k p _ i r are 0.45, 0.5, and 0.55, and the level values of k p _ p l l are 2, 2.5, and 3. A total of 17 sets of experiments are conducted. A quadratic polynomial model is fitted to the experimental data of OD, and the results of ANOVA on the model are shown in Table 7. The constructed statistics are shown in Table 8.

ANOVA results for OD of RSM.

Source Sum of squares df Mean square F-value P-value Significance
Model 1415.17 9 157.24 217.94 <0.0001 ***
A- k p _ d c 675.44 1 675.44 936.17 <0.0001 ***
B- k p _ i r 64.81 1 64.81 89.83 <0.0001 ***
C- k p _ p l l 39.22 1 39.22 54.36 0.0002 **
AB 0.9861 1 0.9861 1.37 0.2806 —
AC 0.0065 1 0.0065 0.009 0.9272 —
BC 6.31 1 6.31 8.75 0.0212 *
A2 0.0664 1 0.0664 0.092 0.7704 —
B2 622.97 1 622.97 863.45 <0.0001 ***
C2 0.4871 1 0.4871 0.6751 0.4384 —
Residual 5.05 7 0.7215 —
Lack of fit 2.87 3 0.9555 1.75 0.2951 —
Pure error 2.18 4 0.546 — — —

P-value is the probability of significance. The smaller the p-value, the more significant the effect of the factor on the response value. p < 0.0001 indicates that the effect of the factor on the response value is extremely significant, denoted by ***; p ≤ 0.01 indicates that the effect is highly significant, denoted by **; 0.01 < p ≤ 0.05 indicates that the effect is significant, denoted by *; p > 0.05 indicates that the effect is not significant.

Statistics for OD of RSM.

C.V.% R 2 Adjusted R 2 Predicted R 2 Adeq precision
2.04 0.9964 0.9919 0.9653 54.9497

As can be seen from Table 7, the p-value of the model is less than 0.0001, indicating that the overall effect of the parameter combination on OD is extremely significant. Linear terms A, B and quadratic term B2 have an extremely significant effect on OD, quadratic term C2 has a highly significant effect on OD, interaction term BC has a significant effect on OD, and the effects of other terms are not significant. The significance of the effect of each factor is A > B > C, that is, k p _ d c > k p _ i r > k p _ p l l . The p-value of the lack of fit in the table is greater than 0.05 (0.2951), indicating that the lack of fit is not significant, so the model is able to fit the data accurately.

The statistics of the experiment data should meet the requirements of C.V.% less than 10, precision greater than 4, R 2, adjusted R 2, and predicted R 2 are close to 0.9, and the difference between adjusted R 2 and predicted R 2 is less than 0.2. From Table 8, it can be found that the statistics of the model all meet the requirements, indicating that the fit and correlation of the model are good and the error is small, so it can also predict OD under different combinations of operating parameters.

A quadratic multiple regression is fitted to a total of 51 data points from 17 sets of experiments to create a regression equation as in (24). The response surface of each factor combination on OD is depicted in Figure 7. O D = 35.72 − 9.19 A + 2.85 B − 2.21 C − 0.4965 A B − 0.0402 A C + 1.26 B C − 0.1256 A 2 + 12.16 B 2 + 0.3401 C 2

Effect of factors on OD. (A) Effect of k p _ d c and k p _ i r on OD; (B) Effect of k p _ d c and k p _ p l l on OD; (C) Effect of k p _ i r and k p _ p l l on OD.

As shown in Figure 7A, OD decreases and then rises as k p _ i r increases from 0.45 to 0.55, and there is a minimal value between 0.47 and 0.51; when k p _ i r is small, OD decreases gradually with the increase of k p _ d c from 45 to 55, and the decrease in OD increases as k p _ i r increases. From Figure 7B, OD tends to decrease as k p _ p l l increases from 2 to 3 and is flatter; OD decreases rapidly as k p _ d c increases. From Figure 7C, OD gradually decreases as k p _ p l l increases, and the effect of k p _ i r on OD is much more significant than that of k p _ p l l . As shown in Figure 7, the interaction between k p _ i r and k p _ p l l is significant, which is consistent with the ANOVA results.

The optimal solution that satisfies the optimization index can be obtained by optimizing OD through (24). The optimal predicted and experimental value can be calculated by substituting the corresponding factor values into the regression equation and simulation model and then the error between them can be obtained. The results are listed in Table 9.

Optimal parameter combination based on RSM.

k p _ d c k p _ i r k p _ p l l Predicted Experimental Error (%)
54.961 0.495 2.979 24.411 24.986 2.3

Noticeably, the OD predicted by the model is 24.411, which is 1.33% lower than the minimum value of the response surface experiment (24.741). The error between the predicted and the experimental value is 2.3% (not more than 5%), indicating that the regression model is accurate and effective for the optimization.

 4.4 Simulation results

To verify the effectiveness of the multi-parameter coupling optimization design strategy in this paper, the strategy combining the OEM and PSO in (Geng et al., 2019) is used for comparative verification. The SCR is changed from 3 to 2 and SCD is changed from 20% to 60% at 2s. The FFT analysis of the grid-connected current and generalized Nyquist curves of the system after 2s are depicted in Figure 8. Figure 9 shows the d-axis stator current waveform.

FFT analysis of the grid-connected current and generalized Nyquist curves of system. (A) FFT analysis of the grid-connected current; (B) Generalized Nyquist curves.

D-axis stator current waveform.

The FFT analysis is used to determine the frequency at which oscillations occur in the grid-connected system. It can be found from Figure 8 that the system appears SSO with a frequency of 39 Hz (lower than power frequency 50 Hz). The generalized Nyquist curve after 2 s changes from enclosing (−1, j0) to not enclosing and gradually moving away when control parameters change from the original parameters to the optimized parameters, indicating that the system changes from unstable to stable. It can be observed from Figure 9 that the system is destabilized after the grid condition changes to very weak grid and high series compensation degree with the original parameters, resulting in divergent SSR. After optimization by the strategy proposed in this paper and the strategy in (Geng et al., 2019), the system oscillations converge. The d-axis stator current returns to stability after about 2 s with the optimized parameters of the OEM-PSO, whereas it only oscillates for about 0.7 s with the optimized parameters of the OEM-RSM. Apparently, both optimization methods are effective in improving the stability of the sub-synchronous state of the system, and the OEM-RSM is more effective. The simulation results are consistent with the theoretical analysis results.

Figure 10 illustrates the DC-link voltage, active power, and grid current angular frequency curves when the SCR is 3 and SCD is 20%, with a three-phase ground short-circuit fault occurring at 2 s and removed after 0.1 s. To facilitate the analysis, Δ U d c , Δ P s , and Δ ω s with 3 parameter combinations are calculated in Table 10.

DC-link voltage, active power, and grid current angular frequency waveform. (A) DC-link voltage; (B) Active power; (C) Grid current angular frequency.

Results of evaluation indexes with 3 parameter combinations.

Evaluation indexes Δ U d c /V Δ P s /MW Δ ω s /pu
Original 67.9229 1.1119 6.3085E-4
OEM-PSO 56.1519 1.0037 4.5033E-4
OEM-RSM 49.8196 0.9692 3.3627E-4

From Figure 10, after optimization of the proposed strategy, the overshoot and settling time of DC-link voltage, active power, and grid current angular frequency with disturbance are much smaller, and the fluctuation is significantly weakened and the anti-interference is strongly enhanced. It can be found from Table 10 that Δ U d c , Δ P s , and Δ ω s with the optimized parameters of the OEM-PSO and the OEM-RSM are both reduced. Visibly, Δ U d c is reduced by 26.65%, Δ P s is reduced by 12.83%, and Δ ω s is reduced by 46.70% compared to the original parameters after optimization by the OEM-RSM. In conclusion, both optimization methods can improve the dynamic characteristics of the system. Moreover, the OEM-RSM is more effective in optimizing both the SSO characteristics and dynamic stability.

 5 Conclusion

Based on the model of DFIG-based wind farm connected to weak grid by series compensation system, this paper proposes a clustering method based on the PCA and WFCM algorithm, with the dominant variables of the system SSO mode as the initial clustering indexes. For the equivalent wind farm after clustering, the DFIG control parameters are taken as the optimization object, and the sub-synchronous state stability and dynamic characteristics of the system are simultaneously optimized as the optimization objective. The optimal design strategy combining the OEM and RSM is adopted to quantitatively evaluate the interaction between controller parameters and the significance of their effects on the SSO and dynamic characteristics of the grid-connected system. The optimal parameter combination is obtained to achieve the effective mitigation of SSO. The highlights are as follows:

1) The proposed equivalence modeling method for DFIG-based wind farm takes into account the correlation between clustering indexes and SSO mode as well as the contribution degree of different clustering indexes to clustering results. It eliminates the redundancy of data and can be applied to study SSO in large-scale wind farm grid-connected systems.

2) Compared with the OEM-PSO, the proposed multi-parameter coupling optimization strategy can effectively improve the anti-interference and system stability in the sub-synchronous state, which can be applied to study the optimization problems of complex systems with multiple parameters and indexes.

3) k p _ d c , k p _ i r , and k p _ p l l have significant effects on SSO characteristics and dynamic characteristics of the system, and there is a strong interaction between k p _ i r and k p _ p l l .

In addition, the current work will also be carried out for better analyzing and solving the SSO problem in actual situations in the near future. Since the control parameters of each cluster are changed synchronously in this paper, the interaction between clusters and the influence of control parameters of different clusters at different locations on SSO characteristics of the system will be fully considered to improve the proposed strategy.

 Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

 Author contributions

QZ put forward the initial concept and gave technical guidance in the whole process. XJ carried out the overall work for the paper and drafted the whole paper. FZ, HY, and BZ contributed to the revision and typesetting of the manuscript. All authors read and approved the final manuscript.

 Funding

This work was supported in part by the State Key Laboratory of HVDC, Electric Power Research Institute, CSG (No. SKLHVDC-2020-KF-14).

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 References Bakir H. Merabet A. Dhar R. K. Kulakistet A. A. (2020). Bacteria foraging optimisation algorithm based optimal control for doubly‐fed induction generator wind energy system. IET Renew. Power Gener. 14 (11), 1850–1859. 10.1049/iet-rpg.2020.0172 Chen A. Xie D. Zhang D. Gu C. Wang K. (2019). PI parameter tuning of converters for sub-synchronous interactions existing in grid-connected DFIG wind turbines. IEEE Trans. Power Electron. 34 (7), 6345–6355. 10.1109/tpel.2018.2875350 Chen Y. Wu H. Shi J. (2018). Application of singular value decomposition algorithm to dimension-reduced clustering analysis of daily load profiles. Automation Electr. Power Syst. 42 (3), 105. 10.7500/AEPS20170309008 Dong X. Hou J. Tian X. Liu Y. (2016). Comparison of the impact level on subsynchronous resonance characteristics of dominant factors by multiple wind farms at different locations. High. Volt. Eng. 42 (07), 2259. 10.13336/j.1003-6520.hve.20160713029 Geng J. Qian X. Zhou S. (2019). New fuzzy clustering validity index. Appl. Res. Comput. 36 (4), 1001. 10.19734/j.issn.1001-3695.2017.10.0991 Han J. Miao S. Li L. Yang W. Li Y. (2020). Wind turbines clustering in wind farm based on multi-view transfer learning and synthetic optimization of parameters in equivalent wind farm. Proc. CSEE 40 (15), 4866. 10.13334/j.0258-8013.pcsee.190952 Han K. L. Huang T. W. Yin L. F. (2021). Quantum parallel multi-layer Monte Carlo optimization algorithm for controller parameters optimization of doubly-fed induction generator-based wind turbines. Appl. Soft Comput. 112, 107813. 10.1016/j.asoc.2021.107813 Hu B. Nian H. Li M. Xu Y. Liao Y. Yang J. (2022). Impedance-based analysis and stability improvement of DFIG system within PLL bandwidth. IEEE Trans. Ind. Electron. 69 (6), 5803–5814. 10.1109/tie.2021.3088369 Huang J. Wu X. Zheng Z. Huang Y. Li W. (2022). Multi-objective optimal operation of combined cascade reservoir and hydrogen system. IEEE Trans. Ind. Appl. 58 (2), 2836–2847. 10.1109/tia.2021.3138949 Lei H. Zhen C. Yue X. Xu G. Bao J. (2013). Mitigating thermal power’s SSR by additional damping controller of DFIG. Electr. Power 46 (4), 63. 10.3969/j.issn.1004-9649.2013.04.014 Leon A. E. Amodeo S. J. Mauricio J. M. (2021). Enhanced compensation filter to mitigate subsynchronous oscillations in series-compensated DFIG-based wind farms. IEEE Trans. Power Deliv. 36 (6), 3805–3814. 10.1109/tpwrd.2021.3049318 Li C. Cai W. Zhao R. Yu C. Zhang Q. (2019). Power consumption behavior analysis of AP cluster users based on optimized SAX and weighted load characteristic index. Trans. China Electrotech. Soc. 34 (S1), 368. 10.19595/j.cnki.1000-6753.tces.l80097 Liu J. Yao W. Wen J. Fang J. Jiang L. He H. (2020). Impact of power grid strength and PLL parameters on stability of grid-connected DFIG wind farm. IEEE Trans. Sustain. Energy 11 (1), 545–557. 10.1109/tste.2019.2897596 Liu J. Yao W. Wen J. (2017). Small signal stability analysis and control of double-fed induction generator considering influence of PLL and power grid strength. Proc. CSEE 37 (11), 3162. 10.13334/j.0258-8013.pcsee.160857 Liu M. Pan W. Zhang Y. Zhao K. Zhang S. Liu T. (2019). A dynamic equivalent model for DFIG-based wind farms. IEEE Access 7, 74931–74940. 10.1109/access.2019.2918359 Ma X. Yang Q. Chen C. Luo H. Zheng D. (2021). Harmonic and interharmonic analysis of mixed dense frequency signals. IEEE Trans. Ind. Electron. 68 (10), 10142–10153. 10.1109/tie.2020.3026288 Miao Z. X. (2012). Impedance-model-based SSR analysis for type 3 wind generator and series-compensated network. IEEE Trans. Energy Convers. 27 (4), 984–991. 10.1109/tec.2012.2211019 Pan X. Zhang C. Ju P. (2015). Coherency-based dynamic equivalence of wind farm composed of doubly fed induction generators. Power Syst. Technol. 39 (03), 621. 10.13335/j.1000-3673.pst.2015.03.006 Shabanikia N. Nia A. A. Tabesh A. Khajehoddin S. A. (2021). Weighted dynamic aggregation modeling of induction machine-based wind farms. IEEE Trans. Sustain. Energy 12 (3), 1604–1614. 10.1109/tste.2021.3057854 Si J. Zhao S. Feng H. Cao R. Hu Y. (2018). Multi-objective optimization of surface-mounted and interior permanent magnet synchronous motor based on taguchi method and response surface method. Chin. J. Electr. Eng. 4 (1), 67. 10.23919/CJEE.2018.8327373 Sun Y. Mao Y. Li Z. Zhang X. Li F. (2021). A comprehensive clustering method of user load characteristics and adjustable potential based on power big data. Proc. CSEE 41 (18), 6259. 10.13334/j.0258-8013.pcsee.201928 Wang Z. Yang C. Oh S. -K. Fu Z. Pedrycz W. (2020). Robust multi-linear fuzzy SVR designed with the aid of fuzzy C-means clustering based on insensitive data information. IEEE Access 8, 184997–185011. 10.1109/access.2020.3030083 Wanigasekara C. Swain A. Almakhles D. Zhou L. (2022). Design of delta–sigma-based PID controller for networked wind energy conversion systems. IEEE Trans. Ind. Appl. 58 (1), 879–889. 10.1109/tia.2021.3126574 Xia Y. Chen Y. Song Y. Strunz K. (2020). Multi-scale modeling and simulation of DFIG-based wind energy conversion system. IEEE Trans. Energy Convers. 35 (1), 560–572. 10.1109/tec.2019.2953893 Zhang F. Yuan H. Zhou B. Zhang Q. Jin X. (2022). “Parameter optimization and SSR mitigation of DFIG considering PLL,” in 2022 Power System and Green Energy Conference (PSGEC), 430 Zhao K. Dai Y. Jia Z. Ji Y. (2022). General fuzzy C-means clustering strategy: Using objective function to control fuzziness of clustering results. IEEE Trans. Fuzzy Syst. 30 (9), 3601–3616. 10.1109/tfuzz.2021.3119240 Zhou Y. Zhao L. Matsuo I. Lee W. (2019). A dynamic weighted aggregation equivalent modeling approach for the DFIG wind farm considering the weibull distribution for fault analysis. IEEE Trans. Ind. Appl. 55 (6), 5514–5523. 10.1109/tia.2019.2929486