## 1. Introduction

Lidar systems are able to provide information about the wind field approaching a wind turbine in advance, which can be used to assist wind turbine control. This field of investigations has increased significantly in recent years, and several controllers for load reduction or energy yield increase have been tested in simulations (see, e.g., Laks et al. 2010; Dunne et al. 2012; Schlipf et al. 2013b; Koerber and King 2011; Henriksen 2011; Kragh et al. 2013). The latest work also considers wind evolution (see Bossanyi 2012; Simley et al. 2012; Laks et al. 2013).

Initial field testing of lidar-assisted feedforward collective pitch control shows that improvement of the control performance can be confirmed in reality and proves that the performance of lidar-assisted control is dependent on the correlation between the lidar system and the turbine (see Schlipf et al. 2012a; Scholbrock et al. 2013).

In this work, a new analytic model of correlation between the rotor effective wind speed of the turbine and the lidar estimate is derived based on models of turbulent wind, lidar measurement, and wind evolution. Previous work is based on simulated or real measurements to the knowledge of the authors (see Schlipf et al. 2010a; Simley et al. 2013). Here, techniques from signal processing are used to build up a modular and extensible framework. This model can be used in the planning of a lidar-assisted control project to determine whether a lidar system is able to provide useful signals or to optimize the configuration or scanning pattern of a lidar system. Furthermore, the model was used successfully in field testing to evaluate the lidar measurements and to design an adaptive filter, the core of the used controller.

This paper is organized as follows. In section 2 the requirements for the feedforward collective pitch control are explained. In sections 3 and 4, the model for a perfect staring and a real scanning lidar is derived, which is evaluated with data from the field testing in section 5. The model is used in section 6 to optimize a lidar system, and conclusions and future work are discussed in section 7.

## 2. Requirements for control

Reducing fatigue and extreme loads of the structure is an important design goal for the control of large wind turbines. Transient events such as gusts represent an unknown disturbance to the control system. Conventional feedback controllers can only provide delayed compensation for such excitations, since the disturbance effects must propagate through the entire wind turbine before showing its effects in the measured outputs. This usually results in additional loads for the wind turbine and requires high actuator rates for the disturbance compensation. These effects can be avoided if the wind ahead of the wind turbine is measured by remote sensing techniques such as lidar and the information is fed to the turbine controller.

The magnitude of load reduction depends on the quality of the wind preview, expressed by the correlation of the lidar measurements and the turbine reaction. In this section the requirements of the correlation for lidar-assisted control are derived by a description of a feedforward collective pitch controller. The same approach can be used for the direct speed controller to optimize the energy yield (see Schlipf et al. 2013a).

The feedforward controller (see Fig. 1) is based on the work in Schlipf et al. (2010a) and combines the baseline feedback controller with a feedforward update. The main control goal of the collective pitch feedback controller Σ_{FB} is to maintain the rated rotor speed Ω_{rated}. The system Σ is disturbed by a wind field _{L} in front of the turbine before reaching the rotor. If the wind does not change on its way (Σ_{E} = 1), then in the case of perfect measurement the measured wind speed *υ*_{0L} and the rotor effective wind speed *υ*_{0} are equal. The disturbance could be perfectly compensated by a feedforward controller

*υ*

_{0}cannot be measured perfectly because of the wind evolution Σ

_{E}and the limitation of the lidar system Σ

_{L}. Therefore, the needed feedforward controller is

_{E}is very complicated. Also, not all information of V can be reconstructed by the inverse of Σ

_{L}. However, if the transfer function

*G*

_{RL}from the measured wind speed to the rotor effective wind speed is known, then it can be used to exploit all information captured by the lidar system:

*G*

_{RL}can be obtained from measurements and approximated by a standard low-pass filter. Therefore, the cutoff frequency (−3 dB) of the corresponding filter can be considered as a quality criterion for the correlation. In the following sections, an analytic way to estimate both the transfer function

*S*

_{LL}from the lidar and

*S*

_{RR}from the rotor, as well as the cross-spectrum

*S*

_{RL}between the rotor and the lidar.

## 3. Correlation of a perfect staring lidar system

*S*

_{LL},

*S*

_{RR}, and

*S*

_{RL}for a perfect staring lidar system is done fully and semianalytically. They are derived based on Kaimal wind spectra, but in principle other wind spectra models can be used. The spectra

*S*

_{ii}_{,h}of the longitudinal (

*h*= 1), lateral (

*h*= 2), and upward (

*h*= 3) velocity component in point

*i*is given in the International Electrotechnical Commission (IEC) standard 61400-1 (IEC 2005) by the equation:

*f*is the frequency;

*σ*and

_{h}*L*are the standard deviation and the integral-scale parameter of the corresponding velocity component, respectively; and

_{h}*r*is defined as

_{ij}*κ*is the frequency-dependent lateral decay parameter and

*S*

_{ij}_{,1}is the cross-spectrum between points

*i*and

*j*, and assuming that the rotor plane is normal to the average wind direction.

*u*

_{HH}at the hub:

*A*with rotor radius

*R*:

*x*,

*y*, and

*z*axes are pointing in the downwind, lateral, and upward directions, respectively. With (6) the explicit solution of (9) can be found as

*I*

_{0}is the modified Bessel function of the first kind. Furthermore, it is assumed that the autospectra of the three wind components in all points are equal to the spectra

*S*

_{HH,h}at the hub.

*S*

_{Hj}

_{,1}and

*γ*

_{Hj}

_{,1}as the cross-spectrum and coherence of the longitudinal wind component between the hub and the point

*j*, respectively, (11) is solved by

*G*

_{RL}and the squared coherence

*k*, independent of the mean wind speed

*n*longitudinal wind components

*u*hitting the rotor plane:

_{i}*G*

_{RL}and

*n*.

## 4. Correlation of a real scanning lidar system

For a real scanning lidar system, the following effects are integrated into the semianalytic model: discrete scanning, wind evolution, spatial averaging, and wind reconstruction. These effects will make the calculation of *υ*_{0L} more complex, but the main idea of the approach can still be used: The measured wind can be considered as a sum of signals and because of the linearity of the Fourier transformation, the spectra can be calculated by a sum of auto- and cross-spectra.

All effects except wind evolution are verified isolated and in combination with a lidar simulator from Schlipf et al. (2009) and a TurbSim wind field simulator (Jonkman and Buhl 2007). The model agrees with the correlation of the rotor effective wind speed from the simulated lidar and the wind field, if the wind field is only evaluated at the grid points, because interpolation in the grid overestimates the correlation.

In the following subsections the effects will be explained separately, giving a simple example. The spectrum of a real scanning lidar system is a combination of all these effects.

### a. Discrete scanning

*T*is represented by

### b. Wind evolution

*d*is given by

_{ij}*α*is the dimensionless longitudinal decay parameter. Other models such as Kristensen (1979) can also be used. If, for example, a second point downwind is added to the perfect staring lidar system (7), then the rotor effective wind speed estimate can be calculated by shifting the measurement downwind in time considering Taylor's frozen turbulence hypothesis:

### c. Wind reconstruction

*i*can be modeled as

*l*,

_{xi}*l*,

_{yi}*l*] is the normalized backscattered laser beam vector and [

_{zi}*u*,

_{i}*υ*,

_{i}*w*] is the wind vector. If the lidar is perfectly aligned with the wind, then the estimated lateral and vertical wind components are assumed to be zero, and the reconstructed longitudinal component

_{i}*l*= 0), the rotor effective wind speed is

_{zi}### d. Spatial averaging

*f*(see Lindelöw 2008). For pulsed systems it can be assumed that

_{L}*f*depends only on the distance

_{L}*a*to the focus point. For continuous wave systems, there is an additional dependency on the distance to the lidar. With the weighting function

*f*(

_{L}*a*), it is possible to calculate the line-of-sight wind speed of each focus point

*i*by

*f*} of the weighting function, depending on the wavenumber

_{L}*k*. But here a discretized approach is used consistent with the overall semianalytic approach. For this purpose, the weighting function is evaluated at discrete distances and the spatial separations are treated similar to the time delays of a temporal filter.

*f*

_{L1}and

*f*

_{L2}at

*a*

_{1}= −

*a*and

*a*

_{2}=

*a*, then the measurement can be modeled by

*i*is the imaginary unit.

## 5. Evaluation with real data

The analytic and the semianalytic model are both evaluated with real data from the field testing of lidar-assisted control. Here, a scanning lidar system from the University of Stuttgart was installed at the National Wind Technology Center on the nacelle of the two-bladed Controls Advanced Research Turbine (CART2) with a rotor diameter of *D* = 42.7 m and a hub height of 36.8 m.

The signals necessary for this correlation study are obtained from lidar and turbine data; see Schlipf et al. (2012a) for more details. The rotor effective wind speed *υ*_{0} is estimated from measured turbine data by a nonlinear estimator. For the lidar system, a circular trajectory with *n*_{fp} = 6 focus points in *n*_{fd} = 5 focus distances equally distributed between 1*D* and 2*D* diameters and a scan time of *T* = 1.33 s were used. Figure 2 shows both signals from the considered 5-min period together with the wind speed measured with the anemometer on the nacelle. The lidar provides a preview signal that coincides with the wind speed estimate from the turbine for low frequencies. The anemometer measures a lower wind speed because of its position behind the rotor, which extracts energy from the wind. With the comparison to the analytic and semianalytic models, it can be confirmed that the lack of correlation for high frequencies is caused by the combination of the effects described in section 4.

Wind speed from the nacelle anemometer (dark gray) and rotor effective wind speed estimates from the CART2 (light gray) and the scanning lidar system (black).

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

Wind speed from the nacelle anemometer (dark gray) and rotor effective wind speed estimates from the CART2 (light gray) and the scanning lidar system (black).

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

Wind speed from the nacelle anemometer (dark gray) and rotor effective wind speed estimates from the CART2 (light gray) and the scanning lidar system (black).

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

### a. Evaluation of the analytic model

The analytic model cannot be evaluated with a real lidar system because wind evolution and volume measurement cannot be neglected. Therefore, the anemometer signal is compared to rotor effective wind speed. Figure 3 shows that the transfer function based on (10) and (12) fits to the data. Confidence bounds are calculated based on Bendat and Piersol (1971). Furthermore, the transfer function based on the semianalytic model (16) with *n* = 357 equally over the rotor area distributed points with 2-m separation is close to the one of the fully analytic model. The separation distance is chosen as a compromise between accuracy and computational effort.

Transfer function of the nacelle anemometer and the CART2: analytic model (black), coinciding semianalytic model (light gray, dotted), measured (dark gray), and ±*σ* confidence bounds (light gray shading). Maximum wavenumber

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

Transfer function of the nacelle anemometer and the CART2: analytic model (black), coinciding semianalytic model (light gray, dotted), measured (dark gray), and ±*σ* confidence bounds (light gray shading). Maximum wavenumber

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

Transfer function of the nacelle anemometer and the CART2: analytic model (black), coinciding semianalytic model (light gray, dotted), measured (dark gray), and ±*σ* confidence bounds (light gray shading). Maximum wavenumber

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

### b. Evaluation of the semianalytic model

*f*(

_{L}*a*) with a full width at half maximum (FWHM) of

*W*= 30 m is used, following the considerations of Lindelöw (2008) and Banakh and Smalikho (1997):

*n*

_{wp}= 5 equally distributed points with a distance of 10 m. In Simley et al. (2013) the volume averaging is considered as a source of error when comparing to a single-point measurement and beneficial when comparing to a blade effective wind speed. The weighting function also is beneficial for the determination of the rotor effective wind speed, because more information is collected over the rotor disk. Furthermore, the dimensionless longitudinal decay parameter

*α*= 0.4 is roughly estimated with a staring lidar system measuring from a wind turbine based on Schlipf et al. (2010b). For the model

*υ*

_{0}is obtained by

*n*= 357 equally over the rotor area distributed points with 2-m separation. The modeled

*υ*

_{0L}is calculated by a sum of 3

*n*

_{fp}

*n*

_{fd}

*n*

_{wp}= 450 signals of the

*u*,

*υ*, and

*w*wind components.

Figure 4 shows good agreement with the data. Based on this model, a filter depending on the mean wind speed

Transfer function of the scanning lidar system and the CART2: analytic (black), measured (dark gray), and ±*σ* confidence bounds (light gray shading). Maximum wavenumber

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

Transfer function of the scanning lidar system and the CART2: analytic (black), measured (dark gray), and ±*σ* confidence bounds (light gray shading). Maximum wavenumber

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

Transfer function of the scanning lidar system and the CART2: analytic (black), measured (dark gray), and ±*σ* confidence bounds (light gray shading). Maximum wavenumber

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

## 6. Lidar system optimization

The proposed model can be used to optimize the configuration or the scanning pattern of a lidar system. As an example the model is used to determine the optimum correlation of a lidar system with three independent beams on a turbine with *D* = 40 m. For this purpose, the measurement distance *x* and the scan radius *r* are varied in a brute force optimization (see Fig. 5). The value of the wavenumber *x* = 40 m and *r* = 12 m.

Scope of a lidar system with three beams.

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

Scope of a lidar system with three beams.

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

Scope of a lidar system with three beams.

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

Maximum coherent wavenumber ^{−1}) of a lidar system with three beams at different distances and radii on a turbine with *D* = 40 m. The optimum is marked with a dot.

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

Maximum coherent wavenumber ^{−1}) of a lidar system with three beams at different distances and radii on a turbine with *D* = 40 m. The optimum is marked with a dot.

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

Maximum coherent wavenumber ^{−1}) of a lidar system with three beams at different distances and radii on a turbine with *D* = 40 m. The optimum is marked with a dot.

Citation: Journal of Atmospheric and Oceanic Technology 30, 10; 10.1175/JTECH-D-13-00077.1

## 7. Conclusions and outlook

This work presents a model to estimate the correlation of the rotor effective wind speed between a lidar system and a wind turbine, considering different rotor diameters, spatial averaging, different scanning patterns, and wind evolution, which—in combination—are responsible for the level of correlation. The derived model is modular, and all changes can be applied in a straightforward manner that makes it usable to optimize the configuration and scan pattern of lidar systems. The model is evaluated with data from a field testing campaign, where the results were used to rate lidar measurements and to design the filter that is the crucial part of the successfully tested feedforward collective pitch controller. In future work the model will be extended to model the correlation of rotor effective shears for lidar-assisted cyclic pitch control.

## Acknowledgments

Thanks to Andreas Rettenmeier, Timo Maul, Davide Trabucchi, Eric Simley, and Ameya Sathe for the great discussions. Thanks to all persons involved in the ongoing control campaigns at NREL. Part of this research is funded by the German Federal Ministry for the Environment, Nature Conservation and Nuclear Safety (BMU) in the framework of the German joint research project “LIDAR II.”

## REFERENCES

Banakh, V. A., and Smalikho I. N. , 1997: Estimation of the turbulence energy dissipation rate from the pulsed Doppler lidar data.

,*Atmos. Oceanic Opt.***10**, 957–965.Bendat, J. S., and Piersol A. G. , 1971:

*Random Data: Analysis and Measurement Procedures*. John Wiley & Sons, 407 pp.Bossanyi, E., 2012: Un-freezing the turbulence: Improved wind field modelling for investigating lidar-assisted wind turbine control.

*Proc. European Wind Energy Association Annual Event,*Copenhagen, Denmark, EWEA. [Available online at http://proceedings.ewea.org/annual2012/allfiles2/912_EWEA2012presentation.pdf.]Dunne, F., Schlipf D. , Pao L. Y. , Wright A. D. , Jonkman B. , Kelley N. , and Simley E. , 2012: Comparison of two independent lidar-based pitch control designs.

*Proc. 50th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition,*Nashville, TN, AIAA, AIAA 2012-1151.Henriksen, L. C., 2011: Model predictive control of wind turbines. Ph.D. thesis, Technical University of Denmark, 132 pp.

IEC, 2005: Wind turbines—Part 1: Design requirements. IEC 61400-1, 92 pp.

Jonkman, B. J., and Buhl M. L. , 2007: TurbSim user's guide. NREL Tech. Rep. NREL/TP-500-41136, 35 pp.

Kapp, S., and Kühn M. , 2012: A five-parameter wind field estimation method based on spherical upwind lidar measurements.

*Proc. Conf. on the Science of Making Torque from Wind,*Oldenburg, Germany, EAWE, in press.Koerber, A., and King R. , 2011: Nonlinear model predictive control for wind turbines.

*Proc. European Wind Energy Association 2011 Annual Event,*Brussels, Belgium, EWEA, 6 pp. [Available online at http://proceedings.ewea.org/annual2011/allfiles2/1105_EWEA2011presentation.pdf.]Kragh, K. A., Hansen M. H. , and Mikkelsen T. , 2013: Precision and shortcomings of yaw error estimation using spinner-based light detection and ranging.

,*Wind Energy***16**, 353–366, doi:10.1002/we.1492.Kristensen, L., 1979: On longitudinal spectral coherence.

,*Bound.-Layer Meteor.***16**, 145–153.Laks, J., Pao L. Y. , Wright A. , Kelley N. , and Jonkman B. , 2010: Blade pitch control with preview wind measurements.

*Proc. 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition,*Orlando, FL, AIAA, AIAA 2010-251.Laks, J., Simley E. , and Pao L. , 2013: Spectral models for evaluating the effect of wind evolution on wind turbine preview control.

*Proc. American Control Conf.,*Washington, DC, IEEE CCS and AACC, TuB18.5.Lindelöw, P., 2008: Fiber based coherent lidars for remote wind sensing. Ph.D. thesis, Technical University of Denmark, 163 pp.

Pielke, R. A., and Panofsky H. A. , 1970: Turbulence characteristics along several towers.

,*Bound.-Layer Meteor.***1**, 115–130.Schlipf, D., and Kühn M. , 2008: Prospects of a collective pitch control by means of predictive disturbance compensation assisted by wind speed measurements.

*Proc. Ninth German Wind Energy Conf.,*Bremen, Germany, German Wind Energy Institute, 4 pp. [Available online at http://elib.uni-stuttgart.de/opus/volltexte/2013/8433/.]Schlipf, D., Trujillo J. J. , Basterra V. , and Kühn M. , 2009: Development of a wind turbine LIDAR simulator.

*Proc. European Wind Energy Conf.,*Marseille, France, EWEA, 6 pp. [Available online at http://proceedings.ewea.org/ewec2009/allfiles2/655_EWEC2009presentation.pdf.]Schlipf, D., Fischer T. , Carcangiu C. E. , Rossetti M. , and Bossanyi E. , 2010a: Load analysis of look-ahead collective pitch control using lidar.

*Proc. 10th German Wind Energy Conf.,*Bremen, Germany, German Wind Energy Institute, 4 pp. [Available online at http://elib.uni-stuttgart.de/opus/volltexte/2013/8432/.]Schlipf, D., and Coauthors, 2010b: Testing of frozen turbulence hypothesis for wind turbine applications with a scanning lidar system.

*Proc. 15th Int. Symp. for the Advancement of Boundary Layer Remote Sensing,*Paris, France, IPSL, O-WIN/02.Schlipf, D., Kapp S. , Anger J. , Bischoff O. , Hofsäß M. , Rettenmeier A. , and Kühn M. , 2011: Prospects of optimization of energy production by lidar assisted control of wind turbines.

*Proc. European Wind Energy Association 2011 Annual Event,*Brussels, Belgium, EWEA, 10 pp. [Available online at http://proceedings.ewea.org/annual2011/allfiles2/914_EWEA2011presentation.pdf.]Schlipf, D., Fleming P. , Haizmann F. , Scholbrock A. , Hofsäß M. , Wright A. , and Cheng P. W. , 2012a: Field testing of feedforward collective pitch control on the CART2 using a nacelle-based lidar scanner.

*Proc. Conf. on the Science of Making Torque from Wind,*Oldenburg, Germany, EAWE, 3-8.Schlipf, D., Rettenmeier A. , Hofsäß M. , Courtney M. , and Cheng P. W. , 2012b: Model based wind vector field reconstruction from lidar data.

*Proc. 11th German Wind Energy Conf.,*Bremen, Germany, German Wind Energy Institute. [Available online at http://elib.uni-stuttgart.de/opus/volltexte/2013/8431/.]Schlipf, D., Fleming P. , Kapp S. , Scholbrock A. , Haizmann F. , Belen F. , Wright A. , and Cheng P. W. , 2013a: Direct speed control using lidar and turbine data.

*Proc. American Control Conf.,*Washington, DC, IEEE CCS and AACC, MoC18.1.Schlipf, D., Schlipf D. J. , and Kühn M. , 2013b: Nonlinear model predictive control of wind turbines using lidar.

, doi:10.1002/we.1533, in press.*Wind Energy*Scholbrock, A., Fleming P. , Fingersh L. , Wright A. , Schlipf D. , Haizmann F. , and Belen F. , 2013: Field testing lidar based feed-forward controls on the NREL controls advanced research turbine.

*Proc. 51st AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition,*Grapevine, TX, AIAA, 2013-0818.Simley, E., Pao L. Y. , Kelley N. , Jonkman B. , and Frehlich R. , 2012: Lidar wind speed measurements of evolving wind fields.

*Proc. 50th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition,*Nashville, TN, AIAA, AIAA-2012-0656 .Simley, E., Pao L. Y. , Frehlich R. , Jonkman B. , and Kelley N. , 2013: Analysis of light detection and ranging wind speed measurements for wind turbine control.

, doi:10.1002/we.1584, in press.*Wind Energy*Taylor, G., 1938: The spectrum of turbulence.

,*Proc. Roy. Soc. London***132A**, 476–490.