Operation of a Wind Turbine-Flywheel Energy Storage System under Conditions of Stochastic Change of Wind Energy

The paper presents the issues of a wind turbine-flywheel energy storage system (WT-FESS) operation under real conditions. Stochastic changes of wind energy in time cause significant fluctuations of the system output power and as a result have a negative impact on the quality of the generated electrical energy. In the author's opinion it is possible to reduce the aforementioned effects by using an energy storage of an appropriate type and capacity. It was assumed that based on the technical parameters of a wind turbine-energy storage system and its geographical location one can determine the boundary capacity of the storage, which helps prevent power cuts to the grid at the assumed probability. Flywheel energy storage was selected due to its characteristics and technical parameters. The storage capacity was determined based on an empirical relationship using the results of the proposed statistical and energetic analysis of the measured wind velocity courses. A detailed algorithm of the WT-FESS with the power grid system was developed, eliminating short-term breaks in the turbine operation and periods when the wind turbine power was below the assumed level.


Introduction
Environmental issues included in long term power strategies of different countries and high accessibility of the renewable sources of energy and their significant potential are the main reasons for an increase in the share of renewable sources of energy in the global generation of electrical energy. When it comes to the widely available solar and wind energy, one should however pay attention to great fluctuations of the converters output power related to a stochastic nature of the irradiation changes and wind velocity V in time. The instability has a negative impact on the cooperation of wind and solar sources with the power grid system [1][2][3]. The issue is important for systems with a high percent share of renewable sources of energy, particularly in those without output power stabilisation [2,4]. The operation of unsustainable sources of energy can cause problems related to stabilisation of a section of a power grid system and generate additional costs related to maintaining the periodically activated conventional sources in a standby mode [5,6].
Energy storage in industrial applications is a current issue and the research in the area led to some practical applications of batteries, artificial and natural compressed air energy storage (CAES), supercapacitors, superconducting magnetic energy storage (SMES), flywheel energy storage, and so forth, [4,[7][8][9][10][11][12][13]. Despite technical sophistication and high costs, their application area in high power systems is gradually extending. There are a growing number of technical devices including energy storage and further recovery as a part of normal operation, for example, emergency supply systems, pumped storage power plants, and hybrid and electrical cars. Due to growing significance of such kind of solutions used in economics of highly developed countries the problem should be considered with regard to widely understood optimization and economic aspects [11,14].
An important issue of the practical application of wind sources is to mitigate the effects of output power fluctuations resulting from the stochastic nature of the wind energy changes in time, when working with a power grid system [15]. Long-lasting (e.g., for several hours) breaks in power generation in wind sources related to the decrease in the wind kinetic energy can be determined with the use of computerassisted systems of the output power prediction [16][17][18]. The situation differs when the breaks are short (up to several minutes) and impossible to predict due to their generation mechanism. Moreover, the frequency and duration of short breaks depend on the parameters of the implemented wind turbine-mainly on its cut-in velocity V cut-in . Figures 1 and 2 present two circadian curves of wind velocity changes obtained by measurements. The measurements were made on 1 and 5 March 2008 in a meteorological station in Strzyżów (South Eastern Poland) at the height of ℎ = 10 m above the ground level.
The obtained values were recalculated for the height of ℎ WT = 60m above the land level, corresponding to the position of the wind wheel hub used for further analyses of Enercon E53 turbine, according to the following relationship: where is the aerodynamic coefficient of terrain roughness, ℎ , ℎ WT are the heights of the wind velocity and wind wheel hub measurement, quoted in reference to the land level, V is the wind velocity at the measurement height ℎ , V WT is the recalculated wind velocity. The model of wind velocity vertical profiles expressed in the formula (1) is simplified but sufficient for the purpose of the study.
The above-mentioned type of wind turbine with nominal power of WTN = 800 kW was used in all calculations and simulations done for the purpose of the study. The horizontal line marked in Figures 1 and 2 stands for the cut-in velocity, which for the reference type of turbine is V cut-in = 2 m/s. The analysis of the curve presented in Figure 1 indicates an almost 5-hour break (long-lasting break) in power generation and many short turbine cut-outs from the power grid system. For the course presented in Figure 2 the average wind velocities are higher, which allows for the uninterrupted operation of the generator 24/7. In practice the nature of the wind velocity changing in time tends to include the features of both courses, and the mean energy additionally depends on the  deterministic components: circadian and annual changes and long term trends.
At the current technological level of energy storage production, solving the problems related to the first type of breaks seems hard and economically not justified. Nevertheless, short breaks in the operation of wind sources can be effectively compensated with energy from appropriately selected energy storage, resulting in a partial stabilisation of the output parameters of a wind power plant connected to the grid, and also contribute to improving the quality of the generated electrical energy [15,20,21].

Cooperation of Energy Storage with
Wind Turbine 2.1. Introduction. With regard to the breaks in the wind turbine operation caused by the stochastic nature of the wind velocity (energy) changes in time, it is necessary to find some engineering solutions preventing the related power cuts to the grid. Considering significant technical difficulty related to eliminating long-lasting power cuts, the methods allowing for preventing power cuts with maximum duration MAX can be considered sufficient but the parameter value is usually determined for a range up to several minutes. One of the practically feasible methods is maintaining the power supplied to the electrical energy system at the assumed level with the value of 3MIN in the reference periods. In order to maintain high quality of energy and stable operation of the system in the connection spot, the change in the power level supplied to the system from the value related to the wind turbine power curve to the value of 3MIN should be as smooth as possible. It is also recommended to implement measures aimed at partial stabilisation of the system output power at the level of 3MIN in periods with reduced wind energy. Such situations occur when the wind velocity value is the velocity corresponding to the power 3MIN + PW ( PW is the house The Scientific World Journal 3 load power of the analysed system). A complete stabilisation of the output power of the WT-ESS for all periods of the turbine operation at a reduced power requires using complex engineering systems and is economically not justified (very high investment expenditure).
The paper assumes that the implementation of the presented measure requires the use of energy storage with appropriate parameters and of appropriate type. The basic aim is to compensate the reduced power supplied by the wind turbine generator in the assumed periods with duration up to MAX . The solution is of particular importance for geographical areas with the values of the average wind velocity not much higher than the cut-in velocity V cut-in of the applied type of turbine. Appropriately selected turbine and energy storage leads to a creation of wind turbine-energy storage (WT-ESS) of a new quality, connected to the power grid, whose features on the one hand result from its being a renewable source of energy but on the other hand are similar to the characteristics of conventional sources [20,22,23].

Characteristics of Energy Storage Systems, Selecting the Energy Storage Type, Advantages and Disadvantages of Kinetic
Storage. Accumulation of energy is a topical and economically expensive problem, of high technological complexity [11,12]. The studies carried out in this field result, among others, from aspiration to improve energetic safety and from the need of long-term accumulation of very large amount of energy. The difficulties in accumulation of electric energy cause the indirect methods are most commonly used. Consequently, it reduces the efficiency of the process. Among the systems that make most often use of the abovementioned method there are accumulator batteries (leadacid and lithium-ion batteries), kinetic storage (flywheels), supercapacitors, superconducting magnetic energy storage (SMES), and compressed air energy storage (CAES) [4,11,12,24].
In case of the storage designed for operation in renewable energy systems the requirements related to their energetic capacity, rated power, charging rate, and durability, and the range of operating temperature results from specific conditions of wind turbine and photovoltaic panel operation, caused directly by weather conditions. The changes in temperature, humidity, pressure, and so forth not only directly affect the equipment but also contribute to stochastic changes of input values delivered by the aforesaid types of the sources.
It was assumed for purposes of the research that functionality of the energy storage in electric power grids is described with the use of a set of parameters including: power and energy densities (W/L and Wh/L, resp.)-determining possible recovery of usable current (power) and energetic capacity, durability (the number of charging-discharging cycles), depth of discharge, the range of operating temperature, discharge rate and transition rate between the operating states, efficiency, unit cost of the equipment, converted to power or unit energy (cost/kW, cost/kWh) and physical dimension of the system. Table 1 presents a comparison of the most important usable parameters of the above mentioned energy storage types [10,12].
In order to carry into effect the algorithm proposed in the paper and aimed at partial stabilization of the power delivered to the system from a wind source a storage is necessary which renders possible a so-called short-term accumulation of energy. It is designed for equalizing the output power of the system in time intervals below 1 h (usually 0,25 h). Such systems are required to deliver the energy to the electric power grid immediately after activation of the storage (with very short delay) and to maintain it at the rated power level in the assumed time [10]. Taking into account a single wind turbine, an energy storage cooperating with it should have the average energetic capacity (usually from tens to hundreds kWh), high charging rate (comparable to discharging ratein the range of minutes), rated power in the range from tens to hundreds kW, very short time of transition between charging and discharging stages (below 1 s) and the range of operating temperature corresponding to yearly temperature variations characteristic for the definite geographic location. Moreover, the storage should be composed of modules allowing for simple development of the system [25].
The SMES storage must be excluded from cooperation with wind turbines due to their low energy density (0,5 Wh/L ÷ 10 Wh/L). Usable current value of a single module reaches even several kA (the superconducting technology and significant reduction of active power loss), nevertheless, the time of cooperation with the system is too short as compared to the one required according to the assumption. Similarly, the CAES storage is excluded too, due to the need of building large systems (pressure vessels) or using a precisely imposed location of the system (natural reservoirs, e.g., old mine excavations, etc.), and relatively poor efficiency of the system. Such a type of the storage is characterized by too long deployment time (from several to 10 minutes) as compared to real dynamics of wind energy variations. On the other hand, high power density is an advantage of this storage type. Nevertheless, in case of the time of energy recovery below one hour this advantage is not decisive for the choice of the storage type [10]. From the group of considered solutions of the problem the supercapacitors must be removed too. This is caused by very low energy density (2 Wh/L ÷ 10 Wh/L) which precludes gaining proper capacity and maintaining the output power at required level within the time from ten to twenty minutes.
Hence, the most important types of energy storage feasible for practical application of the proposed method of equalizing the output power of a wind turbine with stochastic character of the input function are secondary electrochemical cells-accumulators and flywheels [10,12].
Among the advantages of the flywheels as compared to electrochemical cells (lead-acid and lithium-ion batteries) there are constant value of energetic capacity in the whole range of operational temperature (−35 ∘ C to +40 ∘ C), covering yearly variations of weather conditions, very high number of charging and discharging cycles reaching millions (lifetime 15-20 years) and short duration of storage charging (approximating the discharging time with rated power) [10,12,16]. Two first features allow to locate the storage in direct proximity of the turbine and to operate it without any restrictions within the turbine lifetime (15-20 years). High 4 The Scientific World Journal  The Scientific World Journal Power system charging rate [16] enables to use the wind energy even in case of quick variations, without the need of using faster energy storage devices as energetic buffers. Additionally, the kinetic storage is characterized by high efficiency (from 80% to 95%), remarkably higher as compared to lead-acid batteries (75%-80%). For the recent solutions their efficiency is higher even than the one of lithium-ion batteries (83%-86%). It should be noticed that the system occupies relatively small space-a group of modules may be often closed in a container ready for transportation to another location [10,12].
One of the features of the kinetic storage, that might be considered as a fault as compared to accumulator battery is lower energy density (in case of lead-acid battery from 50 Wh/L to 100 Wh/L, while for the lithium-ion one-from 200 Wh/L to 350 Wh/L). Another fault of them is due to high degree of self-discharge (several percent per hour). Nevertheless, the above-mentioned features are not decisive for cooperation between the wind turbine-energy storage system and the electric power grid, since the storage is not required to be characterized by very large energetic capacity and the storage charging and discharging processes last below 1 hour-usually no more than twenty minutes. The investment cost of flywheels, converted to unit power or unit energetic capacity, is several times higher than that of the lead-acid or lithium-ion batteries. Hence, economical aspects of the use of such systems must be considered as their fault worsens appraisal of the technology of kinetic storage [10,12].
Obtaining high energy values requires a high flywheel velocity, which entails the use of modern composite materials. Their density is several times lower than the density of steel, and the boundary strength max related to the presence of high radiation forces is much higher, which results in obtaining the value of characteristic energy several times higher (W/kg). Detailed information on this matter is presented in the paper [26]. Low idle changes and a relatively high total system performance (usually of ca. 86%) are mainly achieved by using magnetic bearings and the rotor operation in a vacuum with the pressure values of about 10 −3 bar [7].
Based on the comparison of technical parameters of the above-mentioned types of energy storage and considering the economic aspects (periodical replacement of batteries) a flywheel type of energy storage was assumed for cooperation with the wind turbine [9].

Algorithm of a Flywheel Energy Storage Cooperation with a Wind Turbine (Farm).
According to the established assumptions, a wind turbine with the nominal power WTN and specific power curve 1 = (V ) working with flywheel energy storage form a complex power system (WT-FESS). Its basic goal is to deliver a relevant level of active power to the power grid system also in the periods when the wind velocity V is below V cut-in . The basic diagram of a flywheelelectrical system is presented in Figure 3. The kinetic energy of wind is transformed in the turbine wheel into the shaft (or gear) and generator rotary motion. According to the turbine power curve, active power 1 ( ) is obtained at the system outlet. The storage operates with the active output power 2 ( ) variable in time; the power can be positive (energy released to the power grid-unloading), negative (energy taken from the generator-loading) or zero energy (idle state of complete unloading of the storage). Hence the storage energy ES ( ) also varies in time and its value ranges from zero to the nominal capacity ES . The current energy value tends to be expressed in the percentage of nominal value, with the use of factor ES% ( ).
Active power 3 ( ), which is an algebraic sum of momentary powers 1 ( ) and 2 ( ) with deducted house load power PW , is released to the system. Due to an automatic change in the WT-FESS configuration, its value also varies in time PW = ( ). According to the assumptions given in Section 2.1 while releasing energy from the storage to the grid, the minimum output power value 3MIN is obtained. However, it covers periods of time with a specific duration (maximum duration MAX ) and depends on meeting several conditions given further on in the algorithm.
The system presented in Figure 3, depending on the momentary value of the wind velocity V ( ) and the energy 6 The Scientific World Journal storage loading ES ( ), can be in one of the four characteristic states: The transition between the above-mentioned states is a continuous and dynamic process, depending on the stochastically changing atmospheric conditions and the current and previous system arrangement. A single continuous operating period of energy collecting from flywheel energy storage is limited with the MAX algorithm parameter.

Selecting the Energy Storage Volume for
Working with a Wind Turbine

Statistical Energy Analysis of the Course of Wind Velocity
Changes V = ( ). Based on theoretical analysis and the conducted tests it was determined that the measurement courses of the wind velocity changes V = ( ) can be used for identifying the minimum capacity of the flywheel energy storage ESMIN that will meet the assumptions of the algorithm of WT-FESS cooperation with the power grid system, according to Section 2. PW and the turbine cut-out velocity V cut-out were marked. This way the course V = ( ) is divided into four areas, where a set of statistical and energy parameters characterising the WT-FESS in the specific geographical location can be determined.
In the area 1 the wind velocity meets the requirement V ( ) < V cut-in , and the generator power is 1 = 0. In practice such periods can last from several seconds to many days. In order to identify the required capacity of a flywheel energy storage ESMIN information about subsequent breaks of the specific type and their average duration is necessary. The parameters proposed and used in further analysis for the area include the average 1AVG and maximum 1MAX duration of power generation breaks (stochastic wind velocity changes) not exceeding the set value of the factor MAX , 1 series coefficient determining the average number of subsequent breaks separated with one turbine operation interval at power guaranteeing the energy storage loading ( 1 > 3MIN + PW ) and the summary turbine operation time in the area 1WT for the assumed period of analysis .
Area 2 covers the wind velocity range meeting the requirement V cut-in < V ≤ V 3MIN . Information concerning the average m 2AVG and maximum 2MAX duration of intervals not exceeding the set value of MAX , the average generator power 1AVG2 and the total turbine operating time in the area 2WT for the assumed period of analysis is determined in the area.
The system operation in area 3 (wind velocity V ≥ V 3MIN ) allows for controlled loading of the storage according to its current energy status ES ( ). The average generator power 1AVG3 and the total turbine operating time in the area 3WT for the assumed period of analysis is determined for the area.
Area 4 covers the turbine cut-out periods due to excess wind velocity V ≥ V cut-out , which can additionally cause mechanical damage. Moreover, the following values of electrical energy generated by the reference type of turbine are determined for the total period and areas 2 and 3: WT , 2WT , and 3WT , respectively.
According to the description above, sets of measurement points whose values constitute the average wind velocity from the period Δ and the duration of 48 seconds are analysed. Hence 1800 measurement points are recorded within 24 hours, and their number amounts to 657 thousand within one year. For high power wind turbines (hundreds kW and more) the moments of inertia of rotating elements are so high that the quoted measurement period is sufficient for the goals presented in the paper. All measurements used in the paper were made with a rotating anemometer placed at 10 m above the land level.
From the point of view of the analysed subject matter it is important to compare the values and relationships between the suggested statistical energy parameters for two characteristic periods of a calendar year: autumn-winter and springsummer. For many geographical locations, including the South Eastern Europe, the autumn-winter period has greater wind energy that the spring-summer one and the differences can be of several dozen percent. Another important element covers determining the impact of the change in the WT-FESS input and output parameters, in particular in the parameter of time MAX

Identifying the Boundary Capacity ESMIN of a Flywheel Energy
Storage. The WT-FESS operation according to the assumptions of the algorithm presented in Section 2.3. requires using a flywheel energy storage with appropriate capacity. The author's research on the analysis of the measurement courses of the wind velocity changes V = ( ) for a period of several years for one geographical location lead to determining an empirical relationship identifying the value of the minimum storage capacity ESMIN that guarantees correct operation of the analysed system. The relationship includes technical parameters of the storage and wind turbine and statistical energy parameters of the measurement courses of the wind velocity changes defined in Section 3.1.
The presented relationship consists of segments corresponding to the turbine operation areas separated in Figure 3. A corrective segment related to the storage additional loading conditions and its ability to use the excess energy generated by the turbine ( 1 > 3MIN ) was also taken into account. Considering these elements in determining the minimum capacity ESMIN of a storage intended for working with a selected type of wind power plant in a specific geographical location, the following relationship was proposed: where 1AVG , 2AVG , 3AVG is the upper ( index) and lower ( index) confidence limit for the subsequent mean time values 1AVG , 2AVG , and 3AVG (Tables 2(a), 2(b), 3(a), and 3(b)), ES % is the idle losses of the flywheel storage expressed in percent of its nominal power ES , AVG is the upper confidence limit of the storage operation on idle gear (the value stands for the mean time between subsequent periods of the storage energy use in areas 1 and 2 whose duration does not exceed the maximum natural unloading time storage ESR ), ME+ , ME− are the flywheel energy storage performance in the loading and unloading process, 2 is the correction factor ( 2 = 0 for 3MIN ≤ 2AVG and 2 = 1 for 3MIN > 2AVG ), 3 is the coefficient of the storage additional loading conditions: identifying the turbine power margin that can be used during the storage additional loading, where 1MIN stands for the minimum turbine power value corresponding with the wind velocity V cut-in , 4 is the ability to use excess power: The other factors and parameters used in the relationship (3) are described in the previous section of the paper. The first three components of the relationship (3) help determine partial capacities related to stabilisation of a power plant output power for areas 1 and 2 at the established maximum continuous duration of the turbine operation with reduced power ( 1 < 3MIN ) and idle loses of the flywheel energy storage Δ ES ( 2 = 0, ES ( ) > 0). The last element is of corrective nature and in special cases reduces the value of the identified capacity. Additionally, it happens that the real capacity of the storage ESMIN must not be lower than the ESMIN determined from the relationship (3) and in practice depends on the nominal data of the modules available for the selected storage type and the possibility of their combining. [Days] [

MWh]
[ %] 1 Jan. The Scientific World Journal 9 The calculations use the power curve and other E53 turbine parameters presented in the manufacturer's technical catalogue [19].

Changes in the Capacity ESMIN in the Function of WT-FESS Parameters.
A computational application was developed with the use of the analysis algorithm of the measurement courses of wind velocity changes V = ( ) proposed in Section 3.1 and empirical relation (3), in the .NET environment (language C#). With regard to a large number of measurement points covering the period of one year and the related long times of statistical analysis, the Task Parallel Library was used, for parallel execution on multicore system which allowed to significantly reduce the total time of calculations.
With the use of the developed application, families of characteristics ESMIN = ( MAX ) and 1 = ( MAX ) were determined for the established set of power values 3MIN and particular geographical location. Based on them it is possible to evaluate the behaviour of the WT-FESS when wind turbines with identical nominal power are used, to differentiate the mounting height of the wind wheel and to analyse the system for different periods of the same year and to compare several years. The above-mentioned families of characteristics were determined separately for two periods of the same year: autumn-winter and spring-summer. The conducted calculations used the values of standard deviations and confidence ranges, assuming the confidence factor of 0,95, which were determined for statistical and power parameters presented in Tables 2(a), 2(b), 3(a), and 3(b).    (Figures 7 and 8). The differences in the characteristics curves 1 = ( MAX ) between the springsummer and autumn-winter period result from different average wind velocity and the dynamics of the wind velocity changes in time. Analysing the obtained characteristics one can note their similarities within the dynamics of the ESMIN storage capacity changes for both analysed periods. The determined capacity ESMIN for the spring-summer period is higher than for the autumn-winter period, which is mainly caused by higher average values of the wind velocity  It was assumed that the correctness of determining the minimum capacity of a flywheel energy storage ESMIN intended for working with a wind turbine is established based on the value of a percentage factor of eliminating the acceptable cutouts % . It is the relationship between the summary working time of a generator with power below 3MIN in unit periods and duration not exceeding MAX , compensated with the flywheel storage energy, and the summary time of all periods of the generator operating at a power not exceeding 3MIN and duration not exceeding MAX (including not compensated periods), in the assumed period of analysis , expressed in percent.

Simulation of WT-FESS Operation under Conditions of Stochastic Wind Energy Change
A set of wind velocity values discrete in time is the simulator input obtained by measurements. According to Section 3.1. of the paper, each measurement point makes the average wind velocity for the period Δ 48 seconds long.
In the numerical algorithm of the simulator, regardless of the energy storage operation state, one should consider idle losses related to mechanical resistance in the system, feeding of magnetic bearings and maintaining the specific vacuum level in the rotating mass housing. If the energy storage is in an idle state they are taken into account as ES % factor. At loading and unloading the idle losses are included in the process efficiency, whereby the efficiency was assumed as identical in both cases and its value is ES .
The momentary power of a wind turbine generator 1 ( ) is determined with the use of the energy curve stored in a discrete form in the database. The values of the generator power are determined for each of the established points separating the time periods Δ ( ) for = 1, 2, . . . , − 1. For the initial ( ) and final ( ) time of the Δ ( ) period, wind velocities amounting to V ( ) and V ( ) respectively and the generator power 1 ( ) and 1 ( ) related to them are determined. The average turbine power in the range Δ ( ) and value 1AVG( ) is used for the calculations made in the WT-FESS operation simulator. The changes in the energy storage power 2 ( ) are established based on the relationships from (2a) to (2d), whereas the output power 3 ( ) of the system is identified based on the determined values of 1 ( ) and 2 ( ) and the house load power PW ( ).
The energy state of the storage in discrete moments of time for = 0, 1, 2, . . . , is determined based on the initial storage loading condition (for = 0, ES ≥ ES0 ≥ 0), previous changes in the storage 2 ( ) and turbine 1 ( ) power, its efficiency, and coefficient of idle losses. The value of energy for discrete time ( = ⋅ Δ ) is determined by adding (considering the sign) the energy gains in all time ranges Δ preceding the point. The storage energy in the moment of time can, thus, be expressed as where is the time step index, is the final time step index used according to the relationship = ⋅ Δ , to determine the time , ES is the nominal power of energy storage, 2( ) is the established value of the energy storage loading or unloading power as the average value for the initial and final point of the time range Δ , , , ∈ {0, 1} are the coefficients from sets , , and , respectively identifying the storage state for the time periods (loading, unloading, idle).
For numerical implementation of proposed model .NET platform, MS Visual C# language and ADO.NET technology for handling the relational database of the wind turbines parameters were used. Elements of object-oriented software were applied for building the programme structures. A library of classes intended for representing the structure and operating principle of the following WT-FESS elements: wind turbine, flywheel energy storage, control system, method of selecting ESMIN storage capacity and identifying the storage energy state at any moment of time were developed. In relation to a very time-consuming nature of the calculations covering a statistical energy analysis of the discrete course of wind velocity changes in time, elements of calculation paralleling were used. That is why Task class was used to divide the calculations onto logical cores of the processor intended for PCs and workstations.

Results of Simulation Analyses.
Simulation tests of a WT-FESS working with the power grid system were carried out for two types of inputs test input V WT = ( ) and real input V = ( ). Two configurations of the system with different nominal power ES , limit capacities ESMIN and initial loading states ES0% of the storage (option I and II- Table 4) were used for the tests. The real input case is covered by parameters presented in Table 4 as option III. ENERCON E 53 turbine with the power of WTN = 810 kW and established generation characteristics was used in all tests.
The first part of the tests was done for the input V WT = ( ), whose curve is presented in Figure 11(a). The analysis covers changes in the wind velocity during 70 minutes, including fluctuations from the cut-in velocity V cut-in , to the velocity V when the turbine reached the nominal power WTN . The velocity changes V WT in time were selected so that in the assumed period of analysis the system WT-FESS reached all working states defined in the defined algorithm (Section 2.3) and shifted between them at diversified dynamics.
The other part of the tests covered a simulation of the investigated system operation for a real input in a form of the curve of wind velocity changes from the one indicated in the geographical location reference for the period between 3 March and 6 March 2008. The nominal (limit) capacity ESMIN of the storage used for the tests was determined for an identical location but using measurement data for the springsummer period in 2010.
According to the assumptions presented in Section 2.3., the numerical simulator model covers four operating states of the system depending on the wind energy, system parameters, and current and previous values of the energy storage. The results of the performed simulations were presented in a form of power curves of the generator 1 ( ), storage 2 ( ) (considering the sign) and the output power of the system  I  200  50  100  1800  100  2  0.5  II  100  0  75  1800  100  2  0.5  III  100  0  150  600  100  2  3 ( ) and a relative percent storage loading ES% ( ) for the assumed period of analysis . Figure 11 shows the results of WT-FESS operation simulation conducted for the test input and two parameter options of the tested system (Table 4). With regard to the short period under analysis and the related high readability in Figures 11(b)-11(d), the curves for the aforementioned parameters are presented simultaneously for two simulation options ( Table 4).
As a result of the wind velocity drop below V cut-in in the period between 37 and 57 minutes, if the turbine works independently, it is disconnected from the power grid system (Figure 11(a)-circled with an intermittent line). However, considering the turbine cooperation with the storage, the break was eliminated thanks to the previously stored energy (Figures 11(b) and 11(c)). For option II, considering the assumption of zero storage energy at the beginning of the analysis period ( ES0 = 0), the stored energy was not sufficient to eliminate the entire break, which resulted in the turbine cut-out after 20 minutes. A similar situation occurred in the first period of the system operation (to ca. minute 4). The enumerated periods are circled with an intermittent line in Figures 11(c) and 11(d). It is the evidence of too low capacity of the applied energy storage, resulting from extremely difficult storage operating conditions not included in the confidence ranges of statistical energy parameters used in the relationship (3). Figure 12 shows the curves of some selected simulator parameters for WT-FESS operation at real input (option III- Table 4).
The analysis of the system operation for a real input covers 50 hours from the period between 3 March 2008 and 6 March 2008, with diversified wind conditions (Figure 12(a)). Next to high wind energy periods (e.g., between the system operation hour 5 and 20) there are periods with boundary energy values from the point of view of the assumed WT-FESS operation parameters (e.g., between hour 20 and 30). This type of periods accumulates breaks in the turbine operation, which are short according to the definition presented in Section 1 of the paper and should be additionally compensated with 14 The Scientific World Journal energy stored in the storage. Furthermore, a period of long-lasting decrease in the wind velocity below the cut-in velocity (between system operation hour 31 and 34) can be additionally seen in Figure 12, whose impact on the system operation will not be analysed in detail.
From the point of view of the developed algorithm, the most important periods are the ones with boundary (limit) values of the wind velocity (energy). The implemented algorithm of WT-FESS cooperation with the power grid system assumes stabilisation of the output power 3 of the system at the assumed level 3MIN , besides eliminating short breaks. It applies to periods where the wind velocity allows for reaching the turbine power 3MIN > 1 > 0 (area 2 in Figure 4) and the assumed duration up to MAX . In the analysed period the greatest number of wind velocity changes corresponding to the transition between areas 1 and 2 ( Figure 4) occurs between hour 15 and 25 of the system operation. This period is circled with an intermittent line in Figures 12(c)-12(e). Unloading of the storage energy is used for eliminating breaks in the turbine operation ( 1 = 0) and equalising the system output power 3 with the value of 3MIN (Table 4, option III) assumed in the algorithm. It is also loaded between the storage unloading periods (positive power 2 ), when the power values 2 are negative (Figure 12(c)).

Comments and Conclusions
Operation of wind sources in geographical locations with moderate wind conditions may generate a number of problems related to their cooperation with the power grid system. The basic reason for such occurrence is stochastically changing kinetic energy of the wind and construction characteristics of the turbines. One of the solutions to mitigate the effect of frequent cut-outs of such sources from the grid is using energy storage. Implementing the proposed algorithm of the wind turbine can control the system operationflywheel energy storage system cooperation with the grid that allows for eliminating a large number of short breaks using the previously stored energy. The author proposed an algorithm using the features of flywheel energy storage, mainly the short period of their loading and shifting between the loading and unloading state, as well as low dependence of the real capacity on temperature. Equalising the active power released to the power grid system at the assumed level 3MIN is done for the breaks in the turbine operation and periods when the turbine reaches the power 1 < 3MIN at maximum duration MAX . The results obtained by simulation (Figures 11 and 12) are the evidence of good efficiency of the developed algorithm and improving the conditions of the wind turbine cooperation with the power grid system. The number of the turbine cut-outs from the grid at appropriately selected flywheel energy storage capacity decreases significantly, which results in an improved quality of electrical energy and the source stability.
Correct operation of the above-mentioned system requires determining the minimum (boundary) capacity ESMIN of the applied energy storage. The process can be conducted in different ways but the author of the paper suggests a proprietary concept based on statistical energy analysis of the measurement time series of changes in the wind velocity in the analysed geographical location for a period of at least one year (Tables 2(a), 2(b), 3(a), and 3(b)). The minimum capacity of the storage ESMIN required for the assumed algorithm at maintaining the specified parameters of cooperation with the power grid system is established based on the empirical relationship (3) connecting the energy storage and wind turbine parameters and states, as well as the results of statistical energy analysis of the measurement curves V ( ). Seasonality of the average wind energy demonstrated based on the tests (Tables 2(a), 2(b), 3(a), and 3(b)) indicated the need to consider this fact in determining the limit storage capacity ESMIN . The simulation results confirm that if this fact is accounted for while establishing the value of ESMIN the real percent index of eliminating the acceptable breaks (duration up to MAX ) is between 75% and 85%. Not meeting this condition results in a significant decrease in the process of eliminating short breaks in the wind turbine operation defined in the paper.
In the author's opinion the statistical energy parameters proposed and determined for the measurement curves can be compared and taken into account while designing WT-FESS systems in various geographical locations. Based on the values of the parameters presented in Tables 2(a), 2(b), 3(a), and 3(b) one can draw more detailed conclusions on the nature of wind conditions in the examined location (energy, dynamics of changes, etc.), similarly to the wind conditions class according to IEC 61400-1. As a result of implementing heuristic methods it is additionally possible to select the optimum components of the WT-FESS (turbine type, tower height, type and size of storage) as regards the unit cost of electrical energy generation.
It was established based on the conducted statistical energy analyses of the curves V = ( ) (Tables 2(a), 2(b), 3(a), and 3(b)) and the tests according to the implemented method of determining the capacity ESMIN that for a specific geographical location, conclusions concerning mutual relations between the parameters characterising the WT-FESS and cooperation with the power grid can be formulated. With this in mind, a series of calculations was made whose results are presented as curves ESMIN = ( MAX ) at 3MIN = const (Figures 4 and 5), and ESMIN = ( MAX ) at ℎ = const ( Figure 6). The coefficient of series 1 has a major impact on the capacity value ESMIN and the shape of the enumerated characteristics. Considering the dependence of the coefficient The families of the aforementioned curves are typical of a particular geographical location, the parameters of the system elements ( WTN , ESN , ℎ TW ) and its cooperation with the power grid ( MAX , 3MIN ). They can be used for an approximate determination of the minimum (limit) capacity of the storage ESMIN when different values of the wind wheel mounting height, power change 3MIN and time of the eliminated breaks MAX are used.
The choice of energy accumulation system in the form of flywheels is an effective solution that enables to fulfill the assumptions formulated for the algorithm of WT-FESS system cooperation with the electric power grid. Exchange of the storage for accumulator batteries would worsen the system properties because of long charging time (the lead-acid batteries), capacity variations (particularly in winter) and shorter lifetime (in higher temperature). On the other hand, the use of supercapacitors would result in significant growth of the cost, since they should be distinguished by high electric capacity. Hence, it appears that despite the disadvantages mentioned in Section 2.2 the kinetic energy storage complies with the largest number of required qualities. Moreover, development of the technology allows forecasting reduction of the kinetic storage prices in the future and their more common use, particularly in the field of renewable power engineering.