COMPLEXITYComplexity1099-05261076-2787Hindawi10.1155/2020/94382859438285Research ArticleThe Small-Signal Stability of Offshore Wind Power Transmission Inspired by Particle Swarm OptimizationLiJiening1HuangHanqi1ChenXiaoning2https://orcid.org/0000-0002-7376-2925PengLingxi13https://orcid.org/0000-0002-1017-0500WangLiang4https://orcid.org/0000-0002-1998-8253LuoPing5HeShuping1School of Mechanical and Electrical EngineeringGuangzhou UniversityGuangzhou 510006Chinagzhu.edu.cn2School of Mathematics and Information ScienceGuangzhou UniversityGuangzhou 510006Chinagzhu.edu.cn3Data Recovery Key Laboratory of Sichuan ProvinceNeijiang Normal UniversitySichuan 641100Chinanjtc.edu.cn4School of Public AdministrationGuangzhou UniversityGuangzhou 510006Chinagzhu.edu.cn5School of Economics and StatisticsGuangzhou UniversityGuangzhou 510006Chinagzhu.edu.cn202015720202020060520202205202015720202020Copyright © 2020 Jiening Li et al.This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Voltage source converter-high-voltage direct current (VSC-HVDC) is the mainstream technology of the offshore wind power transmission, which has been rapidly developed in recent years. The small-signal stability problem is closely related to offshore wind power grid-connected safety, but the present study is relatively small. This paper established a mathematical model of the doubly fed induction generator (DFIG) integrated into the IEEE9 system via VSC-HVDC in detail, and small-signal stability analysis of offshore wind farm (OWF) grid connection is specially studied under different positions and capacities. By selecting two load nodes and two generator nodes in the system for experiments, the optimal location and capacity of offshore wind power connection are obtained by comparing the four schemes. In order to improve the weak damping of the power system, this paper presents a method to determine the parameters of the power system stabilizer (PSS) based on the particle swarm optimization (PSO) algorithm combined with different inertia weight functions. The optimal position of the controller connected to the grid is obtained from the analysis of modal control theory. The results show that, after joining the PSS control, the system damping ratio significantly increases. Finally, the proposed measures are verified by MATLAB/Simulink simulation. The results show that the system oscillation can be significantly reduced by adding PSS, and the small-signal stability of offshore wind power grid connection can be improved.

National Natural Science Foundation of China6177214761100150University Innovation Team Construction Project of Guangdong Province2015KCXTD014Guangdong Province Philosophy and Social Science FoundationGD19CSH03the key Project of Science and Technology Plan of Guangdong Province2020b1010010014Open Research Fund Program of Data Recovery Key Laboratory of Sichuan Province
1. Introduction

Offshore wind power, as a clean and sustainable technology, has been developed rapidly in recent years . At present, the total installed capacity in Europe is increasing every year, among which the UK and Germany dominate the offshore wind power industry . According to Wind Europe’s High Scenario, it is estimated that the offshore wind energy capacity in Europe will reach 99 GW by 2030 . Nowadays, the mainstream transmission technology is VSC-HVDC . Compared to other traditional transmission modes, it has the advantages of independent control of the active and reactive power output, smaller power loss, and lower voltage drop .

Small-signal stability analysis is to study the dynamic response characteristics of the power system after small disturbances (including random fluctuations in power generation or consumption) and to evaluate their ability to suppress oscillations . It is of great significance to promote the development of wind power and improve the safe stability of the system, which needs urgent attention . However, there are few studies on the small-signal stability analysis of the offshore wind power grid connected by VSC-HVDC transmission at present. Based on the small-signal stability problem of the offshore wind power system transmitted by VSC-HVDC, the damping control of the system is increased by introducing PID regulation at the converter station. The modeling method does not consider the aerodynamic model, and the selection of PID parameters based on eigenvalues has limitations . Small-signal stability analysis method from the perspective of voltage and frequency is discussed in [14, 15]. In , the variability of wind power was considered as the interference source of the system and used Prony analysis and swing-based frequency response metric to study the influence of the small signal on the large offshore wind power system. However, the whole system modeling process has not been described in detail but is focused on analysis. In , the model of the offshore AC network for point-to-point VSC-HVDC transmission is deduced in detail, and the droop gain boundary under the stable operation is determined by eigenvalue analysis. A low-pass filter was introduced to improve the system damping, but the effect of different grid-connected positions and capacity on the small signal of the system was not considered. At present, neural networks [16, 17], deep learning, and machine learning methods have been widely used in various fields, and remarkable results have been achieved [18, 19].

To improve the stability of the system under the small signal and large signal, the PSO algorithm is used to optimize the parameters of the DFIG controller. However, the PSO algorithm will fall into local optimum, and more controller equations increase the dimension of the state matrix, which will slow down the fast operation of the system . In these studies, although scholars have studied the small-signal stability of offshore wind power grid connection in different fields, it is still in the exploratory stage to change the grid-connected position and capacity of offshore wind farms for small-signal stability analysis and to apply damping control.

Aiming at the problem of small-signal stability of the offshore wind power grid-connected system, this paper firstly establishes a complete VSC-HVDC transmission offshore wind power grid-connected mathematical model and then analyzes the impact of OWFs on the small-signal stability of the system under different positions and capacities. Secondly, PSS is introduced for damping control. The optimal position of the PSS connected to the power grid is determined by modal control theory, and the parameters of the controller are determined by the PSO algorithm combined with different inertia weight functions. Finally, the correctness of the established model and the effectiveness of the proposed control measures are verified by MATLAB/Simulink simulation.

2. Modeling of Offshore Wind Power Transmitted by VSC-HVDC

This paper uses the typical offshore wind power grid connection system  for reference, and the designed system topology is shown in Figure 1. The OWFs are composed of 10 2-MW DFIGs, and they are equivalent to a DFIG representation according to the aggregated model . The wind turbine (WT) runs on the low-speed (LS) shaft, the generator runs on the high-speed (HS) shaft, and the gearbox connects the two to act as the booster. Each 690 V DFIG is connected to 10 kV bus through the boost transformer XT0 and then sent to 110 kV bus through the boost transformer XT1 and the transmission line XT1. VSC1 and VSC2 are converter stations, which play the roles of the rectifier and the inverter, respectively. The HVDC plays the role of power transmission. In addition, there are resistance-capacitance (RC) filters and phase reactors, whose influence on the system is ignored in this paper. The VSC-HVDC outlet us2 is connected to the 230 kV 3-machine, 9-bus test system  via the boost transformer XT2 and the transmission line XL2, which is represented by IEEE9.

The topology of the whole system.

2.1. Aerodynamic Model

WT is driven by the wind, which is converted into mechanical energy through three blades. The equations of the aerodynamic model  are given by(1)CP=0.22116λi0.4β5e12.5/λi,1λi=1λ+0.08β20.035β3+1,Tt=Ptωt=ρπr2CPv32ωt,where Cp is the coefficient of wind energy utilization; β is the pitch angle; λ is the tip speed ratio (λ = ωtr/ν); ωt is the mechanical angular velocity of the WT; r is the radius of the wind wheel; ν is the wind speed; Pt is the mechanical power output by the WT; Tt is the mechanical torque; and ρ is the air density.

2.2. Shafting Model

In order to ensure the accuracy of calculation, this paper selects the two-mass block shafting model. The mathematical model of the shafting model  is as follows:(2)2Htdωtdt=TtKθDωt1sωs,dθdt=ωbωt1sωs,2Hgωsdsdt=Kθ+Dωt1sωsTe,where Ht and Hg are the inherent inertial time constants of the WT and the generator; Tt and Te are the mechanical torque of the WT and the electromagnetic torque of the generator; K and D are the stiffness and damping coefficients of the shaft; θ is the shaft torsional angle; ωt and ωr are the speeds of the rotor of the WT and the generator; ωs is the stator speed of the generator; and ωb is the base value of rotational speed.

2.3. Pitch Angle Control Model

To ensure the smooth power output of the OWFs, the pitch angle  needs to be controlled. The equations are given by(3)dβdt=1Tββrefβ,βref=Kp0dωtdt+K10Δωt,where Tβ is the inertial time constant of the pitch angle control model; βref and ωt_ref are the reference values of the pitch angle and the WT speed; and Kp0 and KI0 are the proportional and integral coefficients of the controller.

2.4. DFIG Model

Currently, DFIG is the most widely used WT. Ignoring the electromagnetic transient process of stator winding, the mathematical model of the DFIG  is as follows:(4)deddt=1T0ed+sωseqXsXsT0iqsωsLmLrruqr,deqdtsωsed1T0eq+XsXsT0ids+ωsLmLrrudr,udr=Rridr+TbdidrdtsωsTbiqr,uqr=Rriqr+Tbdiqrdt+sωsTauqs+Tbidr,ids=1LssuqsLmLssidr,iqs=LmLssiqr,where T0 = Lrr/Rr; Ta = Lm/Lss; Tb = Lrr − Lm/Lss; Xs = ωsLss; Lss = Ls + Lm; and Lrr = Lr + Lm are the sum of self-inductance of the stator and the rotor; Rs and Rr are the resistances of the stator and the rotor; Lm is the mutual inductance between the stator and the rotor; Xs is the reactance of the stator; ed = −(ωsLm/Lrr)ψqr; eq = (ωsLm/Lrr)ψdr; Xs = ωs/Lrr (LssLrrLm); ed and eq are the d- and q-axis components of the transient voltages; Xs is the transient resistance of the stator; and ψdr and ψqr are the d- and q-axis components of the rotor flux, respectively. ids, iqs, idr, and iqr are the d- and q-axis components of the stator and rotor currents, respectively. Let the stator flux ψs always reunite with the d-axis; then, uds = 0 and uqs = ψs. uds, uqs, udr, and uqr are the voltages of the stator and the rotor.

The active and the reactive power output of OWFs are(5)Pe=32RrPs2+Qs2Ta2us2+2RrQsTa2Lss+1sωPs+Rrus2Lm2,Qe=32sωTbPs2+Qs2Ta2us2+2QsTa2Lss+us2Lm2+1sωsQssωsus2Lss.

2.5. VSC-HVDC Model

According to the relationship between the three-phase stationary coordinate system and the d and q synchronous rotating coordinate system, after the Park transformation, the 7-order mathematical model [13, 21] can be obtained as follows:(6)L1di1ddt=R1i1dω1L1i1q+usd1K1ud1cosδ1,L1di1qdt=R1i1q+ω1L1i1d+usq1K1ud1sinδ1,L2di2ddt=R2i2dω2L2i2q+usd2K2ud2cosδ2,L2di2ddt=R2i2qω2L2i2d+usq2K2ud2sinδ2,Cdud1dt=3K12i1dcosδ1+i1qsinδ1id,Cdud2dt=3K22i2qcosδ2+i2qsinδ2id,Lddidt=ud1ud2Rdid,where the physical quantities corresponding to the rectifier and the inverter station are denoted by subscripts “1” and “2,” respectively. i1d, i1q, i2d, and i2q represent the current components of the AC network; usd1, usq1, usd, and usq2 represent the voltage components of the AC network; ud1, uq1, ud2, and uq2 represent the voltage components between buses of the DC system; and K1 and K2 are the voltage utilization coefficients of the DC system. δ1 and δ2 are the deviation angles between the voltage of the converter station and the bus voltage of the AC system; ω1 and ω2 are the fundamental angular frequencies of the AC system. id is the DC current transmitted by the high-voltage transmission line.

2.6. System Interface Model

Position the stator voltage us on the d-axis, and the voltage vector between the DFIG and the rectifier station is shown in Figure 2.

The voltage vector between the DFIG and the rectifier station.

Its mathematical relationship is(7)usd1usq1=udsuqs+0XTL1XTL10idsiqs,where XTL1 is the total impedance of line 1; XT0 is the equivalent transformer impedance of OWFs; and XT1 and XL1 are the transformer and the circuit impedance of line 1.

Position the terminal voltage u of the power system on the x-axis, and the vector between the inverter station and the power system is shown in Figure 3.

The voltage vector between the inverter station and the power system.

Its mathematical relationship is(8)usd2usq2=cosφsinφu0XTL2XTL20i2di2q.

2.7. Complete Model

All the equations of DFIG and VSC-HVDC are combined and linearized near the stable value:(9)dΔx10dt0=A˜B˜C˜D˜Δx10Δy10,where x10=ωt  θ  s  β  ed  eq  i1d  i1q  i2d  i2q  u1d  u2d  idT and y10=idsiqsT.

The offshore wind power transmitted by VSC-HVDC is integrated into the IEEE9 system, as shown in Figure 4. The small-signal equation of the system is(10)dΔxdt=AΔx,where x=x1x2x3x10T; x1=δ1ω1T; x2=δ2ω2eq2ed2efq2vR2vM2T; and x3=δ3ω3eq3ed3efq3vR3vM3T.

Schematic diagram of the equivalent offshore wind farm connected to the IEEE9 system by VSC-HVDC.

After deducing the complete small-signal stability model of offshore wind power transmitted by VSC-HVDC, it is necessary to do further research on the influence of small-signal stability on the system of offshore wind power in different positions and capacities.

3. Small-Signal Stability Analyses

This section will combine with the concrete example of the offshore wind power grid-connected position and capacity for small-signal stability analysis. Firstly, the parameters of each subsystem are initialized, which are detailed in Appendix and then the state of OWFs is found before incorporation into the system. Finally, two load nodes and two generator nodes are chosen, and the influence of the offshore wind power grid connection is analyzed in detail.

3.1. Eigenvalue Analysis

When the OWFs are added at node 6 and the active power output is 0.5 pu, the eigenvalues of the system are shown in Table 1. It can be seen from Table 1 that the real part of the eigenvalues is all negative, which are on the left side of the imaginary axis, indicating that the original power system is running in a stable state.

Eigenvalues of the system during stable state.

NumberEigenvalueNumberEigenvalue
λ1−99.9360λ2−74.4198
λ3−52.7251λ4−51.3813
λ5−31.1510λ6−28.8327
λ7,8−2.6655 ± 18.2668iλ9−12.9409
λ10,11−0.6959 ± 12.8937iλ12,13−0.1431 ± 8.6461i
λ14,15−4.7160 ± 8.0754iλ16,17−0.1222 ± 4.0611i
λ18,19−4.4655 ± 3.1380iλ20,21−5.4943
λ21−4.1027λ22−0.4476
λ23,24−0.8287 ± 0.9010iλ25,26−0.3560 ± 0.5849i
λ27−0.0290λ28,29−0.0291 ± 0.0502i

λ7,8, λ10,11, λ12,13, and λ16,17 are related to the generator rotational speed (Δω1, Δω2, and Δω3) and the rotational speed difference Δs of wind generator 10, respectively, belonging to the electromechanical oscillation mode.

3.2. Damping Ratio Analysis

Damping ratio  can reflect the speed and the characteristics of oscillation attenuation. The expression is given by(11)ζ=σσ2+ω2,where σ and ω are the real and the imaginary part of the eigenvalue, respectively.

In the actual power system, it is generally required that the damping ratio of the electromechanical oscillation mode should be above 0.05, and then the operating state of the system could be accepted. However, this principle is not unchanged. If the mode of fluctuation is not large when the system operation mode changes, it is acceptable to have such a low damping ratio (for example, 0.03) .

By changing different active power outputs and grid-connected positions of OWFs, the damping ratio curve in various cases can be obtained as shown in Figure 5.

The damping ratio of OWF grid connection.

4. Control Measures

In order to improve the weak damping instability of the system, it is necessary to add a controller to the system for auxiliary regulation. Firstly, PSS needs to be modeled, and the speed difference is used as the feedback signal to be added to the system voltage equation. Secondly, the optimal position of the PSS is determined by the detailed analysis of modal control theory.

Finally, the PSO algorithm is introduced, and the objective function is established according to the research content of this paper. The parameters of PSS controllers are determined with different inertia weight functions.

4.1. PSS Model

The main function of the PSS is to increase damping or suppress low-frequency oscillation of the power system . The mathematical expression of the PSS is(12)dΔV1dt=KgainT6Δω1T6ΔV1,dΔV2dt=KgainT6Δω1T6ΔV11T5ΔV2,dΔV3dt=KgainT1T2T6ΔωT1T2T6ΔV1T1T5T2T5ΔV2ΔV3,dΔVSdt=KgainT1T3T2T4T6ΔωT1T3T2T4T6ΔV1T3T1T5T2T4T5ΔV2T3T2T2T4ΔV31T4ΔVS,where the generator speed difference VIS = ω-ωt is selected as the input signal. VS is the output signal, which serves as the auxiliary input signal of the stator voltage equation Δuqi (i = 1,2,3) of the IEEE9 system. Kgain is the gain of the amplification link; T6 is the time constant of the measurement link. T5 is the time constant of the filtering link. T1, T3, T2, and T4 are the lead and lag time constants of two phase compensation links, respectively.

If the PSS is added to IEEE9 system generator 2, the voltage equation needs to be(13)Δud2=ΔedRa2Δid2+Xq2Δiq2,Δuq2=ΔeqRa2Δiq2+Xd2Δid2+ΔVS.

4.2. PSS Grid-Connected Position Selection

When OWFs integrate into load node 6 and in normal operation, the changes of the electromechanical oscillation mode are shown in Table 2.

Corresponding components of generator speed in the right eigenvector.

Generatorλ7,8λ10,11λ12,13λ16,17
Modulus value (pu)Argument (°)Modulus value (pu)Argument (°)Modulus value (pu)Argument (°)Modulus value (pu)Argument (°)
10.0000−65.650.0002−158.190.0009−108.950.0007−73.93
20.0001−50.390.0012−150.940.002469.350.0006−81.42
30.0003−57.680.004026.210.001372.600.0005−81.49
100.0115−42.980.004623.870.002572.360.0021−78.95

As can be seen from Table 2, in the eigenvectors of the first and fourth modes, the modulus value of Δs is the largest, but the directions of each component are basically the same, indicating that the influence of each generator on the mode is similar. In the second mode, the modulus value of Δω3 and Δs is larger and opposite to the direction of other components (the average argument is about −154.57°), indicating that the oscillation mainly exists in generators 3, 10 and generators 1, 2. The oscillation frequency f = 8.65/(2π) = 1.3767 Hz, belonging to the local oscillation mode. In the third mode, the modulus values of Δω2 and Δs are larger and are basically opposite to Δω1 in the direction (the argument is −108.95°), indicating that the oscillation mainly exists in generators 2, 10 and generator 1. The oscillation frequency f = 12.89/(2π) = 2.0515 Hz, belonging to the local oscillation mode.

In the participation vector, if the generator speed component is positive and the value is large, it indicates that the generator is the best candidate position for installing the PSS, which can significantly increase the damping of the mode. For increasing system damping, it is better to apply control on a larger capacity generator .

The local oscillation modes λ10,11 and λ12,13 are specifically analyzed. The components of the participation vector corresponding to the rotation speed of each generator are shown in Table 3.

Components of the participation vector corresponding to the rotation speed of each generator.

Generator12310
λ10,110.00990.19141.00000.0264
λ12,130.37001.00000.12200.0143

As can be seen from Table 3, for the participation vector of local mode λ10,11, the value of generator 3 is the largest and that of generator 10 is the smallest, indicating that the control is mainly applied in generator 3, and the damping of the system can be significantly increased. Similarly, for the participation vector of local mode λ12,13, the value of generator 2 is the largest and that of generator 10 is the smallest, indicating that fine effects can be achieved by applying control on generator 2.

Furthermore, compared with the capacity of generators 2 and 3, the active power of generator 2 (1.63 pu) is larger than that of generator 3 (0.85 pu), so the damping control of generator 2 can achieve better suppression of the oscillation effect.

4.3. PSO Algorithm Based on Different Inertia Weight Functions

The updating formulas of particle velocity and position are as follows:(14)viDk+1=ωviDk+c1r1piDkxiDk+c2r2pgDkxgDk,xiDk+1=xiDk+viDk+1,where k represents the number of iterations; viDk and xiDk represent the velocity and position of the i-th particle in the D-dimensional space, and the value range is [νlb, νub] and [xlb, xub]; and c1 and c2 refer to the individual learning factors and the group learning factors, and c1 = c2 = 2 is taken in this paper. r1 and r2 are random numbers between 0 and 1; piDk and pgDk represent the individual optimal position and the global optimal position of the i-th and g-th particle in the D-dimensional space; and ω represents the inertia weight.

The inertia weight ω indicates the ability to retain existing velocity. The larger the value is, the stronger the global optimization is. On the contrary, the smaller the value is, the better the local optimization is. In order to weigh the global optimization and local optimization, Shi. Y summarized a new method of calculating weights, linear decreasing inertia weights (LDIW) , which is specifically expressed as(15)ω1k=ωstartωstartωendkTmax,where ωstart is the inertia weight in the initial state; ωend is the inertia weight at the end of the iteration; k is the number of iterations; and Tmax is the maximum number of iterations. It is generally believed that ωstart = 0.9 and ωend = 0.4 in which the algorithm performs best. In addition to LDIW, the common inertia weight functions  are as follows:(16)ω2k=ωstartωstartωendkTmax2,ω3k=ωstart+ωstartωend2kTmaxkTmax2,ω4k=ωendωstartωend1/1+k/Tmax.

4.4. Determination of Optimization Objective Function

In order to improve the weak damping of the system, it is very important to determine the optimal target. As can be seen from Figure 5, when OWFs first join the system, the damping ratio is relatively large. With the increase of the active power output, most damping ratios show a downward trend, and some electromechanical oscillation modes even have weakly damped unstable states. In order to improve the low-frequency oscillation, when generator 10 active output is 1.0 pu, the minimum damping ratio of each generator is selected as the optimal target, and then the maximum value of this target in the iterative process can be obtained by the PSO algorithm. If the optimization target is larger than 0.03, the variation curve of damping ratio can be ensured to be in a stable running state. The mathematical expressions of fitness function are as follows:(17)Dn=minζni,fitness=maxDn,where ξni (i = 1,2,3,4) is the damping ratio of the component related to the rotation speed of generators 1, 2, 3, and 10 during the n-th iteration process.

According to the research of this topic, the steps of the PSO algorithm are as follows:

Initialize particle swarm size, velocity, and position. Particle specifically refers to PSS parameters.

Take the small-signal stability program as a subfunction, and then optimize the designed objective function to obtain the individual and the global optimal fitness value of the particle.

Compare the fitness value of each generation particle. If the current fitness value is better, update the previously recorded optimal fitness value with the current fitness value, and update the previously recorded optimal position with the current position. If not, leave the value unchanged.

Update the position and velocity of each particle.

Determine whether the current fitness value reaches the given standard or reaches the upper limit of the number of iterations. If yes, the program ends. If not, return to (2) and continue to execute the loop body.

4.5. Effect of Applying PSS Control

According to the analysis of Section 4.3, PSS installed in generator 2 works well. Under different inertial weights, the fitness curve can be obtained by solving the function for many times and taking the mean value, as shown in Figure 6.

Fitness curves under different inertia weights.

It can be seen from Figure 6 that the deviation of the optimal fitness obtained under the four inertia weights is small. The optimal fitness of the inertia weight ω2 is 0.034232, the effect of ω1 and ω3 is similar, and the worst is ω4. Based on ω2, the control parameters are deeply optimized, and the optimized parameters of each controller are Kgain = 10, T1 = 0.66, T2 = 0.25, T3 = 0.45, T4 = 0.15, T5 = 4.5, and T6 = 1.25. By substituting the PSS parameters into the original state matrix and changing the position and capacity of the offshore wind power into the power system, the effect can be obtained, as shown in Figure 7.

Curve of damping ratio change after the PSS is added into generator 2.

As can be seen from Figure 7, after adding the PSS to generator 2, each damping ratio is improved correspondingly, and the curves are all above 0.03, indicating that the system works stably. The variation trend of each damping ratio is basically the same as that without the PSS, and the damping ratio of generator 2 rotation speed is significantly increased. Furthermore, it can be verified that the theoretical analysis is reasonable, and generator 2 is the best place to add the PSS.

Through theoretical analysis and programming experiments, the control measures of offshore wind power grid connection have been obtained. However, how the actual power system operates and whether the controller PSS achieves the desired effect need to be verified by simulation. The circuit in Figure 1 was built on MATLAB/Simulink. By applying the small signal (wind speed disturbance and system fault), the proposed control measures are verified and analyzed.

5.1. Wind Speed Disturbance

Assuming the initial wind speed of the OWFs is 10 m/s, a step signal generator disturbs the wind speed to 12 m/s at 1.8 s. The dynamic response effect of the system before and after adding PSS control to generator 2 is shown in Figure 8.

Dynamic response of the system before and after adding PSS2 in the case of wind speed disturbance. (a) The electromagnetic torque of synchronous generator 1. (b) The voltage of synchronous generator 1. (c) The active power of synchronous generator 1. (d) The direct current voltage of HVDC. (e) The active power of VSC1. (f) The reactive power of VSC1.

From the PSS2 control in Figures 8(a) and 8(c), it can be clearly seen that the power of the system fluctuates within a small range after applying the wind speed disturbance, and the system runs stably again at 3.3 s. The process has gone through 1.5 s; it can be seen from Figure 8(a) that, after the addition of PSS2, the electromagnetic torque of synchronous generator 1 is stabilized more quickly, and the oscillation amplitude is effectively suppressed. Figure 8(b) is the voltage curve of synchronous generator 1, and the effect of adding PSS2 is small. It can be seen from Figure 8(c) that the amplitude of the synchronous generator 1 active power is significantly reduced after the addition of PSS2. Figure 8(d) shows the DC curve of HVDC. By comparison, the amplitude of the oscillation is reduced when PSS2 is applied. The amplitude of the voltage stability will drop slightly without PSS2 control, which is also the adverse effect of wind speed disturbance. Figures 8(e) and 8(f) show the active and reactive power of VSC1. It can be seen that, after the addition of the damping controller, the oscillation of the system is better improved. In summary, when the OWFs are disturbed by wind speed, adding PSS2 will improve the small-signal stability of the whole system.

5.2. System Fault

In order to simulate the actual system fault, a three-phase short circuit was applied to the transmission line incorporated into node 8 at 1.5 s, and the fault duration is 0.12 s. The dynamic response of the system controlled by PSS2 is shown in Figure 9.

Dynamic response of the system before and after adding PSS2 in the case of the three-phase short-circuit fault. (a) The electromagnetic torque of synchronous generator 1. (b) The voltage of synchronous generator 1. (c) The active power of synchronous generator 1. (d) The direct current voltage of HVDC. (e) The active power of VSC1. (f) The reactive power of VSC1.

It can be seen from Figures 9(a) and 9(f) that, after applying a three-phase short-circuit fault, the state of the system changes obviously. The system recovers at 2.8 s, and the adjustment process goes through 1.3 s; Figure 9(a) is the electromagnetic torque curve of synchronous generator 1, and the oscillation of the system is significantly reduced after the addition of PSS2. Figure 9(b) shows the voltage curve of synchronous generator 1, and the effect of adding or not adding PSS2 is small. It can be seen from Figure 9(c) that the amplitude of the synchronous generator 1 active power is remarkably suppressed after the addition of PSS2. Figure 9(d)shows the DC curve of HVDC. It can be seen that the DC voltage fluctuation amplitude is reduced, and the system reaches stability more quickly. Figures 9(e) and 9(f) show the active and reactive power of VSC1. It can be seen that, after the damping controller is added, the oscillation of the system is better improved. In summary, when the onshore power system is disturbed by the fault, adding PSS2 improves the small-signal stability of the whole system.

6. Conclusions

Based on the establishment of a complete mathematical model, the influence of offshore wind power on the system at different positions and capacities is analyzed in detail according to the damping characteristics. The experimental results show that most of the damping ratios are decreasing, and some of the electromechanical oscillation modes will show weak damping instability. By comparing the experimental results of two load nodes and two generator nodes, the optimal location and capacity of offshore wind farms are obtained. In order to improve the weak damping of the power system, the optimal positions of the PSS in the power grid are obtained through the modal control. PSO algorithm combining with different inertia weight functions is presented to determine the parameters of the PSS. Finally, the whole system model is built on MATLAB/Simulink platform, and the influence of PSS control on the system is compared and analyzed by applying wind speed disturbance and system fault. The simulation results show that the dynamic response of the system is obviously improved, the oscillation amplitude is significantly reduced, and the accuracy of the proposed control strategy is verified.

Appendix

DFIG parameters of offshore wind farms:

r = 40 m; ρ = 1.242 kg m3; v0 = 10.45 m s−1; Rr = 0.0073 pu; Rs = 0.0076 pu; Lr = 0.0884 pu; Ls = 0.1248 pu; Lm = 1.8365 pu; ωt0 = 1.1055 rad s−1; s0 = −0.003; β0 = 0 rad; Kp0 = 0.1 pu; KI0 = 20 s−1; Tβ = 10 s; Ht = 3 s; Hg = 0.05 s; K = 10 pu.rad−1; and D = 0.5 pu.rad−1.

VSC-HVDC parameters:

R1 = R2 = 0.075 pu; L1 = L2 = 0.016 pu; ω1 = ω2 = 1 pu; δ1 = δ2 = 10°; K1 = K2 = 0.1; Rd = 1.0425 pu; Ld = 0.0119 pu; and C= 0.1733 pu.

System interface parameters:

XT0 = 0.015 pu; XT1 = 0.016 pu; XL1 = 0.02 pu; XT2 = 0.061 pu; and XL2 = 0.014 pu.

Data Availability

All data generated or analyzed during this study are included in this article.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This paper was carefully guided by Professor Ru Yang of Guangzhou University and was supported by the National Natural Science Foundation of China under Grant nos. 61772147 and 61100150, the University Innovation Team Construction Project of Guangdong Province under Grant no. 2015KCXTD014, Guangdong Province Philosophy and Social Science Foundation under Grant no. GD19CSH03, the key Project of Science and Technology Plan of Guangdong Province under Grant No. 2020b1010010014 and Open Research Fund Program of Data Recovery Key Laboratory of Sichuan Province. Here, the authors express their sincere gratitude.

PinedaI.RubyK.HoA.MbistrovaA.CorbettaG.“The European Offshore Wind Industry—Key Trends and Statistics 20152016Brussels, BelgiumEuropean Wind Energy AssociationWalshC.MiróL.SelotF.FraileD.BrindleyG.“Offshore Wind in Europe – Key Trends and Statistics 20182019Brussels, BelgiumWind EuropeTardieuP.VeldeL.NghiemA.PinedaI.Wind Energy in Europe: Scenarios for 20302017Brussels, BelgiumWind EuropeShenL.BarnesM.PreeceR.MilanovićJ.Frequency stabilisation using VSC-HVDCProceedings of the 2016 IEEE Power and Energy Society General Meeting (PESGM)2016Boston, MA, USA15YanL.WeiS.Yong-NingC.Lin-JunW.WeiL.Research on offshore wind farm VSC-HVDC transmission system fault ride through issueProceedings of the 2014 International Conference on Power System Technology2014Chengdu, China21902195AsimeniaK.WuQ.ZhaoH.Review of VSC HVDC connection for offshore wind power integrationRenewable and Sustainable Energy Reviews2016591405141410.1016/j.rser.2016.01.0642-s2.0-84956971553EdrahM.LoK. L.Anaya-LaraO.ElansariA.Impact of DFIG based offshore wind farms connected through VSC-HVDC link on power system stabilityProceedings of the 11th IET International Conference on AC and DC Power Transmission2015Birmingham, UK17FuQ.DuW.CaoJ.WangH. F.“VSC-based HVDC power transmission for the large-scale offshore wind power—a surveyProceedings of the International Conference on Renewable Power Generation (RPG 2015)2015Beijing, China16LuS.XuZ.A novel control strategy of VSC-HVDC for offshore platformsProceedings of the 12th IET International Conference on AC and DC Power Transmission (ACDC 2016)2016Beijing China16GuoY.GaoH.WuQ.ZhaoH.OstergaardJ.ShahidehpourM.Enhanced voltage control of VSC-HVDC-Connected offshore wind farms based on model predictive controlIEEE Transactions on Sustainable Energy20189147448710.1109/tste.2017.27430052-s2.0-85028503042MachowskiJ.BialekJ.BumbyJ.Power System Dynamics: Stability and Control20082ndNew York, NY, USAWileyHuangH.MaoC.LuJ.WangD.Small-signal modelling and analysis of wind turbine with direct drive permanent magnet synchronous generator connected to power gridIET Renewable Power Generation201261485810.1049/iet-rpg.2010.02172-s2.0-84863160194YangY.LiG.The small signal stability control of offshore wind farm based on VSC-HVDCTransactions of China Elect Rotechnical Society20163113101110SajadiA.ZhaoS.ClarkK.LoparoK. A.Small-signal stability analysis of large-scale power systems in response to variability of offshore wind power plantsIEEE Systems Journal201813311010.1109/JSYST.2018.28853022-s2.0-85058878311RazaM.Prieto-AraujoE.Gomis-BellmuntO.Small-signal stability analysis of offshore AC network having multiple VSC-HVDC systemsIEEE Transactions on Power Delivery201833283083910.1109/tpwrd.2017.27259822-s2.0-85023624924WangC.FangH.HeS.Adaptive optimal controller design for a class of LDI-based neural network systems with input time-delaysNeurocomputing202038529229910.1016/j.neucom.2019.12.084HeS.FangH.ZhangM.LiuF.DingZ.Adaptive optimal control for a class of nonlinear systems: the online policy iteration approachIEEE Transactions on Neural Networks and Learning Systems202031254955810.1109/tnnls.2019.2905715HeS.FangH.ZhangM.LiuF.LuanX.DingZ.Online policy iterative-based H∞ optimization algorithm for a class of nonlinear systemsInformation Sciences201949511310.1016/j.ins.2019.04.0272-s2.0-85065212766HeS.ZhangM.FangH.LiuF.LuanX.DingZ.Reinforcement learning and adaptive optimization of a class of Markov jump systems with completely unknown dynamic informationNeural Computing and Applications201910.1007/s00521-019-04180-2WuF.ZhangX.-P.GodfreyK.JuP.Small signal stability analysis and optimal control of a wind turbine with doubly fed induction generatorIET Generation, Transmission & Distribution20071575176010.1049/iet-gtd:200603952-s2.0-34547797926TangG.HVDC Transmission Technology Based on Voltage Source Converter2010Beijing, ChinaChina Electric Power PressQinC.YuY.Small signal stability region of power systems with DFIG in injection spaceJournal of Modern Power Systems and Clean Energy20131212713310.1007/s40565-013-0023-12-s2.0-84902655206LiZ.ZhangX.DengY.HeB.Studies of transient stability of wind power acess grid systemJ of China Three Gorges Univ. (Natural Sciences)20163824750MehtaB.BhattP.PandyaV.Small signal stability analysis of power systems with DFIG based wind power penetrationInternational Journal of Electrical Power & Energy Systems2018586474PrasanthiE.ShubhangaK. N.Stability analysis of a grid connected DFIG based WECS with two-mass shaft modelingProceedings of the 2016 IEEE Annual India Conference (INDICON)2016Bangalore, India16FengS.LiuD.YangX.XiangX.The impact of wind power inverse-peaking characteristics on power system low frequency oscillationProceedings of the 2012 China International Conference on Electricity Distribution2012Shanghai, China14HuangH.Research on Small-Signal Stability of Wind and Photovoltaic Power Generation System2012Wuhan, ChinaDepartment Electronical Engineering, HUST UniversityPh.D. dissertationWangX.FangW.DuZ.Modern Power System Analysis2003Beijing, ChinaScience PressThanpisitK.NgamrooI.Power oscillation damping control by PSS and DFIG wind turbine under multiple operating conditionsTurkish Journal of Electrical Engineering & Computer Sciences2017254354436810.3906/elk-1702-1902-s2.0-85053823292GautamD.VittalV.Impact of DFIG based wind turbine generators on transient and small signal stability of power systems,Proceedings of the 2009 IEEE Power & Energy Society General Meeting2009Calgary, Canada16XinJ.ChenG.HaiY.A particle swarm optimizer with multi-stage linearly-decreasing inertia weight,Proceedings of the 2009 International Joint Conference on Computational Sciences and Optimization2009Sanya, China505508YuL.ShiF.WangH.HuF.Analysis of 30 Cases of Intelligent Algorithms in MATLAB2011Beijing, ChinaBeihang University Press