1. Introduction
Efforts to reduce carbon dioxide (CO
2) emissions have resulted in a focus on road transportation, which accounts for nearly 17% of global emissions [
1,
2]. Electric vehicles (EVs) are a promising solution, but their net zero emissions depend on the energy mix used to generate electricity [
1]. Using renewable energy sources for powering the batteries of EVs can significantly reduce net CO
2 emissions [
1]. It is necessary to increase the proportion of renewable electricity used to power EVs as the propulsion systems improve [
1,
2].
Wind energy has experienced a significant global surge, reaching 906 GW by the end of 2022, thanks to technological advancements and its reliability [
3,
4,
5,
6]. Variable Speed Wind Turbines (VSWTs) utilize both Doubly Fed Induction Generators (DFIG) and Permanent Magnet Synchronous Generators (PMSG). These generators are integrated into the electrical grid through power electronic converters and are referred to as type 3 (DFIG-based VWST) and type 4 (PMSG-based VWST) wind turbines [
7,
8]. PMSGs are preferred over DFIGs due to their lack of gearboxes, grid decoupling, high power density, and simple control mechanisms [
3,
4,
5,
6,
7,
8]. The availability of different PMSG types, including axial and radial flux configurations, allows for customized designs for various applications [
7]. Axial Flux PMSGs are suitable for space-constrained offshore wind farms, while Radial Flux PMSGs are commonly used in onshore direct-driven wind turbines due to their simplicity and reliability [
7].
Efficient wind energy conversion requires a proper VWST selection and an optimal wind farm location [
4]. However, wind power is available primarily remotely, posing various challenges to harvest and integration into the electrical grid [
8]. One challenge is upgrading the grid infrastructure to accommodate increased power flows and ensure reliable electricity delivery [
8,
9,
10,
11]. A solution to increase transmission line capacity is to use series compensation equipment, which may cause Sub-Synchronous Resonance (SSR) due to the interaction between torsional oscillation and resonance frequencies [
8,
9,
10,
11].
SSR is a resonant condition where the electrical and mechanical components of the wind farm exchange energy with the electrical network at frequencies below the fundamental frequency [
12]. If not prevented, these resonances can severely threaten the generating stations, including wind farms [
13,
14,
15].
Firstly, rapid initiation and escalation of SSR can cause over-voltages that damage vital components like series capacitors, crowbar circuitry, and connected loads [
16]. This event was first reported by the Electric Reliability Council of Texas (ERCOT) in 2009 [
12,
17,
18]. A faulted line and subsequent outage in the electrical network caused an interaction between a wind farm and the compensated network. This event damaged the compensation equipment and wind turbines due to the production of SSR [
19]. Investigations concluded that the instability was a consequence of the exchange of electric energy between power electronic converters of the wind turbines and the compensated line [
8,
11], later referred to as Sub-Synchronous Control Interaction (SSCI) [
8,
11].
Secondly, the unstable sub-synchronous oscillations resulting from SSCI can lead to power generation losses [
16]. The over-voltages at wind generators (WG) and the Point of Common Coupling (PCC) trigger the tripping of WGs, significantly reducing the net power output of the wind farm [
16]. In 2017, these incidents were also captured by wind farms connected to the North China Power Grid [
19,
20]. At Guyuan, the power output dropped from 219.5 MW to 74.5 MW within 108 s during an SSCI event [
19]. Three 600 MW steam turbine generators abruptly tripped at Hami due to torsional relays during an SSCI event [
16]. Such reductions in power output undermine the dependability and reliability of the entire power system.
Thirdly, SSCI introduces inter-harmonic and subharmonic components into the system [
16]. The presence of inter-harmonic and subharmonic components can negatively affect power quality and equipment performance [
16]. This includes issues like loud noise and efficiency reduction in grid-connected equipment and electrical loads [
16]. Thus, it is imperative to develop effective SSR countermeasures without delay.
Initially, researchers believed that direct-driven PMSG-based type 4 wind farms were immune to SSR due to grid isolation [
8,
9,
10,
11,
20,
21,
22,
23]. However, recent studies [
20,
22,
23] have shown that the risk of SSR exists in these wind farms. Furthermore, Lui et al. [
20] claims that due to weak grid conditions at the PCC the SSCI can occur. Also, Beza et al. [
22] shows that the closed-loop bandwidths of the Phase Locked Loop (PLL) and the reactive power controller can cause SSCI resulting in system stability. Additionally, Xu et al. [
23] shows that the series compensation levels can lead to oscillatory stability issues under weak grid conditions. In this regard, the control parameters of the power electronic converter on the grid side play a dominant role in system stability. Beza et al. [
22] suggests that the open-loop control strategy should be adopted when encountering SSCIs. However, Lui et al. [
20] and Xu et al. [
23] do not pose any recommendations for the control strategy. These findings clearly emphasize the necessity of extensive analysis and development of techniques to suppress the phenomenon of SSR in type 4 wind farms.
In the field of SSR mitigation for type 4 wind farms [
24], stands as a pioneering work in developing a PLL-based Sub-Synchronous Resonance Damping Controller (SSRDC) [
25,
26,
27]. The innovative approach in [
24] focuses on suppressing both Sub-Synchronous and Super-Synchronous Resonances, contributing to advancing SSR mitigation techniques tailored explicitly for type 4 wind farm configurations. The study demonstrates the critical influence of the PLL bandwidth on system stability by performing a sensitivity analysis of the dominant poles of the impedance closed-loop transfer function. Optimizing the PLL bandwidth can substantially impact mitigating SSR and improving the system’s stability. Similar results regarding the influence of PLL bandwidth on system stability and the risk of SSR in type 4 wind farms were also recorded by Beza et al. [
22].
The PLL is a control system that synchronizes its output frequency and phase with its input using nonlinear negative feedback [
28,
29,
30]. The Synchronous Reference Frame PLL (SRF-PLL) consists of three stages [
28]:
A signal containing the phase error information is created by PD. Driving the phase error to zero is the responsibility of LF. At the same time, the VCO generates a synchronized unit vector that contains phase information at the output [
28]. A good PLL should be computationally light, more resilient to grid disturbances, quick to react dynamically, and stable [
28,
29,
30].
Various researchers have designed numerous upgraded PLLs with improved disturbance rejection capacity [
28,
29,
30]. Most of these PLLs can be viewed as a standard SRF-PLL with extra filters added, which may be placed either inside the SRF-PLL control loop or at its input [
28]. Both [
22,
24] investigate the impact of PLL bandwidth on SSR while using the conventional SRF-PLL and propose that the PLL bandwidth can substantially impact mitigating SSR and improving the system’s overall stability. However, they do not examine the impact of already existing upgraded PLLs on SSR. It is vital to investigate if the risk of SSCI is only triggered by the bandwidth of the PLL or whether architectural variations in the PLL also play a role in developing SSCI. A thorough analysis of improved PLLs is also necessary to understand the predominance of PLL layouts that can reject disturbances in the event of SSCI.
This paper aims to explore the influence of different PLL architectural variations on the development of SSCI in type 4 wind farms, employing a frequency-based analysis. The study investigates a variation of the conventional SRF-PLL by incorporating a Notch Filter (NF) within the Loop Filter (LF) of the PLL. The system’s passivity is then analyzed using eigenvalue analysis, followed by an examination of system stability using the Nyquist Stability Criterion. To gain insights into the influence of PLL architecture on SSCI, the obtained results are compared to those of the conventional SRF-PLL and the SSRDC-based PLL, as proposed by [
22,
24], respectively.
2. System Modelling and Description
Figure 1a shows a typical Permanent Magnet Synchronous Generator (PMSG)-based type 4 wind farm connected to a series compensated grid. In this system, the complete dynamics of the power system includes wind turbine aerodynamics, shaft system, PMSG, Stator Side Converter (SSC) and Grid Side Converter (GSC) converter controllers, DC-link model, Filter, Transformer, and series compensated transmission network dynamics [
22]. However, in regular operation, type 4 wind farms, the generator side is entirely decoupled from the electrical grid side due to the full-power electronic converter.
Therefore, it is convenient to neglect the dynamics of the Stator Side Converter, generator, and wind turbine by replacing it with a variable current source, as shown in
Figure 1b. The GSC is connected to the grid at the PCC via an RL filter having an internal inductance
and an internal resistance
. The signals
,
and
are the space vectors of the GSC output voltage, the filter current and the grid voltage measured at the PCC, respectively. Also, the direct voltage,
, at the input side of the converter is measured across the dc capacitance,
.
Figure 1b shows the system’s control algorithm comprises a cascaded structure with an inner current controller and two outer controllers, a DC voltage controller, and Reactive Power Controller, as shown in
Figure 2. The control system is implemented in the rotating dq frame of reference obtained by aligning with
, the grid voltage vector measured at the PCC. This alignment occurs by tracking the grid voltage angle,
, by the Phase-locked Loop (PLL). The reference voltage vector,
, is the current controller output provided to the PWM modulator.
The following subsections provide a complete design of the control system and the wind farm’s admittance modelling.
2.1. Inner Current Controller
The system equation for the GSC terminal voltage can be written as [
22,
23,
31]:
where,
, is the grid frequency. Then the control law for the inner current controller can be given as [
22]:
where ‘*’ annotates the reference signal.
and
represent the Proportional–Integral (PI) regulator of the current controller and the low-pass feed-forward filter of the grid voltage, respectively.
is the fundamental grid frequency estimated by the PLL. The Proportional and Integral gains
and
can be selected as
and
, respectively, whereas
is the closed-loop bandwidth selected for the current controller. The cut-off frequency of the feed-forward low-pass filter
can be selected as
.
The converter output voltage can be given as [
22]:
where,
is the modulator transfer function given as [
22]:
where,
represents the computation time delay due to the discretization of the output reference voltage,
, of the inner current controller. Then, the closed-loop dynamics of the inner controller in scalar form can be written as:
where,
represents the transfer matrix from reference currents to actual currents while
represents the transfer matrix from reference currents to the input admittance matrix. These matrices are mathematically given as [
22]:
2.2. Outer Controllers
The outer control comprises of DC voltage controller and the Reactive Power Controller. The two controllers are responsible for the power flow to the grid by providing reference current signals to the inner current controller which are described below.
2.2.1. DC Voltage Controller
The dynamics of the direct voltage,
, can be written as [
22,
31]:
where,
is the power from the Stator Side Converter or the power from the variable current source on the dc side of the GSC and can be written as:
While
is the power output of the GSC to the grid and is given as:
Then, the control law governing the dc voltage controller can be written as [
22]:
where ‘*’ annotates the reference signal.
and
represent the PI regulator for dc voltage control and the low pass filter for the magnitude of the grid voltage,
. The magnitude of the grid voltage can be given as:
From the control law, it can be inferred that the power from wind turbine, , acts as a feed-forward term filtered by the filter .
Considering, as the closed loop bandwidth for the DC voltage controller, the Proportional gain, , can be selected as . An Integral with a small integral gain may be added to the system to remove any steady-state errors caused by the mismatch of the analytical model to the actual system. The cut-off frequency of the low pass filter, , is selected as .
2.2.2. Reactive Power Controller
The dynamics of the reactive power injected into the grid at the PCC,
can be written as [
22,
31]:
Then the control law governing the reactive power control can be given as [
22]:
where ‘*’ annotates the reference signal.
and
represent the PI regulator for the reactive power control and the low-pass filter transfer function for grid voltage magnitude. The parameter selection for the PI regulator of the reactive power controller depends on the grid impedance [
22,
32]. However, no design recommendations are provided by [
22] for the reactive power controller.
In this study, only an integral controller is utilized, where the Integrator gain, is selected as . Where is the closed loop bandwidth of the reactive power controller. Also, the cut-off frequency of the low pass filter, is selected as .
2.3. Small Signal Modelling
The small signal model of the outer DC voltage controller and the reactive power controllers can be collectively written as [
22]:
where,
,
and
are considered as the inputs and the transfer matrices
,
,
are given as:
where,
Note that the ‘0′ in the variables’ subscript annotates the respective signals’ steady state values. While the ‘
’ annotates the perturbation in the respective signals.
is the steady state value of the grid voltage and is given as:
Then the combined dynamics of the inner and outer controllers can be written as [
22]:
where, the transfer matrices
and
are given as:
2.4. Phase-Locked Loop
Phase-locked Loop (PLL) is a control system that synchronizes its output to its input [
28,
29,
30] and is shown in
Figure 3. A PLL is a closed-loop control system designed to bring the q component of the grid voltage at the PCC to zero [
28,
29,
30]. This process results in aligning the grid voltage vector to the d axis of the dq frame of reference to realize the voltage-oriented control in the dq coordinate system.
The output of the PLL is the information of the grid angle
which coincides with the true grid angle
when the q component of the grid voltage at the PCC becomes zero. As
defines the converter dq frame and
defines the grid dq frame then the grid voltage vector at the PCC relates to the two frames by the following equation [
22,
32]:
where,
and
are the grid voltages at PCC in the converter dq frame of reference and the actual grid dq frame of reference respectively. It is evident from the above equation that when
coincides with
then
becomes equal to
. This transformation is true for any voltage or current vector of the system. The estimated angle by the PLL can be written as [
22,
32]:
where,
is the open loop transfer function of the PLL and is given as:
where,
and
are the proportional and integral gains of the PLL respectively. If
is selected to be the closed loop bandwidth of the PLL then the gains are selected to be
and
. Then the linearized small signal model of the PLL can be given as [
22]:
where,
Then, the dynamics of the entire system in the grid dq frame can be summarized as [
22]:
where,
is the desired input admittance matrix of the type 4 windfarm and can be written as:
3. System Passivity Analysis
When a system oscillates, it exhibits resistive behavior at its characteristic frequency. This is known as resonance. If the real part of the exhibited resistance is positive, the system dissipates energy, is stable, and is called passive. If the real part of the exhibited resistance is negative, the system accumulates energy and becomes unstable [
29]. Passivity tests can, therefore, be a useful tool for evaluating the stability of our system.
The passivity of the system can be demonstrated if the input admittance matrix,
, of the type 4 wind farm is positive indefinite i.e.,
. This holds true if and only if the eigen values of the
are positive [
22]. This means that if the eigen values of
are positive for all the frequencies below the fundamental frequency of the system then the risk of Sub Synchronous Control Interaction (SSCI) does not exit. On the contrary SSCI can occur if the eigen values of
are negative for all the frequencies below the fundamental frequency.
For the input admittance matrix,
, the parameters are assumed following the criteria presented in [
32] and are given in
Table 1.
While the active power and the reactive power outputs of the wind turbine are se-lected to be 1 p.u. and 0 p.u., respectively, the system eigenvalues are analyzed by in-dividually varying the bandwidths of the inner controller, outer control systems and PLL according to the recommendations provided by [
22,
32]. These values are given in
Table 2.
3.1. Analyzing the Behavior of Eigenvalue for Different Bandwidth Variations
Three different values of 2 p.u, 4 p.u. and 6 p.u. were used for the current controller bandwidth
. While the values of 0.1 p.u., 0.2 p.u. and 0.4 p.u. were used for the bandwidth of each control system in the outer loop i.e., DC Voltage controller, reactive power controller and the PLL as shown in
Table 2. The results obtained are shown in
Figure 4. It can be observed that the eigenvalue
stays positive for all the values of the investigated bandwidth. This shows that
does not influence the system passivity in the concerned range of frequency. Hence,
does not pose any risk of SSR and can be neglected during further investigation and analysis. These findings are consistent with the results reported by [
22].
3.2. Analyzing the Behavior of Eigenvalue for Different Bandwidth Variations
The values of the bandwidths of the control system deployed were varied in the system as shown in
Table 2. The results obtained for the eigenvalue
are shown in
Figure 5. It can be observed that the eigenvalue
becomes negative below the grid frequency. This means that the risk of Sub-Synchronous Control Interactions (SSCI) exists.
It can be observed from
Figure 5 that
has an unnoticeable impact on the behavior of
. This shows that for frequencies higher or lower than the grid frequency, the DC voltage controller has no effect on the system’s passivity. The passivity of the system is less affected by the current controller bandwidth,
, at frequencies in the Sub-Synchronous range. This indicates that the system passivity is slightly impacted by the current controller bandwidth as well. In contrast, the eigenvalue,
, in the sub-synchronous frequency range is more influenced by the reactive power controller bandwidth,
, and the PLL,
. These facts allow it to be confidently concluded that the PLL and the reactive power controller have a negative effect on the system’s passivity. These results are once more consistent with those reported by [
22].
Further investigation of
Figure 5 shows that increasing the PLL bandwidth not only increases the frequency range exhibiting non-dissipative behavior but also the magnitude increases at the saddle points of
. This means the PLL has the most substantial impact on the system passivity in the sub-synchronous frequency range. Hence, the PLL needs to be extensively investigated and improved to enhance its performance in the sub-synchronous frequency range.
4. PLL Enhancement by Sub-/Super-Synchronous Damping Controller
It has been established that PLL is essential for the stability of the Wind Farm [
22,
24]. The analysis presented in the above section shows that a typical PI based PLL depreciates the passivity of the system and needs to be extensively studied to enhance the system passivity. The research in [
24] proposes that an adaptive Sub-/Super-Synchronous Damping Controller (SSRDC) can be utilized to overcome the Sub-/Super Synchronous Oscillations (SSO) that arise within the system. It is suggested in [
24] that when an SSO arises, a direct-driven PMSG-based wind farm is subjected to a super-synchronous and a coupled sub-synchronous frequency component. If
is the Super-Synchronous frequency, then the coupled Sub-Synchronous frequency will be
where
is the fundamental grid frequency. In dq- frame of reference these coupled frequencies, the Super-Synchronous frequency, and the Sub-Synchronous frequency, can be visualized as
and
respectively. These two frequencies serve as the High Cut-off and Low Cut-off frequencies of the ABPF and are annotated as
and
respectively.
The SSRDC based enhanced PLL is shown in
Figure 6, where an Adaptable Bandpass Filter (ABPF) along with a Super-Synchronous Frequency Estimator (SSFE) and a tunable gain is used in parallel to the PI controller of the PLL. The frequency of the SSO i.e.,
, is estimated by the SSFE. Then the transfer function of the PLL can be written as [
24]:
where,
is the Super-Synchronous Oscillation frequency obtained from the SSFE.
is the damping coefficient of the ABPF.
is the gain of the SSRDC and
is the gain of the ABPF. To address the passivity of the system we neglect the adaptability of the SSRDC and vary the parameters of the SSRDC and PLL individually. The parameters of the SSRDC are selected to the following pivot values as per the recommendations provided by [
24] and are given in
Table 3.
Impact of SSRDC Based PLL Parametric Variations on the System Eigenvalues
As, the eigenvalue
showed no impact on the system passivity in
Figure 4. Hence, the eigenvalue
was investigated, and the results are shown in
Figure 7. As per the recommendations provided by [
24] three different values of 1.2 p.u., 1.4 p.u. and 1.8 p.u. were selected for
. Similarly, three different values of 0.2 p.u., 0.4 p.u. and 0.8 p.u. were selected for
and three different values of 0.02 p.u., 0.042 p.u. and 0.078 p.u. were selected for
to investigate the impact of SSO frequency on the eigenvalues of
. It can be observed that none of the three parameters
,
or
depict any noticeable influence on the eigenvalue
. This means that the SSRDC-based PLL does not suppress the risk of Sub-Synchronous Control Interactions (SSCI). Lastly, the PLL bandwidth,
, was varied as per the values given in
Table 2. It was again observed in
Figure 7 that the eigenvalue
exhibited a similar behavior as was seen with the system with typical PLL. Thus, it is imminent that the proposed SSRDC-based PLL in [
24] has no impact on the system passivity, and hence, the risk of SSR still exists.
5. PLL with Enhanced Filtering
Addition of filters to the PLL to enhance its disturbance rejection capability has gained a high importance in the past few years [
28]. Filter addition to the control loop or even before the input have been generally prescribed. Amongst these filters Moving Average Filter and Notch Filters have gained a huge acknowledgment for disturbance rejection due to power quality issues or synchronization challenges [
28]. However, these filters are usually utilized to perform better in rejecting the higher-order harmonics occurring in the system. Thus, here we propose to utilize a notch filter to suppress SSCI.
A notch filter (NF) is a filter that significantly attenuates signals within a narrow band of frequencies, called the rejection band, and passes all other frequencies with negligible attenuation. This makes it an ideal candidate for canceling selective or desired frequency components in the PLL control loop. It can be either adaptive or nonadaptive. To keep the analysis simple, we utilize a nonadaptive NF. The generic transfer function of a NF is given as [
29]:
where,
is the center frequency of the NF and
is the damping coefficient of the NF. Then the transfer function of the PLL with additional filtering capability, shown in
Figure 8, can be written as:
From
Figure 5 it can be analyzed that range of frequencies where the system passivity is jeopardized is up to 45 Hz. Also, from the theory of SSR, we know that these resonances occur between 4 Hz and 46 Hz [
7,
24,
25,
26,
27]. Therefore, we select
= 0.5 p.u., i.e., 25 Hz, while the damping coefficient is
= 0.4 p.u., and the passivity of the system is reanalyzed in the following sections.
5.1. Impact of NF Bandwidth
Three different values of 0.25 p.u., 0.5 p.u. and 0.75 p.u. were selected for
to investigate the impact of NF bandwidth on the eigen values of
. The results are shown in
Figure 9 and
Figure 10. It can be observed that the eigenvalue
shows significant improvement as compared to the system behaviour with typical PLL shown in
Figure 5. It can also be observed that
shows a good reduction in the frequency range with negative eigenvalues. This means that the risk of SSCI does exist but in frequencies below 20Hz. It can also be observed that the increasing bandwidth again increases the frequency range exhibiting non-dissipative behavior; however, magnitude decreases at the saddle points of
.
5.2. Impact of NF Damping Coefficient
Three different values of 0.2 p.u., 0.4 p.u. and 0.6 p.u. were selected for
to investigate the impact of damping coefficient on the eigen values of
. The results are shown in
Figure 9 and
Figure 10. It can be observed that eigenvalue
shows a slight change with the changing values of the damping coefficients. It can be observed that
shows a sleek reduction in the frequency range with increasing coefficient values; however, magnitude increases at the saddle points of
. Again, it is evident that the risk of SSCI does exist but in a thin frequency range.
5.3. Impact of NF on PLL Bandwidth
Three different values of 0.1 p.u., 0.2 p.u. and 0.4 p.u. were selected for
while
and
are set to their pivot values to investigate the impact of PLL bandwidth with NF on the eigen values of
. The results are shown in
Figure 9 and
Figure 10. It can be observed that the eigen value
, shows significant improvement as compared to the system behavior with typical PLL shown in
Figure 5. It can be observed from
Figure 10 that the frequency range with non-dissipative properties is limited to 20Hz. This means that the risk of SSCI exists but in frequencies below 20Hz. It can also be observed that the increasing bandwidth again increases the frequency range exhibiting non-dissipative behavior; however, magnitude increases at the saddle points of
.
8. Conclusions
In this paper, this risk of SSCI in type 4 wind farms connected to series compensated transmission lines due to PLL is explored in detail. In this context, a detailed admittance model of the type 4 wind farm is first provided then the impact of different control parameters on the passivity of the system is explored. It was observed that the bandwidths of the reactive power controller and the PLL show a higher risk of SSCI due to a non-dissipative behavior in the sub-synchronous frequency range. Furthermore, a typical PLL was then modified using an SSRDC, and a Notch Filter was employed separately in the PLL loop to explore their impact on the system passivity. It was observed that SSRDC-based PLL showed no significant improvement in the system passivity in the sub-synchronous frequency range. However, the PLL with Notch Filter showed a considerable improvement in passivity by restricting the non-dissipative nature of the system to less than 20 Hz. Finally, the system’s stability was analyzed for a typical PLL, SSRDC-based PLL and PLL with Notch Filter. It was again observed that the SSRDC-based PLL deteriorated the system stability margins, increasing the risk of SSCI in the sub-synchronous frequency range.
On the contrary, PLL with Notch Filter showed superior performance by increasing the system stability margins. In the future, it is recommended to further investigate the impact of PLL architectural variations in PLL Phase Detector (PD) and Loop Filter (LF) to thoroughly understand the role of PLL in the domain of SSR in wind farms. Also, optimization of PLL is recommended to suppress/mitigate the SSR in grid-integrated, direct-driven PMSG-based wind farms.