1. Introduction
Fatigue deterioration is a common phenomenon in offshore structures such as ships, offshore platforms, offshore wind turbines (OWTs), etc. [
1]. Fatigue cracks usually start at pre-existing defects or stress concentrations and then gradually propagate, eventually leading to structural failure [
2,
3,
4]. This failure may cause great damage to the ocean environment, human life, and property [
5,
6,
7]. Thus, it is important to ensure the reliability of structural systems.
Many classical system reliability analysis methods have been widely adopted in previous studies. One is fault tree analysis (FTA), which focuses on analyzing the causes and paths of system failures. Volkanovski et al. [
8] proposed a novel approach for assessing the reliability of power systems using FTA, which can identify the key elements in the power systems. Ding et al. [
9] established a new method to assess the reliability of residual heat removal systems, and FTA was employed to analyze the logical relationship among events. Another is event tree analysis (ETA), which can analyze all the possible sequences of events and predict the failure probability of the system. Lu et al. [
10] evaluated the system reliability of cable-stayed bridges using ETA incorporating the
-unzipping method. Researchers also have used ETA in reliability analysis considering inspection and maintenance interventions. Kim and Frangopol [
11] quantified the relationship between damage detection delay (associated with reliability) and inspection strategies for fatigue deteriorating structures by ETA. Zou et al. [
1] proposed a holistic approach for inspection and maintenance decision-making and obtaining optimum fatigue reliability using ETA. Thus, ETA is employed as a basic tool in this study.
Inspection and maintenance can gain new information on the state of the system, thereby reducing uncertainties in system reliability assessment. Bayesian inference can update the prior probabilistic model based on the new information and quantify the effect of inspection and maintenance on system reliability [
12]. Classical reliability calculation methods, such as the first-order reliability method (FORM), are efficient in evaluating system reliability but inconvenient for Bayesian updating [
13]. To integrate Bayesian updating into system reliability calculation methods, researchers have proposed many methods. Bayesian network (BN) is one of the most popular models due to its intuitive nature and ability to address correlations among multiple variables. Straub [
14] established a model for computationally efficient and robust reliability analysis of probabilistic deterioration processes based on a dynamic Bayesian network (DBN). Groden and Collette [
15] combined a probabilistic fatigue deterioration model with a permanent set model by BN to update the reliability of marine structures. Hackl and Kohler [
16] proposed a model that combines structural reliability analysis and BN to estimate the reliability of reinforced concrete structures, and this model can be updated based on the outcomes of inspection and monitoring. Another is Markov Chain Monte Carlo (MCMC), aiming at overcoming the inefficiency of Monte Carlo simulation (MCS) in estimating small failure probabilities. Song et al. [
17] proposed an efficient reliability sensitivity algorithm based on the conditional samples generated by MCMC, and this algorithm is suitable for highly nonlinear limit state equations. Papaioannou et al. [
18] introduced a new method for MCMC sampling in standard normal space for subset simulation, and this method is simpler than previous algorithms while maintaining the same accuracy and efficiency.
Inspection and maintenance interventions are common practices to improve the safety and reliability of structural systems [
19,
20]. Maintenance actions depend on inspection outcomes [
11,
21]. If the inspection result is ‘detection’, an immediate repair action will be carried out. By contrast, if the inspection result is ‘no detection’, no action will be performed. Inspection and maintenance should be planned rationally to satisfy the needs of owners. The optimal inspection and maintenance planning is to find a set of rules at each decision node, aiming at optimizing some indicators (e.g., cost minimization).
Numerous algorithms for this optimization problem have been proposed in the past years. In the 1950s, the founding algorithms were created in the field of operations research [
22]. Later, Bellman [
23] developed Markov decision processes (MDPs), providing a strong tool in asset management. However, there are many limitations in MDPs. One is that it assumes inspections always reflect the true state of the system with certainty [
24]. To handle these limitations, partially observable Markov decision processes (POMDPs) were developed and became a more general algorithm. Faddoul et al. [
25] presented generalized POMDPs that combine decision analysis and dynamic programming and applied it to maintenance optimization. Schöbi and Chatzi [
26] described an enhanced variant of POMDPs for the maintenance planning of infrastructure, providing a framework to address non-linearities in action models. MDPs and POMDPs are suitable for simple systems consisting of a limited number of components. However, when the number of components is large, they become impractical because the Markovian transition matrices become extremely large. To this end, Andriotis and Papakonstantinou [
27] proposed a deep reinforcement learning (DRL) algorithm to establish effective life cycle strategies for large multi-component systems.
Regarding the objective function of this optimization problem, most studies focus on cost/ risk minimization. Garbatov and Guedes Soares [
28] presented a maintenance cost minimization method for floating structures and investigated the effect of different parameters on maintenance cost. Straub and Faber [
29] proposed an approach that considers the entire system in risk-based inspection (RBI) planning and discussed the dependencies among system components. Nielsen and Sørensen [
30] evaluated the cost of inspections, repairs, and lost production for a single wind turbine and assessed the impact of parameters such as failure rate and inspection interval on costs. Luque and Straub [
13] presented an integral risk-based optimization method for structural system inspections using the DBN framework, which allows for the rapid calculation of system reliability based on inspection results. Given that risk is the multiplication of cost and corresponding probability, risk optimization is also a kind of cost optimization. Thus, previous studies mainly focus on cost optimization.
For large offshore structures such as offshore platforms, offshore wind farms, etc., reliability performance may be more important than cost. However, a few studies considered reliability growth [
31], and these studies were applied at the component level. In engineering practice, inspection planning is generally aimed at structural systems. Therefore, it is necessary to develop a system-level reliability growth (SRG) model.
In this study, an SRG model is proposed for the first time and applied to efficient inspection of offshore wind farms. This study is novel in the following aspects: (a) SRG is defined as an increase in the percentage of system-level reliability index with and without inspection. This novel definition emphasizes the importance of system-level reliability performance and provides another way of assessing the value of inspections. (b) A framework to quantify the SRG of structural systems is proposed, which also works for a single component. (c) It is interesting to establish the law of diminishing SRG, i.e., the maximum SRG gained from one additional inspection decreases as the number of inspections increases. (d) The characteristics of expected posterior failure probability and the mechanism of optimal inspection time are revealed. In summary, theoretically, the proposed SRG model evaluates system reliability performance from a new perspective, promoting the development of system reliability theory. In practice, this model can be used to develop efficient inspection planning, in order to enhance the reliability and economic performance of offshore wind farms and contribute to the sustainable development of clean energy.
The rest of this study is arranged as follows.
Section 2 introduces the fatigue damage assessment model of the monopile OWT.
Section 3 presents the SRG model based on system-level reliability growth, including prior reliability analysis by S (i.e., stress range)-N (i.e., fatigue life) approach, posterior reliability analysis by Bayesian inference, and the quantification method of SRG.
Section 4 describes the optimization method to obtain the optimal inspection planning based on SRG.
Section 5 gives the results and discussion of the optimal inspection planning for offshore wind farms.
Section 6 shows some important conclusions.
2. Fatigue Damage Assessment Model
In this study, the offshore wind farm is a series system that consists of 15 OWTs (see
Figure 1). These OWTs are subjected to multiple loadings from wind, waves, and currents. The water depth is 25 m, which is relatively shallow, so the monopile foundation is adopted to support each offshore wind turbine. The monopile support structure is composed of a transition piece and a monopile. In the construction of this monopile support structure, steel plates are rolled and then horizontally and vertically welded. Accordingly, a single OWT, which is simplified as a monopile support structure, is defined as a component of the system.
For the
th (
) OWT, when the long-term stress range
is described by a Weibull distribution, and a one-slope S-N curve is used. The accumulated fatigue damage of the monopile support structure
is given by [
32] as follows:
where
is the total number of loading cycles during the design service life;
is the scale parameter of Weibull stress range distribution;
is the shape parameter of Weibull stress range distribution;
is the Gamma function;
is the inverse negative slope of S-N curve;
is the intercept of the S-N curve with the logarithmic N-axis. If the accumulated fatigue damage follows the Palmgren–Miner rule [
33,
34] and the annual number of loading cycles
is a constant over time
years, the accumulated fatigue damage of
th monopile support structure after
years
can be obtained as follows:
4. Optimal Inspection Planning Based on SRG Computation
Following
Section 3.3 for the first inspection, the optimal inspection time can be obtained by maximizing SRG. Then, when the second inspection is considered, the first inspection is conducted at the optimal inspection time, and the first inspection outcome is known as
or
. Based on this inspection outcome, the reliability index can be updated. It is called the posterior reliability index after the first update, which will become the prior reliability index for the second update. In other words, the purpose of the known inspection outcome for the first inspection is to obtain the prior reliability index for the second update. After that, given all the branches (
,
, and
) into consideration, the posterior reliability index can be obtained. Finally, SRG for the second inspection can be obtained, and the optimal second inspection time is when this SRG is maximized.
The procedure mentioned above can be repeated until the
th inspection. Generally, the optimal inspection planning is to solve an optimization problem by maximizing SRG for each inspection, and this problem can be described as follows:
where
is the number of inspections;
is the
th inspection time;
is the SRG corresponding to
;
is the inspection quality;
is the vector of basic parameters in the limit state function (see
Table 1);
is the inspection outcome. Note that
is the design variable,
is the objective function,
,
, and
are deterministic variables, and
is the vector of basic random variables. According to Equations (1)–(18), for given
,
,
, and
,
is a function of a single variable
. Thus, the exhaustive search algorithm is employed to solve this optimization problem, i.e., the
corresponding to each (and all) inspection time
is calculated separately, and the time corresponding to the maximum
is the optimal inspection time.
5. Results and Discussion
Consider two inspections and two different inspection qualities (
and
) are conducted for the offshore wind farm. The repairment effect is described by an ‘as good as new’ model [
42,
43]. It means that when the repair action is performed, the distribution of accumulated fatigue damage will recover to the initial condition. Three different inspection rules are defined as follows:
Inspection rule 1 (IR1). Each inspection is conducted for the same OWT (e.g., No. 1 OWT for two inspections). If the inspection outcome is ‘detection’, all the OWTs will be repaired;
Inspection rule 2 (IR2). Each inspection is conducted for the different OWTs (e.g., No. 1 OWT for the first inspection and No. 2 OWT for the second inspection). If the inspection outcome is ‘detection’, all the OWTs will be repaired;
Inspection rule 3 (IR3). Each inspection is conducted for the different OWTs (e.g., No. 1 OWT for the first inspection and No. 2 OWT for the second inspection). If the inspection outcome is ‘detection’, Only the inspected OWT will be repaired.
MCS has advantages in dealing with nonlinear limit state equations and correlations of random variables, and it has been widely used to calculate system reliability in previous studies [
13,
36,
38]. MCS with 100,000 samples is mainly used to obtain the following results. The maximum system failure probability
is in the order of
, and the coefficient of variation of MCS is in the order of
using Equation (8). Therefore, this sample size is appropriate.
This section is arranged as follows.
Section 5.1 compares the calculations of the reliability index between MCS and FORM to verify the validity of the results of MCS.
Section 5.2 compares the inspection efficiency of the offshore wind farm under different inspection strategies.
Section 5.3 compares the inspection efficiency between a single OWT and the offshore wind farm.
Section 5.4 conducts an in-depth analysis of the optimal inspection time.
Section 5.5 makes some suggestions for future research.
5.1. Verification of Calculation Results by MCS and FORM
In
Section 3.1, MCS and FORM are mentioned to calculate the prior system reliability index.
Figure 4 shows the results of these two methods. It can be seen that the difference between MCS and FORM is sufficiently small. Thus, MCS is a reliable method to obtain the results in the following subsections.
5.2. Comparative Study of SRG under Different Inspection Strategies
5.2.1. SRG under IR1
For the first inspection, the SRG of the offshore wind farm with respect to the inspection time for
and
is shown in
Figure 5. First of all, SRG increases to a maximum and then decreases with the inspection time so that the optimal inspection time can be determined: the 11th year for higher inspection quality and the 12th year for lower inspection quality. Then, the maximum SRG with higher inspection quality is greater than that with lower inspection quality. This is reasonable because higher inspection quality usually means higher inspection costs. If higher inspection quality does not gain more benefits in terms of SRG, the owner will not pay additional costs for higher inspection quality. Finally, it is checked that these two findings do not change for other inspection rules (IR2 and IR3) or the
th inspection (
).
For the second inspection, the initial inspection time is exactly the optimal first inspection time. As mentioned earlier, the first inspection outcome has two possible events: ‘detection’ or ‘no detection’. If the first inspection outcome is ‘detection’, according to the assumption of IR1, the system condition will recover to its original state. SRG over the second inspection time can be obtained by using the same method as the first inspection, and its curve is shown in
Figure 6. The optimal inspection times are the 16th year for higher inspection quality and the 17th year for lower inspection quality. On the other hand, if the first inspection outcome is ‘no detection’, SRG over the second inspection time is shown in
Figure 7. The optimal inspection times are also the 16th year for higher inspection quality and the 17th year for lower inspection quality.
5.2.2. SRG under IR2
Figure 8,
Figure 9 and
Figure 10 show the SRG under IR2. For the first inspection, SRG under IR2 (see
Figure 8) is the same as that under IR1 (see
Figure 5) because of the same inspection object (No. 1 OWT) and same repair action (i.e., if the inspection outcome is ‘detection’, all the OWTs will be repaired).
Consider the second inspection; if the first inspection outcome is ‘detection’, SRG under IR2 (see
Figure 9) is similar to that under IR1 (see
Figure 6) because the system condition recovers to its original state after the first inspection. This means that the system condition is the same under IR1 and IR2 before the second inspection is conducted. Although the inspection object is different (i.e., No.1 OWT under IR1 and No.2 OWT under IR2), the parameters of these two OWTs are identically distributed. As a result, the difference in SRG calculations under IR1 and IR2 is minimal.
However, if the second inspection outcome is ‘no detection’, SRG under IR2 (see
Figure 10) is much different from that under IR1 (see
Figure 7). First, SRG at the initial inspection time is larger than 0 under IR2, while it is equal to 0 under IR1. The reason is that at this time, the first inspection has been conducted on No.1 OWT. If the second inspection is also conducted on No. 1 OWT, no additional information will result in zero SRG. But if the second inspection is conducted on No. 2 OWT at this time, the system will gain new information, so SRG is larger than 0.
Second, comparing
Figure 7 and
Figure 10, the maximum SRG under IR2 (14 for
and 8.3 for
) is more than twice as much as that under IR1 (5.7 for
and 3.4 for
). One possible explanation is that under IR1, two inspections are conducted on the same object (No. 1 OWT), leading to higher information redundancy and lower SRG. This implies that changing the inspection object is a more efficient way to increase the system reliability for the second inspection. It can be expected that the difference in inspection efficiency between IR2 and IR1 grows as the number of inspections increases.
5.2.3. SRG under IR3
Figure 11,
Figure 12 and
Figure 13 show the SRG under IR3. Comparing
Figure 8 and
Figure 11, for the first inspection, the maximum SRG under IR2 is larger than that under IR3, and the difference is huge for the first inspection (e.g., 112 versus 10.8 for
). It is checked that this finding works for the
th (
) inspection. This indicates that the inspection efficiency of IR2 is higher than that of IR3. In other words, when inspecting a ‘detection’ outcome, overall repairing the system is a more efficient way to increase the system reliability than just repairing the inspected object. The reason is that the offshore wind farm is regarded as a series system. Under IR3, if the inspection outcome is ‘detection’, the reliability of the inspected OWT can be enhanced through repair, but the unrepaired OWTs are still the weak links of the system, resulting in a lower SRG.
In summary, from the perspective of gaining the maximum SRG, IR2 has the highest inspection efficiency compared to IR1 and IR3. If the inspection quality is higher, the inspection efficiency is higher. In other words, increasing the number of repair objects when inspecting a ‘detection’ outcome, changing the inspection object for each inspection, and increasing the inspection quality can improve inspection efficiency. This finding is important for structural system reliability engineering. It highlights the need for high-quality, high-coverage inspections and comprehensive repairments to improve the system reliability.
5.2.4. The Maximum SRG with Respect to the Number of Inspections
Table 2 shows the maximum SRG under IR1, IR2, and IR3. For a given inspection quality, the maximum SRG for the first inspection is larger than that for the second inspection. In other words, as the number of inspections increases, the maximum SRG gained from one additional inspection (which can be defined as marginal SRG) decreases. This is in line with the law of diminishing marginal utility in economics and can help find the optimal number of inspections. For the first inspection, the benefits (i.e., maximum SRG) are usually larger than the costs. As the number of inspections increases, the marginal benefits may be smaller than the marginal costs owing to the diminishing marginal benefits, resulting in unprofitability. By weighing the benefits against costs, the optimal number of inspections is the one at which the marginal benefits are equal to the marginal costs.
5.3. Differences between a Single OWT and the Offshore Wind Farm
When the number of components is equal to one, the SRG of the offshore wind farm under IR1 denotes the SRG of a single OWT under IR1. For comparison,
Table 3 presents the maximum SRG for the system (15 OWTs) and a single OWT under IR1. It can be seen that the law of diminishing marginal SRG still works at the component level. The difference is that this diminishing rate of the system is faster than that of a single OWT. One possible explanation is that the system is more difficult to maintain and manage due to complex component interactions and system integration issues, resulting in a sharp reduction in its potential for reliability growth.
5.4. In-Depth Analysis of the Optimal Inspection Time
5.4.1. Characteristics of Expected Posterior Failure Probability
Recall Equation (14): to maximize
,
needs to be maximized, or the expected posterior failure probability
needs to be minimized. Thus, the characteristic of
is discussed in this section. According to the event tree model in
Section 3.2.1,
at the end time of the service life (i.e., the 20th year)
can be expressed as follows:
Note that Equation (19) is a function of inspection time
(
). Consider the first term of the right side of this Equation for a given inspection time because
represents the system condition is ‘failure’ at this time,
, and the system must fail at the 20th-year condition on
, i.e.,
=1. Let
be decomposed into two components
and
, then this Equation can be rewritten as follows:
where
and
are defined as follows:
Clearly, Equation (21) indicates that the first component is equal to the prior failure probability, which reflects the prior physical model, and Equation (22) indicates that the second component is defined as the posterior failure probability considering inspection information. This is exactly the characteristic of the expected posterior failure probability: a probabilistic model that incorporates prior physical model and inspection data. It theoretically explains the nature of the expected posterior failure probability.
Figure 14 shows the trend of
and
under IR2 for the first inspection. As can be seen, the first component
increases as the inspection time increases. This is a negative effect on the expected posterior failure probability, which reflects the natural increase in the failure probability without inspection. However, the second component
decreases as the inspection time increases. This is a positive effect on the expected posterior failure probability, which reflects the reduction effect of inspection on the failure probability. Then,
is smaller than
at the early stage but larger at the later stage. This implies that
plays a dominant role in the early period,
in the later period. As a result, the sum of the two is likely to be smallest in the medium period.
5.4.2. Mechanism of the Optimal Inspection Time
The optimal inspection time is when SRG is maximized, or the expected posterior failure probability is minimized. Therefore, the first-order derivative of Equation (20) should be zero at that time, i.e.,
where
is the increasing rate of the prior failure probability
,
is the decreasing rate of the posterior failure probability considering inspection data
. Note that
and
are both positive, ‘−’ in front of
represents a decreasing trend of
. Accordingly, the optimal inspection time
can be expressed as follows:
Equation (24) reveals the mechanism of the optimal inspection time: it is the time at which the increasing rate of the prior failure probability is equal to the decreasing rate of the posterior failure probability considering inspection data.
To obtain the value of
and
, a central-difference formula is adopted as follows:
Figure 15 shows the trend of
and
under IR2 for the first inspection. There are two critical points,
P1 (horizontal coordinate 11.2) and
P2 (horizontal coordinate 11.7), which lie exactly at the optimal inspection times for
and
, respectively. 11.2 and 11.7 are rounded to be 11 and 12, respectively. This agrees with the results of
Section 5.2.2. Moreover,
shows a monotonically increasing trend, which agrees with the law of fatigue crack propagation: the growth rate of fatigue crack is very small at the early stage and then becomes faster and faster. In addition,
is larger than
at the early stage but smaller at the later stage. For this reason, the intersection of
and
is likely to be at the medium stage, i.e., the optimal inspection time is near the midpoint of the remaining service life.
Note that
Section 5.3 focuses on the first inspection under IR2 (the most efficient inspection rule). It is checked that these results hold for other inspection rules (IR1 and IR3) and the
th (
) inspection.
5.5. Suggestions for Future Research
While this study primarily concentrated on determining the optimal inspection time, it can be extended to consider other performance indicators. For example, researchers can optimize the number of OWTs to be inspected at once. This entails determining, within resource constraints, the ideal quantity of OWTs to maximize the performance and reliability of the entire wind farm. Such studies will contribute to more holistic inspection optimization strategies, enhancing the long-term sustainability and economic viability of offshore wind farms.