Next Article in Journal
Molecular Characterization of the Actin Gene and 5′ Flanking Sequence from Brown Macroalga Saccharina japonica (Laminariales, Phaeophyta)
Next Article in Special Issue
Structural Influence of the Cargo Holds of a 3000 m3 Wellboat on a Double-Bottom Floor
Previous Article in Journal
Ship-Radiated Noise Separation in Underwater Acoustic Environments Using a Deep Time-Domain Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Probabilistic Prediction of Floating Offshore Wind Turbine Platform Motions via Uncertainty Quantification and Information Integration

1
Department of Ocean Science and Engineering, Southern University of Science and Technology, Shenzhen 518055, China
2
School of Civil Engineering, Sun Yat-sen University, Zhuhai 519082, China
3
Department of Offshore Energy, Norwegian Geotechnical Institute, 0484 Oslo, Norway
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(6), 886; https://doi.org/10.3390/jmse12060886
Submission received: 29 April 2024 / Revised: 23 May 2024 / Accepted: 24 May 2024 / Published: 26 May 2024
(This article belongs to the Special Issue Advances in the Performance of Ships and Offshore Structures)

Abstract

:
The accurate prediction of short-term platform motions in a real environment is crucial for the safe design, operation, and maintenance of floating offshore wind turbines (FOWTs). Numerical simulations of motions are typically associated with high uncertainties due to abstracted theoretical models, empirical parameters, initial environment parameters, etc. Therefore, it is necessary to integrate other sources of information associated with less uncertainty, e.g., monitoring data, for accurate predictions. In this paper, we propose a probabilistic prediction based on the Bayesian approach that logically integrates motion monitoring data with simulated motion predictions of FOWTs, considering uncertainties in the environment model, structural properties, motion prediction method, monitoring data, etc. The approach consists of constructing a prior probability density function (PDF) of a random variable (which characterizes the largest value of the initial motion response) via numerical simulations and a likelihood function based on platform motion monitoring data and deriving a posterior PDF of the random variable by Bayesian updating. Then, posterior distributions of short-term extreme motion responses are derived using the posterior PDF of the random variable, representing lower uncertainty and improved accuracy. A Metropolis–Hastings algorithm is adopted to obtain PDFs of complex probability distributions. The effectiveness of the approach is demonstrated on a real FOWT platform in Scotland. The proposed probabilistic prediction approach results in posterior distributions of short-term extreme platform motions associated with less uncertainty and higher accuracy, which is attributed to integrating prior knowledge with monitoring data.

1. Introduction

In recent years, renewable energy has attracted significant attention due to the rising energy demand and growing concerns about environmental degradation. Among renewable energy sources, offshore wind energy holds substantial promise, as evidenced by its rapid growth in the market. Despite this expansion, numerous technical challenges exist within the offshore wind energy industry [1]. For example, considerable efforts have been devoted to understanding their structural response and motion behavior. However, existing studies predominantly focus on analyses under simplified conditions, and in-depth studies on dynamic responses under complex marine conditions (e.g., combining waves, wind, and ocean currents) and real environments, especially extreme environments, are lacking [2]. China is one of the world’s high-incidence countries for typhoons, with several typhoons along the coast each summer, causing significant impacts and losses to wind farms. A review published in Science in 2019 also highlighted that typhoon-induced structural damage to wind turbines poses a significant threat to the operation of wind farms [3]. Therefore, it is crucial to deeply understand and predict the complex dynamic responses of floating offshore wind turbines (FOWTs) under real environments accurately. Reducing uncertainties in simulating dynamic responses of FOWTs is key to accurate predictions.
Guedes Soares [4] classified uncertainty into intrinsic and epistemic types. Epistemic uncertainty arises from assumptions, simplifications, or incomplete information concerning a system or process, which can be reduced by collecting more information about physical quantities and improving measurement methods. Several studies have been conducted to assess the uncertainty in simulation models (e.g., abstracted theoretical models, empirical parameters) of FOWTs and revealed that a significant portion of the discrepancies between simulation results and measured data can be attributed to the uncertainty in simulation models [5,6,7,8,9]. Consequently, some research efforts have been made to optimize simulation models to achieve more accurate results for detecting structural risk probability [10,11,12,13].
On the other hand, structural health monitoring (SHM) serves as another powerful technique to collect real-world data about the platform’s motion responses and monitor structural risks. A variety of classical inference methods based on SHM have been developed to identify and predict structure responses [14,15,16,17,18]. However, structural health monitoring in complex ocean environments is costly and subject to observational noise, particularly in the early stages when monitoring data are scarce or unavailable [19,20]. Frangopol [21] proposed that incorporating prior information about a structural system can help augment limited monitoring data, leading to more reliable structural assessment. Classical statistical methods do not allow for effectively combining the prior knowledge, reflecting experiential judgments, with the likelihood of quantifying the chance of observations, whereas Bayesian methods are better suited for this task [22,23,24,25]. In Bayesian inference, a prior model is a probability distribution that represents the existing state of knowledge about a parameter or property before obtaining additional data, and a posterior model (obtained via Bayesian updating) denotes the probability distribution after obtaining data.
The Bayesian model updating has demonstrated significant computational capabilities in reducing the uncertainty of engineering models and provides valuable insights into model predictions based on observational data [26]. Garbatov and Guedes Soares [27] adopted a Bayesian approach to update parameters of probability distributions (which govern the reliability assessment of maintained floating structures) using inspection data, achieving improved reliability estimation. Okasha et al. [28] proposed a Bayesian updating approach to combine prior design code-based knowledge with SHM data for improved life-cycle performance assessment of ship structures under uncertainty. This approach utilized advanced modeling tools and techniques for lifetime reliability computations, including an optimization-based method for the ultimate failure moment, a progressive collapse method for the first failure moment, and hybrid Latin hypercube sampling. SHM data obtained from a scaled model testing of a joint high-speed sealift ship were used to update the structural life-cycle performance in the Bayesian framework. Okasha and Frangopol [29] extended the approach mentioned earlier into a life-cycle bridge management framework, where prior knowledge and SHM data about an existing bridge located in Wisconsin were integrated for better performance assessment. Zhu and Frangopol [30] further applied this approach to update wave-induced load effects and evaluate the structural performance of a ship. Their study focused on two parameters (location parameter and scale parameter) in the Type I extreme value distribution. It explored three general cases associated with parameter updating (i.e., one parameter, two parameters separately, and two correlated parameters simultaneously). Decò and Frangopol [31] proposed a risk-informed approach for ship hulls, which integrated SHM data via Bayesian updating and presented an approach for real-time optimal routing of ships. In the last five years, there has been growing attention to the use of Bayesian updating methods to reduce the uncertainties of characteristics of structures via response measurements. For example, Fang et al. [32] used displacement influence line data obtained from GPS measurements to update the finite element characteristics of the Tsing Ma Bridge, a long-span cable-stayed bridge in Hong Kong. The results showed that the updated finite element characteristics based on the measured dynamic properties are consistent with the structural characteristics of the bridge. Additionally, Bayesian updating methods have been used to probabilistically update the prior distribution of resistance model uncertainty related to specific combinations of modeling hypotheses in reinforced concrete systems, thereby improving the final model’s accuracy [33]. Kamariotis et al. [34] applied the classical Bayesian updating method to study the deteriorating bridge systems, ultimately achieving the purpose of quantifying the benefit of the availability of long-term SHM vibration data. In order to address the critical issue of premature failure of bridge expansion joints, Ni et al. [35] employed Bayesian updating to update multiple parameters in a linear regression model between bridge temperature and expansion joint displacement, thus considering the uncertainty in the monitoring data and quantifying the model error and the prediction uncertainty.
In the context of wind turbines, Xu and Guedes Soares [36] highlighted Bayesian inference as a robust method for studying extreme mooring tension. They demonstrated the advantages of the mixture gamma-generalized Pareto distribution model in estimating extreme mooring tension for all investigated cases. Ding et al. [37] integrated wind turbine gear physical models and available health condition data via a Bayesian method, predicting wind turbine gearbox fatigue life under instantaneously varying load conditions. In order to accurately predict the performance of unsteady aeroelastic forces in a blade accounting for the uncertainty in aeroelastic characteristics, Sarma et al. [38] proposed a Bayesian method in which experimental data were used to update the structural and rotational parameters of a system of downwind wind turbine blades. Valikhani et al. [39] presented a Bayesian method to update the parameters of a physics-based wind turbine drivetrain model and input loads using measurement data. Moynihan et al. [40] updated the effective stiffness of soil springs for a 6 MW offshore monopile wind turbine and compared the deterministic and probabilistic (Bayesian) methods. Bayesian model updating results successfully estimate the posterior distribution of updating model parameters with an increasing degree of certainty as more data are used. Wang et al. [41] extracted the natural frequencies and mode shapes of offshore jacket platforms from the structural dynamic response of a finite element model and the sensory measurements in the physical system to construct the likelihood function for the Bayesian updating method. The results indicated that the proposed model updating method maintained high effectiveness even with high-level measurement noise and model uncertainty.
The literature review reveals that the Bayesian updating method has shown great potential in reducing structural model uncertainty and could be applied to dynamically and probabilistically predict the extreme value of motion response of FOWTs under real environments. To the best of the authors’ knowledge, there is no such study to date, and a research gap exists. Also, it would be meaningful if we could maximize the value of limited available monitoring data to improve the accuracy of FOWT platform motion predictions. Platform motion monitoring is costly since FOWTs are mostly situated in remote, hash, and complex oceans.
The objective of this study is to bridge the research gap by developing a Bayesian updating approach that integrates prior knowledge (represented by numerical simulation results) with real monitoring data of FOWT motion responses. By probabilistically accounting for uncertainty in both sources of information, this approach leads to more accurate motion response predictions and an improved safety level of FOWTs. The flowchart (Figure 1) of the proposed approach involves constructing prior probability density functions (PDFs) and likelihood functions of model parameters, respectively, from simulation and monitoring data of platform motion responses and obtaining posterior distributions via Bayesian updating, based on which posterior predictions of short-term extreme platform motion responses are obtained and improved in terms of lower uncertainty. Specifically, numerical results of motion responses (performed with the effective knowledge and information of FOWTs) are used to construct the prior PDF of a random variable ( μ n ), which characterizes the largest value of motion response. The likelihood function is established by analyzing platform motion monitoring data, and the posterior probability distribution is derived via Bayesian updating. The posterior prediction results of short-term extreme platform motion responses are obtained based on the posterior probability distribution of μ n . A Markov chain Monte Carlo (MCMC) method coded with a Metropolis–Hastings algorithm is adopted to obtain PDFs of complex probability distributions. Compared to existing studies, the prior distribution model in this paper is more informative, since the distribution is obtained by making full use of all available knowledge and information about FOWTs while performing hydrodynamic simulations. The more informative prior distribution helps to improve the accuracy of posterior distribution with fewer monitoring data. The proposed approach is thus especially useful when monitoring data are limited. A case study on the Hywind wind farm in Scotland demonstrates the effectiveness of the proposed method.
The remainder of this paper is structured as follows: Section 2 presents the prior information and analysis of extreme platform motion responses. Section 3 and Section 4 describe the statistical analysis of platform motion monitoring data and the MCMC-based Bayesian updating approach. A case study of the Hywind wind farm in Scotland is presented in Section 5. Finally, conclusions are drawn in Section 6.

2. Prior Motion Response Analysis

To assess the safety of FOWTs, it is vital to accurately predict the platform motion responses in six degrees of freedom (DOFs). The six DOFs of a rigid floating structure consist of three translation DOFs (surge, sway, and heave) and three rotational DOFs (roll, pitch, and yaw). Physical monitoring and numerical simulation are two common methods for obtaining FOWT motion responses. Platform motion monitoring data are collected with sensors installed on the platform and the tower base. Numerical simulation is often used without monitoring data. A wide range of numerical techniques and software packages (e.g., FAST, HAWC2, AQWA, SESAM, etc.) are used for FOWT motion simulation [42]. Typically, a numerical model of a new design is developed and then validated by model-scale laboratory tests to assess the uncertainties associated with numerical models. In this study, the software package Ansys AQWA 2021 R1 with FAST v7 is utilized to simulate the motion responses of an FOWT subjected to wind, wave, and current loads [43,44]. Considering interactions between the FOWT and the mooring lines, the time-domain equation of motion of FOWTs can be written as follows [45]:
M + A X ¨ + 0 t k t τ X ˙ d τ + C X = F W A + F C U + F M R + F W I
where M is the inertia matrix; A is the added mass matrix at the infinite wave frequency; k is the delay function matrix; C is hydrostatic resilience; and X ¨ ,   X ˙ and X represent the acceleration vector, velocity vector, and displacement vector, respectively. F W A is the wave force, F W I and F C U are the wind load and current load, and F M R is the force passed from the mooring line to the FOWT.
An extreme value distribution of the platform motion responses is adopted as the choice for the prior distribution. Obtaining an accurate extreme value distribution model is a challenging task. One of three asymptotic extreme value distributions is selected to address this challenge based on the tail behavior of the probability distribution of the platform motion responses. Previous studies have indicated that the motion responses of floating structures typically follow the Rayleigh distribution and tend to the Gaussian distribution under severe sea conditions [46]. It is worth noting that exponential tails characterize both the Rayleigh distribution and the Gaussian distribution. Specifically, when the probability distribution of the platform motion responses exhibits an exponentially decaying tail, the Type I extreme value distribution, i.e., the Gumbel distribution, is assigned [47]. Therefore, the Gumbel distribution is usually chosen as an appropriate extreme value distribution for motion response given as follows:
f X x = α n e α n ( x μ n ) e x p [ e α n x μ n ]
where the variable of simulated motion responses is defined as the initial variable X . μ n is the characteristic largest value of the initial variable X , and α n is an inverse measure of the dispersion of the largest value of X . Therefore, the initial values of μ n and α n are related to the mean μ X and standard deviation δ X of X as follows:
α n = π 6 δ X
μ n 0 = μ X γ α n
where γ = 0.577216 is the Euler number. In comparison with the parameter α n , the parameter μ n has a larger effect on the extreme value. In this paper, α n is deterministic, and the initial value of μ n 0 is estimated based on prior knowledge and then updated with platform motion monitoring data in the proposed Bayesian framework. Based on previous studies, μ n follows a Lognormal distribution, which is given by Ang and Tang [47].
π μ n = 1 ζ μ n 2 π e [ 1 2 ln μ n λ ζ 2 ]
where ζ and λ are the hyperparameters of the variable μ n [48] and are related to the mean and standard deviation of the variable μ n according to Equations (6) and (7).
ζ = ln ( 1 + δ 2 μ n 0 2 )
λ = ln μ n 0 1 2 ζ 2
where the mean μ n 0 is obtained from prior knowledge, and δ is set to be the coefficient of variation (COV) of 10% [28].

3. Statistical Analysis of Platform Motion Monitoring Data

In this study, we use monitoring datasets of a real FOWT to construct the likelihood function. We define positive directions as north, east, and above sea level. To obtain the datasets of the extreme value of the six-DOF motion responses, an extreme value extraction algorithm is developed and applied to the signals of the six-DOF motion responses. The signals of six DOFs are represented by X 1 , X 2 , X 3 , X 4 , X 5 , and X 6 . The algorithm is as follows:
  • Let X k = [ x k 1 , x k 2 , , x k n ] be the recorded platform motion monitoring signals of one DOF without duplicate data and calculate the mean value motion of this DOF. (n is the number of non-duplicate platform motion monitoring data, k = 1 , 2 , 3 , 4 , 5 , 6 ).
  • Let the initial i = 1 , starting from x k 1 to x k n ; find the x k ( i + 1 ) that satisfies ( x k ( i + 1 ) x k i ) × ( x k ( i + 2 ) x k ( i + 1 ) ) < 0 .
  • If x k ( i + 1 ) > x k i , plug the absolute value of x k ( i + 1 ) into the P k = m a x ( x k ( i + 1 ) , x k k ( i + 2 ) , , x k ( n 1 ) ) ( P means positive direction) array.
  • If x k ( i + 1 ) < x k i , plug the absolute value of x k ( i + 1 ) into the N k = m a x ( x k ( i + 1 ) , x k k ( i + 2 ) , , x k ( n 1 ) ) ( N means negative direction) array.
  • Calculate the mean value of the cell in the P k array and calculate the mean value of the N k array.
  • Compare the absolute value of these data. The array that has a larger mean value would be implemented in the Bayesian updating.

4. MCMC-Based Bayesian Updating

4.1. Bayesian Framework

The Bayesian theory is described by Gelman et al. [49] and Kruschke [50], which is given as follows:
π θ x = h x , θ m x = p x θ π ( θ ) m x
where θ is the vector of unknown parameters; x are the measured data; π ( θ ) is the prior distribution of the model parameters, and π θ x is the posterior PDF; p x θ is the likelihood function.
Referring to the above equation, a distinctive feature of the Bayesian theory is the capability of combining prior knowledge with the likelihood based on SHM data to obtain the posterior probability distribution of θ . In this paper, the posterior PDF of μ n in the Gumbel distribution of extreme motion responses are updated. The likelihood function is constructed based on the discrete peak samples extracted from platform motion monitoring data, which is given as follows:
p x θ = i = 1 N f X x θ = i = 1 N α n e α n ( x i μ n ) e x p [ e α n x i μ n ]
where f X x θ is the PDF of extreme platform motion responses; X is evaluated based on the platform motion monitoring data value x i , given that the parameter of the PDF is θ . m ( x ) is a normalizing constant required to make π θ x a proper PDF and is calculated independently of θ using Equation (10).
m ( x ) = Θ p x θ π ( θ ) d θ
The new PDF of X is obtained based on the original PDF of X and the posterior PDF of its parameter θ . The updated PDF of X could be described with Equation (11) using the total probability theorem.
f X ( x ) = Θ f X ( x θ ) π θ x d θ
Closed-form solutions for Equation (11) are not always possible. The MCMC method is used, and the algorithm is as follows:
  • Use the MCMC method to sample from the posterior distribution π θ x , obtaining a set of parameter samples θ .
  • For each parameter sample θ i , calculate f X ( x θ i ) .
  • Take the average of f X ( x θ i ) to approximate f X ( x ) .

4.2. Metropolis–Hasting Algorithm

In Bayesian statistical analysis, the evaluation of the PDFs of complex posterior distributions can be challenging. To address this, the MCMC method is utilized. In this method, a Markov chain is constructed with a stationary distribution that is equal to the target sampling distribution, using the chain’s states to generate random samples after an initial burn-in period in which the state distribution converges to the target. Developments on Bayes estimation using MCMC are summarized in detail by Li [51]. Among the multiple available MCMC algorithms, here, we employ the Metropolis–Hasting algorithm. The Metropolis–Hasting algorithm was first introduced by Metropolis et al. [52] and improved by Hastings [53]. It has been extensively employed in Bayesian data analysis. The Metropolis–Hasting algorithm generates samples from a target distribution by simulating an easy Markov chain that has the target density as its stationary distribution [54]. The MATLAB R2022b software package was used for implementation in this study and the Metropolis–Hasting algorithm is as follows:
  • Assume an initial value X ( t ) .
  • Draw a sample y ( t ) from a proposal distribution q ( y | X ( t ) ) .
  • Accept y ( t ) as the next sample X ( t + 1 ) with probability r ( X t , y ( t ) ) and keep X ( t ) as the next sample X ( t + 1 ) with probability 1 r ( X t , y ( t ) ) , where
    r ( x , y ) = m i n f ( y ) f ( x ) q ( x / y ) q ( y / x ) , 1 .
    f ( x ) is the target distribution, also known as the probability density function that we are trying to sample from.
  • Increment t t + 1 and repeat steps 2 and 3 until the desired number of samples is met.
The proposed distribution is crucial for producing high-quality samples effectively using this algorithm. There is currently no consensus in academia for the selection of a burn-in period for Markov chains. In this study, the burn-in period was selected as the initial 10% of samples based on the convergence of the Markov chain shown in the time series plot. To assess convergence, the Markov chain time series plot of three chains was examined, with 200,000 samples generated per chain. Convergence was determined by observing small fluctuations above and below a horizontal line without any trend or period. Additionally, the mean and variance of the samples for each chain were compared. Meanwhile, a sequential updating approach combined with the Metropolis–Hasting algorithm was adopted in this study. In this approach, a posterior distribution is obtained by multiple sequential updates of another update, and 50 SHM samples as one part are used for updating [49].

5. Case Study: Hywind Scotland

5.1. Brief Introduction of the Hywind Wind Farm

The world’s first floating wind farm, the 30 MW Hywind Scotland pilot park, has been in operation since 2017, demonstrating the feasibility of floating wind farms that could be much larger. The platform motion monitoring data of FOWT HS4 were employed in this study. Figure 2 shows the location of the FOWT HS4.

5.2. The Dimensions of Hywind and Some Considerations

Each Hywind unit consists of a turbine SWT-6.0-154 provided by Siemens Gamesa, a tower, and a floater. Referring to Chen et al. [55] and Equinor [56], the main dimensions of the Hywind Spar floater and tower, are shown in Figure 3. In addition, Figure 4 shows the mooring system. The properties of the mooring system are given in Table 1. The mooring system in the initial configuration should provide a vertical load of about 2660 KN. The tension at the upper end of the delta line (bridle) is about 590 KN, while the tension at the upper end of the main line is about 760 KN. Table 2 summarizes the properties of the FOWT.
The turbine of the Hywind is Siemens Gamesa SWT-6.0-154, which is similar to the 5 MW Hywind turbine and tower in OC3 [57]. Relevant information to specify the structural parameters is not available. The parameters of the 5 MW Hywind turbine and tower in OC3 are used to replace the wind turbine of HS4 in this study [57]. The mooring system is also a vital structure that connects a floating foundation to the seabed. The floating foundation is positioned by the restoring force and moment provided by the mooring lines. To enhance the stability of the Hywind floating structure, a delta-connection mooring chain design proposed by Koo et al. [58] was used, comprising three long anchor chains (120 degrees to each other) and six short anchor chains. This design effectively restricts the motion of a floater in the yaw degree of freedom. In this study, the delta-connection of mooring lines was considered to be directly linked to the hull with additional yaw stiffness [55]. These simplifications and alternatives also reflect some real situations in engineering.

5.3. Source of Platform Motion Monitoring Data

The platform motion monitoring data of 11 cases of Hywind in Scotland were provided by Equinor and ORE Catapult [56]. The operation condition (Case 7) and idle condition (Case 10) of the FOWT were chosen to be the representative cases for the FOWT. Wind direction is defined as the direction from which the wind is coming from. Thus, a “northerly wind” is coming from the north, and the wind direction is defined as 0 degrees, while an “easterly wind” has a wind direction of 90 degrees. Wave direction is similarly defined as the direction from which the waves are coming from. Thus, “northerly waves” are coming from the north, and the wave direction is defined as 0 degrees, etc. In contrast, the current direction is defined as the direction toward which the current is going. Thus, a “northerly current” is coming from the south, and the current direction is defined as 0 degrees, while an “easterly current” is coming from the west toward the east, and the current direction is defined as 90 degrees. The wind speed is the undisturbed wind speed measured with the turbine’s anemometer located at elevation +99 relative to sea level [56]. Details of the environmental conditions of Cases 7 and 10 are given in Table 3 and Figure 5.
The displacement of surge and sway was measured by GPS located at elevation +15.3 relative to sea level. The north motion is the positive direction for sway, and the east motion is defined as the positive motion for the surge. It is worth noting that since the GPS is located at the platform level, some of the translational motions may be due to rotations (roll/pitch motions) [56]. The platform motion angles at roll and pitch DOFs were evaluated based on the angles of the tower’s motion measured with sensors. The coordinate system of the tower motion measurements is presented in Figure 6; its origin is set at +16.9 m relative to sea level. This paper concerns the measured data at four DOFs (i.e., surge, sway, roll, and pitch). In this study, any noise treatments to the monitoring data were not applied before updating, because the quality of the monitoring data used was considered acceptable by using fast Fourier transform analysis. Also, filtering techniques for mitigating noise probably introduce subjective judgment biases and increase uncertainty. The Bayesian updating method has an inherent advantage in handling uncertainties, making it reducible and common to use monitoring signals with noise. The posterior distribution can quantify the relative plausibility of different values of the model parameters based on the available incomplete information, such as signals with noise [59].

5.4. Prior Analysis by Numerical Simulation

An aero-hydro-servo-elastic coupling framework for dynamic analysis of FOWTs based on AQWA and FAST v7 (F2A), developed by Yang [44], was utilized to simulate the motion responses of the FOWT subjected to wind, wave, and current loads in an operation condition (Case 7) and an idle condition (Case 10). In the idle condition, the blades were locked at a 90° pitch angle and aligned with the wind direction, while the nacelle yaw was maintained at 0 degrees to mitigate wind loads and motions.
The Blade element momentum (BEM) theory and the potential theory were used for calculating aerodynamics and hydrodynamics, respectively. The Kaimal wind spectrum was adopted in FAST v7. The JONSWAP wave spectrum and current load parameters were defined in the Hydrodynamic Response module of the AQWA Workbench 2021 R1. The full low-frequency quadratic transfer function (QTF) was applied due to the importance of the low-frequency wave loads for the spar buoy. Equation (1) is used to calculate the motion responses of FOWT. Each simulation lasted 3600 s with a time step of 0.005 s. Ten different time-domain realizations were obtained for each condition (of significant wave height, peak-spectral period, wind speed, etc.) by changing random numbers of wave and wind between simulations.

5.4.1. Natural Periods

It is crucial to verify the natural periods of the FOWT and tune damping by free decay tests. In this study, we found that the simulation results and the measurement results from Equinor [56] are in good agreement, with a discrepancy of approximately 5% or less, as shown in Table 4.

5.4.2. Motion Response Results

Figure 7 presents the time series and spectral analysis of surge and roll motion responses for both simulation data and monitoring data of Case 7. Figure 8 presents the time series and spectral analysis of surge and roll motion responses for both simulation data and monitoring data of Case 10. Figure 7b,d and Figure 8b,d show dominant frequencies of motion responses, but low-frequency motion responses are underestimated even with full QTF. Some differences in the mean value of motion response can be observed between simulation data and monitoring data in Figure 7a,c and Figure 8a,c. The reasons can be attributed to three main aspects. Firstly, alternative parameters of the wind turbine and tower in OC3 were used. Secondly, it is challenging to obtain accurate values of drag coefficients. Thirdly, the real ocean conditions may not always follow the wind and wave spectrums we used. The existing differences highlight the value of our study, which aims to reduce the differences and improve predictions by integrating limited monitoring information with numerical simulations via Bayesian updating. At present, the simulated motion responses represent the best of current knowledge about FOWTs and should be utilized to develop a prior model. Table 5 shows the mean and standard deviation values of extreme motion responses for surge and roll in Cases 7 and 10. Table 6 provides the parameters of the corresponding Gumbel and Lognormal distributions for the surge and roll motion responses in Cases 7 and 10. Figure 9 and Figure 10 show the statistical histogram of extreme motion response and fitted Gumbel PDF in Cases 7 and 10.

5.5. Platform Motion Monitoring and Bayesian Updating

The prior distribution model about the extreme motion responses was obtained using AQWA 2021 R1and FAST v7. For the likelihood function, the parameter α n of the Gumbel distribution was derived by calculating the standard deviation δ X of motion responses at each DOF and using Equation (3), while the parameter μ n (as a significant uncertain parameter) was updated by Bayesian updating. Using the mean μ n 0 from the prior analysis, with a COV of 10%, the hyperparameters of the prior PDF of the parameter μ n was calculated using Equations (6) and (7). The platform motion monitoring data were used to update the prior PDF of the parameter μ n via the Metropolis–Hasting algorithm.
Approximately 250 platform motion monitoring samples were applied sequentially in an increment of 50 samples in the Bayesian updating based on the rationale discussed in Section 4.2. The updated PDF of the parameter μ n in the Gumbel distribution at surge and roll DOFs in Cases 7 and 10 are shown in Figure 11 and Figure 12, while Figure 13 and Figure 14 provide a comparison before and after the update of μ n for four DOFs in Cases 7 and 10. The COV of the Gumbel distribution parameter μ n diminishes in both Cases 7 and 10, which means the uncertainties associated with μ n are reduced by incorporating the platform motion monitoring data. As monitoring and updating continue, it is reasonable that uncertainties decrease until they become acceptable, at which point the updating process could be terminated. The decision to stop the updating process was made based on convergence analysis (of the posterior distribution), according to the desired level of reduction in the uncertainty of a specific engineering problem. Meanwhile, Laskey [60] reported that the posterior distribution would be concentrated around the maximum likelihood estimate and relatively insensitive to the prior distribution when the sample size is very large. Finally, although the prior and posterior PDFs of the Gumbel distribution parameter μ n only reflect the reduction in the uncertainty of μ n by Bayesian updating based on the platform motion monitoring data, it is clear that the updated posterior PDFs of the parameter μ n have a favorable impact on the underlying Gumbel distribution of platform motion responses.
Figure 15 and Figure 16 show the original and updated PDFs of the extreme motion responses of surge and roll DOFs in Cases 7 and 10, along with frequency histograms of platform motion monitoring data. It can be observed that the updated distributions of extreme motion responses satisfy the Lognormal distribution and shift toward the platform motion monitoring data. The reduction in the standard deviation in Table 7 is a conclusion to show the reduction in uncertainty. It is clearly observed that the uncertainty of μ n is significantly reduced in most cases after Bayesian updating, with a maximum reduction of up to 91%. Although the uncertainty of variable μ n increases to some extent for roll motion in Case 7, the updated posterior mean value converges toward the monitoring data, outperforming the prior simulation results. The uncertainty of the extreme motion responses X is also slightly reduced due to the favorable impact of μ n . As more data were used, the Bayesian updating method could successfully estimate the posterior distribution of variables with an increasing degree of certainty. This paper presents the framework with two representative load cases, as the primary objective was to develop a Bayesian updating and inference framework leading to improved FOWT motion response probabilistic predictions. The framework can be applied to other load cases or sea states as well since load cases or sea states are just part of the inputs in this framework.

6. Conclusions

In this paper, we introduced an effective Bayesian updating method to dynamically and probabilistically predict the extreme value of motion response of FOWTs by using monitoring data under real sea conditions for the first time. The method provides a theoretical framework for rational and logical integration of monitoring data with prior motion predictions, enhancing the precision and reliability of these predictions under both normal and extreme conditions. The numerical results of motion responses (performed with the effective knowledge and information of FOWTs) were used to construct the prior PDF of a random variable ( μ n ), which characterized the largest value of motion response. The likelihood function was established by analyzing platform motion monitoring data, and the posterior probability distribution was derived via Bayesian updating. The posterior prediction results of short-term extreme platform motion responses were obtained based on the posterior probability distribution of μ n . Finally, accurate posterior predictions play an important role in platform structural design, risk management, and saving costs.
Our proposed approach was demonstrated to be effective in reducing uncertainty in the extreme response prediction of FOWT platforms using the Hywind wind farm in Scotland as a case study. The platform motion response data obtained from Equinor were used. It was clearly observed that the uncertainty of the variables μ n was significantly reduced in most cases after Bayesian updating, with a maximum reduction of up to 91%. Although the uncertainty of variable μ n increased to some extent for roll motion in Case 7, the updated posterior mean value converged toward the monitoring data, outperforming the prior simulation results. The uncertainty of the extreme motion responses X was also slightly reduced due to the favorable impact of μ n . As more data were used, the Bayesian updating method could successfully estimate the posterior distribution of variables with an increasing degree of certainty.

Author Contributions

Conceptualization, G.Z. and N.L.; methodology, G.Z. and N.L.; software, N.L.; validation, N.L. and Y.F.; investigation, N.L., G.Z. and Y.F.; resources, G.Z.; data curation, N.L.; writing—original draft preparation, N.L.; writing—review and editing, G.Z., Y.F., L.A. and N.L.; visualization, N.L.; supervision, G.Z. and Y.F.; funding acquisition, G.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Southern University of Science and Technology, grant number Y01316134.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data utilized in this study were obtained from publicly available sources and previously published datasets, which are appropriately cited in the references. No new data were created for this research. Due to privacy or ethical restrictions, certain data used in this study cannot be shared.

Acknowledgments

Thanks to Equinor for sharing operational data of Hywind Scotland through ORE Catapult.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Esteban, M.D.; Diez, J.J.; López, J.S.; Negro, V. Why offshore wind energy? Renew. Energy 2011, 36, 444–450. [Google Scholar] [CrossRef]
  2. Sun, K.; Xu, Z.; Li, S.; Jin, J.; Wang, P.; Yue, M.; Li, C. Dynamic response analysis of floating wind turbine platform in local fatigue of mooring. Renew. Energy 2023, 204, 733–749. [Google Scholar] [CrossRef]
  3. Veers, P.; Dykes, K.; Lantz, E.; Barth, S.; Bottasso, C.L.; Carlson, O.; Clifton, A.; Green, J.; Green, P.; Holttinen, H. Grand challenges in the science of wind energy. Science 2019, 366, eaau2027. [Google Scholar] [CrossRef] [PubMed]
  4. Guedes Soares, C. Quantification of Model Uncertainty in Structural Reliability; Springer: Berlin/Heidelberg, Germany, 1997; pp. 17–37. [Google Scholar]
  5. Robertson, A.N. Uncertainty analysis of OC5-DeepCwind floating semisubmersible offshore wind test campaign. In Proceedings of the 27th International Ocean and Polar Engineering Conference, San Francisco, CA, USA, 25–30 June 2017. [Google Scholar]
  6. Simos, A.N.; Ruggeri, F.; Watai, R.A.; Souto-Iglesias, A.; Lopez-Pavon, C. Slow-drift of a floating wind turbine: An assessment of frequency-domain methods based on model tests. Renew. Energy 2018, 116, 133–154. [Google Scholar] [CrossRef]
  7. Yang, P.; Li, J.; Zhang, W.; Wu, D.; Gu, X.; Ma, Q. Analysis on statistical uncertainties of wave loads and structural fatigue reliability for a semi-submersible platform. Ocean Eng. 2021, 237, 109609. [Google Scholar] [CrossRef]
  8. Popko, W.; Robertson, A.; Jonkman, J.; Wendt, F.; Thomas, P.; Müller, K.; Kretschmer, M.; Ruud Hagen, T.; Galinos, C.; Le Dreff, J.-B. Validation of numerical models of the offshore wind turbine from the alpha ventus wind farm against full-scale measurements within OC5 phase III. In Proceedings of the International Conference on Offshore Mechanics and Arctic Engineering, Glasgow, UK, 9–14 June 2019. [Google Scholar]
  9. Robertson, A.N.; Gueydon, S.; Bachynski, E.; Wang, L.; Jonkman, J.; Alarcon, D.; Amet, E.; Beardsell, A.; Bonnet, P.; Boudet, B. OC6 Phase I: Investigating the underprediction of low-frequency hydrodynamic loads and responses of a floating wind turbine. J. Phys. Conf. Ser. 2020, 1618, 032033. [Google Scholar] [CrossRef]
  10. Wang, L.; Robertson, A.; Jonkman, J.; Yu, Y.-H. OC6 phase I: Improvements to the OpenFAST predictions of nonlinear, low-frequency responses of a floating offshore wind turbine platform. Renew. Energy 2022, 187, 282–301. [Google Scholar] [CrossRef]
  11. Li, H.; Bachynski-Polić, E.E. Analysis of difference-frequency wave loads and quadratic transfer functions on a restrained semi-submersible floating wind turbine. Ocean Eng. 2021, 232, 109165. [Google Scholar] [CrossRef]
  12. Wang, L.; Robertson, A.; Jonkman, J.; Yu, Y.-H. Uncertainty assessment of CFD investigation of the nonlinear difference-frequency wave loads on a semisubmersible FOWT platform. Sustainability 2020, 13, 64. [Google Scholar] [CrossRef]
  13. Pegalajar-Jurado, A.; Bredmose, H. Reproduction of slow-drift motions of a floating wind turbine using second-order hydrodynamics and Operational Modal Analysis. Mar. Struct. 2019, 66, 178–196. [Google Scholar] [CrossRef]
  14. Mondoro, A.; Soliman, M.; Frangopol, D.M. Prediction of structural response of naval vessels based on available structural health monitoring data. Ocean Eng. 2016, 125, 295–307. [Google Scholar] [CrossRef]
  15. Martinez-Luengo, M.; Shafiee, M.; Kolios, A. Data management for structural integrity assessment of offshore wind turbine support structures: Data cleansing and missing data imputation. Ocean Eng. 2019, 173, 867–883. [Google Scholar] [CrossRef]
  16. Frangopol, D.M.; Strauss, A.; Kim, S. Bridge reliability assessment based on monitoring. J. Bridge Eng. 2008, 13, 258–270. [Google Scholar] [CrossRef]
  17. Strauss, A.; Frangopol, D.M.; Kim, S. Statistical, probabilistic and decision analysis aspects related to the efficient use of structural monitoring systems. Beton-und Stahlbetonbau 2008, 103, 23–28. [Google Scholar] [CrossRef]
  18. Strauss, A.; Frangopol, D.M.; Kim, S. Use of monitoring extreme data for the performance prediction of structures: Bayesian updating. Eng. Struct. 2008, 30, 3654–3666. [Google Scholar] [CrossRef]
  19. Hageman, R.; van der Meulen, F.; Rouhan, A.; Kaminski, M. Quantifying uncertainties for Risk-Based Inspection planning using in-service Hull Structure Monitoring of FPSO hulls. Mar. Struct. 2022, 81, 103100. [Google Scholar] [CrossRef]
  20. Katsoudas, A.S.; Silionis, N.E.; Anyfantis, K.N. Structural health monitoring for corrosion induced thickness loss in marine plates subjected to random loads. Ocean Eng. 2023, 273, 114037. [Google Scholar] [CrossRef]
  21. Frangopol, D.M. Life-cycle performance, management, and optimisation of structural systems under uncertainty: Accomplishments and challenges 1. Struct. Infrastruct. Eng. 2011, 7, 389–413. [Google Scholar] [CrossRef]
  22. Yuen, K.-V.; Kuok, S.-C. Bayesian methods for updating dynamic models. Appl. Mech. Rev. 2011, 64, 010802. [Google Scholar] [CrossRef]
  23. Yuen, K.-V. Bayesian Methods for Structural Dynamics and Civil Engineering; John Wiley & Sons: Hoboken, NJ, USA, 2010. [Google Scholar]
  24. Simoen, E.; Papadimitriou, C.; Lombaert, G. On prediction error correlation in Bayesian model updating. J. Sound Vib. 2013, 332, 4136–4152. [Google Scholar] [CrossRef]
  25. Beck, J.L.; Katafygiotis, L.S. Updating models and their uncertainties. I: Bayesian statistical framework. J. Eng. Mech. 1998, 124, 455–461. [Google Scholar]
  26. Xiao, X.; Li, Q.; Wang, Z. A novel adaptive importance sampling algorithm for Bayesian model updating. Struct. Saf. 2022, 97, 102230. [Google Scholar] [CrossRef]
  27. Garbatov, Y.; Guedes Soares, C. Bayesian updating in the reliability assessment of maintained floating structures. J. Offshore Mech. Arct. Eng. 2002, 124, 139–145. [Google Scholar] [CrossRef]
  28. Okasha, N.M.; Frangopol, D.M.; Decò, A. Integration of structural health monitoring in life-cycle performance assessment of ship structures under uncertainty. Mar. Struct. 2010, 23, 303–321. [Google Scholar] [CrossRef]
  29. Okasha, N.M.; Frangopol, D.M. Integration of structural health monitoring in a system performance based life-cycle bridge management framework. Struct. Infrastruct. Eng. 2012, 8, 999–1016. [Google Scholar]
  30. Zhu, B.; Frangopol, D.M. Reliability assessment of ship structures using Bayesian updating. Eng. Struct. 2013, 56, 1836–1847. [Google Scholar] [CrossRef]
  31. Decò, A.; Frangopol, D.M. Real-time risk of ship structures integrating structural health monitoring data: Application to multi-objective optimal ship routing. Ocean Eng. 2015, 96, 312–329. [Google Scholar] [CrossRef]
  32. Fang, C.; Liu, H.-J.; Lam, H.-F.; Adeagbo, M.O.; Peng, H.-Y. Practical model updating of the Ting Kau Bridge through the MCMC-based Bayesian algorithm utilizing measured modal parameters. Eng. Struct. 2022, 254, 113839. [Google Scholar] [CrossRef]
  33. Castaldo, P.; Gino, D.; Bertagnoli, G.; Mancini, G. Resistance model uncertainty in non-linear finite element analyses of cyclically loaded reinforced concrete systems. Eng. Struct. 2020, 211, 110496. [Google Scholar] [CrossRef]
  34. Kamariotis, A.; Chatzi, E.; Straub, D. Value of information from vibration-based structural health monitoring extracted via Bayesian model updating. Mech. Syst. Signal Process. 2022, 166, 108465. [Google Scholar] [CrossRef]
  35. Ni, Y.; Wang, Y.; Zhang, C. A Bayesian approach for condition assessment and damage alarm of bridge expansion joints using long-term structural health monitoring data. Eng. Struct. 2020, 212, 110520. [Google Scholar] [CrossRef]
  36. Xu, S.; Guedes Soares, C. Bayesian analysis of short term extreme mooring tension for a point absorber with mixture of gamma and generalised pareto distributions. Appl. Ocean Res. 2021, 110, 102556. [Google Scholar] [CrossRef]
  37. Ding, F.; Tian, Z.; Zhao, F.; Xu, H. An integrated approach for wind turbine gearbox fatigue life prediction considering instantaneously varying load conditions. Renew. Energy 2018, 129, 260–270. [Google Scholar] [CrossRef]
  38. Sarma, R.; Dwight, R.P.; Viré, A. Aeroelastic validation and Bayesian updating of a downwind wind turbine. Wind Energy 2020, 23, 864–883. [Google Scholar] [CrossRef]
  39. Valikhani, M.; Jahangiri, V.; Ebrahimian, H.; Moaveni, B.; Liberatore, S.; Hines, E. Inverse modeling of wind turbine drivetrain from numerical data using Bayesian inference. Renew. Sustain. Energy Rev. 2023, 171, 113007. [Google Scholar] [CrossRef]
  40. Moynihan, B.; Mehrjoo, A.; Moaveni, B.; McAdam, R.; Rüdinger, F.; Hines, E. System identification and finite element model updating of a 6 MW offshore wind turbine using vibrational response measurements. Renew. Energy 2023, 219, 119430. [Google Scholar] [CrossRef]
  41. Wang, M.; Leng, J.; Feng, S.; Li, Z.; Incecik, A. Precisely modeling offshore jacket structures considering model parameters uncertainty using Bayesian updating. Ocean Eng. 2022, 258, 111410. [Google Scholar] [CrossRef]
  42. Subbulakshmi, A.; Verma, M.; Keerthana, M.; Sasmal, S.; Harikrishna, P.; Kapuria, S. Recent advances in experimental and numerical methods for dynamic analysis of floating offshore wind turbines—An integrated review. Renew. Sustain. Energy Rev. 2022, 164, 112525. [Google Scholar] [CrossRef]
  43. Yang, Y.; Bashir, M.; Li, C.; Wang, J. Investigation on mooring breakage effects of a 5 MW barge-type floating offshore wind turbine using F2A. Ocean Eng. 2021, 233, 108887. [Google Scholar] [CrossRef]
  44. Yang, Y.; Bashir, M.; Michailides, C.; Li, C.; Wang, J. Development and application of an aero-hydro-servo-elastic coupling framework for analysis of floating offshore wind turbines. Renew. Energy 2020, 161, 606–625. [Google Scholar] [CrossRef]
  45. Chen, D.; Gao, P.; Huang, S.; Fan, K.; Zhuang, N.; Liao, Y. Dynamic response and mooring optimization of spar-type substructure under combined action of wind, wave, and current. J. Renew. Sustain. Energy 2017, 9, 063307. [Google Scholar] [CrossRef]
  46. Song, X.; Wang, S.; Hu, Z.; Li, H. A hybrid Rayleigh and Weibull distribution model for the short-term motion response prediction of moored floating structures. Ocean Eng. 2019, 182, 126–136. [Google Scholar] [CrossRef]
  47. Ang, A.H.-S.; Tang, W.H. Probability Concepts in Engineering: Emphasis on Applications to Civil and Environmental Engineering, 2nd ed.; John Willey & Sons: Hoboken, NJ, USA, 2007; pp. 172–179. [Google Scholar]
  48. Fink, D. A Compendium of Conjugate Priors. Progress Report: Extension and Enhancement of Methods for Setting Data Quality Objectives. DOE Contract 95. 1995. Available online: https://www.researchgate.net/publication/238622435_A_Compendium_of_Conjugate_Priors (accessed on 23 May 2024).
  49. Gelman, A.; Carlin, J.B.; Stern, H.S.; Dunson, D.B.; Vehtari, A.; Rubin, D.B. Bayesian Data Analysis; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  50. Kruschke, J. Doing Bayesian Data Analysis: A Tutorial with R, JAGS, and Stan; Academic Press: Cambridge, MA, USA, 2014. [Google Scholar]
  51. Li, Z. A review of Bayesian posterior distribution based on MCMC methods. In Proceedings of the International Conference on Computing and Data Science, Virtual Event, 12–17 August 2021. [Google Scholar]
  52. Metropolis, N.; Rosenbluth, A.W.; Rosenbluth, M.N.; Teller, A.H.; Teller, E. Equation of state calculations by fast computing machines. J. Chem. Phys. 1953, 21, 1087–1092. [Google Scholar] [CrossRef]
  53. Hastings, W.K. Monte Carlo sampling methods using Markov chains and their applications. Biometrika 1970, 57, 97–109. [Google Scholar] [CrossRef]
  54. Ghosh, J.K.; Delampady, M.; Samanta, T. An Introduction to Bayesian Analysis: Theory and Methods; Springer: Berlin/Heidelberg, Germany, 2006; pp. 218–220. [Google Scholar]
  55. Chen, P.; Jia, C.; Ng, C.; Hu, Z. Application of SADA method on full-scale measurement data for dynamic responses prediction of Hywind floating wind turbines. Ocean Eng. 2021, 239, 109814. [Google Scholar] [CrossRef]
  56. Equinor and ORE Catapult Collaborating to Share Hywind Scotland Operational Data. Available online: https://www.equinor.com/news/archive/2019-11-28-hywind-scotland-data (accessed on 23 May 2024).
  57. Jonkman, J.; Butterfield, S.; Musial, W.; Scott, G. Definition of a 5-MW Reference wind Turbine for Offshore System Development; National Renewable Energy Lab. (NREL): Golden, CO, USA, 2009; pp. 3–13, 15–16. [Google Scholar]
  58. Koo, B.; Goupee, A.J.; Lambrakos, K.; Lim, H.-J. Model test data correlations with fully coupled hull/mooring analysis for a floating wind turbine on a semi-submersible platform. In Proceedings of the International Conference on Offshore Mechanics and Arctic Engineering, San Francisco, CA, USA, 8–13 June 2014. [Google Scholar]
  59. Huang, Y.; Shao, C.; Wu, B.; Beck, J.L.; Li, H. State-of-the-art review on Bayesian inference in structural system identification and damage assessment. Adv. Struct. Eng. 2019, 22, 1329–1351. [Google Scholar] [CrossRef]
  60. Laskey, K.B. Model uncertainty: Theory and practical implications. IEEE Trans. Syst. Man Cybern. Part A Syst. Hum. 1996, 26, 340–348. [Google Scholar] [CrossRef]
Figure 1. The proposed Bayesian updating framework.
Figure 1. The proposed Bayesian updating framework.
Jmse 12 00886 g001
Figure 2. The location of the floating offshore wind turbine (FOWT) HS4.
Figure 2. The location of the floating offshore wind turbine (FOWT) HS4.
Jmse 12 00886 g002
Figure 3. The main dimensions of the Hywind Spar floater and tower.
Figure 3. The main dimensions of the Hywind Spar floater and tower.
Jmse 12 00886 g003
Figure 4. The mooring system.
Figure 4. The mooring system.
Jmse 12 00886 g004
Figure 5. The direction of wave, wind, and current of Cases 7 and 10.
Figure 5. The direction of wave, wind, and current of Cases 7 and 10.
Jmse 12 00886 g005
Figure 6. The coordinate system of the tower motions.
Figure 6. The coordinate system of the tower motions.
Jmse 12 00886 g006
Figure 7. Comparison of surge and roll motion responses predicted by AQWA+FAST and full-scale monitoring in Case 7: (a,c) time series of motion responses; (b,d) spectral analysis results.
Figure 7. Comparison of surge and roll motion responses predicted by AQWA+FAST and full-scale monitoring in Case 7: (a,c) time series of motion responses; (b,d) spectral analysis results.
Jmse 12 00886 g007aJmse 12 00886 g007b
Figure 8. Comparison of sway and pitch motion responses predicted by AQWA+FAST and full-scale monitoring in Case 10: (a,c) time series of motion responses; (b,d) spectral analysis results.
Figure 8. Comparison of sway and pitch motion responses predicted by AQWA+FAST and full-scale monitoring in Case 10: (a,c) time series of motion responses; (b,d) spectral analysis results.
Jmse 12 00886 g008
Figure 9. Statistical histograms of extreme motion responses and fitted Gumbel probability density functions (PDFs) for (a) surge and (b) roll in Case 7.
Figure 9. Statistical histograms of extreme motion responses and fitted Gumbel probability density functions (PDFs) for (a) surge and (b) roll in Case 7.
Jmse 12 00886 g009
Figure 10. Statistical histograms of extreme motion responses and fitted Gumbel PDFs for (a) surge and (b) roll in Case 10.
Figure 10. Statistical histograms of extreme motion responses and fitted Gumbel PDFs for (a) surge and (b) roll in Case 10.
Jmse 12 00886 g010
Figure 11. Posterior histograms of the Gumbel distribution parameter μ n of (a) surge and (b) roll motions in Case 7.
Figure 11. Posterior histograms of the Gumbel distribution parameter μ n of (a) surge and (b) roll motions in Case 7.
Jmse 12 00886 g011
Figure 12. Posterior histograms of the Gumbel distribution parameter μ n of (a) surge and (b) roll motions in Case 10.
Figure 12. Posterior histograms of the Gumbel distribution parameter μ n of (a) surge and (b) roll motions in Case 10.
Jmse 12 00886 g012
Figure 13. Bayesian updating results of the Gumbel distribution parameter μ n for (a) surge and (b) roll motions in Case 7.
Figure 13. Bayesian updating results of the Gumbel distribution parameter μ n for (a) surge and (b) roll motions in Case 7.
Jmse 12 00886 g013
Figure 14. Bayesian updating results of the Gumbel distribution parameter μ n for (a) surge and (b) roll motions in Case 10.
Figure 14. Bayesian updating results of the Gumbel distribution parameter μ n for (a) surge and (b) roll motions in Case 10.
Jmse 12 00886 g014
Figure 15. Original and updated PDFs of the extreme motion responses for (a) surge and (b) roll motions in Case 7.
Figure 15. Original and updated PDFs of the extreme motion responses for (a) surge and (b) roll motions in Case 7.
Jmse 12 00886 g015
Figure 16. Original and updated PDFs of the extreme motion responses for (a) surge and (b) roll motions in Case 10.
Figure 16. Original and updated PDFs of the extreme motion responses for (a) surge and (b) roll motions in Case 10.
Jmse 12 00886 g016
Table 1. The properties of the mooring system.
Table 1. The properties of the mooring system.
PropertyUnitValue
Water depth m100
The radius of the anchorsm640
Draft of mooring pointsm20.6
Number of lines -3
Connection point outside wallm0.8
Segment name (from top to bottom)-BridleMainline
Segment lengthm50609.7
Nominal diametermm132147
Young’s modulusMpa53,94143,624
Segment dry mass per meter kg/m348.5432.2
Weight in water per meterKN/m3.4033.8707
Table 2. A summary of the properties of this FOWT.
Table 2. A summary of the properties of this FOWT.
PropertyUnitValue
Draftm77.6
Displacementtons11,754
Dry masstons11,483
Mooring tensiontons270.9
XYZ
Center of gravitym−0.140−50.03
Center of buoyancym00−42.04
Table 3. The details of environmental conditions of Cases 7 and 10.
Table 3. The details of environmental conditions of Cases 7 and 10.
Case No.Wave HightWave PeriodWave DirectionWind SpeedWind DirectionCurrent SpeedCurrent Direction
74.4 m10.9 s17 deg13.7 m/s11 deg0.21 m/s19 deg
103.9 m8.3 s174 deg30 m/s213 deg0.27 m/s27 deg
Table 4. Comparison of the results from numerical simulations and monitoring data.
Table 4. Comparison of the results from numerical simulations and monitoring data.
Simulated DataMonitoring Data
Natural PeriodDamping RatioNatural PeriodDamping Ratio
Surge91.630.0434960.0445
Sway91.60.0427960.0441
Heave26.280.027225.80.025
Roll34.360.024533.70.0223
Pitch34.590.023933.70.0233
Yaw12.890.0414130.0425
Table 5. The mean and standard deviation of extreme motion responses (surge and roll) in Cases 7 and 10.
Table 5. The mean and standard deviation of extreme motion responses (surge and roll) in Cases 7 and 10.
Monitoring DataMeanStandard Deviation
Surge (Case 7)7.57991.1507
Roll (Case 7)1.23240.3547
Surge (Case 10)2.06990.8446
Roll (Case 10)−0.36010.5436
Simulated DataMeanStandard Deviation
Surge (Case 7)8.10880.6490
Roll (Case 7)0.88400.5725
Surge (Case 10)2.20420.6504
Roll (Case 10)−2.21831.0627
Table 6. The parameters of the Gumbel and Lognormal distributions (surge and roll) in Cases 7 and 10.
Table 6. The parameters of the Gumbel and Lognormal distributions (surge and roll) in Cases 7 and 10.
Likelihood Function α n Distribution Type
Surge (Case 7)1.9763 Gumbel
Roll (Case 7)2.2402 Gumbel
Surge (Case 10)1.9719 Gumbel
Roll (Case 10)1.2374 Gumbel
Prior Function ζ λ Distribution Type
Surge (Case 7)0.09982.0513Lognormal
Roll (Case 7)0.0998−0.4728Lognormal
Surge (Case 10)0.09980.6429Lognormal
Roll (Case 10)0.09980.5628Lognormal
Table 7. Reduction in uncertainty using Bayesian updating in Cases 7 and 10.
Table 7. Reduction in uncertainty using Bayesian updating in Cases 7 and 10.
Case No.Wave HightMotionPrior
Mean
Prior Standard DeviationPosterior
Mean
Posterior Standard DeviationReduction in Standard Deviation
7 Variable _   μ n Surge7.8210.7816.420 0.07091.04%
Roll0.6260.0620.9810.070−12.9%
Extreme motion XSurge8.1090.6496.7120.6401.39%
Roll0.8840.5731.2400.5700.52%
10 Variable _   μ n Surge1.9910.1911.9900.06963.87%
Roll1.7640.1840.9550.05570.11%
Extreme motion XSurge2.2040.6502.2830.6500
Roll−2.2181.063−1.4261.0342.73%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, N.; Zou, G.; Feng, Y.; Ali, L. Probabilistic Prediction of Floating Offshore Wind Turbine Platform Motions via Uncertainty Quantification and Information Integration. J. Mar. Sci. Eng. 2024, 12, 886. https://doi.org/10.3390/jmse12060886

AMA Style

Li N, Zou G, Feng Y, Ali L. Probabilistic Prediction of Floating Offshore Wind Turbine Platform Motions via Uncertainty Quantification and Information Integration. Journal of Marine Science and Engineering. 2024; 12(6):886. https://doi.org/10.3390/jmse12060886

Chicago/Turabian Style

Li, Na, Guang Zou, Yu Feng, and Liaqat Ali. 2024. "Probabilistic Prediction of Floating Offshore Wind Turbine Platform Motions via Uncertainty Quantification and Information Integration" Journal of Marine Science and Engineering 12, no. 6: 886. https://doi.org/10.3390/jmse12060886

APA Style

Li, N., Zou, G., Feng, Y., & Ali, L. (2024). Probabilistic Prediction of Floating Offshore Wind Turbine Platform Motions via Uncertainty Quantification and Information Integration. Journal of Marine Science and Engineering, 12(6), 886. https://doi.org/10.3390/jmse12060886

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