VTT TECHNICAL RESEARCH CENTRE OF FINLAND LTD Estimating icing in Finnish climate conditions Ville Lehtomäki, Timo Karlsson, Simo Rissanen, VTT Teppo Hilakivi, Puhuri Oy Staffan Asplund, Etha Wind Wind Finland 2015 seminar 29 th October 2015, Helsinki
Content Presentation 1) by VTT: Comparing different turbine SCADA to Finnish Icing Atlases Definitions & Introduction to icing Three different Icing Atlases for Finland Case Puhuri Results Conclusions Presentation 2) by Etha Wind: ON ICING: Techno-Commercial Aspects 5.11.15 2
Definitions 3
Meteorological and instrumental icing IEA Wind Task 19 international expert group, goal: Gather and provide information about wind energy in cold climates Download free Best practices for wind energy in cold climates Resource assessment and site classification Task 19 website 5.11.15 4
Ice classes: IEA Ice Classification¹ IEA ice class Duration of Meteorological icing [% of year] Duration of Instrumental icing [% of year] Production loss [% of AEP] 5 >10 >20 >20 4 5-10 10-30 10-25 3 3-5 6-15 3-12 2 0.5-3 1-9 0.5-5 1 0-0.5 <1.5 0-0.5 ¹: IEA Wind Recommended Practices for wind energy projects in cold climates edition 2011, Task 19 5.11.15 5
Introduction
Wind farm development process & lifetime costs Site prospecting & site assessment the most critical decision gate in project development for AEP & risk assessment Topic Today Decision gate Prospecting Planning & Resource assessment Turbine purchase Constructions O&M Decomission X1 x6 x50 x20 x20? 5.11.15 7
Icing effects for wind power in Finland Manual de-icing 5.11.15 8
Icing effects for wind power in Finland AEP losses in Finland are expensive for developer because of FIT = 83.5 /MWh starting 2016-> 3MW turbine, -5 % AEP due to icing = 35 k /turbine/year lost Future issues? Ice throw risks need to be assessed Ice hit not an option Results from Sweden: iced turbine noise 5.11.15 9
Icing atlases for Finland 10
Three icing atlases for Finland Ø Used for site prospecting & initial risk assessment Method Approach Validation Advantage FMI Icing atlas Kjeller Icing atlas VTT WiceAtlas v1.0 - Numerical weather modelling with AROME - Icing hours for vertical Ø3cm, 1m tall cylinder (ISO 12494) - threshold 10g/m/h - 5yr average conditions - Publically available Only one 15 h icing event - Long-term experience of weather modelling in FI - Numerical weather modelling with WRF - Icing hours for vertical Ø3cm, 1m tall cylinder (ISO 12494) - threshold 10g/m/h - Detail simulation 7/2008-7/2009, 20yrs correction - Publically available Multiple validation points in Sweden & Norway - Internationally known & world-class expertize in icing modelling - Cloud and temperature measurements - Point measurements from met stations N=44 - Icing hours when T<0 C and low level clouds present - No threshold - 20yr mean conditions - Interpolate between points -> color map - Non-public available Multiple validation points in Sweden, Canada - No simulations, only measurements 5.11.15 11
Three icing atlases for Finland Icing hours [% of year] FMI Icing atlas 100m agl FI Ave: 1.1 % 5.11.15 Kjeller Icing atlas 140m agl FI Ave: 1.6 % VTT WiceAtlas v1.0 150m agl FI Ave: 6.7 % 12
Case Puhuri 13
Case Puhuri - overview Puhuri Oy is a utility owned developer with portfolio of 4 wind farms in Finland Several projects in pipeline Case wind farm: Turbines A & B (3MW, HH140m, D120m) A & B close to each other 5.11.15 14
Case Puhuri datasets used for analysis SCADA from two turbines: A & B Instrumental Icing (heated/unheated anemometers) on top of turbine A nacelle Meteorological Icing (Labkotec LID) on top of turbine A nacelle 3 x icing atlases: FMI, Kjeller, VTT WiceAtlas 5.11.15 15
Methods used Use IEA Wind Task 19 open source software, download T19IceLossMethod here Input Standard SCADA data Ws, temp, pwr, nacdir, mode Software Calc non-iced, ref power Calc icing events Calc AEP losses due to icing Output Icing timeseries AEP losses Turbine A heated/unheated anemometry as ice detector: if unheated < 20 % of heated, then Instrumental Icing Turbine A meteorological ice detector (Labkotec LID): if warning, then Meteorological Icing 5.11.15 16
Site prospecting analysis: 2, 5 or 10 % losses? FMI met ice Kjeller met ice % of year Wice met ice FMI inst ice FMI AEP loss IEA ice class Duration of Meteorological icing [% of year] Duration of Instrumental icing [% of year] Production loss [% of AEP] Turbine A 3.7 2.0 7.6 30.9 5.70 Turbine B 4.0 2.1 7.6 32.0 5.90 Turbine A Ave Production losses due to icing [% of AEP] (3+12) /2 =7.5 2.8 17.5 >20 Turbine B 7.5 2.8 17.5 >20 Ice meas needed? Yes Yes/No Yes Yes Met icing ave p-loss = 9.3 % AEP Inst icing ave p-loss = > 20 % AEP AEP loss (FMI) = 5.8 % AEP 5 >10 >20 >20 4 5-10 10-30 10-25 3 3-5 6-15 3-12 2 0.5-3 1-9 0.5-5 1 0-0.5 <1.5 0-0.5 Best estimate: (5.8+9.3)/2 = 7.6 % of AEP Decision: perform icing measurements at site! 5.11.15 17
Results 18
Power curve, turbine A 2015 5.11.15 19
Power curve, turbine B 2015 5.11.15 20
Measurement results regarding icing % of year 20 18 16 14 12 10 8 6 4 2 0 Sublimation! 9-12/2013 2014 1-10/2015 A losses B losses Rotor ice Met ice Inst Ice Good correlation to measured icing (met & inst) to AEP losses More icing with warmer winter temperatures? 5.11.15 21
Monthly R^2 correlations Meteorological icing: Wice vs Labko has best correlation Instrumental icing: FMI 5-yr average values Labko Wice FMI Inst meas Labko 1.00 0.79 0.63 Wice 0.79 1.00 0.70 FMI 0.63 0.70 1.00 Inst FMI Inst meas 1.00 0.50 Inst FMI 0.50 1.00 5.11.15 22
Monthly meteorological icing 0,3 May 2013 - May 2015 R^2=0.79 0,3 % of time 0,2 0,2 0,1 0,1 0,0 met ice Wice met ice Wice overpredicts but correlation good -> calibration needed 5.11.15 23
Monthly instrumental icing % of time 0,7 0,6 0,5 0,4 0,3 0,2 0,1 0,0 May 2013 - May 2015 R^2=0.50 FMI large values all winter inst ice FMI inst ice 5.11.15 24
Long-term icing from WiceAtlas timeseries mean Last two winters on same level as long-term average 5.11.15 25
IEA Ice Classification of site Icing atlases Source Kjeller 2 Wice 4 FMI (meteorological) 3 FMI (instrumental) 4 FMI (losses) 3 AVE= 3 Ice class Icing atlases results > measurements Each atlas has pros/cons Instrumental icing atlas an answer? Source Case Met icing (Labko) & Rotor icing Instrumental icing 2-3 A losses 2-3 2 B losses 2 AVE= 2 Ice class Prospecting estimate Mean meas losses - 7.6 % (atlases) -1-5.6 % Ice Class 3 Ice Class 2 5.11.15 26
Conclusions (1/2) Large scatter of results for one site from 2015: Different meteorological icing atlases 2-8 % of year Turbine A Icing 6 % of year Instrumental icing 7-19 % AEP losses due to icing 1-9 % Are long-term effects same? Site prospecting study: Met icing atlases + IEA ice class Good idea to use more than one icing atlas Kjeller ice class close to case study results FMI and VTT need tuning? Long-term WiceAtlas analysis shows similar icing also in future 5.11.15 27
Conclusions (2/2) Different turbines react diffently to same icing conditions Turbine A = -5.6 % mean AEP losses due to icing Turbine B = -1.3 % mean AEP losses due to icing Prospected losses (-7.6 % AEP) were higher than measured (-1-5.6 % AEP) when using atlases with IEA ice classes Next: Etha Wind continues 5.11.15 28
TECHNOLOGY FOR BUSINESS Ville Lehtomäki ville.lehtomaki@vtt.fi +358 50 370 7669
Extra slides
Low temperature & Icing map of world 150m IEA Ice class AEP loss [%] 1 0-0.5 2 0.5-5 3 3-12 4 10-25 5 > 20 5.11.15 31
Cold climate markets Cumula&ve installed capacity by end of 2012 [MW] Forecasted capacity 2013-17 [MW] Low temperature Light icing: safety risk, some economic risk Moderate to heavy icing: economic and safety risk Low temperature Light icing: safety risk, some economic risk Moderate to heavy icing: economic and safety risk 18,945 41,079 11,478 20,025 22,083 8,003 Total 69,000 (*) Total 45,000 50,000 (*) The total capacity is less than the sum of individual capacities because some of the sites have both low temperatures and icing conditions. 30GW of new installagons to icing condigons by 2017 Ø Compare: new offshore 29GW by 2017! source: BTM World market Update 2012 5.11.15 32
Measurement results A lot of short icing events Ave temps rising Global warming? 5.11.15 33
Icing event example Turbine A stops for one week, turbine B only lightly affected 5.11.15 34