1. Introduction
The sea conditions in West Africa are mild all year round, and the combination of floating production storage and offloading (FPSO) + oil offloading lines (OOLs) + catenary anchor leg mooring (CALM) buoy + shuttle tanker has been widely used for deep-water oil transportation in this sea area [
1,
2]. In this combination, OOLs are important structures for oil transportation, with a span of more than 2000 m, and their safety is of great concern [
3]. This paper will study the fatigue characteristics and parameter sensitivity of the OOL according to the environmental statistics in the West Africa Sea area.
Many studies have been published on the fatigue damage of marine pipelines. Liu et al. [
4] studied the coupling effect of direct waves and currents in the context of fatigue damage assessment and proposed a first-order correction factor involving the ratio of current and wave velocities to evaluate the environmental load coupling effect. Arzaghi et al. [
5] proposed a probabilistic method to simulate the pitting and corrosion fatigue degradation process of submarine pipelines by using a dynamic Bayesian network to model the time property and growth rate of change. Zhang et al. [
6] established a numerical model for the fatigue behavior of multiple coplanar cracks in offshore pipelines under cyclic tensile loading and studied the important parameters that affect fatigue crack growth. Al-Janabi et al. [
7] carried out model experiments to study the influence of different parameters on fatigue damage near the touchdown zone of a steel catenary riser at the seabed.
Although these studies mainly focused on the fundamental fatigue characteristics of materials or numerical models on fatigue behavior, most of them are mechanism studies, so the research objects are mostly simple and typical pipeline structures. However, for practical OOL structures in FPSO offloading systems, related fatigue analysis and design optimization research are insufficient. Chaudhury [
8] proposed a “U”-type oil pipeline and made a comparative analysis of the “U”-type oil pipeline and the “W”-type oil pipeline, including strength characteristics, fatigue characteristics, economic advantages, etc. Through this research, it was found that the “U”-type oil pipeline is superior to the “W”-type oil pipeline in fatigue behavior. Odru et al. [
9] analyzed the fatigue characteristics of oil pipelines based on the spectrum analysis method. The results show that the fatigue damage of the pipeline mainly occurs near the buoy. It is difficult for a 22-inch steel pipeline to meet the fatigue requirements, so a buoy submerged in the water is recommended to reduce the wave effect. However, in this analysis, only the first-order wave force is considered. Based on their experience in engineering design in West Africa, Montbarbon et al. [
10] proposed that increasing the effective tension of the oil pipeline is conducive to the improvement of its fatigue life. They also proposed that the fatigue characteristics of the oil pipeline should be considered in the case of shuttle tanker operation. Drumond et al. [
11] conducted a literature review on failure events experienced by the industry concerning pipelines, risers, and umbilical cables and described their causes, consequences, and severity. He et al. [
12] analyzed the effect of pan length on the fatigue life of submarine pipelines, finding that increasing the gap between the suspended pipeline and the seabed is beneficial to enhance the fatigue life of the suspended pipeline.
Due to the fatigue failure of oil pipelines, many accidents such as crude oil leakage have occurred, resulting in serious economic losses and environmental pollution [
13,
14]. Therefore, it is of great importance to improve the fatigue design knowledge of oil pipelines. However, as can be seen from the above overview, research on the fatigue characteristics of oil offloading lines in the FPSO CALM buoy offloading system is insufficient. Although current studies have paid attention to the effects of pipeline layout and buoy position on fatigue damage, the influence of factors such as motion and tension distribution, wave type, wave direction, coverage length of buoyancy material, outer diameter, wall thickness, shuttle tanker, and current loads on OOL fatigue damage is still not fully studied. Therefore, the purpose of this paper is to comprehensively analyze the influence of the above factors on OOL fatigue damage and provide suggestions for relevant design.
The outline of the rest of this paper is as follows: In
Section 2, the methodology for numerical modeling and fatigue calculation is described. In
Section 3, the verification of the numerical model is introduced. In
Section 4, the calculation results are analyzed, specifically the influence of different parameters on OOL fatigue damage.
Section 5 presents the conclusions.
2. Numerical Model of the Offloading System
In this paper, a numerical model of the FPSO CALM buoy offloading system is established according to the principal dimensions of the research object in the experimental investigation previously conducted by the authors’ team [
15]. Firstly, the response amplification operator (RAO), additional mass, and radiation damping of the structures are obtained by the frequency domain analysis using AQWA software (ANSYS, Inc., Canonsburg, PA, USA), and on this basis, the dynamic behavior of the structures are studied using Orcaflex software (Orcina Ltd., Ulverston, Cumbria, UK). Additional damping equal to 5–10% of the critical damping value is added to each degree of freedom of the structures, the specific value of which is determined by comparison with the corresponding free attenuation experimental results, as well as the RAO experimental results, in typical wave angles. Then, the hydrodynamic responses of the structures are calculated, and the accuracy of the numerical model is verified by comparison with the corresponding experimental results.
2.1. Principal Dimensions of the Numerical Model
The schematic diagram of the numerical model is shown in
Figure 1, and the top view of the system is shown in
Figure 2. In the ballast condition, the FPSO and shuttle tanker will have a larger motion amplitude under wave load and have a greater impact on the OOL’s fatigue damage. Therefore, in the fatigue analysis of this paper, the FPSO and shuttle tanker are both in ballast conditions. For principal dimensions and related parameters of the model, refer to the model experiment previously conducted by the authors’ team [
15].
The fatigue analysis in this paper mainly focuses on the OOLs, the specific parameters of which are as follows: The OOL consists of two 20-inch (508 mm) diameter steel pipes, the material of which is Det Norske Veritas (DNV) E460 steel. Buoyancy materials are assembled in the middle of each pipeline to make it form a “W” shape after being subjected to buoyancy, as shown in
Figure 3.
In the spatial arrangement, the two oil pipelines are distributed up and down. The upper OOL is shorter (2140 m) and is connected to the portside of the FPSO and CALM buoy. The lower OOL is longer (2220 m) and is connected to the starboard side of the FPSO and CALM buoy. The detailed attributes of each segment (as shown in
Figure 3) of the OOLs are shown in
Table 1.
The numerical model of the buoy is established by combining the Morison element with the disc element, and the OOLs are connected to the buoy at the top. The numerical model of the buoy and corresponding mooring lines and OOLs is shown in
Figure 4, the detail of which has been introduced in the authors’ previous work [
16]. The other end of the OOLs is connected to the FPSO (see
Table 2 for the specific coordinates of the OOL connecting points on the FPSO and buoy).
2.2. Model of Loads
In the numerical model, the FPSO and buoy are rigid bodies, and their motion equations can be derived from Newton’s second law. The loads applied to the structures are calculated and taken into the motion equations to solve the six-degrees-of-freedom motion of the FPSO and buoy. Refer to Reference [
17] for details. The loads to be calculated mainly include the wave load (both the first-order wave force and the second-order wave force), current load, and the mooring line and OOL’s force.
2.2.1. Mooring Line and OOL’s Force
In this paper, the lumped mass model is adopted for mooring lines and OOLs. The schematic diagram of the line segments is shown in
Figure 5. The line element supports the flexibility of the line to have axial displacement, torsion, tension, and bending. More details can be found in the authors’ previous work [
16].
2.2.2. Wave Load
For the FPSO and shuttle tanker, wave loads are calculated based on potential flow theory [
18], and for the CALM buoy and OOLs, wave loads are calculated by slender-body expressions based on Morrison equations.
2.2.3. Current Load
The current force and moment can be denoted as follows:
where
Fxc,
Fyc. and
Mxyc are longitudinal fluid force, lateral fluid force, and yaw fluid torque.
Cxc,
Cyc, and
Cxyc are longitudinal fluid force coefficient, transverse fluid force coefficient, and yaw fluid torque coefficient.
ρc is the density of seawater (kg/m
3),
T is the average draft, and
VcR is the relative velocity. The coefficients can be determined based on Oil Companies International Marine Forum (OCIMF) [
19].
2.3. Method of Fatigue Analysis
Based on the hydrodynamic coupling calculation results of the system, the S-N curve method, which is widely adopted at present [
20,
21], is used to analyze the fatigue characteristics of the OOLs. According to the requirements of DNV-RP-C203 specification [
22], the DNV-E curve is used for fatigue calculation in this paper. The parameters of the S-N curve used are shown in
Table 3. Based on the above method, the fatigue damage of the OOL is calculated by Orcaflex software. The S-N curve method used in this paper is described as follows:
2.3.1. Calculation of Time-Domain Fatigue Stress
According to the established hydrodynamic numerical simulation model of the FPSO oil transportation system, the time-domain curve of fatigue stress of the OOL is calculated. Assuming that the oil pipeline is isotropic, the fatigue stress can be obtained by the following formula:
where
A is the effective area of the section,
Ixy is the second-order stress moment,
σA is the axial stress,
σM is the bending stress,
T is axial tension,
Cr is the bending moment load factor, and
Mx and
My are the bending moments.
2.3.2. Rain-Flow Counting Method
After obtaining the time-domain calculation results of OOL fatigue stress in the irregular wave, the average stress, stress amplitude, and stress cycle times of each node of the oil pipeline are obtained by the rain-flow counting method.
2.3.3. Stress Correction
The stress obtained by the rain-flow counting method is corrected as follows:
where
Sr is the actual stress,
Sm is the average stress, and
SMTS is the ultimate tensile strength of the material.
The stress concentration factor and wall thickness correction factor are used to consider the influence of welding. The specific formula is:
where
Se is the nominal stress amplitude,
SCF is the stress concentration factor, and (
tfat /
tref)
k is the thickness correction factor, which is applicable to the case that the actual wall thickness
tfat exceeds the reference thickness
tref = 25 mm.
2.3.4. Calculation of Fatigue Damage by S-N Curve
The fatigue damage is calculated according to the following formula:
where
N is the number of stress cycles,
S is the stress amplitude, and
a and
m are material coefficients.
2.3.5. Calculation of Accumulated Fatigue Damage
In this paper, the alternating load of variable amplitude needs to be considered in the fatigue analysis. Therefore, Miner’s linear accumulation criterion is used to calculate the fatigue damage accumulation results under variable amplitude stress cycle. The expression of Miner’s linear accumulation criterion is as follows:
where
n(
Si) is the actual number of stress cycles with amplitude
Si, and
N(
Si) is the total number of stress cycles in the case of fatigue failure of the structure under this stress cycle.
2.4. Environment Condition
The working location of the oil offloading system studied in this paper is the sea area of West Africa. In this sea area, the wave condition is relatively mild, the main wave component is the swell with a long period, and there are wind waves in some parts. Specifically, the wave components in West Africa can be divided into main swell, secondary swell, and wind wave. The Joint North Sea Wave Atmosphere Program (JONSWAP) spectrum [
23] is used to simulate the individual action of swell, and the Ochi Hubble spectrum [
24] is used to simulate the joint action of wind wave and swell. For the single wind wave condition, the wind wave component in the Ochi Hubble spectrum is used for the numerical simulation. The wave angle is considered as a fixed value in this paper.
The working area of the offloading system analyzed in this paper is a certain sea area with a depth of 1500 m in West Africa, and the considered wave distribution is identified according to the observation statistics of the joint wave distribution in West Africa. Three wave conditions need to be considered in this paper: main swell, secondary swell, and wind wave. Their significant wave height (
Hs)—average period (
Tp) distribution characteristics are shown in
Table 4,
Table 5 and
Table 6.
3. Verification of Numerical Model
In order to verify the accuracy of the numerical model, the results of the numerical simulation are compared with the results of the previous experimental investigation conducted by the authors’ team [
15]. Numerical simulation results of CALM buoy motion and portside OOL top tension are compared with the experimental results, including the time-domain curves and frequency domain analysis. The environment condition is the combined action of swell and wind wave, where
Hs (swell) = 2.6 m,
Tp (swell) = 15.9 s;
Hs (wind wave) = 2.05 m, and
Tp (wind wave) = 6.7 s. The wave direction is 0° (as shown in
Figure 1).
The numerical and experimental results of surge, heave, and pitch motions of the CALM buoy, as well as the top tension (in the experiment, a force sensor is set at the connection between OOL and buoy, and the force measured by the force sensor is the top tension), of the portside OOL are shown in
Figure 6,
Figure 7,
Figure 8 and
Figure 9, respectively. As shown in the figures, the numerical simulation results of the CALM buoy’s motion and OOL’s top tension are generally in good agreement with the experimental values, except that the numerical results of pitch motion amplitude of the CALM buoy in
Figure 8 are slightly larger than the experimental results. This should be due to the ignorance of a part of the fluid viscous damping in the hydrodynamic calculation method based on the potential flow theory used in this paper.
Considering that the stress of the OOL is the most important factor affecting its fatigue damage, the numerical simulation results of the OOL
p top tension are compared with the experimental results in detail, as shown in
Table 7. It can be seen that the numerical simulation errors of maximum value, average value, and standard deviation are very small, all less than 5.7%. Only the error of the minimum value is relatively large, reaching 17.1%. Considering that the error of the minimum value is not that relevant because, for design purposes, the minimum top tension is not that important, the overall calculation error is within the acceptable range, proving that the numerical simulation model is reliable.
4. Fatigue Analysis
In this section, the hydrodynamic and fatigue characteristics of the portside OOL are analyzed first, and then the influence of other parameters on the fatigue damage of the OOL is analyzed.
4.1. Hydrodynamic and Fatigue Characteristics of the OOL
In order to study the relationship between the dynamic responses of motion and tension of the OOL and its fatigue damage, the coupling analysis of the system in typical working conditions is carried out. The wave condition is the wind wave, Hs = 1.5 m, Tp = 5.5 s, the wave direction is 0°, and the calculation time is 2000 h.
Figure 10 shows the distribution of maximum motion amplitude of the portside OOL and the fatigue damage, where the abscissa represents the OOL’s segment at the corresponding coordinate position (from FPSO to CALM). As shown in
Figure 10, along the length of the OOL, there are two peaks of fatigue damage. The first peak is located in the middle segments of the OOL, corresponding to the segment covered by buoyancy material (as shown in the red section in
Figure 3). The second peak is near the CALM buoy, where the fatigue damage is much greater than the rest of the OOL. Moreover, there is no obvious correlation between the distribution of motion amplitude and the distribution of fatigue damage. However, in some locations, such as segments covered with buoyancy material or near the CALM buoy, the motion amplitude is smaller due to boundary conditions, but the fatigue damage is greater than the nearby areas.
Figure 11 shows the distribution of tension (variation amplitude and the maximum value) and fatigue damage along the length of the OOL. As can be seen from the figure, the distribution of fatigue damage shows a clear positive correlation with the distribution of tension variation amplitude, but there is no obvious correlation with the distribution of maximum tension, indicating that the tension variation amplitude is a major factor for fatigue damage. Moreover, the growth rate of fatigue damage is obviously greater for segments covered with buoyancy material or located near the CALM buoy. As shown in
Figure 10, both of these two parts of the OOL are subjected to constrained boundary conditions, resulting in a small motion amplitude. This indicates that under similar external force conditions, the constrained boundary conditions may increase the fatigue damage. It should be noted that the maximum stress of the OOL in this case is about 150 MPa, which is far lower than the material’s yield limit of 460 MPa.
4.2. Influence of Wave Parameters on Fatigue Damage
In this part, the influences of the type and direction of wave loads on the fatigue damage of the OOL is analyzed.
4.2.1. Influence of Wave Types
As mentioned above, the wave components in West Africa can be divided into main swell, secondary swell, and wind wave. In order to analyze the different effects of these three wave conditions on the fatigue damage of the OOL, the analysis of parameter sensitivity of different kinds of wave loads on fatigue damage is carried out. It should be noted that the calculation process is to conduct dynamic analysis under each wave condition based on the wave distribution shown in
Table 4,
Table 5 and
Table 6. Then, according to the specific exposure time of each wave condition, the fatigue damage of the OOL under each wave condition is accumulated in proportion.
The fatigue damage of the OOL under three different wave conditions is shown in
Figure 12. It can be found that fatigue damage mainly occurs in the middle parts of the OOL (covered by buoyancy material) and the parts near the CALM buoy. In addition, the fatigue damage of the former is mainly caused by the swell, and the fatigue damage of the latter is mainly caused by the wind wave. Additionally, the fatigue damage caused by the secondary swell is very small.
4.2.2. Influence of Wave Direction
In order to study the influence of wave direction on the fatigue damage of the OOL, the fatigue damage caused by wave load in different directions in the same period is calculated as shown in
Figure 13, where the type of wave load applied is the wind wave, and the calculation method is similar to that in
Section 4.2.1. As shown in
Figure 13, with the increase in wave direction from 0° to 90°, the fatigue damage of the OOL decreases obviously. The fatigue damage at 0° is about 10 times that at 90°, indicating that the surge motion of the CALM buoy caused by the 0° wave loads is the main cause of the fatigue damage of the OOL.
4.3. Influence of Structural Parameters of the OOL on Fatigue Damage
The outer diameter, wall thickness, and coverage area of buoyancy material are important structural parameters of the OOL. This section will analyze the influence of these three parameters on fatigue damage.
4.3.1. Influence of Coverage Length of Buoyancy Material
In order to study the influence of coverage length (
L) of buoyancy material on the fatigue damage of the OOL, three working conditions (
L = 380 m,
L = 420 m, and
L = 460 m) are selected for coupling calculation and fatigue damage analysis. The fatigue damage calculation results are shown in
Figure 14, and the corresponding spatial distribution shapes of the OOL with different coverage lengths of buoyancy material are shown in
Figure 15. The wave load used in this calculation is the 0° wind wave, with a significant wave height of 1.5 m and a period of 5.5 s. The calculation time is 2000 h.
Combined with
Figure 14 and
Figure 15, it can be found that the coverage length of buoyant material can change the spatial distribution shape of the OOL and then change the magnitude and distribution of tension. With the increase in the coverage length of buoyant material, the tension of the OOL will decrease, and the curvature of the overall structure will increase, then the fatigue damage of the structure will increase. Moreover, the variation in fatigue damage is mainly concentrated in the parts near the CALM buoy. Therefore, in order to reduce the fatigue damage caused by wave load, it is suggested that on the premise of meeting the structural strength requirements, the coverage length of buoyancy material can be reduced appropriately.
4.3.2. Influence of Outer Diameter
In order to study the influence of the outer diameter of the OOL on the fatigue damage, keep the wall thickness unchanged in the calculation, and take the outer diameter of
D = 513 mm,
D = 508 mm, and
D = 503 mm, respectively, for coupling calculation and fatigue analysis. The results are shown in
Figure 16. The wave parameters used in the calculation are the same as those in
Section 4.3.1.
As shown in
Figure 16, with the increase in the outer diameter of the OOL (the wall thickness remains unchanged), the fatigue damage of the pipeline shows an increasing trend under the same wave conditions. The variation in fatigue damage caused by the change in outer diameter is mainly concentrated in the area near the CALM buoy, however, it is not very obvious.
4.3.3. Influence of Wall Thickness
In this part, the effect of wall thickness (
WT) on fatigue damage is studied. Keep the outer diameter unchanged in the calculation, and take the wall thickness of
WT = 17 mm,
WT = 18 mm, and
WT = 21 mm, respectively, for coupling calculation and fatigue analysis. The results are shown in
Figure 17.
As shown in
Figure 17, with the increase in the wall thickness of the OOL (the outer diameter remains unchanged), the fatigue damage of the pipeline will be reduced, and the variation in fatigue damage caused by the change in wall thickness is mainly concentrated in the area near the CALM buoy. In addition, compared with the effect of outer diameter (as shown in
Figure 16), the effect of wall thickness on fatigue damage is much more significant.
Considering that the increase in wall thickness will also increase the weight of the OOL and then change its spatial distribution shape, now, the influence of wall thickness (as shown in
Figure 17) and the influence of outer diameter (as shown in
Figure 16) are compared. When the wall thickness increases from 17 mm to 18 mm, the total weight increases about 5%, and the fatigue damage can be reduced by up to 40%. When the outer diameter increases from 503 mm to 513 mm, the total weight increases about 2%, and the fatigue damage of the OOL shows an increasing trend rather than decreasing. It can be inferred that for the FPSO OOL, the influence of wall thickness itself on fatigue damage is much more significant than that of weight and spatial distribution shape variation in the OOL caused by the change in wall thickness. Therefore, in order to reduce the fatigue damage caused by the wave load, it is feasible to increase the wall thickness properly.
4.4. Influence of the Shuttle Tanker
During the service period of the offloading system, about 20–40% of the time, it is moored by the shuttle tanker for oil transportation [
25]. In order to study the effect of the shuttle tanker on fatigue damage of the OOL, coupling calculation and fatigue analysis of two conditions (with shuttle tanker and without shuttle tanker) are carried out, and the results are shown in
Figure 18. The schematic of the shuttle can be seen in
Figure 2, and the wave parameters used in the calculation are the same as those in
Section 4.3.1.
As shown in
Figure 18, although shuttle tankers are not directly in contact with the OOL, by influencing the motion state of the CALM buoy, the fatigue damage of the OOL will increase to about 1.5 times. The object of influence is mainly the part of the OOL near the CALM buoy, and the influence on other parts can be basically ignored.
4.5. Influence of the Flow-Induced Vibration
Vortex-induced vibration is also one of the important factors leading to fatigue damage of offshore pipelines [
26]. In this part, the effect of vortex-induced vibration on the fatigue damage of the OOL will be studied. In this case, the current load is applied. The flow velocity is in the form of shear flow, and the parameters are defined according to the frequently occurring sea conditions in West Africa, as shown in
Figure 19. The vortex force is calculated by the vortex tracking model based on the software Orcaflex, in which the value of the lift coefficient (Cl) adopts the table of the lift coefficient compiled by Norberg [
27].
In order to study the effect of vortex-induced vibration on fatigue damage, the fatigue damage of the OOL considering vortex-induced vibration (
Cl =
Cl (Re)) and not considering vortex-induced vibration (
Cl = 0) are calculated. The fatigue damage distribution of the OOL when the current direction is 45° is shown in
Figure 20.
As shown in
Figure 20, under the action of current loads, the fatigue damage of the OOL mainly occurs near the CALM buoy. Although the fatigue damage curve has local maximum values near the FPSO and the segment covered by buoyancy material, it is far less than the fatigue damage near the CALM buoy. In addition, the vortex force caused by vortex-induced vibration has increased the fatigue damage in the middle part of the OOL by nearly 5–10 times, but it has little influence on the fatigue damage of the parts near the CALM buoy.
The fatigue damage distribution of the OOL when the current load direction is 90° is shown in
Figure 21. Comparing
Figure 20 and
Figure 21, it can be found that the fatigue damage distribution along the OOL is different under different flow load directions, but the vortex-induced vibration can still increase the fatigue damage of the OOL by 5–10 times in the middle part. Meanwhile, the fatigue damage of the parts of the OOL near the CALM buoy is still much larger than other parts, and the effect of vortex-induced vibration is very small in this part. In addition, in the parts near the CALM buoy, the fatigue damage caused by 90° current load is slightly less than that caused by 45° current load.
According to the above analysis, it can be concluded that the fatigue damage caused by current load to the OOL is mainly concentrated in the part near the CALM buoy, and the flow load along the axial direction of the OOL can cause slightly larger fatigue damage. In addition, the effect of vortex-induced vibration is mainly concentrated in the middle part of the OOL, which can increase its fatigue damage by 5–10 times.
5. Conclusions
In this paper, fatigue analysis of oil offloading lines in the FPSO CALM buoy offloading system under wave and current loads in West Africa is carried out. The following conclusions can be drawn:
The accumulated fatigue damage to the OOL is the largest near the area connecting the CALM buoy. Additionally, the section covered by buoyancy material has local maximum accumulated fatigue damage.
Tension variation amplitude is a major factor for fatigue damage. Additionally, constrained boundary conditions may increase fatigue damage.
0° waves cause greater fatigue damage. The fatigue damage of the OOL near the CALM buoy is mainly caused by the wind wave, and the fatigue damage of sections covered by buoyancy material is mainly caused by the swell.
Appropriately reducing the cover length of buoyancy material and increasing the wall thickness can reduce fatigue damage. The outer diameter (with thickness unchanged) has little effect on the fatigue damage of the OOL.
The effect of the shuttle tanker can increase the fatigue damage of the OOL near the CALM buoy by about 1.5 times, and the effect of vortex-induced vibration is mainly concentrated in the middle part of the OOL, which can increase fatigue damage by 5–10 times.
The present study reveals the influence of various factors on the oil offloading line’s fatigue damage, such as motion and tension distribution, wave type, wave direction, coverage length of buoyancy material, outer diameter, wall thickness, shuttle tanker, and current loads, which can provide guidance for the fatigue design and optimization of pipeline layout and parameters in the FPSO oil offloading system. Moreover, the conclusions of this work can also serve as a reference for similar fatigue design in marine renewable energy technologies [
28] or other offshore engineering cases, such as risers, mooring lines, and offshore cables.