Battery Energy Storage System Damper Design for a Microgrid with Wind Generators Participating in Frequency Regulation

: Ancillary frequency control schemes (e.g., droop control) are used in wind farms to improve frequency regulation in grids with substantial renewable energy penetration; however, droop controllers can have negative impacts on the damping of wind turbine torsional mode, thereby reducing the lifespan of the turbine gearbox. This paper presents a battery energy storage system (BESS) damper to improve the damping of torsional vibrations when using doubly fed induction generators (DFIGs) for frequency regulation in a microgrid. We formulated a linearized model comprising diesel generators, a wind turbine with five-mass drivetrain, and BESS. We also designed a feedforward compensator to deal with phase lag between the BESS damper signal and DFIG torque. The proposed BESS damper was shown to improve the torsional mode damping by moving the eigenvalues for torsional mode leftward to desirable locations on the complex plane. Dynamic simulations performed using MATLAB/SIMULINK demonstrated the efficacy of the proposed BESS torsional mode damping scheme in terms of torsional mode 1 damping performance and frequency response.


Background and Literature Review
Ancillary frequency control systems (e.g., droop control) are commonly used in wind farms to enhance the stability of grids with high renewable energy penetration.Unfortunately, droop controllers can have an adverse effect on wind turbine torsional mode [1].In an effort to improve torsional mode damping, researchers have designed resonance dampers for the rotor-side converters used in wind generators [1][2][3].
Ancillary frequency control in wind generation systems is generally performed using droop controllers or inertial controllers [4].The former generates a supplementary frequency regulation signal proportional to deviations in system frequency, whereas the latter generates a frequency regulation signal proportional to the rate at which the system frequency changes.In this paper, we will focus on droop controllers, since it was pointed out in [5] that the optimal inertia gain with the highest frequency nadir was close to zero.Droop gain must be selected carefully to ensure that the frequency nadir can be improved without imposing wind generator stalls or secondary frequency dips, due to a drop in the kinetic energy of the generator.In [6], a dynamic droop controller based on the rate of frequency change was proposed.In [7], droop gain was adjusted according to generator speed.Lee et al. [8] proposed an approach in which droop gain was varied according to the kinetic energy of the wind generator.Researchers have recently developed dynamic droop control methods based on fuzzy logic [9,10].In [11], a frequency control strategy based on coordinated robust dynamic droop power sharing was proposed for an islanded wind-powered microgrid.Researchers have recently developed piecewise linear time-variable droop functions [12,13].Kumar [14] proposed a coordinated droop and inertia control strategy for de-loaded offshore renewable power generators with HVDC-link.The impact of primary and secondary frequency control on the frequency variations caused by wind variability was investigated in [15].In [16], the authors presented a novel frequency control approach that associates the grid-side converter (GSC) with a synchronverter and a BESS.A critical review of synchronverter technology which could increase the inertial response of renewable energy sources (RESs) was presented in [17].
Subsynchronous oscillations have attracted considerable attention for wind generators connected to a weak grid using series capacitors for compensation [18].In [19], the induction generator effect and sub-synchronous control interaction were discussed for two different types of wind turbines: a doubly fed induction generator and a permanent magnet synchronous generator.In [20], the authors presented a novel approach to analyzing damping torque to detect and assess the risk of open-loop modal resonance exciting torsional sub-synchronous oscillations.Li et al. [21] proposed a high-pass filter with a proportional controller to mitigate subsynchronous oscillations.In [22], the authors developed a method that uses a bandpass filter and lead-lag compensator to improve the damping of subsynchronous oscillations.In [23,24], a compensation filter with damping terms was added to the rotor-side converter to mitigate subsynchronous oscillations.Researchers in [25] proposed a nonlinear sliding mode damping controller for subsynchronous oscillations.Two control laws based on stator voltage and current were developed in [26] to regulate active and reactive power with the aim of improving transient response damping.A linear quadratic regulator was proposed in [27] for the damping of subsynchronous interactions in doubly-fed induction generator-based wind farms.Ghaffarzdeh et al. [28] proposed an adaptive supplementary controller based on the multiple-model adaptive control approach to mitigate subsynchronous resonance induced by a type III wind system.Ali et al. [29] investigated sub-synchronous control interactions in doubly-fed induction generators based on eigenvalue analysis.In [30], the authors investigated the impact of wind park operating conditions and DFIG control system parameters on subsynchronous control interactions.
It has recently been reported that ancillary droop controllers can excite torsional vibrations in wind turbines [1].Most ancillary droop controllers take frequency deviation as the input and transmit an output signal to adjust torque commands for the rotor-side converter (RSC) of the DFIG.This implies that oscillations in system frequency imposed by load variations can excite DFIG torsional oscillations via the RSC torque command.In [1], researchers used Bode plots to design a torsional oscillation damper using a two-mass drivetrain model for wind turbines.Mandic et al. [2] designed three bandpass dampers for three torsional modes based on a five-mass drivetrain model.Note, however, that they did not account for the effect of droop control on torsional modes.In [3], researchers developed a sliding mode controller to enhance torsional mode damping.Licari et al. [31] proposed a model-based torsional damper using a Kalman filter.
The ever-increasing penetration of renewable energy into power systems has prompted considerable interest in the use of BESSs for frequency regulation.Teleke et al. [32] developed a model predictive control (MPC) technique based on receding horizon control to manage BESS power in accordance with forecasted wind conditions.Researchers have also developed a fuzzy-logic based frequency controller for wind farms that use a BESS [33].These methods have made it possible to minimize BESS capacity while largely eliminating the need for inflexible wind energy de-loading.Wu et al. [34] constructed a small-scale BESS to avoid the secondary frequency drops during rotor-speed recovery.Researchers have also proposed a proportional integral (PI)-lead and lead-lag controlled BESS to enable frequency, as well as voltage regulation for multimachine power systems [35].Choi et al. [36] proposed a frequency regulation scheme that uses a hybrid strategy involving wind generator kinetic energy control and BESS output power control.In [37], researchers developed a fuzzy-based wind-hybrid energy storage system to reduce frequency variations caused by fluctuations in wind power.In [38], researchers developed an artificial neural network (ANN) controller to regulate the power flow between wind generators and the grid.Shim et al. [39] proposed a droop control and state of charge (SoC) feedback (DaSOF) scheme for BESSs which aimed to improve frequency regulation.Zhao et al. [40] recently constructed a frequency-domain model to enable the pole-based dynamic analysis of control interactions between grid-forming BESSs and offshore wind power plants.
Researchers have also created control systems to regulate the power flow in and out of the BESS with the aim of enhancing the stability of wind-diesel power systems integrated with a BESS [41].Daraiseh [42] assessed the ability of energy storage systems to prevent frequency instabilities while providing a primary frequency response.In [43], Gu et al. developed a mathematical model for a vanadium redox flow battery (VRFB), as well as a scheme by which to smooth fluctuations in wind power.Researchers in [44] evaluated the benefits of demand response control and battery energy storage on frequency regulation in a power system integrated with a wind farm.In [45], researchers proposed a novel control strategy for a hybrid energy storage system comprising a BESS and a supercapacitor energy storage system (SCESS).Their system diverts the low-frequency components of power imbalance to the BESS, while the high-frequency components are diverted to SCESS.Parthasarathy et al. [46] recently proposed a detailed lithium ion battery model for the design of BESS controllers for active network management.In [47], researchers utilized a model predictive controller to enhance the dynamic performance of a system comprising a wind-driven synchronous generator, a fuel cell, and a BESS.In [48], a novel concept of a Smart Battery that brings together batteries with advanced power electronics and artificial intelligence (AI) was proposed for battery lifetime improvement.In [49], researchers reviewed the existing methods of increasing the energy efficiency of electric transport by analyzing and studying the methods of increasing energy storage capacities.In [50], researchers introduced the wear density function to provide a precise lifespan prediction of BESS.Furthermore, a quantitative tradeoff relationship between frequency regulation performance and battery lifecycle was derived using operational data of the actual BESS for frequency regulation.In [51], a frequency control method, in which battery energy storage systems (BESSs) participate in automatic frequency restoration reserve (aFRR) provision through their integration in the AGC of an island system, was presented.In [52], a BESS model with variable efficiency was used to compare the provision of Frequency Containment Reserve (FCR) with different SoC restoration strategies.
In view of widespread use of BESSs in a system with a high penetration of wind generators, we propose a design for a torsional damper on the BESS in this paper.

Contributions
The main contributions of this paper are as follows: (1) We present a complete model for a microgrid comprising diesel generators, wind generators with droop control, and a BESS.(2) We derived a transfer function that relates the BESS damping signal to the DFIG damping torque.We also designed a feedforward compensator for the BESS damper based on the phase lag computed using the transfer function.(3) To improve the damping ratio for the torsional mode, we use the pole assignment method to shift the eigenvalues of poorly damped modes leftward to desired locations.The eigenvalue location with minimal damper gain is chosen in order to minimize power output from the BESS.

Structure of Paper
The remainder of the paper is organized as follows.Section 2 outlines the proposed microgrid with diesel generators, wind generators, and a BESS.A block diagram of the linearized model is also presented.Section 3 presents the eigenvalues, oscillation frequencies, and damping ratios computed for the microgrid in torsional and frequency modes with and without a BESS.Section 4 presents the BESS damper we developed to compensate for the poor damping ratio in torsional mode 1, which is the primary concern in this work.We first derived the transfer function from BESS damper signals to DFIG torque commands for use in developing a feedforward compensator to deal with the phase lag between the BESS damper signal and DFIG torque command.The constants for the BESS damper are determined using the pole assignment method.Section 5 presents MATLAB/SIMULINK simulation results indicating the effectiveness of the proposed BESS damper.As shown in Figure 1, the proposed microgrid comprises five diesel generators (each unit rated at 450 kW), two wind generators (each unit rated at 750 kW), and a BESS (rated at 500 kW). Figure 2 presents the linearized model for the frequency control of the microgrid using the parameters outlined in the Appendix A.  The wind turbine generator is the five-mass drivetrain model described in [2,3].The mechanical torque from the wind turbine is obtained using the method outlined in [53], as follows:

System Model
where ρ is the air density, A indicates the area swept by the blades, V W is wind velocity,  R is the turbine speed, and C P is the power coefficient of the turbine, which can be written as follows: The tip speed ratio λ TSR is defined as where R  indicates the blade radius.
Incremental mechanical torque, ΔT R , in the linearized model (Figure 2) can be derived as follows: where In maximum power point tracking (MPPT) mode, the maximum power generated by the wind turbine generator P wind, opt is related to generator speed ω G in accordance with the following expression: where K opt is a constant.This means that the electromagnetic torque of the wind generator under MPPT operation, T G, opt * , must be set as follows: Incremental electromechanical torque (∆T G, opt

*
) in the linearized model (Figure 2) can be derived as follows: As shown in Figure 2, when the wind turbine generator participates in frequency regulation, an ancillary control signal, ∆u anc , is added to ∆T G, opt * to obtain the desired electromagnetic torque command for the wind generator, ∆T G * , as follows: where K PW is the droop gain of the wind generator.The time constants for the current regulator are far smaller than the wind generator inertia constant H G ; therefore, we assume that the actual electromagnetic torque ΔT G is assumed to be equal to the torque command ∆T G * in (15).
As indicated by the block diagram of the proposed BESS droop controller and BESS damper (Figure 2), the BESS real power command, ∆P BESS * , is obtained as follows: * () where ∆P damper is the damper power and -K BESS ∆f is the power from the BESS droop controller.Details pertaining to the design of the proposed BESS damper are discussed in Section 4. Figure 3 presents the equivalent circuit of the battery [32] and a block diagram of the real and reactive power controllers and d-axis and q-axis current regulators of the BESS [33].
The time constants of the BESS equivalent circuit and associated power controllers and current regulators are far smaller than system inertia constant H T ; therefore, the time delay between BESS real power command P BESS * and its actual real power output P BESS is negligible.Thus, we can assume that ∆P BESS * in Figure 2 is equal to ∆P BESS .Figure 2 also shows that system frequency ∆f is governed by the swing equation of the diesel generators [54,55], as follows: diesel sys wind sys BESS sys load

Modal Analysis of Microgrid with Wind Generator and BESS
The state equation for the linear system in Figure 2 can be written as follows: where ] T is the state vector and ∆d = [ΔP load ΔV W ] T is the disturbance vector.Note that ΔP damper , ΔP damper1 , ΔP damper2 , ΔP damper3 , and ΔP damper4 are the state variables for the BESS damper.The system eigenvalues are computed as the solution to the following equation:

Modal Analysis of Microgrid without BESS
Table 1 lists the eigenvalues for a system without BESS under wind speed of 11 m/s based on the system parameters listed in the Appendix A. Note that only the eigenvalues for the frequency mode and torsional modes 0, 1, 2, 3, and 4 are of major concern in the design of BESS dampers.  1 lists the values for one frequency mode and five torsional modes, identified as 0, 1, 2, 3, 4. Mode 0 refers to a mode without shaft twist characterized by a real eigenvalue.Mode 1 is of particular concern in this work due to the fact that the damping ratio for this mode is poor (0.0319) and the other three modes (modes 2, 3, and 4) can only be excited by wind fluctuations at the frequency specific to that mode [2].By contrast, mode 1 can be excited by a step load disturbance causing a frequency mode oscillation at 0.668 Hz, which is very close to that of mode 1 (2.921Hz), when DFIG participates in ancillary frequency control.

Modal Analysis of Microgrid with BESS but without BESS Damper
Tables 2 and 3 list the eigenvalues, frequencies, and damping ratios, respectively, for torsional mode 1 and frequency mode in a microgrid with BESS but without a BESS damper.As shown in Table 2, the BESS had an insignificant impact on the eigenvalue, frequency, and damping ratio of torsional mode 1.As shown in Table 3, increasing the droop gain for the BESS increased the frequency mode damping ratio and decreased the frequency.Based on the fact that the damping ratio in torsional mode 1 was poor (0.031 for K BESS = 20 ) while the damping ratio in frequency mode was satisfactory (0.751 for K BESS = 20), we developed a BESS damper to improve the damping ratio for torsional mode 1.

BESS Damper Design Based on Modal Control Theory
In the current study, we designed a BESS damper to improve the damping ratio for torsional mode 1.The block diagram in Figure 4 (derived from Figure 2) illustrates how the proposed BESS damper provides damping torque ΔT G that is in-phase with the speed of the generator ∆ω G .As shown in Figures 2 and 4, the BESS damper generates a damping signal ∆P damper with ∆ω G as its input.The transfer function of the proposed BESS damper is as follows: Note that the transfer function in (20) comprises H 1 (s) for a feedforward compensator, H 2 (s) for a washout filter (high-pass filter), and H 3 (s) for a bandpass filter.Details pertaining to the design of the feedforward compensator and bandpass filter are presented in Sections 4.1 and 4.2.
The output from the BESS, ∆P damper , is multiplied by a base conversion factor, , to obtain the damping signal for system base, ∆P damper, sys , where B BESS and B sys , respectively, indicate the BESS base and system base.Transferring ∆P damper, sys from the BESS damper to the microgrid triggers a corresponding change in system frequency, ∆f, as follows: , ( ) ( ) Simple block manipulation in Figure 4 results in the following: The resultant torque ∆T Gdamper provided by the BESS damper is derived as follows: ( ) ( )

Design of Feedforward Compensator
If we let s = jω, then the transfer function for the washout filter H 2 (s) in ( 20) can be expressed as Note that H 2 (jω) ≈ 1 when ω = ω r = ω 1 = 2π × 2.92 = 18.35 since ω ≫ ω CMP = 2 [2].Note also that H 2 (jω) ≈ 0 when ω ≈ 0. Thus, speed signals Δω G of very low frequency are filtered out and mode 1 oscillation signals (with a frequency of ω 1 = 18.35) are allowed to pass.In a similar manner, bandpass filter H 3 (s) in ( 20) can be approximated as If we let s = jω, then we obtain the following: The transfer function at the mode 1 frequency of ω 1 = 18.35 can be written as Note that the torque provided by the BESS damper, ∆T Gdamper , at the mode 1 frequency can be derived using Equations ( 20), (21), and ( 23) as follows: In order for the BESS damper to provide best damping effect for torsional mode 1, it is essential that torque ∆T Gdamper (jω 1 ) be pure damping torque in phase with wind generator speed ∆ω G at mode 1 frequency ω 1 .Since H 2 (jω) ≈ 1 and H 3 (jω) ≈ -g 2 (g 1 ω 1 ) ⁄ when ω = ω 1 , the feedforward compensator must compensate for the phase lag of H f (jω 1 ).Thus, we obtain the following: As shown in Figure 4, we developed a second-order lead-lag compensator for phase compensation.Let T Z = 1, such that the transfer function H 1 (jω 1 ) for the lead-lag network in Figure 4 can be expressed as follows: The desired constant T P for the lead-lag compensator can be computed using Equations (31) and (32), as follows:

Design of Bandpass Filter
Table 2 lists the eigenvalues for torsional mode 1, λ 1 = σ 1 + jω 1 = −0.57± j18.36 as computed using Equation ( 19) for a microgrid with BESS droop control (K BESS = 20) but without a BESS damper (g 2 = 0).The damping ratio for mode 1 is as follows [56]: This eigenvalue λ 1 is very close to the imaginary axis and the damping ratio is not satisfactory.The proposed constant damping ratio method involves moving the unsatisfactory eigenvalue (λ 1 ) to a pre-determined location with the aim of obtaining a constant damping ratio that is better than the original value of 0.031.
The real part of the desired mode 1 eigenvalue, σ 1 ' , is computed for a specific angular frequency ω 1 ' using Equation (34).Substituting the desired eigenvalue λ 1 ' = σ 1 ' + jω 1 ' into Equation ( 19) renders the following: By separating the real part and imaginary part, Equation ( 35) can be rearranged as The desired BESS damper constant can be obtained by solving the following two equations: Note that different eigenvalues yield different damper constants (g 1 and g 2 ) with a corresponding effect on system performance (e.g., frequency nadir after a step load change).
Figure 5 depicts the real part of eigenvalue σ 1 ' , BESS damper constants g 1 and g 2 , and frequency nadir as functions of mode 1 frequency ω 1 ' in the case where B BESS = 0.5 MW, K BESS = 20, and ζ 1 ' = 0.25.As shown in Figure 5c, a frequency range of 11.0≤ω 1 ' ≤17.1 appears to be appropriate, as evidenced by a small BESS damper gain.In the current study, we selected frequency ω 1 ' to yield the smallest damper gain g 2 in order to minimize BESS output power.Figure 5a also shows that varying mode 1 frequency ω 1 ' yielded a corresponding change in the real part of mode 1 eigenvalue σ 1 ' , while the system frequency nadir remained essentially unchanged.The flowchart in Figure 6 highlights the steps of the BESS damper design process.

Simulation Results
Simulations in MATLAB/SIMULINK were conducted to assess the effectiveness of the proposed BESS torsional damper in terms of dynamic system responses to load disturbance and variations wind speed, the results of which are listed in Figure 1.Note that the sampling time (time step) used in the simulation was 0.0001 s.

Effects of Wind Turbine Generator (WTG) Droop Control on Torsional Mode Oscillations: Microgrid without BESS Damper
As mentioned in Section 3.1 (Table 1), underdamped torsional mode 1 can be triggered by load disturbance in a DFIG with droop control.This is due to the fact that the oscillation frequency is close to the system frequency mode when using the system frequency as an input for the droop controller.The effect of droop control on mode 1 oscillation in a microgrid without a BESS damper is illustrated by the simulation results in Figure 7 showing dynamic response curves and FFT spectra as the grid is subjected to a 10% step load increase at t = 1 s.As shown in Figure 7a, the frequency nadir was improved from 59.4 Hz to 59.7 Hz when droop control was employed for the wind turbine generator.This improvement can be attributed to increases in wind generator torque and power, respectively, shown in Figure 7b,c.Note however that the use of droop control triggered mode 1 oscillations, as indicated by generator speed ω G and torsion θ 3G , respectively, shown in Figure 7d,e.The response curves in Figure 7a,f revealed that BESS droop control with a gain of K BESS = 20 increased the BESS output power, P BESS , from zero to 0.01~0.007p.u., as well as the frequency nadir from 59.7 Hz to 59.72 Hz.As shown in Figure 7d,e, mode 1 damping remained essentially unchanged when the BESS droop control gain was increased from K BESS = 0 to K BESS = 20.These results from time domain simulations are consistent with those from modal analysis in Table 2.
In Figure 8b, the two curves revealed an oscillation frequency of roughly 2.9 Hz in the Fast Fourier Transform (FFT) spectra of torsions θ 3G under the effects of wind turbine droop control.The FFT spectra in Figure 8b also revealed that step load disturbances triggered only mode 1 oscillations at a frequency of 2.9 Hz.As shown in Table 1, no excitation was observed in the other torsional modes, including mode 2 (292 Hz), mode 3 (372 Hz), and mode 4 (1974 Hz).In addition, our analysis of system frequency in the FFT spectra revealed that step load disturbance triggered frequency mode oscillations at roughly 0.67 Hz.The fact that step load changes triggered only frequency mode oscillations (0.67 Hz) and torsional mode 1 oscillations (2.9 Hz) (i.e., not torsional modes 2, 3, or 4) can be explained by the FFT spectrum in Figure 9b for the step input in Figure 9a.Here, the magnitudes in frequency mode (0.67 Hz) and torsional mode 1 (2.9 Hz) exceeded 300, whereas the magnitudes in mode 2 (292 Hz), mode 3 (372 Hz), and mode 4 (1974 Hz) were less than 10.This explains the triggering of frequency mode and torsional mode 1 oscillations by step load changes in Figures 7 and 8   The response curves in 10a revealed an increase in the frequency nadir with an increase in BESS droop gain.Increasing BESS droop gain decreased the torque and power of the wind generator (see Figure 10b,c), increased BESS power (see Figure 10f), and had only a small effect on damping for mode 1 (see Figure 10d,e).Taken together, we selected a droop gain of 20 as the optimal trade-off between frequency performance and battery SOH.

Effects of BESS Damper on Torsional Oscillation: Microgrid with a BESS Damper
To improve the damping for torsional mode 1 in Figure 7, the BESS damper designed in Section 4 was applied to the BESS.The constants used in the simulation were as follows: A comparison of response curves with and without the BESS damper in Figure 11bd revealed a notable improvement in damping for mode 1.These results met our goal of increasing the damping ratio for mode 1 from 0.031 to 0.25 (see Section 4). Figure 11a revealed a slight increase in the frequency nadir (from 59.65 Hz to 59.67 Hz) when the BESS damper was applied.We also observed a change in BESS power (P BESS  Figure 11e) following the application of damper signal (P damper in Figure 11f) to the BESS.
The design of BESS damper allows for the selection of damping ratios for torsional mode 1. Table 4 summarizes   As shown in Figure 11b-d, increasing the damping ratio ζ 1 from 0.15 to 0.25 produced a notable improvement in mode 1 damping.

Effects of BESS Capacity on Dynamic Performance
To examine the effect of BESS capacity on the dynamic performance of the microgrid, we designed BESS dampers for three BESS capacities (B BESS = 0.3 MW, B BESS = 0.5 MW, and B BESS = 1.0 MW) using the method outlined in Section 4, the results of which are listed in Table 5. Figure 12 presents the dynamic response curves of the microgrid subjected to a 10% step load change under the three BESS capacities in Table 5.As shown in Figure 12a, using a BESS with large capacity was shown to improve the frequency nadir.Figure 12d shows that better torsional damping can be achieved with larger BESS capacities.Increasing BESS capacity was also shown to increase BESS damper output  damper (Figure 12f) and BESS output power P BESS (Figure 12e).

Comparison with Wind Turbine Generator (WTG) Damper
The dynamic performance of a grid equipped with a BESS damper was compared with that of a grid using a wind turbine generator damper [2] in terms of response curves under the effects of a 10% step load change, the results of which are presented in Figures 13 and 14 and Table 6.As shown in Figure 14 and Table 6, the frequency nadir for the microgrid with a BESS damper (59.6726Hz for ζ 1 = 0.25) was higher than that of the microgrid with a WTG damper (59.6331Hz for ζ 1 = 0.25).In addition, the response curves obtained from the microgrid with a BESS damper (Figure 13b-d) were similar to those from the microgrid with a WTG damper (Figure 11b-d).As shown in Figure 15b-d, the BESS damper improved the damping ratio for torsional mode 1 under the effects of a step change in wind speed.As shown in Figure 15e, the power delivered to the system by the BESS increased with the damping ratio.

Conclusions
This paper presents a BESS damper with the aim of improving the damping ratio for the torsional mode oscillations in a microgrid comprising diesel generators and wind generators.We formulated a complete linearized model that includes diesel generators, a wind turbine with five-mass drivetrain, and a BESS in order to derive the transfer function between the BESS damper signal and DFIG torque.We developed a feedforward mechanism to compensate for the phase lag between the BESS damper signal and DFIG torque.We also designed a bandpass filter to shift the torsional mode eigenvalues leftward to points along a line with the aim of improving the damping ratio on the complex plane.The effectiveness of the proposed BESS damper was assessed in terms of dynamic responses in a microgrid subject to load changes and wind fluctuations using simulations implemented MATLAB/SIMULINK.This analysis led to the following conclusions: (1) Based on FFT spectra of step responses, we determined that only torsional mode 1 (2.92 Hz) was excited by step load changes in the microgrid.Thus, this study focused exclusively on mode 1 in the design of the BESS damper.(2) Providing suitable compensation for the phase lag between the BESS damping signal and DFIG torque made it possible for the BESS damper to provide a damping signal capable of generating damping torque in phase with DFIG speed.(3) The proposed BESS damper improved the damping ratio for torsional mode 1 from 0.0319 to 0.25.(4) When using the feedforward compensator to deal with phase lag, the damping characteristics of the BESS damper were similar to those obtained using a DFIG damper.Note, however, that the frequency nadir obtained using the BESS damper (59.6726Hz) exceeded that obtained using the DFIG damper (59.6331Hz), due to the fact that the damping power was from a BESS instead of a DFIG.(5) The proposed BESS damper improved torsional mode damping in situations involving changes in wind speed, as well as stepped changes in load.(6) In microgrids with low system inertia, frequency deviations are relatively large and the torsional mode damping ratio is relatively small.(7) The proposed feedforward compensator proved effective in dealing with the phase lag between BESS power output and DFIG torque.As a result, the proposed BESS damper provided essentially the same damping effects as those reported in previous works, in which the damper was installed on rotor-side converters (RSCs) while avoiding the degradation of frequency response associated with dampers on the RSC.(8) BESS damper constants have been designed based on a specific operating condition.
Future work will focus on the adaptation of BESS damper constants in accordance with system operating conditions, such as the number of diesel units and wind generators.

Figure 1
Figure1presents a one-line diagram of microgrid with the BESS examined in this study.

Figure 1 .
Figure 1.One-line diagram of microgrid with the BESS examined in this study.

Figure 2 .
Figure 2. Linearized model of frequency control system for microgrid with BESS.As shown in Figure2, the diesel generator is equipped with a proportional controller-K PD ∆f = − 1 R D∆f (droop controller) for primary frequency regulation and an integral con-

Figure 3 .
Figure 3. BESS equivalent circuit of the battery and block diagram of power controllers and current regulators.

Figure 4 .
Figure 4. Block diagram of BESS damper and system frequency response transfer function.

Figure 10 Figure 10 .
Figure 10 presents dynamic response curves of the microgrid subjected to a 10% step load change, illustrating the effects of BESS droop gain (K BESS = 20, 25, and 30) without a BESS damper.

Figure 14 .
Figure 14.Enlargement of frequency response curves in Figures 10a and 12a.

Figure 15 Figure 15 .
Figure15presents the dynamic responses of a microgrid subjected to a stepped decrease in wind speed.

Table 2 .
Torsional mode 1 eigenvalues for the microgrid with BESS but without a BESS damper (K pw = 20, V w = 11 m/s).

Table 3 .
Frequency mode eigenvalues for the microgrid with BESS but without a BESS damper (K pw = 20, V w = 11 m/s).

Table 6 .
Damping ratios and frequency nadirs of microgrid with a BESS damper vs. a WTG damper.