<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JFCMV</journal-id><journal-title-group><journal-title>Journal of Flow Control, Measurement &amp; Visualization</journal-title></journal-title-group><issn pub-type="epub">2329-3322</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jfcmv.2017.51001</article-id><article-id pub-id-type="publisher-id">JFCMV-74028</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  On the Improvement of Wind Power Predictions Based on Terrain Characteristics and Measurements of the Annual Energy Production
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jonathan</surname><given-names>Dorval</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Christian</surname><given-names>Masson</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yves</surname><given-names>Gagnon</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Mechanical Engineering, &amp;amp;EACUTE;cole de technologie sup&amp;amp;ecirc;rieure, Montreal, Canada</addr-line></aff><aff id="aff2"><addr-line>Department of Sciences, Universit&amp;amp;ecirc; de Moncton, Edmundston, Canada</addr-line></aff><pub-date pub-type="epub"><day>09</day><month>02</month><year>2017</year></pub-date><volume>05</volume><issue>01</issue><fpage>1</fpage><lpage>20</lpage><history><date date-type="received"><day>November</day>	<month>17,</month>	<year>2016</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>February</month>	<year>6,</year>	</date><date date-type="accepted"><day>February</day>	<month>9,</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The assessment of the production capacity of wind farms is a crucial step in wind farm design processes, where a poor assessment can cause significant economic losses. Data from Canadian wind farms benefiting from national production incentive programs show that wind farms are typically characterized by an overestimation of the production capacity. In this context, a study has been done to provide insight on the origin of the discrepancies between the energy production estimates and the measured energy generation, and to develop a method to reduce these discrepancies. To this end, the WAsP and MS-Micro models have been studied. Besides the wind speed measurements, topography indices have been developed to identify the influence of the various characteristics of the site on the error in the annual energy production (
  <em>AEP</em>). In addition, roughness classes have been created, including a reference roughness and a roughness complexity. The indices have also allowed establishing correlations and developing equations to evaluate the error based on the site characteristics and the positions of wind turbines on the measured annual energy production. An average reduction of up to 83% on the 
  <em>AEP</em> errors was obtained when the methodology was applied to five wind farms in Canada.
 
</p></abstract><kwd-group><kwd>Wind Energy</kwd><kwd> Wind Farm Assessment</kwd><kwd> WAsP</kwd><kwd> MS-Micro</kwd><kwd> Roughness Classes</kwd><kwd>  Roughness Complexity</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Similar to the rest of the world, the wind power industry in Canada has experienced significant growth over the last decade. The total Canadian wind capacity surpassed the 1000 MW mark in 2006 and it reached more than 11,000 MW in 2015 [<xref ref-type="bibr" rid="scirp.74028-ref1">1</xref>] . In order to assist in the growth of wind energy, the Government of Canada took several initiatives to accelerate the development of renewable energy in the country. One of the government’s initiatives has been in the form of wind power production incentive programs, notably the Wind Power Production Incentive program (WPPI) and the EcoEnergy program (EcoERP) [<xref ref-type="bibr" rid="scirp.74028-ref2">2</xref>] . These programs subsidize wind energy projects based on their estimated annual energy production.</p><p>A study from Natural Resources Canada on the annual energy production (AEP) of 59 projects funded by the WPPI and EcoERP programs [<xref ref-type="bibr" rid="scirp.74028-ref3">3</xref>] has shown that, in 2011, while a capacity factor of 34.2% was expected for the wind farms benefitting from the programs, the effective capacity factor was only 29.9%. This difference means an underperformance of about 12.5% and represents significant losses in economic terms. This behavior of underperformance of wind farms has been observed year after year.</p><p>Besides the daily and seasonal variations, the wind resource on a specific site is influenced by the physical characteristics of the terrain. Carvalho et al. [<xref ref-type="bibr" rid="scirp.74028-ref4">4</xref>] showed that these characteristics are significantly enhanced when dealing with a complex terrain, where the wind predictions are more difficult to achieve. Thus, the estimation of the wind resource, and more specifically the annual energy production (AEP) using the wind speed distribution of the IEC 61400-12-1 standard, might lead to significant uncertainties [<xref ref-type="bibr" rid="scirp.74028-ref5">5</xref>] .</p><p>Studies have been published in the literature to assess the feasibility of wind energy in various conditions. For example, Rehman et al. [<xref ref-type="bibr" rid="scirp.74028-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.74028-ref7">7</xref>] presented wind resource assessments, along with the design and the economic assessment of a 20 MW wind power plant, using 2 MW turbines, located in the Eastern region of Saudi Arabia. Specific to wind farms, Touani et al. [<xref ref-type="bibr" rid="scirp.74028-ref8">8</xref>] developed of a Wind to Power Model for wind farm power production forecasting, where wind farm power curves are used to convert the predicted wind speeds to predicted power productions of wind farms. In their work, they have however not considered the site characteristics of the wind farms.</p><p>Kwon [<xref ref-type="bibr" rid="scirp.74028-ref9">9</xref>] proposed a methodology to better assess the AEP of wind farms. While making some simplifying assumptions, this methodology introduces a new approach for long-term corrections, the vertical extrapolation of the wind distribution and the power performance of the wind turbines by taking into account the uncertainties. However, the effects of the terrain on the AEP are not considered.</p><p>It is known that the land cover (roughness) and terrain height variations (orography) influence the wind. When these effects are not well interpreted, the horizontal and vertical extrapolation of the wind speed can have significant uncertainties, thus providing erroneous AEP estimates.</p><p>The present study aims to better understand the effect of orography and roughness on the production capacity of wind farms and to develop a methodology to reduce the gap between the energy production estimates and the actual measured energy generation, and to develop a method to reduce these discrepancies.</p><p>The models studied are WAsP [<xref ref-type="bibr" rid="scirp.74028-ref10">10</xref>] and MS-Micro [<xref ref-type="bibr" rid="scirp.74028-ref11">11</xref>] , which are linear models based on the analysis of the theory of flow over hills by Jackson and Hunt [<xref ref-type="bibr" rid="scirp.74028-ref12">12</xref>] . The simplifications of the Navier-Stoke equations in the models generate uncertainties, notably when flow separation occurs.</p><p>Bowen and Mortensen [<xref ref-type="bibr" rid="scirp.74028-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.74028-ref14">14</xref>] developed a methodology to evaluate the accuracy of the WAsP model in complex terrains. It has been shown [<xref ref-type="bibr" rid="scirp.74028-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.74028-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.74028-ref15">15</xref>] that a relationship, called the ruggedness index (ΔRIX), exists between the mean wind speed errors and a flow separation indicator.</p><p>Mortensen et al. [<xref ref-type="bibr" rid="scirp.74028-ref16">16</xref>] investigated the optimum configuration of the parameters for the RIX calculation and showed that there is also a relationship between the AEP errors and ΔRIX. A coefficient of determination (R<sup>2</sup>) of 97.7% was obtained between the ΔRIX and the AEP errors on five sites located in rugged mountainous areas of northern Portugal.</p><p>Berge et al. [<xref ref-type="bibr" rid="scirp.74028-ref17">17</xref>] also obtained good correlations between the ΔRIX and the mean wind speed errors on five mountainous sites in Norway by using WAsP. Moreover, they investigated the relationship between the wind speed errors and the differences in heights and horizontal distances. No correlation was visible between the mean wind speed errors and the horizontal distances and a weak correlation was observed between wind speed errors and the difference in heights.</p><p>Feng et al. [<xref ref-type="bibr" rid="scirp.74028-ref18">18</xref>] validated the ruggedness index methodology based on energy production with data from an existing wind farm. The energy production corrected was more accurate; however, AEP errors up to 9.4% were still observed. For their part, Jimenez et al. [<xref ref-type="bibr" rid="scirp.74028-ref19">19</xref>] investigated WAsP on two offshore, one onshore and three island sites located in northern Germany. Four sites correlated well with each other, with mean wind speed errors up to 2%, while the two other sites gave significant mean wind speed errors.</p><p>These studies show that methodologies to increase the accuracy of the WAsP model in complex terrains can give relatively good results. However, these methodologies do not remove all the errors in complex terrains. Moreover, it has been observed that significant AEP errors occur even on simple terrains, where the methodologies cannot be used.</p><p>Based on these studies, several indices have been developed to better identify the effect of the terrain on the AEP errors. According to these indices, a new methodology has been investigated which uses a multi-variable correlations to give more accurate results. Finally, the efficiency of this methodology will also be applied to the MS-Micro model.</p><p>Following a presentation of the measured production data available and the quality control performed on these data, the quantification of the AEP errors are presented. Thereafter, the indices created to characterize the sites in terms of roughness, slopes and the position of the wind turbines relative to the meteorological mast are introduced; this is followed by the correlations, based on the indices proposed, used to relate the AEP error and the site characteristics. Finally, the results of the analysis without improvement are presented and compared to the improvements obtained due to the methodology proposed in this work. A short conclusion and recommendations for future work ends the paper.</p></sec><sec id="s2"><title>2. Measured Production Data</title><p>The production data measured from five Canadian wind farms, with a global installed capacity of 350 MW, have been used to perform this study. For confidentiality reasons, the information that identifies these wind farms is not divulged.</p><p>Wind Farm I is characterized as a very simple site, with a relatively flat terrain and low roughness classes. Wind Farms II and III are complex sites with steep terrains and high roughness classes (notably forest). Finally, Wind Farms IV and V are simple sites with small slopes and high roughness classes.</p><p>The measured meteorological data, contour maps, roughness maps and wind farm layouts have been provided by the owners of the wind farms.</p><p>The meteorological data are ten minute averages measured at one 80 m meteorological (met) mast for every wind farm, while the corresponding production data cover one year of energy production for every wind farm, except for Wind Farm I where two years of meteorological and production data are available.</p><p>A quality control has been performed on the met mast measured data according to three main steps. First, the missing values and the values in the sectors that are disturbed by the tower shadow and by the wake of wind turbines in proximity to the met masts have been removed. The disturbed sector calculations have been done according to the relevant IEC standard [<xref ref-type="bibr" rid="scirp.74028-ref20">20</xref>] . Then, following the approach of the Wind Resource Assessment Handbook [<xref ref-type="bibr" rid="scirp.74028-ref21">21</xref>] , data which did not satisfy one or more criteria (see <xref ref-type="table" rid="table1">Table 1</xref>) have been flagged. Finally, a visual inspection of the flagged data has been done to identify and remove time series where abnormal operations of the instruments have occurred.</p><p>A quality control has also been performed on wind turbine data. The first step in the quality control of wind turbine data is to eliminate the missing values, the values in the perturbed sectors by the wake of the wind turbines in proximity to the met mast, and the time series which are not available at the met mast. Then, the time series resulting in active power values and/or reactive power values equal to or below zero have been removed.</p><p>Since it is not common to use a met mast as far as 8 km from any wind turbine position when performing wind assessments, this study is limited to wind turbines located within a radius of 8 km from the corresponding wind farm met mast. Moreover, it has been observed that wind turbines located at a distance greater than 8 km have significant AEP errors (see Section 6.3).</p></sec><sec id="s3"><title>3. Quantification of the AEP Errors</title><p>In order to have the largest amount of unperturbed data in the various analyses, each wind turbine has been analyzed individually. Thus 197 analyses have been</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Met mast quality control criteria for the measured data</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Categories</th><th align="center" valign="middle" >Parameters</th><th align="center" valign="middle" >Validation Criteria</th></tr></thead><tr><td align="center" valign="middle" >Icing</td><td align="center" valign="middle" >Temperature and, Relative Humidity and, Wind direction standard deviation or Mean Wind Speed</td><td align="center" valign="middle" >&lt;2˚C &gt;80% =0˚ =0 m/s</td></tr><tr><td align="center" valign="middle"  rowspan="5"  >Range test criteria at 80 m</td><td align="center" valign="middle" >Mean wind speed</td><td align="center" valign="middle" >0 m/s &lt; value &lt; 25 m/s</td></tr><tr><td align="center" valign="middle" >Wind speed standard deviation</td><td align="center" valign="middle" >0 m/s &lt; value &lt; 3 m/s</td></tr><tr><td align="center" valign="middle" >Maximum gust of wind speed</td><td align="center" valign="middle" >0 m/s &lt; value &lt; 30 m/s</td></tr><tr><td align="center" valign="middle" >Wind direction standard deviation</td><td align="center" valign="middle" >1˚ &lt; value &lt; 75</td></tr><tr><td align="center" valign="middle" >Max gust vs. mean wind speed</td><td align="center" valign="middle" >Max gust ≤2.5* Mean WS</td></tr><tr><td align="center" valign="middle"  rowspan="3"  >Relational test criteria between 80 m and 50 m</td><td align="center" valign="middle" >Δ Mean wind speed</td><td align="center" valign="middle" >≤2 m/s</td></tr><tr><td align="center" valign="middle" >Δ Mean wind speed by day</td><td align="center" valign="middle" >≤5 m/s</td></tr><tr><td align="center" valign="middle" >Δ Mean wind direction</td><td align="center" valign="middle" >≤20˚</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Relational test criteria between 80 m and 30 m</td><td align="center" valign="middle" >Δ Mean wind speed</td><td align="center" valign="middle" >≤4 m/s</td></tr><tr><td align="center" valign="middle" >Δ Mean wind speed by day</td><td align="center" valign="middle" >≤7.5 m/s</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Trend test criteria</td><td align="center" valign="middle" >Mean wind speed</td><td align="center" valign="middle" >&lt;5 m/s</td></tr><tr><td align="center" valign="middle" >Temperature</td><td align="center" valign="middle" >≤5˚C</td></tr></tbody></table></table-wrap><p>done for every numerical model (WAsP and Ms-Micro), which is the number of wind turbines available for this study.</p><p>To evaluate the accuracy of the analyses, the AEP errors between the measured values and the estimated values from the models have been calculated for each wind turbine:</p><disp-formula id="scirp.74028-formula33"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x2.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x3.png" xlink:type="simple"/></inline-formula> is the AEP calculated with the models and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x4.png" xlink:type="simple"/></inline-formula> is the AEP measured as the output energy of the wind turbine.</p><p>Since the AEP estimated by the models do not take into account the energy losses, most of these loss events (availability losses, wake losses and icing losses) have been removed from the database when calculating the<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x5.png" xlink:type="simple"/></inline-formula>. Also, for the same reasons, the data from wind turbines that were not operating in normal conditions (e.g. failure of a major component of the wind turbine) have been removed from the analysis. The method used to identify these wind turbines is presented in Section 5. The effects of the electric losses in the analysis have been minimized by taking the energy measurements directly at the wind turbine output. Finally, it is assumed that any other losses, which are not considered in this study, are negligible and constant for all the wind turbines analyzed.</p></sec><sec id="s4"><title>4. Index Classification</title><p>In order to evaluate the effect of the different site characteristics on the AEP, an index classification has been developed. The indices have been divided in three categories: an index related to the position of the wind turbines relative to the met mast, an index related to the orography characteristics and an index related to the roughness of the terrain.</p><p>The site characteristics (orography, roughness) have been taken from the vector maps provided by the wind farm owners. Thus, it is assumed that the maps provided are representatives of the site characteristics of the wind farms studied.</p><p>The position of the wind turbines relative to the met mast has been characterized by the horizontal distance between the wind turbines and the met mast, along with the difference of height between the hub of the turbines and the top measurement height at the met mast of the wind farm.</p><p>The terrain characteristics have been evaluated on ten degree sectors over a radius corresponding to 20 times the rotor diameter (D) and it includes two ruggedness parameters (RIX). By definition, the RIX is the fractional extent of the surrounding terrain which is steeper than a certain critical slope [<xref ref-type="bibr" rid="scirp.74028-ref16">16</xref>] .</p><p>The first parameter is called the drop index, RIX<sub>10</sub>, where the critical slope is 10%. This index has been based on a classification of the IEC [<xref ref-type="bibr" rid="scirp.74028-ref22">22</xref>] . <xref ref-type="table" rid="table2">Table 2</xref> shows the different classes of the drop index.</p><p>The second parameter is the flow separation index, which is already used by WAsP to correct the AEP [<xref ref-type="bibr" rid="scirp.74028-ref23">23</xref>] . This method is based on the fact that a linear relationship exists between the difference of flow separation between the met mast area and the wind turbine area, and the AEP error [<xref ref-type="bibr" rid="scirp.74028-ref16">16</xref>] .</p><p>The minimum critical slope before separation of the flow has been the subject of several studies. Wood [<xref ref-type="bibr" rid="scirp.74028-ref24">24</xref>] summarized the critical slopes given by Grant and Mason [<xref ref-type="bibr" rid="scirp.74028-ref25">25</xref>] from various numerical studies of two dimensional flows. <xref ref-type="table" rid="table3">Table 3</xref> provides a summary.</p><p>Taking into account the different values for the minimum critical slope for the separation of the flow [<xref ref-type="bibr" rid="scirp.74028-ref25">25</xref>] , a critical slope of 40% has been used in this work.</p><p>Using the roughness classes in the European Wind Atlas [<xref ref-type="bibr" rid="scirp.74028-ref26">26</xref>] , the land surrounding wind turbines can contain different types of roughness. However, the roughness complexity, denoted as the variation of roughness over an area of interest, cannot be properly quantified by the current roughness classes. To this</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Classes for the drop index RIX<sub>10</sub></title></caption><table><tbody><thead><tr><th align="center" valign="middle" >RIX<sub>10</sub> (%)</th><th align="center" valign="middle" >Terrain class</th></tr></thead><tr><td align="center" valign="middle" >0 ≤ RIX<sub>10</sub> &lt; 8</td><td align="center" valign="middle" >1</td></tr><tr><td align="center" valign="middle" >8 ≤ RIX<sub>10</sub> &lt; 16</td><td align="center" valign="middle" >2</td></tr><tr><td align="center" valign="middle" >16 ≤ RIX<sub>10</sub> &lt; 24</td><td align="center" valign="middle" >3</td></tr><tr><td align="center" valign="middle" >24 ≤ RIX<sub>10</sub></td><td align="center" valign="middle" >4</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Critical slope before separation of the flow</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Study number</th><th align="center" valign="middle" >Source</th><th align="center" valign="middle" >Critical Slope</th></tr></thead><tr><td align="center" valign="middle" >I</td><td align="center" valign="middle" >Newley [<xref ref-type="bibr" rid="scirp.74028-ref27">27</xref>]</td><td align="center" valign="middle" >43%</td></tr><tr><td align="center" valign="middle" >II</td><td align="center" valign="middle" >Newley [<xref ref-type="bibr" rid="scirp.74028-ref27">27</xref>]</td><td align="center" valign="middle" >40%</td></tr><tr><td align="center" valign="middle" >III</td><td align="center" valign="middle" >Taylor et al. [<xref ref-type="bibr" rid="scirp.74028-ref28">28</xref>]</td><td align="center" valign="middle" >31%</td></tr><tr><td align="center" valign="middle" >IV</td><td align="center" valign="middle" >Mason and King [<xref ref-type="bibr" rid="scirp.74028-ref29">29</xref>]</td><td align="center" valign="middle" >47%</td></tr></tbody></table></table-wrap><p>end, a new method of classification of the roughness has been developed to take into account the different roughness lengths and roughness variations over an area surrounding the wind turbines. This new classification has been described by the reference roughness length and the roughness complexity.</p><p>The reference roughness length is the average roughness length of the site and is obtained using the reference friction coefficient. From the equation of velocity distribution in a neutral boundary layer:</p><disp-formula id="scirp.74028-formula34"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x6.png"  xlink:type="simple"/></disp-formula><p>where κ is the Von Karman constant, U is the wind speed, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x7.png" xlink:type="simple"/></inline-formula>is the friction velocity, z<sub>0</sub> is the roughness length and z is the reference height. Considering the surface shear stress:</p><disp-formula id="scirp.74028-formula35"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x8.png"  xlink:type="simple"/></disp-formula><p>the friction coefficient is:</p><disp-formula id="scirp.74028-formula36"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x9.png"  xlink:type="simple"/></disp-formula><p>Currently, one of the most common hub height of wind turbines is approximately 80 meters, which is much higher than the blending height [<xref ref-type="bibr" rid="scirp.74028-ref30">30</xref>] . Thus, the reference height (z) for the evaluation of the friction coefficient has been fixed at 80 meters.</p><p>The roughness in the vicinity, but outside of a sector, may influence the fluid inside the sector, particularly near the wind turbine of interest. Thus, the reference roughness cannot only be defined from the roughness values within the sector delimitation. The considered area is taken as a rectangle with a certain width, as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>, which shows the area used to calculate the roughness influence in a specific sector.</p><p>The equation of the reference friction coefficient thus becomes:</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Area used to evaluate the roughness influence in a specific sector</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-2760114x10.png"/></fig><disp-formula id="scirp.74028-formula37"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x11.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x12.png" xlink:type="simple"/></inline-formula> is the surface length at the surface n, R is the radius of the sector, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x13.png" xlink:type="simple"/></inline-formula>is the friction coefficient at the surface n and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x14.png" xlink:type="simple"/></inline-formula> is the reference friction coefficient. <xref ref-type="fig" rid="fig2">Figure 2</xref> shows the friction coefficients over the considered area.</p><p>From Equation (4), the reference roughness length is:</p><disp-formula id="scirp.74028-formula38"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x15.png"  xlink:type="simple"/></disp-formula><p>Then, the roughness length is converted to roughness classes according to the European Wind Atlas [<xref ref-type="bibr" rid="scirp.74028-ref26">26</xref>] :</p><p>For<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x16.png" xlink:type="simple"/></inline-formula>:</p><disp-formula id="scirp.74028-formula39"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x17.png"  xlink:type="simple"/></disp-formula><p>For<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x18.png" xlink:type="simple"/></inline-formula>:</p><disp-formula id="scirp.74028-formula40"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x19.png"  xlink:type="simple"/></disp-formula><p>The roughness variations have been quantified as the sum of the differences of the friction coefficients on the site:</p><disp-formula id="scirp.74028-formula41"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x20.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x21.png" xlink:type="simple"/></inline-formula> is the roughness complexity index and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x22.png" xlink:type="simple"/></inline-formula> is defined as:</p><disp-formula id="scirp.74028-formula42"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x23.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x24.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x25.png" xlink:type="simple"/></inline-formula> are respectively the drag coefficients before and after the roughness change n.</p><p>Subsequently, the roughness complexity index <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x26.png" xlink:type="simple"/></inline-formula> has been divided in classes. The development of these classes has been done as follows: First, the roughness complexity index has been calculated for each ten degree sector of the wind turbines studied. Then, the 8028 roughness complexity indices have been analyzed statistically. <xref ref-type="table" rid="table4">Table 4</xref> shows the resulting statistical parameters. Assuming</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Friction coefficients over the considered area</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-2760114x27.png"/></fig><p>that the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x28.png" xlink:type="simple"/></inline-formula> indices analyzed represent the range of possible <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x29.png" xlink:type="simple"/></inline-formula> in a wind farm area, the classes have been selected by quartiles (see <xref ref-type="table" rid="table5">Table 5</xref>).</p></sec><sec id="s5"><title>5. Correlation Parameters of the Proposed Indices</title><p>To be able to evaluate the effect of the different proposed indices on the AEP accuracy, a multiple linear correlation, using the least square method, between the AEP error and the indices has been done.</p><p>The indices have been transformed mathematically to obtain a relationship between terrain characteristics at the met mast position and at the wind turbine position. <xref ref-type="table" rid="table6">Table 6</xref> presents these parameters. On this basis, the correlation equation proposed in this work is:</p><disp-formula id="scirp.74028-formula43"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x30.png"  xlink:type="simple"/></disp-formula><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Statistical values of the C<sub>R</sub> analyzed</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Statistical parameter</th><th align="center" valign="middle" >Value</th></tr></thead><tr><td align="center" valign="middle" >Number of C<sub>R</sub><sub> </sub></td><td align="center" valign="middle" >8028</td></tr><tr><td align="center" valign="middle" >Average</td><td align="center" valign="middle" >0.0154</td></tr><tr><td align="center" valign="middle" >Standard deviation</td><td align="center" valign="middle" >0.0148</td></tr><tr><td align="center" valign="middle" >Minimum</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >Q1</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >Median (Q2)</td><td align="center" valign="middle" >0.0124</td></tr><tr><td align="center" valign="middle" >Q3</td><td align="center" valign="middle" >0.0258</td></tr><tr><td align="center" valign="middle" >Maximum</td><td align="center" valign="middle" >0.0818</td></tr></tbody></table></table-wrap><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> Classes for the roughness complexity index C<sub>R</sub></title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="2"  ></th><th align="center" valign="middle" >C<sub>R</sub> Values</th></tr></thead><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >From</td><td align="center" valign="middle" >To</td></tr><tr><td align="center" valign="middle" >Class 1</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >Class 2</td><td align="center" valign="middle" >&gt;0</td><td align="center" valign="middle" >0.0124</td></tr><tr><td align="center" valign="middle" >Class 3</td><td align="center" valign="middle" >&gt;0.0124</td><td align="center" valign="middle" >0.0258</td></tr><tr><td align="center" valign="middle" >Class 4</td><td align="center" valign="middle" >&gt;0.0258</td><td align="center" valign="middle" >0.0818</td></tr></tbody></table></table-wrap><table-wrap id="table6" ><label><xref ref-type="table" rid="table6">Table 6</xref></label><caption><title> Correlation parameters</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Index</th><th align="center" valign="middle" >Variable</th><th align="center" valign="middle" >Correlation parameters</th></tr></thead><tr><td align="center" valign="middle" >Drop</td><td align="center" valign="middle" >RIX<sub>10</sub></td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x31.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" >Flow separation</td><td align="center" valign="middle" >RIX<sub>40</sub></td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x32.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" >Reference roughness</td><td align="center" valign="middle" >z<sub>0ref</sub></td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x33.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" >Roughness complexity</td><td align="center" valign="middle" >C<sub>R</sub></td><td align="center" valign="middle" ><sub><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x34.png" xlink:type="simple"/></inline-formula> </sub></td></tr><tr><td align="center" valign="middle" >Distance between the turbine and the met mast</td><td align="center" valign="middle" >l</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x35.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" >Height difference between the turbine hub height and met mast top measurement height</td><td align="center" valign="middle" >Δh</td><td align="center" valign="middle" >Δh</td></tr></tbody></table></table-wrap><p>where a, b, c, d, e and f are the variables to be solved.</p><p>In order to avoid using suspicious data, usually resulting from situations which cannot be identified only by looking at the power curve, the wind turbines which have a different trend than most of the other wind turbines of the wind farms have been removed from the correlation. The methodology used to remove these wind turbines is based on the Peirce criterion [<xref ref-type="bibr" rid="scirp.74028-ref31">31</xref>] .</p><p>The Peirce criterion is a simple method to eliminate the outlier data which are not representative of the sample of data set. The suspicious measurements are eliminated when the following condition is fulfilled:</p><disp-formula id="scirp.74028-formula44"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x36.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x37.png" xlink:type="simple"/></inline-formula> is a measured data value, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x38.png" xlink:type="simple"/></inline-formula>is the mean of the data set, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x39.png" xlink:type="simple"/></inline-formula>is the “sample” standard deviation of the complete data set and P is the Pierce criterion. Ross [<xref ref-type="bibr" rid="scirp.74028-ref31">31</xref>] denoted the Pierce criterion as the representation of the ratio of the maximum allowable deviation of a measured value from the data mean to the standard deviation, namely:</p><disp-formula id="scirp.74028-formula45"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x40.png"  xlink:type="simple"/></disp-formula></sec><sec id="s6"><title>6. Results</title><sec id="s6_1"><title>6.1. Average Indices</title><p>The average indices obtained for each wind farm are shown in <xref ref-type="table" rid="table7">Table 7</xref>. The average indices give a good overview of the site complexity. It also shows that the proposed methodology is applicable to different categories of sites.</p><p>As mentioned in Section 2, Wind Farm I is characterized by a very simple terrain with a reference roughness of Class 1, a roughness complexity of Class 2 and a drop index of Class 1. Wind Farms IV and V are characterized by simple terrains with a reference roughness of Class 1, the roughness complexity is relatively high (Class 3) and the drop index is low (Class 2).</p><p>Wind Farms II and III have a relatively low reference roughness (Class 2) and a low roughness complexity (Class 2); however the drop index is high (Class 4)</p><table-wrap id="table7" ><label><xref ref-type="table" rid="table7">Table 7</xref></label><caption><title> Average indices for each wind farm</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Indices</th><th align="center" valign="middle" >Wind Farm I</th><th align="center" valign="middle" >Wind Farm II</th><th align="center" valign="middle" >Wind Farm III</th><th align="center" valign="middle" >Wind Farm IV</th><th align="center" valign="middle" >Wind Farm V</th></tr></thead><tr><td align="center" valign="middle" >Reference roughness</td><td align="center" valign="middle" >0.0027 (Class 1)</td><td align="center" valign="middle" >0.6899 (Class 2)</td><td align="center" valign="middle" >0.6366 (Class 2)</td><td align="center" valign="middle" >0.1474 (Class 1)</td><td align="center" valign="middle" >0.1521 (Class 1)</td></tr><tr><td align="center" valign="middle" >Roughness complexity</td><td align="center" valign="middle" >0.0011 (Class 2)</td><td align="center" valign="middle" >0.0051 (Class 2)</td><td align="center" valign="middle" >0.0043 (Class 2)</td><td align="center" valign="middle" >0.0104 (Class 3)</td><td align="center" valign="middle" >0.0111 (Class 3)</td></tr><tr><td align="center" valign="middle" >Drop</td><td align="center" valign="middle" >0 (Class 1)</td><td align="center" valign="middle" >36 (Class 4)</td><td align="center" valign="middle" >45 (Class 4)</td><td align="center" valign="middle" >3 (Class 2)</td><td align="center" valign="middle" >2 (Class 2)</td></tr><tr><td align="center" valign="middle" >Flow separation</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >Height difference (m)</td><td align="center" valign="middle" >−7</td><td align="center" valign="middle" >−2</td><td align="center" valign="middle" >−3</td><td align="center" valign="middle" >29</td><td align="center" valign="middle" >24</td></tr><tr><td align="center" valign="middle" >Distance (m)</td><td align="center" valign="middle" >1529</td><td align="center" valign="middle" >4608</td><td align="center" valign="middle" >3753</td><td align="center" valign="middle" >3095</td><td align="center" valign="middle" >3686</td></tr></tbody></table></table-wrap><p>and flow separation is observed in these area. Thus, these wind farms are considered as complex terrains.</p></sec><sec id="s6_2"><title>6.2. Annual Energy Production Errors</title><p><xref ref-type="fig" rid="fig3">Figure 3</xref> shows the average AEP errors and the average absolute AEP errors for the wind turbines, along with the global AEP errors for the wind farms. From <xref ref-type="fig" rid="fig3">Figure 3</xref>, it can be seen that for Wind Farm I, WAsP underestimates the AEP while MS-Micro gives better results but nonetheless overestimates the AEP.</p><p>Regarding the wind farms located in complex sites, both models underestimate the AEP of Wind Farm II, with an AEP error being relatively high by exceeding 10%. However, WAsP gives the best results for this wind farm. For its part, Wind Farm III is characterized by the largest AEP errors, where the models overestimate the AEP of Wind Farm III, with WAsP having the smallest AEP errors.</p><p>For the wind farms located in simple sites, the AEP of Wind Farm IV is overestimated by WAsP and MS-Micro; however, the AEP errors are relatively low, which is expected for a simple terrain wind farm. For this wind farm, MS-Micro gives the best results, albeit by a relatively small margin.</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Average AEP errors between the measured values and the estimated values for both numerical models</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-2760114x41.png"/></fig><p>Finally, Wind Farm V has high average absolute AEP errors for the wind turbines, and these AEP errors are in the same order as the AEP errors for the wind farms in complex terrains. This result suggests that errors, or a component of the errors, are induced by characteristics that are not related to the terrain. For this wind farm, MS-Micro gives the best results.</p><p>As expected, these results show that the largest AEP errors occur for wind farms in complex terrains. Furthermore, it appears from these results, and for the wind farms studied, that MS-Micro gives better results than WAsP for wind farms characterized by simple sites, while, in contrast, WAsP gives better results for wind farms characterized by complex sites. This is only an observation since the numerous calculations done were not clearly indicating any physical reason for this behavior.</p></sec><sec id="s6_3"><title>6.3. Example of Correlation Results</title><p>This section shows an example of detailed correlation results. For this purpose, Wind Farm III has been analyzed with MS-Micro; <xref ref-type="table" rid="table8">Table 8</xref> shows the initial AEP errors.</p><p>As expected and as shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>, especially in complex terrains, the distance of the wind turbines from the met mast has a significant influence on the AEP<sub>error</sub>.</p><p>It has been observed that above 8 km, the influence of the distance on the AEP<sub>error</sub> is so high that the correlation with the other parameters, enumerated in</p><table-wrap id="table8" ><label><xref ref-type="table" rid="table8">Table 8</xref></label><caption><title> Initial AEP errors―wind farm III―MS-Micro model</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >Turbines available</th><th align="center" valign="middle" >Turbines ≤ 8 km</th></tr></thead><tr><td align="center" valign="middle" >Number</td><td align="center" valign="middle" >30</td><td align="center" valign="middle" >26</td></tr><tr><td align="center" valign="middle" >AEP error―average by turbines</td><td align="center" valign="middle" >18.5%</td><td align="center" valign="middle" >15.3%</td></tr><tr><td align="center" valign="middle" >AEP error―average absolute by turbine</td><td align="center" valign="middle" >18.5%</td><td align="center" valign="middle" >15.3%</td></tr><tr><td align="center" valign="middle" >AEP error―wind farm</td><td align="center" valign="middle" >17.1%</td><td align="center" valign="middle" >14.4%</td></tr></tbody></table></table-wrap><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Influence of the distance on the AEP<sub>error</sub>―wind farm III―MS-Micro model</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-2760114x42.png"/></fig><p><xref ref-type="table" rid="table6">Table 6</xref>, is not representative. According to this observation, and in order to evaluate the influence of all the parameters, the wind turbines located at a distance above 8 km were removed from the analyses.</p><p><xref ref-type="fig" rid="fig5">Figure 5</xref> presents the correlation curve obtained for Wind Farm III using the parameters of Equation (11). The values of the parameters are given in Equation (14).</p><disp-formula id="scirp.74028-formula46"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x43.png"  xlink:type="simple"/></disp-formula><p>This correlation has an R<sup>2</sup> of 96.8% and a mean absolute error of 2.2%. The AEP errors corrected and the improvement of this methodology are presented in <xref ref-type="table" rid="table9">Table 9</xref>.</p><p>Equation (15) shows the formulation for the AEP errors corrected, namely:</p><disp-formula id="scirp.74028-formula47"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x44.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x45.png" xlink:type="simple"/></inline-formula> is the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x46.png" xlink:type="simple"/></inline-formula> calculated with the correlation equation.</p><p>Then, Equation (16) shows the improvement calculation:</p><disp-formula id="scirp.74028-formula48"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-2760114x47.png"  xlink:type="simple"/></disp-formula></sec><sec id="s6_4"><title>6.4. Overall Results</title><p><xref ref-type="fig" rid="fig6">Figure 6</xref> and <xref ref-type="fig" rid="fig7">Figure 7</xref> show, for each wind farm, the AEP error corrected and</p><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Correlation results―wind farm III―MS-Micro model</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-2760114x48.png"/></fig><table-wrap id="table9" ><label><xref ref-type="table" rid="table9">Table 9</xref></label><caption><title> Correction results―wind farm III―MS-Micro model</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >AEP errors</th><th align="center" valign="middle" >AEP errors corrected</th><th align="center" valign="middle" >Improvement</th></tr></thead><tr><td align="center" valign="middle" >Average by turbines</td><td align="center" valign="middle" >15.3%</td><td align="center" valign="middle" >−0.8%</td><td align="center" valign="middle" >94.9%</td></tr><tr><td align="center" valign="middle" >Average absolute by turbine</td><td align="center" valign="middle" >15.3%</td><td align="center" valign="middle" >2.7%</td><td align="center" valign="middle" >82.2%</td></tr><tr><td align="center" valign="middle" >Wind farm</td><td align="center" valign="middle" >14.4%</td><td align="center" valign="middle" >−1.9%</td><td align="center" valign="middle" >87.1%</td></tr></tbody></table></table-wrap><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> AEP errors corrected for each wind farm</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-2760114x49.png"/></fig><p>the improvement obtained in predicting the AEP using the proposed methodology, respectively.</p><p>The AEP errors for the wind farms, initially up to 17%, were reduced to values below 5% when using the methodology based on the proposed indices. Similar observations are applicable for the AEP errors and the absolute average AEP errors for the wind turbines, where the corrected errors are respectively below 2% and 5%, with the exception of Wind Farm V.</p><p>It has been seen that Wind Farm V has relatively high AEP errors for a simple terrain wind farm. This observation shows that the AEP errors could be influenced by other parameters that this study does not take into account.</p><p>Nevertheless, the methodology proposed shows a significant improvement of the AEP errors, up to 98% for Wind Farm III. The average improvement is 83% for the average AEP errors for the wind turbines and 74% for the average AEP error of the wind farms.</p></sec></sec><sec id="s7"><title>7. Conclusion and Future Work</title><p>A methodology has been developed to reduce the errors in the prediction of the</p><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> AEP error improvement for each wind farm</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-2760114x50.png"/></fig><p>annual energy production (AEP) of wind farm power production. Various indices have been proposed to better identify the characteristics of wind farm sites. Notably, a roughness classification has been developed, which includes a reference roughness index and a roughness complexity index.</p><p>Moreover, a drop index, a flow separation index, the horizontal and vertical distance between the met mast and the wind turbines have also been proposed and used. These indices have been transformed mathematically to obtain relationships between terrain characteristics at the met mast locations and the wind turbine locations. Then, multiple linear correlations have been performed in order to develop equations of the predicted AEP error as a function of these indices.</p><p>The efficiency of this methodology has been proven by reducing the AEP errors by an average of 83% for the wind turbines and by an average of 74% for the wind farms.</p><p>The methodology could also be applied to the wind speed instead of the energy production. If the methodology is validated with the wind speed as the reference variable, it would be possible to use this methodology in the wind condition assessment studies before the implantation of wind turbines. Thus, the errors on the wind condition at the wind turbine locations and on the energy production assessment could be considerably reduced.</p><p>Some additional analysis could also be done to improve the classifications; notably, the variation of the radius of analysis, defined initially by 20 D, could be optimized.</p><p>Moreover, the terrain and the roughness indices could be weighted according to their distance from the turbine, while other indices could be included in the correlation (e.g. flow stability, turbulence, etc.).</p></sec><sec id="s8"><title>Acknowledgements</title><p>This work has been done within the Wind Energy Strategic Network (WESNet), a Canadian research network supported by the Natural Sciences and Engineering Research Council (NSERC) of Canada and by industry. The authors are also very grateful to the PEI Energy Corporation and an undisclosed Canadian wind energy independent power producer for providing production data of their wind farms.</p></sec><sec id="s9"><title>Cite this paper</title><p>Dorval, J., Masson, C. and Gagnon, Y. (2017) On the Improvement of Wind Power Predictions Based on Terrain Characteristics and Measurements of the Annual Energy Production. Journal of Flow Control, Measurement &amp; Visualization, 5, 1-20. https://doi.org/10.4236/jfcmv.2017.51001</p></sec><sec id="s10"><title>Nomenclature</title><p>a Correlation coefficient (see Equation (11))</p><p>b Correlation coefficient (see Equation (11))</p><p>c Correlation coefficient (see Equation (11))</p><p>d Correlation coefficient (see Equation (11))</p><p>e Correlation coefficient (see Equation (11))</p><p>f Correlation coefficient (see Equation (11))</p><p>κ von Karman constant</p><p>l Distance between the turbine and the mast [m]</p><p>σ Sample standard deviation</p><p>ρ Air density [kg/m<sup>3</sup>]</p><p>τ<sub>w</sub><sub> </sub>Surface shear stress [N]</p><p>u<sub>* </sub><sub> </sub>Friction velocity [m/s]</p><p>x Data value</p><p>z Elevation with respect to the ground [m]</p><p>z<sub>o</sub> Roughness length [m]</p><p>AEP Annual energy production [kWh]</p><p>C<sub>d </sub>Friction coefficient</p><p>Class Roughness class</p><p>C<sub>R</sub><sub> </sub>Roughness complexity</p><p>Improvement Improvement calculation (see Equation (16))</p><p>P Pierce criterion</p><p>R Radius of the sector [m]</p><p>RIX Ruggedness parameter<sub> </sub></p><p>U Wind speed [m/s]</p><p>ΔC<sub>d</sub> Variation in friction coefficient</p><p>Δh Difference in turbine and mast heights [m]</p><p>Δr Surface length [m]</p></sec><sec id="s11"><title>Subscript</title><p>corrected related to corrected values</p><p>correlated related to the correlated values</p><p>error related to the relative error</p><p>estimated related to the estimation values</p><p>i related to a given value</p><p>measured related to the measurement values</p><p>n related to the n surface</p><p>ref related to reference value</p><p>Mast related to the mast position</p><p>Turbine related to the turbine position</p><p>1 related to the upstream roughness change</p><p>2 related to the downstream roughness change</p><p>10 computed values using 10% slope</p><p>40 computed values using 40% slope</p></sec><sec id="s12"><title>Superscript</title><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2760114x51.png" xlink:type="simple"/></inline-formula> averaging operator</p><disp-formula id="scirp.74028-formula49"><graphic  xlink:href="http://html.scirp.org/file/1-2760114x52.png"  xlink:type="simple"/></disp-formula><p>Submit or recommend next manuscript to SCIRP and we will provide best service for you:</p><p>Accepting pre-submission inquiries through Email, Facebook, LinkedIn, Twitter, etc.</p><p>A wide selection of journals (inclusive of 9 subjects, more than 200 journals)</p><p>Providing 24-hour high-quality service</p><p>User-friendly online submission system</p><p>Fair and swift peer-review system</p><p>Efficient typesetting and proofreading procedure</p><p>Display of the result of downloads and visits, as well as the number of cited articles</p><p>Maximum dissemination of your research work</p><p>Submit your manuscript at: http://papersubmission.scirp.org/</p><p>Or contact jfcmv@scirp.org</p></sec></body><back><ref-list><title>References</title><ref id="scirp.74028-ref1"><label>1</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Ross</surname><given-names> S.M. </given-names></name>,<etal>et al</etal>. (<year>2003</year>)<article-title>Peirce’s Criterion for the Elimination of Suspect Experimental Data</article-title><source> Journal of Engineering Technology</source><volume> 20</volume>,<fpage> 38</fpage>-<lpage>41</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.74028-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Mason, P.J. (1988) The Formation of a Really-Averaged Roughness Lengths. Quarterly Journal of the Royal Meteorological Society, 114, 399-420.  
https://doi.org/10.1002/qj.49711448007</mixed-citation></ref><ref id="scirp.74028-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Mason, P.J. and King, J.C. (1984) Atmospheric Flow over a Succession of Nearly Two-Dimensional Ridges and Valleys. Quarterly Journal of the Royal Meteorological Society, 110, 821-845. https://doi.org/10.1002/qj.49711046604</mixed-citation></ref><ref id="scirp.74028-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Taylor, P.A., Gent, P.R. and Keen, J.M. (1976) Some Numerical Solutions for Turbulent Boundary Layer Flow above Fixed Rough, Wavy Surfaces. Geophysical Journal International, 44, 177-201. https://doi.org/10.1111/j.1365-246X.1976.tb00281.x</mixed-citation></ref><ref id="scirp.74028-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Newley, T.J. (1985) Turbulent Airflow over Hills. Cambridge University, Cambridge.</mixed-citation></ref><ref id="scirp.74028-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Troen, I. and Petersen, E.L. (1989) European Wind Atlas. Ris&amp;oslash; National Laboratory, Roskilde.</mixed-citation></ref><ref id="scirp.74028-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Grant, A.L.M. and Mason, P.J. (1990) Observations of Boundary-Layer Structure over Complex Terrain. Quarterly Journal of the Royal Meteorological Society, 116, 159-186. https://doi.org/10.1002/qj.49711649107</mixed-citation></ref><ref id="scirp.74028-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Wood, N. (1995) The Onset of Separation in Neutral, Turbulent Flow over Hills. Boundary-Layer Meteorology, 76, 137-164. https://doi.org/10.1007/BF00710894</mixed-citation></ref><ref id="scirp.74028-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Mortensen, N.G., Heathfield, D.N., Myllerup, L., Landberg, L. and Rathmann, O. (2007) Getting Started with WAsP 9. Riso National Laboratory, Roskilde.</mixed-citation></ref><ref id="scirp.74028-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">IEC (2013) Power Performance of Electricity Producing Wind Turbines Based on Nacelle Anemometry. 61400-12-2, IEC, Geneva.</mixed-citation></ref><ref id="scirp.74028-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">AWS SI (1997) Wind Resource Assessment Handbook. National Renewable Energy Laboratory, Golden.</mixed-citation></ref><ref id="scirp.74028-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">IEC (2005) Power Performance Measurements of Electricity Producing Wind Turbines. 61400-12-1, IEC, Geneva.</mixed-citation></ref><ref id="scirp.74028-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Jimenez, B., Durante, F., Lange, B., Kreutzer, T. and Tambke, J. (2007) Offshore Wind Resource Assessment with WAsP and MM5: Comparative Study for the German Bight. Wind Energy, 10, 121-134. https://doi.org/10.1002/we.212</mixed-citation></ref><ref id="scirp.74028-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Feng, S.-L., Wang, W.-S., Liu, C. and Dai, H.-Z. (2009) Study on the Error Evaluation Method of Calculated Energy Production by WAsP When the Terrain Gets Complex. International Conference on Sustainable Power Generation and Supply, Nanjing, 6-7 April 2009, 1-4. https://doi.org/10.1109/SUPERGEN.2009.5348148</mixed-citation></ref><ref id="scirp.74028-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Berge, E., Nyhammer, F., Tallhaug, L. and Jacobsen, O. (2006) An Evaluation of the WAsP Model at a Coastal Mountainous Site in Norway. Wind Energy, 9, 131-140.  
https://doi.org/10.1002/we.191</mixed-citation></ref><ref id="scirp.74028-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Mortensen, N.G., Bowen, A.J. and Antoniou, I. (2006) Improving WAsP Predictions in (too) Complex Terrain. Proceedings of 2006 European Wind Energy Conference and Exhibition, Athens, 27 February-2 March 2006, 175-183.</mixed-citation></ref><ref id="scirp.74028-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Rathmann, O., Mortensen, N.G., Landberg, L. and Bowen, A. (1997) Assessing the Accuracy of WAsP in Non-Simple Terrain. Proceedings of BWEA Conference, Exeter, 25-27 September 1996, 413-418.</mixed-citation></ref><ref id="scirp.74028-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Bowen, A.J. and Mortensen, N.G. (2004) WAsP Prediction Errors Due to Site Orography. Ris&amp;oslash; National Laboratory, Roskilde, 65.</mixed-citation></ref><ref id="scirp.74028-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Bowen, A.J. and Mortensen, N.G. (1996) Exploring the Limits of WAsP the Wind Atlas Analysis and Application Program. Proceedings of 1996 European Union Wind Energy Conference, G&amp;ouml;teborg, 20-24 May 1996, 584-587.</mixed-citation></ref><ref id="scirp.74028-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Jackson, P.S. and Hunt, J.C.R. (1975) Turbulent Wind Flow over a Low Hill. Quarterly Journal of the Royal Meteorological Society, 101, 929-955.  
https://doi.org/10.1002/qj.49710143015</mixed-citation></ref><ref id="scirp.74028-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Zephyr North (2010) MS-Micro Version 3.</mixed-citation></ref><ref id="scirp.74028-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Ris&amp;oslash; National Laboratory, WAsP 10 (2010).</mixed-citation></ref><ref id="scirp.74028-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Kwon, S.-D. (2010) Uncertainty Analysis of Wind Energy Potential Assessment. Applied Energy, 87, 856-865. https://doi.org/10.1016/j.apenergy.2009.08.038</mixed-citation></ref><ref id="scirp.74028-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Touani, I., Yu, W., Gagnon, Y. and Masson, C. (2013) Development of a Wind to Power Model for Wind Farm Power Production Forecasting. Wind Engineering, 37, 347-368. https://doi.org/10.1260/0309-524X.37.4.347</mixed-citation></ref><ref id="scirp.74028-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Rehman, S., Aftab, A. and Al-Hadhrami, L.M. (2011) Development and Economic Assessment of a Grid Connected 20 MW Installed Capacity Wind Farm. Renewable and Sustainable Energy Reviews, 15, 833-838.  
https://doi.org/10.1016/j.rser.2010.09.005</mixed-citation></ref><ref id="scirp.74028-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Rehman, S., Mahbub, A.M., Meyer, J.P. and Al-Hadhrami, L.M. (2012) Wind Speed Characteristics and Resource Assessment Using Weibull Parameters. International Journal of Green Energy, 9, 800-814. https://doi.org/10.1080/15435075.2011.641700</mixed-citation></ref><ref id="scirp.74028-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Jung, S., Vanli, O.A. and Kwon, S.D. (2013) Wind Energy Potential Assessment Considering the Uncertainties Due to Limited Data. Applied Energy, 102, 1492-1503. https://doi.org/10.1016/j.apenergy.2012.09.011</mixed-citation></ref><ref id="scirp.74028-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Carvalho, D., Rocha, A., Santos, C.S. and Pereira, R. (2013) Wind Resource Modeling in Complex Terrain using Different Mesoscale-Microscale Coupling Techniques. Applied Energy, 108, 493-504.  
https://doi.org/10.1016/j.apenergy.2013.03.074</mixed-citation></ref><ref id="scirp.74028-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Royer, J. (2011) Getting the Facts on Wind Energy in Canada: Notable Information Derived from the Federal Wind Programs. CanWEA 2011, Vancouver, Canada.</mixed-citation></ref><ref id="scirp.74028-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">ecoERP (2016) ecoEnergy for Renewable Power, Natural Resources Canada.  
www.nrcan.gc.ca/ecoaction/14145</mixed-citation></ref><ref id="scirp.74028-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">CanWEA, Canadian Wind Energy Association (2016) Installed Capacity.  
http://canwea.ca/wind-energy/installed-capacity/</mixed-citation></ref></ref-list></back></article>