Offshore Structural Reliability Assessment by Probabilistic Procedures—A Review

: Offshore installations must be built to resist fatigue as well as extreme forces caused by severe environmental conditions. The structural reliability analysis is the popular practise to assess a variety of natural waves determined by the long ‐ term probability distribution of wave heights and corresponding periods on the site. In truth, however, these structures are subjected to arbitrary wave ‐ induced forces in the open ocean. Hence, it is much more reasonable to account for the changed loading characteristics by determining the probabilistic characteristics of the random loads and outcomes responses. The key challenges are uncertainties and the non ‐ linearity of Morison’s drag element, which results in non ‐ Gaussian loading and response distributions. This study would analyze advances achieved to date in a comprehensive probabilistic review of offshore fixed jacket ‐ type platforms.


Introduction
Sea surface elevation is always random phenomena that continuously change in space and time [1].An analytical mathematical model is not reasonable to describe the actual condition of sea surface elevation, but stochastic procedures can make an illustration of the sea surface elevation [2].For practices related to sea state statistics, the sea surface is usually described by at least two parameters from either wind speed, wave height, wave period, peak period, or wavelength etc. [3,4].In this study, Met-Ocean historical (or hind-cast) data of wave height (H) and wave period (T) is considered.
As stated by Vanem, et al. [5], through statistical wave analysis, it is possible to describe the sea surface elevation using the aforementioned parameters.Once the sea surface elevation is known, the wave-induced loads and its responses of offshore structures can be simulated regarding the schematic of ocean wave mechanics [6,7], as described in Section 3.3.In order to obtain the extreme structural responses induced by wave loads, various approaches have been proposed for defining the maximum values according to its return period [8,9].
Therefore, the prediction of extreme response values by long-term analysis related to a specific return period is utilized as design waves [10,11].Usually, the design wave based on API standards consider a return period of 25 years, 50 years or 100 years for extreme occurrences, which refers to structural responses (i.e., base shear, overturning moment etc.) that are used in the offshore structural assessment [12].Specifically, the structural assessment can be implemented either by employing a deterministic or probabilistic approach.This approach, associated with the calculated design wave, is described in Section 2.1 below.

Overview of the Offshore Structural Assessment
A fixed offshore platform is exposed to different types of loads in the ocean environment.Usually, the loads acting on offshore structures can be classified into dead (permanent) loads, live (operating) loads, environmental loads, deformation loads, installation (construction) loads, and accidental loads [13][14][15][16].In this study, the harsh environmental load conditions such as wind, wave, current including marine growth effects are considered as the influential loadings in the design of offshore structures [17,18], in which the wind-induced wave load impact is the most dominant loads [19].Thus, the estimated structural response is required to be used in the structural reliability assessment in order to check the structure integrity of platforms.
Structural integrity is a part of design engineering, which deals with the capability of a structure to sustain a designed structural capacity (weight, force, etc.) without breaking structural components [20,21].Structural failure denotes the deterioration (loss) of structural integrity or the loss of the load-carrying ability of a structural component.Ensuring sustainable structure lifespan, predicting wave loads and structural responses on the fixed offshore platform are required for assessing the structural strength capability [22,23].The structural strength (resistance) of the offshore platform is a significant part of reaching a reliable and vigorous design exceptionally capable of withstanding the extreme environmental loadings [21].
In order to achieve safe conditions, one of the plans to succeed in this task is by implementing a structural reliability analysis [24].Structural reliability analysis (SRA) is implemented to guarantee the required design loads (structural responses) and its equivalent resistance load for the offshore platforms [25].Moan [26] presented that offshore platforms are designed and assessed by reliability-based assessment to prevent any structure collapses based on the safety and structure integrity specifications.

Structural Reliability Analysis for Fixed Offshore Platforms
The idea of safety may be described as the non-existence of accidents, disasters, or failures.Safety and reliability are associated with the design stages of the structure's longevity (lifetime) [27][28][29].As explained by San Tint [30], a profound level of safety is the probability that the offshore structures or other components will remain safe for the specified period.The safety assessment is carried out by performing probabilistic-based structural reliability analysis, which requires the coupled loads of design (stress) and resistance (strength) as main models [31].Such probability models will be then assessed in limit state function to achieve the platform's probability of failure (POF) in specific comparative periods [32].
In practice, the design of the structure is determined related to the return period (e.g., 100-year event) [33].In order to ensure the structure design meets the crucial prerequisite of acceptable SRA, the safety performance will be analyzed with regard to the POF associated with the reliability index [34,35].Thus, the SRA can be described as the probability of the structure having proper appearance throughout its lifetime [36,37].Figure 1 illustrates how the probabilistic approach is used in the structural assessment, which is the calculation of the extreme structural responses adopted to assess the structure failure (POF of offshore structures).The probabilistic approach was preferred in this study to account for all sorts of uncertainties in the stochastic load and resistance modelling, whereas the deterministic approach inaccurately assesses when dealing with random processes.For performing the structural assessment, the computer numerical analysis is initiated by calculating the offshore hydrodynamics (wave motion, acceleration, loads), which are imposed on the platform model based on Met-Ocean data provided [38,39].Based on the probabilistic approach [40], the measurement of safe structures between design and resistance loads could be determined.From this gathered data collection, the offshore structures that depend on the data records of stress (design) loads caused by environmental loads (i.e., wind-induced wave and currents) is calculated first [41].Then, the strength (resistance) loads that occur due to the environmental load, gravity load, material, and geometric variability is defined below ultimate limit state [42,43].
To ensure the structural integrity can cope with the extreme waves, the design load (stress model) should not exceed the resistance (strength model) of the platform structures.As displayed in Figure 2, once the stress is greater than the strength, the possibility of the structure's collapses, or POF fails, results in higher chances of failure.A decrease in the POF can be assumed equal to the increase of safety [44].The approximated POF depends on the value of reliability index, which is the space from the origin (mean model distribution) to the limit state function [33,45].Thus, a fair safety margin must be designated with the POF that can regulate the overlap red area scale.Once the longest distance of reliability index is obtained [46], a red zone could be minimized, which decrease the POF directly.In order to evaluate POF of the platform, a simple design procedure based on the failure threshold (return period) of the probability distribution of structural responses is considered [36].The procedure is on the basis of the structural reliability method following the given Equation ( 1), in which the POF function is described accordingly [20]: where POF is the probability of failure, P is the probability distribution (model), R is the system's "resistance" (or ultimate limit state/strength), L is the system's "loads" (or assessed structural response distribution from environmental load effects), and G is the function of failure, which could be represented as G = R − L.
In reliability assessment, the POF was developed in reliability engineering to analyze structural safety and execute reliability analysis [47,48].This means that when the POF number declines, the offshore structure and other structures become more reliable.In this context, reliability may be described as a 1 in 100 or 1 failure in 100 probability [46].Hence, a static pushover analysis or a nonlinear collapse analysis such as SACS and USFOS is extensively employed by various oil and gas organizations in order to obtain the RSR and POF values based on ultimate base shear over a design return time [37,49].In this analysis, the probabilistic approach-based probability, frequency, and time domains are discussed in more detail.Regarding load analysis procedure, the wave load-induced responses are then simulated based on numerical wave modelling

Determination on the Wave Force Procedure
In order to calculate the extreme wave and its related current forces for static analysis in correct sequences, the typical requirement is provided by the American Petroleum Institute (API) standard according to the recommended practice (RP) of 2A (API RP 2A).Since the Met-Ocean data is associated with this standard [50], the API RP-2A procedure is utilized, as illustrated in Figure 3. Thus, this wave force procedure including the twodimensional wave theory, the current kinematics, the surface roughness led by marine growth is referred to in calculating the wave force in Morison's equation.All related topics will be discussed in Section 3. Based on this procedure, it is universal enough to be applicable in all structures for the estimation of extreme wave forces.Due to the impact of the random nature of waves, the calculation of extreme responses could be realized by stochastic processes using specific ocean-wave modelling and its statistical analysis [51].Before proceeding into the random processes, each stage of calculations should be followed orderly related to the API RP-2A guidelines, as explained in Section 3.1.

Analysing Hydrodynamic Wave Load-Induced Offshore Structural Responses
As a practice, the sea surface elevations (surfaces waves) could be acquired by three methods: theoretical analysis, either field experiments or laboratory experiments and digital simulations [52][53][54].In this study, the analysis focuses on the digital (numerical) simulations, which include hybrid time-frequency domain simulation as the basis.Under these domains as per Kareem and Kijewski [55], the Fourier transform method is used to transform the time domain into the frequency domain.Besides, the frequency domain (wave spectrum) could be transformed into the time domain (broadband form) using the inverse Fourier transform.This subsequent section explained briefly the numerical process of calculations related to statistical analysis.

Numerical Simulation Analysis
In the 1950s, Phillips' [56] effort was based on dimensional considerations and on the fact that due to the wave steepness, the high-frequency part of the spectrum is limited by wave breaking.Energy inputted by wind in these frequencies is lost by wave breaking.For that reason, Phillips elaborated the equilibrium range of spectrum in the portion of higher energy (in the band of frequencies from 1.5×fp to 3.0×fp) and is independent of the wind speed and fetch.Phillips' equation is not used in practice, but the form became the basis for most subsequent wave spectra.
According to wave spectral formulations, a spectral analysis based on the fast Fourier transform will be carried out on all the times series for each wave record analyzed [7].Conforming to Faltinsen [57], the wave elevation for a wide area where the sea is traveling along the x-axis, can be described as a summation of sinusoidal wave components.By the fundamental theory of waves [58][59][60], numerical analysis can generate the sea-wave surface elevation of irregular sea states until obtaining the wave loads and offshore structural responses.This is the key concept for describing the waves of the ocean regarding the ocean wave mechanics scheme [6,61,62], as illustrated in Figure 4.

Short-Term and Long-Term Analysis
Analyzing the hydrodynamic flow field (wind, wave, and current) exerted on structures that the wave-induced load varies time in nature [63].It is important to differentiate between short-and long-term definitions in the characterization of the sea climate.In the short-term definition, the state of the sea, which is called stationary, is described for a short period of time.In the long-term definition, the aim is to explain the transition in short term characterization over a long period of time, e.g., over the lifespan of the structure [64,65].
The short-term definition of the offshore platform architecture is used explicitly for the estimation of the wave loading.The long-term definition can be used to decide the most important definition of the short term sea condition in an investigation of the first excursion loss or in fatigue analysis to determine statistical variance in the different shorter sea states encountered by a platform [66].Sea waves arising in a short time scale more closely follow the response time of fixed offshore structures and are described by deterministic definitions.The longer time frame correlates mostly to the structure's fatigue fracture and is defined as probabilistic.
In term of structure design, the short-term analysis is utilized as a first step to evaluate the probability of failure for a particular sea state in a specified timeline.Meanwhile, the long-term analyses take more realistic account of the total sea state condition (i.e., met-ocean data based on scatter diagram).The probabilistic approach could be used to calculate the design responses based on several domains.

Probabilistic Approaches for Fixed Offshore Structural Response Analyses
A systematic technique to calculate all statistical data associated to uncertainties is a probabilistic analysis method [67][68][69][70].Because of its capability on uncertainty, probabilistic analysis is suitable for simulating wave characteristics related to randomness mechanisms such as geometry, material properties, and measurement (i.e., wave load-induced responses) [71][72][73].As reported by Karadeniz [74],Reid and Naess [75], Moon, et al. [76], uncertainty is a significant aspect that influences the accuracy of simulations in order to estimate the numerous simulated data of extreme structural responses.
Therefore, a probabilistic approach is the most encouraging way to predict sporadic wave-inducing forces in ocean environments [77].The probabilistic method is more reliable in simulations [78], as it can quantify arbitrary wave load-induced responses in the ocean environment.Several different techniques have been presented for response assessments.In general, probabilistic structural analysis can be established in three modes; probability, frequency, or time domains [31,[79][80][81].Sections 3.1-3.3give a detailed background of domains under probabilistic structural response analysis.

Probability Domain Methods
The basic concept of the probability domain used to calculate the responses was based on the assumption of wave forces from the theoretical Pierson-Holmes distribution [82].In line with [83,84], the hypothesis began with the mathematical model like wave spectrum, linear wave theory, water particle, wave load, and responses had been expressed into a simplified term of the variable used in the probability distribution function (PDF).By the probability and stochastic process, the hydrodynamic response was simulated employing the Gaussian technique based on the revised PDF [85].Numerous simulated data are then tabulated into the probability density histogram.In order to determine extreme response values, the simple statistical techniques of second and fourth moments were analyzed [86].
In the probability domain, Burrows [87,88] revealed that the Second-Fourth Statistical Moment (SFSM) approach alongside the Principal Component Technique (PCT) had used a linear wave theory, but the Simple Sampling Technique (SST) used the non-linear wave theory.All these approaches were tested on the structures by quasi-static responses.In terms of statistical analysis, only the SFSM method had performed both short-term and long-term analysis.However, the fourth statistical moment approach is extremely computationally challenging in assessing the structural responses.Therefore, Najafian,et al. [89] implemented the development of the random sampling method, allowing the scattered wave loadings to be idealized from a wide range of point loads.The technique is, however, effective only when used in high seas.At low seas, the unrelated water particle kinematic on the surfaces is assessed with a greater sample size accurately.
On the other side, Najafian, et al. [90] presented the SST technique in order to perform effectively in both high and low sea condition.A 'Gaussian Booster' enhances the sampling process by taking benefits from an auxiliary variable in the Gaussian response.Without taking into account the current impacts, load intermittency in the splash zone, and wave directionality in the process, this approach managed to significantly reduce computational effort for quasi-static responses consistently with a slight variability in the kurtosis study.For quasi-static reaction, this technique has proved accuracy and highly effective.
None of the methods considered the current effect in their analyses, Najafian [91] introduced the PCT method, which was tested in the presence and absence of current propagation including load intermittence.Researchers have suggested a simpler way to minimize complexity with the use of the PCT for various random variables.A set of uncorrelated variables known as principal components that are linear to the initial variables characterize much of the correlated (clustered) random variables.It will then produce an individual response that is counter to the associated time simulation and results in a smaller sampling variability.Yet, the development considered for quasi-static structural analysis is neither short-term nor long-term analysis.Currently, under this domain, as presented in Table 1, the PCT system was the latest method for quasi-static assessment.

•
The application of PCT approach is mostly significant for the assessment of quasi-static structural responses

Frequency Domain Methods
Stochastic processes in the frequency domain had been previously practiced in the offshore oil and gas industry [92,93].On the other hand, the frequency domain becomes another option to minimize computational demands, although it is less precise than the time domain [94,95].This domain has been widely assigned to the marine structures, mainly offshore manufacturers.Marine designers in the offshore industry applied one such Volterra series-based model.
This series-based model enabled an estimation of offshore structural responses correlated with cumulant spectra (e.g., unispectrum, bispectrum, and trispectrum) [96][97][98].The application of the Volterra model requires that polynomials are used to approximate the nonlinear Morison's drag term for defining the transfer functions.The model has two shortcomings that drastically decrease its applicability amid rising polynomial degrees: exciting increment of complexity in obtaining higher-order spectra [99] and immediately expanding numerical problems in computing the higher-order statistics [100].
In 1997, Naess and Pisano [101] performed research on dynamic structural analysis by applying a stochastic method.In his study, the current propagation is inconsiderable, and the splash zone is neglected.The recommended approach presented a correct forecast of the response, although it was restricted to linear structures.Kareem, et al. [102] explained that the nonlinearities' behavior forms an essential issue in the probabilistic method by using Stokes' perturbation expansion solution.Nevertheless, the investigation found that the platform is not subjected to the higher responses.
Therefore, it can be stated that the Volterra series can quickly estimate the nonlinearities of an uncomplicated system.Zheng and Liaw [103] also worked on the nonlinearity issues of wave forces and structural responses by a frequency-domain cumulant spectral analysis method.Ersdal [104],Puskar, et al. [105] first discussed a major problem on wave loading evaluation applying the Ultimate Force for Framed Offshore Structures (USFOS).This approach increased its efficiency in calculating the structural failure using USFOS nonlinear analysis.On the other hand, Marzocca, et al. [106] found Volterra functional series-spectrum analysis to solve the higher order in complex adaptive (dynamic).This theoretical approach is a better guide to help, for example (1) how non-linearities relate to the dynamics of the non-linear system, (2) how the computational methods of prediction are derived, (3) how the non-linear system analysis techniques are oriented, and (4) details on bispectrum responsiveness and the reliance on amplitude of such nonlinearities.
Carassale and Kareem [107] improved the technique utilizing the 5th Volterra model to determine the higher order of the nonlinear system.Even so, the model was restricted at a specific condition that affected reliability assessment.Agarwal and Manuel [108] presented the theory for a nonlinear irregular model and use it as part of the widely used wind turbine simulation program, FAST (Facilitated Application Specification Technique software).A powerful algorithm is used for nonlinear wave elevation and kinematic estimation, allowing for simple performance of large numbers of time-domain simulations for long-term load prediction with statistical extrapolation.The analysis focuses on higher estimation on the use of second-order wave theory than the first-order extreme offshore response in structures.
For reliability evaluation, Norouzi and Nikolaidis [109] introduced an improvised Probabilistic Re-analysis method.This method is based on the idea that a small failure probability can be calculated as a product of larger conditional probabilities of intermediate events.Consequently, this method is more efficient because it is much faster to calculate several large probabilities than a single low probability.Subsequently, Ebrahimian, et al. [110] presented a more detailed review of influencing factors that have an effect of uncertainty in the assessment process.However, this PIWA technique-based frequency domain has lower accuracy in relation to a complex nonlinear system such as more severe environmental conditions [110].
On the other side, when the wave height is equal to the diameter of the cylindrical structure, nonlinear effects could be excited and influence the steadiness oscillations by exhibiting steady or transient response.Furthermore, first-order wave theory, which uses linear constraints, does not explain wave kinematics in the wave field and then leaves out loads in the wave zone.Solving these issues, Chen [111] introduced the frequency range to be defined by the higher order terms added by the procedure that can best be calculated.Better estimation for extreme offshore structural responses in relation to first order for second-order theory.
The simulations based on linear wave theory are currently considerably less predictive than nonlinear analysis [110,112,113].While non-linear wave analysis still predicts offshore structural responses as a result of the wave load, it is more practical compared to linear wave analysis, as seen in the measurements and laboratories [114].At least the addition of the second nonlinear term is necessary to provide the stochastic characteristic of the actual sea state, thus producing a better outcome for predicting the precision of extreme value responses [112,115].In this respect, the contribution of higher-order Stokes modes to the linear approximation makes it more complicated and difficult to solve the equation.
Kim [115] developed a NARX (Non-linear Autoregressive with Exogenous Input) method.With this approach, a different higher-order frequency component is precisely calculated in the time series of the structural scheme responses, even the presence of random excitation.Following Mérigaud and Ringwood [116], an NLFD framework has the potential to be implemented on quadratization and cubicization of the nonlinear wave loads that facilitate great computational conservations in the estimation.In probabilistic based on frequency domains, Table 2 listed all the approaches with its respective results.

•
PIWA technique by the application of a Collapse Wave Height (CWH)-based approach, which is suggested to measure the mean annual frequency (MAF) above the Collapse Prevention (CP) limit condition.
• Decreased precision when faced with a dynamic non-linear system, such as more severe ocean conditions.

•
The frequency spectrum to measure the structure can be better described with higher order terminology added in wave analysis Kühn [117] proposed a simplified method for fatigue evaluation.Based on the frequency domain as a basis, the calculations of aerodynamic wave load on the structure could be reached via superimposing time-domain simulations.Later, Van Der Tempel [118] employed frequency domain analysis to achieve structure design, where a judgment of time-domain outcomes proved to be consistent and accurate.Recently, Seidel [119] proclaimed a profoundly uncomplicated procedure to determine if structural responses induced wave loads on the structures.

Time Domain Methods
In general, frequency-domain methods are efficient, but not as accurate as timedomain simulations.Time-domain simulation is applied to facilitate nonlinear analysis with linked simulations of wind and wave [120].The time-domain analysis operated as the most accurate and regularly managed for verification purposes.Higher accuracy in results is an essential part of the analysis, particularly for proving ] time-domain methods.However, time-domain analysis demands an excessive volume of computer run-time due to the nonlinearity consideration in calculations [121].
Previous researchers have revealed multiple time-domain methods, including the Monte Carlo Time Simulation (MCTS), NewWave theory (NWT), Constrained New Wave (CNW), Finite-Memory Nonlinear System (FMNS), Efficient Time simulation (ETS), and Efficient Threshold Upcrossing (ETU).Beginning with Orkin, et al. [122], this study describes the way that the MCTS procedure is used to help design an offshore terminal system in a very hostile environment.By structural analysis, nonlinear excitation under random responses needs the MCTS method with hundreds of thousands of simulated records [123].In their study, the MCTS method accurately evaluates the responses of fixed platforms, which are used as benchmarks in deriving the short-term and long-term probability distribution of extreme responses.
Further research work is also required in the model development to improve the accuracy towards MCTS methods.In addition, many researchers tried to improve the level of efficiency of simulation-based methods.Tromans, et al. [124] were the first to present the NWT in a time-domain method.This scheme involves the derivation of the probability distribution of extreme responses.In this method, the separation of the NewWave crest elevation is computed, employing the auto-correlation function of a sea state.The rationale of autocorrelation is to interpret the Fourier transform of the frequency spectrum.On the other side, Taylor et al., (1997) introduced a technique by the addition of two non-random functions.This technique considers short constrained simulations in order to determine the largest response correlated with a specified crest.
In order to overcome this previous problem, the NWT was expanded to a second-order approach that analyses the non-linearity of random waves [125][126][127].After that, an idea of CNW theory is developed by [128], which is produced by mathematically constraining a NewWave within random time series, allowing the irregularity of ocean waves to be considered.These wave theories have been extended in this thesis to include 2nd-order and directionality effects, and their formulations have been written into a new Fortran code for calculation of the water surface and water particle kinematics.The effects of irregularity, directionality, and nonlinearity of ocean waves on dynamically sensitive structures are then shown for examples for the mobile jack-up drilling platform.
Cassidy, et al. [129] developed a probabilistic model to calculate the statistical extreme response and failure of jack-up structures considering different sources of uncertainties in loading and structural and foundation modelling and demonstrated that the uncertainties in the values of some of the essential variables could significantly affect the extreme response statistics.CNW theory and an advanced model to consider spud can-soil interactions were used to perform a dynamic evaluation on sample jack-up offshore platforms.In further studies, an association MCTS method is linked with a resulted linear extrapolation procedure in a central processing unit (CPU) for time-saving.
Naess, et al. [130] introduced a simple and robust technique in order to compute the extreme offshore structural responses without any simplifying exceeding those indirectly expressed in any computers.Linear extrapolation technique applied MCTS approach and optimized it with linear extrapolation for more accurate prediction.This method is an extension of MCTS method utilizing linear extrapolation procedure followed in the time-saving CPU to determine extreme response statistics.It is robust and straightforward, which counts all nonlinear influences without any oversimplification exceeding these embedded in any computer models.The outcome was authenticated by utilizing the MCTS method by varying sea states for the 100 years.This technique was used as a structural analysis for jacket platforms, which is more accurate, while at the same time, diminished sampling variability relies upon the same number of simulated records as that of MCTS method.
Najafian [131,132] later applied the same approach to predict the extreme response values and introduced a finite-memory nonlinear system designed for a further efficient procedure.Najafian and Zaki [133] announced a new procedure called the FMNS method.Note that the proposed technique offered an efficient procedure in simulating the structure response and provided good accuracy on the prediction of extreme response values.This procedure comprises three major systems, which are system A (linear system), system B (zero-memory nonlinear system), and system C (linear system).The result of FMNS method demonstrated to be more efficient than the MCTS procedure by a scale of 35 times.However, the outcomes indicate the level of accuracy reduced at low wave height.
In order to improve the FMNS technique, Najafian and Zaki [133] and Mohd Zaki, et al. [134] proposed a new approach using a modified finite-memory nonlinear system (MFMNS).This was done by separating the structure into two zones and applying the FMNS technique separately based on the zone, and the responses obtained were added to be employed to generate the probability distribution for the purpose of predicting an extreme response value.This was required to ensure a good correlation of water particle kinematics at the corresponding zone due to the shorter wavelength.However, the accuracy reduced once the current was considered in the analysis.The relationship of responses based on the corresponding technique was demonstrated.
On the other hand, Abu Husain and Najafian [135] pioneered an ETS method to resolve the sampling variability obstacle in the MCTS method.The ETS method uses the superiority of the excellent relationship between two variables of extreme surface elevation and extreme structural responses.To end this, the ETS method was compared with the benchmark of the MCTS method by deriving them into the probability distribution of extreme responses.The ETS method was validated with the results of the MCTS method.Thus, the ETS method analyzed two versions of the wave conditions.As a result, this ETS method was revealed to be a more efficient approach for determining high sea state intensity than the low sea state.The ETS method offered a more precise result at high sea state compared to the low sea state.
The unique approach merges the speed of the ETS method with excellent accuracy based on threshold up crossings identified as Efficient Threshold Upcrossing (ETU) by [136].This approach is related to the ETS system, but it removes more numerous sample data records.The result of the ETU method was insignificantly more constant than the ETS method.However, the ETU method is restricted to the quasi-static structural responses without consideration of wave-current and intermittency load around the splash zone.
In a previous study, Mallahzadeh, et al. [137] conceived an extremely efficient simulation technique known as the Relationship Time Simulation (RTS) procedure.At the beginning of its appearance, this RTS technique is established by taking advantage of their relationship from ETS methods [138].Based on this relationship, an empirical model was derived from the given sample data to be used in the probability distribution to predict the extreme structural responses.Mallahzadeh's ideas were proposed as a new concept from efficient time simulation (ETS) procedure in order to extemporize current ETS methods [139].
This concept was known as an Efficient Time Simulation-Relationship (ETS-RTS) method employed to improve accuracy and efficiency from MCTS method.Due to the study limitations, Mallahzadeh's study was focused only on the single-legged structure problem, and the method was only validated for high sea state intensity.Thus, the pro-posed approach used in this study was from this conceptual study, which expended more details on cases of wave conditions including criteria consideration.Table 3 outlined the entire procedure with its associated performance in the probabilistic analysis based on time domains.• ETU method has an ability to be up to 3 grades of magnitude better than the MCTS method.ETU method was slightly more stable than the advanced ETS method.As noticed, an estimation approach based on the MCTS method was introduced recently in [140] with the intention of minimizing the computational time.For the application as an example, buckling collapse strength, Gaspar, et al. [141] proposed an MCTS method associated with a model of structural finite element process.With a nonlinear structural finite element process, the redundant computational time in order to assess the strength of such a vessel hull girder could be reduced.Regarding excessive samples of the wave, Aarland [142] proposes an alternative to FFT by presenting equal-area principle.This approach was implemented in order to prevent repeating the wave history with several thousands of wave components.
Johari [143] introduced a technique for an optimization of the ETS method.The optimization ETS method has been validated with the MCTS method by comparing the short-term probability distribution of extreme responses at several sea states with the appearance of current impacts.This study provides a better understanding of increasing the level of accuracy and efficiency for estimating the extreme offshore structural response.In order to get result with higher accuracy, a large number of groups (simulated response records) will systematically distribute probability of occurrence at the higher segment.Indeed, several enhanced simulated response records will decrease the level of efficiency of simulations.
Mukhlas [144] introduces a new model by combining the FMNS approach with a non-linear wave analysis for a more accurate outcome.This model was developed based on the advantage of the efficiency of the FMNS method in deriving the non-Gaussian responses.This invented FMNS method was known as the FMNS NL model.For the recent findings, an ETS-Reg procedure was conceived to fulfil limitations of the ETS method, including to evaluate the level of efficiency.It is a combination of two methods of the conceptual ETS-RTS method and optimization of ETS method.Regarding conceptual ETS-RTS method, the proposed model was developed based on the excellent correlation found between the extreme values of surface elevation (input) and their corresponding responses (output) [138,139].Meanwhile, the suggested number of simulations equal to N = 260 from the optimization of the ETS method would be used in the model development [143].Hence, a simplified analysis-based time domain is introduced to overcome this issue by applying direct calculation procedures.

Overall Review of the Probabilistic Domains on Offshore Structural Assessment
Based on [20], the structural reliability analysis (SRA) and the quantitative risk analysis (QRA) are two categories of reliability assessments.The SRA intends to describe the POF due to the degradation mechanism (i.e., coating damage, corrosion, cracking, deformations, and changes in material properties) [145,146].Also, the degradation process connected with the uncertainties should be considered in structural analysis.Meanwhile, the QRA is conducted based on decision-making to form a framework [147].This framework is, along with approximations on the possibility of risks, endangering civilised (human) lives and health, assets, and commercial businesses [148].
Since the SRA approach is concerned with structure analysis, the study emphasizes that the probabilistic domain is a method for estimating the extreme structural responses as the POF value is calculated [149].Under the SRA approach, several methods have been used to calculate the load impact on structures.Wave loading can be performed using two approaches, either applying a deterministic or probabilistic technique [32].The deterministic approach would not take uncertainty into consideration during calculation implementation and corresponds only to normal (regular) waves [150].Due to this restriction, another approach is introduced using probabilistic modelling on computing the wave estimation, which is to be more practical than the deterministic process when determining the hydrodynamic loadings of the structures [151].
As described in the previous Sections 3.1-3.3,three main types of domains such as probability, frequency, and time are categorized under the probabilistic approach.Each domain may have its benefits and drawbacks depending on the performance of method potential.Table 4 has summarized all the methods in probabilistic analysis with its corresponding output.The application of those methods was taken into account regarding the type of wave model, system structure be it evaluation either quasi-static or dynamic structures, statistical analysis (i.e., short-term or long-term), and the wave condition of current impacts.
Generally, the probability domain was used as a fundamental concept in the range of the first to the fourth statistical moment of their calculations [86,87].Essentially, the probability domain is quite similar to the frequency domain in terms of efficiency, but the limited study shows that the level of accuracy is obviously far behind the frequency domain and time domain [152].For that reason, there was not much effort that was expended in this domain to improve the previous methods.On the other hand, the frequency domain is what the industry is most looking toward, because it gives an accurate prediction with the efficient outcome that considers real structures [153,154].In comparison to the time domain, the frequency domain does not take a long process of calculation in order to complete the full offshore structure assessment [155].
However, the frequency domain approach does not consider the nonlinearities in its calculation [102].Due to that reason, the prediction of a specific critical condition was sometimes under-predicted or over-predicted, which gave inconsistent actual values [156,157].Thus, to improve these prediction values, the probabilistic-based time domain was preferred.This is because the time domain considers all sort of nonlinearities in the calculation.Also, it can consider the randomness process from the wave, loads, and the effect of current and velocity kinematics, including the intermittency wave loadings around splash zones [158][159][160][161]. Thus, the time domain was a more accurate method compared to the probability and frequency domains [161][162][163].

•
The fundamental principle of the design and analysis of offshore structures is that it enables their stability over design life to be maintained.It follows from the non-linear behaviourof fixed offshore platforms that they involve a considerable number of major uncertainties pertaining to the techniques of structural analysis, the mathematical models for the structural behaviour, the loading mechanism, and the response to the parameters concerned.

•
The probabilistic approach has been shown to have several benefits compared to the deterministic approach.It offers, in particular, a clear framework for calibrating code design practices like the evolving load and resistance factor design (LRFD) methods.This will explain a decrease in the safety factors associated with the individual components through a clear reflection of the distinctive aspect of the wave environments and hydrodynamic response.

•
There are three fields in which the response properties can be built on the basis of the probabilistic method;the probability, frequency, and time domain.The principal challenge in probabilistic analysis is the randomness of stochastic waves.An exact model for predicting extreme responses is important as awareness is necessary if the structural design of the offshore is to be optimised.

•
With regard to design structural analysis, establishing a practicable technique that combines the latest short-term and long-term non-linear analysis methods with the conventional design wave methods seems desirable in a manner compatible with the reliability analysis.

•
The long-term probability distribution, which is then required to analyse fatigue concerning the structural extreme during its service life, is needed to analyse the first excursion failure.It can be accomplished by transforming their respective short-term distributions into a long-term wave environment.

Figure 1 .
Figure 1.Determination of the probability of failure using the structural reliability analysis approach.

Figure 2 .
Figure 2. Probability of failure of load-resistance normal assumption.

Figure 4 .
Figure 4.The numerical procedure and its statistical analysis of the structural response design based on the conventional methods in determining return periods.

Table 1 .
Probability domain method in evaluating extreme responses.

Table 2 .
Frequency domain method in evaluating extreme responses.

Table 3 .
Time domain method in evaluating extreme responses.
Author Contributions: S.Z.A.S.A. compiled the methods from previous researchers, analyzed the findings, and wrote the original draft paper; M.K.A.H. and N.I.M.Z.contributed to the review, discussion and verification on the results; N.A.M. extracted several methods based on the nonlinear wave modelling and also reviewed the paper; E.M.S. and N.U.A. shared about reliability analysis with industry experience perspective; G.N. contributed as a supervision.All authors have read and agreed to the published version of the manuscript.The paper is financially supported by the Ministry of Higher Education under FRGS (Registration Proposal No: FRGS/1/2018/TK01/UTM/02/15; FRGS/1/2019/TK01/UTM/02/17) and Universiti Teknologi Malaysia (Registration Proposal No: Q.K130000.2456.09G36;Q.K130000.2456.09G41)which is gratefully acknowledged.
Funding:Informed Consent Statement: Not applicable.