Journal topic
Wind Energ. Sci., 4, 677–692, 2019
https://doi.org/10.5194/wes-4-677-2019
Wind Energ. Sci., 4, 677–692, 2019
https://doi.org/10.5194/wes-4-677-2019

Research article 18 Dec 2019

Research article | 18 Dec 2019

# Uncertainty identification of blade-mounted lidar-based inflow wind speed measurements for robust feedback–feedforward control synthesis

Uncertainty identification of blade-mounted lidar-based inflow wind speed measurements for robust feedback–feedforward control synthesis
Róbert Ungurán1, Vlaho Petrović1, Lucy Y. Pao2, and Martin Kühn1 Róbert Ungurán et al.
• 1ForWind, Institute of Physics, University of Oldenburg, Küpkersweg 70, 26129 Oldenburg, Germany
• 2Department of Electrical, Computer, and Energy Engineering, University of Colorado, Boulder, USA

Correspondence: Róbert Ungurán (robert.unguran@uni-oldenburg.de)

Abstract

1 Introduction

The ongoing trend of steadily growing rotor diameters of wind turbines results in dynamic loads across the rotor swept area which are becoming more uneven. Due to the so-called rotational sampling or eddy slicing effect, the blade samples the inhomogeneous wind field with frequencies determined by the rotor speed. Hence, the dynamic blade loads are concentrated at the multiples of the rotational frequency, i.e. 1P, 2P, 3P, …, nP .

The scope of this paper is particularly geared to the relevance of three aspects of recent developments in controls to mitigate such loading. First, the control surfaces on the rotor are becoming more localized and consequently in addition to individual (blade) pitch control, local active or passive blade load mitigation concepts (e.g. trailing edge flaps) have been researched for several years. Second, in addition to the proven feedback control based on rotor speed and individual blade root bending moment measurements, feedforward control using either observer techniques or lidar-assisted preview information of the inflow has been investigated for collective and individual pitch as well as trailing edge flap control. Third, there are methods that can be applied in the feedback–feedforward controller design to guarantee robust stability and performance in the presence of inherent uncertainties in the lidar measurement.

For obvious reasons, it is necessary to consider the uncertainties in the lidar measurements to achieve robust stability and performance of the feedback–feedforward controller. Furthermore, the source of such uncertainties must be identified and modelled, which can then be incorporated into the design and analysis of the controller, to ensure performance even for uncertain lidar measurements. Several authors have already addressed this problem, e.g. , , and with their numerical investigations. performed field tests to assess the influence of the “cyclops dilemma”, spatial averaging error, induction zone, and wind evolution on a hub-mounted lidar measurement. used a hub-mounted continuous-wave (CW) lidar to investigate the effect of the cyclops dilemma and concluded that a compromise in the preview distance existed. Spatial averaging increases with increasing distance from the rotor plane, leading to correlation attenuation between the rotor-effective wind speed and the lidar-estimated inflow wind speed, with increasing frequency. As measurements are taken closer to the rotor plane, the contribution of the lateral and vertical wind components to the line-of-sight lidar measurements also increases. Thus, it is not possible to accurately reconstruct the longitudinal wind component from a single hub-mounted lidar system, which results in over- or underestimation of the rotor effective wind speed. investigated how wind evolution affects controller performance; they used a single-point measurement, without spatial averaging, in front of the wind turbine blade as a feedforward IPC input. Using the feedback–feedforward IPC, they acquired the highest load reduction at the blade root bending moment at a preview time of only 0.2 s. The further the measurement was taken from the rotor plane, the more the wind evolved at high frequencies (i.e. the so-called “wind evolution”), leading to overactuation by the feedforward IPC. It should be noted that the required preview time depends on many factors, e.g. wind turbine size, 1P frequency, inflow wind speed, and induced phase shift by the feedforward controller and blade pitch actuators.

The blade-mounted lidar system is a novel technique that enables us to sample the wind component parallel to the rotor shaft axis around the swept area (Bossanyi2013) and has been demonstrated to be technologically viable . Such a feature of the system enables addressing the mitigation of higher-harmonic dynamic blade loads through feedback–feedforward individual pitch and trailing edge flap controllers , while simultaneously posing challenges with the presence of the induction zone. The closer the lidar measurement is taken to the rotor plane, the higher the deficit between the measured inflow and free flow wind speeds. Additionally, this deficit depends on where the lidar is mounted along the blade radius, which shows the importance of analysing how the blade-mounted lidar measurement is affected by wind evolution, the induction zone, and the assumptions made during the inflow wind speed reconstruction.

Therefore, in this study, our objective is to identify the nominal measurement transfer functions and model the uncertainties of the blade-mounted lidar measurement as a frequency-dependent uncertain weight for inclusion in the feedback–feedforward individual pitch and trailing edge flap control development and to analyse the impact of the induction zone effect on the preview time.

The rest of the paper is organized as follows. Sect. 2 provides a description of the framework and methods we use for identifying the uncertainties and preview time of the blade-mounted lidar measurement, after an introduction of the blade-mounted lidar-based simulation set-up in Sect. 2.1. In Sect. 2.2 we describe the method we use to estimate the inflow wind speed. The method we employed for determining the blade effective wind speed to assess the efficiency of the blade-mounted lidar-based inflow wind speed measurement is discussed in Sect. 2.3. Section 2.4 describes the general control implementation and presents the multiblade coordinate transformation and its importance in the controller design, while Sect. 2.5 details how the lidar-based measurement uncertainty is considered in control development and analysis. Section 2.6 proposes a method to identify the uncertainties of the blade-mounted lidar measurement as a frequency-dependent uncertainty weight, Sect. 2.7 presents the method we apply for estimating the preview time, and Sect. 2.8 introduces a cost function which we use to evaluate the initially selected lidar and telescope parameters. The results of a reference case are presented in Sect. 3, where in Sect. 3.1 we analyse the effect of the multiblade coordinate transformation on the measurement. The simulation set-up is established in Sect. 3.2, and we systematically analyse the uncertainties of various telescope and control parameters in Sect. 3.3. The results are discussed in Sect. 4 prior to the conclusions in Sect. 5.

2 Methodology

## 2.1 Blade-mounted lidar

A telescope is mounted on each blade and is connected to a hub-based continuous-wave lidar with fibre optical cables. The lidar samples the inflow wind speed in front of the rotor plane at a rate of 5 Hz, and we intend to use the lidar measurements for control purposes. The lidar measurements are integrated into the system model according to Fig. 1, and we use a combination of large-eddy simulations and an aeroelastic simulation code to simulate and evaluate the lidar-based inflow measurements. Thus, lidar measurements are simulated in a realistic environment, where the effect of the induction zone and wind evolution, as well as the dynamic behaviour of the wind turbine, are taken into account. Moreover, the lidar simulator considers volumetric measurement, dynamics of the blade and tower, i.e. displacement, rotation, and linear velocity in 3-D space, and blade-rotation-induced velocity. Nevertheless, the rotational effect of the blade is not accounted for during the accumulation of a single measurement.

Figure 1Block diagram of the blade-mounted lidar-based simulation set-up. LOS corresponds to line of sight.

Figure 2Configuration of the lidar measurement system, with a telescope mounted on each blade and connected to a continuous-wave lidar in the hub via fibre optics. The line-of-sight wind speed is computed on the basis of a weighting function (W(Fξ)), which is dependent on the focus distance (F) and the range along the beam (ξ).

Figure 2 illustrates the coordinate systems and the telescope orientation. Here, the line-of-sight (LOS) wind speed measurement from blade i (ulos,i) is defined as

$\begin{array}{}\text{(1)}& {u}_{\mathrm{los},i}=\frac{\underset{{\mathit{\xi }}_{\mathrm{min}}}{\overset{{\mathit{\xi }}_{\mathrm{max}}}{\int }}W\left(F,\mathit{\xi }\right){V}_{i}\left(\mathit{\xi }\right)\mathrm{d}\mathit{\xi }}{\underset{{\mathit{\xi }}_{\mathrm{min}}}{\overset{{\mathit{\xi }}_{\mathrm{max}}}{\int }}W\left(F,\mathit{\xi }\right)\mathrm{d}\mathit{\xi }},\end{array}$

where Vi(ξ) is defined in Eq. (3) and W(F, ξ) is the lidar's weighting function, defined according to as

$\begin{array}{}\text{(2)}& W\left(F,\mathit{\xi }\right)=\frac{\mathrm{1}}{{\mathit{\xi }}^{\mathrm{2}}+{\left(\mathrm{1}-\frac{\mathit{\xi }}{F}\right)}^{\mathrm{2}}{R}_{\mathrm{R}}^{\mathrm{2}}},\end{array}$

where RR is the Rayleigh range, set at 1573 m herein, as proposed by ; F is the focus distance and ξ is the range along the beam. The limits ξmin and ξmax, introduced in Eq. (1), refer to the minimum and maximum range, respectively, along the beam. For a practical implementation of the lidar simulator, these values are chosen such that $\frac{W\left(F,\mathit{\xi }\right)}{W\left(F,F\right)}$ equals 0.02 at these limits. During discretization of Eq. (1), the spatial resolution is set empirically at Δξ=0.1 m. A single-point measurement is given by

$\begin{array}{}\text{(3)}& {\mathbf{V}}_{i}\left(\mathit{\xi }\right)={\left(\left[\begin{array}{c}{u}_{\mathrm{h},i}\left(\mathit{\xi }\right)\\ {v}_{\mathrm{h},i}\left(\mathit{\xi }\right)\\ {w}_{\mathrm{h},i}\left(\mathit{\xi }\right)\end{array}\right]-\left[\begin{array}{c}{\stackrel{\mathrm{˙}}{x}}_{\mathrm{t},\mathrm{h},i}\\ {\stackrel{\mathrm{˙}}{y}}_{\mathrm{t},\mathrm{h},i}\\ {\stackrel{\mathrm{˙}}{z}}_{\mathrm{t},\mathrm{h},i}\end{array}\right]\right)}^{T}\left[\begin{array}{c}{\mathrm{\ell }}_{x,\mathrm{h},i}\\ {\mathrm{\ell }}_{y,\mathrm{h},i}\\ {\mathrm{\ell }}_{z,\mathrm{h},i}\end{array}\right],\end{array}$

where $\left[{u}_{\mathrm{h},i}{v}_{\mathrm{h},i}{w}_{\mathrm{h},i}{\right]}^{T}$ is the wind speed vector along the laser beam expressed in the rotating hub coordinate system; $\left[{\stackrel{\mathrm{˙}}{x}}_{\mathrm{t},\mathrm{h},i}{\stackrel{\mathrm{˙}}{y}}_{\mathrm{t},\mathrm{h},i}{\stackrel{\mathrm{˙}}{z}}_{\mathrm{t},\mathrm{h},i}{\right]}^{T}$ is the linear velocity vector of the blade segment where the telescope is mounted, expressed in the rotating hub frame of reference; and $\left[{\mathrm{\ell }}_{x,\mathrm{h},i}{\mathrm{\ell }}_{y,\mathrm{h},i}{\mathrm{\ell }}_{z,\mathrm{h},i}{\right]}^{T}$ is the unit vector of the laser beam in the rotating hub coordinate system. The aeroelastic simulation tool is capable of providing full kinematics information, i.e. positions, orientations, and linear and angular velocities, of any blade segment in the hub coordinate system.

## 2.2 Wind speed estimation

During the inflow wind speed estimation, the velocity, displacement, and rotation of the blade segment are assumed to be known; therefore, the wind speed component parallel to the rotor shaft axis can be reconstructed as indicated in Eq. (4). Without loss of generality, in the wind speed estimation, the weighting function of W(F,ξ) from Eq. (1) is neglected, and two assumptions are made: (1) the vh,i and wh,i components are zero, and (2) the mean wind velocity is parallel to the rotor axis; i.e. no tilt and no yaw misalignments are considered. Consequently, an estimate of the wind speed parallel to the rotor shaft axis (${u}_{\mathrm{h},\mathrm{est},i}$) is

$\begin{array}{}\text{(4)}& {u}_{\mathrm{h},\mathrm{est},i}\approx \frac{{u}_{\mathrm{los},i}+{\stackrel{\mathrm{˙}}{y}}_{\mathrm{t},\mathrm{h},i}{\mathrm{\ell }}_{y,\mathrm{h},i}+{\stackrel{\mathrm{˙}}{z}}_{\mathrm{t},\mathrm{h},i}{\mathrm{\ell }}_{z,\mathrm{h},i}}{{\mathrm{\ell }}_{x,\mathrm{h},i}}+{\stackrel{\mathrm{˙}}{x}}_{\mathrm{t},\mathrm{h},i}.\end{array}$

Nevertheless, such assumptions introduce errors into the lidar measurement that are presumed to exist in the identified uncertainty weight and, thus, are consequently considered during the controller development.

Figure 3Normalized longitudinal inflow wind speed ($\frac{u\left(F,R\right)}{{u}_{\mathrm{0}}}$) as a function of focus distance (F) and blade span position (R), with an undisturbed inflow wind speed u0=13 m s−1.

Figure 3 illustrates the induction zone effect for laminar inflow. Note that the lidar measurement is affected by the rotor induction. The reduction depends on the position of the telescope along the blade radius (R) and the focus distance of the laser beam (F), where the wind speed measurement takes place. To account for this effect in the lidar-based inflow wind speed measurement, we construct a second-order polynomial function (f), whose inputs are chosen as rotor speed (ωr), blade pitch angle (βi), and blade root flapwise and edgewise moments (Mfw,i, Mew,i). Rotor speed and blade pitch angles are easily measured, and we assume that blade root flapwise and edgewise moment sensors are also available for implementing this method. Therefore, the estimated wind speed parallel to the rotor shaft axis (${u}_{\mathrm{h},\mathrm{est},i}$) is corrected as

$\begin{array}{}\text{(5)}& {u}_{\mathrm{cor},i}={u}_{\mathrm{h},\mathrm{est},i}+\mathrm{\Delta }{u}_{\mathrm{est},i},\end{array}$

where

$\begin{array}{}\text{(6)}& {u}_{\mathrm{0}}-u\left(F,R\right)\approx \mathrm{\Delta }{u}_{\mathrm{est},i}=f\left(F,R,{\mathit{\omega }}_{\mathrm{r}},{\mathit{\beta }}_{i},{M}_{\mathrm{fw},i},{M}_{\mathrm{ew},i}\right).\end{array}$

The second-order polynomial function (f) is fitted to the data extracted from 10 min large-eddy simulations with laminar inflow for mean wind speeds between 4 and 25 m s−1. u(FR) is the wind speed at an upstream distance F from the blade and at a blade radial position of R, and u0 is taken from the same blade radial position of R but at an upstream distance of 3 times the rotor diameter (3-D).

## 2.3 Blade effective wind speed

To assess the performance efficiency of the blade-mounted lidar-based inflow wind speed measurement, we introduce a new signal called the blade-effective wind speed (ubeff,i), which is determined as the contribution of the inflow wind speed on each blade segment ui(r) to the flapwise blade root bending moment; the inflow wind speed refers to the longitudinal wind speed in the rotor axis direction. The contribution depends on the radial distance (r) and the local thrust coefficient (CT) of the blade segment as expressed by

$\begin{array}{}\text{(7)}& {u}_{\mathrm{beff},i}=\sqrt{\frac{\underset{{R}_{\mathrm{hub}}}{\overset{{R}_{\mathrm{tip}}}{\int }}{C}_{\mathrm{T}}\left(r,{u}_{i}\left(r\right)\right){r}^{\mathrm{2}}{u}_{i}^{\mathrm{2}}\left(r\right)\mathrm{d}r}{\underset{{R}_{\mathrm{hub}}}{\overset{{R}_{\mathrm{tip}}}{\int }}{C}_{\mathrm{T}}\left(r,{u}_{i}\left(r\right)\right){r}^{\mathrm{2}}\mathrm{d}r}}.\end{array}$

The local thrust coefficients are resolved from steady-state simulations for each blade segment from cut-in to cut-out wind speeds.

## 2.4 Multiblade coordinate transformation (MBC)

In the subsequent step, we introduce the multiblade coordinate transformation (MBC) that simplifies the controller design by transforming a time-varying system into a time-invariant system and decouples the individual pitch from the collective pitch control.

Figure 4Implementation of the feedforward collective and individual pitch control, where the inputs (ucor,1, ucor,2, and ucor,3) are the estimated wind speeds parallel to the rotor shaft axis and the outputs are the blade pitch angles (β1, β2, and β3). The feedforward controller (Kff,f) is implemented in the non-rotating (fixed) frame of reference and is, therefore, denoted with an extra index “f”. Further, the multiblade coordinate transformation (Tmbc) is applied to the inputs, and the pseudo-inverse transformation (${\mathbf{T}}_{\mathrm{mbc}}^{+}$) is applied to the outputs.

Figure 4 demonstrates the manner in which the feedforward controller is implemented. First, the measured inflow wind speed is transformed to the non-rotating frame of reference by applying MBC transformation (Tmbc(θ+ϕ)) in accordance with Eq. (8), where θ denotes the azimuth angle.

$\begin{array}{}\text{(8)}& \left[\begin{array}{l}{u}_{\mathrm{cor},\mathrm{col}}\\ {u}_{\mathrm{cor},\mathrm{yaw}}\\ {u}_{\mathrm{cor},\mathrm{tilt}}\end{array}\right]={\mathbf{T}}_{\mathrm{mbc}}\left(\mathit{\theta }+\mathit{\varphi }\right)\left[\begin{array}{c}{u}_{\mathrm{cor},\mathrm{2}}\\ {u}_{\mathrm{cor},\mathrm{3}}\\ {u}_{\mathrm{cor},\mathrm{1}}\end{array}\right],\end{array}$

where

$\begin{array}{}\text{(9)}& {\mathbf{T}}_{\mathrm{mbc}}\left(\mathit{\theta }\right)=\left[\begin{array}{ccc}\frac{\mathrm{1}}{\mathrm{3}}& \frac{\mathrm{1}}{\mathrm{3}}& \frac{\mathrm{1}}{\mathrm{3}}\\ \frac{\mathrm{2}}{\mathrm{3}}\mathrm{cos}\left({n}_{\mathrm{h}}\mathit{\theta }\right)& \frac{\mathrm{2}}{\mathrm{3}}\mathrm{cos}\left({n}_{\mathrm{h}}\left[\mathit{\theta }+\frac{\mathrm{2}\mathit{\pi }}{\mathrm{3}}\right]\right)& \frac{\mathrm{2}}{\mathrm{3}}\mathrm{cos}\left({n}_{\mathrm{h}}\left[\mathit{\theta }+\frac{\mathrm{4}\mathit{\pi }}{\mathrm{3}}\right]\right)\\ \frac{\mathrm{2}}{\mathrm{3}}\mathrm{sin}\left({n}_{\mathrm{h}}\mathit{\theta }\right)& \frac{\mathrm{2}}{\mathrm{3}}\mathrm{sin}\left({n}_{\mathrm{h}}\left[\mathit{\theta }+\frac{\mathrm{2}\mathit{\pi }}{\mathrm{3}}\right]\right)& \frac{\mathrm{2}}{\mathrm{3}}\mathrm{sin}\left({n}_{\mathrm{h}}\left[\mathit{\theta }+\frac{\mathrm{4}\mathit{\pi }}{\mathrm{3}}\right]\right)\end{array}\right].\end{array}$

A phase shift (ϕ) is introduced into the transformation to consider that the measured inflow wind speed hits the wind turbine blade after this azimuth angle change. This value varies with respect to several parameters, including the selected focus distance, inflow wind speed, and rotor speed. The estimated wind speed parallel to the rotor shaft axis from blade 1 is used to determine the blade pitch control at blade 1; hence, the order of the estimated wind speeds parallel to the rotor shaft axis has changed to ucor,2, ucor,3, and ucor,1. Further, the control signals or the blade pitch angles (βcol, βyaw, βtilt) are determined by the feedforward controller (Kff,f). If the preview time provided by the lidar is greater than the time delay induced by the feedforward controller, an additional time delay (${e}^{-s{T}_{\mathrm{id}}}$) is introduced into the system. Finally, the delayed control signals (βcol,d, βyaw,d, and βtilt,d) are transformed to the rotating frame of reference using the pseudo-inverse of the MBC transformation (${\mathbf{T}}_{\mathrm{mbc}}^{+}\left(\mathit{\theta }\right)$). The main structure of the feedforward individual pitch controller in Fig. 4 can be used in the feedforward trailing edge flap controller as well.

The MBC transformation plays a considerably important role because it can transform a frequency component of interest, such as 1P, 2P, or 3P , to a low-frequency component, named 0P. It is dependent on the selected value of nh in Eq. (9). For example, 1P will be transformed to 0P when nh is specified as 1, and 2P will be transformed to 0P when nh is specified as 2.

In this study, we focus on identifying the uncertainty weight that can be used during the feedback–feedforward individual and collective pitch control development with an objective to mitigate the 1P loads at the flapwise blade root bending moments and to enhance the rotor speed tracking. This indicates that by considering nh as 1, the measured inflow wind speeds are transformed to the non-rotating frame of reference in Eq. (9), where the uncertainty weight identification is conducted. Further, the same methodology can be applied to identify the uncertainty weight for higher-harmonics control by selecting a larger integer value of nh.

We have already mentioned that the measured inflow wind speeds are transformed to the non-rotating frame of reference by applying the MBC transformation. In order to assess the performance efficiency of the blade-mounted lidar-based inflow wind speed measurement, the blade effective wind speeds are also transferred to the non-rotating frame using the MBC transformation as follows

$\begin{array}{}\text{(10)}& \left[\begin{array}{l}{u}_{\mathrm{bef},\mathrm{col}}\\ {u}_{\mathrm{bef},\mathrm{yaw}}\\ {u}_{\mathrm{bef},\mathrm{tilt}}\end{array}\right]={\mathbf{T}}_{\mathrm{mbc}}\left(\mathit{\theta }\right)\left[\begin{array}{c}{u}_{\mathrm{bef},\mathrm{1}}\\ {u}_{\mathrm{bef},\mathrm{2}}\\ {u}_{\mathrm{bef},\mathrm{3}}\end{array}\right],\end{array}$

where Tmbc(θ) is defined in Eq. (9).

## 2.5 System modelling with uncertain lidar measurements

We use the blade-mounted telescopes to measure the disturbance, or the inflow wind speed in this case. Afterward, the three measurements are transformed into the non-rotating frame of reference where they are used as inputs to the feedforward individual and collective pitch controllers. Figure 5 illustrates the disturbance rejection controller set-up with uncertainty. Each block in the figure represents a three-input and three-output system with a three-by-three matrix transfer function.

Figure 5Block diagram of the disturbance rejection control design with performance weight and uncertain input measurement. Kfb,f and Kff,f are the feedback and feedforward controllers, Gwt,f is the wind turbine model from the control input to output, Gd,f is the wind turbine model from the disturbance to the output, Gn,f is the nominal disturbance measurement model, Δ is the uncertainty, W is the measurement uncertainty weight, and Wp is the performance weight. The “f” in the index refers to the non-rotating (fixed) frame of reference.

The control development is aimed at achieving disturbance rejection up to a certain frequency with measurement uncertainties. In other words, we want to find a controller that satisfies Eq. (11) for a chosen performance weight Wp.

$\begin{array}{}\text{(11)}& \parallel {\mathbf{W}}_{\mathrm{p}}{\mathbf{S}}_{\mathrm{fb}}{\mathbf{S}}_{\mathrm{ff},\mathrm{p}}{\parallel }_{\mathrm{\infty }}<\mathrm{1},\end{array}$

where the frequency-dependent feedback (Sfb) and feedforward sensitivity (Sff,p) functions with additive uncertainty are given by

$\begin{array}{}\text{(12)}& & {\mathbf{S}}_{\mathrm{fb}}={\left(\mathbf{I}+{\mathbf{G}}_{\mathrm{wt},\mathrm{f}}{\mathbf{K}}_{\mathrm{fb},\mathrm{f}}\right)}^{-\mathrm{1}},\text{(13)}& & {\mathbf{S}}_{\mathrm{ff},\mathrm{p}}=\mathbf{I}+{\mathbf{G}}_{\mathrm{wt},\mathrm{f}}{\mathbf{K}}_{\mathrm{ff},\mathrm{f}}\left({\mathbf{G}}_{\mathrm{n},\mathrm{f}}+{\mathbf{\Delta }}_{\mathrm{\ell }}{\mathbf{W}}_{\mathrm{\ell }}\right){\mathbf{G}}_{\mathrm{d},\mathrm{f}}^{-\mathrm{1}},\end{array}$

and

$\begin{array}{}\text{(14)}& {\mathbf{\Delta }}_{\mathrm{\ell }}=\left[\begin{array}{ccc}{\mathit{\delta }}_{\mathrm{\ell },\mathrm{1}}& \mathrm{0}& \mathrm{0}\\ \mathrm{0}& {\mathit{\delta }}_{\mathrm{\ell },\mathrm{2}}& \mathrm{0}\\ \mathrm{0}& \mathrm{0}& {\mathit{\delta }}_{\mathrm{\ell },\mathrm{3}}\end{array}\right]\phantom{\rule{0.25em}{0ex}}\in {\mathbb{C}}^{\mathrm{3}×\mathrm{3}},\phantom{\rule{0.25em}{0ex}}\parallel {\mathrm{\Delta }}_{\mathrm{\ell }}{\parallel }_{\mathrm{\infty }}\le \mathrm{1}.\end{array}$

This equation highlights the importance of knowing the frequency-dependent uncertainty weight W(jω) in advance, so as to ensure that the closed-loop system is stable and that the objective in Eq. (11) is satisfied for all perturbations ($||{\mathbf{\Delta }}_{\mathrm{\ell }}|{|}_{\mathrm{\infty }}\le \mathrm{1}$). For control development, the frequency dependent uncertainty weight of W(jω) and the nominal disturbance measurement model of Gn,f(jω) are missing; we discuss how they can be identified in the next subsections and later illustrate the process for the reference cases in Sect. 3.3.

Remark: only one objective is introduced in Eq. (11); nevertheless, other objectives can be added, such as penalizing the control signal magnitude at high frequencies .

## 2.6 Uncertainty modelling for control development

We employ black box system identification to establish the transfer functions (G) from the blade effective wind speeds (ubeff) to the corrected lidar-based inflow wind speeds (ucor) in the non-rotating (fixed) frame of reference

$\begin{array}{}\text{(15)}& {u}_{\mathrm{cor},\mathrm{f}}={\mathbf{G}}_{\mathrm{\ell }}{u}_{\mathrm{beff},\mathrm{f}}\end{array}$

with

$\begin{array}{}\text{(16)}& {\mathbf{G}}_{\mathrm{\ell }}=\left[\begin{array}{ccc}{G}_{\mathrm{\ell },\mathrm{col}}& \mathrm{0}& \mathrm{0}\\ \mathrm{0}& {G}_{\mathrm{\ell },\mathrm{yaw}}& \mathrm{0}\\ \mathrm{0}& \mathrm{0}& {G}_{\mathrm{\ell },\mathrm{tilt}}\end{array}\right]\phantom{\rule{0.25em}{0ex}}\in {\mathbb{C}}^{\mathrm{3}×\mathrm{3}}.\end{array}$

The system identification is performed via the ssest function from with a 15th-order state-space model, which can capture all the relevant information. The order of the state-space model is found empirically through analysis of the Hankel singular values.

We separately identify the nominal disturbance measurement model (Gn,k(jω)) and the uncertainty weight (wℓ,k(jω)), where k{col, yaw, tilt}, as a fifth-order minimum phase filter for each of the inputs in such a way as to satisfy the following inequalities

$\begin{array}{}\text{(17)}& \left|{G}_{\mathrm{n},k}\left(j\mathit{\omega }\right)\right|<\left|{G}_{\mathrm{\ell },k}\left(j\mathit{\omega }\right)\right|,\phantom{\rule{0.25em}{0ex}}\forall \mathit{\omega }\end{array}$

and

$\begin{array}{}\text{(18)}& \left|{G}_{\mathrm{n},k}\left(j\mathit{\omega }\right)+{w}_{\mathrm{\ell },k}\left(j\mathit{\omega }\right)\right|>\left|{G}_{\mathrm{\ell },k}\left(j\mathit{\omega }\right)\right|,\phantom{\rule{0.25em}{0ex}}\forall \mathit{\omega },\end{array}$

leading to the diagonal nominal disturbance measurement model matrix of

$\begin{array}{}\text{(19)}& {\mathbf{G}}_{\mathrm{n}}=\left[\begin{array}{ccc}{G}_{\mathrm{n},\mathrm{col}}& \mathrm{0}& \mathrm{0}\\ \mathrm{0}& {G}_{\mathrm{n},\mathrm{yaw}}& \mathrm{0}\\ \mathrm{0}& \mathrm{0}& {G}_{\mathrm{n},\mathrm{tilt}}\end{array}\right]\end{array}$

and uncertainty weight matrix of

$\begin{array}{}\text{(20)}& {\mathbf{W}}_{\mathrm{\ell }}=\left[\begin{array}{ccc}{w}_{\mathrm{\ell },\mathrm{col}}& \mathrm{0}& \mathrm{0}\\ \mathrm{0}& {w}_{\mathrm{\ell },\mathrm{yaw}}& \mathrm{0}\\ \mathrm{0}& \mathrm{0}& {w}_{\mathrm{\ell },\mathrm{tilt}}\end{array}\right].\end{array}$

The order of the transfer functions are determined empirically during the analysis of the data. Lower orders could be selected as well; however, these would lead to higher uncertainties at high frequency.

The ideal case would be to measure, with a telescope, the exact inflow wind speed hitting the rotor blades to result in a nominal disturbance measurement transfer function with a gain of 1 over the entire frequency range. However, not only the inflow condition, but also the telescope parameters influence the nominal disturbance measurement model and the measurement uncertainty weight. In Sect. 3.3 we identify these transfer functions (Gn,k and wℓ,k), which can then be used for control development and analysis. Furthermore, we analyse how much the low-frequency gains of G deviate from 1 for several cases.

We neglect the cross-coupling between the yaw and tilt components in the system identification, but these are considered in the wind turbine and disturbance transfer functions in line with , so that the cross-coupling between the yaw and tilt components is included in the controller development.

## 2.7 Preview time estimation

Preview time plays an important role in the development of feedforward control. It must be larger than or equal to the time delay introduced by the feedforward controller and actuator dynamics. It is preferable for it to be equal, but a larger value is acceptable, as an additional time delay can be easily introduced into the feedforward controller, as shown in Fig. 4. To determine the optimal preview time for a given focus distance, we evaluate the cross-correlation between the blade effective (ubeff,k) and the corrected inflow (ucor,k) wind speeds, with k{col, yaw, tilt}, and we choose the index of the peak value as the available preview time.

## 2.8 Telescope parameter estimation

We introduce a cost function which is based only on the coherence (${\mathit{\gamma }}_{k}^{\mathrm{2}}$) between the blade effective (ubeff,k) and the corrected inflow (ucor,k) wind speeds, with k{col, yaw, tilt}:

$\begin{array}{}\text{(21)}& {J}_{\mathrm{lp}}\left(f\right)=\sum _{k}{J}_{\mathrm{lp},k}\left(f\right)=\sum _{k}{\mathit{\gamma }}_{k}^{\mathrm{2}}\left(f\right).\end{array}$

By evaluating Jlp for the discrete set of sampled lidar and telescope parameters, the maximum of the objective function results in the optimal telescope parameters within the discrete set of sampled lidar and telescope parameters. In this way, we are able to judge the initially chosen telescope parameters.

3 Results

## 3.1 Multiblade coordinate transformation effect on the blade-mounted lidar measurement

To perform an analysis of the MBC transformation, we create three generic wind speed measurement signals with

$\begin{array}{}\text{(22)}& {u}_{\mathrm{cor},i}={u}_{\mathrm{0}}+\sum _{j=\mathrm{1}}^{\mathrm{6}}\frac{\mathrm{1}}{{j}^{\mathrm{3}}}\mathrm{sin}\left(j\left[\mathrm{2}\mathit{\pi }{f}_{\mathrm{0}}t+\left(i-\mathrm{1}\right)\frac{\mathrm{2}\mathit{\pi }}{\mathrm{3}}+\frac{\mathit{\pi }}{\mathrm{6}}\right]\right),\end{array}$

where u0, i, f0, and t are the offset or undisturbed inflow wind speed, blade index, 1P frequency, and time, respectively. Here, we considered harmonics of up to 6P (j=1 … 6). Figure 6 shows an example time series of the generated signals.

Figure 6Time series of three generic wind speed measurements at the same amplitude, used for analysing the impact of the multiblade coordinate transformation. The first, second, and third signals have a phase shift of 30, 150, and 270, respectively. The signals are constructed to include harmonics up to 6P.

Figure 7Power spectral densities of the generic signals in the rotating (ucor,1) and non-rotating (ucol, uyaw, utilt) frames of reference during the application of the multiblade coordinate transformation.

Figure 7 presents the power spectral densities of the wind speed measurement obtained from the first blade (ucor,1) and the collective (ucol), yaw (uyaw), and tilt (utilt) components after the MBC transformation, which is applied to the generic wind speed measurement signals (ucor,1, ucor,2, ucor,3). The figure highlights the MBC transformation keeping only 0P, 3P, and multiples of 3P. As describe, the frequency (f) in the non-rotating frame of reference arises from f±f0 from the rotating frame of reference; e.g. the 3P in the non-rotating frame of reference arises from the 2P and 4P contributions in the rotating frame of reference.

Several cases may illustrate the transfer of the measurement errors from the rotating to the non-rotating reference frame. First, we should consider the effect of over- or underestimation of the measured wind speed with one of the blade-mounted lidar systems, due to, e.g., different radial positions of the telescope along the blade radii or one of the telescopes having a different orientation, which reduces the DC offset (u0 in Eq. 22) for one of the three generic signals. Next, the signals are transformed into the non-rotating frame of reference, which can be compared to the case where all the DC offsets are maintained for each of the three signals at the same level.

Figure 8Power spectral densities of collective, yaw, and tilt components of the generic signals with partial DC offset. The expression u…,e indicates the case where the DC offset (u0 in Eq. 22) of one of the signals differs from the other two in the rotating frame of reference.

As Fig. 8 highlights, an undesired peak appears at 1P in the yaw and tilt components in the non-rotating frame of reference, due to the presence of asymmetries in the signals in the rotating frame of reference .

Figure 9Power spectral densities of collective, yaw, and tilt components of the generic signals with partial DC offset and phase shift. The expression u…,e indicates the case where a different DC offset is set and a 1 phase shift is added to the 1P harmonics of one of the blade signals in the rotating frame of reference.

Second, in addition to the reduction in the DC offset for one of the signals, a 1 phase shift is added to the 1P harmonic to the same signal in the rotating frame of reference, which represents the case, for example, where one of the blade-mounted lidar focus distances differs from the other two. Figure 9 reveals that after applying the MBC transformation to the three generic signals, undesired higher-harmonic peaks arise in the non-rotating frame of reference. Interestingly, the phase shift that is introduced to one of the signals in the rotating frame of reference results in different higher harmonics in the components in the non-rotating frame of reference, e.g. a peak observed at 1P of the collective component and both 1P at 2P of the tilt and yaw components.

## 3.2 Simulation set-up

The reference case we use in this investigation is based on the NREL 5 MW generic wind turbine . We use an actuator line model through the coupling between the FASTv7 aeroelastic simulation code and PALM (Parallelized Large-Eddy Simulation Model) as explained by . The operating conditions correspond to a hub-height mean wind speed of 13.06 m s−1, which is above the rated value of 11.4 m s−1. Furthermore, the 10 min simulation results in a turbulence intensity of 8.5 % and a wind shear corresponding to a power law description with an exponent of approximately 0.12. The baseline controller of the wind turbine ensures that the generator speed is kept at 1173.7 rpm , thereby resulting in a mean rotor speed (ωr) of 11.74 rpm and further leading to a 1P frequency of f0=0.195 Hz.

For an analysis of the induction zone effect, we set the range of the focus distance and telescope position along the blade radius at F∈[10, 40 m] and R∈[20, 60 m], based on a previous investigation . The range of the other input variables are determined by the results of simulations with laminar inflow and power law wind shear with coefficients of 0.1, 0.2, and 0.3. An approximation of the induction zone effect introduces some uncertainties into the measurement, but they are included in the identified uncertainty weight.

## 3.3 Nominal plants and uncertainty weights identification

stress that an elevated peak around the crossover frequency (just below the 1P frequency) of the feedback–feedforward controller sensitivity function leads to increased loads across the wind turbine components. Here, the crossover frequency of the controller is defined where the sensitivity function first crosses −3 dB from below. Uncertainties impose limitations on the achievable performance ; e.g. the peak of the sensitivity function may increase due to uncertainties in the system. Therefore, it is important to analyse how the lidar measurement uncertainty is affected by, e.g., mounting misalignment of the telescope on the blade or cases where the focus distance or position of the telescope along the blade span differs from the optimal parameters. Identifying the lidar measurement uncertainty as a frequency-dependent minimum-phase filter enables the inclusion of such parameters in the control development, allowing an analysis of its impact on the stability and performance of the closed-loop system. As we explain in detail in Sect. 3.3.1, a straightforward solution to determine the telescope and lidar parameters, such as focus distance, telescope position along the blade radius, and telescope orientation on the blade, is to assume that the blades are rigid, that the rotor speed and pitch angle are constant, and that Taylor's frozen turbulence hypothesis (Taylor1938) holds . We perform large-eddy simulation (LES) in the subsequent sections to examine the usefulness and limitations of these assumptions and further analyse the uncertainties in the blade-mounted lidar measurement as well as the measurement sensitivity with respect to lidar and telescope parameter changes. The investigated cases are described in Sect. 3.3.13.3.5 and summarized in Table 1. Section 3.3.6 describes how the measurement uncertainties are affected when one or two telescopes are aligned differently than the others. First, we assume that the orientation angle misalignment is unknown. Second, we assume that this orientation angle misalignment can be identified, so that the lidar-based inflow wind speed measurement can be corrected.

Table 1The cases investigated in this study, along with the lidar and telescope parameters for each case. If one or more parameters in the third column are not specified, then the parameters defined in the first case are used. F is the focus length, R is the radial position of the telescope along the blade, and Φℓ,i and Γℓ,i are the orientation angles of the telescope (see Fig. 2).

For each case, first the transfer functions (Gℓ,k) from the blade effective wind speeds (ubeff,k) to the corrected lidar-based inflow wind speeds (ucor,k) are identified. Next, the nominal disturbance measurement models (Gn,k) and the uncertainty weights (wℓ,k) for each of the inputs are estimated to satisfy Eqs. (17) and (18). Figure 10 provides a summary of the identified DC gain upper (${G}_{\mathrm{n},k}+{w}_{\mathrm{\ell },k}$) and lower (Gn,k) bounds of the transfer functions (Gℓ,k) from the blade effective wind speeds (ubeff,k) to the corrected lidar-based inflow wind speeds (ucor,k).

Figure 10Identified DC gain upper (${G}_{\mathrm{n},k}+{w}_{\mathrm{\ell },k}$) and lower (Gn,k) bounds of the transfer functions (Gℓ,k) from the blade effective wind speeds (ubeff,k) to the corrected lidar-based inflow wind speeds (ucor,k), with k{col, yaw, tilt}; C1C5 represent the investigated cases (outlined in Table 1).

We would like to act only below the 1P (0.195 Hz) frequency; therefore, below this frequency, it is desired that the gain of Gn,k is 1 and that the measurement uncertainty is small but still covers the worst case. A higher percentage of measurement uncertainty can be tolerated at frequencies above 1P by designing the feedforward controller accordingly, e.g. a model inversion-based feedforward controller with a low-pass filter with a crossover frequency below 1P. With Fig. 10, we show how wide variation in the DC gain of Gℓ,k agreed with the identified nominal disturbance measurement models and the additive uncertainty weights.

### 3.3.1 Telescope parameters for no-induction case (C1)

The basic concept of the feedforward controller is the use of measured inflow wind speed from blade i to control the blade and trailing edge flap angles at blade i−1. Assuming rigid blades, constant rotor speed and pitch angle and that Taylor's frozen turbulence hypothesis (Taylor1938) holds, it is easy to compute the minimum preview time of 1.7 s ($=\frac{\mathrm{2}\mathit{\pi }}{\mathrm{3}}\frac{\mathrm{30}}{\mathit{\pi }{\mathit{\omega }}_{\mathrm{r}}}$, ωr=11.74 rpm), which is the time needed for blade i−1 to reach the position of blade i, i.e. a 120 azimuth angle change. The simulation set-up presented in Sect. 3.2 results in a hub-height mean wind speed of 13.06 m s−1. The assumption that the wind evolves according to Taylor's frozen turbulence hypothesis and with the induction zone effect being negligible, a focus distance of 22.2 m (=1.7 s  13.06 m s−1) is determined. In accordance with and , the inflow at 70 % (≈44 m) of the blade radius can be assumed as most representative of the blade effective wind speed; hence, the telescope is located at this radial position. The telescope orientation angles Φℓ,i and Γℓ,i are found through aeroelastic simulation where laminar inflow is considered. The telescope orientation angles are the counter rotation of the blade segment angular orientation so that the lidar beam becomes parallel to the rotor shaft axis (see Fig. 2).

Figure 11A selected time series of the blade effective wind speed from blade 1 (ubeff,1) and the estimated (${u}_{h,\mathrm{est},\mathrm{2}}$) and corrected (ucor,2) inflow wind speeds from blade 2 in the rotating frame of reference shown in (a). The power spectral densities (PSDs) of the three signals are displayed in the lower plot.

Figure 11a shows a selected time series of the blade effective wind speed from blade 1 (ubeff,1), as well as the estimated (${u}_{\mathrm{h},\mathrm{est},\mathrm{2}}$) and corrected (ucor,2) inflow wind speeds from blade 2. The three signals are in the rotating frame of reference. The lower plot displays the power spectral densities (PSDs) of the three signals. The dominant frequencies are clearly visible, as a result of the rotational sampling of the inflow wind speed by the blade-mounted telescope. The PSD analysis highlights these dominant frequencies as 1P, 2P, and 3P. Moreover, the plot reveals a good match at 1P between ubeff,1 and ucor,2, although ucor,2 is slightly underestimated at higher harmonics.

Figure 12Power spectral densities of the blade effective wind speeds (ubeff,k) and the corrected inflow wind speeds (ucor,k) in the non-rotating frame of reference, with k{col, yaw, tilt}.

We transform the different blade effective and corrected inflow wind speeds from the rotating to the non-rotating frame of reference via the multiblade coordinate transformation (Tmbc(θ)) as discussed in Sect. 2.4. Afterward, we evaluate the PSD for the collective, yaw, and tilt components of the signals, and the results are displayed in Sect. 12. The plot highlights the absence of 1P and 2P components (as observed in the rotating frame of reference; see Fig. 11) in the non-rotating frame of reference, in line with Sect. 2.4. Below 0.1 Hz, a good match between the collective and tilt components is observed, but the yaw component of the corrected inflow wind speed (ucor,yaw) is slightly underestimated. Furthermore, the 3P component of ucor,k (with k{col, yaw, tilt}) in the non-rotating frame of reference, which is the contribution of 2P and 4P from the rotating frame of reference, is likewise underestimated in all three components.

Figure 13Coherences (γ2) between the blade effective wind speeds (ubeff,k) and the corrected inflow wind speeds (ucor,k) in the non-rotating frame of reference, with k{col, yaw, tilt}.

Figure 13 reveals a good coherence at the frequencies where the power is concentrated, i.e. below 0.1 Hz, and at 3P and 6P. Additionally, the plots disclose the declining coherence with increasing frequency; i.e. higher coherence is achieved at 0P than at 3P; the same could be implied between 3P and 6P. With Fig. 12 highlighting the low-power content of the signals between 0P and 3P and between 3P and 6P, low coherences are similarly seen at the same frequencies in Fig. 13.

Figure 14The identified disturbance measurement transfer functions (Gℓ,k(jω)). The dashed–dotted lines indicate the estimated nominal disturbance measurement models (Gn,k(jω)). The dashed lines show the sum of the estimated nominal disturbance measurement models and uncertainty weights (${G}_{\mathrm{n},k}\left(j\mathit{\omega }\right)+{w}_{\mathrm{\ell },k}\left(j\phantom{\rule{0.125em}{0ex}}\mathit{\omega }\right)$), where k{col, yaw, tilt}.

Furthermore, we determine the disturbance measurement models (Gℓ,k(jω)), the nominal disturbance measurement models (Gn,k(jω)), and the measurement uncertainty weights (wℓ,k(jω)), shown in Fig. 14, which can be incorporated into the feedback–feedforward individual pitch control development and analysis. This case is labelled as C1 in Fig. 10. Figure 14 shows that this case only covers very small gain variations. The figure highlights that the mean value of the corrected inflow wind speed measurement is slightly underpredicted on the collective and yaw components, where the low-frequency gain is below 1, and is slightly overpredicted on the tilt component, where the low-frequency gain is above 1.

### 3.3.2 Uncertainties around the no-induction telescope parameters (C2)

In this section, we investigate the impact on the uncertainty weights when the telescope parameters cannot be selected as defined for the no-induction case, but are close to these values. We carried out simulations involving a discrete set of sampled values for the focus distance, radial position of the telescope along the blade radii, and orientation angles of the telescope. The identified disturbance measurement transfer functions (Gℓ,k(jω)) for the discrete set of sampled values are shown as overlapping blue lines in Fig. 15. The plot underscores that the disturbance measurement transfer functions are influenced by the telescope parameters. The low-frequency gain variation is different at each of the three components, which is also seen in Fig. 10, where it is labelled as C2. The highest low-frequency gain variation is observed on the tilt component.

Figure 15The identified disturbance measurement transfer functions (Gℓ,k(jω)) for a discrete set of sampled telescope parameters. The dashed–dotted lines indicate the estimated nominal disturbance measurement models (Gn,k(jω)). The dashed lines show the sum of the estimated nominal disturbance measurement models and uncertainty weights (${G}_{\mathrm{n},k}\left(j\mathit{\omega }\right)+{w}_{\mathrm{\ell },k}\left(j\mathit{\omega }\right)$), where k{col, yaw, tilt}.

### 3.3.3 Optimal focus distance and available preview time (C3)

We determine the preview time in accordance with Sect. 2.7. We keep the telescope parameters constant as defined in Sect. 3.3.1, except for the focus distance, which is allowed to vary between 20.2 and 30 m. We determine a preview time of 1.9 s for all the focus distances, which is slightly higher than the initially calculated value of 1.7 s in Sect. 3.3.1.

This case is denoted as C3 in Fig. 10, and that figure highlights that there is a smaller low-frequency gain variation for this case compared to the previous case (C2).

### 3.3.4 Telescope position along the blade span (C4)

proposed that a blade-mounted lidar placed at 70 % of the blade radius is most suitable for feedforward control input. We assess in this subsection whether placing the blade-mounted lidar at 70 % (≈44 m) of the blade radius would result in the maximum of the objective function in Eq. (21). We set fmax in Eq. (21) as 0.1 Hz, while we maintain a focus distance of 22.2 m.

We find that the telescope placed at a radial position of 46 m leads to the maximum value of the objective function in Eq. (21); in other words, the telescope positioned at a radial position of 46 m results in the highest coherence between the blade effective and the corrected inflow (ucor,k) wind speeds. This corresponds to 73 % of the blade span. The value found is quite close to the findings of . Varying the telescope radial position in a fairly small range (42–47 m) results in a higher low-frequency uncertainty in the tilt component than in the collective and yaw components (see C4 in Fig. 10). In this case, at the yaw and tilt components, the low-frequency gain variation is higher than in C3 but still smaller than in C2.

### 3.3.5 Telescope orientation (C5)

In this section, we evaluate whether the initially selected telescope orientation angles (Φℓ,i and Γℓ,i, with i=1, 2, 3) would result in a maximized objective function in Eq. (21). For this purpose, we fixed the telescope parameters as described in Sect. 3.3.1, with the exception of the orientation angles (Φℓ,i and Γℓ,i). The two angles are changed around the initially selected values. We simulate the lidar measurements with each new set of parameters.

Figure 16Optimal angular orientation of the telescope. The maximum frequency (fmax) in the objective function (Jlp) is set at 0.1 Hz. The blue diamond marks the parameters initially chosen; the red dot marks the maximum point of the function from Eq. (21).

We determine the optimal orientation of the telescope in Fig. 16 based on the objective function in Eq. (21). In the plot, the blue diamond marks the initial telescope orientation based on the no-induction calculation, where ${\mathrm{\Phi }}_{\mathrm{\ell },i}=-\mathrm{3.7}$ and ${\mathrm{\Gamma }}_{\mathrm{\ell },i}=\mathrm{7.0}$. The red dot indicates the obtained optimal value, where ${\mathrm{\Phi }}_{\mathrm{\ell },i}=-\mathrm{2.6}$ and ${\mathrm{\Gamma }}_{\mathrm{\ell },i}=\mathrm{8.1}$, which is only marginally different from the no-induction values.

The identified Gℓ,k low-frequency (DC) gain upper   (${G}_{\mathrm{n},k}\left(j\phantom{\rule{0.125em}{0ex}}\mathit{\omega }\right)+{w}_{\mathrm{\ell },k}\left(j\mathit{\omega }\right)$) and lower (Gn,k) bounds are labelled as C5 in Fig. 10. This case results in similar low-frequency gain variations as C2 at the yaw and tilt components; however, it has a smaller gain variation in the collective component than C2.

### 3.3.6 Telescope orientation misalignment (C6)

In this subsection, transfer functions (Gℓ,k) from the blade effective wind speeds (ubeff,k) to the corrected lidar-based inflow wind speeds (ucor,k) are identified for the cases where one or two of the telescopes have been aligned differently relative to their values in the no-induction case (C1). Such cases could occur, for example, during telescope installation. Initially, we assume this misalignment is unknown but detectable to allow for a correction of the lidar-based inflow wind speed measurement. To simulate these cases, we fixed the telescope parameters as described in Sect. 3.3.1, except for the orientation angles of Φℓ,i and Γℓ,i of the telescopes mounted on the second and third blades. The angular values are changed around the no-induction values by ±5 (${\mathrm{\Phi }}_{\mathrm{\ell },i}={\mathrm{\Phi }}_{\mathrm{\ell },\mathrm{1}}±\mathrm{5}$ and ${\mathrm{\Gamma }}_{\mathrm{\ell },i}={\mathrm{\Gamma }}_{\mathrm{\ell },\mathrm{1}}±\mathrm{5}$ for i=2, 3) as follows. First, the values are changed only for the telescope mounted on the second blade, then for the telescopes mounted on both the second and third blades.

Figure 17Identified transfer functions (Gℓ,k) from the blade effective wind speeds (ubeff,k) to the corrected lidar-based inflow wind speeds (ucor,k) for the discrete set of sampled telescope parameters with unknown and known telescope orientation misalignment, where k{col, yaw, tilt}.

We evaluate such a set-up via simulations. Figure 17a displays the identified transfer functions (Gℓ,k) from the blade effective wind speeds (ubeff,k) to the corrected lidar-based inflow wind speeds (ucor,k). Figure 17b reveals a 1P peak at the collective component and 1P and 2P peaks at the yaw and tilt components. As shown in Sect. 3.1, adding a phase shift of 1 to the 1P harmonic and reducing the DC offset for one of the signals in the rotating frame of reference results in such undesired higher-harmonic peaks at the collective, yaw, and tilt components in the non-rotating frame of reference. Figure 17b underlines that, by assuming that the misalignment angles are identifiable and that the lidar-based inflow wind speed measurement is corrected accordingly, the undesired peak at 1P is reduced by a factor of 10, although existent on all the components.

4 Discussion

We have shown that the determined telescope parameters with assumptions of rigid blades, constant rotor speed and pitch angle, absence of induction, and Taylor's frozen turbulence hypothesis provide a good trade-off between simplicity and accuracy (see C1 in Fig. 10). First, the low-frequency gains of the identified disturbance measurement models (Gℓ,k(jω)) only have small absolute deviations from 1, which are found to be 3 %, 2.5 %, and 3.1 % for collective, yaw, and tilt components, respectively. Second, the optimal telescope parameters in C4 and C5, which maximize a cost function based on the coherence between the blade effective (ubeff,k) and the corrected inflow (ucor,k) wind speeds, are close to the telescope parameters in C1. Such a small deviation is expected with respect to the assumptions we made during the calculation of the values for the no-induction case (see Sect. 3.3.1).

By evaluating the cross-correlation between the blade effective (ubeff,k) and the corrected inflow (ucor,k) wind speeds for a discrete set of sampled values of the focus distance in Sect. 3.3.3, we found that the preview time is constant for all the selected focus distances. It is closely coupled to the time needed for blade i−1 to reach the position of blade i, i.e. 120 azimuth angle change. For example, by considering laminar inflow with wind shear, no matter what the focus distance is, the delay time between the corrected inflow wind speed from blade 1 and the blade effective wind speed from blade 3 will always be the same, which is the time needed for blade i−1 to reach the position of blade i. If the focus distance has changed, the ϕ in the MBC transformation also has to be changed; furthermore, the control signal should be delayed accordingly. Note that control development must proceed with sufficient attention so as to ensure that the feedforward controller does not result in a higher time delay than the available preview time. For example, a feedforward controller with a crossover frequency of 0.1 Hz may result in higher time delay compared to that with a crossover frequency of 0.2 Hz . With this, we want to point out that the feedforward controller crossover frequency and the focus distance are coupled. Hence, defining the former typically leads to a minimal selectable focus distance.

As stated above, the lidar and telescope parameters based on the assumptions we made in Sect. 3.3.1 provide a good trade-off between simplicity and accuracy. They are close to the optimal parameters we found for the discrete set of sampled values of the focus distance, the radial position of the telescope along the blade, and the orientation angles of the telescope in Sect. 3.3.33.3.5. Nevertheless, this is not the case for the preview time; the available preview time is slightly increased from 1.7 to 1.9 s, as we demonstrated in Sect. 3.3.3. This could be due to the assumptions we made: (a) the blades are rigid; (b) constant rotor speed and blade pitch angle; (c) Taylor's frozen turbulence hypothesis holds; and (d) the induction effect is absent during our calculation in Sect. 3.3.1. Furthermore, the signals were sampled with a sampling time of 0.2 s, which also imposes limitations on the resolution of the preview time. Note that LES simulations with lower sampling time are resource and time expensive. The crossover frequency of the feedforward controller affects the time delay. With a higher preview time available, we can select a lower crossover frequency. This understanding gives us more room during the feedback–feedforward control development. The available preview time could be determined online in field tests and used to delay the feedforward control signal accordingly. This can be done online by, for example, storing 10 min of blade effective (ubeff,k) and corrected inflow (ucor,k) wind speed measurements and evaluating the cross-correlation between them.

We found that the blade-mounted lidar placed at the 73 % span of the blade radius results in the best coherence between the corrected inflow wind speed and the blade effective wind speed. This finding is close to the value (70 % of the blade radius) found by for a blade-mounted lidar and for a hub-mounted lidar system.

Any unknown orientation angle misalignment for one of the telescopes leads to an unknown contribution of the rotational speed to the lidar-based line-of-sight wind speed measurement. This is the reason why the low-frequency gain can vary between 0.7 and 1.2. Nevertheless, this can be reduced to a low-frequency gain variation between 0.96 and 0.98 by assuming that we are able to detect the angular offset. By detecting the angular offset, we are able to better estimate what is the mean value of the blade effective wind speed, and the resulting Gℓ,k low-frequency gain lower bound is 0.96, which is very close to 1. In addition, an undetected misalignment of the telescope orientation angle results in a phase shift of the 1P harmonic and a reduction or increase in the DC offset of the signal in the rotating frame of reference. This subsequently leads to undesired peaks at 1P and 2P frequencies at the collective, yaw, and tilt components in the non-rotating frame of reference. Assuming the angular offsets to be known, we can reduce 1P and 2P peaks by a factor of 10 for both yaw and tilt components, but we cannot completely eliminate them. Thus, the question as to whether robust stability and performance can be ensured with such peaks still remains. To avoid such a peak, the telescopes need to be well aligned with each other and the blade segment orientation angles and linear velocities should be measured well. We showed in Sects. 3.1 and 3.3.6 that an unknown orientation angle misalignment leads to a 1P peak at the yaw and tilt components in the frequency domain. Therefore, the orientation angles of the telescope can be identified by formulating an optimization problem, whose main objective is to minimize the 1P peaks at the yaw and tilt components, with the orientation angles of the telescopes as the decision variables.

The nominal measurement transfer functions and uncertainty weights identified for the no-induction case can be directly included in robust feedback–feedforward individual pitch and trailing edge flap control development to guarantee robust stability and performance. However, this would be a very optimistic approach, as we considered only one reference wind turbine with a single inflow wind condition, and we need to assess how the measurement uncertainties change for other wind turbines with different wind speeds, turbulence intensities, yaw misalignments, etc. The nominal measurement transfer functions and uncertainty weights found in Sect. 3.3.2 might cover these cases and may be better for robust control development. C2 covers a wide range of telescope parameter variations; hence, if for some reason one or more lidar and telescope parameters cannot be selected as for the no-induction case but are close to these values, the established transfer functions from C2 can be used for robust feedback–feedforward control development. In addition, C2 also covers the situations where the mean blade pitch angle is increased or decreased due to the wind turbine operating at a different point. The final selection of the nominal measurement transfer functions and uncertainty weights depends on whether the lidar and telescope parameters are varied dynamically with the operating points of the wind turbine and the wind speed or are kept constant over the entire operating range. Nevertheless, this may result in a conservative feedforward controller, thus limiting the benefits of the lidar system.

The methodology we presented in this paper can be applied in identifying the uncertainty weight for higher-harmonics control development; i.e. selecting nh as 2 in Eq. (8) can be used to identify the uncertainty weight for the controller to mitigate 2P dynamic blade loads.

5 Conclusions

Our paper has aimed to identify the nominal measurement transfer functions and model the uncertainties in blade-mounted lidar measurements as a frequency-dependent uncertainty weight for inclusion into the feedback–feedforward individual pitch and trailing edge flap control development.

We found that the preview time with the lidar mounted on the blade is more linked to the time it takes the previous blade to reach the position of the blade from which the measurement took place rather than the focus distance. For a given focus distance, the preview time can be estimated online; hence, the feedforward control signal can be delayed accordingly. While the selected focus distance should provide sufficient preview time, it is desirable that the time delay introduced by the feedforward controller and actuators be eliminated. This sets the lower limit for the selectable focus distance.

Accordingly, we introduced a simple method, based on steady-state data, to calculate the telescope and lidar parameters. Nevertheless, we showed in a large-eddy simulation that such an approach provides a good trade-off between an efficient determination of the telescope parameters and accurate inflow wind speed measurement. The low-frequency gains of identified disturbance measurement transfer functions had small absolute deviations from 1, which were due to wind evolution, the cyclops dilemma, using a single-point measurement to estimate the blade effective wind speed, and the assumptions we made to correct the measurements. The nominal measurement transfer functions and uncertainty weights, as we have identified in this paper for several cases, can be directly included in the robust feedback–feedforward individual pitch and trailing edge flap control development to ensure robust stability and performance. However, to prevent the transfer functions (G) from the blade effective wind speeds (ubeff) to the corrected lidar-based inflow wind speeds (ucor) from having a large high-frequency gain at 1P and 2P in the non-rotating frame of reference, the telescopes must be well aligned with each other and the blade segment orientation angles and linear velocities should be measured well.

Data availability
Data availability.

The presented results are archived by the University of Oldenburg and can be obtained by contacting the corresponding author.

Author contributions
Author contributions.

RU carried out the research and prepared the article. VP, LYP, and MK gave guidelines for the development of the methodology and had a supervisory function. All authors discussed the results and helped to improve the article.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

The support of a fellowship from the Hanse–Wissenschaftskolleg in Delmenhorst (HWK) is likewise gratefully acknowledged.

Financial support
Financial support.

This work has been partly funded by the Federal Ministry for Economic Affairs and Energy according to a resolution by the German Federal Parliament (projects DFWind 0325936C and SmartBlades2 0324032D).

Review statement
Review statement.

This paper was edited by Jakob Mann and reviewed by two anonymous referees.

References

Barlas, T. K., van Wingerden, J.-W., Hulskamp, A. W., van Kuik, G. A. M., and Bersee, H. E. N.: Smart dynamic rotor control using active flaps on a small-scale wind turbine: aeroelastic modeling and comparison with wind tunnel measurements, Wind Energy, 16, 1287–1301, https://doi.org/10.1002/we.1560, 2013. a

Berg, J. C., Resor, B. R., Paquette, J. A., and White, J. R.: SMART wind turbine rotor. Design and field test, Tech. rep., Sandia National Laboratories, Sandia National Laboratories, Albuquerque, New Mexico, USA, available at: https://www.energy.gov/sites/prod/files/2014/03/f14/smart_wind_turbine_design.pdf (last access: 29 November 2019), 2014. a

Bergami, L. and Poulsen, N. K.: A smart rotor configuration with linear quadratic control of adaptive trailing edge flaps for active load alleviation, Wind Energy, 18, 625–641, https://doi.org/10.1002/we.1716, 2015. a

Bossanyi, E. A.: Individual blade pitch control for load reduction, Wind Energy, 6, 119–128, https://doi.org/10.1002/we.76, 2003. a, b, c

Bossanyi, E. A.: Further load reductions with individual pitch control, Wind Energy, 8, 481–485, https://doi.org/10.1002/we.166, 2005. a

Bossanyi, E. A.: Un-freezing the turbulence: application to LiDAR-assisted wind turbine control, IET Renewable Power Generation, 7, 321–329, https://doi.org/10.1049/iet-rpg.2012.0260, 2013. a, b, c, d, e, f

Bossanyi, E. A., Kumar, A., and Hugues-Salas, O.: Wind turbine control applications of turbine-mounted LIDAR, J. Phys.: Conf. Ser., 555, 012011, https://doi.org/10.1088/1742-6596/555/1/012011, 2014. a, b

Bromm, M., Vollmer, L., and Kühn, M.: Numerical investigation of wind turbine wake development in directionally sheared inflow, Wind Energy, 20, 381–395, https://doi.org/10.1002/we.2010, 2017. a

Castaignet, D., Barlas, T., Buhl, T., Poulsen, N. K., Wedel-Heinen, J. J., Olesen, N. A., Bak, C., and Kim, T.: Full-scale test of trailing edge flaps on a Vestas V27 wind turbine: active load reduction and system identification, Wind Energy, 17, 549–564, https://doi.org/10.1002/we.1589, 2014. a, b

Dunne, F. and Pao, L. Y.: Optimal blade pitch control with realistic preview wind measurements, Wind Energy, 19, 2153–2169, https://doi.org/10.1002/we.1973, 2016. a

He, K., Qi, L., Zheng, L., and Chen, Y.: Combined pitch and trailing edge flap control for load mitigation of wind turbines, Energies, 11, 2519, https://doi.org/10.3390/en11102519, 2018. a

Jonkman, J. M. and Buhl, M. L.: FAST user's guide, Tech. rep., National Renewable Energy Laboratory, Golden, Colorado, 2005. a

Jonkman, J. M., Butterfield, S., Musial, W., and Scott, G.: Definition of a 5-MW reference wind turbine for offshore system development, Tech. rep., National Renewable Energy Laboratory, Golden, Colorado, 2009. a, b

Kapp, S.: Lidar-based reconstruction of wind fields and application for wind turbine control, PhD thesis, Carl von Ossietzky University of Oldenburg, Oldenburg, available at: http://oops.uni-oldenburg.de/3210/1/kaplid17.pdf (last access: 29 November 2019), 2017. a, b

Laks, J., Simley, E., and Pao, L.: A spectral model for evaluating the effect of wind evolution on wind turbine preview control, in: American Control Conference (ACC), 17–19 June 2013, Washington, DC, USA, 3673–3679, https://doi.org/10.1109/ACC.2013.6580400, 2013. a, b

Lu, Q., Bowyer, R., and Jones, B.: Analysis and design of Coleman transform-based individual pitch controllers for wind-turbine load reduction, Wind Energy, 18, 1451–1468, https://doi.org/10.1002/we.1769, 2015. a, b

Maronga, B., Gryschka, M., Heinze, R., Hoffmann, F., Kanani-Sühring, F., Keck, M., Ketelsen, K., Letzel, M. O., Sühring, M., and Raasch, S.: The Parallelized Large-Eddy Simulation Model (PALM) version 4.0 for atmospheric and oceanic flows: Model formulation, recent developments, and future perspectives, Geosci. Model Dev. 8, 2515–2551, https://doi.org/10.5194/gmd-8-2515-2015, 2015. a

Marten, D., Bartholomay, S., Pechlivanoglou, G., Nayeri, C., Paschereit, C. O., Fischer, A., and Lutz, T.: Numerical and experimental investigation of trailing edge flap performance on a model wind turbine, in: 2018 Wind Energy Symposium, AIAA SciTech Forum, American Institute of Aeronautics and Astronautics, 8–12 January 2018, Kissimmee, Florida, USA, https://doi.org/10.2514/6.2018-1246, 2018. a

MATLAB: version 9.4.0 (R2018a), MathWorks, available at: https://www.mathworks.com/products/matlab.html (last access: 1 April 2019), 2018. a

Mikkelsen, T., Siggaard Knudsen, S., Sjöholm, M., Angelou, N., and Tegtmeier, A.: WindScanner.eu – a new remote sensing research infrastructure for on- and offshore wind energy, in: International Conference on Wind Energy: Materials, Engineering, and Policies, 22–23 November 2012, Hyderabad Campus, India, 2012. a

Pechlivanoglou, G.: Passive and active flow control solutions for wind turbine blades, PhD thesis, Technical University of Berlin, Berlin, https://doi.org/10.14279/depositonce-3487, 2013. a

Petrović, V., Jelavić, M., and Baotić, M.: Advanced control algorithms for reduction of wind turbine structural loads, Renewable Energy, 76, 418–431, https://doi.org/10.1016/j.renene.2014.11.051, 2015. a

Schlipf, D., Schlipf, D. J., and Kühn, M.: Nonlinear model predictive control of wind turbines using LIDAR, Wind Energy, 16, 1107–1129, https://doi.org/10.1002/we.1533, 2013. a

Simley, E., Pao, L. Y., Frehlich, R., Jonkman, B., and Kelley, N.: Analysis of light detection and ranging wind speed measurements for wind turbine control, Wind Energy, 17, 413–433, https://doi.org/10.1002/we.1584, 2014a.  a, b, c, d, e, f

Simley, E., Pao, L. Y., Gebraad, P., and Churchfield, M.: Investigation of the impact of the upstream induction zone on LiDAR measurement accuracy for wind turbine control applications using Large-Eddy simulation, J. Phys.: Conf. Ser., 524, 12003, https://doi.org/10.1088/1742-6596/524/1/012003, 2014b. a

Simley, E., Angelou, N., Mikkelsen, T., Sjöholm, M., Mann, J., and Pao, L. Y.: Characterization of wind velocities in the upstream induction zone of a wind turbine using scanning continuous-wave lidars, J. Renew. Sustain. Energy, 8, 013301, https://doi.org/10.1063/1.4940025, 2016. a

Skogestad, S. and Postlethwaite, I.: Multivariable feedback control: Analysis and design, 2nd Edn., John Wiley & Sons, Chichester, UK, 2005. a

Taylor, G. I.: The spectrum of turbulence, P. Roy. Soc. A, 164, 476–490, https://doi.org/10.1098/rspa.1938.0032, 1938. a, b

Ungurán, R. and Kühn, M.: Combined individual pitch and trailing edge flap control for structural load alleviation of wind turbines, in: American Control Conference, 6–8 July 2016, Boston, USA, 2307–2313, https://doi.org/10.1109/ACC.2016.7525262, 2016. a

Ungurán, R., Petrović, V., Pao, L. Y., and Kühn, M.: Performance evaluation of a blade-mounted LiDAR with dynamic versus fixed parameters through feedback-feedforward individual pitch and trailing edge flap control, J. Phys.: Conf. Ser., 1037, 032004, https://doi.org/10.1088/1742-6596/1037/3/032004, 2018. a, b, c

Ungurán, R., Boersma, S., Petrović, V., van Wingerden, J.-W., Pao, L. Y., and Kühn, M.: Feedback-feedforward individual pitch control design for wind turbines with uncertain measurements, in: American Control Conference, 10–12 July 2019, Philadelphia, USA, 2019. a, b, c, d, e

van Engelen, T. G.: Design model and load reduction assessment for multi-rotational mode individual pitch control (higher harmonics control), Tech. Rep. ECN-RX–06-068, Energy Research Centre of the Netherlands, Petten, the Netherlands, 2006. a, b, c

van Wingerden, J.-W., Hulskamp, A., Barlas, T., Houtzager, I., Bersee, H., van Kuik, G., and Verhaegen, M.: Two-degree-of-freedom active vibration control of a prototyped “smart” rotor, IEEE T. Control Syst. Technol., 19, 284–296, https://doi.org/10.1109/TCST.2010.2051810, 2011. a

Zhang, W., Bai, X., Wang, Y., Han, Y., and Hu, Y.: Optimization of sizing parameters and multi-objective control of trailing edge flaps on a smart rotor, Renew. Energy, 129, 75–91, https://doi.org/10.1016/j.renene.2018.05.091, 2018. a