17.11.2014 Views

Evaluation of AQ500 sodar performance, Sweden

Evaluation of AQ500 sodar performance, Sweden

Evaluation of AQ500 sodar performance, Sweden

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!