Validierungsbeispiele und Modellprobleme: Difference between revisions

From Wiki-WindPRO
Jump to navigation Jump to search
No edit summary
No edit summary
 
(4 intermediate revisions by the same user not shown)
Line 1: Line 1:
#REDIRECT [[:Category:Validierung]]
{{Hauptlink}}
{{Hauptlink}}
===Meso data long term consistency===
[[File:DE_Energy_Validation(1).png|800px]]
'''Image:Long-term consistency using Meso data'''
This screen shot shows how well the meso data based calculation matches measured production on a time scale of 18 years. The Tunø Knob offshore wind farm has been operating very well with few problems and a high availability during all 18 years. Thereby, it is a good validation case. It is seen that year-by-year model errors are within +/-5% and month by month within +/- 10%, apart from few months where the problems probably were related to availability issues.
===Deep array wake loss model===
====Horns rev 1, Danish Offshore project====
A more comprehensive test of the wake models, especially with deep array settings, has been performed on the Horns Rev 1 wind farm, where EMD has access to wind measurements as well as production data at a 10-minute level for one year. A calculation was set up and data loaded in the Performance Check module for detailed filtering and evaluations of the measurements against calculations on different aggregation levels. Note the deep array settings are only available for the time-step based PARK calculation, for which only the NO2005 model is available.
[[File:DE_Energy_Validation(2).png|300px]][[File:DE_Energy_Validation(3).png|220px]]
'''Image: Wind farm layout and the direction distribution of energy'''
Note the turbine orientations give the densest spacing in north, east, south and west sectors. The results are evaluated in 30-degree sectors, but calculated in 1-degree steps. The first is shown for the 360 degree calculation, meaning a full year with all data.
For the total wind farm, with 360 degrees and all wind speeds, the wake loss is calculated as 10.8% with the original N.O.Jensen model, exactly as measured, but with smaller deviations by turbine row. The "measured" wake losses are found by calculating with 1 year of 10-minute data and finding the ratios of measured/calculated as averages for all non-error time stamps by use of the windPRO Performance Check module. Based on the prediction error and the calculated wake loss, the measured wake loss can then be established for the sample.
[[File:DE_Energy_Validation(4).png|450px]]
'''Image:NO2005 variant and original N.O. Jensen model'''
Inclusion of mirror turbine wake (light blue line), which is new in 3.1 for NO2005 model, has almost no impact for this wind farm size. Orange line is without mirror wake.
The NO2005 variant is seen to calculate less wake loss than the original N.O.Jensen model, and too low compared to measurements.
There is a possibility to improve the NO2005 variant by using the deep array settings when calculating by time step:
[[File:DE_Energy_Validation(5).png|550px]]
'''Image:The deep array loss settings for time step calculations'''
The linear weight of 35% is recommended. The additional option “Use downwind change of WDC by number of wake turbines are only for going down at a very detailed level. This is now in a version 2, that work somewhat differently and require new parameters. This can be used for “ultimate tuning” when very good data are available to calibrate against. See later validation tests.
[[File:DE_Energy_Validation(6).png|450px]]
'''Image:Using 35% linear weight in deep array model setting solve'''
Looking at a more detailed analysis by wind direction sector, two deep array settings are tested and compared to a baseline:
:0. NO2005 standard, WDC 0.04 - the baseline.
:1. Linear weight 35%, WDC 0.04.
:2. Linear weight 100%, but with increased WDC from 0.04 to 0.07.
Setting 2 matches measurements almost perfect at 360 degrees, but looking at individual directions (see below), it is seen that this might not be a robust solution for all projects. For the directions where turbines are most densely spaced (north, east, south and west), an under-prediction is observed for the most wake-affected turbines. Moreover, it can be tricky to “guess” which WDC shall be used.  Setting number 1 handles all directions somewhat better (i.e. there is no systematic bias), and the default WDC is used. This is likely to be a better solution when the wind farm performance is not already known.
The graphs below shows direction by direction in 30 degree sectors and illustrates how the measured/calculated ratio comes out across the average row by row, where all rows are arranged so the upwind row is row 1. The goal for a good wake model is horizontal lines in such graphs – the absolute values are less important, as they can reflect a wind speed bias.
[[File:DE_Energy_Validation(7).png|650px]]
'''Image: Setting 0: Calculation with default NO2005 model, WDC 0.04'''
Row 1, leftmost, is upwind row. Wake losses are under estimated the more downwind in the row.
[[File:DE_Energy_Validation(8).png|650px]]
'''Image: Setting 1: Calculation with 35% linear weight'''
Reasonably straight lines for all direction sectors, meaning the wake model settings work well. (almost no data from NNE and ENE may explain poor behaviour here).
[[File:DE_Energy_Validation(9).png|650px]]
'''Image: Setting 2: Calculation with 100% Linear weight and WDC 0.07'''
Wake losses are over-predicted in back rows for the direction sectors where turbines are oriented along the sector.
As observed, the more upwind the turbines (higher row numbers), the more energy under-prediction (wake loss over prediction) is seen for the east, south and west directions with 100% linear weight – this is somewhat critical. There are, due to the nature of the data, no perfect matches for any row. The north sector is not included due to very few data points in this sector. '''The presented case study can be concluded with no doubt that setting 1; the NO2005 with 35% linear weight is the best alternative.'''
The original N.O. Jensen model is not tested by direction, as the time-step based calculation concept cannot use this model and the very detailed analyses are more difficult to perform.
In addition, and partly for evaluation of the mesoscale time consistence, a longer period is calculated, WITHOUT filtering poor performance problems out. The results can be seen in the figure below:
[[File:DE_Energy_Validation(10).png|600px]]
'''Image:Monthly measured and calculated'''
The graph above shows partly the measured results (including availability loss and internal cabling grid loss) and partly the calculated results based on the standard WDC of 0.04. The deep array loss model alone and the deep array loss model combined with time step turbulence from mesoscale data (the last one is only from 2013, since meso data does not include turbulence before 2013). As can be seen, the difference in the three calculation models are hardly visible in the graph, but there are differences, which can be seen in the graph below:
[[File:DE_Energy_Validation(11).png|400px]]
'''Image: Summarised result with different corrections'''
Summarizing the results for the different calculation setups since March 15, 2013. While the losses are assumed to be in the order of 5% (2% grid and 3% availability loss), it is assumed the most refined calculation setup provides the best result. It is noted that use of the deep array loss model calculates 0.8% higher loss and the WDC changed by time step adds a further 0.7% loss. In combination, that is an addition of 1.5% loss. This seems very little, but, with 10% wake loss with a standard calculation, it is an increase in calculated wake loss of 15%.
[[File:DE_Energy_Validation(12).png|600px]]
'''Image: Monthly ratios measured/calculated'''
The results shown in the above graph are shown as ratios. The reason some months have higher measured values than calculated values is because mesoscale data is not that accurate. Note that only 6 months out of 125 exhibits higher measured values than calculated. For the months with lower measured values than calculated values, it is most likely due to availability problems. Here is seen how the added wake model corrections give slightly lower calculation results (higher ratios).
The conclusion of the test with HR-1 calculations is due partly to the Mesoscale wind data performing “very precisely” (no post scaling) and partly to the deep array model modifies well. That the WDC is decided by TI adds further improvement. All in all, we recommend including all corrections. However, we must set a “warning flag” regarding the use of the meso data turbulence. Usage of the meso data in this manner has proven to provide correct results at this and a few other tested sites, but no comprehensive validation on many sites has been performed so far.
====Lillgrund, Sweden offshore project====
This project is special due to the dense spacing - around 3.2x the rotor diameter (RD). The main wind direction is from WSW, along the row orientation.
[[File:DE_Energy_Validation(14).png|700px]]
'''Image: Measured and calculated from Performance Check module'''
Shown above is the output from the PERFORMANCE CHECK module, where both measured and calculated turbine by turbine can be seen. Data is filtered by taking out larger down times in both measured and calculated values. Monthly production data for each turbine for 5 years (from December 2008) is used. Calculation is based on Mesoscale data. In general, it is a good match, but zooming in based on the ratios shows deviations. This is especially evident when based on the “standard” wake settings. Using the new deep array wake loss options help.
[[File:DE_Energy_Validation(15).png|700px]]
'''Image:Wind farm layout and ratios measured/calculated, Lillgrund offshore'''
As seen by Horns Rev, the new deep wake default settings improve the results quite a bit. The back row calculation performs, essentially, better. It is also observed how the 25% linear weight, by itself, does most of the modification. However, it overestimates the wake losses for the inner turbines, and it still underestimates the back row wake loss to a certain degree. These two problems are both improved by including the WDC reduction by number of up wind turbines. Note, although, that this option does not change results drastically - it just fine-tunes a bit. But, it is a good change, overall, in that larger windfarms with more rows will have a greater improvement in the wake decay results than previously.
===Test of model vs wind farm size===
[[File:DE_Energy_Validation(16).png|550px]]
'''Image: The development of calculated wake loss by wind farm size'''
Using a logarithmic x-axis, the calculated wake loss increases almost linearly with wind farm size when spacing is kept constant.
For small wind farms (< 20 turbines), the three variants calculate almost identical results. The deviations increase as project size increases.
By the time the project reaches 100 turbines, the NO2005 calculates around 2% more AEP than the original N.O. Jensen model. At 250 turbines, this increases to 3%.
Using a linear weight of 35% in the NO2005 combination model, brings the result closer to the original N.O. Jensen model and AEP deviations are less than +/- 1%. Mirror wake is used in NO2005 in the figure above.
===Test of reduced WDC by upwind turbines on Horns Rev-1===
[[File:DE_Energy_Validation(17).png|550px]]
'''Image: Wake loss test. Blue graphs show the two WDC reduction methods'''
The presented graph shows the variation in calculated wake loss by row for several methods.
Green line represents the observed wake loss, the red line represents the wake loss calculated with the original N.O. Jensen model. '''The two blue lines represent the wake calculated based on time varying data with NO 2005 method with 35% linear weight and WDC reduction by using method 1 and 2'''. Finally, the orange line represents the wake calculated with NO2005 with 35% linear weight in combination model and mirroring effect (for more details see [http://help.emd.dk/knowledgebase/content/TechNotes/TechNote_5_Park%20model%20revision.pdf Park model revision paper]).
Although the deviations between all methods are small, it can be seen that the results from the time varying calculation with NO 2005 method and WDC reduction gets closer to the observations than the rest of the methods. Zooming in, a clearer picture is seen:
[[File:DE_Energy_Validation(18).png|550px]]
'''Image: Ratios measured/calculated wake loss illustrates the differences more precise'''
Part of the reason for the lack of full match can be explained by the data. But the trend is clear, the reduced WDC by upwind turbines straighten out the lines by calculating less wake loss in the “inner rows” while the outmost rows are calculated roughly as without WDC reduction. If the wind farm were larger, the reduced WDC would result in higher calculated wake losses compared to not using reduced WDC.
[[File:DE_Energy_Validation(19).png|550px]]
'''Image: Comparison to the Original N.O.Jensen model and NO2005'''
Compared to original N.O.Jensen model it is seen that the NO2005 with 35% Linear weight’ calculates slightly higher wake losses (lower Meas/calc wake loss ratio). Compensating with reduced WDC, slightly lower at the rows which get too high loss calculation by org. NO model. But the most important NO2005 modification is the 35% Linear weight in the combination model. Without this, the NO2005 (grey line), calculates round 15% too low wake loss, which in this example mean about 1.5% too high AEP.
===Test of reduced WDC by upwind turbines on Egypt wind farm===
Here the distances between the rows are essentially larger, while the in row distance is much smaller compared to HR1. Main wind direction is almost uniform, and row 1 is the upwind row. The Figure below presents the wind farm design.
[[File:DE_Energy_Validation(20).png|400px]]
'''Image: Wind farm layout, main wind direction is from NW'''
Approximate RD spacing: in row: 3, between rows: 14.
[[File:DE_Energy_Validation(21).png|550px]]
'''Image: Egypt wind farm with 100 turbines in 7 rows'''
It is seen how the NO2005 model under predicts wake losses remarkable for the back rows (the measured/ calculated AEP is low for back rows). The “base” WDC is here set differently by wind direction based on measured turbulence.
With the applied 35% linear weight problem is almost solved. Including the WDC reduction by upwind turbines improves slightly for rows 2-3, but for the back rows 6-7 it is arguable if an improvement is seen, however other reasons could be behind this.
[[File:DE_Energy_Validation(22).png|550px]]
'''Image: Turbine by turbine measured/calculated energy production'''
Quite some deviations within the rows are seen, which can be partly explained by topographic*) issues, partly by operational issues. It is therefore difficult to judge if the row average values give a correct picture of which model variant is the best. A reason for the small under prediction of the back rows with the WDC reduction model, can here also be the use of too low WDC. Setting this slightly up, would lower the measured/calculated values especially for the back rows.
The good thing is that using the same settings for the WDC reduction as for Horns Rev example, still work well on a fully different layout.
*) It shall be noted that later calculations on the wind farm above show that with more accurate elevation data the most remarkable outliers are brought in line. It can thereby for sure be said that the outlier problems mainly are related to lack of elevation data quality.
===Test of WDC reduction on varying wind farm sizes===
[[File:DE_Energy_Validation(23).png|550px]]
'''Image: WDC reductions by upwind turbines tested on varying wind farm sizes'''
As seen from the graph, the calculated wake loss with NO 2005 35% linear weight model (grey line), as the major model parameter, and the two applied WDC reduction models (yellow and green), shows quite similar values, by looking at the average for the entire wind farm. The benefit from the WDC reduction is that it compensates for the over prediction of the wake losses typically for front rows (2 to aprox.5) and compensates for the under prediction of the wake losses for the downwind rows.
===WAsP versions modifications===
The WAsP model itself will not be explained in detail here (see Risø/DTU WAsP manual), only the changes of high importance for the user in the more recent versions. From the very first versions, until and including ver. 9, the model itself has only changed marginally, and the calculation results, thereby, also. The major improvement during the earlier WAsP model changes is the capability to handle more map file points.  However, from version 10.0, model modifications have been made. These mainly relate to stability correction handling, especially for offshore and near shore, but, also, the roughness map interpretation has been improved. The corrections were partly included in ver. 10.0, but first fully implemented in ver. 10.1 and 10.2. We, therefore, do '''not''' recommend using ver. 10.0 in offshore or coastal regions. The corrections relate to the default heat flux parameters - the way roughness in a coastal zone is interpreted and to formula modifications. The result of the corrections is a smoother change between on shore and offshore stability correction. The windstatistics for offshore is different if it is generated from WAsP 9 or 10.2+.
'''This means that an offshore or near shore windstatistics made from WAsP 9 SHOULD NOT be used in any version from WAsP 10.2 or vice versa.'''
This can best be illustrated by an example: 4 windstatistics are generated from the same time series data:
[[File:DE_Energy_Validation(24).png|275px]][[File:DE_Energy_Validation(25).png|275px]]
[[File:DE_Energy_Validation(26).png|275px]]
[[File:DE_Energy_Validation(27).png|275px]]
'''Image: Four windstatistic results for Lillgrund offshore, different WAsP versions'''
Top row is using WAsP 9 and bottom row is using WAsP 10.2. Left column is assuming a class 0 site and the right column is assuming a class 1 site.
In the left column, the class 0 data is almost the same using the two WAsP versions. The highest level is at 200m, although the result is slightly lower with WAsP 10.2 – a part of the formula modification. BUT, for the onshore classes, the wind speeds are essentially higher with WAsP 10.2 - around 0.2 m/s.
'''This means that a WAsP 10.2 calculation could calculate up to around 10% higher AEP than WAsP 9 at an onshore site, if the data basis were offshore.'''
In the right column, the onshore site, the onshore class data 1, 2 and 3 are almost identical from the two versions. Again 200m has slightly lower results with WAsP 10.2 – a part of the formula modification. BUT the class 0 data is around 0.3 m/s lower when calculated with 10.2.
'''This means that a WAsP 10.2 calculation would calculate up to around 10% lower AEP at an offshore site, if the data basis were onshore.'''
The changes are, of course, not serious as long data is used in the same environment it comes from (e.g., offshore for offshore). But, it should be mentioned that, in the coastal region, where there is part water and part land, quite unpredictable changes can be seen. Here, it will be especially important to use the same WAsP version for generation of windstatistics.
From the tables, it is also seen that, with large hub heights (>100m), WAsP 10.2 calculates slightly lower wind speeds than WAsP 9 – this is so both on and offshore have a general modification of the stability correction model. Based on quality test data in an offshore environment and from tall masts, this correction appears to give a better reproduction of the measurements.
[[File:DE_Energy_Validation(28).png|250px]][[File:DE_Energy_Validation(29).png|370px]]
'''Image: Test case calculations showing the WAsP stability model shift'''
The example above illustrates the improvement regarding smoothing the stability model shift. Two turbines just 200m apart, both around 4 km offshore, are calculated. In sector 5, there is a roughness change from class 0 to class 0.4 at a 10 km distance. The graphs show the ratio of the calculated AEP between WAsP9/WAsP10.2. Looking at the rose, it is seen that, in sector 5, what happens for the two turbines is very different. This is where the distance to the roughness change in sector 5 is just around 10 km. A change from 0 to a higher class decides that WASP 9 shall change between on- and offshore stability in that direction. The parameter: “Width of coastal zone” decides this and can be changed, but the default is 10 km. It’s an obvious inconsistence, which is in WAsP 9, but not in WAsP 10.2. But, the graph also shows how the calculation results are smaller in general with WAsP 9 when the wind come from land, while the 9 and 10.2 results are the same when wind comes from the open sea.
[[File:DE_Energy_Validation(30).png|300px]][[File:DE_Energy_Validation(31).png|300px]]
'''Image: Wind profiles, measured and calculated'''
A test on reproducing measured profiles in large heights is shown above. It is seen that WAsP 10.2 predicts the vertical profile slightly better than WAsP 9. The 160m points (purple squares) are from another mast nearby and, therefore, not fully comparable to the measured (red triangles). Based on both masts measured at 100m, the 160m point is scaled (yellow triangle), and it is observed that the WAsP calculated profile matches measurements well up to 160m, with a small advantage for WAsP 10.2 relative to WAsP 9.
====Tests of calculation in a coastal region with different hub heights====
The calculation setup can be seen in the figure on the next page:
[[Category:PARK (DE)]]

Latest revision as of 14:41, 10 January 2020