Wind Power Integration in Hybrid Power System ...

3 downloads 80683 Views 262KB Size Report
properties of available technologies can be spread on ..... are then stored into classes (usually 20 classes are used) of .... Automotive Technology, Vol. 9, No.
March 26-29

Wind Power Integration in Hybrid Power System Active Energy Management M. A. Tankari, M.B.Camara, B. Dakyo, C. Nichita GREAH Laboratory, University of Le Havre, Le Havre, France, 25 rue Philippe Lebon, BP 540, 76058 LE HAVRE Cedex France E-mail: [email protected] E-mail: [email protected] Abstract: The wind speed is assumed to be a sum of slow and turbulence components containing the total energy. The caption of such large frequency band of energy will increase the profitability of a hybrid wind-diesel power source. Indeed, to absorb all fluctuations of the wind generator induced current, it is necessary to oversize diesel generator, or to connect storage devices for medium and high frequencies components absorption. Inserting a flywheel to absorb the high frequencies would increase significantly the battery lifetime. An original strategy presented in this paper is based on the frequency share of diesel and storage behaviour ensuring a better lifetime of these sources and meeting the demand of the load. The effectiveness of this strategy lies in its ability to control the frequency bands allocated to the diesel engine and to extend the batteries lifetime by individual charge/discharge cycles moderation. The aim of this paper is to analyze the degree of selectivity and the configurations of filters through a specific study of the spectral density of signals. The batteries lifetime is then evaluated using the rainflow cycle counting algorithm for effectiveness comparison. Keywords: Hybrid energy system, Energy storage, battery, wind turbine, flywheel, power spectral density, diesel engine. 1. Introduction 2. Wind speed model Wind diesel hybrid system is subject to disturbances due to wind speed fluctuations. Long life and low capital cost of diesel engine is guaranteed when storages devices are used. According to energy storage issue, the properties of available technologies can be spread on dual energy-power definition related to time horizon of effectiveness. This paper deals an energy transfer management strategy in hybrid system including wind turbine, diesel engine, flywheel, battery, capacitor, and emulated load management. In first case, the importance of taking into account of wind speed high fluctuations is introduced. In second case, the hybrid system management is presented through first and second order low pass filters studies. Finally, some simulations results are presented and analyzed.

1

The wind speed spectrum shown in Fig. 1 was constructed by Van Der Hoven (1957) [1] from long- and mean-term records at Brookhaven, New York. Few days’ variations (synoptic), variations with the time of day (diurnal) and shorter time-scales of minutes down to seconds or less wind-speed variations (turbulence) are observed. The phenomenon called ‘spectral gap’ occurs between the diurnal and turbulent peaks. It shows that the synoptic and diurnal variations can be treated as quite distinct from the higherfrequency fluctuations of turbulence. There is very little energy in the spectrum in the region between 2 h and 10 min. These turbulent fluctuations then have a zero mean when averaged over about 10 min to 2h. . This phenomenon is not observed in data analysed in [2] that show a smooth behaviour through de concerned region. The ‘spectral gap’ may

occur when measures at various altitudes or different sampling periods are combined [2]. One of the important information from Fig. 1, which is correlated by [2], is that high frequency fluctuations contain considerable power and it can have a very significant effect on the design and performance of the individual wind turbines, as well as on the quality of power delivered to the network and to the consumers. In [2], acquired data of 104 turbines located at 4 sites were studied. It is shown that wind speed frequency bandwidth is from 2.5*10-7 Hz to 1.4*10-4 Hz when measurements are carried out with a sampling period of 1 hour on 366 days, while it lies between 9.2 * 10-9 Hz and 0.5 Hz for a sample of 1s on 10 days. The significant part of the generator power fluctuations is located in the low (below 0.01 Hz) and medium (between 0.01 and 1 Hz) frequency regions. Most of the high frequency fluctuations (above 1 Hz) are effectively damped out by the large inertia of turbine generator; their magnitudes are insignificant [3].

are required to match medium and high fluctuations. The lifetime of the battery would be multiplied by four when a flywheel is used to absorb fluctuations in high frequencies instead of using only the battery [6]. The hybrid system presented in Fig. 2 is DC bus linked.

Fig. 2. Hybrid system with flywheel

The energy transfer between different sources and load is based on filtering wind power to generate the current reference of storage devices. The batteries allocated frequency band is 1.3mHz to 20 mHz. The flywheel absorbs high frequency currents (over 20mHz) and the diesel engine is used to regulate the DC bus voltage continuously (below 1.3mHz). 4. Filtring Principle

Fig. 1. Wind Spectrum Farm Brookhaven Based on Work by van der Hoven (1957)

So, the wind speed v(t) can be expressed as the sum of the harmonics characterised by the magnitudes Ai, the pulsation wi and the phase ϕi generated randomly [4][5]. v(t) = vl (t) + vt (t) =

2

Nl

2

N

∑ A cos(wt +ϕ ) + π ∑ A cos(wt +ϕ ) π i =0

i

i

i

i

i

i

Nt

(1) Nl are Samples for the slow component vl(t) and N-Nt are Samples for the component of turbulence vt(t). 3. Hybrid system management The energy capacity of diesel generator operating in slower fluctuations would be double sized to fill in fast fluctuations with amplitudes of 1% of the maximum fluctuation [2]. The fast-ramp-rate energy storage systems such as flywheel, batteries or supercapacitors

2

The hybrid system management is based on filtering of wind turbine current by two lowpass filters. The optimal current Iwind_opt provided by wind turbine to the DC bus is measured and filtered through the low pass filter HLPF1(p) whose output is subtracted from Iwind_opt to get the reference current IFly_ref of the flywheel regulator (Fig.3). The filter HLPF1(p)cut-off frequency fc1 corresponds to the battery allowed maximum frequency (20mHz in this case). As shown by fig.4, resulting signal from subtraction of measured flywheel actual current IFly from wind turbine current Iwind_opt is filtered by HLPF2(p). The filter output is subtracted from its entrance, to generate the battery reference current Ibat_ref. The filter HLPF2(p)cut-off frequency fc2 corresponds to the battery allowed minimum frequency (1.3mHz in this case). Using this strategy, the diesel generator current is expected to be very low frequency.

The effectiveness of this method relies on the ability of the flywheel to absorb high fluctuations currents. That is why we will focus on the control of the flywheel in the next section.

Bode Diagram 20

Cut-off Frequency : 20mHz Magnitude (dB)

0 1st Order

-20

2nd Order, Q=0.05 -40

2nd Order, Q=0.5 2nd Order, Q=0.707

-60

2nd Order, Q=5 2nd Order, Q=20

-80 0

Phase (deg)

-45 -90 -135 -180

Fig. 3: Flywheel reference current

-5

10

-4

10

-3

10

10

-2

10

-1

0

1

10

10

Frequency (rad/sec)

Fig. 6: 1er and 2nd Low Pass Filters (Fc=20mHz)

5. Simulations results The wind diesel hybrid system is simulated in a runtime of 7000s in various configurations based on low pass filters of first and second order. The effects of three characteristics quality factors (0.5, 0707, and 20) of secondorder filters are analyzed and compared to those of the first order filter of the same frequency. The profile of the battery current generated in the case of a 2nd order filter quality factors of 0.5 is quite close to that with a quality factor of 0707 (Fig.7). When Q varies from 0.5 to 20, it can be observed that current profile of the battery becomes more fluctuant (fig.7 and fig.8), which increases the number of charges / discharges partial cycles and reduce the battery’s lifetime.

Fig. 4: Battery reference current

A First order Low pass filter The first order low pass filter transfer function is expressed by (2).With fc=1/(2.π.τc) the filter cut-off frequency and τ c the time constant [7] [8].

1 H LPF ( p )= 1+ pτ c

(2)

B Second order Low pass filter The second-order filters attenuate more the signal’s magnitude beyond the cutoff frequency. Second order low pass filter transfer function is expressed by (3).

τc p

Q

+ (τ c p )

30

(3) 2

20 Current (A)

1+

Low Pass Second Order Q=0.707 Low Pass Second Order Q=0.5

1

H LPF ( p ) =

Battery Currents 40

With Q the quality factor which effects are illustrated by fig. 5 and fig.5; The signal near cut-off frequency, is amplified when Q is higher than 0.707 and is reduced when Q is lower than 0.5.

10

0

-10

1400

1600

1800 Times (s)

2000

2200

2400

Fig. 7 Zoom of instantaneous battery currents

Bode Diagram 20

Cut-off Frequency : 1.3mHz Magnitude (dB)

0

Fig. 9 and 10 suggest that it is more efficient to use a filter of first order than a 2nd order filter with quality factor equal to 20. As shown by the current spectral density (fig. 10) the quality factor of 20 led to a considerable amplification of the battery current value beyond the cutoff frequency.

1st Order

-20

2nd Order, Q=0.05 -40

2nd Order, Q=0.5 2nd Order, Q=0.707

-60

2nd Order, Q=5 2nd Order, Q=20

-80 0

Phase (deg)

-45 -90 -135 -180 -5

-4

10

10

-3

10

-2

10

-1

10

0

10

Frequency (rad/sec)

Fig. 5: 1er and 2nd Low Pass Filters (Fc=1.3mHz)

3

1

10

Battery Currents

Battery Currents spectrum 100

30

90

20

80 70 Magnitude

10 Current (A)

1st order LPF 2nd order LPF(Q=0.5)

0 -10 -20

60 50 40 30

-30

20

Low Pass Second Order Q=20 Low Pass Second Order Q=0.707

-40 1200

1300

1400

1500

10

1600 1700 Times (s)

1800

1900

2000

2100

0 -4 10

-3

10

-2

10

Fig. 8 Zoom of instantaneous battery currents

-1

10 Frequency (Hz)

0

10

1

10

2

10

Fig. 11: Battery currents spectrum Battery Currents

Battery Currents spectrum 100

60 Low Pass Second Order Q=20

2nd order LPF(Q=0.707) 2nd order LPF(Q=0.5)

90

Low Pass First Order

40

80

30

70

20

60

Magnitude

Current (A)

50

10 0 -10

50 40 30

1900

1950

2000 2050 Times (s)

2100

20

2150

10

Fig. 9 Zoom of instantaneous battery currents

0 -4 10

10

-3

-2

10

-1

10 Frequency (Hz)

0

10

1

10

2

10

Battery Currents spectrum 200

Fig. 12: Battery currents spectrum

1st order LPF 2nd order LPF(Q=20)

160

Magnitude

140 120 100 80 60 40 20 0 -4 10

-3

10

-2

10

-1

10 Frequency (Hz)

0

10

1

10

2

10

Fig. 10: Battery currents spectrum

Fig. 11 to 13 show good attenuation at frequencies above 20mHz in the case of first order filters and second-order filters with factor of quality between 0.5 and 0.707 have similar spectra at frequencies above 1.3mHz and good attenuation at frequencies above 20mHz. Against by the magnitude spectrum is amplified at frequencies below 1.3mHz by second-order filters. The amplification increases with the decrease of the quality factor from 0.707 to 0.5 (fig.11 and 12). The flywheel power spectrums (fig. 13 and fig.14) show the complexity to analyze the filtering effect because of the high density and low amplitudes of the wind turbine fluctuating currents. It is observed that the magnitudes of the spectrums frequencies below 20mHz are higher than those of the flywheel allocated frequency band.

Indeed, it is finding that the magnitudes of the battery currents are greatly larger and less dense than those of flywheel currents at frequencies below 20mHz. Beyond this frequency, the flywheel currents are larger and denser. The filtering principle has therefore effect to allocate to the flywheel the densest and lowest magnitudes signals along the frequency band of wind currents, thereby reducing the number of battery partial cycles (fig.15). The flywheel spectrum does not show much difference in the case of the first filter and second order filter with factor equal to 0.5 and 0707. Exchanged currents between wind turbine, diesel engine and load with first order filter control are illustrated by fig.16. Flywheel Currents spectrum 30 1st order LPF 2nd order LPF(Q=0.5) 25

20 Magnitude

180

15

10

5

0 -4 10

-3

10

-2

10

-1

10 Frequency (Hz)

0

10

1

10

Fig. 13: Flywheel currents spectrum

4

2

10

Flywheel Currents spectrum

complete discharge and recharge of a fully battery [9][10]. From data provided by the batteries manufacturers, the approximated expression of the curve of numbers of cycles depending on the amplitude of discharge current is calculated (eq.4) [10]

30 2nd order LPF(Q=0.707) 2nd order LPF(Q=0.5) 25

15

10

5

C F = a1 + a2 e − a3 R + a4 e − a5 R -3

-2

10

-1

10

0

10 Frequency (Hz)

10

1

10

10

Fig. 14: Flywheel currents spectrum Flywheel and Battery Currents 60 Flyw heel (2nd Order LPF Q=0.5)

50

Battery (2nd Order LPF Q=0.5)

40

Current (A)

30 20 10 0 -10 -20 -30 -40

0

1000

2000

3000 4000 Times (s)

5000

6000

7000

Fig. 15: Flywheel and Battery currents Currents 700 Wint turbine current Load current Diesel engine current

600

Magnitude (A)

400

300

200

100

1000

2000

3000 4000 Times (s)

Avec, CF : cycles to failure ai : constants to determine R : cycles magnitude (fraction of discharge depth); The lifetime is then calculated using the Palmgren-Miner rule [10] [12] considering that the fraction of life consumed during a given cycle is 1/CF. When the sum of the numbers of cycles multiplied by their corresponding fraction exceeds 1, it is considered that the battery is dead and must be replaced. Considering Ni cycles of 20 fractions of discharge depths classified from 0.05 to 1. (i.e. 5% to 100%) the total damage, D, over a given time, is expressed by (5). 20 1 (5) D = ∑ Ni CF ,i

i =1

500

0 0

5000

6000

7000

Fig. 16: Wind turbine, Diesel and Load currents

As shown by simulation results, the battery charge state typically does not follow a regular pattern of cycle and spectral profiles analysis are not sufficient to conclude on the effectiveness of a filter in relation to others. The study will be completed by calculating the lifetime of the battery. Method used in this paper supposes that the number of cycles a battery can tolerate is only function of the discharge current depths.

With CF,i the number of cycles to failure according with class i. Assuming that the battery returns to its initial state of charge after two hours of the hybrid system’s operation, simulation results can be extrapolated to the years to determine the battery life. The software KiBaM elaborated by RERL [13] is used to determine the numbers of cycles according to their depths for each type of filter (fig. 17 to 20). Table 1 presents the OPzS battery lifetime with different filter types. A second order filter with quality factor equal to 0.5 tends to ensure better life of the battery. Battery Current (2nd Order LPF Q=0.5)

Battery Current (2nd Order LPF Q=0.707)

6

6

5

5

4

3

2

1

The counting algorithm called rainflow cycle counting is applied to the time series of the battery currents to find the individual charges/discharges cycles, which are then stored into classes (usually 20 classes are used) of equal size. The classes correspond to different depths of discharge, and the last class corresponds to the 5

(4)

2

No. of Cycles (cycles)

0 -4 10

No. of Cycles (cycles)

Magnitude

20

0

4

3

2

1

0

0.01

0.02

0.03

0.04 0.05 0.06 Cycle Depth

0.07

0.08

0.09

0.1

0

0

0.01

0.02

0.03

0.04 0.05 0.06 Cycle Depth

0.07

0.08

0.09

0.1

Fig.17 Number of cycles Fig.18 Number of cycles

Battery Current (2nd Order LPF Q=20) 12

5

10

No. of Cycles (cycles)

No. of Cycles (cycles)

Battery Current (1st Order LPF) 6

4

3

2

1

0

8

6

4

2

0

0.01

0.02

0.03

0.04 0.05 0.06 Cycle Depth

0.07

0.08

0.09

0.1

0

0

0.01

0.02

0.03

0.04 0.05 0.06 Cycle Depth

0.07

0.08

0.09

0.1

Fig.19 Number of cycles Fig.20 Number of cycles Table 1. Battery Lifetime’s estimation for different filter types. Low 1st 2nd 2nd 2nd Pass Order Order Order Order Q=0.5 Q=0.707 Q=20 Filter Lifetime (years)

1.8

3.1

2.2

1

6. Conclusion In this paper, it is showed that wind power due to high frequencies is considerable and that taking it into account increases wind-diesel hybrid system performance. As against, while the diesel should compensate all wind current fluctuations it should be double sized compared to where it should carry only lowfrequency wind currents. In addition, a diesel operating at high dynamics lifetime is considerably reduced. Big savings will be made so by using storage devices to mitigate the wind currents of medium and high frequency but much lower amplitude compared to the low frequencies current. The original strategy proposed is based on frequency allocation of diesel, battery and flywheel missions by cascading two first order filters. Second order filters quality factors and first order filters effects are studied. Spectral density analysis is completed by battery lifetime estimation to compare the effectiveness of each filter configurations. Simulations results highlight the interest to operate battery with second order filter. References [1] T. Burton, D. Sharpe, E. Bossanyi, “Wind Energy Handbook”, John Wiley & Sons, Ltd, 2001 [2] J. Apt, “The spectrum of power from wind turbines”, Journal of Power Sources 169 (2007) 369-374; [3] W. Li, G. Joós, and C. Abbey, “Attenuation of Wind Power Fluctuations in Wind Turbine Generators using a DC Bus Capacitor Based Filtering Control 6

Scheme” Industry Applications Conference, 2006. 41st IAS Annual Meeting. Conference Record of the 2006. IEEE Vol. 1, Oct. 2006 Page(s):216 - 221 [4] C. Nichita, D. Luca, B. Dakyo and E. Ceanga, “Large band simulation of the wind speed for real time wind turbine simulators”, IEEE Transactions on energy conversion, Vol 17, No. 4, Dec 2002; [5] B. G. Rawn, P. W. Lehn, M. Maggiore, “Control Methodology to Mitigate the Grid Impact of Wind Turbines”, IEEE Transaction on Energy Conversion, Vol 22, N°2, June 2007 [6] M. A. Tankari, B. Dakyo, C. Nichita, “Improved sizing method of storage units for hybrid wind-diesel powered system”, Power Electronics and Motion Control Conference, 2008. EPE-PEMC 2008. 13th, 1-3 Sept. 2008 Page(s) : 1911 – 1917 [7] J. Max, “Pratique du filtrage, filtrage analogique”, Techniques de l’ingénieur, R1102 [8] G. Lissorgues, “Filtres actifs Synthèse et réalisation”, Techniques de l’ingénieur, E 115 [9] H.Bindner and al., “Lifetime Modelling of Lead Acid Batteries”, Risø National Laboratory, Roskilde Denmark, April 2005. [10] A.D. Hansen and al.; “Models for a standalone PV system”; Riso National Laboratory, Roskilde December 2000 [11] T. Christen, Martin W. Carlen, “Theory of Ragone plots”, Journal of Power Sources 91, pp210–216, March 2000 [12] S. H. Baek and al., “Fatigue life prediction based on the rainflow cycle counting method for the end beam of a freight car bogie”, International Journal of Automotive Technology, Vol. 9, No. 1, pp. 95-101 (2008) [13] Http://www.ceere.org/rerl/