Next Article in Journal
Policy Analysis of the Challenges to an Effective Switch to Low-Carbon Energy in the Economic Community of West African States
Previous Article in Journal
Architecture for Co-Simulation of Transportation and Distribution Systems with Electric Vehicle Charging at Scale in the San Francisco Bay Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applicability of WorldCover in Wind Power Engineering: Application Research of Coupled Wake Model Based on Practical Project

1
Beijing RETEC New Energy Technology Co., Ltd., Beijing 100079, China
2
Department of Electrical Engineering, School of Electrical and Electronic Engineering, North China Electric Power University, Beijing 100096, China
3
CSSC Windpower Development Co., Ltd., Beijing 100097, China
4
Department of Atmospheric Science, School of Atmospheric Science, Nanjing University of Information Science and Technology, Nanjing 210044, China
5
Institute of Engineering Thermophysics, Chinese Academy of Sciences, Beijing 100190, China
*
Authors to whom correspondence should be addressed.
Submission received: 2 February 2023 / Revised: 20 February 2023 / Accepted: 22 February 2023 / Published: 24 February 2023

Abstract

:
This paper discusses how the incorporation of high-resolution ground coverage dataset ESA WorldCover into a wind flow field and wake simulation calculation, as well as the use of the coupled wake model for wind farm output simulation, can improve the accuracy of wind resource assessment using engineering examples. In the actual case of grid-connected wind farms in central China, SCADA wind speed data is reconstructed to the free flow wind speed in front of the wind turbine impeller using the transfer function of the nacelle, and the wind farm is modeled using OpenWind software, simulating the wind speed at the height of each wind turbine hub and each wind turbine output. The results show that when other initial data are consistent, using ESA’s high-precision land cover dataset WorldCover 10 m to make roughness lengths which improves the wind farm output simulation accuracy by 8.91%, showing that it is worth trying to apply WorldCover 10 m to the wind farm simulation design. At the same time, this case is used to compare and analyze the application of the Eddy-Viscosity wake model and the two coupled wake models based on the Eddy-Viscosity wake model. The results show that the coupled wake model will have higher accuracy than the Deep Array Eddy Viscosity wake model and it is 1.24% more accurate than the Eddy Viscosity wake model, and the ASM Eddy Viscosity wake model is 5.21% more accurate than the Eddy Viscosity wake model.

1. Introduction

As people become more concerned about environmental protection and green energy, the wind power industry continues to thrive. Wind energy consumption increased by 88% between 2002 and 2017 in the world’s 17 major wind energy-consuming countries, including China. During this time, China’s renewable energy consumption increased by 39%, India by 23%, Ireland by 21%, and the United Kingdom by 20% [1]. China commissioned 17 GW offshore wind power capacity in 2021, accounting for 80% of the global offshore growth, making China’s total offshore wind capacity 27.7 GW [2], rapidly surpassing the UK and Germany to take the lead globally in terms of offshore wind energy capacity with a 44.8% market share [3]. China promised in 2020 that it would achieve carbon neutrality by 2060. Wind energy’s carbon dioxide emission intensity in China is more than 98% lower than that of traditional fossil fuels, according to studies, and the emission reduction effect can reach 84–98% [4]. The world’s three major economies, the European Union, China, and the United States, have all implemented effective incentives to promote the development of the wind power industry and increase the proportion of green power electricity [5].
The evaluation of Chinese wind resources has gradually progressed to fine assessment, and a more accurate and reliable assessment of resources will help to develop wind resources more reasonably, save costs, and increase output. Some Chinese researchers have emphasized that, despite becoming the world’s largest wind power installer in the last two decades, the quality of installed capacity has been subpar [6]. As a result, it is both meaningful and necessary to improve the accuracy and reliability of resource assessment. Engineers must use advanced technology and resources to obtain high-resolution and high-precision wind resource maps in wind farms to support the construction of wind power projects.
The surface cover’s most direct influence on near-surface airflow is to change the near-surface flow field, thereby changing the near-surface wind speed profile. Many scholars have already conducted extensive research in this area. Neff et al. (1998) discovered that densely vegetated ridges and hills can cause nonlinear interactions [7]. Jacobs et al. (1995) measured the vertical and horizontal changes in wind speed in the maize canopy and discovered a 20–30% difference between the average wind speed and the horizontal average [8]. This research is extremely practical. Corn is a popular economic crop in northern China. According to the findings by Fu et al. (2020), the influence of changes in soil geology on surface wind is limited, changes in the surface cover have a significant impact on simulated wind speed, and roughness length is the most important factor influencing wind speed [9]. Wen et al. (2014) demonstrate that improving the accuracy of ground cover can improve the accuracy of wind speed simulation [10]. According to the findings of Frank Baier et al. (2022), the higher the resolution of the ground cover data used, the smaller the error of wind speed simulation, especially near the surface [11]. Matthew J.C. et al. (2012) discussed the effect of surface roughness on wind turbine wake turbulence at a given atmospheric stability in their study [12]. D.L. Elliott and J.C. Barnard (1990) investigated the variability of wind speed and turbulence intensity as a function of wind direction and surface roughness across the site [13].
According to the findings, the surface cover has a significant impact on wind resources and wind turbine wake simulation, and changes in the surface cover and data resolution are also important factors. To improve the accuracy and reliability of the simulation, high-resolution surface cover data should be utilized in wind resource assessment.
With the rapid development of wind power projects in recent years, researchers have not stopped researching and developing wake models, shifting from a single wake model to a coupled wake model [14]. Another question that wind resource engineers have is whether the coupled wake model is more exact and can help improve wind resource assessment accuracy.
At the moment, the most commonly used land cover products in wind resource assessment are NASA GeoCoverLC90m resolution dataset (2009) [15], European Space Agency (ESA) Copernicus Global Land Service (CGLS) land cover 100 m data set (2015–2019) [16], and China National basic Geographic Information Center global surface cover data product GlobeLand30 (2020) [17]. This paper will use the ESA’s 10 m resolution global surface cover data released in 2020 to make roughness length as the initial data for wind farm software simulation design, and use different coupled wake models to calculate the wind power output. The experimental project selects a grid-connected wind farm in central China, and the experimental results are compared to the actual data. The study of this paper is expected to provide a reference for accurate wind resource assessment.

2. Background

2.1. Data Source

The surface cover data used in this paper are from the European Space Agency’s WorldCover project. The main outcome of this project is the release in 2020 of a free global land cover (GLC) product with a resolution of 10 m based on Sentinel-1 and Sentinel-2 data. It includes 11 different types of land cover and has been independently verified with a global accuracy of 74.4%. In Asia, the overall accuracy is 80.7%. The accuracy of the categories of tree cover, snow and ice, farmland, water, and bare/sparse vegetation is higher, exceeding 80%. The accuracy of grassland and architecture is moderate, while shrubs, wetlands, and moss/lichens are less accurate. Moss, lichens, and grasslands are overestimated, while trees are underestimated [18].
The 10 m resolution surface cover data products, which are also based on Sentinel satellite data, include Google’s DynamicWorld product coverage [19]. Water, trees, cultivated land, and other coverage types are divided into nine categories, with an overall accuracy of 71.3% Landcover coverage [20]. With an overall accuracy of 85%, the coverage type includes nine categories, such as water, grassland, ice, and snow. Although LandCover appears to have higher overall accuracy, studies have shown that WorldCover has higher accuracy and LandCover accuracy is lower at lower cell resolution (100 m2) [21,22,23].
The National Basic Geographic Information Center of China’s land cover data GlobeLand30 (2020) is widely used in the current wind resource assessment work. GlobeLand30 products cover an area of 80° S–80° N with a resolution of 30 m. Ground cover is classified into ten types, including cultivated land, forest, water, and artificial cover. Among them, the minimum classification accuracy of six types, including cultivated land and forest, is controlled at or above 70%, while the lowest classification accuracy of four types, including water bodies and artificial cover, is controlled at or above 80%. In China, the precision of the data set is 80.4% [17,24,25].

2.2. The Context of the Experimental Project

The experimental project Yingchen wind farm is located in Qingling Mountain in central China, in the north temperate zone, and forest coverage in this area can reach 37.1%. The natural woodland is dominated by deciduous broad-leaved forests such as Robinia pseudoacacia, willow, Toona sinensis, paulownia, and so on, while the economic woodland is dominated by walnut, persimmon, peach, and other fruit forests, and the crops are primarily wheat, corn, sweet potato, and so on. The wind farm has a capacity of 48 MW and employs 146–3000 wind turbines from a Chinese manufacturer. The wheels are 90 m tall and were connected to the grid in 2020.
At the project site, obtain WorldCover (2020), LandCover (2020), and GlobeLand30 (2020) ground coverage data. Figure 1 and Figure 2 depict the surface coverage of wind power sites in the area and its surrounding areas.
As shown in Figure 2, WorldCover, LandCover, and GlobeLand30 are basically the same in terms of overall surface cover classification, but LandCover data lacks woodland cover data; WorldCover and GlobeLand have a better approximation in the distribution of all types of surface cover, while WorldCover shows more ground details, which is confirmed by the GoogleEarth satellite map.

2.3. Making the Roughness Length

To use the ESA WorldCover surface cover data set in wind resource assessment, the mapping relationship between surface cover data and surface roughness must be established. The mapping relationship between ground coverage and roughness length is established using the European Wind Atlas [26] and the New European Wind Atlas [27], and the roughness length is adjusted based on the local actual situation (Table 1) [28,29,30,31,32]. Expand the 10 km range from the site location to the east, west, north, and south to account for the impact of upstream surface cover on the experiment [33]. When the wind blows over the top of the mountain, a pressure difference exists between the windward and leeward sides of the mountain, increasing friction loss and effective roughness length [34]. Because there is no dense vegetation near the wind measurement site, displacement height is not required. The EWA roughness length value and the improved roughness length value were used to assign the land cover data to generate the roughness length map that the program could read (Figure 3).

3. Verification Mode

Accurate energy calculations are critical in wind power engineering. Improving the precision of output estimation can help to improve the precision of investment return expectations. Wind resource engineers must select a suitable wake model based on the characteristics of different wind farm terrain environments when calculating energy, and when the wake model performs simulation calculations of energy, ground roughness is a very dependent parameter. As a result, in the same simulation environment, a high-precision land cover map can improve the wake model’s calculation accuracy.
The Yingchen wind farm, which has been connected to the grid, is chosen as the actual measurement project in this paper. According to this project’s feasibility study report (FSR), the Openwind [35] wind farm design software is also used for modeling, the WindMap [36] wind flow model is used to calculate the wind flow field, and the Eddy-Viscosity wake model is used to calculate energy. The initial data uses the same wind measurement data, surveying and measured topographic map, and the same uncertainty table. On this basis, the roughness length is calculated using EAS land cover data, and the energy is calculated using the Deep Array EV wake model and the ASM-EV wake model. It should be noted that the wake model’s general parameters use the same settings, while the non-general parameters use the default parameters without modification. The experiment’s energy data is compared to the wind farm feasibility study report data and the wind farm’s actual operation data. On one hand, it confirms that WorldCover data can be used in the wind power industry. On the other hand, it compares the errors of different wake models. The flowchart process will be used for verification in this paper (Figure 4).

3.1. Wake Model

The Eddy-Viscosity wake model (EV) [37] is a linear CFD calculation that calculates the wake of an axisymmetric structure wind turbine using the time-averaged Navier–Stokes equations (RANS). The model employs cylindrical coordinates and assumes incompressible airflow.
The Eddy-Viscosity wake model automatically observes mass and momentum conservation in the wake. The average Eddy-Viscosity of each downstream wake segment is used to calculate the wind-turbulence interaction. The wake in the Eddy-Viscosity model has a Gaussian distribution, and its recovery depends on the turbulence intensity. The more turbulence there is, the more the wake wind mixes with the free wind around it, and the wake travels and recovers faster [33]. Its equation for calculation is as follows:
1 v v 0 = D M   e x p { 3.56 ( r b ) 2 }
where v 0 represents the average free stream wind speed, v represents the wind speed at a distance r from the wake centerline, D M   represents the initial wind speed attenuation at the wake centerline, and b represents the wake width parameter, and the equation is as follows:
b = 3.56 C T 8 D M ( 1 0.5 D M )  
where D M is related to thrust coefficient C T and turbulence intensity, the expression is as follows:
D M = C T 0.05 ( 16 C T 0.5 ) I A 1000  
where I A is the ambient turbulence intensity (%).
The Eddy-Viscosity wake model was used to calculate the energy in the feasibility study. Figure 5 depicts the Eddy-Viscosity wake model simulating the wake behind a single wind turbine, which employs a Gaussian distribution to improve the prediction of the wake velocity deficit, and is run by Openwind. The X-axis in the figure is a multiple of the wind turbine impeller’s diameter(D), and the background free wind speed is 10 m/s, which is displayed in red. The legend assumes that the background wind speed is 1(100%), and different colors correspond to different multiples of that speed.
The Deep Array Eddy-Viscosity wake model (DAWM-EV) [38] is a modification and improvement of the Eddy-Viscosity wake mode, a coupling model based on Sten Frandsen’s theory. Deep Array and Eddy-Viscosity are separately modeled in the DAWM model and combined by taking the maximum value of roughness influence and basic wake influence (background roughness is obtained by reading the imported rough map).
An infinite array of wind turbines is represented as a uniform region of high surface roughness in Frandsen’s theory. The roughness drags on the atmosphere, causing changes in the structure of the planetary boundary layer (PBL) downstream, most notably a decrease in free-stream wind speed at turbine hub height. According to this theory, the equivalent roughness of the wind farm is:
z 00 = h H e x p ( k c t + ( k / ln ( h H / z 0 ) ) 2 )
where h H is the hub height, k is the von Karman constant (about 0.4), z 0 is the ambient roughness between turbines, and c t is the distributed thrust coefficient,
c t = π 8 S d S c C T
where C T is the thrust coefficient, S d is the mean downwind spacings in rotor diameters and S c is the mean crosswind spacings in rotor diameters.
Assume that each wind turbine occupies a distinct area that contributes to increased surface roughness. Increased roughness creates the internal boundary layer as the wind reaches the wind turbine. The wind width line or shear in this internal boundary layer is defined by the wind turbine roughness rather than the ambient roughness, and the velocity at the top of the internal boundary layer must immediately match the velocity above it. After the wind has passed through the wind turbine, a second internal boundary layer will form to indicate the transition back to the environment’s surface. Both internal boundary layers grow with downstream distance, according to the following equation.
h i b l [ ln ( h i b l z 0 ) 1 ] = ( x z 0 1 ) z 0  
In this equation, h i b l represents the height of the internal boundary layer, x represents the distance generated by the internal boundary layer, and z 0 represents the downstream roughness (wind turbine roughness for the first internal boundary layer and ambient roughness for the second internal boundary layer). The downstream wind turbine produces its own internal boundary layer under the upstream wind turbine.
Once the equivalent roughness is defined, the meteorological theory is used to estimate the impact on the hub-height wind speed deep within the array (i.e., where the PBL has reached equilibrium with the array roughness) under the assumption of a constant geostrophic wind speed G and a neutral logarithmic profile throughout the PBL. Taking only the first turbine’s IBL pair and assuming that both IBLs have grown to exceed hub height, the equation for hub-height speed is as follows:
v H v H = ln ( h 1 / z 0 ) ln ( h 2 / z 00 ) ln ( h 2 / z 0 ) ln ( h 1 / z 00 )
v H is the wind speed of the wind turbine hub height downstream of the wind farm, v H is the wind speed of the upstream wind turbine hub height, h1 and h2 are the heights of the first and second internal boundary layers, respectively, and z 0 and z 00 are the roughness of the wind turbine and the environment (background).
Figure 6 depicts the wake behind a single wind turbine as simulated by the Deep Array EV wake model, which predicts wake velocity deficits using a Gaussian distribution in the near-wake region and a linear extension in the far-wake region, and it is run by Openwind. The X-axis in the figure is a multiple of the wind turbine impeller’s diameter(D), and the background-free wind speed is 10 m/s, which is displayed in red. The legend assumes that the background wind speed is 1(100%), and different colors correspond to different multiples of that speed.
The Area Slowdown Eddy-Viscosity wake model (ASM-EV) [39] considers the wind farm as a whole. The wind farm is related to the inversion layer at the top of the atmospheric boundary layer as additional surface roughness and as a gravity wave generator, which means that even though the wind speed of the atmospheric boundary layer above the wind field is uniform, it is very sensitive to the decrease of the lower wind speed and the decrease of energy, so there will be changes in the pressure gradient and the generation of gravity waves. The model calculates the interaction force between the atmosphere and wind farm based on the conservation of kinetic energy and the above-influencing factors, which is a top-down simulation method of fluid kinetic energy. The basic equation is:
v ( t ) v 0 = 1 + ( v v 0 1 ) e x p ( α t )
where v 0 is the free wind speed of the upstream wind turbine at hub height, v is the wind speed of the downstream wind turbine at hub height.
α t = k u * z / Δ z 2  
where z is PBL’s height, h is wind turbine hub height, Δ z = z − h, u * is the friction velocity.
The ASM EV wake model’s single wake wind speed is depicted in the figure below (Figure 7), which ran by Openwind. The X-axis in the figure is a multiple of the wind turbine impeller’s diameter(D), and the background free wind speed is 10 m/s, which is displayed in red. The legend assumes that the background wind speed is 1 (100%), and different colors correspond to different multiples of that speed.

3.2. Nacelle Transfer Function

Wind farm operational data are derived from the SCADA (Supervisory control and data acquisition) system. The wind speed data in SCADA is the nacelle anemometer measurement data. The anemometer is easily affected by the wind rotor and the nacelle because it is installed behind the wind rotor on the top of the nacelle. As a result, the nacelle transfer function must be used to correct it to the free stream wind speed at hub height in front of the impeller. The nacelle transfer function establishes the relationship between the nacelle wind speed and the wind speed measured in the wind measuring tower [40].
In this paper, the wind speed of the anemometer tower and the wind speed of the nacelle are divided into 0.5 m/s intervals in this study [41], and the wind speed in front of the nacelle impeller is calculated using the following equation [42]:
v H = v H , i + 1 v H , i V N , i + 1 v N , i × ( v N v N , i ) + v H , i
where v H is the wind speed at the hub height of the wind turbine, v H , i , and v H , i + 1 is the bin-average wind speeds in bin i and bin i + 1, v N , i , and v N , i + 1 is the wind speed measured by the nacelle anemometer in the bin i and bin i + 1.

3.3. The Sector Affected by the Wake

According to the calculation formula of the influence sector given in the IEC 61400-12-1: 2017 standard, determine the influence sector range of the first row of the wind turbine on the rear row of wind turbines [43].
The calculation equation is as follows:
θ = 1.3 a r c t a n ( 2.5 D n l n + 0.15 ) + 10  
where D n is the diameter of the wind turbine impeller, l n is the distance between wind turbines.

4. Results and Analysis

The table below shows the results of the estimated energy in the Yingchen wind farm feasibility study report, the estimated energy in this experiment, and the actual energy of the wind farm (Table 2). Except for the roughness length data, which is made using WorldCover 10 m, all initial input data in this verification use data from the feasibility study stage, and the values of the uncertainty parameters are all set by the feasibility study stage. The Eddy-Viscosity wake model uses the same settings as the feasibility study stage. According to the results in Table 2, the deviation between the calculated and actual power generated by the Eddy-Viscosity wake model in this verification is 6.76%, which is 8.91% higher than the feasibility study EV model. As a result, the high-accuracy land cover dataset, WorldCover 10 m can be considered to improve measurement accuracy.
In this verification, the energy is calculated using the Eddy-Viscosity wake model and two basic coupled models of the model. According to the results, the deviation between the Deep Array Eddy-Viscosity wake model and actual energy is 5.52%, which is 10.15% greater than the feasibility study report (FSR) and 1.24% greater than the Eddy-Viscosity wake model. The difference between the ASM Eddy-Viscosity wake model and actual energy is 1.55%, which is 14.12% greater than the feasibility study report and 5.21% greater than the Eddy-Viscosity wake model.
To further clarify the differences between the three wake models used in this experiment, a ridge in the wind farm area is chosen and the wake attenuation is analyzed one by one. Although the arrangement is more regular and the number of wind turbines is greater on the southwest side, the influence of the topography is greater, which is not conducive to the analysis of the influence of the wake, and the elevation drop of the ridge on the northeast side of the wind farm is limited. The wind farm ridge wind turbines on the northeast side are used in this experiment to study wake attenuation.
The wind turbines on the northeast ridge are 3# to 8#. Because the 3# wind turbines are not on the same side of the mountain as the other wind turbines, and their position is relatively offset to the east relative to other wind turbines on the same ridge, they are not included in this analysis. Figure 8 depicts the ranking positions of 4#–8# wind turbines. The wake of the upstream wind turbine will affect the downstream wind turbine, and the influence range of the upstream wind turbine on downstream wind turbines can be calculated using the wake influence sector equation, and the north is 0°. Table 3 shows the sector range of the downstream wind turbine affected by the wake of the upstream fan when 4# is the first upstream wind turbine, using the formula given above with reference to the IEC 61400-12-1.
This study reviews the sectors where the downstream wind turbines are impacted by the upstream 4# wind turbine. On 3 August, 2020, 0:45 to 2:45 lasted 2 h, and the wind direction is 135° to 160°, which is covered by the range in Table 3, and it is relatively heavily affected by the wake.
Table 4 displays the average wind speed calculated by different wake models for each wind turbine in this experiment during the above time period and wind direction. According to the trend, the free wind loses a lot of energy after passing through the 4# wind turbine, and the wind speed is greatly reduced when it reaches the 5# wind turbine. The wind blows through the 5# wind turbine, losing some of its energy. The wind is supplemented with energy from the surrounding environment as it passes through the 6# wind turbine and loses some of its energy to reach 7#. At the same time, because 7# is not completely parallel to the first three wind turbines, the wake flow received by the 7# wind turbine is reduced. In addition, the impact is minor. The 8# wind turbine has a greater offset to the east, and the distance from the upstream 7# wind turbine is 367.5 m, while the distance from the upstream 4# wind turbine has increased to 1267.5 m, indicating that the wind has recovered enough energy from the environment.
The wind blows from the 4# wind turbine to the subsequent 5#–8# wind turbine. Combining with Figure 9, it is also discovered that the altitude of the 8# wind turbine has a large drop from the wind turbine arranged in front of it, and the altitude difference from the 7# is 51.38 m, the altitude difference with the 4# wind turbine is 91.96 m, and the wind will eventually reach the 8# wind turbine through continuous climbing. If the free flow wind passes through the 4# position and reaches the 8# position, the wind speed should be continuously increased by the terrain.
As a result of the environment, terrain, and wind turbine arrangement, wind speeds of 4#, 5#, and 6# wind turbines decreased row by row, while wind speeds of 7# and 8# wind turbines recovered noticeably. Table 4 displays the results.
The turbulence at each wind turbine shown in Table 5, when combined with the turbulence tables generated by the wake model, confirms the energy loss and recovery at the wind turbine hub.
The power generation calculated by the ASM Eddy-Viscosity wake model is the closest to the actual energy, as shown by the summary of the calculation results in Table 6. In the third row, the difference between the Deep Array wake model and the Eddy-Viscosity wake model becomes apparent. The Eddy-Viscosity wake model’s simulated wake flow has begun to recover greater than the loss at the fourth wind turbine 7#, and power generation has begun to rise. 8# The wind turbines’ comprehensive geographical location and simulation calculation of the wake flow show that the power generation has a higher recovery than the previous three wind turbines 5#, 6#, and 7#. The Deep Array wake model simulation also shows energy recovery for the 7# wind turbine. Greater than the loss, but recovery takes longer. The ASM Eddy-Viscosity wake model simulation data trend is consistent with the previous two models. Because of its top-down simulation method, the energy lost in the wind farm can only be supplemented by the turbulent kinetic energy above, reducing wind resistance and increasing turbulence. reflected more clearly. The energy of the 5# and 6# wind turbines is decreasing, the energy of the 7# wind turbine is being restored, the energy of the 8# wind turbine is being less affected by the wake, and the energy supplemented from the environment is becoming more sufficient. Figure 10 depicts the average wind speed and wind turbine output power during the analysis period.

5. Conclusions

The following conclusions can be drawn from the preceding analysis:
When other initial data are consistent, using ESA’s high-precision land cover dataset WorldCover 10 m to make roughness lengths improve simulation accuracy by 8.91%, indicating that it is worthwhile to try to apply WorldCover 10 m to wind farm simulation design. When conducting wind resource assessments, it is recommended to use a higher-resolution ground cover dataset and to adjust the roughness length appropriately based on the actual local vegetation coverage.
The coupling wake model outperforms the basic wake model in terms of accuracy. The Deep Array Eddy-Viscosity wake model is 1.24% more accurate than the Eddy-Viscosity wake model, and the ASM Eddy-Viscosity wake model is 5.21% more accurate than the Eddy-Viscosity wake model. Encourage the use of coupled wakes in wind resource assessment work and help improve wake models.
Although the Deep Array Eddy-Viscosity wake model performed less well than the ASM Eddy-Viscosity wake model in this verification, one reason is that it is a wake model designed for large wind farms, whereas the wind farm used in this verification was on a smaller scale. Follow-up studies are needed to determine which of the two performs better in large-scale wind farms.
We use the nacelle in transfer function in this paper to reconstruct the wind speed of the wind turbine, and we also advocate the construction of production wind measuring towers in wind farms to improve the accuracy of wind speed reconstruction, and thus the accuracy of wind power prediction.
For verification, this study employs an actual grid-connected wind farm. Although a relatively ideal experimental result was obtained, more experimental items and data must be used to form a statistically significant result that can be used to fine-tune wind power projects. Wind resource evaluation provides new ideas to help improve wind resource evaluation accuracy, allowing for more accurate budgets, more valuable investment suggestions, and more convincing returns.

Author Contributions

Formal analysis, Z.Y. and Y.C.; resources, J.Y.; writing—original draft preparation, J.Z.; writing—review and editing, H.L.; supervision, Q.L.; Investigation, J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to project data is confidential data within the company.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sadorsky, P. Wind energy for sustainable development: Driving factors and future outlook. J. Clean. Prod. 2021, 289, 125779. [Google Scholar] [CrossRef]
  2. Lee, J.; Zhao, F. GWEC Global Wind Report 2022; Global Wind Energy Council: Brussels, Belgium, 2022. [Google Scholar]
  3. Musial, W.; Spitsen, P.; Duffy, P.; Beiter, P.; Marquis, M.; Hammond, R.; Shields, M. Offshore Wind Market Report: 2022 Edition; National Renewable Energy Lab (NREL): Golden, CO, USA, 2022. [Google Scholar]
  4. Xu, K.; Chang, J.; Zhou, W.; Li, S.; Shi, Z.; Zhu, H.; Guo, K. A comprehensive estimate of life cycle greenhouse gas emissions from onshore wind energy in China. J. Clean. Prod. 2022, 338, 130683. [Google Scholar] [CrossRef]
  5. DeCastro, M.; Salvador, S.; Gómez-Gesteira, M.; Costoya, X.; Carvalho, D.; Sanz-Larruga, F.J.; Gimeno, L. Europe, China and the United States: Three different approaches to the development of offshore wind energy. Renew. Sustain. Energy Rev. 2019, 109, 55–70. [Google Scholar] [CrossRef]
  6. Dong, C.; Qi, Y.; Dong, W.; Lu, X.; Liu, T.; Qian, S. Decomposing driving factors for wind curtailment under economic new normal in China. Appl. Energy 2018, 217, 178–188. [Google Scholar] [CrossRef]
  7. Neff, D.E.; Meroney, R.N. Wind-tunnel modeling of hill and vegetation influence on wind power availability. J. Wind Eng. Ind. Aerodyn. 1998, 74, 335–343. [Google Scholar] [CrossRef]
  8. Jacobs, A.F.G.; Van Boxel, J.H.; El-Kilani, R.M.M. Vertical and horizontal distribution of wind speed and air temperature in a dense vegetation canopy. J. Hydrol. 1995, 166, 313–326. [Google Scholar] [CrossRef] [Green Version]
  9. Fu, D.; Liu, Y.; Li, H.; Liu, S.; Li, B.; Thapa, S.; Tian, C. Evaluating the impacts of land cover and soil texture changes on simulated surface wind and temperature. Earth Space Sci. 2020, 7, e2020EA001173. [Google Scholar] [CrossRef]
  10. Wen, X.; Liao, X.; Yuan, W.; Yan, X.; Wei, Z.; Liu, H.; Dong, W. Numerical simulation and data assimilation of the water-energy cycle over semiarid northeastern China. Sci. China Earth Sci. 2014, 57, 2340–2356. [Google Scholar] [CrossRef]
  11. Baier, F.; Metz-Marconcini, A.; Esch, T.; Schroedter-Homscheidt, M. Impact of higher-resolved satellite-based land cover classification on near surface wind speed forecasts. Meteorol. Z. 2022, 31, 101–116. [Google Scholar] [CrossRef]
  12. Churchfield, M.J.; Lee, S.; Michalakes, J.; Moriarty, P.J. A numerical study of the effects of atmospheric and wake turbulence on wind turbine dynamics. J. Turbul. 2012, 13, N14. [Google Scholar] [CrossRef]
  13. Elliott, D.L.; Barnard, J.C. Observations of wind turbine wakes and surface roughness effects on wind flow variability. Sol. Energy 1990, 45, 265–283. [Google Scholar] [CrossRef]
  14. Stevens, R.J.; Gayme, D.F.; Meneveau, C. Coupled wake boundary layer model of wind-farms. J. Renew. Sustain. Energy 2015, 7, 023115. [Google Scholar] [CrossRef] [Green Version]
  15. Cunningham, D.; Melican, J.E.; Wemmelmann, E.; Jones, T.B. GeoCover LC-A moderate resolution global land cover database. In Proceedings of the ESRI International User Conference, San Diego, CA, USA, 8–12 July 2002. [Google Scholar]
  16. Buchhorn, M.; Lesiv, M.; Tsendbazar, N.E.; Herold, M.; Bertels, L.; Smets, B. Copernicus global land cover layers—Collection 2. Remote Sens. 2020, 12, 1044. [Google Scholar] [CrossRef] [Green Version]
  17. National Geomatics Center of China. 30-Meter Global Land Cover Dataset (GlobeLand30) Product Description; National Geomatics Center of China: Beijing, China, 2020.
  18. ESA. Chapter 3.1 Statistical Accuracy Assessment. In WorldCover Product Validation Report v1.0; ESA: Paris, France, 2021. [Google Scholar]
  19. Brown, C.F.; Brumby, S.P.; Guzder-Williams, B.; Birch, T.; Hyde, S.B.; Mazzariello, J.; Czerwinski, W.; Pasquarella, V.J.; Haertel, R.; Ilyushchenko, S.; et al. Dynamic World, Near real-time global 10 m land use land cover mapping. Sci. Data 2022, 9, 251. [Google Scholar] [CrossRef]
  20. Karra, K.; Kontgis, C.; Statman-Weil, Z.; Mazzariello, J.C.; Mathis, M.; Brumby, S.P. Global land use/land cover with Sentinel 2 and deep learning. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Brussels, Belgium, 11–16 July 2021; pp. 4704–4707. [Google Scholar]
  21. Venter, Z.S.; Barton, D.N.; Chakraborty, T.; Simensen, T.; Singh, G. Global 10 m Land Use Land Cover Datasets: A Comparison of Dynamic World, World Cover and Esri Land Cover. Remote Sens. 2022, 14, 4101. [Google Scholar] [CrossRef]
  22. Wang, J.; Yang, X.; Wang, Z.; Cheng, H.; Kang, J.; Tang, H.; Li, Y.; Bian, Z.; Bai, Z. Consistency Analysis and Accuracy Assessment of Three Global Ten-Meter Land Cover Products in Rocky Desertification Region—A Case Study of Southwest China. ISPRS Int. J. Geo-Inf. 2022, 11, 202. [Google Scholar] [CrossRef]
  23. Kang, J.; Yang, X.; Wang, Z.; Cheng, H.; Wang, J.; Tang, H.; Li, Y.; Bian, Z.; Bai, Z. Comparison of Three Ten Meter Land Cover Products in a Drought Region: A Case Study in Northwestern China. Land 2022, 11, 427. [Google Scholar] [CrossRef]
  24. Chen, J.; Chen, J.; Liao, A.; Cao, X.; Chen, L.; Chen, X.; He, C.; Han, G.; Peng, S.; Lu, M.; et al. Global land cover mapping at 30 m resolution: A POK-based operational approach. ISPRS J. Photogramm. Remote Sens. 2015, 103, 7–27. [Google Scholar] [CrossRef] [Green Version]
  25. Xie, H.; Wang, F.; Gong, Y.; Tong, X.; Jin, Y.; Zhao, A.; Wei, C.; Zhang, X.; Liao, S. Spatially Balanced Sampling for Validation of GlobeLand30 Using Landscape Pattern-Based Inclusion Probability. Sustainability 2022, 14, 2479. [Google Scholar] [CrossRef]
  26. Troen, I.E.L.P.; Lundtang Petersen, E. European Wind Atlas; Commission of the European Communities: Brussels, Belgium, 1989. [Google Scholar]
  27. Dörenkämper, M.; Olsen, B.T.; Witha, B.; Hahmann, A.N.; Davis, N.N.; Barcons, J.; Ezber, Y.; García-Bustamante, E.; González-Rouco, J.F.; Navarro, J.; et al. The making of the new european wind atlas–part 2: Production and evaluation. Geosci. Model Dev. 2020, 13, 5079–5102. [Google Scholar] [CrossRef]
  28. De Bruin, H.A.; Moore, C.J. Zero-plane displacement and roughness length for tall vegetation, derived from a simple mass conservation hypothesis. Bound.-Layer Meteorol. 1985, 31, 39–49. [Google Scholar] [CrossRef]
  29. Garratt, J.R. The atmospheric boundary layer. Earth-Sci. Rev. 1994, 37, 89–134. [Google Scholar] [CrossRef]
  30. Davenport, A.G.; Grimmond, C.S.; Oke, T.R.; Wieringa, J. Estimating the roughness of cities and sheltered country. In Proceedings of the 12th Conference on Applied Climatology (Asheville Meeting), Asheville, NC, USA, 8–11 May 2000; Volume 96, p. 99. [Google Scholar]
  31. Anderson, J.R. A Land Use and Land Cover Classification System for Use with Remote Sensor Data; US Government Printing Office: Washington, DC, USA, 1976.
  32. Karagali, I.; Mann, J.; Dellwik, E.; Vasiljević, N. New European wind atlas: The Østerild balconies experiment. J. Phys. Conf. Ser. 2018, 1037, 052029. [Google Scholar] [CrossRef] [Green Version]
  33. Mason, P.J. The formation of areally-averaged roughness lengths. Q. J. R. Meteorol. Soc. 1988, 114, 399–420. [Google Scholar] [CrossRef]
  34. Bottema, M. Landscape roughness parameters for Sherwood Forest–Experimental results. Bound.-Layer Meteorol. 1998, 89, 285–316. [Google Scholar] [CrossRef]
  35. OpenWind Online Help. Available online: https://openwind.ul-renewables.com/index.html (accessed on 2 December 2022).
  36. Phillips, G.T.; Traci, R.M. A Preliminary Users Guide for the NOABL Objective Analysis Code; Science Applications, Inc.: La Jolla, CA, USA, 1979. [Google Scholar]
  37. Ainslie, J.F. Calculating the flowfield in the wake of wind turbines. J. Wind Eng. Ind. Aerodyn. 1988, 27, 213–224. [Google Scholar] [CrossRef]
  38. Brower, M.C.; Robinson, N.M. The OpenWind Deep-Array Wake Model: Development and Validation; AWS Truepower: New York, NY, USA, 2012. [Google Scholar]
  39. Emeis, S. Wind Energy Meteorology: Atmospheric Physics for Wind Power Generation; Springer: Berlin/Heidelberg, Germany, 2018. [Google Scholar]
  40. St Martin, C.M.; Lundquist, J.K.; Clifton, A.; Poulos, G.S.; Schreck, S.J. Atmospheric turbulence affects wind turbine nacelle transfer functions. Wind Energy Sci. 2017, 2, 295–306. [Google Scholar] [CrossRef] [Green Version]
  41. IEC 61400-12-2; Wind Turbines-Part 12-2: Power Performance of Electricity Producing Wind Turbines Based on Nacelle Anemometry. International Electrotechnical Commission: Geneva, Switzerland, 2013.
  42. Shin, D.; Ko, K. Application of the nacelle transfer function by a nacelle-mounted light detection and ranging system to wind turbine power performance measurement. Energies 2019, 12, 1087. [Google Scholar] [CrossRef] [Green Version]
  43. IEC 61400-12-1; Wind Turbines-Part 12-1: Power Performance Measurements of Electricity Producing Wind Turbines. International Electrotechnical Commission: Geneva, Switzerland, 2017.
Figure 1. Comparison of remote Angle of View on different Surface coverage data of Wind Farm.
Figure 1. Comparison of remote Angle of View on different Surface coverage data of Wind Farm.
Energies 16 02193 g001
Figure 2. Close-angle comparison of different surface coverage data of wind farms.
Figure 2. Close-angle comparison of different surface coverage data of wind farms.
Energies 16 02193 g002
Figure 3. The roughness length map created by land cover data.
Figure 3. The roughness length map created by land cover data.
Energies 16 02193 g003
Figure 4. Verification flowchart.
Figure 4. Verification flowchart.
Energies 16 02193 g004
Figure 5. Single wake wind speed by Eddy-Viscosity wake model.
Figure 5. Single wake wind speed by Eddy-Viscosity wake model.
Energies 16 02193 g005
Figure 6. Single-wake wind speed by Deep Array EV wake model.
Figure 6. Single-wake wind speed by Deep Array EV wake model.
Energies 16 02193 g006
Figure 7. Single wake wind speed by ASM-EV wake model.
Figure 7. Single wake wind speed by ASM-EV wake model.
Energies 16 02193 g007
Figure 8. Schematic diagram of wind turbine layout.
Figure 8. Schematic diagram of wind turbine layout.
Energies 16 02193 g008
Figure 9. Profile of ridge in northeast region of wind farm.
Figure 9. Profile of ridge in northeast region of wind farm.
Energies 16 02193 g009
Figure 10. Average wind speed and average power during the analysis period.
Figure 10. Average wind speed and average power during the analysis period.
Energies 16 02193 g010
Table 1. Roughness length.
Table 1. Roughness length.
CodeLand Use Land Cover ClassEWA Roughness ClassificationEWA Roughness Length (m)Modified Roughness Length (m)
10Tree cover2.50.650.7
20Shrubland1.50.150.2
30Grassland10.030.05
40Cropland10.050.03
50Built-up310.55
60Bare10.010.01
70Snow and ice00.0010.001
80Water00.00010.0001
90Wetland10.030.02
95Mangroves20.150.15
100Moss and lichen10.010.01
Table 2. Comparison of actual energy and calculated energy.
Table 2. Comparison of actual energy and calculated energy.
FSREVDAWM-EVASM-EVActual
Mean Speed [m/s]5.496 5.265 5.265 5.265 5.291
Gross Energy [GWh]127.644 125.611 125.611 125.611 123.644
Waked Energy [GWh]114.107 105.322 104.097 97.119 96.411
Wake Loss [%]10.6316.13 17.10 22.66 22.03
Net Energy [GWh]91.445 84.405 83.423 77.831 79.057
Table 3. The influence sector of upstream wind turbine on downstream wind turbine.
Table 3. The influence sector of upstream wind turbine on downstream wind turbine.
Down5#6#7#8#
Up
4#103.3°~173.3°115.5°~168.1°123.5°~167.9°130.3°~169.3°
5# 101.7°~177.7°120.3°~176.5°129.3°~174.5°
6# 118.1°~196.9°131.7°~186.7°
7# 123.6°~196.6°
Table 4. Wind speed comparison table.
Table 4. Wind speed comparison table.
4#5#6#7#8#
Actual Waked Speed [m/s]8.403 5.587 4.815 5.294 6.967
EV Waked Speed [m/s]8.379 5.274 5.134 5.177 7.092
DAWM-EV Waked Speed [m/s]8.379 5.274 4.519 4.338 7.034
ASM-EV Waked Speed [m/s]8.379 5.494 4.868 4.822 6.741
Table 5. Simulated turbulence table.
Table 5. Simulated turbulence table.
4#5#6#7#8#
Ambient TI [%]6.66 6.52 6.51 6.27 5.77
EV Total TI [%]6.66 21.16 22.28 21.48 15.68
DAWM-EV Total TI [%]6.66 20.15 25.86 27.03 16.66
ASM-EV Total TI [%]15.80 34.89 43.34 48.38 35.17
Table 6. Array yield comparison.
Table 6. Array yield comparison.
4#5#6#7#8#
Actual Array Yield [kWh]4320.583 1566.864 1011.900 1240.540 2786.203
EV Array Yield [kWh]4619.868 1372.439 1214.295 1482.740 3143.623
DAWM-EV Array Yield [kWh]4619.868 1372.439 806.209 902.434 3071.128
ASM-EV Array Yield [kWh]4070.235 1534.638 1040.528 1155.667 2803.804
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, J.; Chen, J.; Liu, H.; Chen, Y.; Yang, J.; Yuan, Z.; Li, Q. Applicability of WorldCover in Wind Power Engineering: Application Research of Coupled Wake Model Based on Practical Project. Energies 2023, 16, 2193. https://0-doi-org.brum.beds.ac.uk/10.3390/en16052193

AMA Style

Zhang J, Chen J, Liu H, Chen Y, Yang J, Yuan Z, Li Q. Applicability of WorldCover in Wind Power Engineering: Application Research of Coupled Wake Model Based on Practical Project. Energies. 2023; 16(5):2193. https://0-doi-org.brum.beds.ac.uk/10.3390/en16052193

Chicago/Turabian Style

Zhang, Jing, Jixing Chen, Hao Liu, Yining Chen, Jingwen Yang, Zongtao Yuan, and Qingan Li. 2023. "Applicability of WorldCover in Wind Power Engineering: Application Research of Coupled Wake Model Based on Practical Project" Energies 16, no. 5: 2193. https://0-doi-org.brum.beds.ac.uk/10.3390/en16052193

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop