Analysis of Various Energy Storage Systems for Variable Speed Wind Turbines

Abstract

During the high penetration of wind power, wind turbines can affect power quality directly due to an unstable and intermittency source. Voltage fluctuations, harmonics, and voltage drops might be factors in this environment. Energy storage systems (ESSs) with variable speed wind turbines (VSWTs) as a permanent magnetic synchronous generator (PMSG) and a doubly fed induction generator (DFIG) could be a solution to improve the power quality from the “variability” of wind power. This chapter investigates the proposed system, which comprises a hybrid ESS for the VSWT. It analyzes the ability of various ESSs (B, SC, and EDLC) based on VSWTs with various ESSs for power quality in terms of average THD (%) specified in reference to IEEE std-519-1992 and IEC 61400-21-Ed.2.0. In addition, this chapter investigates the DFIG with hybrid energy storage systems (Li ion battery and super capacitor ESS) for the economic evaluation in terms of payback time. The simulation results have been verified by a power system computer-aided design/electromagnetic transients direct current (PSCAD/EMTDC) to demonstrate the system performance under different scenarios.

1. Introduction

The DOE’s 2009 Annual Energy Outlook projected energy accumulation in the USA. The report describes that the future growth of domestic electric energy will increase by 26% for electricity sector by 2030. This might not only require additional generation capacity of 259 GW but also forecast to increase as 1.7% of total installed capacity until 2030 [1, 2]. Variable speed wind turbine prefers to the conventional operation to extract the maximum power from the unstable wind generation. Since speed of wind turbine is nonconstant, the generator should be controlled by the power electronic circuit. Two variable speed wind turbines (VSWTs) on a principle of power electronics could be classified as DFIGs and PMSGs [3, 4].

DFIG-based wind turbine is a wound rotor induction generator. It consists of three-phase windings on the rotor and stator. The stator winding of the wind turbine is directly connected to the grid winding, while the rotor winding of the generator is fed by variable frequency bidirectional back-to-back PWM based on the voltage source converter, which consists of rotor side converter (RSC) and grid side converter (GSC). It has typically about 30% of nominal generator power [3, 4]. The synchronous machine has the ability to provide its own excitation on the rotor. Such excitation may be obtained by means of either a current carrying winding or permanent magnets (PMs) [5]. PM excitation avoids the field current supply or reactive power compensation facilities [2, 6, 7, 8]. PMSG is a permanent magnet synchronous machine with its stator windings connected to the grid through a frequency converter. Signal frequency is generated via the pulse width modulation (PWM) with the DC link of back-to-back voltage source converters (VSCs) consisting of machine side converter and grid side converter [9, 10].

There are several state-of-the-art technologies on wind power and energy storage to improve power quality and stability by optimal control of ESS [11, 12, 13, 14, 15, 16, 17], hybrid ESS [18, 19], and renewable forecasting modeling [20, 21]. It is widely popularized that an elaborately optimal control strategy for ESS is to smooth wind power fluctuations for the renewable energy [11, 12, 13, 14]. Two configurations of DFIG-BESS for internal and external controllers have improved the capability integrated to the grid [15]. In [16], an optimal control scheduling was considered on the variable smoothing time constant and charging/discharging of power limits, which is to mitigate the wind power fluctuations, while extending the battery life cycle of the BESS. A control strategy of the ESS in the wind farm focused an open-loop optimal control scheme is to incorporate the operating limits of BESS based on the forecasted wind condition [17]. The author [18] introduced the hybrid energy storage system (HESS) to overcome the fast PV-wind power generation fluctuations by smart scheduling, which has a statistical approach for the capacity distribution of the HESS. However, they did not demonstrate real-time simulation with few limited cases. The control and energy management of the hybrid ESS (battery and super capacitor) with DFIG has coordinated the power flows and load demand [19]. Prediction scheme of wind power generation connected to battery ESS developed by using numeric weather prediction model which required to the input as detailed wind information [20]. However, they did not smooth to power variation scheme in this chapter. Energy management system of flywheel ESS used to fuzzy logic for two optimized models as constraint condition and determine optimization objective to control of FESS [21].

2. Configuration of the hybrid ESS of the variable speed wind turbines

The configuration of the HYESS-DFIG is illustrated in Figure 1(a). DFIG is a variable speed wind turbine with a partial scale power converter in the rotor circuit. The main objective of the RSC controller is to regulate active power (Ps) and reactive power (Qs) of the stator side independently, while the GSC controller is to keep the constant DC voltage (V_dc) to adjust reactive power (Qg) of the GSC, which transfers from the grid side [22]. Figure 1(b) describes that PMSG which is a variable speed wind turbine with a direct-drive generator connected to the grid through a full-scale power converter [2, 6]. Three different ESS such as battery (B), super capacitor (SC), and electrical dual layer capacitor (EDLC) were introduced in Figure 1. Each ESS controller comprised bank, inductance and two-quadrant DC/DC converter connected to the DC link as illustrated in Figure 1. Two DC/DC converter design includes two insulated gate Bi-polar junction transistors (IGBTs) diode switches as S1 and S2. To compare two different VSWTs, that is, DFIG and PMSG, three-phase voltage sources are employed at 20 kV, 60 Hz, 0.04 H, and 2.5 Ω, and PSCAD parameter of the wind turbine given in Table 1.

Figure 1.

Block diagram of the hybrid energy storage system connected to the configuration variable speed wind turbines: (a) DFIG (left) and (b) PMSG (right).

Generated Rated MVA	2 [MVA]
Machine rated angular mechanical speed	125.667 [rad/s]
Rotor radius	44 [m]
Rotor area	2124 [m^2]
Air density	1.229 [kg/m^3]
Gear box efficiency	0.97 [P.U]
Gear ratio-machine/Turbine	60
Equation for power coefficient	Mode 2

Table 1.

PSCAD/EMTDC parameters of the wind turbine for DFIG, PMSG in the proposed systems [23].

2.1. Model of the wind turbines

Variable speed wind turbine could be mathematically described by the follows (1), (2), (3). A wind turbine extracts kinetic energy (P_wind(kin)) from the swept area (Area) of the blades. The power of the airflow ($P_{\mathit{\text{airflow}}})$ is given by [3, 24]:

$${\mathrm{P}}_{\mathrm{air}}=\frac{1}{2}\bullet \mathit{\rho }\bullet \text{Area}\bullet {\mathrm{v}}^3$$E1

where $\rho ,\mathit{\text{Area}},v$ is the symbol of air density (approximately 1.225 $\mathrm{kg}/{\mathrm{m}}^3$), swept area of rotor $\left(m^2\right)$, upwind free wind speed (m/s), respectively.

The wind power (${\mathrm{P}}_{\text{wind}\left(\mathrm{kin}\right)})$ transferred to the rotor of VSWT could be expressed by

$${\mathrm{P}}_{\text{wind}\left(\mathrm{kin}\right)}={\mathrm{C}}_{\mathrm{pt}}\left(\mathrm{\lambda },\mathrm{\beta }\right)\bullet \frac{1}{2}\bullet \mathrm{\rho }\bullet \text{Area}\bullet {\mathrm{v}}^3$$E2

where ${\mathrm{C}}_{\mathrm{pt}}\left(\mathrm{\lambda },\mathrm{\beta }\right)$ is a power coefficient, $\mathrm{\lambda }$ is a function of tip speed ratio $,$ and $\mathrm{\beta }$ is pitch angle. Wind turbines are characterized by their aerodynamic torque or turbine torque as Ta [3, 24, 25]:

$$\mathrm{Ta}=\frac{{\mathrm{C}}_{\mathrm{pt}}\left(\mathrm{\lambda },\mathrm{\beta }\right)\frac{1}{2}\bullet \mathrm{\rho }\bullet \text{Area}\bullet {\mathrm{v}}^3}{{\mathrm{\Omega }}_t}$$E3

where A is the maximum value of C_pt defined by the Betz limit, which a turbine can never extract more than 59.3% of the power from an air stream [21, 25]. In reality, wind turbine rotors have maximum C_p values in the range 25–45% [7].

2.2. Park’s model of the VSWT

For our proposed system, the Park’s model was used in VSWT (DFIG and PMSG) [3]. The voltage equations in the d, q axes-frame are given as follows [4, 5, 22, 26]:

$$V_d=\frac{d{\varphi }_d}{\mathit{dt}}-{\omega }_r{\varphi }_q+R_a{I}_d$$E4

$$V_q=\frac{d{\varphi }_q}{\mathit{dt}}+{\omega }_r{\varphi }_d+R_a{I}_q$$E5

$${\omega }_r=\frac{\mathit{d\theta }}{\mathit{dt}}$$E6

$$\theta ={\int }_0^t{\omega }_r\mathit{dt}+{\theta }_0$$E7

${\varphi }_d,{\varphi }_q,I_d,I_q,V_d,V_q$ are the stator fluxes, current and voltage of the d and q components; $R_a\text{is the stator resistance}$; ${\omega }_r$ is the rotor speed in electrical [radians/s]. The general Park’s model is used to explain about the control scheme of an induction machine. The equation of voltage vector in the DFIGs is expressed as follows [4, 5, 22, 26]:

$$V_s=R_s{i}_s+\frac{\mathit{d\varphi }}{\mathit{dt}}+j{\omega }_s{\varnothing }_s$$E8

$$V_r=R_r{i}_r+\frac{\mathit{d\varphi }}{\mathit{dt}}+j{\omega }_r{\varnothing }_r$$E9

$${\varnothing }_{\mathrm{s}}={\mathrm{L}}_{\mathrm{s}}{\mathrm{i}}_{\mathrm{s}}+{\mathrm{L}}_{\mathrm{m}}{\mathrm{i}}_{\mathrm{r}}$$E10

$${\varnothing }_r=L_m{i}_s+L_r{i}_r$$E11

The equivalent two-phase model of the symmetrical variable speed wind turbine (dq frame) used to synchronously rotating reference frame (dq frame) as follows [7]:

$$v_{\mathit{sd}}=R_s\bullet i_{\mathit{sd}}+\frac{d{\varphi }_{\mathit{sd}}}{\mathit{dt}}-{\omega }_s{\varphi }_{\mathit{sq}}$$E12

$$v_{\mathit{sq}}=R_s\bullet i_{\mathit{sq}}+\frac{d{\varphi }_{\mathit{sq}}}{\mathit{dt}}+{\omega }_s{\varphi }_{\mathit{sd}}$$E13

$$v_{\mathit{rd}}=R_r\bullet i_{\mathit{rd}}+\frac{d{\varphi }_{\mathit{rd}}}{\mathit{dt}}-{\omega }_r{\varphi }_{\mathit{rq}}$$E14

$$v_{\mathit{rq}}=R_r\bullet i_{\mathit{rq}}+\frac{d{\varphi }_{\mathit{rq}}}{\mathit{dt}}+{\omega }_r{\varphi }_{\mathit{rq}}$$E15

The stator and rotor fluxes using the synchronously rotating reference frame (d,q frame) are given as follows [5, 8, 9, 22]:

$$\left\{\begin{array}{c}{\mathrm{\varphi }}_{\mathrm{sd}}={\mathrm{L}}_{\mathrm{s}}{\mathrm{i}}_{\mathrm{sd}}+{\mathrm{L}}_{\mathrm{m}}{\mathrm{i}}_{\mathrm{rd}}\\ {\mathrm{\varphi }}_{\mathrm{sq}}={\mathrm{L}}_{\mathrm{s}}{\mathrm{i}}_{\mathrm{sq}}+{\mathrm{L}}_{\mathrm{m}}{\mathrm{i}}_{\mathrm{rq}}\end{array}\right.$$E16

$$\left\{\begin{array}{c}{\varphi }_{\mathit{rd}}=L_r{i}_{\mathit{rd}}+L_m{i}_{\mathit{sd}}\\ {\varphi }_{\mathit{rq}}=L_r{i}_{\mathit{rq}}+L_m{i}_{\mathit{sq}}\end{array}\right.$$E17

where $R_s,R_r$ indicate equivalent resistance of stator and rotor windings. $L_s,L_r,L_m$ indicate self and mutual inductances of stator and rotor windings, respectively. Ignoring the power losses in the stator and rotor resistances, the active and reactive powers from the stator are given by

$$P_S=\frac{3}{2}\left(v_{\mathit{sd}}\bullet i_{\mathit{sd}}+v_{\mathit{sq}}{i}_{\mathit{sq}}\right)$$E18

$$Q_S=\frac{3}{2}\left(v_{\mathit{sd}}\bullet i_{\mathit{sd}}-v_{\mathit{sq}}{i}_{\mathit{sq}}\right)$$E19

$$P_r=\frac{3}{2}\left(v_{\mathit{rd}}\bullet i_{\mathit{rd}}+v_{\mathit{rq}}{i}_{\mathit{rq}}\right)$$E20

$$Q_r=\frac{3}{2}\left(v_{\mathit{rd}}\bullet i_{\mathit{rd}}+v_{\mathit{rq}}{i}_{\mathit{rq}}\right)$$E21

The electromagnetic torque is given by [4, 5, 22, 26]:

$$T_e=\frac{3}{2}p{\varphi }_s{i}_{\mathit{sq}}=-\frac{3}{2}p{\varphi }_s\frac{L_m}{L_s}{i}_{\mathit{rq}}=-\frac{3}{2}p\frac{V_{\mathit{sq}}}{{\omega }_s}\frac{L_m}{L_s}{i}_{\mathit{rq}}$$E22

The synchronous model is also expressed in the (d,q) synchronous Park’s model, and the voltage equations of the PMSG are represented as follows [7]:

$$\frac{d{i}_{\mathit{sd}}}{\mathit{dt}}=\frac{1}{L_d}\left(-R_s{i}_{\mathit{sd}}+\omega L_q{i}_{\mathit{sq}}-V_{\mathit{sd}}\right)$$E23

$$\frac{d{i}_{\mathit{sq}}}{\mathit{dt}}=\frac{1}{L_q}\left(-R_s{i}_{\mathit{sq}}-{\omega }_e{L}_d{i}_{\mathit{sd}}-V_{\mathit{sq}}+{\omega }_{\mathit{PM}}{\varnothing }_f\right)$$E24

where V_sd, V_sq, i_sd and i_sq are voltage and current for the d,q axis of the stator side. ∅_f is the magnitude of the flux linkages by using the permanent magnetic flux [Web]; ω_PM is the value as rotating speed [rad/s] of the PMSG. The electromagnetic (EM) torque in the rotor could be expressed as follows [7]:

$$T_e=\frac{3}{2}p{\varnothing }_f{i}_{\mathit{sq}}$$E25

p is the number of pole pairs in the PMSG.

3. Modeling and control of the grid connected to various ESSs

Energy storage systems (ESSs) have a function of converting electrical energy from a power system network into a form that can be stored for converting back to electrical energy when needed [8, 9]. ESS has numerous applications including portable devices, smart grid, building integration, energy efficiency, transport vehicles, and stationary renewable energy resources [27]. In this chapter, only three different distributed ESSs for renewable generation systems were introduced as battery, super capacitor, and electrical dual layer capacitor [23]. Table 2 summarizes the characteristics of three different ESSs.

Table 2.

Properties of three different ESSs [23].

3.1. Modeling and control of energy storage systems

Figure 2 shows that each E-ES, for this study, consists of an energy source bank, an inductance and a two-quadrant DC/DC converter connected to the DC link. It also describes how the controller of the DC/DC buck-boost mode generates the gate signals for gate 1 and gate 2 [23].

Figure 2.

Block diagram of the ESS in PSCAD/EMTDC: modeling of ESS (upper) [14] and controller of ESS (lower) [29].

The reference value of the i_ess can be made by using a voltage feedback control between V_dc-ref and V_dc with PI controller of the voltage. It can be expressed as follows [2, 25, 28]:

$$i_{\mathit{ess}₋\mathit{ref}}=\left[k_{\mathit{vp}}+\frac{k_{\mathit{vi}}}{s}\right]\left(v_{\mathit{dc}₋\mathit{ref}}-v_{\mathit{dc}}\right)$$E26

The inner current control loop using a PI controller of current can be produced by using the duty ratio to generate the gate-signal $i_{\mathit{ess}₋\mathit{ref}}$, as given by Eq. (27) [2, 25, 28]:

$$D=\left[k_{\mathit{bp}}+\frac{k_{\mathit{bi}}}{s}\right]\left(i_{\mathit{ess}₋\mathit{ref}}-i_{\mathit{ess}}\right)$$E27

Finally, the gate signal is generated by comparing the difference between the duty ratio and the carrier frequency.

3.2. Energy flow and Total harmonic distortion of various energy storage systems

The HYESS which is applied to the VSWT system uses DC-DC converters with controller of the ESS illustrated in Figure 2. When active power of DFIGs and PMSGs in the rotor side (Pr) or stator side (Ps) is greater than active power of DFIG and PMSG on the grid side (Pg), the ESS charges the energy flow from the DC bus to the ESS Bank through the S1 switch and S2 diode [26, 29, 30].

Therefore, the source of ESS works to absorb active power from the DC voltage, while acts as a step-down converter when V_ess boosts. The DC-DC converter design used as a buck converter circuit, and then D1 (duty ratio) of S1 in the buck mode can be computed by (28) [31]:

$$\text{D1}=\left[{\mathrm{V}}_{\mathrm{ess}}/{\mathrm{V}}_{\mathrm{dc}}\right]$$E28

When Pg is bigger than Pr, the E-ES energy discharges through S1 and S2, and energy flows to the DC bus (V_dc) [27, 28, 29]. In this case, the converter acts as a boost converter mode [27, 28, 29]. The ESS bank serves as a source to supply active power, which results in the decrease of the voltage V_ess. The duty ratio D2 of S2 in the boost mode can be expressed as (29) [25]:

$${\mathrm{D}}_2=1-{\mathrm{D}}_1$$E29

Table 3 summarizes the energy flow in the two modes, that is, buck and boost modes.

Energy flow	Charge (Pr or Ps > Pg): From DC bus to the ESS	Discharge (Pr or Ps < Pg): From ESS to DC Bus
Bank	Sink as absorb active power	Source as supply active power
Vess	Increase	Reduce
Switch direction of two quadratic converters	S1 switch and S2 diode	S2 switch and S2 Diode
Converter act as	Buck converter as Vdc→reduce	Boost converter as Vdc→increase

Table 3.

Energy flow of the energy storage systems.

The harmonic distortion of the voltage and current waveforms is generally expressed in terms of the fundamental frequency [32]. Power injection from the HY ESS affects the power quality, which is described as voltage and current THD. In Ref. [30], the practices and requirements for harmonic control in electrical power systems are established. Figure 3 shows the PSCAD/EMTDC modeling for the voltage (Vg), current (Ig), and active power (Pg) of the THD measurement at the point of common connection (PCC) in the proposed system. With the increasing use of the nonlinear devices, harmonic distortion of the voltage waveform is a problem which is receiving the considerable attention. Any periodic waveform of nonsinusoidal form can be synthesized by expressing it as the sum of a series of harmonics of the fundamental frequency by using Fourier analysis (30) [31]:

$$\mathrm{f}\left(t\right)=a_0+a_{1}cos\mathit{\omega t}+a_{2}cos2\mathit{\omega t}+a_{3}cos3\mathit{\omega t}+\dots +b_{1}cos\mathit{\omega t}+b_{2}cos2\mathit{\omega t}+b_{3}cos3\mathit{\omega t}+\dots$$E30

where

$$a_0=\frac{1}{T}{\int }_0^{T}f\left(t\right)\mathit{dt},a_n=\frac{2}{T}{\int }_0^{T}f\left(t\right)cosn\bullet \mathit{\omega t}\bullet \mathit{dt},b_n=\frac{2}{T}\bullet {\int }_0^{T}f\left(t\right)sinn\bullet \mathit{\omega t}\bullet \mathit{dt}$$

Figure 3.

PSCAD modeling of 7th order 60 Hz for THD of active power, current, voltage at the PCC.

Total harmonic distortion (THD) is usually expressed as a percentage of fundamental voltage by the expression as (31):

$$\sqrt{\sum _{n=2}^{n=7}{\left(\frac{V_n}{V_1}\right)}^2\times 100\%}$$E31

where V1 = fundamental frequency voltage component, Vn = nth harmonic voltage component.

3.3. Battery energy storage systems (B-ESS)

Batteries have been widely used in many fields, which are referred to as electrical energy storage system, which accumulate electric energy in electrochemical form and delivers direct current (DC). The promising battery which is used as energy storage devices is lithium ion (Li-ion) battery. Due to high electrical potential and energy density of battery, Li-ion batteries are one of the promising solutions for the storage compared to other battery options. In addition, Li-ion batteries do not have poisonous metals (lead, mercury, or cadmium) and memory effect. The main disadvantage of Li-ion battery is the required high production cost [33]. The energy of battery could be calculated by (32):

$$E=\mathit{IVh}$$E32

where I (A) and V(V) are current and voltage of battery and h (hours) is the charging time. For this study, a Panasonic LJ-SK84A Li-ion storage battery system [34] is used to design about the modeling for the PSCAD battery as 8 KWh from 110 to 165 V in Figure 4.

Figure 4.

Modeling of PSCAD/EMTDC: BESS.

3.3.1. State of charge

In the Li-ion battery, over charging or over discharging could result in reducing the battery stack life and indirectly increasing the cost. However, accurate calculation of the state of charge (SOC: %) for battery ESS is essential for the battery outputs in powered system that aims at maximizing energy storage system’s performance, extending battery life, and realizing the safe operation of the systems including EV, renewable generation, and building integration.

The SOC of the battery is defined as the ratio of the remaining capacity ($Q_{\mathit{\text{present}}}$) to the nominal capacity of the cell ($Q_{\mathit{\text{nominal}}})$. It could be described as (33)

$$\mathrm{SOC}\left(\%\right)=\frac{Q{₋}_{\mathit{\text{present}}}}{Q₋\mathit{\text{nominal}}}$$E33

Q is usually measured in the unit of Ah or mAh.

The battery state of charge (SOC) estimation could be calculated by the Coulomb counting (CC) method through the integration of measured battery current or open circuit voltage (OCV) [35, 36, 37]. However, unknown initial SOC which caused by the measurement errors and noise in real practice can affect performance, which means unreliable SOC value results in reduced performance and potential risk to the battery system. Therefore, it is critical to develop algorithms that can estimate accurate SOC to improve the performance for the battery energy storage. There are several advanced methods which can be closer to true values over time to compensate for nonlinearity such as extended Kalman filter (EKF) and sigma-point (unscented) Kalman filter (SP/UKF) [35, 36, 37].

3.4. Super capacitor ESS

Figure 5 describes about the super capacitors (SCs), which have higher power density, higher round trip efficiency, longer life cycle life, and lower capital cost per cycle than batteries [9, 37]. Therefore, SC is a good candidate for short-term storage (i.e., seconds to minutes). Assuming the initial voltage across the super capacitor (C_sc) is V_sc (0) after charging for a period t, the instantaneous voltage across the capacitor V_sc is given by (34)

$${\mathrm{V}}_{\mathrm{sc}}\left(\mathrm{t}\right)={\mathrm{V}}_{\mathrm{sc}}\left(0\right)+\frac{I_L\bullet t}{C_{\mathit{sc}}}$$E34

where $I_L$ is the load current entering the capacitor from the unregulated power supply.

Figure 5.

PSCAD modeling of the super capacitor as 2F and EDLC [9, 39].

3.5. Electrical dual layer capacitor ESS

Capacitors composed of two conducting plates, which are separated by an insulating material [8, 29]. Conventional capacitors store little energy due to the limited charge storage areas and the separation distance between the two charged areas and the separation distance between the two charged plates [29, 38]. However, super capacitors based on the EDL mechanism can store significantly more energy because of the large interface area and the atomic range of charge separation distances [8]. A simple resistive capacitive equivalent circuit of the EDLC-cell is designed in PSCAD/EMTDC as shown in Figure 5.

4. Simulation results and discussion

Case I—Analysis of the Variable Speed Wind Turbines with Various ESSs for Power Quality in terms of THD (%).

The simulation results in Case I verify the performance of the hybrid ESS-VSWT by using the PSCAD design tool in terms of the THD (%). The mean wind speed is to make realistic results by using the PSCAD/EMTDC with random noise components during 60 s in Figure 6.

Figure 6.

PSCAD modeling and noise component for the mean wind speed (10 m/s).

Power quality issue could be verified by the THD which was monitored at the point of common coupling (PCC) in the VSWT-HYESS as illustrated in Figure 1. The data (THD: %) in Table 4 were extracted from using excel file and then taking an average value (THD: %) during 1 min.

Comparison of THD (%) with the HY ESS-VSWT (2 MW, 60 Hz, 10 m/s with random noise) during 60 s
	DFIG			PMSG
	Avg THD for P	Avg THD for I	Avg THD for V	Avg THD for P	Avg THD for I	Avg THD for V
No ESS	3.61	6.90	9.73	14.12	4.38	7.37
B- ESS	2.90	4.44	4.58	11.31	2.78	4.94
SC-ESS	2.77	5.68	4.65	10.97	2.30	3.73
E- ESS	2.83	5.65	4.7	11.14	2.78	5.00
B-SC ESS	2.90	5.74	4.43	11.18	1.82	4.96
B- E ESS	2.89	5.75	4.69	11.14	2.79	4.93
SC-E ESS	3.02	5.73	4.60	11.23	2.77	3.71
B-SC-E ESS	3.05	6.06	4.60	1.42	1.82	3.71

Table 4.

Comparison of THD at the PCC.

4.1. The case of the DFIG-HYESS

Figure 7 describes average THD (%) of a DFIG at the PCC for P, I, and V during 60 s. Table 5 shows that average THD (%) of the DFIG is a greater than other DFIG-HYESS. This result showed that the proposed HY ESS-DFIG could improve the power quality by reducing the THD. SC-ESS (2.77%), B-ESS (4.4%), and B-SC ESS (4.43%) are the best choices to enhance power quality in the system, as shown in Figure 8.

Figure 7.

THD (%) of only DFIG during 60 second: active power THD (Pthd: upper), current THD (ithd: middle), voltage THD (Vthd: lower).

Item	Cost ($)
Overnight cost of DFIG	2949000.00
Pb battery ESS	16862.66
Li-ion battery ESS	11992.85
Super capacitor ESS	5574125.98

Table 5.

Cost of 1.5 MW Wind Turbine [9].

Figure 8.

THD (%) of HY ESS DFIG during 60 secs: active power THD of SCESS, current THD of BESS, voltage THD of B-SC-E ESS.

4.2. In the HYESS-PMSG

Table 5 shows that the average THD of PMSGs is bigger than HY ESS-PMSG. These results indicate that an ESS which reduces THD (%) is to be improved. SC-E-B ESS is the most effective option in reducing THD for the voltage (3.71%), current (1.82%), and active power (1.42%), as shown in Figures 9 and 10.

Figure 9.

THD (%) of only PMSG during 60 s: active power THD (Pthd: upper), current THD (ithd: middle), and voltage THD (Vthd: lower).

Figure 10.

THD (%) of HYESS PMSG during 60 sec: active power THD of B-SC-E ESS, current THD of B-SC-E ESS, voltage THD of B-SC-E ESS.

Case II– Financial analysis of the DFIG with two ESSs in terms of the payback time ($)

Wind energy is free. However, financial analysis can be used to assess wind project investments such as how to screen cost and the benefits of the project. Regarding the financial analysis, the payback time depends on the amount of time to evaluate the performance of different ESSs such as super capacitor ESS (SCESS), battery ESS (BESS). Payback time is the time in which investors can recover all their investment. It depends on several factors such as O&M, inflation, and depreciation. It will show whether a project is worth establishing the investment to escape the risk. Shorter payback time indicates a more economical project. A 1.5 MW wind generator, which sells the electricity generated to utilities, can serve as an example. There are several manufacturers who provide large capacity DFIG. Generally, the overnight cost of a wind generator is $1966/kW. The Cost of O&M per year was assumed as $30.98/kW-Year for this work to find payback time from the [39]. A 83.33 kWh super capacitor ESS and Li Ion battery are chosen [40, 41]. The discount rate is also assumed as 9%. The capital cost of a wind generator is provided in Table 5.

Payback time in terms of accumulated cash flow could be expressed as (35)

$${\left(\frac{1}{1+\mathit{\text{Discount\hspace{0.17em}Rate}}\left(0.9\%\right)}\right)}^{\mathit{\text{year}}}=\text{Capital\hspace{0.17em}Cost}+\mathrm{Net}\text{Annual Saving}\left(\text{\$}1,032,838.80\right)$$E35

where capital cost is total overnight cost of DFIG ($) + cost of Li ion or SC ESS. The payback time for DFIG with Li ion ESS is 3.47 years; however, the payback time for DFIG with super-capacitor is 15.76 year as shown in Figure 11. If a hybrid ESS (B-SC) ESS is chosen, the payback time will be between 3.47 and 15.76 years, according to the proportion of BESS and SCESS given in this system.

Figure 11.

Accumulated cash flow.

5. Conclusion

This section analyzes the ability of various HY ESSs based on VSWTs to improve power quality in terms of THD under the same situations: three-phase voltage sources used at 20 kV, 60 Hz, 0.04 H, and 2.5 $\mathrm{\Omega }$. From the previous results, HYESS-VSWT has been verified to increase the power quality by reducing THD. SC-ESS (2.77%), B-ESS (4.4%), and B-SC ESS (4.43%) of the DFIG’s case would be a better option to improve power quality for active power, voltage, and current, respectively, compared to other cases. In the case of the PMSGs, SC-E-B ESS indicated as the best option to reduce THD for voltage (3.71%), current (1.82%), and active power (1.42%). Compared to average active power THD of a DFIG as 3.61% and PMSG as 14.42% without ESS, a THD (%) for PMSG has lower power quality difference output as 10.81% compared to THD (%) of DFIG. This might be a converter configuration issue because DFIGs are partially controlled by the back-to-back converter such as RSC and GSC, while PMSGs are directly connected through the back-to-back converter to the grid, which make 100% power rated. Regarding average THD, the PMSG-B E SC-ESS as 1.4% has the most outstanding output results compared to other options. Furthermore, the best option’s values from the simulation results satisfy the IEEE Standard requirements, which should be less than 4% for current THD and 5% for voltage THD [30].

In addition, a DFIG with various ESSs is a financially attractive option if a chemical battery ESS is used. The payback time is acceptable, and it still has a benefit for the environment. If the new energy storage technology as a super capacitor is considered, the payback time right now is relatively long. However, if there is a technology breakthrough, the cost of a super capacitor will probably decrease and finally be financially attractive. Presently, a hybrid system will still be relatively affordable. This financial forecast helps achieve the high-power quality and reasonable payback time during the short period.