EP4127465A1 - Procédé de surveillance prédictive de l'état d'éoliennes - Google Patents

Procédé de surveillance prédictive de l'état d'éoliennes

Info

Publication number
EP4127465A1
EP4127465A1 EP20737569.2A EP20737569A EP4127465A1 EP 4127465 A1 EP4127465 A1 EP 4127465A1 EP 20737569 A EP20737569 A EP 20737569A EP 4127465 A1 EP4127465 A1 EP 4127465A1
Authority
EP
European Patent Office
Prior art keywords
time period
data
monitoring
wind turbine
wind
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
EP20737569.2A
Other languages
German (de)
English (en)
Other versions
EP4127465C0 (fr
EP4127465B1 (fr
Inventor
Gianmarco PIZZA
Eskil JARLSKOG
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fluence Energy LLC
Original Assignee
Fluence Energy LLC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fluence Energy LLC filed Critical Fluence Energy LLC
Publication of EP4127465A1 publication Critical patent/EP4127465A1/fr
Application granted granted Critical
Publication of EP4127465C0 publication Critical patent/EP4127465C0/fr
Publication of EP4127465B1 publication Critical patent/EP4127465B1/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D17/00Monitoring or testing of wind motors, e.g. diagnostics
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00Control
    • F05B2270/30Control parameters, e.g. input parameters
    • F05B2270/303Temperature
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00Control
    • F05B2270/40Type of control system
    • F05B2270/404Type of control system active, predictive, or anticipative
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction

Definitions

  • the present invention relates to a method for predictive monitoring of the condition of power generating plants, and in particular, it relates to the condition monitoring of wind turbines.
  • a method for predictive monitoring of the condition of wind turbines comprising the steps of: - selecting at least one wind turbine inside a wind farm and at least one component of the wind turbine;
  • SCADA supervisory control and data acquisition
  • - processing SCADA data comprising calculating differential data, wherein the differential data is a difference between the temperature values of the selected turbine component of the selected wind turbine and an average temperature of the selected wind turbine component in at least two wind turbines in the wind farm;
  • the extraction of the features comprises calculating at least one predetermined statistic of the differential data during the monitoring time period, and saving the predetermined statistics as a monitoring feature;
  • the method is characterized in that the step of the extracting the features further comprises calculating at least one predetermined statistic of the differential data during a reference time period, and that the reference time period is at least partially overlapping the monitoring time period and saving the predetermined statistic as a reference feature.
  • the threshold value is a function of at least one reference feature.
  • the method comprises the step of sending an alarm in the case that at least one monitoring feature exceeds the corresponding threshold value.
  • the predetermined statistics comprise one or the combination of linear interpolation, average, and standard deviation of the differential data.
  • the testing step further comprises calculating the total number of the differential data points during the monitoring time period and/or the reference time period.
  • the method comprises the step of filtering SCADA data to eliminate low-power data, wherein the low-power data is determined based on a preselected power threshold.
  • the step of processing SCADA data further comprises averaging of the SCADA data over a predefined period of time.
  • the averaging is performed over the period ranging from one hour to one year.
  • the average temperature of the wind turbine component is calculated over all the wind turbines in the entire wind farm that are of the same model.
  • the wind turbine component is one of main bearing, generator, hydraulic oil system, inverter, transformer and gearbox.
  • the preselected time period over which the SCADA data is collected is longer than one year.
  • the monitoring feature is one- day average or one-month average of the differential data.
  • the monitoring feature is one- week average or one-month average of the differential data
  • the reference feature is one-week or one-month average of the differential data.
  • the monitoring feature and/or the reference feature is the slope of the linear interpolation of the differential data.
  • the slope is calculated over at least one of the time periods with a duration of half a month, one month, three months, six months and nine months.
  • a computer- readable storage medium comprising instructions, which when executed by a computer cause the computer to carry out the steps of the method of any of embodiments described above.
  • Fig. 1 shows a flow diagram of the method for condition monitoring in accordance with one embodiment of the invention
  • Fig. 2 shows a flow diagram of the selection step in accordance with one embodiment of the invention
  • Fig. 3 shows a flow diagram of the processing step in accordance with one embodiment of the invention
  • Fig. 4 shows a flow diagram of the filtering of data in the processing step in accordance with one embodiment of the invention
  • Fig. 5 shows a flow diagram of the example of feature extraction in accordance with one embodiment of the invention
  • Fig. 6A and Fig. 6B show a schematic illustration of various time periods.
  • FIG. 1 shows a flow diagram of the method for condition monitoring in accordance with one embodiment of the invention.
  • the first step in the method may be a selection step 11 as shown in more detail in Figure 2.
  • a specific wind farm is selected.
  • at least one component of the wind turbines in the wind farm is selected.
  • this component may be, but not limited to: main bearing, generator, hydraulic oil system and gearbox of the wind turbine.
  • the final step 23 of the selection process 11 at least one wind turbine to be monitored is selected.
  • the method allows also to monitor multiple turbines and components of the wind turbines in parallel.
  • step 12 of the method comprises acquiring supervisory control and data acquisition (SCADA) data of the wind farm, or of the one part of the wind farm comprising a minimum of two wind turbines.
  • SCADA supervisory control and data acquisition
  • the SCADA data across whole farm is collected.
  • the SCADA data contains operation data of the wind turbine or wind farm.
  • the SCADA system is a control system of a plant or turbine used for high-level process supervisory management.
  • the SCADA system performs a supervisory operation over a variety of the devices and components of the controlled system.
  • SCADA data contains vast information collected by the sensors inside the wind turbine.
  • the SCADA data comprises temperature values of the at least one component of the wind turbine during a preselected time period.
  • the SCADA data comprises the temperature values of all the critical components of the wind turbine.
  • the SCADA data is collected in time periods which may range from few minutes to few years.
  • the time intervals between the measurement points in SCADA data varies depending on the required accuracy, for example, the intervals could range from a few seconds or multiple minutes.
  • SCADA data may record the series of information as a tuple (T, t, P, Temp), where T is name of the turbine, t is the time of the measurement, P is power of the turbine at time t, and Temp is the temperature of the selected component at time t.
  • SCADA data is normally available to the users, so advantageously the method may not need any additional measurements or installation of new devices.
  • step 13 the data is processed in step 13 of the method.
  • the SCADA data may go through the three steps of: filtering low-power data 31, time averaging of the data 32, and calculation of differential data 33.
  • the first two steps, 31 and 32, may serve to improve the quality and to make processing faster.
  • step 13 is consisting only of step 33.
  • the step of filtering out low power data is shown in more detail in flow chart in Figure 4.
  • power level is tested in step 41. In the case that the power level is lower than the certain percentage of maximum power the data point is discarded.
  • the certain percentage of the maximum power may be for example in the range of 0 to 90%.
  • Step 42 checks if the turbine has been consistently producing power (producing more than a certain threshold) during a certain, so-called, consistency time period prior to the time of the measurement of the that particular data point. If this is not the case the data point is discarded. Otherwise the data is saved as filtered SCADA data in step 43.
  • step 32 time-averaging, which may have the following steps, is performed: defining duration and resolution of time intervals; calculation and storage of the time intervals; for each time interval the average of all data points in the interval is calculated and saved.
  • the SCADA data points originally recorded every 10 minutes are averaged to daily resolution.
  • the resolution parameters that may be used are time unit (hour or day), the number of time units (positive integer) and moving average (positive integer). The process of averaging may be repeated for every turbine of the wind farm that is monitored.
  • Step 33 which is the calculation of the differential data, may comprise the calculation of the difference between temperature of the selected wind turbine component of the selected wind turbine and an average temperature of the selected wind turbine component in at least two wind turbines in the wind farm.
  • the average temperature of the wind turbine component is calculated over all the wind turbines in the wind farm that are of the same model.
  • the differential data point T t - T t is calculated, where t t - denotes the time corresponding to the i:th data point; T t - denotes the component temperature corresponding to the i:th data point; T t - denotes the average farm component temperature of the turbines which have the same model as the one that is described by the data point.
  • the process of calculating the differential data may be repeated for every turbine in the wind farm that is analyzed.
  • Step 14 shown in Figure 1 shows the feature extraction step.
  • the extraction of the features comprises calculating at least one predetermined statistic of the differential data during a monitoring time period and saving the predetermined statistics as a monitoring feature.
  • the monitoring time period may be selected in the selection step 11, or it may be selected at a later stage, but before the step of the feature extraction.
  • the monitoring time period is one day or one month.
  • the step of the feature extraction further comprises calculating at least one predetermined statistic of the differential data during a reference time period, wherein the reference time period is at least partially overlapping the monitoring time period and saving the predetermined statistic as a reference feature.
  • the predetermined statistic may comprise one or the combination of: linear interpolation, average and standard deviation of the differential data.
  • the predetermined time period is usually a long period lasting months or years.
  • the monitoring time period is usually short period lasting, for example, days or months.
  • the reference time period is at least partially overlapping the monitoring time period, and it can be, for example, days or months.
  • the examples of the durations are provided for illustration purposes, and they are not limiting.
  • step 51 the monitoring time period MTP is selected, which may be for example one week as shown in Figure 6A.
  • step 52 five reference time periods (RTP1, RTP2, RTP3, RTP4, and RTP5) are selected in step 52.
  • the selected reference time periods have durations of 0.5, 1, 3, 6 and 9 months as shown in Figure 6B.
  • the monitoring time period and the reference time periods end at the same time ts as shown in Figure 6B.
  • Linear interpolation, average and standard deviation are calculated and saved for each reference time period in step 53.
  • step 54 the linear combinations of the statistics is calculated and saved as the reference features in final step 55. Using the monitoring time period similar monitoring features may be calculated.
  • the features may be categorized in the three categories:
  • - Short term increase features for example, one-day average and one- month average;
  • - Slope features for example, slopes of the interpolated lines of the periods with a duration of 0.5, 1, 3, 6 and 9 months;
  • the feature may be the linear combination of at least two of the short term increase feature, the slope feature and the consistent long-term increase feature.
  • testing step 15 is testing if the at least one monitoring feature exceeds a corresponding threshold value.
  • the threshold value is a function of the at least one reference feature.
  • all the thresholds are based on the standard deviation of the pre-processed time series containing the temperature of the component which has been monitored. For example, an alarm is triggered if the following condition is satisfied: mean_day > mean_month + x*std_month where mean_day is a monitoring feature corresponding to statistic of average value of the component temperature during the monitoring time of one day; mean_month is reference feature corresponding to the statistic of mean value temperature during the period of one month before the monitoring time; std_month is reference feature corresponding to the statistic of standard deviation of the temperature during the one month before the monitoring time; and x is an adjustable parameter, which is dependent on the component.
  • test may be performed as follows:
  • mean(T, ts, n) denotes the mean of the pre-processed time series, before the alerting stage, starting at the timestamp corresponding to n days before ts (exclusive) and ending at ts (inclusive).
  • mean(T, 2019-01-10, 7) is the mean of the time series over the days 2019-01-04, 2019-01-05, ..., 2019-01-10 which is 7 days in total
  • slope(T, ts, n) denotes the slope of the linear interpolation of the pre- processed time series, before the alerting stage, starting at the timestamp corresponding to n days before ts (exclusive) and ending at ts (inclusive);
  • slope(T, 2019-01-10, 7) is the slope of the linear interpolation of the time series over the days 2019-01-04, 2019-01-05, ..., 2019-01-10 which is 7 days in total. If the time series would have the same values for all of these days then the slope would be equal to 0; std(T, ts, n) denotes the standard deviation of the pre-processed time series for all turbines of the same model as T in the wind farm, before the alerting stage, starting at the timestamp corresponding to x days before ts (exclusive) and ending at ts (inclusive); count(T, ts, n) denotes the total number of data points that is in of the pre- processed time series, before the alerting stage, starting at the timestamp corresponding to x days before ts (exclusive) and ending at ts (inclusive).
  • count(T, ts, n) denotes the total number of data points that is in of the pre- processed time series, before the alerting
  • the method according to the invention offers several advantages compared to the traditional methods for the power generation assets.

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Wind Motors (AREA)

Abstract

L'invention concerne un procédé de surveillance prédictive de l'état d'éoliennes. Ledit procédé comprend les étapes : la sélection d'au moins une éolienne à l'intérieur d'un parc éolien et au moins un composant de l'éolienne ; l'acquisition de données SCADA comprenant des données opérationnelles du parc éolien pendant une période de temps présélectionnée, les données SCADA comprenant des valeurs de température du ou des composants de l'éolienne pendant la période de temps présélectionnée ; le traitement des données SCADA incluant le calcul de données différentielles, les données différentielles étant une différence entre les valeurs de température du composant d'éolienne sélectionnée de l'éolienne sélectionnée et une température moyenne du composant d'éolienne sélectionné dans au moins deux éoliennes dans le parc éolien ; la définition d'une période de temps de surveillance pour surveiller le composant ; l'extraction des caractéristiques incluant le calcul d'au moins une statistique prédéterminée des données différentielles pendant la période de temps de surveillance, et la sauvegarde de la statistique prédéterminée en tant que caractéristique de surveillance ; l'essai si au moins une caractéristique de surveillance dépasse une valeur seuil.
EP20737569.2A 2020-06-30 2020-06-30 Procédé de surveillance prédictive de l'état d'éoliennes Active EP4127465B1 (fr)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/IB2020/056187 WO2022003397A1 (fr) 2020-06-30 2020-06-30 Procédé de surveillance prédictive de l'état d'éoliennes

Publications (3)

Publication Number Publication Date
EP4127465A1 true EP4127465A1 (fr) 2023-02-08
EP4127465C0 EP4127465C0 (fr) 2024-02-07
EP4127465B1 EP4127465B1 (fr) 2024-02-07

Family

ID=71527848

Family Applications (1)

Application Number Title Priority Date Filing Date
EP20737569.2A Active EP4127465B1 (fr) 2020-06-30 2020-06-30 Procédé de surveillance prédictive de l'état d'éoliennes

Country Status (6)

Country Link
US (1) US20230184223A1 (fr)
EP (1) EP4127465B1 (fr)
AU (1) AU2020455928A1 (fr)
CA (1) CA3197065A1 (fr)
IL (1) IL301469A (fr)
WO (1) WO2022003397A1 (fr)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116203910B (zh) * 2023-04-27 2023-07-07 三峡智控科技有限公司 一种基于异构同源的风机状态映射与判断系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2223048A4 (fr) * 2007-12-11 2014-12-03 Vestas Wind Sys As Système et procédé permettant la détection de performance
DE102009004385B4 (de) * 2009-01-12 2010-11-25 Repower Systems Ag Verfahren und Anordnung zum Überwachen einer Windenergieanlage
EP2267305B1 (fr) * 2009-06-24 2016-01-13 Vestas Wind Systems A/S Procédé et système pour contrôler le fonctionnement d'une éolienne
US8190394B2 (en) * 2011-05-31 2012-05-29 General Electric Company System and methods for monitoring oil conditions of a wind turbine gearbox
DE102016117190A1 (de) * 2016-09-13 2018-03-15 fos4X GmbH Verfahren und Vorrichtung zum Überwachen eines Zustands wenigstens einer Windkraftanlage und Computerprogrammprodukt
CN109779846B (zh) * 2019-01-11 2020-04-14 北京京运通科技股份有限公司 基于风电机组温度的故障预警方法

Also Published As

Publication number Publication date
CA3197065A1 (fr) 2022-01-06
AU2020455928A1 (en) 2023-08-24
IL301469A (en) 2023-07-01
EP4127465C0 (fr) 2024-02-07
EP4127465B1 (fr) 2024-02-07
US20230184223A1 (en) 2023-06-15
WO2022003397A1 (fr) 2022-01-06

Similar Documents

Publication Publication Date Title
Schlechtingen et al. Wind turbine condition monitoring based on SCADA data using normal behavior models. Part 2: Application examples
CN101555806B (zh) 发电厂生产实时参数分类报警和识别辅助方法
CN103912448B (zh) 一种区域风电场机组功率特性监测方法
CN108074197B (zh) 风机故障数据分析系统的控制方法
Cui et al. An anomaly detection approach using wavelet transform and artificial neural networks for condition monitoring of wind turbines' gearboxes
CN108072524B (zh) 一种风电机组齿轮箱轴承故障预警方法
WO2009016020A1 (fr) Système de surveillance d'éolienne
Butler et al. A feasibility study into prognostics for the main bearing of a wind turbine
CN107654342A (zh) 一种考虑湍流的风电机组功率异常的检测方法
EP4127465B1 (fr) Procédé de surveillance prédictive de l'état d'éoliennes
de Andrade Vieira et al. Wind turbine condition assessment based on changes observed in its power curve
CN112633614A (zh) 一种基于特征提取的实时故障程度诊断系统及方法
CN115406483A (zh) 一种水电机组故障识别方法、装置、设备及可读存储介质
CN111103136A (zh) 基于scada数据分析的风机齿轮箱故障检测方法
CN111176226A (zh) 一种基于运行工况的设备特征参数报警阈值自动分析方法
CN117119783B (zh) 一种模块待机功耗的控制方法
US20210182749A1 (en) Method of predicting component failure in drive train assembly of wind turbines
Sarma et al. Early life failure modes and downtime analysis of onshore type-III wind turbines in Turkey
Lima et al. Wind turbine failure prediction using SCADA data
CN111879522A (zh) 基于时序概率的汽轮机运行监测和故障判别方法及系统
Osadciw et al. Wind turbine diagnostics based on power curve using particle swarm optimization
CN106650154B (zh) 一种基于热力性能指标的核电常规岛设备可靠性监测方法
McKinnon et al. Investigation of anomaly detection technique for wind turbine pitch systems
CN114323260A (zh) 汽轮机组全工况振动故障自动诊断方法及计算机可读介质
Chacón et al. SCADA data analytics for fault detection and diagnosis of wind turbines

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: UNKNOWN

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20221024

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)
INTG Intention to grant announced

Effective date: 20230927

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602020025333

Country of ref document: DE

U01 Request for unitary effect filed

Effective date: 20240207

U07 Unitary effect registered

Designated state(s): AT BE BG DE DK EE FI FR IT LT LU LV MT NL PT SE SI

Effective date: 20240212