Evaluation of AQ500 sodar performance, Sweden
Evaluation of AQ500 sodar performance, Sweden
Evaluation of AQ500 sodar performance, Sweden
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
<strong>Evaluation</strong> <strong>of</strong> <strong>AQ500</strong> <strong>sodar</strong><br />
<strong>performance</strong>, <strong>Sweden</strong><br />
Report number: KVT/YY/2011/043 – Rev. 1
KVT/YY/2011/R043<br />
Content<br />
1 INTRODUCTION .................................................................................................. 3<br />
2 SITE INFORMATION ............................................................................................. 4<br />
3 SODAR DATA ..................................................................................................... 6<br />
3.1 AVAILABILITY OF STANDARD DEVIATION MEASUREMENTS 6<br />
3.2 AVAILABILITY OF WIND SPEED MEASUREMENTS WITH VARYING WIND SPEED 7<br />
3.3 OPERATION OF THE SODAR 8<br />
4 COMPARISON BETWEEN SODAR AND MAST ................................................................ 9<br />
4.1 FILTERING OF MEASUREMENT DATA 9<br />
4.2 COMPARISON OF MEASURED WIND SPEED AND DIRECTION 10<br />
4.3 COMPARISON OF MEASURED TURBULENCE INTENSITY (TI) 12<br />
4.4 COMPARISON OF MEASURED WIND SHEAR 13<br />
5 SUMMARY AND CONCLUSIONS ............................................................................... 15<br />
6 REFERENCES .................................................................................................... 16<br />
APPENDIX A ................................................................................................................. 17<br />
APPENDIX B ................................................................................................................. 19<br />
2/19
KVT/YY/2011/R043<br />
1 Introduction<br />
In the recent years the use <strong>of</strong> remote sensing equipment for assessing the wind resources when<br />
developing wind energy projects has increased significantly. For Swedish wind farm projects, in<br />
particular, the deployment <strong>of</strong> the <strong>sodar</strong> <strong>AQ500</strong> manufactured by AQ System AS has become<br />
widespread. In this report Kjeller Vindteknikk (KVT) evaluates the <strong>performance</strong> <strong>of</strong> one <strong>AQ500</strong><br />
<strong>sodar</strong> deployed for 6 months (mid-October to mid-April) at a forested site in <strong>Sweden</strong>. The<br />
measurements from the <strong>sodar</strong> are compared to those <strong>of</strong> a nearby 100m meteorological mast<br />
equipped with cup- and ultrasonic anemometers. The work presented in this report is ordered<br />
by Statkraft Development AS.<br />
This analysis <strong>of</strong>fers a comparison <strong>of</strong> the basic wind energy climate parameters measured with<br />
<strong>sodar</strong> and anemometry. The scope <strong>of</strong> this work is to assess the accuracy <strong>of</strong> the data output from<br />
the <strong>AQ500</strong> <strong>sodar</strong> presented in the result files. We have, therefore, not applied any additional<br />
filtering routines to the <strong>sodar</strong> data nor investigated the underlying raw data (that is usually not<br />
accessible for the average <strong>sodar</strong> operator). There are a wide range <strong>of</strong> dependencies and<br />
questions to be answered on the subject <strong>of</strong> <strong>sodar</strong> measurements. With limited time available,<br />
we have prioritised to investigate the relations and parameters we believe to be <strong>of</strong> most<br />
importance for wind resource assessments carried out for similar sites today.<br />
The <strong>sodar</strong> site, henceforward simply called Sodar, is situated in a clear-fell area. The site can<br />
be defined as complex with regard to the forest. However the area where the mast and <strong>sodar</strong><br />
are located is relatively flat. The elevation <strong>of</strong> the area in the vicinity <strong>of</strong> the site is between 490<br />
and 500 m.a.s.l. The height variations in the area around the <strong>sodar</strong> are displayed in Figure 1-1.<br />
Figure 1-1: Map with height variations at the Sodar site and the Mast site.<br />
3/19
KVT/YY/2011/R043<br />
2 Site information<br />
The <strong>sodar</strong> measurement campaign at the Sodar site extends from late October 2010 to mid-<br />
April 2011. Statkraft has carried out measurements in a 100 m lattice mast in the vicinity <strong>of</strong> the<br />
Sodar. There is mature forest <strong>of</strong> about 10-15 m height a couple <strong>of</strong> hundred meters away from<br />
the <strong>sodar</strong> and the mast in several directions. The <strong>sodar</strong> site is a typical Swedish forested area<br />
with a mix <strong>of</strong> clear-fell areas, young forest, and mature forest. With respect to the surface<br />
roughness (forest) the site can be characterized as somewhat complex, however the area is<br />
relatively flat with no hills. Aerial images with an overview <strong>of</strong> the <strong>sodar</strong> site and the exact<br />
position <strong>of</strong> the Sodar and the Mast are shown in Figure 2-1 and Figure 2-2 respectively.<br />
Figure 2-1: Aerial photo <strong>of</strong> the area surrounding the Sodar site. The black frame indicates the borders <strong>of</strong> the cut<br />
out shown in Figure 2-2. The width <strong>of</strong> the black square is approximately 650 m.<br />
Figure 2-2: Aerial photo with the position <strong>of</strong> the Sodar site indicated with a blue dot and the position <strong>of</strong> the<br />
Mast site with a red dot. The distance between the <strong>sodar</strong> and mast is 200 m.<br />
4/19
KVT/YY/2011/R043<br />
The Mast is located 200 meters from the Sodar in the direction 154° (from <strong>sodar</strong> to mast where<br />
north is 0°). The 100 m lattice mast is equipped with cup-anemometers in several heights and a<br />
heated ice-free ultrasonic sensor close to the top mounted cup-anemometers. This makes the<br />
measurements from the mast well suited for comparing the wind speeds measured in 100 m<br />
height, and the wind shear measured by the <strong>sodar</strong> and mast.<br />
The details regarding the measurement period, heights, and placing <strong>of</strong> the <strong>sodar</strong> and mast are<br />
outlined in Table 2-1.<br />
Table 2-1: Details regarding the <strong>sodar</strong> and mast measurements.<br />
Campaign name Meas. heights [m] Campaign length<br />
Site elevation<br />
[m.a.s.l]<br />
Sodar [50,55,...,200] 19.10.2010-18.04.2011 ~495<br />
Mast 100,99, 96, 72, 47 9.4.2010 - present ~496<br />
The technical details regarding the type and settings <strong>of</strong> the <strong>sodar</strong> can be found in Appendix A<br />
together with an instrumentation overview for the Mast.<br />
5/19
KVT/YY/2011/R043<br />
3 Sodar data<br />
The average wind speed measured by the <strong>sodar</strong> at different heights throughout the<br />
measurement campaign is shown in Table 3-1. The underlying dataset is only filtered for<br />
erroneous (9999) values in the wind speed. Other than that the data is used as output from the<br />
<strong>sodar</strong>. The availability <strong>of</strong> the wind speed measurements between 50 and 140 meters is very<br />
satisfying. One should also observe that the <strong>sodar</strong> was deployed during the winter season when<br />
we generally expect the availability <strong>of</strong> measurements campaigns to be lowest.<br />
Table 3-1: Measured wind speed and data availability for various heights at the Sodar site. All valid 10-min wind<br />
speed measurements from the entire <strong>sodar</strong> campaign are included, i.e. only wind speed records with 9999 are<br />
removed.<br />
Meas. Height [m] 50 60 70 80 90 100 110 120 130 140 150 160 180 200<br />
Wind speed [m/s] 6.8 7.4 7.9 8.3 8.8 9.1 9.4 9.7 10.1 10.3 10.5 10.9 11.4 11.9<br />
Availability [%] 98.6 98.6 98.5 98.0 97.6 97.0 96.2 94.9 92.9 90.1 86.7 82.3 73.1 63.2<br />
3.1 Availability <strong>of</strong> standard deviation measurements<br />
The <strong>AQ500</strong> <strong>sodar</strong> can output valid wind speed measurements while the corresponding standard<br />
deviation (STD) measurement is rendered erroneous (9999). The latter is possible because there<br />
is a stricter signal requirement (sound-to-noise ratio) in the <strong>sodar</strong> s<strong>of</strong>tware for standard<br />
deviation values compared to the requirement for a wind speed measurement. The reason for<br />
this is that the <strong>sodar</strong> uses vector average measurements from the entire 10-min period to<br />
calculate the average wind speed, while the standard deviation is calculated every<br />
measurement cycle.<br />
For the average wind speed values presented in Table 3-2 we have removed all wind speed<br />
measurements where the wind speed and/or the corresponding STD measurement are<br />
erroneous. This obviously leads to a lower availability, but it also increases the wind speed with<br />
3.2 % in average over all heights. The increase in average wind speed ranges from 1.7 % in 50 m<br />
height to 4.1 % in 200 m height. In 100 m height the increase is 3.2 %. This result suggests that<br />
the probability <strong>of</strong> the <strong>sodar</strong> to acquire STD measurements is not independent <strong>of</strong> wind speed,<br />
but that there is an increasing probability for loosing STD data with decreasing wind speed. One<br />
should be careful not to exclude wind speed measurements based on the availability <strong>of</strong> the<br />
corresponding STD measurement. The latter makes it, for instance, difficult to calculate the<br />
average wind speed distribution and turbulence intensity from the same dataset without<br />
introducing a bias in one <strong>of</strong> the measures.<br />
Table 3-2: Measured wind speed and data availability for various heights at the Sodar site. Only 10-min wind<br />
records that have both valid wind speed and standard deviation measurements are included, i.e. 10-min records<br />
with either 9999 in the wind speed and/or the standard deviation are excluded.<br />
Meas. height 50 60 70 80 90 100 110 120 130 140 150 160 180 200<br />
Wind speed 6.9 7.5 8.1 8.6 9.0 9.4 9.7 10.1 10.4 10.7 11.0 11.3 11.8 12.5<br />
Availability 85.9 95.5 83.2 80.1 77.3 74.2 71.0 68.0 64.3 60.7 56.8 52.4 43.7 35.3<br />
6/19
KVT/YY/2011/R043<br />
3.2 Availability <strong>of</strong> wind speed measurements with varying wind speed<br />
We have in Figure 3-1 plotted the availability in 200 m as a function <strong>of</strong> the wind speed<br />
measured in 50 m height. The figure indicates how the availability <strong>of</strong> the wind speed<br />
measurement varies with wind speed. The figure shows that failure <strong>of</strong> retrieving valid wind<br />
speed measurements (for this example in 200 m) is more likely to occur at low and high wind<br />
speeds compared to medium wind speeds (6-11 m/s) . The trend is also the same for lower<br />
heights. For lower heights, however, the availability is significantly better than in 200 m and<br />
the problem is therefore <strong>of</strong> less importance. The result shows that there probably is introduced<br />
a bias in the wind speed measurements because the availability depends on wind speed. Note<br />
that we here only consider errors in the mean wind speed measurements, and not the influence<br />
<strong>of</strong> the more frequent errors in the STD measurements.<br />
Figure 3-1: The availability in 200 meters as a function <strong>of</strong> the wind speed measured in 50 m height.<br />
The average wind speeds calculated in 50 m height when concurrent measurements in different<br />
heights are erroneous is shown in Table 3-3. The average wind speed in 50 is lower when there<br />
are errors at higher heights. And correspondingly we see that the average wind speed in 50 m is<br />
higher when only measurements that have valid concurrent measurements at higher heights are<br />
included. The effect <strong>of</strong> this error diminishes with decreasing height since the availability then<br />
improves.<br />
7/19
KVT/YY/2011/R043<br />
Table 3-3: Average wind speed measured in 50 m height for different subsets <strong>of</strong> measurements.<br />
Wind speed in 50 m height based on: Wind speed [m/s] # measurements<br />
all measurements in 50 m : 6.8 25783<br />
measurements when erroneous measurement in 200 m: 6.5 9290<br />
measurements when valid measurement in 200 m: 7.0 16493<br />
measurements when erroneous measurement in 150 m: 6.0 3163<br />
measurements when valid measurement in 150 m: 6.9 22620<br />
measurements when erroneous measurement in 100 m: 4.8 441<br />
measurements when valid measurement in 100 m: 6.8 25342<br />
The crude analysis shown here suggests that the wind speed can be overestimated because the<br />
failure to acquire measurements occurs for climate conditions with low wind speeds. When<br />
computing the wind shear, the problem is to some degree evaded since only concurrent wind<br />
speed measurements from heights the wind shear is computed between are used. If the wind<br />
shear at the site is dependent <strong>of</strong> the wind speed, which is <strong>of</strong>ten the case, on must be careful<br />
not to use data for heights with reduced availability.<br />
3.3 Operation <strong>of</strong> the Sodar<br />
The Sodar was deployed 19.10.2010. About a month later, 16.11.2010, AQ System installed<br />
their cold climate option on the <strong>sodar</strong>. This mainly consists <strong>of</strong> a heating fan for the diesel<br />
generator that is also used for directing hot air on the parabolic dishes in the <strong>sodar</strong> antenna<br />
unit. There have been episodes when the data availability has been low/zero. This is probably<br />
caused by snow on the reflective dishes. The system has, however, autonomously removed the<br />
snow by melting, and returned to normal operation again. The detailed logbook for the Sodar<br />
with an overview <strong>of</strong> the availability each day and all site visits is supplied in Appendix B. All<br />
personnel involved in the maintenance and operation <strong>of</strong> the <strong>sodar</strong> has been instructed to report<br />
their site visits to Kjeller Vindteknikk. We believe that there also has been removed snow from<br />
the unit during the site visits. Unfortunately this preventive action is not specified in the <strong>sodar</strong><br />
logbook, nor has it been reported to us. During the 6 months the campaign lasted the system<br />
has been refilled with diesel 3 times.<br />
The internal clock in the <strong>AQ500</strong> <strong>sodar</strong> has been drifting during the measurement campaign. This<br />
is a known problem with the <strong>AQ500</strong> <strong>sodar</strong> system that we have reported to the manufacturer<br />
several times in the past. Typically the clock on an <strong>AQ500</strong> <strong>sodar</strong> drifts about 5 minutes per<br />
month. We have mended the issue by adjusting the <strong>sodar</strong> time via the manual remote<br />
connection s<strong>of</strong>tware SODWIN.NET whenever the discrepancy between the <strong>sodar</strong> and our server<br />
time is 5 minutes or larger. A more appropriate solution to the problem would have been to<br />
automatically synchronize the <strong>sodar</strong> clock through the automatic s<strong>of</strong>tware SODWIN-MCOM that<br />
normally communicates and downloads data from the <strong>sodar</strong> once every 24 hours. Currently the<br />
latter program does not have this option.<br />
8/19
KVT/YY/2011/R043<br />
4 Comparison between <strong>sodar</strong> and mast<br />
Measurement data from the Mast has been used for comparison with measurements from the<br />
Sodar. The wind speed measurements in 100 m from the <strong>sodar</strong> have been compared with the<br />
top mounted Thies First Class Advanced cup-anemometer at 100.0 m in the mast. The three<br />
side mounted Thies First Class cup-anemometers at heights 96.2 m, 71.6 m, and 47.1 m have<br />
been utilized for comparing the wind shear measured in the mast and by the <strong>sodar</strong>.<br />
4.1 Filtering <strong>of</strong> measurement data<br />
The internal filtering routines developed by AQ system aims at only releasing quality controlled<br />
data in the <strong>sodar</strong> output files and their customers are encouraged to use the data as it appears<br />
and only screen the dataset for obvious errors. The <strong>sodar</strong> data is, therefore, screened and used<br />
as it is reported in the raw files from the <strong>sodar</strong> without further data conditioning. We have only<br />
removed invalid measurements recorded as 9999. We have attempted to use the data both with<br />
and without removing wind speed measurements with corresponding errors in the STD<br />
measurements. Our manual screening <strong>of</strong> the <strong>sodar</strong> data, including the comparison with the<br />
mast data, reveals no obvious errors in the dataset recorded by the <strong>sodar</strong>.<br />
All the cup-anemometer data used in this analysis has been filtered for ice influenced periods<br />
by comparing the measurements to those <strong>of</strong> the Gill heated ultrasonic anemometer mounted in<br />
the mast and the <strong>sodar</strong> data. For the particular campaign at the Mast site the heated ultrasonic<br />
sensor was used for identifying the ice periods through an automatic routine. We emphasise,<br />
however, that ice filtering the cup-anemometers based on deviation from the <strong>sodar</strong><br />
measurements would be possible as the <strong>sodar</strong> and ultrasonic anemometer agrees very well. It is<br />
integral to identify and remove ice influenced measurements from the dataset used for<br />
calculating the long term wind speed. Figure 4-1 shows a typical example from the<br />
measurement time series where there is ice on the cup-anemometer while the agreement<br />
between the ice and snow free <strong>sodar</strong> and ultrasonic anemometer measurements is very good.<br />
The <strong>sodar</strong> and ultrasonic anemometer measurements are both valuable as they verify that ice<br />
and snow do not influence either <strong>of</strong> the instruments.<br />
Figure 4-1: Example <strong>of</strong> an icing episode on the unheated Thies cup anemometer (red). The <strong>sodar</strong> (green) and<br />
the heated sonic anemometer in the mast (blue) continue to measure.<br />
We have filtered out all anemometer measurements when the wind comes from sector 1 – 6<br />
(346°-165°). The latter filtration is carried out to ensure minimal mast influence on the top<br />
mounted cup-anemometer used in the comparison with the <strong>sodar</strong>. The mast has two cup-<br />
9/19
KVT/YY/2011/R043<br />
anemometers mounted in 100 m height, denoted Ano_Thies_A and Ano_Thies_B in the<br />
instrumentation overview supplied in Appendix A. The two anemometers show excellent<br />
correlation (10-min) and have a discrepancy <strong>of</strong> less than 0.1 % for the data used in the<br />
comparison between the mast and <strong>sodar</strong>.<br />
Only concurrent data from the Sodar and the Mast have been used in the comparison presented<br />
in the report. The comparison is carried out based on 10-min averaged wind speed and STD<br />
measurements in 100 m. The direction measurements used in the analysis are recorded in 95 m<br />
and 97.4 m height for the <strong>sodar</strong> and mast respectively.<br />
4.2 Comparison <strong>of</strong> measured wind speed and direction<br />
The measured wind speed and correlation between the anemometers and the <strong>sodar</strong> is shown in<br />
Table 4-1. The availability <strong>of</strong> concurrent measurements after ice and sector filtering is about<br />
50% (or 12969 10-min records) for the 6 months long <strong>sodar</strong> campaign. The table shows that the<br />
average wind speed measured by the <strong>sodar</strong> and in the mast deviate with 0.3 m/s or 3.5 %,<br />
where the anemometers in the Mast display the lower wind speed.<br />
Table 4-1: Mean wind speed measured by the Sodar and in the Mast. All measurements are concurrent,<br />
amounting to 12969 measurements (total availability <strong>of</strong> about 50 %).<br />
Measurement<br />
Meas.<br />
height [m]<br />
Mean wind<br />
speed [m/s] Correlation: Sodar Ano A Ano B<br />
Sodar 100 8.89<br />
Mast Ano_Thies_A 100 8.59 0.98<br />
Mast Ano_thies_B 100 8.60 0.98 1.00<br />
The monthly average wind speeds measured in 100 m height by the Sodar and in the Mast are<br />
shown in Table 4-2. As was the case for the mean <strong>of</strong> all measurements, the discrepancy<br />
between the monthly average wind speeds are also about 0.3 m/s.<br />
Table 4-2: Monthly mean wind speed measured in 100 m height at the Sodar site and the two top mounted<br />
anemometers at the Mast site.<br />
Month Availability (%) Sodar [m/s] Mast Ano A [m/s] Mast Ano B [m/s]<br />
2010 October 33 8.84 8.56 8.56<br />
2010 November 27 8.67 8.36 8.39<br />
2010 December 22 8.26 7.66 7.67<br />
2011 January 38 9.88 9.56 9.61<br />
2011 February 38 8.51 8.19 8.19<br />
2011 March 86 9.24 9.01 9.02<br />
2011 April 52 8.26 7.96 7.94<br />
Mean <strong>of</strong> months: 54 8.81 8.47 8.48<br />
We have also carried out the above comparison for all wind speeds above 2.0 m/s (filtering out<br />
all measurements when either the anemometers or the <strong>sodar</strong> recorded wind speeds below 2.0<br />
m/s). This did not change the outcome, the discrepancy <strong>of</strong> 0.3 m/s persisted.<br />
We do not expect a difference in the wind field <strong>of</strong> this magnitude over the 200 m distance<br />
between the <strong>sodar</strong> and mast position, but some variation can be expected between the sites.<br />
We have preformed a cross-prediction check with the micro scale model WAsP 10 between the<br />
<strong>sodar</strong> and mast site. The model predicts that the average wind speeds at the two locations are<br />
identical. Our qualitative assessment <strong>of</strong> the sites and the WAsP model both indicate that the<br />
difference in measured wind speed cannot be ascribed to the terrain variation alone. We<br />
cannot, however, expect a perfect agreement in the absolute wind speed measured with the<br />
10/19
KVT/YY/2011/R043<br />
the Sodar and the Mast measurements. There are uncertainties in both the calibration <strong>of</strong> the<br />
anemometers (estimated to about 1%) and in the calibration <strong>of</strong> the <strong>sodar</strong>. Further the <strong>sodar</strong><br />
measures in a large volume over the <strong>sodar</strong> site while the anemometer measures in a single<br />
point. We have taken great care to avoid influence by configuration and ice on the<br />
anemometers, but we cannot rule out that a slight bias is still present in the dataset from the<br />
Mast. We have here compared state <strong>of</strong> the art measurements carried out with two different<br />
techniques, and the result shows that one must accept and expect a discrepancy in the<br />
absolute wind speed.<br />
Difference between vector (SODARS) and scalar (LIDARS and anemometers) averaging<br />
It is <strong>of</strong>ten expected that <strong>sodar</strong> data recorded in modest to very complex terrain will yield a<br />
lower wind speed than measurements made by an anemometer in the same wind field. One<br />
reason for this is the difference in the averaging technique between the two instruments<br />
discussed in the following. Note that correcting for the different averaging techniques in this<br />
case will make the discrepancy larger.<br />
The <strong>sodar</strong> computes the 10-minute wind speed based on the 10-minute average radial wind<br />
speed sampled in three directions separated by an azimuth angle <strong>of</strong> 120°, i.e. averaging<br />
vectors. Common practice with other wind speed measuring techniques, such as LIDAR and<br />
anemometer, is to calculate the wind speed (scalar) every time a sample is recorded and then<br />
average this result to form the 10-minute average wind speed. The latter calculation routine is<br />
referred to as scalar average and the former one (which the <strong>AQ500</strong> <strong>sodar</strong> utilize) is called<br />
vector averaging. When the same homogenous wind field is measured, it is theoretical<br />
impossible for the vector average to be larger than the scalar average. They can be equal if the<br />
wind direction in the averaged interval is constant. If the direction however varies, the scalar<br />
mean will be larger than the vector mean. This is why it is commonly expected to measure a<br />
lower wind speed with the <strong>sodar</strong> compared to other measuring techniques. The discrepancy<br />
arising from the different averaging is <strong>of</strong>ten referred to as DP-error (“data processing-error”) in<br />
the literature. The DP-error is dependent on the fluctuations in the wind direction, which in<br />
turn is dependent on the turbulence in the atmosphere. A typical value <strong>of</strong> 2.5% for the DP-error<br />
is reported in (Ioannis Antoniou (ed.) 2003).<br />
In an attempt to reproduce results that are in better agreement with anemometer<br />
measurements one can employ a “vector-to-scalar” transformation on the SODAR data. The<br />
vector-to-scalar transformation <strong>of</strong> the volume measured 10-minute average wind speed, U vector ,<br />
is carried out with the equation,<br />
, where<br />
U scalar is the calculated point wind speed corresponding more closely to what is measured by an<br />
anemometer. The height is denoted by z and the turbulence intensity by TI (10-min standard<br />
deviation <strong>of</strong> wind speed/10-min mean wind speed). Note that TI is calculated from the <strong>sodar</strong><br />
measurements and therefore will contain the bias we are trying to correct for. This equation is<br />
today the best Kjeller Vindteknikk knows for reducing the discrepancies between vector and<br />
scalar averaged measurements. The equation is based on the findings in the article (Ioannis<br />
Antoniou (ed.) 2003). The <strong>sodar</strong> data from the Sodar site that has undergone the above<br />
transformation yield a 0.6 % higher average wind speed than the <strong>sodar</strong> data that has not been<br />
transformed. The transformation from vector averaged to scalar averaged wind speed will<br />
always yield an increase in the wind speed. Correcting for the different averaging techniques<br />
11/19
KVT/YY/2011/R043<br />
for the campaigns in question can only increase the discrepancy in wind speed measured by the<br />
Sodar and the Mast.<br />
The time series with concurrent direction measurements from the Sodar and Mast agrees very<br />
well at 95 m height. The scatter plot in Figure 4-2 shows the agreement between the direction<br />
measurements made in 95 m in the mast and <strong>sodar</strong>.<br />
Figure 4-2: Scatter plot showing the agreement between the concurrent direction measurements in 95 m height<br />
for the <strong>sodar</strong> and in the mast.<br />
4.3 Comparison <strong>of</strong> measured turbulence intensity (TI)<br />
The turbulence conditions at a wind turbine site is normally described through the turbulence<br />
intensity (TI) defined as the 10-min standard deviation divided with the 10-min average wind<br />
speed. The mean and 90 th percentile <strong>of</strong> the TI as a function <strong>of</strong> wind speed measured by the<br />
Sodar and the Mast in 100 m height is presented in Figure 4-3 and Table 4-3. Only<br />
simultaneously recorded measurements are used in the comparison. The figures and table show<br />
that the TI measured by the <strong>sodar</strong> and mast between about 2.5 m/s up to about 13.5 m/s are<br />
very similar. Above this wind speed the TI recorded with the <strong>sodar</strong> and anemometer starts to<br />
deviate. The deviation can possibly be explained by the more frequent occurrence <strong>of</strong> smaller<br />
sized turbulence structures at higher wind speeds. The temporal and the spatial resolution <strong>of</strong><br />
the <strong>sodar</strong> and anemometer measurements are very different, and it is plausible that the<br />
instruments ability to record turbulence varies with the size <strong>of</strong> the turbulence structures. One<br />
should note that there are relatively few concurrent measurements available at the higher wind<br />
speeds. The result shows that at the particular site in question the TI is described very similar<br />
by the two measuring techniques at moderate wind speeds.<br />
12/19
KVT/YY/2011/R043<br />
Figure 4-3: Turbulence intensity with wind speed recorded in 100 m height with the Sodar (left) and Mast<br />
(right).<br />
Table 4-3: The table that shows the measured turbulence intensity in 100 m height measured by the Sodar and<br />
the Mast. Only concurrent measurements from the mast are used in the calculations.<br />
Velocity class [m/s]<br />
Sodar<br />
# observations TI 90 TI mean<br />
Mast<br />
#observations TI 90 TI mean<br />
[2.5 – 3.5> 207 0.21 0.11 223 0.21 0.12<br />
[3.5 – 4.5> 375 0.19 0.11 463 0.20 0.11<br />
[4.5 – 5.5> 501 0.19 0.10 562 0.19 0.10<br />
[5.5 – 6.5> 702 0.16 0.09 757 0.17 0.09<br />
[6.5 – 7.5> 928 0.16 0.09 1077 0.16 0.09<br />
[7.5 – 8.5> 1402 0.15 0.09 1467 0.16 0.09<br />
[8.5 – 9.5> 1686 0.15 0.09 1660 0.16 0.09<br />
[9.5 – 10.5> 1547 0.15 0.10 1587 0.15 0.09<br />
[10.5 – 11.5> 1359 0.15 0.10 1216 0.15 0.09<br />
[11.5 – 12.5> 849 0.15 0.10 701 0.16 0.10<br />
[12.5 – 13.5> 456 0.16 0.11 345 0.17 0.11<br />
[13.5 – 14.5> 254 0.15 0.12 235 0.18 0.13<br />
[14.5 – 15.5> 158 0.15 0.12 142 0.17 0.13<br />
[15.5 – 16.5> 92 0.14 0.12 68 0.18 0.14<br />
[16.5 – 17.5> 33 0.14 0.11 25 0.18 0.14<br />
[17.5 – 18.5> 7 0.12 0.10 9 0.17 0.15<br />
[18.5 – 19.5> 1 0.07 0.07 4 0.16 0.15<br />
4.4 Comparison <strong>of</strong> measured wind shear<br />
The velocity pr<strong>of</strong>ile within a wind farm, i.e. the wind speed change with height, is <strong>of</strong>ten<br />
expressed through the dimensionless wind shear coefficient α. The wind shear coefficient can<br />
be estimated with measurement data through the following expression,<br />
<br />
V<br />
z<br />
V<br />
r<br />
z d <br />
<br />
<br />
zr<br />
d <br />
,<br />
where z is the height above ground level, z r is the reference height, d is the local zero-plane<br />
displacement height, V is the horizontal velocity, and V r is the velocity in the reference height.<br />
13/19
KVT/YY/2011/R043<br />
By rearranging the expression above we can express α through the wind speed measured at<br />
different heights in the mast. We have estimated α using the wind pr<strong>of</strong>ile measured by the<br />
anemometers mounted alongside the Mast and by the Sodar. The wind shear coefficients<br />
calculated between various heights are presented in Table 4-4.<br />
Table 4-4: Comparison <strong>of</strong> measured wind shear by the Sodar and in the Mast. The calculated coefficients are<br />
based on 3487 concurrent 10-min wind speed measurements from the <strong>sodar</strong> and mast.<br />
Wind shear coefficient<br />
measured between:<br />
50m-70m with <strong>sodar</strong><br />
47.1m-71.6m in mast<br />
50m-95m with <strong>sodar</strong><br />
47.1m-96.2m in mast<br />
70m-95m with <strong>sodar</strong><br />
71.6m-96.2m in mast<br />
Sodar 0.40 0.39 0.38<br />
Mast 0.38 0.37 0.35<br />
The agreement between the wind shear measurements made in the mast and by <strong>sodar</strong> is good.<br />
The wind shear measured in the Mast is systematically somewhat lower than the wind shear<br />
measured by the <strong>sodar</strong>. This could be attributed to the differences in the roughness and terrain<br />
around the instruments. In addition, for sites with high wind shear, an overestimation <strong>of</strong> the<br />
wind shear measured is expected because <strong>of</strong> the <strong>sodar</strong> measurement principle. The <strong>sodar</strong><br />
averages the wind speed over a volume with a vertical extent <strong>of</strong> about 20 m. As a result <strong>of</strong> the<br />
high wind shear the wind speed measured by the <strong>sodar</strong> over the volume will be lower than a<br />
point measurement in the middle <strong>of</strong> the volume. This effect will drop <strong>of</strong> with decreasing wind<br />
shear at higher heights, which leads to a higher calculated wind shear coefficient for the <strong>sodar</strong><br />
compared to anemometer measurements in the same heights. An example given by the IEA<br />
Topical Experts Sodar Recommended Practises Group states that a <strong>sodar</strong> will overestimate the<br />
wind shear between 30/50 m with 5 % (Kathleen Moore 2010). This example supposes a volume<br />
average 20 m in depth (same as <strong>AQ500</strong>) and a wind shear coefficient (point measurement) <strong>of</strong><br />
0.40. Another factor that can lead to an overestimation <strong>of</strong> the wind shear measured with the<br />
<strong>sodar</strong> is the effect <strong>of</strong> the DP-error (vector vs. scalar averaging). The effect <strong>of</strong> the DP-error, that<br />
leads to a lower measured wind speed, is reduced with decreasing turbulence and thus with<br />
increasing height.<br />
14/19
KVT/YY/2011/R043<br />
5 Summary and conclusions<br />
We have in this report assessed the <strong>performance</strong> <strong>of</strong> the <strong>AQ500</strong> <strong>sodar</strong> and compared its<br />
measurements with those made with cup-anemometry in a 100 m tall lattice mast located 200<br />
m from the <strong>sodar</strong>.<br />
We have seen that the correlation between the time series from the Sodar and the<br />
anemometers in the Mast is excellent. The strong covariance with the mast, and high<br />
availability <strong>of</strong> the <strong>sodar</strong> measurements, makes the <strong>sodar</strong> suited for identifying ice periods on<br />
the cup-anemometers in the mast with high accuracy.<br />
We have found a discrepancy <strong>of</strong> about 3.5 % in the absolute wind speed measured by the <strong>sodar</strong><br />
and the cup anemometers in the mast. We conclude that a disagreement <strong>of</strong> this order can be<br />
expected between <strong>sodar</strong> and anemometer measurements because <strong>of</strong> the intrinsic uncertainties<br />
and differences in the measuring techniques. We have found that the availability <strong>of</strong> valid wind<br />
speed measurements is correlated with wind speed. This is important to account for when<br />
analysing the data from the <strong>AQ500</strong> <strong>sodar</strong>.<br />
The turbulence intensity (TI) measured with the <strong>sodar</strong> and in the mast agrees well for wind<br />
speeds between about 3 – 13 m/s, and we can conclude that the <strong>sodar</strong> is able to describe the TI<br />
at these wind speeds at the considered site. For wind speeds above this the <strong>sodar</strong> might have<br />
problems resolving the turbulence structures properly. It is uncertain if this actually is the case<br />
since there are few concurrent measurements from the <strong>sodar</strong> and mast at higher wind speeds.<br />
The agreement in the measured wind shear with the two techniques is satisfying, and we<br />
believe that the observed discrepancies can be explained by differences in the measuring<br />
techniques.<br />
We regard the <strong>AQ500</strong> <strong>sodar</strong>, on background <strong>of</strong> our findings, as a valuable contribution to the<br />
ongoing measurement campaign in the Mast. The measurements from the instrument will lead<br />
to a reduction <strong>of</strong> the uncertainty in the predicted wind climate at the site and give valuable<br />
information regarding the wind shear at the site.<br />
15/19
KVT/YY/2011/R043<br />
6 References<br />
Ioannis Antoniou (ed.), Hans E. Jørgensen (ed.) et al. On the Theory <strong>of</strong> SODAR Measurement.<br />
Final reporting on WP1, EU WISE project NNE5-2001-297, Risø-R-1410(EN): Risø National<br />
Laboratory, Roskilde, Denmark, 2003.<br />
Kathleen Moore, et Al. Recomended Practices for the Use <strong>of</strong> Sodar in Wind Energy Resource<br />
Assesment. Expert recomendations (draft), Roskilde: IEA Topical Expert Group, 2010.<br />
16/19
KVT/YY/2011/R043<br />
Appendix A<br />
Table 0-1: Deployment and system details for the Sodar site<br />
Deployment details<br />
Installation Crew :<br />
Max Gråhed (AQ Systems), Peter Rothmund (KVT)<br />
Location coordinates and elevation :<br />
RT90: N XX, E XX. ~497masl<br />
System orientation (North = 0°) : 334° (measured 330° + magnetic declination 4°)<br />
System details<br />
System type<br />
<strong>AQ500</strong>C<br />
System name<br />
AQS3<br />
Trailer reg. number, system serial number<br />
XX<br />
Sodar/control unit phone number<br />
XX<br />
Alarm number<br />
+47 XX<br />
SODAR s<strong>of</strong>tware version<br />
RE01<br />
Measurement height 50-200m<br />
Control unit s<strong>of</strong>tware version 4.00<br />
Clock settings (measurements)<br />
UTC (GMT+0)<br />
17/19
KVT/YY/2011/R043<br />
Table 0-2: Instrumentation overview for the Mast site.<br />
18/19
KVT/YY/2011/R043<br />
Appendix B<br />
Log book for the Sodar site<br />
Site: XXXX AQS3 Location: Sverige, XXXX Customer: Statkraft Development<br />
* ok . Limited data Quality - No data High wind or snow is "ok" data, but can be low quality<br />
Date/Sign. Year Month 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31<br />
Sodar mounted 19.October 2010 by AQ systems / Max Gråhed and KVT / Peter<br />
2010 October * * * * * * * * . * . * *<br />
2010 Nevember * * * * * * * . * . . * * * * . * * * * . . * * * * * * * *<br />
2010 December * * . * * . * . * * . * . * * * . * * * * * . * * * * . * . *<br />
2011 January * . . * * * . * * . * * . . * * * * * * * * . . . . * * * * *<br />
2011 February * * * * * * * * * * * * . * . * * . . * * . * * * * * *<br />
2011 March * * * * * * * * * * * * * * * * * * * * * . . * * * * * . . .<br />
2011 April * * * * . * . . * *<br />
Site visits and notes<br />
19.10.2010 YY Sodar deployed<br />
16.11.2010 10 nov installert dieselvarmere av Max.<br />
29.11.2010 AJ Message sms Low fuel level 27.11.2010<br />
06.12.2010 AJ 1.12.2010 filled up 160 l fuel. 40 l reserve.<br />
13.12.2010 AJ Bad spectrum and small amount <strong>of</strong> data. Snow on membrans?<br />
15.12.2010 AJ Looking good again.See other side logg. It seems that the snow has been melted by the heat produced in the agg unit.<br />
27.12.2010 AJ Sms Low fuel level 23.12.2010. SMS Fredrik regard. assist. 27.12. Filled up 29.12.2010 with 150 l <strong>sodar</strong> and app 70 l res.<br />
05.01.2011 AJ AQ system at location 28.12.2010. Shut down between 19:00 and 24:00.<br />
26.01.2011 AJ Changed <strong>sodar</strong> time from 9:30 to 9:35<br />
22.02.2011 AJ SMS Low Fuel Level yesterday. Informed Fredrik. He will take move out today or Wednesday<br />
24.02.2011 AJ 23.02: Fredrik filled up 154 liters <strong>sodar</strong> and 70 l agg.<br />
11.03.2011 AJ Changed <strong>sodar</strong> time from 15:33 to 15:39<br />
22.03.2011 YY Changed <strong>sodar</strong> time from 13:53 to 13:55<br />
02.05.2011 AJ Dismounted Sodar 20.04.2011<br />
19/19