ENERGY 3D Nacelle Mounted Lidar in Complex Terrain PCWG Hamburg, Germany Paul Lawson 25.03.2015 1 DNV GL 125.03.2015 SAFER, SMARTER, GREENER
Agenda Introduction and Project Background Lidar Specifications Wind Speed Derivation Met Mast and Lidar Comparisons Lidar Uncertainty from Terrain Concluding Remarks 2
Project Motivation Demonstrate the capability of 3D nacelle mounted lidar to produce an accurate power curve in complex terrain Potentials Reduce installation time, campaign duration, and overall cost More measurement points across the rotor area to better estimate power performance Track the wind with the yaw position Perform site calibration of sites where it is not feasible or possible Multiple measurement distances could provide an alternative to site calibration Challenges Must make several assumptions to resolve horizontal wind Complex terrain results in variable beam heights above ground Measures a volume rather than a fixed point An accurate power curve assessment revolves around an accurate site calibration 3
Project Background Campaign Details Previous site calibration DNVGL 2009 in accordance with IEC standards Site complexity characterization IEC 61400-12-1 Cold climate, winter in north-eastern United States Duration of current campaign: November 5, 2014 February 14, 2015 Two valid 10 sectors: 265-275 275-285 Met Mast Anemometry Met mast located at 240m and 265 from turbine 4
Terrain 5
Lidar Specifications Cooperative Project Lidar Specifications Manufacturer Model Type Beam configuration Measurement ranges [m] Outputs Measurement Width Refresh Rate Operational Availability Avent Prototype Pulsed 4 beams with rectangular arrangement. Horizontal separation of 30 Vertical separation of 10 0.5D, 0.8D, 1.3D, 1.7D, 2.1D, 2.5D, 2.9D, 3.3D, 3.8D, 4.2D Horizontal wind speed, wind direction, shear and veer exponents, REWS, TI ±14.5m 1 second No recorded downtime due to environmental effects Beta test for prototype lidar based off Wind Iris platform 6
0.85D 0.8D Nacelle Lidar in Complex Terrain: Beam Positioning Top View Side View 1.63D 0.5D 0.8D 1.3D 1.7D 2.1D 2.5D 0.5D 0.8D 1.3D 1.7D 2.1D 2.5D First 6/10 measurement distances shown for the top and side views Met mast located at 2.5D The lidar beam exits the projected swept area [---] at 1.63D 7
Lidar Hub Height Horizontal Wind Speeds Two horizontal wind speeds are derived at different heights (from the upper and lower pairs of beams) The power law can be incorporated to reconstruct a wind profile and interpolate for hub height wind speeds H Ref 1 H Ref 2 2.5D V Hub Height = V Ref 1,2 H Hub Height H Ref 1,2 ln V Ref 1 V Ref 2 ln H Ref 1 H Ref 2 8
Site Calibration and Wind Speed Interpolation Wind speed reconstitution is the same process, only turbine conditions (RPM speed) differ 9
What is a Met Mast Site Calibration? A ratio of wind speeds used to describe far away wind speeds and transpose them to the turbine position Previous site calibration completed in 2009 in accordance with IEC standards Chose future turbine installation Install Mast A and Mast B Measure hub height wind speeds Calculate R=V A /V B Remove Mast A for turbine installation Transpose future V B by multiplying by R 10
What is a Lidar Site Calibration? A ratio of wind speeds used to describe far away wind speeds and transpose them to the turbine position V 0.5D V 2.5D Nacelle lidar installation Turbine idling or RPM <1 Derive upper and lower wind speeds at 0.5D and 2.5D Interpolate height wind speeds at 0.5D and 2.5D Calculate R=V 0.5D /V 2.5D Transpose future V 2.5D by multiplying by R 11
Wind Speed Correlations Lidar 2.5D [m/s] 20 18 16 14 12 10 8 6 4 2 Hub Height Wind Speeds from Lidar 2.5D and Met Mast Wind Speeds Linear Fit Line 0 0 2 4 6 8 10 12 14 16 18 20 Met Mast Anemometer [m/s] Important Filters Minimum 3.5m/s for lidar and met mast No Maximum wind speed Wind shear 0.0<α<0.5 Availability >0.40 Bin Center Slope R 2 270 0.996 0.986 280 0.996 0.988 The compared wind speeds here are the raw measurements before any site calibration factors have been applied 12
HH HH Wind Wind Lidar Lidar Speeds Speeds Normalized to to 240m 240m Wind Wind Speeds Speeds Site Calibration Turbine OFF 1.02 1.02 1 1 0.98 0.98 Elevation of Terrain Range Gate Width Normalized Hub Height Wind Speeds data4 data5 data6 Terrain Terrain Effect Effect Ratio Ratio for Yaw for Position Yaw Position 265-285 265-285 0.96 0 50 100 150 200 250 300 350 400 610 0.96 Distance from Turbine 0 0.5D 50 1.0D 100 1.6D 150 2.1D 200 2.6D 250 3.1D 300 3.6D 350 4.2D 400 610 [m] Distance from Turbine 630 625 620 615 Terrain Elevation Above Sea Level [m] 630 625 620 615 Terrain Elevation Above Sea Level [m] Turbine RPM must be off or idling 265 <yaw<285 Wind direction 265 <WD<285 Min wind 4 m/s Max wind 16 m/s Lidar hub height wind speeds normalized to those measured at 2.5D Red lines are the measured range gates ±( 14.5m) Comparison to IEC approved site calibration from 2009 (pre-turbine installation) Bin Center IEC - 2009 Lidar- V 0.5D /V 2.5D 270 0.980 0.970 280 0.983 0.964 13
Cal Factor i / Cal Factor mean Convergence of the Terrain Correction Factor 1.03 1.02 1.01 1 0.99 0.98 Convergence of Moving Average Site Cal. Factor onto Average Sector Center 270 Sector Center 280 Total 265-285 Turbine RPM must be <1 265 <yaw<285 Turbine not intentionally shut down 0.97 0 2 4 6 Time 8 10 12 [hours] The moving average of V 50m /V 240m normalized to the final, overall average outlined Initial stages of convergence (analysis) of site calibration factor in accordance with IEC 61400-12-1, Annex C, Section C.3 After 7.5 hours (45 samples) the 270 centred sector is trending at 0.9999, 0.9999, 1.000, 1.0004 Does not yet meet full requirements of 24hrs of data per 10 sector 14
Power Curves Comparison Electrical Power (kw) 2400 1600 800 Electrical Power (kw) 3000 2000 1000 0.6 0.4 0.2 C p Reference Met Mast Lidar data4 0 0 2 4 6 8 10 0 12 14 16 18 20 0 Binned Wind Speeds [m/s] 0 2 4 6 8 10 12 14 Binned Wind Speeds [m/s] Respective site calibration factors have been applied Lidar power curve is shifted to the left as a result of lower measured wind speeds and lower site calibration factors Close proximity of lidar power curves, suitable for real world application! 15
Shear Exponent Comparison Lidar 2.5D 0.5 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 Linear Regression of Shear Exponent Slope of 0.849 for a linear line forced through 0 R 2 value of 0.5497 10 minute averages 0 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 Anemometer A1 Modelling of the wind profile using different reference elevations Lidar shear was calculated with the upper and lower beam pairs while the met mast used two cup anemometers at 80.2m and 51.7m above ground level Lidar measures above and below hub height - better representation of the profile 16
Start of Explanation for Lidar-Mast Shear Differences Lower wind shear above hub height Higher wind shear below hub height Lidar measures above and below hub height - better representation of the profile The cup anemometer measures wind speeds at and below hub height Different measurement points along the wind profile (i.e. mast and lidar) can result in wind shear exponent differences 17
Shear Exponent Comparison 0.35 Binned Shear Exponent of the Lidar Lidar 240m 0.3 Lidar 240m 0.25 0.2 0.15 0.1 2 4 6 8 10 12 Binned Wind Speed 14 16 18 [m/s] Wind shear values binned into 0.5m/s bins Very low sample count in the 4m/s bin, not sufficient to be statistically representative (a total of 9 samples) 18
Turbulence Intensity 0.22 0.2 Anemometer D1 50.7m Anemometer A1 80.2m Lower Lidar Beams ~62m Upper Lidar Beams ~104m Turbulence Intensity for Binned Wind Speeds Binned wind speed data (0.5m/s bins) Turbulence Intensity 0.18 0.16 0.14 0.12 0.1 2 4 6 8 10 12 14 16 18 Binned Wind Speed [m/s] A good approximation exists and the expected trend of decreasing TI as height increases is true. Following of the extreme peaks and drops is a promising result Lidar TI is calculated from the RWS rather than the actual horizontal wind speed 19
Uncertainty Lidar Specific Challenges Power law is incorporated, how do uncertainties in beam elevation translate to uncertainties in wind speed? Four main sources are identified in the elevation uncertainty: Dynamic tilting of the turbine during operation causing the height beams to change (tilt of turbine binned to the 0.5m/s bins) In calculating the average height of the terrain for the left and the right beams Inherent uncertainty of the lidar position due to uncertainty in the lidar s internal inclinometer Elevation uncertainty from the terrain source and roughness Not the final solutions but a first step! 20
Uncertainty due to Terrain Lidar 2.5D -Anemometer 0.6 0.5 0.4 0.3 0.2 0.1 0-0.1-0.2-0.3-0.4-0.5 translated elevation Difference in Binned Wind Speeds Count of Samples per Bin 100 90 80 70 60 50 40 30 20 10 Number of Samples Only uncertainty in wind speed due to uncertainty in height of beams is displayed (red) Site specific! -0.6 2 4 6 8 10 12 14 16 18 0 Wind Speed [m/s] The four main sources of elevation uncertainty have been translated from measurement elevation above ground (meters) to wind speed (m/s) There exists other lidar uncertainties which are not shown These uncertainties would remain within the bounds of cup anemometer uncertainty 21
Blockage Effect HH Wind Lidar Speeds Normalized to 240m Wind Speeds 1.01 1 0.99 0.98 0.97 0.96 Blockage Effect for Yaw Position 265-285 Range Gate Width Normalized Hub Height Wind Speeds Turbine RPM > 8 Min wind speed 3.5m/s 0.95 0 0.5D 1.0D 1.6D 2.1D 2.6D 3.1D 3.6D 4.2D 450 Distance from Turbine [m] Lidar hub height wind speeds normalized to those measured at 2.5D Site calibration factors have been applied Beams exit the swept area at 1.63D, no blockage at 2.1D which is less than 2.5D Lidar beams exit the swept area of the turbine here 22
Deviations? V 1.9D The compression zone is the lidar measurement range(s) which are encompassed with the swept area During terrain effect calculations we pretend this affected area doesn t exist. V 2.5D -0.5 y = - 0.00267*x - 0.693 Tilt During Turbine OFF Linear Fit Line Tilt During Turbine ON There is tilting of the turbine when during (terrain effect) measurements. Passive blockage? Tilt Angle [deg] -1-1.5 0 2 4 6 8 10 Lidar 240m 12 14 16 18 20 [m/s] 23
Conclusions Positive outcomes Beta test in complex terrain show a high correlation between cup anemometer and lidar measured wind speeds at turbine hub height Successful wind speed retrieval at hub height using power law interpolation method and turbine tilting compensation technique. Encouraging results in shear and TI High convergence on lidar measured site calibration factors with a measurable difference from met mast site calibration factors No lidar downtime due to environmental conditions Future Plans Trial at a site with high terrain slopes to attempt vertical wind speed estimations Investigate rotor equivalent wind speed measurements (REWS) and TI renormalization Utilize the lidar measurements as the turbine yaw for additional sector site calibrations 24
Special Thanks A collaborative and international project with promising initial results could not have been done with out all the support from my advisors at: Avent Lidar Technology DNVGL Forwind SunEdison 25
Thank you for your attention! Paul Lawson paul.lawson@dnvgl.com OR paul.lawson.eng@gmail.com +49(0)1511 6536623 www.dnvgl.com SAFER, SMARTER, GREENER Luke Simmons luke.simmons@dnvgl.com Erik Tüxen Erik.Tuexen@dnvgl.com Vineet Parkhe Vineet.Parkhe@dnvgl.com Matthieu Boquet mboquet@leosphere.com Florian Rebeyrat frebeyrat@leosphere-avent.com 26