# Wind Characteristics of Three Meteorological Stations in China

Authors: Yang Yang; Yao Gang; Wang Rong; Wang Hengyu
Journal: Journal of Control Science and Engineering
 (2015)

---

## Abstract

With rapid economic development of China, demand for energy is growing rapidly. Many experts have begun to pay attention on exploiting wind energy. Wind characteristics of three meteorological stations in China were analyzed to find out if or not it is possible to build a wind farm in this paper. First of all, studies about the wind characteristics and potential wind energy were summarized. Then ways of collecting and manipulating wind data were introduced. Wind-generation potential was assessed by the method of Weibull distribution. Wind shear exponent, extreme wind speed in 50 years, and turbulence intensity were calculated. The wind characteristics were summarized and assessment of wind-generation potential was given. At last, the wind was simulated with autoregressive method by Matlab software.


---

## Body

## 1. Introduction

With rapid economic development of China, demand for energy is growing rapidly [1]. Electricity produced by consuming coal had taken up of 76% in the whole electricity. This led to the spring of 2013 when the fog hazy weather lasted for a long time in China.China is one of the countries which have a large number of wind resources and make use of the wind resource in the early days. According to statistics [2], density of wind energy in China is about 100W/m2 on average, and the total wind energy amount is 3226GW. The available development and utilization wind energy is about 253GW on land. The wind energy utilization has been growing at a fast pace around all over the world. But the utilization is held back by the shortage of wind data which are sufficient and credible. The wind speeds’ distribution is not only significant to the wind farm designers but also to other power generators. Nowadays, vast quantities of information is offered to researchers by research of wind characteristic, investors, and wind farm planners that pay close attention to the renewable energy’s evolution.Gansu province which is the fifth on the potential wind energy in China has 237 million KW potential wind energy in theory [3]. Jiuquan region in Gansu is rich of wind energy, and its annual average wind density is above 150W/m2. Nowadays, energy policy of China is committed to developing a kind of new energy which is abundant and free pollution with low cost. So considering the exploitation of wind energy resources in Jiuquan is required.In the last five years, lots of researches had been done in China on potential wind energy and wind characteristics. Song et al. [4] investigated major measures of observation data and computational precision based on field measurement and analyses of a number of wind power facilities, in particular over complex terrain. Du et al. [5, 6] proposed an improved method which had the advantage of less uncertainty and less error in computing the wind speed of the representative year in wind farm and presented changes of the wind shear index basing on the measured data of three 70 meter-high wind towers of two consecutive years in Inner Mongolia. Peng et al. [7] studied the wind shear exponent for wind resource assessment and found that the wind shear exponent in different vertical layers calculated at different speed bins had the highest precision. Chen et al. [8] studied the wind characteristics in northwest and southeast of China by analyzing the observational data obtained from towers which were representative in the two areas. Fang et al. [9] analyzed the wind power in northeast of China. The method and model for wind energy resource investigation in cold area had been established in his study. Wen et al. [10] analyzed the reserves and distribution characteristics of wind speed and wind energy at 70m height with statistical method, basing on observational data at 70m height from 18 wind towers in Fujian coastal areas. Tong [11] put forward a new method of wind speed forecasting-wavelet decomposition autoregressive combined with wind speed change parting forecast method by bringing wind speed change characteristics types into wind speed prediction. Shang and Xie [12] developed logistic models according to the development characteristics of each province to predict development trend of wind power in 2013~2020, and predicting the main problems and measures should be taken by different provinces in the future. Bassyouni et al. [13] collected eleven years’ wind data in Jeddah to determine the wind characteristics including the daily, monthly, and annual wind speed, wind probability density distribution, and shape (k) and scale (c) parameters at 10m height. Alam et al. [14] studied the variability of fluctuating nature of the wind both in time and spatial domain. By utilizing daily mean values of wind speed from different meteorological stations, dynamic nature of the wind at nine stations was analyzed in conjunction with wavelet transform and fast Fourier transform power spectrum techniques. Wu et al. [15] developed a prediction system to predict wind power with a method of combining statistical model and physical model. The inlet condition of the wind farm is forecasted by the autoregressive model in this system. The proposed prediction system was tested by the data from Wattle Point Wind Farm in Australia, and this system was effective for power output prediction of wind farm. Santamaría-Bonfil et al. [16] proposed a hybrid methodology based on Support Vector Regression for wind speed forecasting with autoregressive model, and results showed that forecasts made with this method were more accurate for medium short term wind power forecasting and wind speed forecasting than other models.China has sufficient potential wind energy especially in the northwest region. But development of wind resource is under the expected standard until now. This paper aims to make sure of the wind’s characteristics and simulate the wind with autoregressive method.


## 2. Wind Speed Data Analysis

### 2.1. Description of Wind Data Collection System

Guazhou meteorological station which belonged to national principal station was built in 1951 with a 1171m altitude above sea level. The wind speed data from 1981 to 2010 was collected at Guazhou meteorological station. Measurement data was kept with EL and EL15-2/2A anemometer.Wind farm was about 67km away the northeast of Guazhou meteorological station, which had a similar geographic and geomorphic condition with Guazhou meteorological station. Therefore, quantitative data of Guazhou meteorological station could be used in analyzing the characteristics of wind.The annual average wind speed and wind direction during 1981–2010 year in Guazhou was given in Figure1. It can be calculated that 2.55m/s was the average speed in 30 years. Average wind speed kept stability and east wind direction was the prevailing wind direction. The month average wind speed during 1981–2010 year in Guazhou was given in Figure 2. The highest average speed was 3.4m/s in March and the smallest was 2.0m/s in October. It evidenced that wind in Guazhou kept stability in each month. Therefore it is possible to build a wind farm in this region.Figure 1
Annual average wind speed and wind direction during 1981–2010 year in Guazhou meteorological station.
(a)
Wind speed
(b)
Wind directionFigure 2
Month average wind speed during 1981–2010 year in Guazhou meteorological station.More than fifty meteorological stations with 70m height had been built in Jiuquan region to develop the wind recourse. There were three meteorological stations in the survey region, and hence all of them were selected to analyze the wind characteristics of the wind farm. The numbers of meteorological stations were 8158#, 7411#, and 8601#. Table 1 gave the basic condition of the three meteorological stations.Table 1
Site specific information of meteorological stations considered in this paper.
Number
H (m)
Period
Coordinate
Elevation (m)
Wind tower configuration
Instrument device
8158#
70
2008.09.01–2011.04.28
N: 40°48′39′′ E: 96°31′20′′
1566
Wind speed: 10, 30, 50, 70 direction: 10, 70
NRG
7411#
70
2008.09.01–2011.04.27
N: 40°44′28′′ E: 96°31′18′′
1494
Wind speed: 10, 30, 50, 70 direction: 10, 70
NRG
8601#
70
2008.08.25–2011.04.28
N: 40°45′33′′ E: 96°24′46′′
1500
Wind speed: 10, 30, 50, 70 direction: 10, 70
NRG


### 2.2. Manipulation of Wind Data

In order to get the assessment which was most close to the actual wind power generation potential in the wind farm, the collected wind data were separately checked by all of wind data and the wind data which were measured in 2009 year. NRG anemoscopes were used in the meteorological stations, and sampling frequency of wind speed was 1Hz. Data of wind speed and wind direction was output each 10min. Sampling frequency of extreme wind speed was 2Hz. Collection of wind direction kept the same pace with the wind speed, and all of wind direction records were moment sampling values. The anemoscopes could store the wind data collected in three months completely, and the staff got the wind data every two months. According to [17], the wind data were checked with data integrity, data scope, data dependency, and trend of wind speed. The unreasonable and missing data were deleted. At last, all the unreasonable data were collected and the rationality of data was judged again. The valid data were picked out and returned to the original data set. Table 2 showed that the entire rate of valid data was beyond 80% during observation period.Table 2
Valid data’s rate of three meteorological stations.
Time
2008.9.01–2011.4.28
2008.9.01–2008.12.31
2009.1.01–2009.12.31
2010.1.01–2010.12.31
2011.1.01–2011.4.28
8158#
93.64%
95.16%
95.83%
89.76%
97.31%
7411#
96.47%
97.20%
96.71%
95.28%
98.68%
8601#
89.32%
95.33%
80.09%
95.59%
96.74%In order to check the correlation between wind speed and wind direction, wind data collected from different meteorological stations were selected. Correlation between wind speed and wind direction was analyzed between 8601# and 8158#, 8601# and 7411#, and 8158# and 7411# meteorological stations. Table3 showed the related equation and correlation coefficient between 8601# and 8158#, 8601# and 7411#, and 8158# and 7411# at 70m height in 2009 year. The wind data between 8601# and 7411# was taken as an example, the total correlation coefficient was 0.859, and the related equation was Y=0.869X+0.623. The 8601# and 7411# mainly wind correlation coefficient was above 0.75. The wind speed had a good correlation coefficient in three meteorological stations. It can be seen that wind power generation potential in the wind farm was almost the same.Table 3
Correlation coefficient of three meteorological stations.
(a)
Correlation coefficient of 8601# and 8158# meteorological stations at 70 m height
The total correlation equation:Y = 0.884X + 0.556
Correlation coefficient: 0.859
Quadrant
Correlation equation
Correlation coefficient
N
Y = 0.605X + 1.824
0.627
NNE
Y = 0.776X + 0.414
0.879
NE
Y = 0.906X + 0.132
0.8
ENE
Y = 0.956X + 0.099
0.843
E
Y = 0.862X + 1.130
0.824
ESE
Y = 0.797X + 1.485
0.746
SE
Y = 0.787X + 1.096
0.8
SSE
Y = 0.605X + 1.370
0.664
S
Y = 0.532X + 1.369
0.603
SSW
Y = 0.481X + 1.371
0.553
SW
Y = 0.714X + 0.907
0.676
WSW
Y = 0.869X + 0.512
0.829
W
Y = 0.893X + 0.365
0.871
WNW
Y = 0.860X + 0.328
0.855
NW
Y = 0.834X + 0.497
0.846
NNW
Y= 0.840X + 1.090
0.772
(b)
Correlation coefficient of 8601# and 7411# meteorological stations at 70 m height
The total correlation equation:Y = 0.869X + 0.623
Correlation coefficient: 0.859
Quadrant
Correlation equation
Correlation coefficient
N
Y = 0.607X + 2.020
0.635
NNE
Y = 0.796X + 0.487
0.882
NE
Y = 0.820X + 0.731
0.75
ENE
Y = 0.917X + 0.173
0.819
E
Y = 0.876X + 0.623
0.798
ESE
Y = 0.785X + 0.967
0.766
SE
Y = 0.840X + 0.623
0.849
SSE
Y = 0.717X + 0.946
0.769
S
Y = 0.665X + 1.041
0.721
SSW
Y = 0.576X + 1.219
0.626
SW
Y = 0.800X + 0.719
0.798
WSW
Y = 0.902X + 0.441
0.885
W
Y = 0.918X + 0.429
0.898
WNW
Y = 0.910X + 0.561
0.876
NW
Y = 0.860X + 0.902
0.846
NNW
Y = 0.807X + 1.597
0.772
(c)
Correlation coefficient of 8158# and 7411# meteorological stations at 70 m height
The total correlation equation:Y = 0.918X + 0.682
Correlation coefficient: 0.932
Quadrant
Correlation equation
Correlation coefficient
N
Y = 0.786X + 0.671
0.802
NNE
Y = 0.959X + 0.114
0.964
NE
Y = 1.000X + 0.207
0.933
ENE
Y = 0.974X + 0.575
0.949
E
Y = 0.907X + 1.112
0.943
ESE
Y = 0.958X + 1.081
0.926
SE
Y = 0.903X + 0.800
0.875
SSE
Y = 0.743X + 0.983
0.735
S
Y = 0.669X + 0.878
0.695
SSW
Y = 0.573X + 1.129
0.642
SW
Y = 0.746X + 0.687
0.737
WSW
Y = 0.886X + 0.514
0.858
W
Y = 0.905X + 0.482
0.909
WNW
Y = 0.841X + 0.658
0.873
NW
Y = 0.878X + 0.377
0.898
NNW
Y = 0.921X + 0.261
0.875From the checked results of wind data, the 8601#, 8158#, and 7411# meteorological stations can meet the requirement of design with more than 95% effective data. The monthly average wind speed data in 2009 collected from three meteorological stations at different heights were compared in Figure3. It showed that wind speed data in 2009 collected from three meteorological stations at the same height had a good agreement, and hence wind speed data in 2009 had a certain representation. From Figure 3, it also can be seen that the wind speed was concentrated mainly on 5m/s–9m/s. Invalid wind speed under 3m/s was rare, and there was no destructive wind all over the year. With the small change of the wind speed, this wind farm can generate electricity through the whole year.Figure 3
Wind speed of three meteorological stations at different heights in 2009 year.


## 2.1. Description of Wind Data Collection System

Guazhou meteorological station which belonged to national principal station was built in 1951 with a 1171m altitude above sea level. The wind speed data from 1981 to 2010 was collected at Guazhou meteorological station. Measurement data was kept with EL and EL15-2/2A anemometer.Wind farm was about 67km away the northeast of Guazhou meteorological station, which had a similar geographic and geomorphic condition with Guazhou meteorological station. Therefore, quantitative data of Guazhou meteorological station could be used in analyzing the characteristics of wind.The annual average wind speed and wind direction during 1981–2010 year in Guazhou was given in Figure1. It can be calculated that 2.55m/s was the average speed in 30 years. Average wind speed kept stability and east wind direction was the prevailing wind direction. The month average wind speed during 1981–2010 year in Guazhou was given in Figure 2. The highest average speed was 3.4m/s in March and the smallest was 2.0m/s in October. It evidenced that wind in Guazhou kept stability in each month. Therefore it is possible to build a wind farm in this region.Figure 1
Annual average wind speed and wind direction during 1981–2010 year in Guazhou meteorological station.
(a)
Wind speed
(b)
Wind directionFigure 2
Month average wind speed during 1981–2010 year in Guazhou meteorological station.More than fifty meteorological stations with 70m height had been built in Jiuquan region to develop the wind recourse. There were three meteorological stations in the survey region, and hence all of them were selected to analyze the wind characteristics of the wind farm. The numbers of meteorological stations were 8158#, 7411#, and 8601#. Table 1 gave the basic condition of the three meteorological stations.Table 1
Site specific information of meteorological stations considered in this paper.
Number
H (m)
Period
Coordinate
Elevation (m)
Wind tower configuration
Instrument device
8158#
70
2008.09.01–2011.04.28
N: 40°48′39′′ E: 96°31′20′′
1566
Wind speed: 10, 30, 50, 70 direction: 10, 70
NRG
7411#
70
2008.09.01–2011.04.27
N: 40°44′28′′ E: 96°31′18′′
1494
Wind speed: 10, 30, 50, 70 direction: 10, 70
NRG
8601#
70
2008.08.25–2011.04.28
N: 40°45′33′′ E: 96°24′46′′
1500
Wind speed: 10, 30, 50, 70 direction: 10, 70
NRG


## 2.2. Manipulation of Wind Data

In order to get the assessment which was most close to the actual wind power generation potential in the wind farm, the collected wind data were separately checked by all of wind data and the wind data which were measured in 2009 year. NRG anemoscopes were used in the meteorological stations, and sampling frequency of wind speed was 1Hz. Data of wind speed and wind direction was output each 10min. Sampling frequency of extreme wind speed was 2Hz. Collection of wind direction kept the same pace with the wind speed, and all of wind direction records were moment sampling values. The anemoscopes could store the wind data collected in three months completely, and the staff got the wind data every two months. According to [17], the wind data were checked with data integrity, data scope, data dependency, and trend of wind speed. The unreasonable and missing data were deleted. At last, all the unreasonable data were collected and the rationality of data was judged again. The valid data were picked out and returned to the original data set. Table 2 showed that the entire rate of valid data was beyond 80% during observation period.Table 2
Valid data’s rate of three meteorological stations.
Time
2008.9.01–2011.4.28
2008.9.01–2008.12.31
2009.1.01–2009.12.31
2010.1.01–2010.12.31
2011.1.01–2011.4.28
8158#
93.64%
95.16%
95.83%
89.76%
97.31%
7411#
96.47%
97.20%
96.71%
95.28%
98.68%
8601#
89.32%
95.33%
80.09%
95.59%
96.74%In order to check the correlation between wind speed and wind direction, wind data collected from different meteorological stations were selected. Correlation between wind speed and wind direction was analyzed between 8601# and 8158#, 8601# and 7411#, and 8158# and 7411# meteorological stations. Table3 showed the related equation and correlation coefficient between 8601# and 8158#, 8601# and 7411#, and 8158# and 7411# at 70m height in 2009 year. The wind data between 8601# and 7411# was taken as an example, the total correlation coefficient was 0.859, and the related equation was Y=0.869X+0.623. The 8601# and 7411# mainly wind correlation coefficient was above 0.75. The wind speed had a good correlation coefficient in three meteorological stations. It can be seen that wind power generation potential in the wind farm was almost the same.Table 3
Correlation coefficient of three meteorological stations.
(a)
Correlation coefficient of 8601# and 8158# meteorological stations at 70 m height
The total correlation equation:Y = 0.884X + 0.556
Correlation coefficient: 0.859
Quadrant
Correlation equation
Correlation coefficient
N
Y = 0.605X + 1.824
0.627
NNE
Y = 0.776X + 0.414
0.879
NE
Y = 0.906X + 0.132
0.8
ENE
Y = 0.956X + 0.099
0.843
E
Y = 0.862X + 1.130
0.824
ESE
Y = 0.797X + 1.485
0.746
SE
Y = 0.787X + 1.096
0.8
SSE
Y = 0.605X + 1.370
0.664
S
Y = 0.532X + 1.369
0.603
SSW
Y = 0.481X + 1.371
0.553
SW
Y = 0.714X + 0.907
0.676
WSW
Y = 0.869X + 0.512
0.829
W
Y = 0.893X + 0.365
0.871
WNW
Y = 0.860X + 0.328
0.855
NW
Y = 0.834X + 0.497
0.846
NNW
Y= 0.840X + 1.090
0.772
(b)
Correlation coefficient of 8601# and 7411# meteorological stations at 70 m height
The total correlation equation:Y = 0.869X + 0.623
Correlation coefficient: 0.859
Quadrant
Correlation equation
Correlation coefficient
N
Y = 0.607X + 2.020
0.635
NNE
Y = 0.796X + 0.487
0.882
NE
Y = 0.820X + 0.731
0.75
ENE
Y = 0.917X + 0.173
0.819
E
Y = 0.876X + 0.623
0.798
ESE
Y = 0.785X + 0.967
0.766
SE
Y = 0.840X + 0.623
0.849
SSE
Y = 0.717X + 0.946
0.769
S
Y = 0.665X + 1.041
0.721
SSW
Y = 0.576X + 1.219
0.626
SW
Y = 0.800X + 0.719
0.798
WSW
Y = 0.902X + 0.441
0.885
W
Y = 0.918X + 0.429
0.898
WNW
Y = 0.910X + 0.561
0.876
NW
Y = 0.860X + 0.902
0.846
NNW
Y = 0.807X + 1.597
0.772
(c)
Correlation coefficient of 8158# and 7411# meteorological stations at 70 m height
The total correlation equation:Y = 0.918X + 0.682
Correlation coefficient: 0.932
Quadrant
Correlation equation
Correlation coefficient
N
Y = 0.786X + 0.671
0.802
NNE
Y = 0.959X + 0.114
0.964
NE
Y = 1.000X + 0.207
0.933
ENE
Y = 0.974X + 0.575
0.949
E
Y = 0.907X + 1.112
0.943
ESE
Y = 0.958X + 1.081
0.926
SE
Y = 0.903X + 0.800
0.875
SSE
Y = 0.743X + 0.983
0.735
S
Y = 0.669X + 0.878
0.695
SSW
Y = 0.573X + 1.129
0.642
SW
Y = 0.746X + 0.687
0.737
WSW
Y = 0.886X + 0.514
0.858
W
Y = 0.905X + 0.482
0.909
WNW
Y = 0.841X + 0.658
0.873
NW
Y = 0.878X + 0.377
0.898
NNW
Y = 0.921X + 0.261
0.875From the checked results of wind data, the 8601#, 8158#, and 7411# meteorological stations can meet the requirement of design with more than 95% effective data. The monthly average wind speed data in 2009 collected from three meteorological stations at different heights were compared in Figure3. It showed that wind speed data in 2009 collected from three meteorological stations at the same height had a good agreement, and hence wind speed data in 2009 had a certain representation. From Figure 3, it also can be seen that the wind speed was concentrated mainly on 5m/s–9m/s. Invalid wind speed under 3m/s was rare, and there was no destructive wind all over the year. With the small change of the wind speed, this wind farm can generate electricity through the whole year.Figure 3
Wind speed of three meteorological stations at different heights in 2009 year.


## 3. Wind Characteristic Analysis

### 3.1. Weibull Parameter of Wind Speed

In the studied areas, frequency distribution of wind speed is a principal element to evaluate potential wind power. Potential wind power in this wind farm must be determined by the Weibull distribution. The range of wind data is wide. As a consequence, keeping a few important parameters which could certify the behavior of wind speed data with a wide range is requisite. Using a distribution function is the easiest method for the procedure, and the most widely used function is Weibull function. The two-parameter distribution Weibull function could be showed in what follows:(1)fV=kcVck-1×e-V/ck,where k is a dimensionless parameter of Weibull shape, c is a Weibull scale parameter, and V is the wind speed, m/s. For the sake of assessing k and c, different kinds of ways have been come up with. Wind power can be assessed with the following:(2)W=12ρv3,where W is the wind power density, ρ is average air density, and v is the average speed at 10m height. Air density is calculated by(3)ρ=1.2761+0.00366t×p-0.378e1000,where t is temperature on average, p is the air pressure on average, and e is the average vapor pressure. The air density in the wind farm is 1.045kg/m3.The Weibull parameters were determined by using the software WASP10.0 in this paper. Table4 showed the results of two Weibull parameters k and c at different height and different meteorological stations. Weibull parameters k and c ranged from 1.85 to 2.16, 6.4 to 9.1 separately. Both of k and c increased with the height under 50m and kept nearly steady above the 50m height. The wind speed frequency distribution Weibull fit curves at different height were described in order to check the calculate results. Table 5 showed density of wind power and wind speed at different height. Wind speeds in three meteorological stations were almost the same and had an error of no more than 4.6% at the same height. Average power density of wind ranged from 374W/m2 to 503W/m2 when the height of meteorological station changed from 50m to 90m height and had an error of no more than 9.9% at the same height. According to the standard [18], the grade of this wind farm was beyond three, and wind energy resources were abundant and steady.Table 4
The Weibull parameters at different heights.
Number
Parameter
10 m
30 m
50 m
65 m
70 m
75 m
80 m
90 m
8158#
c
6.4
7.9
8.6
8.8
8.8
8.9
9.0
9.1
k
1.85
2.2
2.24
2.16
2.16
2.16
2.16
2.16
7411#
c
6.6
7.8
8.5
8.8
8.9
9.0
9.1
9.2
k
1.92
2.14
2.21
2.15
2.15
2.16
2.16
2.16
8601#
c
6.5
7.8
8.3
8.4
8.5
8.5
8.6
8.8
k
1.86
2.16
2.15
2.05
2.05
2.05
2.05
2.06Table 5
Power density of wind and wind speed of three meteorological stations.
Number
Average wind speed
(m/s)
Average power density of
wind (W/m2)
50 m
75 m
90 m
50 m
75 m
90 m
8158#
7.60
7.9
8.08
396
453
487
7411#
7.55
7.98
8.16
389
465
503
8601#
7.38
7.61
7.79
374
422
453


### 3.2. Calculation of Wind Shear Exponent

Wind shear exponent is a key factor on wind turbine safety in wind resource assessment. Wind shear exponent in different vertical layers calculated in different speed bins has the highest precision in all methods, while the wind shear exponent calculated by the exponent law has the characteristic of high credibility and high stability. Thus exponent law was used in this paper. Wind shear exponent can be calculated by(4)a=lg⁡v2/v1lg⁡z2/z1,where v2 is wind speed at z2 height and m/s, v1 is speed at z1 height, m/s.Table6 showed the results of wind shear at different heights. When the height is above the 70m height, the wind shear exponent is constant and independent of height. Therefore the wind shear exponent was assigned a value of 0.13 when the height was above 70m height.Table 6
Results of wind shear at different heights.
Number
8158#
7411#
8601#
H (m)
10
30
50
70
10
30
50
70
10
30
50
70
α
0.146
0.147
0.132
0.131
0.138
0.139
0.139
0.130
0.142
0.142
0.135
0.129


### 3.3. Extreme Wind Speed in 50 Years

Wind speed and turbulence intensity determine the grade of wind turbine generator system. The wind farm can be divided into three grades by extreme wind speed in 50 years, associating with three grades of wind turbine generator system. Extreme wind speed in 50 years determines the security and rationality of the wind farm. Table7 showed the grade of wind farm.Table 7
Security grade of the wind farm.
IEC grade
I
II
III
Extreme wind speed in 50 years (m/s)
<70
<59.5
<52.5Type I extreme value distribution method shown in (5) was used to assess the extreme wind speed in 50 years:(5)V50-max=u-1ωln⁡ln⁡250249,where u is distribution location parameter, u=17.04. ω is distribution scale parameter, and ω=0.4892.Air density in this wind farm is 1.045kg/m3. Considering the change of air density, the extreme wind speed in 50 years should be modified with the condition of standard air density. Figure 4 showed the extreme wind speed in 50 years should be modified with the condition of standard air density.Figure 4
Comparison of wind speed between standard and actual air density.


### 3.4. Turbulence Intensity at Different Heights

Turbulence intensity function at the speed of 15m/s could be calculated from(6)IT=δV,where V is the average speed between 14.5m/s and 15.5m/s; σ is corresponding wind speed standard deviation. Table 8 showed the turbulence intensity at different heights. Turbulence intensity was decreased with the increasing of height. Range of turbulence intensity variation is 0.108 to 0.059. The C standard wind farm turbulence intensity is 0.12, so the turbulence intensity in this wind farm can meet the requirement [18].Table 8
Turbulence intensity at different heights.
No.
70 m
50 m
30 m
10 m
8158#
0.059
0.073
0.093
0.109
7411#
0.059
0.074
0.093
0.109
8601#
0.062
0.076
0.089
0.108


## 3.1. Weibull Parameter of Wind Speed

In the studied areas, frequency distribution of wind speed is a principal element to evaluate potential wind power. Potential wind power in this wind farm must be determined by the Weibull distribution. The range of wind data is wide. As a consequence, keeping a few important parameters which could certify the behavior of wind speed data with a wide range is requisite. Using a distribution function is the easiest method for the procedure, and the most widely used function is Weibull function. The two-parameter distribution Weibull function could be showed in what follows:(1)fV=kcVck-1×e-V/ck,where k is a dimensionless parameter of Weibull shape, c is a Weibull scale parameter, and V is the wind speed, m/s. For the sake of assessing k and c, different kinds of ways have been come up with. Wind power can be assessed with the following:(2)W=12ρv3,where W is the wind power density, ρ is average air density, and v is the average speed at 10m height. Air density is calculated by(3)ρ=1.2761+0.00366t×p-0.378e1000,where t is temperature on average, p is the air pressure on average, and e is the average vapor pressure. The air density in the wind farm is 1.045kg/m3.The Weibull parameters were determined by using the software WASP10.0 in this paper. Table4 showed the results of two Weibull parameters k and c at different height and different meteorological stations. Weibull parameters k and c ranged from 1.85 to 2.16, 6.4 to 9.1 separately. Both of k and c increased with the height under 50m and kept nearly steady above the 50m height. The wind speed frequency distribution Weibull fit curves at different height were described in order to check the calculate results. Table 5 showed density of wind power and wind speed at different height. Wind speeds in three meteorological stations were almost the same and had an error of no more than 4.6% at the same height. Average power density of wind ranged from 374W/m2 to 503W/m2 when the height of meteorological station changed from 50m to 90m height and had an error of no more than 9.9% at the same height. According to the standard [18], the grade of this wind farm was beyond three, and wind energy resources were abundant and steady.Table 4
The Weibull parameters at different heights.
Number
Parameter
10 m
30 m
50 m
65 m
70 m
75 m
80 m
90 m
8158#
c
6.4
7.9
8.6
8.8
8.8
8.9
9.0
9.1
k
1.85
2.2
2.24
2.16
2.16
2.16
2.16
2.16
7411#
c
6.6
7.8
8.5
8.8
8.9
9.0
9.1
9.2
k
1.92
2.14
2.21
2.15
2.15
2.16
2.16
2.16
8601#
c
6.5
7.8
8.3
8.4
8.5
8.5
8.6
8.8
k
1.86
2.16
2.15
2.05
2.05
2.05
2.05
2.06Table 5
Power density of wind and wind speed of three meteorological stations.
Number
Average wind speed
(m/s)
Average power density of
wind (W/m2)
50 m
75 m
90 m
50 m
75 m
90 m
8158#
7.60
7.9
8.08
396
453
487
7411#
7.55
7.98
8.16
389
465
503
8601#
7.38
7.61
7.79
374
422
453


## 3.2. Calculation of Wind Shear Exponent

Wind shear exponent is a key factor on wind turbine safety in wind resource assessment. Wind shear exponent in different vertical layers calculated in different speed bins has the highest precision in all methods, while the wind shear exponent calculated by the exponent law has the characteristic of high credibility and high stability. Thus exponent law was used in this paper. Wind shear exponent can be calculated by(4)a=lg⁡v2/v1lg⁡z2/z1,where v2 is wind speed at z2 height and m/s, v1 is speed at z1 height, m/s.Table6 showed the results of wind shear at different heights. When the height is above the 70m height, the wind shear exponent is constant and independent of height. Therefore the wind shear exponent was assigned a value of 0.13 when the height was above 70m height.Table 6
Results of wind shear at different heights.
Number
8158#
7411#
8601#
H (m)
10
30
50
70
10
30
50
70
10
30
50
70
α
0.146
0.147
0.132
0.131
0.138
0.139
0.139
0.130
0.142
0.142
0.135
0.129


## 3.3. Extreme Wind Speed in 50 Years

Wind speed and turbulence intensity determine the grade of wind turbine generator system. The wind farm can be divided into three grades by extreme wind speed in 50 years, associating with three grades of wind turbine generator system. Extreme wind speed in 50 years determines the security and rationality of the wind farm. Table7 showed the grade of wind farm.Table 7
Security grade of the wind farm.
IEC grade
I
II
III
Extreme wind speed in 50 years (m/s)
<70
<59.5
<52.5Type I extreme value distribution method shown in (5) was used to assess the extreme wind speed in 50 years:(5)V50-max=u-1ωln⁡ln⁡250249,where u is distribution location parameter, u=17.04. ω is distribution scale parameter, and ω=0.4892.Air density in this wind farm is 1.045kg/m3. Considering the change of air density, the extreme wind speed in 50 years should be modified with the condition of standard air density. Figure 4 showed the extreme wind speed in 50 years should be modified with the condition of standard air density.Figure 4
Comparison of wind speed between standard and actual air density.


## 3.4. Turbulence Intensity at Different Heights

Turbulence intensity function at the speed of 15m/s could be calculated from(6)IT=δV,where V is the average speed between 14.5m/s and 15.5m/s; σ is corresponding wind speed standard deviation. Table 8 showed the turbulence intensity at different heights. Turbulence intensity was decreased with the increasing of height. Range of turbulence intensity variation is 0.108 to 0.059. The C standard wind farm turbulence intensity is 0.12, so the turbulence intensity in this wind farm can meet the requirement [18].Table 8
Turbulence intensity at different heights.
No.
70 m
50 m
30 m
10 m
8158#
0.059
0.073
0.093
0.109
7411#
0.059
0.074
0.093
0.109
8601#
0.062
0.076
0.089
0.108


## 4. Simulation of Wind Speed with Autoregressive Method

Simulation of wind speed is particularly important for wind farms because of cost-related issues, dispatch planning, and energy markets operations. Typically, wind farm energy production is estimated using a fixed weighted measure of the wind farm’s nominal power and forecasts from historical atmospheric data. Further, it has been stated that wind speed is one of most important variables related to wind power generation. Autoregressive model is commonly used for time series forecasting wind speed since they are able to capture persistence in a time series. Therefore, wind speed was simulated with autoregressive method in this paper. Principle of autoregressive model can be expressed from(7)Vt=v-z+vt,v1j·tv2j·t⋮vMj·t=-∑k=1pψk·v1j-k·tv2j-k·t⋮vMj-k·t+N1j·tN2j·t⋮NMj·t,j·t=0,1,2,3,…,T,k≤j,v-z=v-10z10α,where V(t) is simulation wind, v-(z) is average wind, v(t) is turbulent wind, and ψk is regression coefficient. N(t) is a Gaussian function, p is the order of autoregressive mode, α can be found in Table 6, z is the height of wind, and v-(10) is the average wind at the height of 10m.Simulation of wind speed with autoregressive method was realized by Matlab software. Time step was defined as 0.1, calculation step was defined as 1024, and operation order was 4. Wind in April 2009 at 70m height was simulated according to davenport spectrum with autoregressive method. Figures 5 and 6 gave the turbulent and average wind in 400 seconds. From Figure 5, it can be seen that turbulent wind ranged from −1m/s to 2m/s at 70m height, and it fluctuated around 0m/s. From Figure 6, it can be seen that average wind ranged from 6.5m/s to 9.5m/s at 70m height, and it fluctuated around 7.4m/s which kept a good agreement with the measured wind speed data.Figure 5
Turbulent wind in 400 seconds.Figure 6
Average wind in 400 seconds.Simulated power spectrum was compared with the standard one in Figure7, the simulated power spectrum agreed well with the standard power spectrum. With the data collected, turbulent wind and average wind at different heights can be simulated accurately with autoregressive method.Figure 7
Simulated and standard power spectrum.


## 5. Conclusions

(1)
East wind is the prevailing wind direction in this wind farm. The wind speed is strong in winter and spring. By analyzing the rate of valid data and comprising wind speed data between meteorological stations, wind data collected in 2009 were used to assess wind-generation potential. The wind frequency mainly concentrated on 5m/s–9m/s. Invalid wind speed under 3m/s was rare, and there was no destructive wind. With the small change of the wind speed, this wind farm can generate electricity through the whole year.
(2)
Average wind power density at different heights was calculated. Average wind speed at 50m height was 7.51m/s and the density of wind power on average was 387W/m2. The grade of wind farm was beyond three.
(3)
Under the condition of standard air density the maximum wind speed in 50 years is 36.4m/s and 37.3m/s at the height of 75m and 90m. Both of them are less than 37.5m/s which is the reference wind speed in third grade of wind farm. And hence it is possible to build a wind farm of third grade in this region.
(4)
Range of turbulence intensity variation is 0.108 to 0.059. TheC standard wind farm turbulence intensity is 0.12, so the turbulence intensity in this wind farm can meet the requirement.
(5)
With data collected and davenport spectrum, both wind speed at different heights and wind speed in this wind farm can be simulated accurately with autoregressive method.



---

*Source: 325435-2015-10-22.xml*