The parameters required for building a wind power plant have been calculated using the fuzzy logic method by means of Wind Atlas Analysis and Application Program (WAsP) in this study. Overall objectives of the program include analysis of raw data, evaluation of wind and climate, construction of a wind atlas, and estimation of wind power potential. With the analysis performed in the application, the average wind velocity, average power density, energy potential from micrositing, capacity factor, unit cost price, and period of redemption have been calculated, which are needed by the project developer during the decision-making stage and intended to be used as the input unit in the fuzzy logic-based system designed. It is aimed at processing the parameters calculated by the designed fuzzy logic-based decision-making system at the rule base and generating a compatibility factor that will allow for making the final decision in building wind power plants.
Various methods are being used for the determination of wind power potentials. One of the most important of these is WAsP (Wind Atlas Analysis and Application Program), which is made in the Denmark Riso National Laboratory and used to generate the wind atlas of the European continent (EWA) [
WAsP is being used for the determination of regional wind atlas statistics and energy yield potential by using the wind speed and direction information as well as the obstacles around the wind observation station, land surface roughness, and land topographical characteristics.
Fuzzy logic is a multiple-logic system that is developed against the binary-logic system and assigns membership grades to variables used in everyday life and determines the proportions at which events occur. Everything is right or wrong in two-valued logic. In fuzzy logic, everything, including truth, is only a matter of degree.
Wind is formed by heating and cooling different surfaces at different speeds on the earth from the sunrise until sunset. Kinetic energy of the air in motion is called as wind power. Since the humans are being concerned with the environment, the interest in utilisation of renewable power is increasing. Wind power among the renewable energy is being used widely [
In this study, analysis of wind measurement data obtained at 10 m height and 10-minute interval has been carried out with WAsP, and the wind speed, average power density, form parameter, scale parameter, dominant wind direction, and the Hellmann coefficient have been calculated in accordance with the Weibull distribution. While planning to install the wind turbines to the plant, the required regional wind atlas has been generated. Meteorological data from weather stations outside urban areas are sometimes used, where wind measurements are not available [
In this micrositing study, a wind plant with 12 MW in the region has been planned to be established, and the use of a 2 MW V80-type turbine belonging to the VESTAS company available in the WAsP library has been planned for this establishment. Annual energy generation values of the turbines have been calculated as a result of the micrositing study performed in this way.
The capacity factor, which is the most important parameter during the definition of wind energy potential of one region, is identified as the proportion of energy generated by a wind power plant to the energy that has to be generated at nominal power [
WAsP has been used during the analysis of wind data and energy generation potential. WAsP is a computer program developed in the Denmark Riso National Meteorology Laboratory which is performing analysis with the assumption that wind speed data comply with the Weibull distribution with 2 parameters.
WAsP carries out some analyses by assessing four various basic data in its submodels. Hourly wind data, roughness data of the region, obstacle data in the near vicinity, and topographical data of the region are the basic information used by the program. For the wind farm location, it is also necessary to enter data such as the power curve, thrust coefficient curve, turbine hub height, and turbine rotor diameter in the program [
It is possible to group them under four main titles of general purpose of WAsP so as to analyze the raw data, generate the wind atlas, and assess the wind climate and wind energy potential.
WAsP is performing the time-series analysis by adjusting the raw meteorological data including wind direction and speed data. Weibull parameters are being calculated with the performed analysis [
It is possible to carry out wind climate assessment in any region by using the wind atlas generated by WAsP or by using the data series or data from other reliable sources. Wind climate is being assessed with Weibull parameters and regional distribution of the wind [
WAsP uses some submodels while determining the wind potential. These submodels are obstacle screening, orographic, and roughness exchange models.
The logarithmic wind profile is only applicable where the surface is homogeneous. Pressure changes as per the average surface tensions and surface conditions of surface wind speed until where the gradient force is equal to the fraction force. For the roughness length values belonging to two different surfaces, the following equation may be written for the increase of the limit layer height values:
The wind profile is deformed below the level
In this equation, friction speed at the point taken into consideration is the surface friction speed.
The wind profile is deformed at close distances and at lower parts of the flow of wind coming from the turbine. For this, these objects creating an obstacle effect should be separately studied. Generally, there should be a distance of two objects’ height on the upper wind part of the obstacle and five objects’ height on the lower wind part of the obstacle. Since the meteorological stations are being impacted by the obstacles in the vicinity, this model corrects the fault caused by the obstacle on the data.
The results of wind tunnel studies on simple two-dimensional objects have the following equation:
In this equation,
If there are more than one object, they have been assessed with sectors formed with angles of 30°.
In this model, on the measured wind, data faults of local impacts caused by the topographical structure are being corrected [
For a given radius
For the specific problems, coefficients shall be calculated from surface kinematical limit conditions:
Functions are as Fourier–Bessel series.
In order to make the wind power more competitive against other power generation methods, bigger wind turbines are being designed and established on the clusters especially on the overseas regions called as the wind plant [
To be able to establish a wind power plant, both meteorological and financial parameters are needed. This also shows that one wind plant establishment actually requires to think and use more than one discipline at the same time. For this reason, during the analysis, WAsP has been used which has been already used for the preparation of the European Wind Atlas and Turkish Wind Atlas.
This study is aimed at finding out the average wind speed, average power density, energy yield potential obtained as a result of micrositing, capacity factor, amortization period, and unit cost price required for the establishment of the wind power plant as a result of the performed analysis. These obtained values are being planned to be used to generate a compliance factor for the establishment of the wind power plant by scaling in the rule basis prepared in the fuzzy logic method. In Figure
Matlab/Simulink model of the designed system.
Mathematical location is between 39°05′ N latitude and 38°30′ E longitude. We can examine the province in three parts in terms of surface features: Mountainous areas in the western and northern parts of the province The section east of the district center Medium-height and slightly rugged Arapgir rubble forming the southern part [
As a result of continental climate, most part of the areas of the province are coated with steppe. This means that forest land is very small. Existing forest lands are not in good quality. These are mostly very small oak forests. Satellite view taken from Google Earth of Arapgir Province is shown in Figure
Satellite view of Arapgir Province.
In this study, wind speed and wind direction data belonging to 10 m of the height taken from OMGI (Automatic Meteorological Observation Station) are being used between January 2014 and December 2015.
Monthly average of the measurements taken at 10 m of the height in the region is shown in Table
Monthly average wind speed measured at 10 m height.
Month | Average wind speed at 10 m height (m/s) |
---|---|
January | 1.6 |
February | 2 |
March | 2.6 |
April | 2.8 |
May | 3.1 |
June | 3.8 |
July | 4.1 |
August | 3.6 |
September | 2.9 |
October | 1.9 |
November | 2.5 |
December | 2.2 |
In Table
Average wind severity and average power density of the region.
Calculated values | |
---|---|
Average wind severity (m/s) | 2.75 |
Average power density (W/m2) | 53 |
Twelve sectors have been used in Figure
Wind rose belonging to Arapgir Province.
Most frequently, the distribution used for the calculation of wind power potential is the Weibull distribution. This distribution has been found by the Swedish physicist Waloddi Weibull. This distribution is considerably flexible and simple and also complies with the real data. In other words, since the Weibull distribution is in compliance with wind speed data, it is generally accepted in wind power analysis [
In order to find Weibull parameters (
After examining the wind data in detail and performing the above-given calculations, according to the Weibull distribution in Figure
The Weibull distribution belonging to Arapgir Province.
Wind speed measurements are generally being made at different heights than at the tower heights of wind turbines. For this reason, these measured wind speeds are being extrapolated to turbine tower heights by using the formula known as 1/7 wind power law [
In this study, wind speed values obtained at 10 m in Arapgir with ten minutes of interval have been moved to 80 m which is the turbine hub height with the calculated Hellmann coefficient.
The Hellmann coefficient can be taken as 1/7 in the most general case. However, the other ways can be used to determine it more accurately [
Calculation in terms of roughness length:
Calculation in terms of speed and height:
Calculation of roughness length and speed:
Equation (
Change of the Hellmann coefficient as per wind direction.
In this study, wind atlas statistics have been prepared belonging to the plant area with the help of WAsP software by using the frequency distribution table obtained from wind speed and direction data belonging to 10 m of the height, obstacles in near vicinity, roughness data, and digitized map with a scale of 1/25.000 representing the region topography.
Data imported into WAsP are being analyzed, and the average wind severity map of the region in Figure
Average wind severity map calculated in the area.
Power density map of the area.
While planning the installation of wind turbines in the site, wind atlas is being used. By using WAsP and wind atlas statistics, locations where the power generation amount will be high on the digital map may be defined from colour distribution.
In this study, it has been planned to establish a wind plant of 12 MW in the area, and for this establishment, it has been planned to use 2 MW of the V80-type turbine belonging to the VESTAS company within the WAsP library. In order to minimize the interference of turbines with each other, the micrositing study has been performed.
During this micrositing study, design has been made so as to keep the track area losses of the turbines on each other in the minimum level. Furthermore, while designing the optimum turbine placement so as to obtain a maximum yield, wind turbines have been placed on the dominant wind direction not only taking into consideration the generation amount but also the operability limits of the turbines in technical terms. In Figure
Micrositing process of turbines.
Following this micrositing study with 12 MW installed power where 6 VESTAS V80 wind turbines are used, electricity generation of annual gross 4.1316 GWh and net 38.352 GWh has been calculated. As a result of the interference of turbines with each other, it has been understood that there will be a 7.17% loss. Annual generation values have been shown in Table
Annual generation values of 6 V80 wind turbines.
Parameter | Total | Min | Max |
---|---|---|---|
Net annual generation (GWh) | 38.352 | 5.713 | 7.082 |
Gross annual generation (GWh) | 41.316 | 6.034 | 7.897 |
Wake (track) loss (%) | 7.173 | — | — |
In Table
Annual generation and loss values of 6 V80 wind turbines.
Turbine number | Hub height (m) | Net annual generation (GWh) | Wake (track) loss |
---|---|---|---|
Turbine 1 | 80 | 6.789 | 7.09 |
Turbine 2 | 80 | 6.145 | 6.18 |
Turbine 3 | 80 | 5.824 | 9.62 |
Turbine 4 | 80 | 7.082 | 5.34 |
Turbine 5 | 80 | 6.799 | 8.29 |
Turbine 6 | 80 | 5.713 | 9.79 |
As can be shown in Table
The capacity factor, which is an important parameter and has to be known both by the generators and consumers, is the division of energy generated in a specific time frame to the maximum energy that can be generated at that specific time frame [
In this study, the capacity factor has been calculated for the turbine used during micrositing and turbine selection studies and generated micrositings:
In the micrositing study, the annual calculated generation amount for six 2 MW V80 wind turbines, maximum amount that can be generated from turbines, and capacity factor are shown in Table
Capacity factor of the area.
Turbine model | Turbine number | Installed power | Annual calculated production amount (GWh) | Nominal power to be produced (GWh) | Capacity factor (%) |
---|---|---|---|---|---|
VESTAS V80 | 6 | 12 | 38.352 | 105.12 | 36.48 |
In case of operation under nominal power of six 2 MW turbines during 8760 hours in one year, they may generate maximum 8760 × 12 = 105.120 GWh power in one year. As a result of the calculations made in WAsP, it has been determined that the wind power generated by the wind power plant formed by six 2 MW turbines shall generate 38.352 GWh power in one year, and the capacity factor shall be realized as 38.352/105.120 = 0.3648, in other words 36.48%. The capacity factor which is calculated at and above 25% for the wind turbines deemed appropriate for the investment in Turkey [
One of the most important studies that have to be carried out while establishing a wind turbine to a region is the calculation of kWh power cost. Generally, the cost of one wind power project per kWh is found by proportioning the annual total cost to the annual power generation amount.
The annual power generation amount changes depending on the parameters such as the hub height of turbine, rotor diameter, average wind severity of the area, and annual cost may be correlated with turbine price, turbine foundations, inner site road construction costs, investment costs, and project lifetime. In Figure
Cost structure per kWh for wind power projects.
In this study, within the first establishment, turbine cost, network connection, foundation and establishment costs, electricity installation, and road construction and control systems have been considered as mentioned in the report of the European Wind Power Union, and land rental, consultancy, maintenance and repair, and yearly licence costs have been included into the analysis as permanently paid costs within years (Table
Turbine establishment cost [
Establishment parameters | Cost (€) |
---|---|
Turbine and installation | 9.000.000 |
Project preparation costs (permission and licences) | 650.000 |
Land rental cost | 100.000 |
Network connections (cable, transformer, and communication) | 1.800.000 |
Foundation and installation (foundation, road, and site) | 1.200.000 |
Finance (consultancy, insurance, and bank) | 450.000 |
Total | 13.200.000 |
Annual 38.352 GWh power generation has been calculated for the wind power plant formed with six 2 MW VESTAS V80 wind turbines. In this case, annual income has been calculated as 2.109.360 € with 5.5 €cent/kWh electricity sales price. Emission sales income for the 38.352 GWh power generation has been calculated as 82.608 € annually. As a result of the calculations made, annual income and expenses as well as the annual margin are shown in Table
Annual income-expense table [
€/year | |
---|---|
Electricity sales income | 2.109.360 |
Emission sales income | 82.608 |
|
|
|
|
Rental amount | 100.000 |
Maintenance-repair insurance | 130.000 |
Labour | 30.000 |
Annual licence cost | 1.000 |
General administrative costs | 5.000 |
|
|
|
|
|
|
|
|
In this study, operational expenses, maintenance, repair, insurance, rental, and licence costs, and some values changing as per year such as interest have been reflected as average values in the income-expense table by taking into consideration the project lifetime. In case of selling the electricity at 5.5 €cent/kWh at the wind power plant formed by 6 VESTAS V80 wind turbines, it is seen that the cash flow will pass to the positive state at 7 years and the project returns at year 7.
In this study, cost per kWh for the capacity factor obtained from selected turbines has been calculated. By using the capacity factor as 36.48% calculated for the VESTAS V80 turbine used in the established system, the annual generation amount of the plant in kWh has been found.
In order to calculate the unit cost of electricity power obtained from wind turbines, it is necessary to know the regaining factor of investment. The regaining factor of capital (
Generation cost is as follows:
The power to be generated by the wind turbine (kWh/year) in the region (
In the calculations made, unit cost values to be used as the input unit in the fuzzy logic system are found as 4.26 €c/kWh.
In this study, average wind speed and wind direction data of Arapgir province measured at 10 m height and 10-minute intervals and taken from OMGI (Automatic Meteorological Observation Station) were analyzed by WAsP software. In the analysis made, the average wind speed is found as 2.75 m/s, the average power density is found as 53 W/m2, the figure parameter is found as 1.51, and the scale parameter is found as 4.2 m/s obtained at 10 m according to the Weibull distribution.
As a result of the calculations made in WAsP, 12 sectors have been used in the wind mill and dominant wind direction has been observed as SE at 120° sector. The change of the Hellmann coefficient as per the direction by using existing data in addition to the analysis made as per the wind blowing has been examined, and the average Hellmann coefficient has been calculated as 0.34.
As a result of this micrositing study with 12 MW installed power where 6 VESTAS V80 wind turbines are used, with the wind power plant to be established in the region, it has been calculated that annual gross 41.316 GWh and net 38.352 GWh electricity shall be generated. It has been understood that 7.17% loss shall occur as a result of the interference of wind turbines with each other.
In case six 2 MW turbines operate under nominal power during 8760 hours, a maximum of 105.120 GWh energy may be generated. As a result of the calculations made in WAsP, it has been determined that the wind power plant formed by six 2 MW turbines shall generate 38.352 GWh energy in one year, and the capacity factor shall be realized as 36.48%.
Annual generation values of the plant have been found as 4.26 €c/kWh by using 36.48% capacity factor calculated for the VESTAS V80 turbine used in the established system. The amortization period is calculated as 7 years. Unpredictable expenditure was selected as 2% of the initial installation cost. The height of the tower of the wind turbine was determined as 80 m. The wind speed at 80 m height was calculated as 7.1 m/s. The power density at the height of 80 m was calculated as 309. In Figure
Display of the result values generated by the fuzzy-based decision-making system designed according to the input values calculated by the WAsP.
The data used to support the findings of this study are available from the corresponding author upon request.
The authors certify that there are no actual or potential conflicts of interest in relation to this article.