1. Introduction
With the characteristics of stable wind energy resources, large power generation, and easy consumption, offshore wind energy has developed rapidly in recent years and has a broad market prospect [
1,
2]. Offshore wind energy has significant advantages compared with onshore. (1) Energy reserves are more than three times those of onshore wind. (2) Average available hours of offshore wind power are more than 3000 h per year, which is higher than the 2000 h of onshore wind power [
3]. This happens because the offshore wind is more stable and the average wind speed is higher, which can meet the minimum wind speed requirements for longer time than onshore wind power. (3) The offshore wind turbine (OWT) has little impact on human life and is close to the coastal electricity concentration area, which is convenient for the supply of wind energy.
Despite the great potential of offshore wind power, its high risk of failure and serious consequences of failure pose challenges to safe and efficient energy development [
4,
5]. As an important wind energy collection and conversion equipment, OWTs are in a harsh climate, and the waves, tides, and other factors are complex and variable, resulting in the risk of failure changing with the increase of running time [
6,
7]. The design life of OWT is 20 years, and its cost of the operation and maintenance (O&M) occupy approximately 25–30% of the 75–90% of the investment costs or whole energy generation cost [
8,
9]. It can be explained by the following: (1) insufficient accessibility due to the unpredictable weather and the remote location; (2) compared with onshore ones, more failures owing to the oceanic environment, for example, storm surge, sea wave, sea ice; (3) extra inventory expense, specific vessels, and technicians are needed.
Risk control measures should analyze the influences of various factors such as safety, economy, and cost [
10]. Under the traditional periodic maintenance and time-based maintenance methods, if the system conditions are acceptable when performing preventive maintenance work, it will cause a waste of the remaining life of the system, and if the system deteriorates faster than expected, a failure may occur. Recently, researches related to data-based wind turbines condition monitoring are proposed. Pandit et al. [
11] presented a reference power curve using Supervisory Control And Data Acquisition (SCADA) datasets from a healthy turbine, which is developed by using a Gaussian Process and then was compared with a power curve from an unhealthy turbine. Yang et al. [
12] studied a new condition monitoring method with the help of the concept of the transmissibility of Frequency Response Functions. Experiment of verification showed that the new technique is effective not only in damage detection but in damage location in certain conditions. With the continuous maturity of state monitoring, storage and analysis technologies, condition-based maintenance (CBM) is a solution of this problem and has become a hot issue in the field of risk control [
13,
14]. The CBM method focuses on the specific operating conditions of the system, evaluates the true state of the monitored equipment, overcomes the blindness of maintenance, and can effectively reduce accidental failures caused by insufficient maintenance and waste of resources due to over-maintenance. Simulation methods are usually used in studying the effect of CBM in any type of system [
15]. Scheu et al. [
16] presented a transparent risk-based methodology for the prioritization of OWT systems toward the application of condition monitoring systems. Calculated results contained information related to various kinds of wind turbines and substructure concepts. This work shows details of paths that leads to critical failure modes, which is the basis of designing a condition monitoring system. Shafiee et al. [
17] proposed a CBM optimization model for a multi-blade fan system. The calculation results show that the CBM method can effectively reduce the maintenance cost. Ghamlouch et al. [
18] proposed a method of condition-based maintenance strategy mainly adopting option theory, which solved the maintenance interval optimization problem of complex systems considering the uncertainty of production and deterioration. Verbert et al. [
19] analyzed the multi-component problem of the research object and proposed a new CBM optimization method. Based on the analysis of economic correlation and structural correlation between components, the rational planning of maintenance tasks can be made to achieve system-level maintenance program optimization.
In terms of multi-component systems such as OWTs, whenever a component of the OWT system needs to be repaired, other components will be given the opportunity to repair in advance and form a certain number of group maintenance plans [
20]. Opportunistic maintenance, often abbreviated as OM, has superiority in the following aspects: (1) it is critical for OWTs that joint maintenance processes only spend one portion of the fixed cost; (2) vessels, technicians, and tools are all maintenance resources that can also be in common use. Whereas, the preventative maintenance (usually abbreviated as PM) should not be conducted when the previous maintenance date results in the loss of assemblies. Regarding Remaining Useful Life (RUL), by comparing the maintenance costs corresponding to different grouping schemes, the best repair solution could be determined. Song et al. [
21] introduced a two-stage framework to optimize the operation of offshore wind farms. Computational models were developed for opportunistic condition-based maintenance. Minimum maintenance costs are calculated and used to determine the schedule of periodic inspections. Pandit et al. [
22] presented a support vector regression-based pitch curve, which was adopted in anomaly detection exploration. The comparative results in consideration of a binned pitch curve showed that the blade pitch curve closely follows the binned pitch curve, but above the rated wind speed. Zhou and Yin [
23] formulated a dynamic opportunistic condition-based maintenance strategy. The varying maintenance lead-time effect on the maintenance decision was analyzed, which remarkably affects the maintenance cost. The comparative analysis illustrated the capability of the proposed strategy. Lu et al. [
24] proposed a CBM optimization method for OWT systems in consideration of opportunities. An artificial neural network (ANN) was employed to analyze the RUL based on the system monitoring data. A comparative study was performed, and the feasibility of the presented approach was proved.
The objective of this paper is to achieve an optimized economy and availability by an opportunistic maintenance policy. First, we calculate the preventive repair threshold and the opportunity maintenance threshold from the parameters such as the lost shutdown loss, maintenance cost, and failure probability of the component. By comparing the support vector machine (SVM) fault prediction results with the maintenance thresholds, the components that need to be repaired at that time are identified, and finally the maintenance operation optimization scheme in the entire calculation interval is formed. The core idea of this method is to obtain the optimal repair time and maintenance combination by analyzing the correlation between components, so as to reduce the number of repairs and reduce the repair cost.
The following sections of the paper are organized as this:
Section 2 presents the OWT system description. The fault diagnosis model based on SVM is introduced in
Section 3. A case study is performed in
Section 4, including the discussion of the results. Conclusions are addressed in the last section.
2. System Description of a Generic Offshore Wind Turbine
The basic function of the OWT is to use the rotor blade system to absorb wind energy and convert it into mechanical energy; then, transfer the mechanical energy to the generator system through the transmission system; then, convert it into electric energy by the generator; and finally, output the electric energy through the grid to complete the whole process. In the entire process of converting to electricity, OWTs can be divided into a certain amount of assemblies according to various functions. In this paper, the OWTs system [
25] is divided into four main assemblies according to the various functions (
Figure 1), which are the more relevant ones from a failure rate point of view [
26].
A rotor can be divided into three components, including rotor blades, rotor bearings, and a rotor hub. A rotor system functions to absorb and transmit the wind energy. A generator is installed inside the nacelle. This equipment is used to convert mechanical energy to electrical energy and adapt the output energy from the wind turbine to the grid. The gearbox functions to transform high-torque to low-torque and transform the low speed of the main shaft to the high speed of the generator. A pitch system is a mechanism that turns the blade, or part of the blade, in order to adjust the angle of attack of the wind.
Carroll et al. [
27] addressed three parameters of offshore wind turbines including maintenance cost, repair time, and failure rate. The results indicate that the pitch system and generator are the largest contributors to failure rates, and the gearbox has the highest average cost per failure. The downtime of the gearbox and rotor are the longest, which can be explained by the high installation position and weight.
3. Support Vector Machine Classification Algorithm
Fault analysis based on big data is the frontier method for the risk identification of offshore wind turbines. Support vector machine (SVM) classifiers is an effective data processing technology, which is based on the structural risk minimization principle solving a quadratic programming problem [
28]. By applying the kernel function, the method can map the data into a higher dimensional input space and construct an optimal hyperplane. In this study, “ν”-Soft Margin Support Vector Classifiers (SVC), a class of ν-SVM, is applied, which can predetermine the fraction of training sample that are support vectors [
29].
Given a training dataset as
, the goal is to find a function
to construct the classifier, taking the form
where
is the weight vector;
is the nonlinear function; and
is the bias term. For the original Support Vector Classifier (SVC) algorithm, the optimization process can be described as follows
subject to
The regularization constant
determines the trade-off of a large margin and noise tolerance, which is a hyperparameter that needs to be determined. Data points that are closer to the hyperplane and affect the position and orientation of the hyperplane are called support vectors (as shown in
Figure 2). In the
-SVC algorithm, the hyperparameter
is replaced by the formulation of constant
to control the number of margin errors and the support vectors. The optimization problem becomes
subject to
Note that
does not appear in the objective function. Instead, a parameter
and an additional variable
are added to be optimized. The new hyperparameter
is an upper bound on the fraction of margin errors and a lower bound on the fraction of Support Vectors (SV) [
29]. To derive the dual, the Lagrangian can be constructed as
where
are Lagrange multipliers. The optimal solution is given by the saddle point of the Lagrangian. Thereby, the conditions are
Substituting Equation (7) into Equation (6), the
-SVC is obtained in a Lagrange dual form
subject to
In the above equations, the dot product of
is substituted by kernel functions. In this study, the Gaussian kernel is used to map the data into a higher dimensional input space, meanwhile avoiding the dimensionality curse
is the hyperparameter that determines the width of the Gaussian kernel. Then, the regression function Equation (1) can be rewritten as follows
Only the training data belonging to the Support Vectors will affect the classifier, and the number of Support Vectors is controlled by the hyperparameters . It should be noted that the selection of hyperparameters and is significant to avoid overfitting. They can be selected by the cross-validation procedure, finding the support vector classifiers with a good performance on data not yet observed.
4. Opportunistic CBM Optimization for OWTs
The downtime of OWT brings an opportunity for joint maintenance activities. This section proposed a maintenance method that aims to take advantage of the opportunity to determine the optimized maintenance plan aiming to obtain the minimum long-term operational expenditure.
The Opportunistic CBM method presented in this study is based on an assumption that an OWT is regarded as a multi-component system consisting of n independent units. A binary model is applied to describe the condition of each unit, which means that the components are either in a functional or failure state.
4.1. Reliability Threshold Calculation
According to the Opportunistic CBM strategy, the maintenance decision for each component can be concluded.
The probability density
f(t) can be described by
where
is the corresponding survivor function. The failure rate
can be derived as
Since
R(0) = 1, then
where
expresses the cumulative failure risk of component
i in maintenance cycle
j, implying that the number of corrective maintenance (CM) for component
i in each maintenance cycle is equal to
, which is considered as the PM threshold.
Since the number of CM during Δ
ti,j is
, assuming that the maintenance times for component
i during the operation period is
Mi, the mean maintenance cost per unit time of this component is defined as follows
where
τi,j is the PM duration and
τ’i,j is the CM duration.
τi,j consists of the repair time and travel time, and
τ’i,j is usually longer than
τi,j because CM requires extra logistic time, resulting in more downtime losses. The influence of changeable marine environment is considered in the presented model. In Equation (15),
Wt is the waiting time for an appropriate sea condition, which satisfied the requirement of offshore operation.
CCM(i) is the CM cost, including the recovery expense and the downtime losses.
C′PM(i) and
CD(i) represent the imperfect PM expense and unit downtime losses, respectively.
C′PM(i) is related to age reduction factor
σ i,j and PM cost
CPM(i), which can be evaluated as
The individual optimal Ri,j(t) can be obtained by minimizing . Obtaining threshold Ri,j(t) is the basis of the presented approach.
4.2. The Calculation of Assemblies’ Remaining Life
In this paper, the SVM method proposed in
Section 3 is implemented to address the collected OWT monitoring information. A dataset with 80% of the historical data is used as incentive data to complete the training of the model. The remaining 20% of the historical data are used as verification data to prove the accuracy of the model prediction. In practice, temperature and vibration are closely related to system failure. Therefore, this research selects these two feature quantities as the regressors. The trained SVM model is employed to identify the real-time status of the components. If the threshold
Ri,j(t) is reached, a PM activity should be performed.
In the modeling process, “kernel trick” is used to approximate the upper bound of the generalization error by minimizing the norm of the weights in the feature space. The Gaussian kernel function is used to determine the nonlinear mapping of regressors. Compared to the polynomial kernels, the Gaussian kernel has less numerical difficulty, and only one hyperparameter needs to be determined, which is easier to choose. In this way, the model structure is chosen with a trade-off between empirical errors in training data and model complexity. By using the kernel dot product trick, the curse of dimensionality, which usually happens in other black-box modeling methods such as ANN, can be avoided, and a unique global solution can be solved through a convex optimization problem.
4.3. Determination of the Maintenance Group
For an OWT system, when units
k and
i are jointly repaired, the corresponding cost preservation can be evaluated as
where
is the downtime cost preservation resulting from the joint maintenance, which can be defined as
When the component
i is recovered in advance, the unexpected failures would be reduced. The cost saving for the decline of random malfunctions is
, which can be expressed as
After one PM activity is performed in advance, every planned maintenance time should be rearranged. Assuming
represents the previous maintenance time and
represents the latest one, the overall time switch should be
in which
represents the originally planned number of PM, and
is the number of PM activities after the update.
However, the substantial RUL is wasted when the unit is repaired in advance and the penalty cost is
The cost savings can be achieved based on Equations (17)–(21).
The maintenance activities of components can be regarded as a finite set Φ in terms of a multi-component system. When a maintenance opportunity appears, each maintenance combination Φ
1, Φ
2, …, Φ
l is a subcollection of Φ, satisfying
The size of Φ explodes as the amount of element rises, making the solution process extremely complicated. When a CM or planned PM is performed on component
, all the opportunistic maintenance (OM) combination can be identified, and the expense preservation should be
The OM combination corresponding to the maximum C(Φl) can be identified as the optimized solution.
4.4. Rolling-Horizon Update
If
ti,1 is the
j-th maintenance execution time of unit
i (
i = 1, 2, …,
N), thus, the following equations can be established:
and
where
tbegin is the initial time, which is generally set as 0.
Δti,1 is the interval of the first PM activity.
Wt is the waiting time for a weather window. After completing a maintenance activity, the failure rate of the corresponding component is renewed. If one preventative maintenance or failure replacement is conducted for unit
I, a new life period need to be started.
4.5. Maintenance Schedule Determination and Maintenance Cost Calculation
The total maintenance expenditure of OWT can be derived as follows when the last inspection is accomplished.
represents the inspection interval, and
represents the inspection time of the OWT operation process. The entire cost at inspection time
is denoted as
Ct, which can be calculated by
where N represents the amount of OWT assemblies. The complete calculation process of opportunistic CBM optimization based on SVM is presented in
Figure 3.