Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH

Similar documents
Kernel Methods and SVMs Extension

Statistical Evaluation of WATFLOOD

IRO0140 Advanced space time-frequency signal processing

This column is a continuation of our previous column

Chapter 3 Describing Data Using Numerical Measures

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors

Regulation No. 117 (Tyres rolling noise and wet grip adhesion) Proposal for amendments to ECE/TRANS/WP.29/GRB/2010/3

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

Negative Binomial Regression

[The following data appear in Wooldridge Q2.3.] The table below contains the ACT score and college GPA for eight college students.

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA

Statistics Chapter 4

Chapter 9: Statistical Inference and the Relationship between Two Variables

Cokriging Partial Grades - Application to Block Modeling of Copper Deposits

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6

SIMPLE LINEAR REGRESSION

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation

Statistics MINITAB - Lab 2

A Robust Method for Calculating the Correlation Coefficient

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Chapter 13: Multiple Regression

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

/ n ) are compared. The logic is: if the two

Chapter - 2. Distribution System Power Flow Analysis

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Lecture Notes on Linear Regression

A strong wind prediction model for the railway operation based on the onsite measurement and the numerical simulation

Grover s Algorithm + Quantum Zeno Effect + Vaidman

Lecture 10 Support Vector Machines II

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

Chapter 14 Simple Linear Regression

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

RELIABILITY ASSESSMENT

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

2016 Wiley. Study Session 2: Ethical and Professional Standards Application

Chapter 11: Simple Linear Regression and Correlation

UNR Joint Economics Working Paper Series Working Paper No Further Analysis of the Zipf Law: Does the Rank-Size Rule Really Exist?

Physics 181. Particle Systems

Statistics for Economics & Business

Linear Regression Analysis: Terminology and Notation

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Global Sensitivity. Tuesday 20 th February, 2018

where I = (n x n) diagonal identity matrix with diagonal elements = 1 and off-diagonal elements = 0; and σ 2 e = variance of (Y X).

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

UNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours

Chapter 8 Indicator Variables

Uncertainty and auto-correlation in. Measurement

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Definition. Measures of Dispersion. Measures of Dispersion. Definition. The Range. Measures of Dispersion 3/24/2014

FUZZY FINITE ELEMENT METHOD

Probability Theory (revisited)

Numerical Heat and Mass Transfer

Statistics for Business and Economics

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Comparison of Regression Lines

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

Introduction to Regression

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing

Modal Strain Energy Decomposition Method for Damage Detection of an Offshore Structure Using Modal Testing Information

Econ Statistical Properties of the OLS estimator. Sanjaya DeSilva

The Study of Teaching-learning-based Optimization Algorithm

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

Basic Statistical Analysis and Yield Calculations

International Power, Electronics and Materials Engineering Conference (IPEMEC 2015)

Statistical analysis using matlab. HY 439 Presented by: George Fortetsanakis

Evaluation of Validation Metrics. O. Polach Final Meeting Frankfurt am Main, September 27, 2013

CALCULUS CLASSROOM CAPSULES

Module 9. Lecture 6. Duality in Assignment Problems

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M

Testing for seasonal unit roots in heterogeneous panels

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

NUMERICAL DIFFERENTIATION

The internal structure of natural numbers and one method for the definition of large prime numbers

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Lecture 3: Probability Distributions

Uncertainty on Fatigue Damage Accumulation for Composite Materials Toft, Henrik Stensgaard; Sørensen, John Dalsgaard

Uncertainty as the Overlap of Alternate Conditional Distributions

A new Approach for Solving Linear Ordinary Differential Equations

Pulse Coded Modulation

Durban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

STAT 3008 Applied Regression Analysis

Chapter 1. Probability

Finite Element Modelling of truss/cable structures

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Research Article Green s Theorem for Sign Data

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Correlation and Regression. Correlation 9.1. Correlation. Chapter 9

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Assessing inter-annual and seasonal variability Least square fitting with Matlab: Application to SSTs in the vicinity of Cape Town

REAL-TIME DETERMINATION OF INDOOR CONTAMINANT SOURCE LOCATION AND STRENGTH, PART II: WITH TWO SENSORS. Beijing , China,

Statistical Inference. 2.3 Summary Statistics Measures of Center and Spread. parameters ( population characteristics )

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands

Transcription:

Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant attrbutes of a wnd turbne to consder, s the fatgue of the wnd turbne s components. The fatgue helps to understand the magntude of the varaton on the loads and how they are affectng the turbne, ndcatng therefore the expected lfetme of the turbne. Both a measurement campagn and a smulatons on the turbne have results of the fatgue analyzes of wnd turbne s components. The problem begns when both fatgue results are beng compared. Dfferent atmospherc condtons are not always consdered n smulatons. The duraton of a load measurement campagn s often not more than 3 to 6 months and these condtons are dependent on the season, the ste and the nfluence of the locaton on the wnd turbne then the results can be obvously mslead. In ths paper a new classfcaton method s presented whch helps to compare measurement wth smulaton results. Ths s an approach not to change the certfcaton procedures, but to mprove the current classfcaton procedures. The current classfcaton method Currently the fatgue of wnd turbne components s obtaned Capture matrx Wndturbne: Example by analyzng (measurement or Wnd speed bn sze (x axs): 1 m/s 2% Turbulence bn sze (y axs) : smulaton) data wth certan 0 <3 2 16 21 57 49 49 23 1 6 2 2 4 3 5 3 5 16 36 62 117 133 105 53 32 27 11 4 6 3 1 wnd and turbulence condtons. Ths s done Turbulence 5 7 3 24 37 45 64 103 127 81 56 46 28 16 12 16 4 3 1 7 9 7 19 31 30 42 84 89 63 41 39 28 21 16 7 6 2 9 11 2 19 27 32 50 59 83 49 27 16 17 11 4 1 1 11 13 4 9 19 29 50 47 51 30 17 4 11 8 1 1 13 15 2 4 11 36 39 26 24 10 7 3 4 1 15 17 1 4 14 35 32 20 8 4 1 17 19 1 8 19 23 7 4 1 1 19 21 1 15 16 12 2 1 accordng to the normatve IEC 21 23 1 8 6 5 23 25 5 7 1 25 27 1 4 1 61400-13 [1] for mechancal 27 29 1 1 1 >29 1 load measurement on wnd Datasets : 3156 Turbulence at 15m/s: 8,18 turbne. The normatve gves a Fatgue capture matrx whch classfes a database nto wnd and Fgure 1 Capture matrx turbulence bns. Fgure 1 shows ths capture matrx where the wnd speed s ndcated on the horzontal cells n bns of 1m/s and the turbulence n vertcal cells wth bn wdths of 2%. From the matrx a sub-database wth one or more turbulences are selected for the fatgue analyss. The approach for a new classfcaton method The new method conssts of three dfferent phases. The database for the evaluatons s prepared accordng to the IEC 61400-13. Devatng from ths, 1 Mnute (and not 10 Mnute) tme seres are used n the new method to gather a sutable bg database n a reasonable measurement tme. Fgure 2 gves an overvew of the new method. Ths frst phase s dvded nto two parts. The frst part s the varables defnton, where for meteorologcal mast and turbne s quanttes several parameters are obtaned. Furthermore the fatgue of turbne s calculated from measured load quanttes. In the second part the correlaton between the parameters and the fatgue s nvestgated. The last phase s the reproducblty and the conclusons. Phase 1A Varables Wnd speed V(m/s) 0 1,5 2,5 3,5 4,5 5,5 6,5 7,5 8,5 9,5 10,5 11,5 12,5 13,5 14,5 15,5 16,5 17,5 18,5 19,5 20,5 21,5 22,5 23,5 >24.5 I(% ) 1,5 2,5 3,5 4,5 5,5 6,5 7,5 8,5 9,5 10,5 11,5 12,5 13,5 14,5 15,5 16,5 17,5 18,5 19,5 20,5 21,5 22,5 23,5 24,5 V out 0 0 21 101 232 340 493 530 541 314 181 142 101 64 43 31 16 5 0 1 0 0 0 0 0 For ths paper, the followng quanttes from a meteorologcal mast are used: Mean Wnd Speed Standard Devaton of Wnd speed

Phase 0: Select the database NOP1 1 mnute No IEC-61400-13? Yes Met Mast Turbne Fatgue IEC Phase 1A: Varables Wnd, Atmosphere Varables Loads, Wnd, Control Classfcaton Range Par Spectra EQL Met Mast: Wnd speed + Turbulence Phase 1B: Compare Varables Wnd, Atmosphere versus Fatgue Loads, Wnd, Control versus Fatgue Phase 2: Classfcaton & Results % varables wth nfluence n IEC classfcaton Learn whch condtons affects to the turbne Fatgue: turbulences + wnd speeds Phase 3: Reproducblty & Conclusons Recombne the database or change condtons n model Conclusons Fgure 2 Scheme of the new classfcaton method Turbulence Gust factor for wnd speed Trend of the wnd speed The gust factor s calculated accordng to the equatons 1 to 2 Gust = val{ max( TaL( f )) ( Envelope x ) Whereas = Equatons 1 and 2 t The followng quanttes from the nacelle measurement are used: Mean Wnd Speed Standard Devaton of Wnd speed Man frequences of the nacelle anemometer Integraton of the FFT of the nacelle anemometer accordng to equaton 3 f Percentage = value Ampltude =0 x = = Ampltude x Equaton 3

For an example the fatgue s obtaned by usng the blade root moment out of plane n every dataset. The followng fatgue quanttes are used: M40 M20 M05 Counts As the fgure 3 llustrates the parameters: M05, M10, M20 untl M60 The equvalent loads Ampltude Fgure 3 Fatgue Indexes Phase 1b Comparson of the varables. In ths phase the fatgue and the varables are compared n order to fnd: a) The man correlatons between the met mast varables and the fatgue of turbne components b) The man correlatons between the turbne varables and the fatgue of turbne components There are several methods to fnd these correlatons. To name some: lnear or polynomal regressons, statstcal approaches, neural networks, or methods based on the convoluton product (Laplacan). 1.30 1.20 1.10 1.00 0.90 Fgure 4 Correlaton between M20 and the frequency of the nacelle anemometer The correlatons can be used to estmate the fatgue: a) on the turbne s ste wth not measured load condtons. b) of same model turbnes on dfferent stes c) wth the worst case scenaro. [ MNm ] 0.70 0.50 0.30 0.10 0.00 Fgure 4 shows a correlaton to obtan M20 as equaton 4 ndcates: [ Hz ] 1.00 [ m/s ] 0.00 M 2 20 = a vmean + b f medan + c f medan Equaton 3 Phase 2 Classfcaton and Results At last the man correlatons are ndcated n the new capture matrx (fgure 5). Every cell of the capture matrx s splt n two dmensons: the horzontal have the varables found before, and the vertcal contans the bn analyss of the varables. In ths paper for clarty just the gust s ncluded although more parameters are possble.

Wndturbne: Example Wnd speed bn sze (x axs): 1 m/s Turbulence bn sze (y axs) : 2% Capture matrx V(m/s) 0 1,5 2,5 3,5 4,5 5,5 6,5 7,5 8,5 9,5 10,5 11,5 12,5 13,5 14,5 15,5 16,5 17,5 I(%) 1,5 2,5 3,5 4,5 5,5 6,5 7,5 8,5 9,5 10,5 11,5 12,5 13,5 14,5 15,5 16,5 17,5 18,5 gust gust gust gust gust gust gust gust gust gust gust gust 0 <3 3 5 5 7 7 9 9 total class 1 3 81 326 1 483 599 5 518 2 344 7 180 9 66 7 25 6 8 2 3 - class 2 2 14 45 98 114 9 49 16 6 2 1 11 class 3 3 72 262 364 40 296 174 82 35 6 3 1 class 4 7 49 74 94 102 71 37 6 6 1 1 class 5 2 4 2 3 2 1 11 13 13 15 15 17 17 19 19 21 21 23 23 25 25 27 27 29 >29 0 0 3 81 326 483 599 518 344 180 66 50 16 6 0 0 0 0 Datasets : 2672 Mean Turbulence 10,00 gust mnmum 0,82 class sze 0,29 Fgure 5 New capture matrx Ths representaton helps: a) If the model cannot be smulated wth specfc atmospherc condtons, then the capture matrx can be fltered wth the possble ones b) to study the nfluence of the atmospherc effects on the fatgue c) to understand the statstcal dstrbuton of the atmospherc condtons at the ste In the followng example (fgure 6) for two turbnes wth dfferent szes, locatons and weather condtons several varables are fltered. The parameters are gust and frequency on the nacelle anemometer (hgh or low). The fatgue of the moment out of plane s evaluated for ther correspondent webull dstrbuton, extrapolated over 20 years, the wnd speed range s n all cases from 4 m/s to 14 m/s. Dfferent wnd condtons alter the result as t can be seen of the fatgue analyss.

counts 1.00 10 8 Normal: Mbop_4_RP No_gust: Mbop_4_RP low_medan: Mbop_4_RP hgh_medan: Mbop_4_RP counts 1.00 10 8 normal: Mbop_4_RP no_gust: Mbop_4_RP low_medan: Mbop_4_RP hgh_medan: Mbop_4_RP 1.00 10 7 1.00 10 7 1.00 10 6 1.00 10 6 1.00 10 5 1.00 10 5 1.00 10 4 1.00 10 4 1.00 10 3 20.0 40.0 60.0 80.0 100.0 1.00 10 3 20.0 40.0 60.0 80.0 100.0 %MNm %MNm Fgure 6 Fatgue results Conclusons Applyng the descrbed method allows to adapt the gven measurement condtons to the smulaton and so obtan a better understandng of the wnd turbne and the effects of dfferent behavors on the loads of the turbne. Ths helps to estmate the fatgue n any turbne and to understand the dfferences between the smulatons and the measurements. References [1] Techncal Specfcaton TS IEC 61400-13: Wnd Turbne Generator Systems, Part 13: Measurement of Mechancal Loads. 1. Edton 2001-06. Internatonal Electro-techncal Commsson. Geneva, Swtzerland. 2001.