Harmonic Resonance Analysis for Wind Integrated Power System and Optimized Filter Design

As the contribution of renewable energy sources is increasing year over year, the effect of harmonics on power system becomes important, and it requires special attention. In conventional power sources, the harmonics is not generated at the source side; only load side is contributing in the harmonics. But renewable energy sources, particularly wind and solar, are based on power electronic devices, so it generates harmonics. This harmonics may have an adverse effect on the system. Harmonic resonance is one of the phenomena, due to which the harmonics are amplified and give rise to several trivial issues. Various methods are used to control the harmonics in the system. Harmonic filter is one of the simple ways to absorb the harmonics generated at load and generation side. Various filter designs have been found in literature as well as in the field. The filters are classified according to their design, construction and operation method. There are two main categories, active filters and passive filters. The passive filters are widely used due to its simplicity and lesser cost. However, to achieve the better performance, it is also used with active filters, and this combination is known as hybrid filter. The response of filters is modified as per the system requirement using various techniques. In this work, the impedance characteristics of various filters are discussed and analysed. Also, how the control structure of power electronic devices affects or modifies the output impedance of converter is also discussed.


Introduction
Voltage source converter (VSC) is widely used in industrial, commercial and renewable power generation applications. Out of many possible configurations, the three-phase and three-wire is widely used. During the last decade, the penetration of renewable energy sources has increased around the world. This is due to the increased concern worldwide about the carbon emission [1,2]. The utilisation of energy from uncertain and variable sources has become possible with the use of these converters only.
For the purpose of this study, the analysis of VSCs is carried out assuming the ideal operating conditions of grid. But VSCs are never operating under such conditions in practical. In a relatively weak system, VSCs are subject to various power quality disturbances, such as unbalance voltage, voltage swell and swag, notches, etc. [3]. The occurrence of such disturbances causes various problems like ripple in torque of generator and motor, increased losses, abnormal tripping of protective devices, malfunction of sophisticated control system, reduction in the expected lifetime of equipment, etc. [4].
There are two types of harmonics generated by VSCs. One is characteristics harmonics, which are related with the switching operations of the IGBTs inside the VSC. And second is non-characteristics harmonics. The voltage ripple on DC side of VSC generates harmonics on its AC side current. According to [5], the noncharacteristics harmonics are generated by the unbalanced voltage in the AC side. However the quantification of magnitude of such harmonics is not simple and cannot be done with deterministic method. The non-characteristics harmonics are considered as the steady-state low-frequency components which would not appear if the grid voltage is balanced. The unbalance grid voltage has fundamental frequency negative sequence component and third-order positive sequence component.
Though it is possible to eliminate zero sequence third harmonic component using transformer of proper vector group, the non-characteristics third-order positive sequence harmonics cannot eliminate transformers with delta-connected winding.
In [6], the author has proposed DC voltage control to eliminate DC oscillating voltage when AC side is unbalanced. To achieve this, VSC has to operate with constant AC power control. However, the effect of control on AC current is not discussed. It is important to analyse the distorted and unbalanced AC side current, when such control is implemented. In this case, the currents of AC side of converter contain non-characteristics low-frequency component such as fundamental negative sequence and third harmonic components of positive and negative sequence.

Modelling of the voltage source converter
The total DC current I dc is flowing to the DC side of VSC. V dc is the instantaneous voltage across the capacitor, and id is the instantaneous DC current. The quality factor of DC capacitor is assumed to be high, so the series resistance is neglected. The instantaneous power is supplied by the renewable energy source, i.e. wind turbine. The instantaneous current I s is supplied by the external source. It is equal to zero when VSC operates as a reactive power compensator.
The block diagram of VSC control structure is depicted in Figure 1. The control system consists of (i) voltage control, (ii) a phase lock loop (PLL), (iii) current reference calculation block and (iv) current control. The real power reference is calculated by the PI controller, which considers the DC voltage and desired active power through the VSC to the grid. The instantaneous reactive power is calculated by the separate loop, which may consider either the desired power factor or the reference voltage. The dq frame is synchronised with the positive sequence fundamental voltage of the grid at PCC with the help of PLL. It converts three phase voltages V a , V b and V c into V d + and V q +, which are converted in the alpha-beta reference frame voltages V α and V β . Then current reference is obtained by the α-β to abc transformation. These reference currents are compared with actual current, and modulating signals m d and m q are generated. These signals are finally transformed into m a , m b and m c to generate pulse width modulated (PWM) signal.
Both switching frequency and grid voltage distortion can cause poor power quality. A filter design is a subject that requires trade-off between filter performance and the control bandwidth. Filters are required to meet power quality standard, avoid parallel resonance and improve power quality.
Inverter for grid interfacing will need to incorporate interface filters to attenuate the injection of current harmonics.

Impedance-based stability analysis
This method was proposed in [7]. The system impedance is partitioned into source and grid impedance. The source impedance is either represented by Thevenin's equivalent circuit or Norton equivalent circuit (Figures 2 and 3). Thevenin's equivalent circuit consists of ideal voltage source in series with the series impedance (Z s ), whereas the load impedance is modelled by series impedance (Z l ). Since the converter circuit is non-linear, it is represented by the small signal circuit. This linear representation of circuit is valid only for the small perturbation of signal. With this assumption, the current (I) flowing from source to load is given by  Which is rearranged as, System stability analysis is based on the assumption that the source voltage and load impedance remain stable. So, V(s) and 1/Z(s) are stable. So, the stability depends on the extreme right-hand side of Eq. (2). It is given by The close observation of Eq. (3) reveals the important characteristics. It is a transfer function with unity gain and feedback equal to Z s (s)/Z l (s). According to the linear control theory, the H(s) is stable, if and if only, when ration Z s (s)/Z l (s) meets the requirement of Nyquist stability criterion [7].
In the above analysis for stability, it is assumed that the source is ideal voltage source and it remains stable under unloaded condition. However, the gridconnected inverters are usually current controlled. Hence above analysis is not much useful. So, the source should be represented as current source. To arrive at the equivalent current source, the same small signal voltage source is modified. The voltage across load is given by Similar to the above analysis, the current source is assumed to be stable under unloaded condition. The load is stable when connected to ideal current source. Under this condition, I(s) and 1/Y l (s) are stable. Under this condition, stability of V(s) depends on stability of second term of Eq. (5). This again resembles the closedloop transfer function with negative feedback. The gain is unity and the feedback factor is Y s (s)/Y l (s). Therefore the system is stable, if, along with above conditions, it meets the Nyquist criterion. In Eq. (5) admittances are used instead of impedance, though the analysis still can be carried out in terms of the impedances. In this case, Eq. (5) becomes It is important to note the requirement of stability. In the voltage source model, the output impedance of source should be as low as possible (ideally zero); whereas in the current source model, the output impedance should be as high as possible (ideally infinite).

Grid-connected inverters
The modelling of impedance of grid-connected VSC has very important use in analysis of stability and resonance phenomena when converter is integrated into the grid [5]. The grid-connected converter used in renewable sources is modelled as a current source in parallel with an impedance, i.e. Norton's equivalent circuit [6]. The stability of grid-connected inverter can be determined by Nyquist criterion [8]. The control structure of most of the VSCs is developed in the rotating dq reference frame [7]. The phase lock loop is used to synchronise converter with grid [9]. The use of complex multiple-loop control structure introduces non-linearities. These are generally overlooked in the simplified low-order modelling [10]. On the flip side, the detailed model introduces complexity and cross-couplings between various terms, which makes the determination of output impedance cumbersome. The trade-off way suggested in some literatures is to linearise the model by small signal analysis technique.
The impedance of converter-interconnected generator is affected by various factors such as control parameters, PLL, switching delays and converter harmonic filters. The converter is basically controlled by output current signal. In Figure 4, i 1 is the converter current and i 2 is the grid current. The converter is controlled either by i 1 or i 2 . If the grid current is the control variable, then the current control loop is where Y 21 is the forward trans-admittance of the filter, G PI s ð Þ is the proportional-integral-type current controller and G D s ð Þ is the switching delay. Here, the converter output voltage is considered as pure sinusoidal; if there is a noise in the voltage, then it needs to be considered as a disturbance signal. If, the converter is controlled by taking converter current i 1 , then Y 21 is replaced by output conductance Y 11 . The forward transconductance Y 21 is given by The simple transfer function of current controller is given by Higher-order controller also can be used for current control, but here simple PI controller is considered for the sake of simplicity. The delay, estimated by Pade's approximation technique, is given by The frequency response of output impedance is plotted here with the line impedance. The phase at the intersection of two curves gives the phase margin (PM).

Phase lock loop
The phase lock loop is required for PI-based controller. So, the structure of PLL should be thoroughly assessed for stability. Proper structure of PLL and control scheme helps to mitigate adverse impact on the system. The purpose of PLL is to synchronise converter with grid. Harmonics in the grid may penetrate to the converter through the PI controller. If control parameters are not properly chosen, this may produce harmonics through circular effect. Figure 5 shows the block diagram of PLL. The detail block diagram is given in Figure 6. Figure 6 is the traditional second-order generalised integrator-quadrature signal generator (SOGI-QSG) PLL. It can filter out higher-order harmonics, where u i is the input signal; u i 0 and qu i 0 are two output signals, which are in quadrature; k is the damping coefficient; and ω 0 is the output angular frequency of PLL. Eqs. (14) and (15) show the transfer function of PLL.
Simplifying above equation gives Similarly, for quadrature output transfer function is  In traditional SOGI-QSG PLL, the DC component in the input signal is not suppressed by the PLL. To overcome this problem, a minor modification is made in above PLL. The modified PLL is shown in Figure 7.
Transfer functions of modified PLL are given in Eqs. (16)-(18). This structure reduces the tracking error of PLL.
The bode plot of one of the PLL used in [12] is given here. The bandwidth of the PLL is 33 Hz. It attenuates harmonics of 1 kHz to À30 dB. However, the effect of PLL with other controllers and output filter needs to be investigated for crucial stability analysis (Figure 8).

Damped passive filter topologies
Different types of passive filters are used in converter-based renewable generation sources. The effectiveness of filter, particularly passive type, primarily depends on the grid strength and variation of grid impedance with time. LCL is the common type of filter used widely. The variation in grid impedance affects the performance  of LCL filter. Hence, the design of LCL filter is a trade-off between robustness and damping of resonance. The effective impedance with simple PI controller is explained here with RL filter topology.

R-L filter
Applying KVL at the converter output gives Converter output is controlled by output current With simple PI controller, the output becomes Putting Eq. (20) into Eq. (18) Equation (25) shows that there are two zeros and one pole. First, the output impedance decreases at À20 dB per decade up to first zero at α. At α, the impedance response becomes flat, and at β the impedance starts increasing at 20 dB per decade ( Figure 9). So, the response of integrated filter becomes similar to that of series resonance filter. The selection of α and β depends on the parameter selection of L fc , R fc , k P and k i . Bode plot for RL filter without and with PI current controller is given in Figures 10 and 11, respectively.

L-C filter
LC filter is used in converters for industrial applications like variable frequency drive (VFD) and uninterrupted power supply (UPS). It is simple in construction and relatively less costly. The analysis of LC filter with PI controller is given here ( Figure 12).
As per the standard practice, the value of L is selected such that its impedance at fundamental frequency should not drop by more than 3% of rated voltage. The capacitive reactance offers 1/5th of the fundamental inductive reactance at switching frequency of converter (around 3-4 kHz) to absorb harmonics effectively. Based on these criteria, LC filter is widely designed. The frequency response of LC filter with and without PI controller is given in Figures 13 and 14,   respectively. It is clear from the difference in bode plot that PI controller changes the frequency response.

L-C-L filter
Initially LC filter was used for converter applications, but grid-connected inverter has unique requirements that LC filter may not provide. Properly designed LCL filter may overcome the drawbacks of LC filter ( Figure 15).
Applying KVL,  The converter output voltage V c is a function of the grid current I g . Then Putting this in Eq. (26) and simplifying further give Figure 15. L-C-L filter.

Figure 16.
Frequency response of LCL filter with PI controller. Figure 16 shows the frequency response of LCL filter together with PI controller. It is clear from Eq. (29) and from Figure 16 that the PI controller increases the order of filter, so the frequency response of passive filter gets changed by the controller action.

Optimised filter design
The design criterion for filter design should comply with the regulatory requirement. As per IEEE 519-1992 standard, the current harmonics for weak grid condition (I sc /I L ) should be less than 0.3%. The ripple is caused by pulse width modulated signal. Output voltage varies from zero level to DC voltage level (V dc ). The modulated wave causes ripple in the current, which can be reduced by proper selection of output filter parameters. Typical L-C-L filter is used in most of the inverter. The L1-C-L2 filter has three unknowns. The selection of these parameters depends on various factors. Grid condition is one of them. The strength of the grid decides the effectiveness of filters. Typical grid impedance varies from 5 to 8% with X/R ratio in the range of 7-10 [13].
The switching frequency decides the ripple level and ripple frequency. Generally, inverter switching frequency remains in the range of 3-5 kHz. The ripple current is reduced either by increasing switching frequency or by using passive filter at the inverter/converter output. Higher switching frequency is selected to reduce the ripple current at the generation point, but it adversely affects the converter (IGBT) losses [14]. The second option is to use large inductor at the output of converter, but this not only incurs high cost but also increases the core losses. Typically 20% ripple current is expected in the output current. Keeping this in consideration, the inductor L1 is given by [7] The capacitor rating is selected such that the reactive power of capacitor is neither too high nor too low. Higher reactive power demands more power from converter, which causes more loss in reactor L1 and also more loss in converter switches. A lower value of capacitor will increase the inductor size. So, the capacitor is selected so that the reactive power should be in the range of 15-20% of the rated power.

Type I filter
In LCL filter, the damping can be achieved by simply adding series resistance in series with capacitor C. It is obvious that large value of damping resistance (Rd) gives large damping. But damping is effective around the resonance point only [15]. Above the resonance point, damping weakens. Also, the large value of resistance causes higher losses (Figure 17).
Impedance of Type I LCL filter with damping resistance is given by Frequency response of LCL filter with damping resistor is given in Figure 18. The grid parameters are Lg = 0.212 mH and Rg = 0.0095 Ohm, whereas the filter parameters are L 1 = 0.450 mH, L 2 = 0.300 mH and C = 270 uF. The damping resistance is varied from 0.05 to 0.8 Ohm. It is clear from the plot that the response of impedance is similar to inductive impedance. The notch is observed at the resonant frequency, which can be dampened by resistor in series with capacitor.  The gain margin is À40 dB and phase margin is around 120°. So, as per the Nyquist criterion, the filter response is very stable.

Type II filter
In Type II LCL filter (Figure 19), the damping can be achieved by simply adding a parallel combination of resistor and inductor in series with capacitor C. The effect of inductor is investigated using Nyquist plot. The value of R f = 10 Ohm and L f is varied from 50 to 500 μH. The effect of increase in L f observed using Nyquist plot is given in Figure 20.
Impedance of Type II LCL filter with damping inductor and resistance is given by Figure 19.
Nyquist plot of Type II LCL filter.
From the bode graph, it is observed that the gain margin is 70-100 dB for various values of L f . Similarly, the phase margin is 270°. As the value of L f is increased, there will be reduction in the gain margin.

Type III filter
In Type III LCL filter, the damping can be achieved by simply adding a parallel combination of resistor and inductor in series with capacitor C. The effect of inductor is investigated using Nyquist plot (Figure 21).
Impedance of Type III LCL filter with damping inductor and resistance is given by Impedance shows two resonance points: first is parallel resonance and second is series resonance (Figure 22). At first resonance point, the output impedance increases, and at series resistance it is at the minimum value. The phase margin of filter is around 210°. This ensures the stability of filter. The resonant frequency will shift with change in grid resistance, so the damping effect is difficult to predict.

Type IV filter
In Type IV LCL filter, the damping can be achieved by simply adding a parallel combination of resistor and capacitor Cf. The effect of resistor value is investigated using Nyquist plot (Figure 23).
Impedance of Type IV LCL filter with damping resistor R f in series with capacitor C f is given by The response of Type IV filter is similar to the tuned filter. The notch in frequency response of impedance is observed at the resonant frequency ( Figure 24). This notch can be damped by putting higher value of resistor in series with the capacitor C f . The phase of impedance sharply changes from À90 to 90°, which means the nature of impedance turns from capacitive to inductive. The ratio of C f /C decides the damping effectiveness. The larger the C f /C ratio, the larger will be the damping. Also, with increasing value of C f /C ratio, the loss in resistor R f is also increases as more and more current tends to flow in it [16][17][18].

Type V filter
In Type V LCL filter, the damping can be achieved by simply adding a parallel combination of resistor and inductor with capacitor C. The effect of inductor is investigated using Nyquist plot ( Figure 25).
Impedance of Type V LCL filter is given by Figure 24.
Nyquist plot of Type IV LCL filter.
Type V filter response shows two resonance frequency ( Figure 26). The first resonance creates high voltage distortion, while the second resonance point gives rise to current distortion. If converter generates harmonics equal to this resonance, then there will be high voltage distortion and current distortion. However, with lower value of resistance, damping can be achieved. An additional inductor L f can reduce the resistive loss.

Active damping of filter
Grid-connected converters may not function stably under the harmonic resonance condition. Harmonic resonance occurs when the converter impedance and grid impedance becomes equal in magnitude and 180°out of phase. The converter is generally connected to grid through filters, which is mostly LCL type. So, the parameters of LCL filters play an important role in keeping the successful functioning of converter. Also, the control structure of the converter shapes its output impedance. So, the control parameter should be selected such that it keeps the converter in the safe zone at all frequency. This is explained further here with analysis.  The inverter output characteristics depend on several factors like parameter of LCL filter, grid impedance, type of controller and control parameters. By selecting the appropriate parameter, the inverter can be operated with resistive output impedance, inductive output impedance and capacitive output impedance. Most of the inverters are operated with inductive (L-type) output impedance. This inverter is known as the L-type inverter. Here it is explained how an inverter impedance can be made capacitive (C-type inverter) with selection of virtual current control loop.
The output impedance of inverter is controlled by the virtual impedance. The single-line diagram of inverter with control loop is shown in Figures 27 and 28. The control is implemented with current loop and voltage loop, which gives good tracking behaviour and good output voltage. The feedback of current is taken from branch between the two capacitors. The first capacitor has a value of βC, and the second capacitor has a value of (1-β)C. Thus, the overall value of capacitor is C. The output impedance of system after adding virtual branch G V is given by The value of virtual impedance G V decides the overall impedance characteristics of output impedance. If G V is selected as per Eq. (38),  With the above value of G V , the output impedance becomes capacitive. So, by selection of proper control loop and control parameter, the shaping of inverter output impedance can be effectively done. Also, the desired value of damping can be achieved. This method is known as active damping method. Like this, the inverter can be made either L-type or R-type.

Conclusion
The characteristics and response of impedance of grid-connected inverter depend on various factors like selection of output filter and its parameter, controller type and its parameters, structure of PLL, the delay in switching of converter, etc. In this chapter, using first-order PI controller and different controller types, it is explained in a simple way. The complexity of impedance increases with the order of the controller and its structure. Also, the reshaping of inverter impedance can be done by selection of suitable control structure. By adding virtual impedance, the output impedance can be made inductive, capacitive or reactive without adding any additional hardware. In this work different variants of passive filter are explained, and how their output impedance behaves at different frequency is also explained with the use of Nyquist plot. The active damping method is explained, whereby damping is achieved without using resistor in the filter, which is considered as a very energy-efficient method of achieving damping. The output impedance can be reshaped either by external hardware in the form of filter or by adding virtual control loop in the controller. The stability of converter depends on how converter output impedance interacts with the grid impedance. Numerous efforts have been made and can be found in literatures to achieve the stability only by reshaping the impedance. However, the stability region can be expanded by reshaping of impedance up to a certain extent only. Only reshaping of impedance does not guarantee the converter stability with changing grid condition. Also, the converter is more likely to fall in an unstable region with weak grid condition, but reshaping of impedance helps to prevent converter-based generation from going unstable. By analysing the output impedance in the frequency domain, the controller parameters can be adjusted to reshape the output impedance. By increasing the proportional gain of the PI controller, the magnitude of the output impedance can be increased and can enhance the ability of harmonic rejection. By increasing the integral gain of the PI controller, the phase of the output impedance can be increased to improve the stability of the system [19,20] (Table 1).