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.
Power plants of different types have in common the requirements for the maintenance of various critical parts such as turbines, pumps, gearboxes, generators, etc. The requirements for such kind of services and maintenance are demanding with the need for high reliability with low downtime. Traditionally, the preventive type of maintenance is applied where the service and control are performed periodically regardless of the state of the power generating device or the specific component. This approach has the advantage of simplicity and ease of implementation. However, there are various disadvantages related to the preventive maintenance approach such as the unnecessary servicing of certain parts and the inability of the prediction of the need for servicing other parts more frequently. These disadvantages could drive the total cost of the maintenance, and in some cases, they can lead to failures of the components and costly downtime.
In recent decades, wind power has become increasingly popular as a way of power generation, and it is already replacing a significant part of the traditional power plants such as coal plants and gas plants. At the same time, the maintenance of wind turbines is becoming a critical issue.
It is an aim of the present invention to mitigate at least some of the disadvantages of the condition monitoring and maintenance currently applied for power generation plants and in particular wind farms.
According to the present invention, there is provided a method for predictive monitoring of the condition of wind turbines, the method comprising the steps of:
In an embodiment the method comprises the step of sending an alarm in the case that at least one monitoring feature exceeds the corresponding threshold value.
In an embodiment of the method the predetermined statistics comprise one or the combination of linear interpolation, average, and standard deviation of the differential data.
In an embodiment of the method 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.
In an embodiment 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.
In an embodiment of the method the step of processing SCADA data further comprises averaging of the SCADA data over a predefined period of time. Preferably, the averaging is performed over the period ranging from one hour to one year.
In an embodiment of the method 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.
In an embodiment of the method the wind turbine component is one of main bearing, generator, hydraulic oil system, inverter, transformer and gearbox.
In an embodiment of the method the preselected time period over which the SCADA data is collected is longer than one year.
In an embodiment of the method the monitoring feature is one-day average or one-month average of the differential data.
In an embodiment of the method the monitoring feature is one-week average or one-month average of the differential data, and/or the reference feature is one-week or one-month average of the differential data.
In an embodiment of the method the monitoring feature and/or the reference feature is the slope of the linear interpolation of the differential data.
In an embodiment of the method 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.
In another aspect of the invention, there is provided 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.
The invention will be better understood with the aid of the description of embodiments given by way of example only and illustrated by the figures, in which:
As shown in
After acquiring the SCADA data, the data is processed in step 13 of the method. In one preferred embodiment shown in
In one embodiment, 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. In one embodiment, the SCADA data points originally recorded every 10 minutes are averaged to daily resolution. In general, 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. In one preferred embodiment, 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. In one embodiment, for each pre-processed or raw SCADA data point di=(ti, Ti) the differential data point δi=Ti−
Step 14, shown in
Various time periods as a function of time are visualized in
The example of the embodiment of the feature extraction is shown in the flow diagram in
In one non-limiting example, the features may be categorized in the three categories:
After the feature extraction is performed, there is the testing step 15 as shown in
In one preferred embodiment 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.
In one embodiment, the test may be performed as follows:
Further examples are provided in the following paragraphs using definitions:
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). For example, 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); For instance 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 Tin 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). There might for various reasons be no data point in the time series for a date, this because all points might have been filtered out in the power filtering stage so that there are no data points left for the entire day.
Below is a list of condition of alerts. If any of these conditions are satisfied then the alarm is sent.
mean(T,ts,7)>mean(T,ts,30)+x_week*std(T,ts,10) and count(T,ts,7)>1, where x_week=0.9. Short term condition 1:
mean(T,ts,1)>mean(T,ts,30)+x_day*std(T,ts,10), where x_day=1.9. Short term condition 2:
slope(T,ts,14)>x_2 week*std(T,ts,10) and count(T,ts,14)>3, where x_2 week=0.07. Slope condition 2 weeks:
slope(T,ts,30)>x_month*std(T,ts,10) and count(T,ts,30)>8, where x_2month=0.05. Slope condition 1 month:
slope(T,ts,90)>x_3month*std(T,ts,10) and count(T,ts,90)>count(T,ts,30)+10 and slope(T,ts,30)>0, where x_3month=0.026. Slope condition 3 months:
slope(T,ts,180)>x_6month*std(T,ts,10) and count(T,ts,180)>count(T,ts,90)+20 and slope(T,ts,30)>0 and slope(T,ts,90)>0, where x_6month=0.013. Slope condition 6 months:
slope(T,ts,270)>x_9month*std(T,ts,10) and count(T,ts,270)>count(T,ts,180) and slope(T,ts,30)>0 and slope(T,ts,90)>0 and slope(T,ts,180)>0, where x_9month=0.004. Slope condition 9 months:
mean(T,ts,1)>mean (T,ts,30) and mean (T,ts,30)−mean (T,ts,90)>y_month*std(T,ts,10) and mean (T,ts,90)−mean (T,ts,180)>y_3month*std(T,ts,10), where y_month=0.21 and y_3month=0.15. Long term condition:
The method according to the invention offers several advantages compared to the traditional methods for the power generation assets.
Using the method according to the invention it is possible to identify the failure patterns and to release alarms at an early stage, which allows for the more efficient risk management of the critical components' failures and increase of the energy output by minimizing costly downtimes.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/IB2020/056187 | 6/30/2020 | WO |