Next Article in Journal
A Simple and Robust Shock-Capturing Approach for Discontinuous Galerkin Discretizations
Next Article in Special Issue
Risk Assessment and Value of Action Analysis for Icing Conditions of Wind Turbines Close to Highways
Previous Article in Journal
Detecting West Texas Intermediate (WTI) Prices’ Bubble Periods
Previous Article in Special Issue
Response-Based Assessment of Operational Limits for Mating Blades on Monopile-Type Offshore Wind Turbines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Maintenance Optimization of Offshore Wind Turbines Based on an Opportunistic Maintenance Strategy

School of Energy Power and Mechanical Engineering, North China Electric Power University, Beijing 102206, China
*
Author to whom correspondence should be addressed.
Energies 2019, 12(14), 2650; https://doi.org/10.3390/en12142650
Submission received: 19 May 2019 / Revised: 28 June 2019 / Accepted: 8 July 2019 / Published: 10 July 2019

Abstract

:
Owing to the late development of offshore wind power in China, operational data and maintenance experience are relatively scarce. Due to the harsh environmental conditions, a reliability analysis based on limited sample fault data has been regarded as an effective way to investigate maintenance optimization for offshore wind farms. The chief aim of the present work is to develop an effective strategy to reduce the maintenance costs of offshore wind turbines in consideration of their accessibility. The three-parameter Weibull distribution method was applied to failure rate estimation based on limited data. Moreover, considering the impacts of weather conditions on the marine maintenance activities, the Markov method and dynamic time window were used to depict the weather conditions. The opportunistic maintenance strategy was introduced to cut down on the maintenance costs through optimization of the preventive maintenance age and opportunistic maintenance age. The simulation analysis we have performed showed that the maintenance costs of the opportunistic maintenance strategy were 10% lower than those of the preventive maintenance strategy, verifying the effectiveness of the proposed maintenance strategy.

1. Introduction

In recent years, the wind power industry has made rapid progress globally. Due to abundant offshore wind energy resources, the installed capacity and power generation of offshore wind power are growing rapidly [1]. However, considering the large investment and high risk involved in offshore wind power, and the influence of weather conditions, the requirements for offshore wind farms in terms of operation and maintenance vessels, spare-parts management, and other influential factors are more strict. The special operating environment of offshore wind turbines, in which the equipment is affected by natural conditions, such as typhoons, tides, waves, can accelerate the failure of unit components, and the increased failure rate of electrical and mechanical systems finally leads to lower reliability. Thus, the reliability of wind turbines has become an important issue.
The Weibull distribution is an important statistical measure in reliability engineering, known to have good adaptability and availability for various forms of failure rate simulation in relation to mechanical and electrical products. In this paper, the construction of a Weibull equation is pursued. Based on the two-parameter Weibull equation, a location parameter was introduced to create what is known as a three-parameter Weibull equation. Many researchers are interested in three-parameter Weibull equations. In [2], a three-parameter Weibull distribution was used in the asphalt‒concrete fatigue assessment of aging. In [3], optimization of distribution was studied for an inventory system established by Weibull. In [4], a three-parameter Weibull distribution was proposed in a study on the mechanical sector. In [5], the three-parameter Weibull distribution equation was applied to predict wind power density. In [6], a mixed Weibull distribution parameter estimation method was discussed. In [7], the Weibull distribution was studied by using the three best linear unbiased parameter estimation based on previous unbiased estimation and quasi-optimal linearity. This literature review reveals that the three-parameter Weibull equation has been used in many fields, including machinery, architecture, and aerospace; however, so far few applications of the three-parameter Weibull distribution have been used in fault prediction for offshore wind turbines, especially for those with limited fault data.
Moreover, access to offshore wind farms is greatly restricted by wind and wave conditions, which hinders the operation and maintenance of the unit, thus resulting in an increase in outage time, a decrease in the availability of units, and an increase in operation and maintenance costs. It is necessary to consider the impacts of sea weather conditions related to unit operation and maintenance. The impact of weather conditions on wind turbine maintenance strategies has been investigated in previous works. In China, in [8], a wind turbine reliability evaluation method was proposed. Considering the special operating environment of offshore wind turbines, in [9], the impacts of weather, accessibility, and maintenance time on the maintenance strategy of offshore wind turbines were analyzed. In [10], an offshore wind turbine maintenance strategy with minimum cost and maximum reliability was established, considering the influence of wind speed and maintenance waiting time. In [11,12], a normal behavior model was proposed to further investigate wind turbine vibration and fatigue load by using neural networks and a stochastic approach considering the influence of wind speed. In [13], a survey of stochastic models for sea state and different offshore wind farms was made according to the different maintenance schemes of offshore wind turbines, combined with the failure rate and maintenance time for various components of the unit. In [14], a novel method for simulating wind and wave conditions for offshore sites was revealed, and the result indicated that the persistence of weather windows for significant wave height values can be captured by this approach. All these tests show promising results, and the weather model has been deemed accurate enough for simulations of offshore wind parks. Moreover, in [15], a Markov model based on a statistical analysis of wind speed and wave height was established, and finally, the predictions were made for maintenance waiting time. Although there have been some related studies on the effect of weather conditions on offshore wind turbine maintenance strategies, wind speed and wave height have not yet been simultaneously considered to evaluate the accessibility of offshore wind turbine maintenance.
Meanwhile, many researchers have theoretically investigated the opportunistic maintenance strategy from a different perspective. In [16], a two-level maintenance threshold strategy for wind farms was developed, considering opportunistic maintenance and imperfect maintenance based on reliability. In [17], opportunistic condition-based maintenance for systems subjected to degradation and shocks was proposed to determine an optimal maintenance policy for multi-bladed offshore wind turbines. In [18], a new bi-objective opportunistic maintenance optimization model was proposed, and a three-phase discrete event simulation was used to evaluate the performance measures considering the stochastic behavior of wind and limited maintenance capacity. In [19], an opportunistic maintenance approach for wind farms was developed to take advantage of the maintenance opportunities, considering imperfect maintenance actions. In [20], the trade-off between wind farm configuration and the maintenance strategy was investigated by a new bi-objective redundancy and maintenance optimization model.
The primary goal of this research is to carry out studies in an effort to reduce the operating costs of offshore wind turbines. In this paper, the impacts of weather conditions on the maintenance activities are considered, the Markov chain method and dynamic time window are introduced to represent the weather conditions, and a maintenance waiting time is proposed for offshore wind turbines. In addition, the opportunistic maintenance strategy was used to optimize the key components of the maintenance of the offshore wind turbines. Furthermore, the minimum maintenance cost within the maintenance duration was deemed the optimal objective, and the preventive maintenance time and opportunistic maintenance time have been optimized for the main components of wind turbines. In particular, the novelties can be reflected in the following aspects: (1)The three-parameter Weibull method, instead of the two-parameter Weibull method, was applied to the reliability analysis of offshore wind turbines with limited fault data; (2) the Markov chain method was adopted to evaluate the prediction of maintenance waiting time, considering the wind speed and wave height simultaneously; (3) the preventive opportunistic maintenance strategy has proven to be one of the methods available to generate maintenance cost efficiency, instead of the preventive maintenance strategy.
The remainder of this paper is organized as follows. Section 2 introduces the two-parameter Weibull model of onshore wind turbines by the least squares and maximum likelihood methods, and the three-parameter Weibull model of offshore wind turbines by the correlation coefficient, probability-weighted moment, and bilinear regression methods. In Section 3, we describe the construction of a Weibull equation for offshore wind turbines based on scant sample fault data. Section 4 describes the Markov prediction method. In Section 5, we develop an opportunistic maintenance strategy based on accessibility evaluation. Section 6 gives the opportunistic maintenance simulation based on maintenance waiting time. In Section 7, we propose conclusions and recommendations from our study and suggest areas for further research.

2. Reliability Analysis Model for Wind Turbines

For wind turbines, the commonly used reliability analysis model is the Weibull distribution. This distribution is derived from the analysis of material strength by Weibull. It is a very important distribution form in reliability engineering [21]. The probability density function of the Weibull model obtained from the weakest link theory is
f ( t ) = { β α ( t γ α ) β 1 exp [ ( t γ α ) β ] ( t γ ) 0 ( t < γ )
Its distribution function is given by
F ( t ) = { 0 t f ( t ) d t = 1 exp [ ( t γ α ) β ] ( t γ ) 0 ( t < γ )
Its inefficiency function is given by
λ ( t ) = f ( t ) 1 F ( t )
In Equations (1) and (2), α is the scale parameter, β is the shape parameter, β > 0 , and γ is the position parameter, γ > 0 . γ is the threshold value for failure, β is used to describe the dispersion of the measured values, and α is related to the average measured values.
In the failure analysis of the product, β is associated with the failure mechanism of the product, and the different values of β are accompanied by different fault mechanisms. When β < 1 , the failure-rate function is a decreasing function, showing the life distribution of products under the wear-out failure period. When β = 1 , the failure-rate function is constant, representing the products in the life distribution of the random failure period. When β > 1 , the failure-rate function is an increasing function, showing the distribution of products in the period of life [22].

2.1. Weibull Model of Onshore Wind Turbines

For onshore wind turbines, because of their simple operating conditions and abundant operation and fault data, in practical applications, the common assumption is that the device fails at time γ = 0 , and the corresponding expressions, Equations (1) and (2), respectively, can be simplified as follows:
f ( t ) = β α ( t α ) β 1 exp [ ( t α ) β ] ( t 0 ) ,
F ( t ) = 1 exp [ ( t α ) β ] ( t 0 ) .
At present, the two-parameter Weibull analysis method is roughly divided into several categories, i.e., least-squares estimation, median rank regression, maximum-likelihood estimation, and method of moments. For engineering purposes, maximum likelihood is usually considered first due to its high accuracy, as compared to method of moment, which tends to be adopted only if the maximum likelihood function will be difficult to construct. In this paper, we used maximum likelihood rather than method of moment.
Practically, median-rank regression can perform well in the case of extreme values occurring when the number of failures is small. Generally, based on the least squares method, the mean variable will be replaced by the median-rank variable so as to alleviate the estimated deviation caused by the extreme value. In this study, the maintenance optimization was based on a small sample, and the extreme value was carved out; therefore, the least squares method, rather than median-rank regression, was pursued, as follows.

2.1.1. Least-Squares Estimation

By the principle of extremum, the equation can be defined as follows:
{ F a = 2 i = 0 n x i ( y i a x i b ) = 0 F b = 2 i = 0 n ( y i a x i b ) = 0
The solution equation is shown as follows:
{ a = n i = 0 n x i y i i = 0 n x i i = 0 n y i n i = 0 n x i 2 ( i = 0 n x i ) 2 b = 1 n i = 1 n y i a n i = 1 n x i
The two-logarithm transformation of Equation (5) is
ln ( ln ( 1 1 F ( t ) ) ) = β ln ( t ) β ln ( α )
Comparison with the one-element linear-regression equation gives:
{ x = ln ( t ) y = ln ( ln ( 1 1 F ( t ) ) ) a = β b = β ln ( α )
Parameters α and β can be obtained from Equation (7).

2.1.2. Maximum-Likelihood Estimation

Maximum-likelihood estimation is an effective method of parameter estimation. It has strong applicability, especially in the case of incomplete life. The basic idea is to select the undetermined parameters first so that the probability of samples appearing in the field of observation is greatest; then the sample value is taken as the unknown parameter’s point estimated value [23].
The likelihood function is
L ( t 1 , t 2 , , t n , α , β ) = i = 1 n β α ( t i α ) β 1 exp [ ( t i α ) β ]
For the convenience of calculation, the logarithm of Equation (10) is taken as a natural logarithm:
{ ln L ( t 1 , t 2 , , t n , α , β ) = n ( ln β ln α ) + ( β 1 ) i = 1 n ln ( t i α ) i = 1 n ( t i α ) β ln L α = ln L β = 0
Then, the likelihood equation is given by
{ β α i = 1 n ( t i α ) β n β α = 0 n β + i = 1 n ln ( t i α ) i = 1 n ln ( t i α ) β ln ( t i α ) = 0 .
Solving the equation gives the parameters α and β .

2.2. Weibull Model of Offshore Wind Turbines

The three-parameter Weibull distribution is used in the life tests of materials with a low stress level, which have better characteristics of curve fitting as compared to the two-parameter Weibull distribution. In the Weibull distribution function, the positional parameter is the threshold value for the failure time of targeted devices and is used to estimate the earliest failure time. Changes in positional parameters will affect the curve displacement of the probability density curve. Generally, when the unit operation time exceeds the positional parameter, the unit starts to fail.
In contrast with onshore wind farms, the conditions of operation and maintenance of offshore wind farms are extremely complicated due to tidewater, typhoons, and corrosion, which result in an increasing probability of the unit components’ failure as well as increased maintenance costs. In addition, due to the particularity of the environment, when a unit failure occurs at an offshore wind farm, maintenance personnel may not be able to reach the site for several months due to harsh weather conditions, which results in difficulty in obtaining component failure data and failure models. Therefore, it is crucial to carry out reliability prediction as soon as possible.
As part of our research, we performed a comparative analysis of the two-parameter Weibull distribution and the three-parameter Weibull distribution. The Weibull three-parameter equation prevails in short-term failure prediction and can be used effectively before t = γ. With regard to short-term failure rate prediction, an abrupt and abnormal rise in the failure rate has taken place under the two-parameter Weibull distribution. On the contrary, the failure rate curve under the three-parameter Weibull distribution displays a rational and smooth rising trend.
However, three parameters exist at the same time, so estimating the parameters of the three-parameter Weibull distribution is difficult. The least-squares method and the maximum-likelihood method are no longer applicable to the three-parameter Weibull distribution. MATLAB is a mathematics software produced by MathWorks. It is a high-level technical computing language and interactive environment for algorithm development, data visualization, data analysis, and numerical calculation, mainly including MATLAB and Simulink. The paper applied MATLAB software for data calculation and simulation.

2.2.1. Correlation Coefficient Method

The two-logarithm transformations of Equation (2) are obtained by
ln ( ln ( 1 1 F ( t ) ) ) = β ln ( t γ ) β ln ( α ) .
Equation (13) is changed into
{ a = β b = β ln ( α ) X = ln ( t γ ) Y = β ln ( t γ ) β ln ( α ) = a X + b
The relationship between X and Y is a linear relationship, as shown in Equation (14). Equation (14) shows that when the estimate of γ is correct, a linear relationship exists between X and Y, i.e., the correlation coefficient between X and Y is a maximum. The correlation coefficient r between X and Y is obtained by
r = C o v ( X , Y ) V a r [ X ] V a r [ Y ]
According to the flowchart shown in Figure 1, column vectors containing (k − 1), correlation coefficients are obtained. The step size p can be adjusted according to the actual situation and is selected as 0.01 in the paper. The correlation coefficient r initially increases and then decreases. The position of the corresponding maximum correlation coefficient is found as N, and the optimal estimation value of the position parameter is recorded as
{ t i = t i ( N 2 N ) p / 2 γ ^ = ( N 2 N ) p / 2
After obtaining the estimated value of the position parameter, we can obtain the estimated value of the parameters by least-squares or maximum-likelihood estimation.

2.2.2. Probability-Weighted Moment Method

For small samples, probability-weighted moments are unbiased and have small errors; this is suitable for the reliability assessment of offshore wind turbines with a small amount of sample data. Therefore, the probability-weighted moments of the sample are used to estimate the distribution parameters, and the accuracy is high.
The inverse probability distribution function for the random variables X of the Weibull three-parameter distribution, which is shown in Equation (17), is obtained using Equations (1) and (2).
x = γ + α [ ln ( 1 F ) ] 1 / β
From this equation, the following can be obtained:
M 1 , 0 , k = γ 1 + k + α Γ ( 1 + 1 / β ) ( 1 + k ) ( 1 + 1 / β )
Let M 1 , 0 , k = M ( k ) ; then
{ M ( 0 ) = γ + α Γ ( 1 + 1 / β ) M ( 1 ) = γ 2 + α Γ ( 1 + 1 / β ) 2 ( 1 + 1 / β ) M ( 3 ) = γ 4 + α Γ ( 1 + 1 / β ) 4 ( 1 + 1 / β )
The equation is solved, and the estimated values of parameters α , β , γ can be acquired.
{ α ^ = M ( 0 ) γ Γ { ln [ ( M ( 0 ) 2 M ( 1 ) ) ( M ( 1 ) 2 M ( 3 ) ) ] / ln ( 2 ) } β ^ = ln ( 2 ) ln [ ( M ( 0 ) 2 M ( 1 ) ) ( 2 M ( 1 ) 4 M ( 3 ) ) ] γ ^ = 4 [ M ( 3 ) M ( 0 ) M ( 1 ) 2 ] ( 4 M ( 3 ) + M ( 0 ) 4 M ( 1 ) )

2.2.3. Bilinear Regression Method

Least-squares estimation is a linear-regression analysis method for a linear equation. It can only be used to solve two-parameter estimation problems. For the estimation of three parameters, two linear equations must be combined.
Let α = t 0 1 / β ; then, Equation (2) can be transformed as follows:
ln ( ln ( 1 1 F ( t ) ) ) = β ln ( t γ ) ln ( t 0 )
( ln ( 1 1 F ( t ) ) ) 1 / β = t α γ α
Equations (21) and (22) are linearly independent; they are linear equations that can be expressed as y k = a k x k + b k ,   k = 1 , 2 and x i , y i can be depicted, respectively, as
{ y 1 = ln ( ln ( 1 1 F ( t ) ) ) a 1 = β x 1 = ln ( t γ ) b 1 = ln ( t 0 ) y 2 = ( ln ( 1 1 F ( t ) ) ) 1 / β a 2 = 1 α x 2 = t b 2 = γ α
According to the least-squares method, Equations (21) and (22) can be substituted into Equation (7) separately, and the combination can be obtained:
{ γ = t ¯ [ t 2 ¯ ( t ¯ ) 2 ] ( ln ( 1 1 F ( t ) ) ) 1 / β ¯ t ( ln ( 1 1 F ( t ) ) ) 1 / β ¯ t ¯ ( ln ( 1 1 F ( t ) ) ) 1 / β ¯ α = [ exp ( β ln ( t γ ) ¯ ln ( ln ( 1 1 F ( t ) ) ) ¯ ) ] 1 / β β = ln ( t γ ) ln ( ln ( 1 1 F ( t ) ) ) ¯ ( ln ( t γ ) ) 2 ¯ ( ln ( t γ ) ¯ ) 2 ln ( t γ ) ¯ ln ( ln ( 1 1 F ( t ) ) ) ¯ ( ln ( t γ ) ) 2 ¯ ( ln ( t γ ) ¯ ) 2
As long as the accuracy is given, the values of β and γ can be estimated iteratively, and Equation (24) is substituted to obtain the estimated value of α .

3. Analysis of Results

3.1. Verification of Reliability Analysis Method for Wind Turbines

Taking the failure shutdown data of a 3 MW wind turbine at an offshore wind farm in Jiangsu province in China as an example, the method of parameter estimation is verified in this work. The wind farm has 19 wind turbine units installed, each operating at 3 MW, and first operated in June 2017. Typical fault data were collected from the Supervisory Control and Data Acquisition system (SCADA) and selected from July 2017 to July 2018, as shown in Table 1.
The method of Weibull parameter estimation proposed in this paper is used to estimate the parameters and failure-rate function of the fault data, as shown in Table 2. The failure-rate curves of two- and three-parameter Weibull distribution are shown in Figure 2.
According to the failure-rate function estimates given in Table 2, the reliability of the wind turbine is predicted from the last fault time, as shown in Table 3.
Comparative analyses are given as follows:
1)
As shown in Table 2, for the two-parameter Weibull distribution, the least-squares, and maximum-likelihood methods are suitable to obtain similar results. For the three-parameter Weibull distribution, the correlation coefficient method, bilinear regression method, and probability-weighted moment method produce similar results.
2)
As shown in Figure 2, given the two-parameter Weibull failure-rate curve from the beginning (t = 0), and the three-parameter Weibull distribution of the failure-rate curve from the beginning (t = γ), the Weibull three-parameter equation prevails for short-term failure prediction of key defects occurring in newly operated wind turbines and can be used to effectively avoid the impact of short-term reliability prediction of atypical failure data before t = γ.

3.2. Simulation Analysis of Maintenance Costs of Wind Turbines

As shown in Figure 3, according to the wind turbine preventive maintenance strategy, the maintenance cost of component i is
C i ( T p ( i ) ) = C i c F i ( T p ( i ) ) + C i p [ 1 F i ( T p ( i ) ) ] + C N 0 T p ( i ) t f ( t ) d t + T p ( i ) T p ( i ) f ( t ) d t .
In this equation, C i c is the failure maintenance cost of component i ; C i p is the preventive maintenance cost of component i ; T p ( i ) is the preventive maintenance service life of component i ; C N is the fixed maintenance cost; C i ( T p ( i ) ) is the expected cost rate of long-term operation of component i (Chinese ¥/day).
The key components of the wind turbine are subject to Weibull distribution; the distribution parameters are shown in Table 4. The maintenance costs are shown in Table 5.
Using the two- and three-parameter Weibull reliability analysis methods, the maintenance costs are seen to vary over time, as illustrated in Figure 4, Figure 5 and Figure 6.
Using the two- and three-parameter Weibull reliability analysis methods, the optimal preventive maintenance time and maintenance cost of key components of the wind turbine can be determined; these are shown in Table 6.
It can be seen from Table 6 that the maintenance cost of the unit can be reduced by 8% using the reliability analysis method of the three-parameter Weibull model.

4. Markov Prediction Method

This method uses a mathematical model to analyze the evolution of objects [24]. In the Markov process, the time series is regarded as a stochastic process. By studying the initial probability of different states of things and the state transition matrix, the state change trend is determined, and the future state of things can be predicted. This method has been widely used in communications, biology, the social sciences, and other fields [25].

4.1. Markov Theory

4.1.1. Markov Process

If the random process { X ( t ) , t T } satisfies the following condition [26]:
(1)
If the state of the stochastic process { X ( t ) , t T } at time t is known, and the state at time t + 1 is only related to the state at time t and independent of the state before time t, then it is considered that the stochastic process { X ( t ) , t T } has Markov property.
(2)
When the state space of the stochastic process { X ( t ) , t T } is S, if for any n ≥ 2, and any t 1 < t 2 < < t n T in the condition X ( t i ) = x i , x i S , i = 1 , 2 , , n 1 , the conditional probability distribution function of X ( t n ) is equal to its probability distribution function under the condition X ( t n 1 ) = x n 1 , which is
P ( X ( t n ) x n | X ( t 1 ) = x 1 , X ( t 2 ) = x 2 , , X ( t n 1 ) = x n 1 ) = P ( X ( t n ) x n | X ( t n 1 ) = x n 1 )

4.1.2. Markov Chain

The Markov chain is a discrete-time stochastic process with Markov property, and the model is usually expressed as λ = ( S , P , Q ) , where the meaning of each element is as follows [27]:
1)
S represents a non-empty set of all possible states in a random process. State is the result of a random process occurring at a certain moment. P is a one-step state transition probability matrix. Conditional probability can be expressed as follows:
P { X ( t + 1 ) = j | X ( t ) = i } = P i j
2)
Q is the initial state probability distribution vector, Let Q = [ q 1 , q 2 , , q n ] denote the probability that the stochastic process is in state xi at time t = 0 is qi, then
i = 1 n q i = 1

4.2. Data Processing Based on Markov Method

When using the Markov chain model for predictive analysis, the random variables need to be processed as follows:
(1)
State division of random variables;
(2)
Calculation of state transition probability;
(3)
Test of “Markov property”.

4.2.1. State Division of Random Variables

For the state division of random variables, the following methods are possibilities [28]:
(1)
Sample mean value-mean variance partition method;
(2)
Ordered clustering method;
(3)
Fuzzy clustering method.
In this paper, only the first method is used and will be briefly described below.
Suppose a random observation set of random variable X is x 1 , x 2 , , x n , the sample mean is x ¯ , and the sample mean variance is s. If the absolute value | r | of the auto-correlation coefficients of this set is ≤0.30, it can be approximated as an independent and identically distributed sequence. According to the central limit theorem:
{ P { x ¯ s x ¯ x ¯ + s } = 2 Φ ( 1.0 ) 1 0.68 P { x ¯ 1.5 s x ¯ x ¯ + 1.5 s } = 2 Φ ( 1.5 ) 1 0.87

4.2.2. Calculation of State Transition Probability

Before using the Markov chain for prediction, it is necessary to estimate the transition probability of the Markov chain from the historical observation sequence.
Suppose a set of observation sequences in a random variable X is x 1 , x 2 , , x m , which contains n states, namely the state space S = { 1 , 2 , , n } . Use fij to indicate the frequency of the observation sequence from state i to state j through one step, i , j S . A matrix ( f i j ) i , j S composed of f i j ( i , j S ) is called a transfer frequency matrix. The value obtained by dividing the element of the i-th row and the j-th column in the transfer frequency matrix by the sum of the elements of the i-th row is called the maximum likelihood estimation of the transition probability, P ^ i j , i , j S , as follows:
P ^ i j = f i j j = 1 n f i j

4.2.3. Test of “Markov Property”

It is necessary to test whether the random variable sequence has “Markov property” before applying the Markov chain model to solve practical problems. Usually, a discrete sequence of the Markov chain can be tested using the χ 2 statistic.
The marginal probability is denoted as P · j , as follows:
P · j = i = 1 n f i j i = 1 n j = 1 n f i j
Then, when n is sufficiently large, the statistics obeys the χ 2 distribution with a degree of freedom of ( n 1 ) 2 , where Pij is the transition probability.
χ 2 = 2 i = 1 n j = 1 n f i j | log P i j P · j |
Given the significance level α , we can get the value of χ α 2 ( n 1 ) 2 from the lookup table, and finally calculate the value of the available statistic χ α 2 . If χ α 2 > χ α 2 ( n 1 ) 2 , it can be considered that the random sequences x 1 , x 2 , , x m have “Markov property” [29].

4.3. Markov Prediction Model

The operation and maintenance of offshore wind turbines are affected by weather conditions, such as sea wind speed and wind direction, and sea conditions, such as wave height and cycle. Among them, wind speed and wave height are the two most important factors affecting the accessibility of wind turbines [30]. As a consequence, in this paper, the accessibility of offshore wind turbines is studied mainly by considering the wind speed and wave height.
Offshore wind speed and wave height are characterized by time variation and randomness. Commonly used research methods include Gaussian statistics, the Markov method, and the auto-regressive moving average method. The Markov method was used to build the model since the method can accurately describe the long-term and seasonal wind speed and wave height distribution [31].
There are three main methods for Markov chain prediction [32]:
(1)
Markov chain prediction based on absolute distribution;
(2)
Superimposed Markov chain prediction;
(3)
Weighted Markov chain prediction.
In the case of unrestricted visual conditions, the operation and maintenance of offshore wind farms only need to consider the effects of wind speed and wave height. On the other hand, according to the monitoring of coastal sea conditions by the National Ocean Forecasting Station, the relationship between the average wind speed at sea and the effective wave height can be expressed as in Equation (33) [33]:
H s = 0.038 V 2.0 + 0.1950
where Hs is the effective wave height, and V represents the wind speed. Based on this, the paper simplifies the weather conditions at sea, and the wind speed prediction value can be obtained according to Equation (33).
Due to the fact that the wind speed and wave height at sea have obvious periodicity and seasonality and the state transition probability satisfies the “time homogeneity” characteristic, this paper uses the Markov chain prediction method based on absolute distribution to model the weather of offshore wind farms. The detailed method and steps for modeling are as follows:
(1)
State division of historical observation data of wind speed and wave height. In this paper, ∆Hs = 0.1 m, ∆V = 1 m/s or ∆Hs = 0.4 m, ∆V = 1 m/s were used as the standard interval grouping method to classify the state of the sea wind and waves.
(2)
According to the state group established by step (1), a frequency transfer matrix ( f i j ) i , j S and a Markov chain state transition probability matrix ( P ^ i j ) i , j S having a step size of one can be obtained.
(3)
Test of “Markov property”.
(4)
If the state at time t is known at which the wind speed and the wave height is i, the initial distribution can be considered to be as follows:
P ( 0 ) = ( 0 , , 0 , 1 , 0 , , 0 )
The absolute distribution of wind speed and wave height at t + 1 is
P ( 1 ) = P ( 0 ) P = ( p i 1 , p i 2 , , p i n )
The predicted state j at time t + 1 satisfies Equation (36):
p i j = max { p i j , j S } S = { 1 , 2 , , n }
Then, the wind speed and wave height at time t + k can be obtained.
(5)
Duration prediction of accessibility and inaccessible state. According to the current ship performance and sea level in China, referring to the domestic and foreign related literature [34], and combined with the weather and ocean conditions in the sea area where the offshore wind farm is located, the operation and maintenance of offshore wind farms are generally operated at wind speed V ≤ 10 m/s and wave height Hs ≤ 2 m. According to Equation (32), when the wave height is less than 2 m, the wind speed can meet the ship’s sea conditions. Therefore, the paper only needs to consider the wave height.
Suppose H is a set of accessibility states that contains multiple wind speeds and wave height states; according to Equation (36), the one-step state transition probability matrix P ˜ for H and its complement S H (denoted as H ˜ ) can be obtained. The steady-state distribution vector is π ˜ = { π ˜ H , π ˜ H ˜ } , where π ˜ H , π ˜ H ˜ represent the probability of state set H and H ˜ , respectively.
P ˜ = [ p ˜ H H p ˜ H H ˜ p ˜ H ˜ H p ˜ H ˜ H ˜ ]
Assuming that τ ( τ = 1 , 2 , , l ) is the duration of the state set H, the probability that the state set H lasts for l durations is
p { τ = l } = p { X 1 H , X 2 H , , X l H | X 0 H ˜ , X l + 1 H ˜ } = π ˜ 2 p ˜ 21 ( p ˜ 11 ) l 1 p ˜ 12
where X 1 , X 2 , X s , X l represent the states of the first, second, s-th, and l-th periods, respectively. Then, the average duration E(τ) can be obtained according to Equation (39):
E ( τ ) = i = 1 p { τ = i } i
Similarly, the probability p ˜ { τ = l } of the state H ˜ lasting for l periods of time and the average duration E ˜ ( τ ) are
p ˜ { τ = l } = π ˜ 1 p ˜ 12 ( p ˜ 22 ) l 1 p ˜ 21
E ˜ ( τ ) = i = 1 p ˜ { τ = i } i
(6)
Maintenance waiting time prediction analysis. The operation and maintenance of offshore wind turbines are affected by climate and tides, resulting in increased waiting time for maintenance and limited actual maintenance time. Assuming that personnel, vessels, spare parts, etc., are adequately prepared, the waiting time for maintenance depends only on the weather conditions.
Assuming that the unit components fail at time t and the total time required for failure repair and transportation is t n , the maintenance waiting time can be estimated using the dynamic time window. The specific steps are as follows:
1)
Determine the season w at time t , and then calculate the cumulative probability distribution based on the probability ( π w ) of each state in the season. A random number between (0, 1) is generated, and the cumulative probability distribution interval in which the random number is located is determined; thereby, a state X t , w of the wind speed and the wave height at the same moment can be obtained.
2)
If X t , w H , then initialize the duration t l o = Δ t of the reachable state and the maintenance waiting time t w a = 0 ; if not, let t l o = 0 , t w a = Δ t , where ∆t is the Markov model step size.
3)
According to the one-step state transition matrix P w , the probability vector of the state X t , w converted to the next moment state is P X , w , and X t , w = P X , w . The next time state X t + 1 , w is obtained according to the method in step (1).
4)
If X t + 1 , w H , update the reachable state duration t l o = t l o + Δ t , t = t + 1 ; otherwise, t l o = 0 , t w a = t l o + Δ t , t = t + 1 .
5)
Repeat steps (3) and (4) until t l o t n , where t w a is the waiting time for this maintenance activity.

5. An Opportunistic Maintenance Strategy Based on Accessibility Evaluation

A wind turbine is a complex multi-component system. Once a component fails, the functions of other components will also be affected, resulting in overall shutdown and the need for maintenance of the unit. The paper intends to study a maintenance strategy for offshore wind turbines in combination with equipment fault repair and preventive maintenance.

5.1. Opportunistic Maintenance Strategy

The opportunistic maintenance strategy refers to the fact that, when a component in the system is repaired or replaced due to a malfunction or other reason, it also provides an opportunity for preventive maintenance or replacement of other components, considering the fact that the maintenance of multiple parts at the same time can effectively save fixed maintenance costs and increase availability.
The principle of the opportunistic maintenance strategy for failure repair combined with preventive maintenance of components is shown in Figure 7.

5.2. Opportunistic Maintenance Model

5.2.1. Opportunistic Maintenance Probability Density Function

According to the opportunistic maintenance strategy in Figure 7, the expected replacement rate of the component consists of three parts, complete replacement, incomplete repair, and minimum maintenance. This paper only considers the opportunistic maintenance strategy for a complete replacement. According to the update process theory, when the component is repaired for a period of time, the maintenance rate is close to constant, which is the reciprocal of the expected replacement cycle of the component. According to the opportunistic maintenance strategy shown in Figure 7, the replacement rate λ i of one component consists of three parts, as shown in Equation (42):
λ i = 1 / E i ( T ) = λ c i + λ o i + λ p i
where λ i = P i ( f ) × λ i ; λ o i = P i ( o ) × λ i ; λ p i = P i ( p ) × λ i are the failure replacement rate, opportunistic replacement rate, and preventive replacement rate of one component, respectively. P i is the replacement probability and E i ( T ) is the expected life of the component.
Since the opportunistic replacement of component i is caused by the failure replacement or preventive replacement of other components, the opportunistic replacement probability density function of the component i can be approximated as an exponential distribution function, as follows:
g i ( t , T o ( i ) ) = δ i exp [ δ i ( t T o ( i ) ) ]
δ i = j i N ( λ c i + λ p i )
where g i ( t , T o ( i ) ) represents the opportunistic replacement probability density function of the component i , δ i is the exponential distribution function, and N is the number of key components of the unit.

5.2.2. Key Component Replacement Probability of the Wind Turbine

According to the opportunistic maintenance strategy shown in Figure 7, the probability of each replacement is as follows:
P i ( c ) = 0 T o ( i ) f i ( t ) d t + T o ( i ) T p ( i ) f i ( t ) × G ¯ i ( t , T o ( i ) ) d t ,
P i ( o ) = T o ( i ) T p ( i ) F ¯ i ( t ) × g ¯ i ( t , T o ( i ) ) d t ,
P i ( p ) = F ¯ i ( T p ( i ) ) × G ¯ i ( T p ( i ) , T o ( i ) ) ,
F ¯ i ( t ) = 1 0 t f i ( u ) d u ,
G ¯ i ( t , T o ( i ) ) = 1 T o i t g i ( u , T o ( i ) ) d u ,
where f i ( t ) is the failure probability density function of component i, which obeys the Weibull distribution.

5.2.3. Opportunistic Maintenance Optimization Model

During one replacement cycle of a wind turbine component, the maintenance cost of component i consists of replacement costs, downtime losses, and rental vessel costs. The expected repair cost of component i during a replacement cycle can be expressed as follows:
E i ( C ) = [ C c i + ( t w a + t n ) C e ] P i ( c ) + C o i P i ( o ) + [ C p i + ( t w a + t n ) C e ] P i ( p ) + C N
where C c i , C p i , C o i represent the costs of failure maintenance, preventive maintenance, and opportunistic maintenance, respectively, of component i; C e is the loss per unit downtime, and C N is the rental vessel cost.
During one replacement cycle of a wind turbine component, the expected life of component i can be expressed as follows:
E i ( T ) = 0 T o ( i ) t × f i ( t ) d t + T p ( i ) × F i ¯ ( T p ( i ) ) × G i ¯ ( T p ( i ) , T o ( i ) ) + T o ( i ) T p ( i ) t [ f i ( t ) × G i ¯ ( t , T o ( i ) ) + F i ¯ ( t ) × g t ( t , T o ( i ) ) ] d t
According to the above analysis, by optimizing the opportunistic maintenance age and preventive maintenance service age of each key component of the unit to minimize the average cost rate, the following optimization model is established:
min Z ( T o , T p ) = i = 1 N E i ( C ) E i ( T ) s . t .    0 < T o < T p t l o t n
where Z ( T o , T p ) is the average cost rate of the unit.

6. Simulation Analysis

6.1. Accessibility Evaluation Simulation Analysis

Due to a lack of weather and ocean data for offshore wind farms, this paper takes the historical observation data from an observation station in the Yangtze River estuary (30° N, 128° E) in the East China Sea as an example to model the sea weather.
The data on wind speed and wave height in the article were recorded every half hour from 0:00 on 1 January 2014. The 17,306 data points recorded in 2014 were used to make a forecast analysis of wind speed and wave height in 2015.

6.1.1. Wind Speed and Wave Height State Prediction

Using the Markov method, the wind speed and wave height can be divided into n states, and the state set is S = {1,2,⋯,n}. In this paper, H s = 0.1 m, 0.2 m, 0.3 m, 0.4 m are used as grouping criteria to predict the 100 data points of wind speed and wave height in January 2015. The results are shown in Figure 8, Figure 9, Figure 10 and Figure 11.
The difference between the predicted data and the observed data and the correlation coefficient are used as criteria of the quality of the predicted data. The correlation coefficients of the observed and predicted data of average wind speed and average wave height are calculated separately, as shown in Table 7.
In Table 7, the data of the wind speed and wave height when Δ H s = 0 are observed data; r is the correlation coefficients of the observed and predicted data; d is the absolute error of prediction data and observation data of the wave height, expressed as in Equation (52):
d = i = 1 100 | ( H s ) p r e d i c t i v e   v a l u e ( H s ) o b s e r v e d   v a l u e |
Taking the average error of the wave height prediction data as a standard to measure the advantages and disadvantages of different state grouping methods, the average error, corresponding to H s = 0.1 m, 0.2 m, 0.3 m, 0.4 m, is shown in Figure 12.
From Figure 12, when Δ H s = 0.1 m, the average error of wave height prediction data and observation data reached the minimum, which means the state grouping method of Δ H s = 0.1 m is superior to the other three groups. As a consequence, this state grouping method is used as a grouping criterion for the state in the following.

6.1.2. Duration Prediction of Wind Speed and Wave Height

Markov models are established for each of the four seasons. The occurrence probability of the accessible state H and the inaccessible state H ˜ in different seasons is shown in Figure 13.
The duration of each state under different seasons is shown in Table 8. One duration is the step size of the Markov model, which is 0.5 h.
Take the inaccessible state H ˜ in January 2015 as an example, the observed and predicted values of the probability distribution and the cumulative distribution of the state duration are presented in Figure 14, which shows the accuracy of the Markov chain model.
According to Figure 13 and Table 8, the occurrence probability of accessibility in summer is larger than in any other seasons, and the average duration is the longest, indicating that the wind speed and wave height in summer are more suitable for sea. On the contrary, the occurrence probability of the inaccessible state in winter is the largest, indicating that the wind speed and wave height are not suitable, and the bad weather will last for a long time.

6.1.3. Maintenance Waiting Time Estimate

Taking the maintenance waiting time in January 2015 as an example, the relationship between the maintenance waiting time and the state of the reachable state is shown in Figure 13. One duration is set to be 0.5 h. It can be seen that the waiting time for maintenance increases with the duration of the accessible state.
As seen in Figure 15, at the initial moment, there is a large difference between the predicted value and the observed value when Δ H s = 0.1 m. As the duration increases, the observations gradually approach the predicted value, and finally, only a small gap is observed. When Δ H s = 0.4 m, the predicted value is far from the observed value at the initial moment; as the duration increases, the gap decreases, but in the end, it still maintains a large gap.

6.2. Simulation Analysis of the Opportunistic Maintenance Strategy

In this paper, a 3 MW wind turbine in an offshore wind farm in China is used as an example to optimize the opportunistic maintenance strategy for the critical components of the unit, including the gearbox, the main bearing, and the generator. According to the analysis of the statistical data of failure record of the unit, the failure time of the component follows the Weibull distribution, and the specific parameters are shown in Table 9 below, where 1, 2, 3, and 4 represent the blade, gearbox, main bearing, and generator, respectively; α and β are the shape parameters and scale parameters of the components, respectively; unit downtime loss C e = 720 yuan/h; and the ship chartering fee C N = 150,000 yuan/day.
When the traditional age-preventive maintenance strategy is adopted, the preventive maintenance age of the critical components of the unit is optimized to minimize the average cost rate. The results after optimization are shown in Table 10.
When using a multi-component opportunistic maintenance strategy, the opportunistic maintenance age, and preventive maintenance age of the components are optimized to minimize the average cost rate. The results are shown in Table 11.
From the comparison of the results in Table 10 and Table 11, it can be seen that the maintenance cost of the wind turbine using the opportunistic maintenance strategy is 10% lower than that of the preventive maintenance strategy.
Although the proposed model gains competitiveness in terms of the unit maintenance cost, its limitations are as follows:
(1)
The research on the maintenance cost was aimed at the key components of a single wind turbine rather than multiple wind turbines.
(2)
Influential factors, such as the availability of vessels, the adequacy of the spare parts, and the meteorological conditions, were linked to maintenance waiting time. In this paper, wind speed and wave height were considered in the determination of maintenance waiting time, assuming vessels and spare parts were available.

7. Conclusions

Before concluding, one additional point needs to be discussed. Because the operating environment of offshore wind turbines is complex, wind is a crucial factor that should be considered in optimal maintenance decisions. Moreover, the randomness of wind impacts on the power generated by wind turbines, which contributes to unsatisfactory energy costs. In this paper, the influence that fluctuations in wind speed had on the costs of offshore wind farming was considered in the determination of maintenance waiting time. From the perspective of reducing the entire-life costs of an offshore wind farm, we paid attention to how much maintenance costs will be reduced with the introduction of the opportunistic maintenance strategy. Actually, in this sense, how the wind fluctuation will act on the unit cost is not the research priority in this study. Indeed, the influence of wind fluctuation and grid impact on the unit is worth intensively study in further research.
In this paper, we attempted to analyze the characteristics of operation data of offshore wind turbines and investigated reliability analysis methods for offshore wind turbines based on limited fault data. Considering the influence of weather factors, such as wind speed and wave height, we studied maintenance waiting time prediction methods for offshore wind turbines. Combining failure maintenance and preventive maintenance, we proposed an opportunity-based offshore wind turbine maintenance strategy. The main study results are as follows:
(1)
The construction of a Weibull equation for offshore wind turbines was based on a small amount of sample fault data. Different to [8], based on the construction of a two-parameter Weibull equation, a three-parameter Weibull equation was proposed. The results show that the maintenance costs can be reduced by 8% with the adoption of a three-parameter Weibull model, and the fitting curve and failure rate short-term prediction of the three-parameter Weibull distribution is superior to the two-parameter Weibull distribution where there are limited fault data;
(2)
A maintenance waiting time prediction method was introduced for offshore wind turbines. The Markov chain method and dynamic time window were used to describe wind speed and wave height, and a maintenance waiting time prediction model was established. Different to [9], the impacts on maintenance waiting time arising from wind speed and wave height were considered in this paper. The results show that the deviation of the predicted value of the wave height obtained by grouping interval with ∆Hs = 0.1 m was the smallest, which was close to the true value;
(3)
Combining failure maintenance and preventive maintenance, an opportunistic maintenance strategy was presented for offshore wind turbines. The minimum expected maintenance cost was regarded as an objective function to optimize the opportunistic maintenance time and preventive maintenance time. Compared with [10], the opportunistic maintenance strategy reduces maintenance duration and decreases the maintenance waiting time and downtime, thereby reducing maintenance costs. The results show that the maintenance cost was reduced by 10% under the opportunistic maintenance strategy for offshore wind turbine maintenance, which verified the effectiveness and superiority of the opportunistic maintenance strategy for offshore wind turbine maintenance.

Author Contributions

Data curation, X.H.; Investigation, L.X.; Methodology, S.L.; Supervision, X.R.; Visualization, S.L.; Writing—original draft, L.X.

Funding

This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

Acknowledgments

We would like to thank AIP Author Services for English language editing.

Conflicts of Interest

The authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

References

  1. Mahmood, S. Maintenance logistics organization for offshore wind energy: Current progress and future perspectives. Renew. Energy 2015, 77, 182–193. [Google Scholar]
  2. Li, X.; Liang, N.; Cheng, Z.J. Estimation of fatigue failure of asphalt concrete based on three-parameter Weibull distribution. Chongqing Jiaotong Univ. (Nat. Sci.) 2017, 1, 178. [Google Scholar]
  3. Sanni, S.S.; Chukwu, W.I.E. Inventory model with three-parameter Weibull deterioration, quadratic demand rate and shortages. Am. J. Math. Manag. Sci. 2016, 35, 159. [Google Scholar] [CrossRef]
  4. Deng, B.; Jiang, D.; Gong, J. Is a three-parameter Weibull function really necessary for the characterization of the statistical variation of the strength of brittle ceramics? Eur. Ceram. Soc. 2017, 20, 123. [Google Scholar] [CrossRef]
  5. Chaurasiya, P.K.; Ahmed, S.; Warudkar, V. Study of different parameters estimation methods of Weibull distribution to determine wind power density using ground based Doppler SODAR instrument. Alex. Eng. J. 2017, 37, 140. [Google Scholar] [CrossRef]
  6. Chai, H.; Shin, J.; Lim, T.; Trans, J. Failure rate calculation using the mixture Weibull distribution. Kim. Korean Inst. Electr. Eng. 2017, 66, 500. [Google Scholar]
  7. Xu, X.L. Statistical inference for incomplete data three-parameter Weibull distribution. Shanghai Teachers Univ. (Nat. Sci.) 2002, 31, 22. [Google Scholar]
  8. Fu, Y.; Xu, W.X.; Liu, L.J. An opportunistic maintenance strategy for offshore wind turbine based on accessibility evaluation. electric power. Wind Eng. 2016, 49, 74–80. [Google Scholar]
  9. Liu, L.J.; Fu, Y.; Ma, S.W. Preventive maintenance strategy for offshore wind turbine based on reliability and maintenance priority. Proc. CSEE 2016, 36, 5732–5740. [Google Scholar]
  10. Zheng, X.X.; Zhao, H.; Liu, L.J. A Combined maintenance strategy for offshore wind turbine considering accessibility. Power Syst. Technol. 2014, 38, 3030–3036. [Google Scholar]
  11. Pedro, G.L.; Luis, V.T.; Matthias, W. Normal behaviour models for wind turbine vibrations: Comparison of neural networks and a stochastic approach. Energies 2017, 10, 1944–1957. [Google Scholar]
  12. Pedro, G.L.; Iván, H.; Matthias, W. Fatigue load estimation through a simple stochastic model. Energies 2014, 7, 8279–8293. [Google Scholar]
  13. Monbet, V.; Aillot, P.; Prevosto, M. Survey of stochastic models for wind and sea state time series. Probab. Eng. Mech. 2007, 22, 113–126. [Google Scholar] [CrossRef]
  14. Matha, D.; Scheu, M.; Muskulus, M. Validation of a Markov-based weather model for simulation of O&M for offshore wind farms. In Proceedings of the 22th International Offshore and Polar Engineering Conference, Rhodes, Greece, 17–22 June 2012; pp. 463–468. [Google Scholar]
  15. Scheu, M.; Matha, D.; Hofmann, M.; Muskulus, M. Maintenance strategies for large offshore wind farms. Eng. Proc. 2012, 24, 281–288. [Google Scholar] [CrossRef]
  16. Zhang, C.; Gao, W.; Guo, S.; Li, Y.L.; Yang, T. Opportunistic maintenance for wind turbines considering imperfect, reliability-based maintenance. Renew. Energy 2017, 103, 606–612. [Google Scholar] [CrossRef]
  17. Mahmood, S.; Maxim, F.; Christophe, B. An opportunistic condition-based maintenance policy for offshore wind turbine blades subjected to degradation and environmental shocks. Reliab. Eng. Syst. Saf. 2015, 142, 463–471. [Google Scholar]
  18. Hadi, A.; Karim, A.; Morteza, A. Multi-objective opportunistic maintenance optimization of a wind farm considering limited number of maintenance groups. Renew. Energy 2016, 88, 247–261. [Google Scholar]
  19. Ding, F.F.; Tian, Z.G. Opportunistic maintenance for wind farms considering multi-level imperfect maintenance thresholds. Renew. Energy 2012, 45, 175–182. [Google Scholar] [CrossRef]
  20. Atashgar, K.; Abdollahzadeh, H. Reliability optimization of wind farms considering redundancy and opportunistic maintenance strategy. Energy Convers. Manag. 2016, 112, 445–458. [Google Scholar] [CrossRef]
  21. Hallinan, A.J. A review of the Weibull distribution. Qual. J. Technol. 1993, 25, 85–93. [Google Scholar] [CrossRef]
  22. Zhu, W.Y. Design of Mechanical Reliability; Press of Shanghai Jiaotong University: Shanghai, China, 1992; pp. 62–66. [Google Scholar]
  23. Zhu, Y.H.; Tai, S.C.; Sun, Y.Y. Applied Mathematical Statistics; Press of Wuhan University of Hydraulic and Electrical Engineering: Wuhan, China, 1999; pp. 68–72. [Google Scholar]
  24. Lu, D.J. Stochastic Processes and Applications; Tsinghua Press: Beijing, China, 1986; pp. 100–115. [Google Scholar]
  25. Liu, C.H. Stochastic Processes; Press of Huazhong University of Science and Technology: Wuhan, China, 2008; pp. 42–45. [Google Scholar]
  26. Liu, Y.F. Research on Reliability of Software System Based on Markov Chain Method; Changchun University of Science and Technology: Changchun, China, 2005. [Google Scholar]
  27. Sheldon, M.R. Stochastic Processes, 2nd ed.; John Wiley and Sons: Hoboken, NJ, USA, 1996. [Google Scholar]
  28. Qian, M.P. Application on Stochastic Processes; Peking University Press: Beijing, China, 1998. [Google Scholar]
  29. Luo, J.Y. Economic Statistical Analysis Methods and Forecasts; Tsinghua Press: Beijing, China, 1987; pp. 347–348. [Google Scholar]
  30. Van Bussel, G.J.W.; Zaaijer, M.B. Reliability, availability and maintenance aspects of large-scale offshore wind farms, a concepts study. In Proceedings of the MAREC, Newcastle, UK, 17–18 April 2001; pp. 5–9. [Google Scholar]
  31. Hagen, B.; Simonsen, I.; Hofmann, M.; Muskulus, M. A multivariate Markov weather model for O&M simulation of offshore wind parks. Energy Procedia 2013, 35, 137–147. [Google Scholar]
  32. Lyding, P.; Faulstich, S.; Hahn, B. Reliability of the electrical parts of wind energy systems-a statistical evaluation of practical experiences. In Proceedings of the EPE Wind Energy Chapter Symposium, Staffordshire University, Stafford, UK, 15–16 April 2010. [Google Scholar]
  33. Feng, W.B.; Peng, X.L.; Zhang, S.L. Analysis on wave condition of south ocean of China: No 14. In Proceedings of the China Ocean Science and Technology Symposium, Hohhot, China, 16–17 May 2009. [Google Scholar]
  34. Sinha, Y.; Steel, J.A. A progressive study into offshore wind farm maintenance optimization using risk based failure analysis. Renew. Sust. Energ. Rev. 2015, 42, 735–742. [Google Scholar] [CrossRef]
Figure 1. Correlation coefficient method.
Figure 1. Correlation coefficient method.
Energies 12 02650 g001
Figure 2. Failure-rate curves of Weibull distributions. (a) Comparison of failure rate for three- and two-parameter Weibull distributions. (b) Partial enlarged drawing of maximum likelihood estimation and least-squares estimation.
Figure 2. Failure-rate curves of Weibull distributions. (a) Comparison of failure rate for three- and two-parameter Weibull distributions. (b) Partial enlarged drawing of maximum likelihood estimation and least-squares estimation.
Energies 12 02650 g002
Figure 3. Preventive maintenance strategy based on age.
Figure 3. Preventive maintenance strategy based on age.
Energies 12 02650 g003
Figure 4. The cost rate of blade maintenance based on two- and three-parameter Weibull: the left-hand side (a) shows the blade maintenance cost rate of the two-parameter Weibull, and the right-hand side (b) shows the three-parameter Weibull.
Figure 4. The cost rate of blade maintenance based on two- and three-parameter Weibull: the left-hand side (a) shows the blade maintenance cost rate of the two-parameter Weibull, and the right-hand side (b) shows the three-parameter Weibull.
Energies 12 02650 g004
Figure 5. The cost rate of gearbox maintenance based on two- and three-parameter Weibull: the left-hand side (a) shows the gearbox maintenance cost rate of the two-parameter Weibull, and the right-hand side (b) shows the three-parameter Weibull.
Figure 5. The cost rate of gearbox maintenance based on two- and three-parameter Weibull: the left-hand side (a) shows the gearbox maintenance cost rate of the two-parameter Weibull, and the right-hand side (b) shows the three-parameter Weibull.
Energies 12 02650 g005
Figure 6. The cost rate of generator maintenance based on two- and three-parameter Weibull, the left figure (a) shows the generator maintenance cost rate of the two-parameter Weibull, and the right figure (b) shows the three-parameter Weibull.
Figure 6. The cost rate of generator maintenance based on two- and three-parameter Weibull, the left figure (a) shows the generator maintenance cost rate of the two-parameter Weibull, and the right figure (b) shows the three-parameter Weibull.
Energies 12 02650 g006
Figure 7. Opportunistic maintenance strategy.
Figure 7. Opportunistic maintenance strategy.
Energies 12 02650 g007
Figure 8. The comparison of predicted and observed values of wave height, wave height prediction curve with ∆Hs = 0.1 m.
Figure 8. The comparison of predicted and observed values of wave height, wave height prediction curve with ∆Hs = 0.1 m.
Energies 12 02650 g008
Figure 9. The comparison of predicted and observed values of wave height, wave height prediction curve with ∆Hs = 0.2 m.
Figure 9. The comparison of predicted and observed values of wave height, wave height prediction curve with ∆Hs = 0.2 m.
Energies 12 02650 g009
Figure 10. The comparison of predicted and observed values of wave height, wave height prediction curve with ∆Hs = 0.3 m.
Figure 10. The comparison of predicted and observed values of wave height, wave height prediction curve with ∆Hs = 0.3 m.
Energies 12 02650 g010
Figure 11. The comparison of predicted and observed values of wave height, wave height prediction curve with ∆Hs = 0.4 m.
Figure 11. The comparison of predicted and observed values of wave height, wave height prediction curve with ∆Hs = 0.4 m.
Energies 12 02650 g011
Figure 12. The average error of the wave height prediction data with ∆Hs = 0.1 m, 0.2 m, 0.3 m, 0.4 m.
Figure 12. The average error of the wave height prediction data with ∆Hs = 0.1 m, 0.2 m, 0.3 m, 0.4 m.
Energies 12 02650 g012
Figure 13. The probability of each state occurring in different seasons.
Figure 13. The probability of each state occurring in different seasons.
Energies 12 02650 g013
Figure 14. The duration prediction of state H ˜ .
Figure 14. The duration prediction of state H ˜ .
Energies 12 02650 g014
Figure 15. Maintenance waiting time and duration diagram.
Figure 15. Maintenance waiting time and duration diagram.
Energies 12 02650 g015
Table 1. Failure times.
Table 1. Failure times.
Sequence12345678910111213
Failure time (h)4560456846604756487948994904496749884990499549964998
Sequence14151617181920212223242526
Failure time (h)4999503450965137513852785289536754675567567857785879
Sequence27282930313233343536373839
Failure time (h)5889590060676078617862786378647865786678678968796888
Sequence40414243444546474849505152
Failure time (h)6978698869907089711572457345735674567468756875697709
Sequence53545556
Failure time (h)7809788979007908
Table 2. Parameters estimated.
Table 2. Parameters estimated.
MethodWeibull Distribution ParametersFailure-Rate Functions
αβγ
Least-squares estimation4.9979 × 10318.012-4.7388 × 10−66t17.023
Maximum-likelihood estimation4.9941 × 10317.963-7.2376 × 10−66t16.976
Correlation coefficient method4.298 × 1021.48744621.8379 × 10−4 (t − 4457)0.484
Probability-weighted moment4.069 × 1021.38344723.2298 × 10−4 (t − 4486)0.393
Bilinear regression method3.859 × 1021.14445131.4007 × 10−3 (t − 4543)0.124
Table 3. Short period prediction on failure rates.
Table 3. Short period prediction on failure rates.
Failure RateTwo-ParameterThree-Parameter
Least-Squares EstimationMaximum-Likelihood EstimationCorrelation Coefficient MethodProbability-Weighted Moment Method
λ ( 5378 ) 0.01290.01310.00500.0047
λ ( 8760 ) 51.772651.14630.01050.0086
Table 4. Weibull distribution parameters for major components of wind turbines.
Table 4. Weibull distribution parameters for major components of wind turbines.
Component α /Day β
Blade30002
Gearbox24003
Generator33003
Table 5. Maintenance costs for the major components of wind turbines.
Table 5. Maintenance costs for the major components of wind turbines.
Component C i c / ¥ C i p / ¥ C N / ¥
Blade112,00028,00035,000
Gearbox152,00038,000
Generator100,00025,000
Table 6. Optimal maintenance time and maintenance cost of unit components.
Table 6. Optimal maintenance time and maintenance cost of unit components.
ComponentTwo-ParameterThree-Parameter
T p ( i ) / Day C i ( T p ( i ) ) T p ( i ) / Day C i ( T p ( i ) )
Blade277051.7344249049.8622
Gearbox166468.4917150563.5657
Generator247038.2659228431.0392
Table 7. Comparison of observation and prediction data.
Table 7. Comparison of observation and prediction data.
Δ H s ( m ) V(m/s) H s ( m ) d(m) r
09.93004.561801.0000
0.19.27563.936016.07640.8742
0.29.77624.196017.72520.8235
0.38.11543.684018.63770.8712
0.48.76243.524019.37320.8525
Table 8. The duration of each state under each season.
Table 8. The duration of each state under each season.
StatePersistence
SpringSummerAutumnWinter
H0.6240.7860.5700.455
H ˜ 0.3760.2140.4300.545
Table 9. Weibull distribution parameters and maintenance costs of critical components.
Table 9. Weibull distribution parameters and maintenance costs of critical components.
Component Number α i / Day β i γ i / Day C 1 i /Yuan C 2 i /Yuan
125292.8465672,000168,000
222282.5136912,000228,000
340362.591360,00090,000
431542.236600,000150,000
Table 10. Optimization results of age maintenance strategy.
Table 10. Optimization results of age maintenance strategy.
ComponentPreventive Maintenance Age Tp/DayMinimum Cost Rate/¥ Day−2Total Minimum Cost Rate/¥ Day−1
Blade1816184.8898782.6744
Gearbox1423300.5747
Main bearing2822108.3843
Generator2284188.8256
Table 11. Optimization results of the opportunistic maintenance strategy.
Table 11. Optimization results of the opportunistic maintenance strategy.
ComponentOpportunistic Maintenance Age To/DayPreventive Maintenance Age Tp/DayTotal Minimum Cost Rate/¥ Day−1
Blade7212709702.5834
Gearbox6391812
Main bearing10013000
Generator8972430

Share and Cite

MDPI and ACS Style

Xie, L.; Rui, X.; Li, S.; Hu, X. Maintenance Optimization of Offshore Wind Turbines Based on an Opportunistic Maintenance Strategy. Energies 2019, 12, 2650. https://doi.org/10.3390/en12142650

AMA Style

Xie L, Rui X, Li S, Hu X. Maintenance Optimization of Offshore Wind Turbines Based on an Opportunistic Maintenance Strategy. Energies. 2019; 12(14):2650. https://doi.org/10.3390/en12142650

Chicago/Turabian Style

Xie, Lubing, Xiaoming Rui, Shuai Li, and Xin Hu. 2019. "Maintenance Optimization of Offshore Wind Turbines Based on an Opportunistic Maintenance Strategy" Energies 12, no. 14: 2650. https://doi.org/10.3390/en12142650

APA Style

Xie, L., Rui, X., Li, S., & Hu, X. (2019). Maintenance Optimization of Offshore Wind Turbines Based on an Opportunistic Maintenance Strategy. Energies, 12(14), 2650. https://doi.org/10.3390/en12142650

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop