The Optimization of a Subsea Pipeline Installation Configuration Using a Genetic Algorithm

: The most commonly used subsea pipeline installation method is the S-Lay method. A very important and complex task in an S-Lay installation engineering analysis is to find the optimal pipelay vessel installation configuration for every distinctive pipeline route section. Installation loads in the pipeline are very sensitive to small changes in the configuration of the pipeline supports during laying and other influential parameters, such as the tensioner force, stinger angle, trim and draft of the pipelay vessel. Therefore, the process of an engineering installation analysis is very demanding, and there is a need for an automated optimization process. For that purpose, installation engineering methodology criteria and requirements are formalized into a nonlinear optimization problem with mixed continuous and discrete variables. A special tailored multi-objective genetic algorithm is developed that can be adjusted to any desired combination of criteria and offshore standards’ requirements. The optimization algorithm is applied to the representative test cases. The optimization procedure efficiency and quality of the achieved solution prove that the developed genetic algorithm operators and the whole optimization approach are adequate for the presented application.


Introduction
Subsea pipelines are an important means of continuous oil and gas transport over long distances, with great economic and geopolitical significance.The safety of the hydrocarbon supply has paramount importance nowadays.Consequently, the vulnerability and resilience of subsea pipelines to the various factors affecting their integrity and operability are receiving increased attention and attracting significant research focus [1,2].Subsea pipeline resilience can be defined as the pipeline's capacity to withstand manmade and natural hazards, and this resilience plays a significant role in the design process [3].These hazards could include accidental pulling with fishing gear, damages caused by terrorist activity, earthquakes, extreme climate and environmental conditions and unforeseen installation or operational loads [4][5][6].The increasing demand and growth of energy consumption affect the design and installation of ever-longer subsea pipelines at increasing depths.In such conditions, the requirements of the safety, reliability and integrity of subsea pipelines become more complex.
Subsea pipeline installation is a rather complex, demanding and expensive offshore operation.Related engineering projects are driven by contradictory requirements whose optimal combination should be found.The most important are the minimization of pipeline installation loads and fatigue; minimization of residual loads in the installed pipeline; and demand for as-low-as-possible installation operation costs.Residual tension in the installed subsea pipeline is directly proportional to the applied tensioner force during installation.Minimization of residual tension is an important factor from the standpoint of increasing pipeline resilience to hazards like accidental pulling with fishing gear or the operational hazard of damage due to vortex-induced vibrations (VIVs) in the free-span areas.That demand is, at the same time, in contradiction with the aim of minimizing installation loads, because a higher applied tensioner force decreases installation loads and vice versa.These complex requirements raise a need for an efficient optimization method of finding subsea pipeline installation configuration parameters.
A pipelay installation engineering analysis plays a pivotal role during installation project preparation.The process of laying subsea pipelines on the seabed is a very demanding operation dependent on environmental conditions.It involves careful planning and a detailed analysis in the design phase.The most commonly used subsea pipeline installation method is the S-Lay method.It is performed with S-Lay pipelay vessels; during the procedure, the pipe is controlled to bend on the seabed in the shape of the letter S (Figure 1).During the S-Lay installation, the pipeline is connected by welding on the deck of the laying ship (tubular) and is continuously lowered from the deck all the way to the seabed.The pipeline consists of a series of individual pipe elements (pipe joints) with a standard length of 12.2 m or 24.4 m, the ends of which are connected on the production line.The pipeline production line (firing line) usually contains a pipe beveling station, a pipe centering station, a few welding stations, a nondestructive testing station and at minimum one tensioning device (tensioner).The average pipeline laying capacity for most pipelay vessels is 4.5 km/day, although some vessels can set up to 7 km/day.operational hazard of damage due to vortex-induced vibrations (VIVs) in the free-span areas.That demand is, at the same time, in contradiction with the aim of minimizing installation loads, because a higher applied tensioner force decreases installation loads and vice versa.These complex requirements raise a need for an efficient optimization method of finding subsea pipeline installation configuration parameters.
A pipelay installation engineering analysis plays a pivotal role during installation project preparation.The process of laying subsea pipelines on the seabed is a very demanding operation dependent on environmental conditions.It involves careful planning and a detailed analysis in the design phase.The most commonly used subsea pipeline installation method is the S-Lay method.It is performed with S-Lay pipelay vessels; during the procedure, the pipe is controlled to bend on the seabed in the shape of the letter S (Figure 1).During the S-Lay installation, the pipeline is connected by welding on the deck of the laying ship (tubular) and is continuously lowered from the deck all the way to the seabed.The pipeline consists of a series of individual pipe elements (pipe joints) with a standard length of 12.2 m or 24.4 m, the ends of which are connected on the production line.The pipeline production line (firing line) usually contains a pipe beveling station, a pipe centering station, a few welding stations, a nondestructive testing station and at minimum one tensioning device (tensioner).The average pipeline laying capacity for most pipelay vessels is 4.5 km/day, although some vessels can set up to 7 km/day.On the pipelay vessel's deck, the pipeline rests on a series of discontinuous supports, which can be rollers or track guides.Except in extremely shallow waters, as a rule, a long fixed ramp or an immersed floating strut is used (stinger), which serves to support the pipeline behind the stern of the ship as it descends into the water.A stinger is a lattice structure mounted on the ship's stern that can be attached to the vessel via a fixed or articulated link.The floating stinger is connected by a hinge and can freely float in the water.In floating freely, the stinger supports the pipe with its buoyancy force and opposes its weight.The stinger's buoyancy force, and hence its angle, can be controlled by regulating the amount of water in the stinger ballast tanks.Similar to the deck, cylindrical supports on the stinger provide a platform for the pipeline to rest upon.
After the point of separation from the stinger, the pipe extends in a long free span all the way to the point of contact with the seabed, called the touch-down point (TDP).The free span length in the S-Lay procedure typically ranges between two and three times the water depth.
To reduce stresses in the pipeline to acceptable levels, one or more tensioners are employed, transmitting a substantial tensile force to the pipeline during assembly.This tensioning action straightens the pipeline, reducing its curvature and extending the free span length.When laying in deep waters, tensile force is crucial to minimize the angle at which the pipe detaches from the stinger (Figure 1) and ensure the pipe's separation from On the pipelay vessel's deck, the pipeline rests on a series of discontinuous supports, which can be rollers or track guides.Except in extremely shallow waters, as a rule, a long fixed ramp or an immersed floating strut is used (stinger), which serves to support the pipeline behind the stern of the ship as it descends into the water.A stinger is a lattice structure mounted on the ship's stern that can be attached to the vessel via a fixed or articulated link.The floating stinger is connected by a hinge and can freely float in the water.In floating freely, the stinger supports the pipe with its buoyancy force and opposes its weight.The stinger's buoyancy force, and hence its angle, can be controlled by regulating the amount of water in the stinger ballast tanks.Similar to the deck, cylindrical supports on the stinger provide a platform for the pipeline to rest upon.
After the point of separation from the stinger, the pipe extends in a long free span all the way to the point of contact with the seabed, called the touch-down point (TDP).The free span length in the S-Lay procedure typically ranges between two and three times the water depth.
To reduce stresses in the pipeline to acceptable levels, one or more tensioners are employed, transmitting a substantial tensile force to the pipeline during assembly.This tensioning action straightens the pipeline, reducing its curvature and extending the free span length.When laying in deep waters, tensile force is crucial to minimize the angle at which the pipe detaches from the stinger (Figure 1) and ensure the pipe's separation from the end of the stinger, i.e., the last support.The maximum capacity of the tensioners ranges from 50 to 750 ton-force (tf).
During pipeline installation, pipelay vessels maintain their position using a radial anchor system, typically comprising four to twelve anchor lines with winches and anchors arranged radially around the vessel.As new pipe segments are added to the assembly line, the ship advances forward by the length of the added segment.The forward movement is achieved by shortening (hauling in) anchor lines on the bow and releasing (paying out) anchor lines at the stern of the vessel.After advancing several hundred meters, it becomes necessary to reposition the anchors, which is accomplished by specially equipped boats known as anchor tugs.The anchor lines counteract the tensile force in the pipeline, keeping the ship in the desired position and mitigating environmental forces such as waves, currents and wind.Fourth-generation pipeline ships employ dynamic propulsors and a dynamic positioning system, which effectively replaces the role of anchor lines.Dynamic positioning enables pipeline installation at deeper depths where traditional anchor lines cannot be used.
Pipeline loads are highly sensitive to minor alterations in the configuration of the pipeline supports during installation.Based on the analysis results, the engineer often needs to position more than ten supports.Additionally, other critical parameters such as tensioner force, the angles of one or more stingers, pipelay vessel trim and draft must be varied.Consequently, the process of an engineering installation analysis is quite demanding, necessitating the development of an automated optimization process.
The optimization of subsea pipeline installation analysis parameters has been addressed in a limited number of research papers.The first work on this topic was Daley's paper [7], describing the graphical method for determining the optimal tensioner force based on the requirement of minimizing stinger length or on the requirement to achieve the maximum radius of curvature of the stinger supports.Maier et al. [8] optimized the articulated stinger geometry (a stinger composed of multiple floating parts jointly connected to each other) using the nonlinear programming method.Bhavikatti et al. [9] developed a method for optimizing tensioner force and free span length with the goal of minimizing the maximum bending moment in the pipeline's "S" curve.An enhanced sequential linear programming moving boundary method was utilized for optimization.Zhu and Cheung [10] presented an analytical method for analyzing a simplified S-Lay installation model with an articulated stinger, employing a singular perturbation technique.The derived analytical solution was employed to determine the optimal buoyancy combination for articulated stinger elements, albeit using an unspecified nonlinear programming deterministic method.More recent work by Ivić et al. [11] applied particle swarm multi-objective optimization to the S-Lay method installation configuration, including buoyancy tank distribution, for the installation of heavier pipelines.
In the context of applying genetic algorithms (GAs) to structural optimization problems of a similar nature, Wang and Chen [12] employed GAs to identify optimal locations for elastic vibrating rod supports, aiming to minimize vibration.Chiba et al. [13] optimized the arrangement of supports for a pipe system subjected to dynamic loads.Tabakov in [14] presents an overview of multi-criteria optimization of layered structures using GA.Vieira et al. [15] applied GA to optimize the configuration of vertical submarine pipes (risers) and in [16] they compared this application of GA to other optimization algorithms inspired by biological processes.Gantovnik [17] provides an overview of GAs adapted for the optimization of composite structures.
A significant number of studies have explored the application of GA optimization in offshore engineering.Shafieefar and Rezvani [18] demonstrate the optimization of anchor lines for floating oil platforms utilizing GA, while Boulougouris and Papanikolaou [19] apply GA to multi-criteria optimization of floating LNG terminals.Clauss and Birk investigated the optimization of the hydrodynamic form of general large marine structures using the tangent search method in [20,21], and then compared the efficiency of optimization with the sequential quadratic programming deterministic method on a similar problem.In [22], Lee and Clauss present a method of automated design of the shape of a floating marine structure with the objective of minimizing motion induced by sea waves.They performed shape optimization using the adaptive simulated annealing method.
Subsea pipeline design is inextricably linked to and dependent upon subsea pipeline installation engineering.A series of studies have been conducted on the topic of subsea pipeline route optimization utilizing GA and the development of associated software tools-de Lima Jr. et al. [23], Baioco et al. [24], de Lucena et al. [25] and Baioco et al. [26].
Among recent studies addressing offshore engineering problem optimization with GA are Kim et al. [27], who explored the design of an underwater chain trencher based on GA optimization.Liu et al. [28] employed hybrid algorithms to optimize the offshore wind turbine jacket substructure.Ghigo et al. [29] conducted platform optimization and a cost analysis in a floating offshore wind farm.Yin et al. [30] optimized static and discrete berth allocation for large-scale marine-loading problems using an iterative variable grouping GA.Zhang et al. [31] employed GA optimization to optimize configuration parameters of a new energy hybrid system.Sun Q. et al. [32] performed optimization of the anchor chain design of a catenary anchor leg mooring system.Sun H. et al. [33] optimized the number, hub height and layout of offshore wind turbines using GA.
The first research objective of this study was to systematically analyze and standardize the process of a subsea pipeline installation analysis.To achieve this goal, installation engineering methodology criteria and requirements were formalized into a nonlinear optimization problem involving both continuous and discrete variables.
The second research objective was to develop a complete optimization procedure based on a genetic algorithm (GA) that would enable the automated discovery of the optimal combination of all crucial parameters of subsea pipeline installation configuration, characteristic for the S-Lay method.Optimization using a genetic algorithm specifically designed for this problem should meet all the requirements that are set when analyzing and determining the configuration of the pipelaying procedure in practical project implementation.
A specifically designed multi-objective genetic algorithm is created that can be customized to suit any prescribed combination of criteria and offshore standards' requirements for subsea pipeline installation with the S-Lay method.The optimization algorithm is applied to the representative test cases.The effectiveness of the optimization procedure and quality of the obtained solution demonstrate that the developed genetic algorithm operators and whole optimization approach are suitable for the given optimization problem and application.

Materials and Methods
The development of multi-objective GA for S-Lay installation configuration optimization (Section 3) was founded on the following four methodological components, which are detailed in this section: • An S-lay installation method analysis; • Subsea pipeline installation analysis methodology; • Optimization problem formulation; • Genetic algorithm optimization.

S-Lay Installation Method Analysis
On the production line, where pipeline segments are connected, the pipeline layout can be considered nearly horizontal and parallel to the water surface (Figure 1).At the stern, where the pipe transitions into the water, the deck and pipe supports descend at a slope to ensure that the pipe enters the water at a defined angle.Pipe supports are typically configured to allow the pipe to bend downwards with a specific curvature radius.This applies to supports on both floating stinger and fixed ramps.Due to this curvature, the pipe adopts a specific separation angle from the stinger.Consequently, the stinger angle should be at least equal to the separation angle to prevent excessive pipe bending over the last support on the stinger, leading to excessive stresses at that point.
The "S" curve features an inflection point dividing it into two curvature regions: • Overbend-the initial section of the "S" curve, spanning the supports on board and the stinger.The pipeline curve y = f (x) is concave (y ′′ < 0); • Sagbend-the terminal section of the "S" curve extending to the seabed.The pipeline curve is convex (y ′′ > 0).
The inflection point is usually found immediately after the stinger's end, specifically after the last stinger support.At this point, the bending moment is zero and undergoes a sign change, as does the curvature, which, according to the Euler-Bernoulli beam theory, is linearly proportional to the bending moment.
In relatively shallower water depths and with larger-diameter pipes, the shape of the "S" curve of the pipeline is relatively gentle and is principally determined by the pipe's bending stiffness.The behavior of the pipe in the free span, that is, in the sagbend, resembles that of an elastically supported linear rod.Conversely, in deeper water depths and with smaller-diameter pipes that exhibit relatively low bending stiffness, the pipeline's behavior approximates that of a catenary model.The pipe's curvature is more pronounced and is predominantly determined by the tensioner's force.
Upon reaching the end of the sagbend region, the pipe comes into contact with the seabed at the TDP.If the seabed is relatively flat, the pipe can be regarded as lying on the seabed in an undeformed state, provided that residual tensile stress caused by the tensioner's force and bending near the contact point is disregarded.At a specific distance of approximately 100 m from the TDP, there is a point of apparent fixation of the pipeline, beyond which the pipeline's condition can be considered to be consistent and unaffected by the installation process.
Static stresses in the pipe within the overbend region are primarily influenced by four key parameters:

•
Common radius of curvature corresponding to the supports' configuration; • Pipe free span weight between the supports; • Bending torque local increase over supports; • Axial pipe tensile force generated by the tensioner force.These stresses are highly sensitive to changes in the height of the roller supports.The relatively short distance between the supports amplifies the impact of even small changes in support height.
The primary factors affecting static stresses in the pipe in the sagbend region are • Tensioner force; • Pipe weight in the free span; • Free span length.
Since the primary influencing factors that determine static stresses in the pipeline's overbend and sagbend regions operate in the vertical plane, the static analysis of laying pipelines can be simplified to a two-dimensional (2D) problem.A three-dimensional (3D) analysis is typically required only when laying pipelines along a winding route, where a detailed analysis of the ship's curved laying path is necessary.In practice, the pipe bends more significantly over the supports than between them, leading to a typical increase in strain over the supports and decreases in strain between them.
The most suitable model for analyzing static pipe stresses during S-Lay pipelay operation is a simplified one-dimensional homogeneous tensed rod subjected to its own weight and large deflections (geometrical nonlinearity), described using Euler-Bernoulli beam Equation (1) [52][53][54][55][56][57]: where T(s)-Tensile (axial) force component in an arbitrary pipe section; s-Curvilinear coordinate; θ-Slope of the pipe curve; E-Young modulus of elasticity; I-Moment of inertia of the pipe cross section; W(s)-weight of the pipe in an arbitrary pipe section.
Upon introducing Equation (2) for the horizontal force component H(s) in an arbitrary pipe section, where F(s)-Shear force in an arbitrary pipe section.
Nonlinear beam bending Equation (3) can be derived in a form commonly found in relevant literature [49,51,52]: Equation ( 3) is applicable to both shallow and deep waters and for small and large displacements.It is a third-order nonlinear differential equation with an unknown pipe length in the free span and an unknown sea bottom reaction force at the free end.This distinctive characteristic makes pipelaying problems more challenging to solve than classical nonlinear beam models and allows the formulation of the problem as either an initial value problem (IVP) or a boundary value problem (BVP).Equation (3) can also be expressed in terms of the function y(s) instead of θ(s).
Some earlier research works introduced the solution of Model (3) as an IVP [34,36], where initial values of θ(s) or y(s) and their derivatives are specified at one fixed pipe end, typically the tensioner end or TDP.An advantage of the IVP approach is that it eliminates the need for an additional boundary equation to determine the unknown pipe span length.However, a major drawback of the IVP approach is its inability to converge for nonlinear problems in general.This has led to the adoption of the more robust BVP approach in recent models.The BVP solution necessitates three boundary conditions (BCs) to be defined.Assuming a fixed pipelay vessel, typically at the tensioner end (s = 0), two BCs are established, zero rotation θ(0) = 0 or vertical displacement y(0) = 0, along with tensioner force T(0) = T 0 .The third BC is specified at the free pipe end, where a vertical displacement of y(L) = −D is defined, corresponding to the depth of the seabed D. In addition to the three BC, an additional boundary equation is required to determine the unknown elastic line free span length L. Furthermore, the arbitrarily supported pipe in the overbend and at the TDP end can be generally described using a supporting vertical displacement constraint function, y svdc (s).The specific form of this function depends on the model, vessel and stinger geometry.The conventional approach is to assume perpendicular contact reactions between the pipe and the supports.Various combinations of BCs can be employed to close the BVP.
Nonlinear beam bending Equation (3) can only be solved analytically with certain simplifying assumptions [42,47].Numerical methods are necessary to solve the full Equation (3) in conjunction with the general supporting vertical displacement constraint function y svdc (s).The most accurate results are obtained using the finite element method (FEM) with beam elements specifically designed to handle large displacements (geometrical nonlinearity) [36].Dedicated finite element software packages are employed for a subsea pipeline installation engineering analysis in the offshore construction industry.These packages can effectively model and analyze problems involving various components, including rigid, slender and flexible elements such as pipes, tubular members, flexible risers and cables, under static and dynamic loads.Offpipe software [58] stands as a prominent industry standard for static and dynamic pipelay analyses, having been extensively used in numerous worldwide projects over the years.It is based on the finite element method described in [40].Offpipe employs specialized finite element models of the entire pipeline system, encompassing the pipelay vessel, pipe supports, stinger and seabed.The software can handle both 2D and 3D analyses and perform a dynamic analysis for single or multiple simultaneous sea states, including regular and random waves.Other widely used general offshore construction analysis software packages for a pipelay installation analysis include Orcaflex by Orcina Ltd. [59] and Flexcom by Wood Group Kenny [60].
The same model can be analyzed and solved using general-purpose nonlinear finite element analysis (FEA) computer-aided design (CAD) software packages such as SIMULIA Abaqus by Dassault Systèmes [61], Ansys Mechanical by Ansys Inc. [62] and Autodesk Inventor Nastran [63].However, specialized software solutions specifically designed for a pipeline installation analysis are more efficient in handling the large number of parameters involved in a pipelay installation analysis and provide faster solver speeds due to their reduced model complexity.Ivić et al. [64] implemented and solved an S-Lay pipelaying model using SIMULIA Abaqus 6.11 [61] and validated the results with standard pipelay analysis software Offpipe v2.05.The static subsea pipeline installation methodology and optimization method described in the following sections are independent of the specific pipelay installation analysis software used.In this research, Offpipe software v2.07 [58] was employed as the S-Lay method static analysis solver for the implemented genetic algorithm optimization method.

Subsea Pipeline Installation Analysis Methodology
During the static installation analysis of the S-Lay pipelay method, a large number of parameters need to be optimized to determine their optimal combination that enables safe laying operation with minimal execution time, ultimately minimizing the entire operation costs.The engineering analysis of S-Lay pipeline installation typically involves evaluating static pipeline loads along the "S" curve, including tensile force, bending, pressure, contact forces on supports and contact forces on the seabed.This analysis aims to determine the optimal S-Lay static laying configuration.
The overbend curvature of the "S" curve is influenced by the stinger angle and the configuration (heights) of the supports on the pipelay vessel and stinger.The tensioner force affects the sagbend curvature of the "S" curve and the bending moment at the stinger tip.The required tensioner force is dependent on the water depth, immersed pipe weight, permissible overbend radius of curvature, pipe-stinger separation angle and permissible sagbend radius of curvature.
Key parameters of the S-lay pipelay method configuration that are adjusted during an installation analysis are depicted in Figure 2. Detailed parameter descriptions are provided in Section 2.3.The water depth and pipe weight are determined by pipeline design and therefore are not considered as optimizable parameters during an installation analysis.
A primary objective of pipelay configuration parameter optimization is to fulfill the project criteria.These criteria are established to safeguard the pipeline's integrity during the laying process and ensure that the installation equipment adheres to its operational limitations as defined for the specific project.
Design criteria for ensuring pipeline integrity are established by industry regulations.The DNV-ST-F101 standard "Submarine pipeline systems" [65] is the most widely used and relevant design criterion for pipeline installation projects, making it the basis for defining, testing and calibrating the genetic algorithm optimization.Other criteria essentially involve similar combinations of the permissible bending moment or equivalent strain and permissible deformation in the overbend and sagbend.Therefore, the flexible ability to define the project criteria enables optimization based on any of these criteria as described in the following sections.To accommodate the variations in overbend and sagbend criteria, a flexible approach has been adopted that allows setting both the permissible static bending moment and the permissible deformation in both the overbend and sagbend.A primary objective of pipelay configuration parameter optimization is to fulfill the project criteria.These criteria are established to safeguard the pipeline's integrity during the laying process and ensure that the installation equipment adheres to its operational limitations as defined for the specific project.
Design criteria for ensuring pipeline integrity are established by industry regulations.The DNV-ST-F101 standard "Submarine pipeline systems" [65] is the most widely used and relevant design criterion for pipeline installation projects, making it the basis for defining, testing and calibrating the genetic algorithm optimization.Other criteria essentially involve similar combinations of the permissible bending moment or equivalent strain and permissible deformation in the overbend and sagbend.Therefore, the flexible ability to define the project criteria enables optimization based on any of these criteria as described in the following sections.To accommodate the variations in overbend and sagbend criteria, a flexible approach has been adopted that allows setting both the permissible static bending moment and the permissible deformation in both the overbend and sagbend.
Operational limitations of installation equipment include permissible loads (reactions) on pipe supports (vessel and stinger), maximum applicable tensioner force and stinger ramp angle restrictions.
The DNV-ST-F101 standard [65] design criterion establishes distinct criteria for permissible combined loads based on the division into two regions: 1. Displacement controlled condition (DCC) region-pipeline deformation and the position of the "S" curve are primarily determined by the placement of supports and stingers, making them geometrically conditioned.The criterion prescribes a maximum allowable total strain, ε ≤ εallow, at all evaluated points (nodes) within the overbend.This calculated allowable installation strain represents unique input data for each pipeline section.The standard value is εallow = 0.002 (0.2%).While this region is often incorrectly referred to as the overbend in engineering practice, this terminology will be retained for consistency.However, the boundary between the overbend and sagbend is defined by the inflection point of the "S" curve (see Figure 1); 2. Load controlled condition (LCC) region-pipeline deformation and position of the "S" curve are primarily determined by the tensioner force; therefore, they are conditioned by the loads on the pipeline.The criterion prescribes a maximum allowable bending moment, Mb ≤ Mallow, at all evaluated points (nodes).This region is also often imprecisely labeled as the sagbend, although it encompasses the area above the inflection point.Therefore, this work will further distinguish between the sagbend(+) region above the inflection point of the "S" curve with a positive bending moment and the sagbend(−) area below the inflection point with a negative bending moment.
The boundary between the DCC and LCC regions is typically situated before the last support on the pipelay vessel or the fixed ramp preceding the articulated connection with the floating stinger.Consequently, the exact position of this boundary varies depending on the type of the stinger ramp attached, rigid or floating.For the pipelay vessel depicted Operational limitations of installation equipment include permissible loads (reactions) on pipe supports (vessel and stinger), maximum applicable tensioner force and stinger ramp angle restrictions.
The DNV-ST-F101 standard [65] design criterion establishes distinct criteria for permissible combined loads based on the division into two regions: 1.
Displacement controlled condition (DCC) region-pipeline deformation and the position of the "S" curve are primarily determined by the placement of supports and stingers, making them geometrically conditioned.The criterion prescribes a maximum allowable total strain, ε ≤ ε allow , at all evaluated points (nodes) within the overbend.This calculated allowable installation strain represents unique input data for each pipeline section.The standard value is ε allow = 0.002 (0.2%).While this region is often incorrectly referred to as the overbend in engineering practice, this terminology will be retained for consistency.However, the boundary between the overbend and sagbend is defined by the inflection point of the "S" curve (see Figure 1); 2.
Load controlled condition (LCC) region-pipeline deformation and position of the "S" curve are primarily determined by the tensioner force; therefore, they are conditioned by the loads on the pipeline.The criterion prescribes a maximum allowable bending moment, M b ≤ M allow , at all evaluated points (nodes).This region is also often imprecisely labeled as the sagbend, although it encompasses the area above the inflection point.Therefore, this work will further distinguish between the sagbend(+) region above the inflection point of the "S" curve with a positive bending moment and the sagbend(−) area below the inflection point with a negative bending moment.
The boundary between the DCC and LCC regions is typically situated before the last support on the pipelay vessel or the fixed ramp preceding the articulated connection with the floating stinger.Consequently, the exact position of this boundary varies depending on the type of the stinger ramp attached, rigid or floating.For the pipelay vessel depicted in Figure 2, if the stinger is rigid, the last support within the DCC region would be SR4.In the case of a floating stinger, the last support within the DCC region would be VR3.The bending moment allowance criterion for the LCC region is generally more conservative and is sometimes applied in the sagbend region.

Optimization Problem Formulation
The objective of static laying analysis optimization is to determine the optimal combination of parameters that define the S-Lay installation method configuration, ensuring that all project criteria are satisfied.The entire configuration can be characterized by a set of invariant (fixed) parameters and a set of variable parameters.
Invariant parameters are unique to each section along the pipeline route and establish the optimization environment as a set of variable states during the laying process.This group comprises parameters that serve as input variables for static pipelay analysis software and are derived from pipeline properties along the route: Water depth; 2.
Pipe weight in air; 5.
Pipe equivalent Young modulus of elasticity (accounting for concrete coating if applicable).
Variable parameters are those whose values need to be determined during the optimization (optimized parameters) and can be categorized as follows (Figure 2 Stinger angle (φ s ).
Optimization conditions are classified into two categories based on the requirement to fulfill them at the end of the optimization process: mandatory and additional optimization conditions.
Mandatory optimization conditions-compliance with these conditions is mandatory and serves as the foundation for establishing optimization constraint functions and defining the overall optimization process.These mandatory optimization conditions are derived from the previously described static laying analysis criteria:

•
Meeting the pipeline integrity criteria;

•
Meeting the installation equipment operational limiting criteria.
Additional optimization conditions-the optimization process should aim to fulfill these conditions to the greatest extent possible, but they are not mandatory.Most of these conditions serve as the foundation for defining an optimization goal function.When multiple goals are defined, the optimization problem becomes multi-objective and each goal (criterion) should have a corresponding goal weighting factor that reflects its relative importance compared to others.Optimization based on these conditions can be customized based on the specified weight inputs.Additional optimization conditions include

•
Minimization of maximum values for defined pipeline integrity criteria, e.g., the maximum value of bending, stress or deformation moment at any finite element node within overbend or sagbend regions; • Minimization of maximum values for defined operational limiting criteria of installation equipment, e.g., the maximum reaction force value at any support; • Uniform distribution of support reaction forces: minimizing the deviation of supports' reaction forces from their mean value while avoiding situations where supports are not in contact with the pipe; • Uniform distribution of load on supports-minimizing of load deviation at individual supports from their mean value; • Minimum required tensioner force application.
The optimization of a subsea pipeline S-Lay installation configuration falls under the category of general nonlinear optimization problems (Rao [30]).This problem is characterized by a combination of continuous and discrete variables, a nonlinear objective function and nonlinear constraint functions.Both the objective function and constraint functions are not expressible in an analytical form but rather represent the outcomes of a numerical simulation of the laying process.
The constrained nonlinear optimization problem can be formulated as follows: Find the vectors X * = (x 1 * , x 2 * , . .., x m * ) T and Y * = (y 1 * , y 2 * , . .., y n * ) T that minimize the objective function subject to the constraints and the given upper and lower bounds for the optimization variables where f (X, Y, P, R)-objective function; X = (x 1 , x 2 , . .., x m ) T , X ∈ R m -the vector of real (continuous) optimization variables; Y = (y 1 , y 2 , . .., y n ) T , Y ∈ Z n -the vector of integer (discontinuous) optimization variables; P = (p 1 , p 2 , . .., p p ) T , P ∈ R p -the vector of invariant real parameters; R = (r 1 , r 2 , . .., r q ) T , R ∈ Z q -the vector of invariant integer parameters; g k (X, Y, P, R) ≥ 0-inequality constraint functions; h l (X, Y, P, R) = 0-equality constraint functions; m-the number of real optimization variables; n-the number of integer optimization variables; r-the number of constraint functions in the form of inequality; s-the number of constraint functions in the form of equality; p-the number of invariant real parameters; q-the number of invariant integer parameters.
Each equality constraint function can be replaced by two inequality constraint functions, thereby eliminating the need for Equations ( 5) and ( 6) and replacing them with Equation (9):

Genetic Algorithm Optimization
The optimization of a subsea pipeline S-Lay installation configuration falls into the category of challenging nonlinear optimization problems.The nonlinear objective function and constraint functions are derived from numerical simulation of the S-Lay installation method, making them unavailable in analytical form.Consequently, first-and second-order derivatives, which are crucial for gradient-based optimization methods, must be estimated using approximate techniques, such as the finite difference method.This necessitates at least two pipelay simulations for each partial derivative of the objective function and each optimization parameter.
Standard nonlinear programming methods (gradient and others) would be inefficient for this problem, demanding excessive computational resources and often finding local optima near the starting point.Heuristic optimization methods are better suited for such nonlinear problems, as they employ a directed random search to find the global optimum.Genetic algorithms have demonstrated efficacy in solving such problems, often identifying globally optimal solutions with high probability.
A genetic algorithm (GA) is an algorithm that utilizes a directed heuristic, stochastic search method inspired by natural evolution (refs.[66][67][68][69][70][71][72]).Natural evolution relies on a selection process that prioritizes the survival of the most exceptional (best) individuals within a particular species.The operators of a genetic algorithm operate on individuals within a population over multiple generations, striving to progressively enhance their quality (fitness).Individuals that represent potential solutions are frequently compared to chromosomes and are represented using strings or binary numbers.Similarly to other heuristic methods (e.g., simulated annealing), a genetic algorithm should find a global minimum even when the objective function has multiple extremes, including local maxima and minima.A general genetic algorithm flowchart is depicted in Figure 3.
cient for this problem, demanding excessive computational resources and often finding local optima near the starting point.Heuristic optimization methods are better suited for such nonlinear problems, as they employ a directed random search to find the global optimum.Genetic algorithms have demonstrated efficacy in solving such problems, often identifying globally optimal solutions with high probability.
A genetic algorithm (GA) is an algorithm that utilizes a directed heuristic, stochastic search method inspired by natural evolution (refs.[66][67][68][69][70][71][72]).Natural evolution relies on a selection process that prioritizes the survival of the most exceptional (best) individuals within a particular species.The operators of a genetic algorithm operate on individuals within a population over multiple generations, striving to progressively enhance their quality (fitness).Individuals that represent potential solutions are frequently compared to chromosomes and are represented using strings or binary numbers.Similarly to other heuristic methods (e.g., simulated annealing), a genetic algorithm should find a global minimum even when the objective function has multiple extremes, including local maxima and minima.A general genetic algorithm flowchart is depicted in Figure 3.The stop condition for GA optimization can be based on the maximum number of generations, the fitness function limit, the change in the best fitness function tolerance or other criteria.The stop condition is typically defined as a combination of these criteria.The stop condition for GA optimization can be based on the maximum number of generations, the fitness function limit, the change in the best fitness function tolerance or other criteria.The stop condition is typically defined as a combination of these criteria.
A chromosome represents an individual within a genetic algorithm, encapsulating the values of optimization variables.To simplify representation, variables are not encoded in their decimal form but rather as integer values that represent discretized real values between the minimum and maximum for a specified discretization step.An optimization variable is often referred to as a gene.
A population is a set of solutions representing the current iteration step (generation) of a genetic algorithm.Population size is a crucial parameter in GAs and its accurate determination is essential.If the population size is too small, the genetic algorithm may converge prematurely; if it is too large, the genetic algorithm unnecessarily consumes computing resources, and the time required to improve the solution may be unnecessarily prolonged.Two critical concepts emerge in the evolutionary process and genetic search: population diversity and selective pressure.Both factors are significantly influenced by population size.Population size flexibility can be relatively easily implemented within a GA.
At each successive generation, a portion of the existing population is selected to create a new generation.Individual solutions are chosen through a fitness function process, favoring those with higher fitness values.Some selection methods evaluate the quality of each solution and select the best solutions based on direct comparison.Others perform assessments on a random sample of the entire population to expedite the process.
Many selection mechanisms are stochastic and are designed to select a small percentage of inferior solutions.This practice helps maintain population diversity, preventing premature convergence toward inferior solutions.Among the most renowned selection algorithms are simple selection (roulette wheel), proportional selection and tournament selection.
The most crucial role of crossover is to combine distinct chromosomes (individuals) and transmit their genes to a new population [66].For each new solution to be generated, a pair of (parental) solutions are chosen from a pre-selected set for reproduction.The creation of a "child" (or "offspring") employing the crossbreeding and mutation yields a novel solution that typically incorporates novel characteristics of its "parents".New parents are selected for each child and the cycle continues until a new population of appropriately sized solutions is created.
Mutation is essential because the population cannot encompass all possible genes and there must be a mechanism for introducing new genes that could potentially represent optimal solutions or serve as valuable stepping stones toward the optimum.In a genetic algorithm, mutation serves as a mechanism for maintaining population diversity.
Similar to crossover, mutations are not applied with every iteration.The probability of mutation is typically set to a low value, often below 0.05 for binary-encoded chromosomes.However, for integer coding and floating-point encoding, the mutation probability can be significantly higher, ranging up to 0.8.

Multi-Objective Genetic Algorithm for S-Lay Installation Configuration Optimization
A multi-objective GA is developed to achieve the goal of creating a flexible automated algorithm that adheres to the guidelines of the subsea pipeline installation analysis methodology outlined in Section 2.2.
The GA's stop condition combines two termination criteria and the algorithm terminates when either of the following criteria is met: No improvement in the best individual's fitness value for a predetermined number of generations, with a specified fitness function value tolerance; 2.
The maximum number of generations has been reached.

Chromosome Encoding
A chromosome represents an individual in the genetic algorithm and encapsulates the values of optimization variables.For simplicity, variables are not encoded in decimal form but are represented as integer values that discretize the real-valued range between the minimum and maximum for the specified discretization step.Optimization variables are commonly referred to as genes.The mathematical representation of the general form of the solution (chromosome) is a vector denoted by X = (x 1 , x 2 , . .., x m ) T .Figure 4 depicts the chromosome encoding employed for pipelay configuration optimization, corresponding to the pipelay vessel defined for the test cases presented in Section 3.

Fitness Function
Fitness represents a relative assessment (relative to the average or best chromosome) of chromosomes.The fitness function serves as a mechanism for evaluating chromosome quality.To optimize pipeline laying, a multi-objective function is employed and constrained GA optimization is implemented using the penalty function approach [73].Constraints are defined as permissible limits on maximum stresses in the pipeline (bending moment, relative stress or deformation) and/or permissible reaction forces on supports.The fitness function for constrained multi-objective optimization with the application of the penalty function is defined as where O(X)-Multi-objective goal function to be minimized; P(X)-Penalty function.

Chromosome Encoding
A chromosome represents an individual in the genetic algorithm and encapsulates the values of optimization variables.For simplicity, variables are not encoded in decimal form but are represented as integer values that discretize the real-valued range between the minimum and maximum for the specified discretization step.Optimization variables are commonly referred to as genes.The mathematical representation of the general form of the solution (chromosome) is a vector denoted by X = (x1, x2,..., xm) T .Figure 4 depicts the chromosome encoding employed for pipelay configuration optimization, corresponding to the pipelay vessel defined for the test cases presented in Section 3.

Fitness Function
Fitness represents a relative assessment (relative to the average or best chromosome) of chromosomes.The fitness function serves as a mechanism for evaluating chromosome quality.To optimize pipeline laying, a multi-objective function is employed and constrained GA optimization is implemented using the penalty function approach [73].Constraints are defined as permissible limits on maximum stresses in the pipeline (bending moment, relative stress or deformation) and/or permissible reaction forces on supports.The fitness function for constrained multi-objective optimization with the application of the penalty function is defined as where O(X)-Multi-objective goal function to be minimized; P(X)-Penalty function.
To streamline selection method implementation, the fitness function is transformed into an absolute fitness function that needs to be maximized: The nonlinear multi-objective goal function and the nonlinear penalty function are calculated based on the results of the pipeline S-Lay installation method static analysis To streamline selection method implementation, the fitness function is transformed into an absolute fitness function that needs to be maximized: The nonlinear multi-objective goal function and the nonlinear penalty function are calculated based on the results of the pipeline S-Lay installation method static analysis simulation, making them unavailable in analytical form.A flowchart outlining the genetic algorithm fitness function evaluation process is presented in Figure 5.

Objective Functions
Optimization objectives are related to the loading condition of the pipe (or support) and can be applied to the entire length of the pipe or specified for distinctive regions of the pipe (overbend, sagbend).The multi-objective goal function is formulated as a sum of weighted objective functions: where

Objective Functions
Optimization objectives are related to the loading condition of the pipe (or support) and can be applied to the entire length of the pipe or specified for distinctive regions of the pipe (overbend, sagbend).The multi-objective goal function is formulated as a sum of weighted objective functions: where N O -number of objectives; o i (X)-objective function; ow i -i-th objective function weight factor, where The weighting factor approach enables flexible customization of pipelay analysis criteria tailored to specific offshore construction standards' requirements.Unnecessary criteria can be readily disabled by setting their weighting factor as equal to zero while adjusting the remaining factors to ensure compliance with Equation (13).
Pipelay analysis criteria that can be independently configured for the overbend and sagbend include 1.
The objective function that minimizes the maximum overbend bending moment is where M b_max_OV -maximum overbend bending moment; M b_allow_OV -allowable overbend bending moment.
The objective function that minimizes maximum sagbend relative stress is where σ rel _ max_SAG -maximum sagbend relative stress defined with σ e _ max_SAG -maximum sagbend equivalent Von Mises stress; σ rel _ allow_SAG -allowable sagbend relative stress.The objective function that minimizes maximum sagbend strain is where ε max_SAG -maximum sagbend strain; ε allow_SAG -allowable sagbend strain.At least one of these criteria must be included for optimization in the overbend and one in the sagbend, and all three can be included simultaneously, resulting in a minimum of two and a maximum of six criteria.All included criteria have equal importance in either the overbend or sagbend, so they have equal weight factors.Since maximum loads typically occur at the pipe support location in the model, the distribution of weight factors between overbend and sagbend regions depends on the relative number of supports in each region.Additionally, an FT objective weight scaling factor is introduced to prioritize tensioner force minimization over other objectives when all constraints are met.The following general expression for objective function weighting factors addresses these optimization requirements for defined overbend objectives (criteria): and for defined sagbend objectives, where N S_OV -number of supports in overbend area; N S_SAG -number of supports in sagbend area; N S = N S_OV + N S_SAG -total number of supports; N CRIT_OV = 3-number of overbend objectives (criteria); N CRIT_SAG = 3-number of sagbend objectives (criteria); N CRIT = N CRIT_OV + N CRIT_SAG = 6-total number of objectives (criteria) in the overbend and sagbend; oa i -involvement (activity) of objective (criteria) in the overbend and sagbend (0-inactive| 1-active); FT = 0.1-objectives' weight scaling factor compared to the tensioner force minimization priority.
A scaling factor of 0.7 is applied to the overall fitness function value when all weighted objective functions are summed together.A scaling factor of 0.2 is applied for the weighting factor of the objective function that minimizes the deviation of the criteria on the supports (Equation ( 26)), and a scaling factor of 0.1 is used for the weighting factor of the objective function that minimizes the cumulative distance of the supports from the pipe (Equation ( 28)).
The objective of minimizing the deviation of criteria on supports aims to equalize loads on individual supports and achieve a uniform load distribution.The standard deviation σ STD_σ_rel of relative stress on all supports in the overbend and sagbend is calculated using The objective function that minimizes the deviation of the criteria on the supports is scaled to the interval [0,1] using the arctan function The corresponding weighting factor is calculated based on the previously defined objectives' weight scaling factor compared to the tensioner force minimization priority FT.
Minimizing the cumulative distance of the supports from the pipe aims to bring non-contacting supports closer to the pipe and achieve a configuration where all supports are in contact and have positive reaction forces.That is achieved using the following objective function: where d i is the distance of the i-th support from the pipe and the corresponding weighting factor is calculated based on the previously defined objective weight scaling factor and the tensioner force minimization priority FT.
The primary goal of tensioner force optimization is to apply the minimum necessary force while ensuring that all constraints are met.Initially, the primary focus is to achieve a configuration that satisfies all constraints, regardless of the tensioner force magnitude.Once all constraints are satisfied, the tensioner force minimization becomes the highest priority.When all constraints are satisfied, the quality function takes a value of f (x) < 1, and the absolute quality function takes a value of F(x) > 1.
The objective function and weighting factor defined as follows achieve this goal:

Penalty Functions
The penalty function method is widely used to transform constrained optimization problems into unconstrained problems suitable for genetic algorithm optimization [73].
A cost or penalty is associated with the fitness function evaluation for each constraint violation.The penalty function for constraint violation has the following general form: where NC represents the number of constraints.Each component of the penalty function is typically defined as a linear function with a penalty jump: where cw j -the constraint component weighting factor.In this implementation, for all constraints, it is set to one: cw j = 1, j ∈ {1, . .., NC}; PJ j -the penalty jump for a specific constraint.In this implementation, for all constraints, it is set to one: PJ j = 1, j ∈ {1, . .., NC}; C max (X)-the maximum value of an arbitrary criterion within the observed region; C allow -the allowable value of an arbitrary criterion within the observed region.
The penalty function for exceeding the allowable overbend bending moment is The penalty function for exceeding the allowable overbend relative stress is The penalty function for exceeding the allowable overbend strain is The penalty function for exceeding the allowable sagbend bending moment is The penalty function for exceeding the allowable sagbend relative stress is The penalty function for exceeding the allowable sagbend strain is The number of penalty functions for exceeding allowable support reactions is equal to the total number of supports, N S .Each support can have a unique allowable reaction force.
where R i (X)-reaction force on the i-th support; R allow,i -allowable reaction force on the i-th support.

Crossover
Due to the variability in optimized variable types, a flexible crossover procedure is selected as the most appropriate approach.The developed flexible crossover procedure combines the following standard basic crossover operators: Each basic crossover operator has a corresponding probability assigned to it.This probability determines the likelihood of the individual operator being activated when a crossover operation takes place.These three individual probabilities sum up to one.The total crossover probability defines the likelihood of a crossover operation occurring at all.

Mutation
A complex mutation operator application procedure has been developed to adapt GA behavior to the specific requirements of installation configuration optimization, ultimately achieving maximum optimization efficiency.In addition to the standard general mutation operators (strong and weak mutation), special mutation operators tailored to the specific S-Lay installation method configuration optimization problem have been developed.These operators are grounded in the geometry and physics of the S-Lay method model.Each mutation operator from the following list has a corresponding mutation operator probability assigned.Similar to the crossover operator probabilities, each mutation operator probability determines the likelihood of activating the particular operator when mutation operation takes place.Moreover, these individual probabilities sum up to one.The total mutation probability determines the likelihood of mutation occurring.The following is a list of developed mutation operators with their descriptions:

•
Strong mutation-each gene can randomly change its value within the entire range of the associated optimization variable.This corresponds to the standard GA mutation operator; • Weak mutation-each gene can randomly change its value, but only to the nearest integer value of the associated optimization variable.This type of mutation operator is well suited for shape optimization problems and should have a significantly higher probability of being applied compared to the strong mutation operator; • Support-to-pipe mutation-this mutation simultaneously modifies multiple genes associated with supports' heights, raising each support until it touches the pipe if there is a gap between them.The gap between the support and the pipe is calculated based on their vertical distance and the corresponding support height is increased to close the gap or to the maximum allowable support height if necessary; • Stinger angle mutation-this mutation combines a change in the stinger angle with adjustments to the supports, similar to the support-to-pipe mutation operator.Since a small change in the stinger angle can affect the contact of the stinger rollers with the pipe, it is necessary to correct the rollers' positions to maintain proper configuration and contact with the pipe without overstressing the pipe due to excessive support heights; • Reaction-based support mutation-this mutation adjusts multiple genes associated with support heights to achieve more uniform distribution of support reactions.It compares the reactions at all vessel or stinger supports, lowering the heights of supports with higher reactions and raising the heights of supports with lower reactions than the average.The magnitude of the adjustments is proportional to the deviation of the corresponding supports' reactions from the average value; • Tension mutation-due to the inherent geometry of S-Lay installation configuration, small changes in tensioner force can significantly impact bending loads and support reactions.This can lead to infeasible configurations and fitness functions' drops due to the constraint violation if tensioner force is adjusted without corresponding adjustments to support heights.The tension mutation operator addresses this issue by simultaneously altering tensioner force and supports' heights.When tensioner force is increased, the heights of more distant vessel and stinger supports are proportionally increased in order to maintain a consistent radius of curvature and vice versa; Bending moment-based support mutation-this operator is similar to the reactionbased support mutation operator.Multiple genes associated with support heights are adjusted to achieve a more uniform distribution of bending moments at pipe nodes located above supports; • Extreme support mutation-this operator corrects the heights of supports that have extreme values of constraint functions.The support height is adjusted toward the heights of its two neighboring supports.

Results and Discussion
The developed multi-objective genetic algorithm was applied and tested on two distinctive and representative test cases: 1.
A 32-inch outer diameter (OD) pipeline with concrete weight coating (CWC) installation at a moderate water depth of 50 m; 2.
A 16-inch pipeline without CWC installation at a higher water depth of 400 m.
Both test cases adhere to prescribed mixed pipe integrity criteria based on the DNVGL-ST-F101 standard [65].In the overbend region, the DCC-combined loading criterion with its defined limiting total strain was applied, while in the sagbend region, the LCC-combined loading criterion with its defined limiting bending moment was applied.Both test cases utilize the same pipelay vessel and stinger (Section 4.1) and employ the same GA optimization parameters (Section 4.2).
Due to the stochastic nature of GA optimization, each optimization run on the same test case will exhibit a unique course but should converge to optimal values for the main optimized parameters-the tensioner force, maximum bending moment and strain.A key characteristic of S-Lay installation configuration optimization is that there exist multiple combinations of the support heights, stinger angle and pipelay vessel trim and draft that can produce the same pipe string geometry and installation loads for the same tensioner force.To assess the effectiveness of the developed GA, both test case optimizations were conducted five times and primary optimization results are presented for each run, along with the mean value of the observed optimization parameter.The best optimal configurations are chosen and presented based on the maximum fitness function value achieved across all optimization runs.

Test Pipelay Vessel and Stinger
Tables 1-3 summarize the pipelay vessel and stinger geometry employed for both test cases involving an S-Lay method installation configuration analysis and optimization.The pipelay vessel incorporates four adjustable deck supports and is equipped with a rigid stinger of an approximate length of 50 m with five adjustable supports.Both test cases employ the same two installation analysis criteria, as defined in [60].The allowable overbend strain limit triggers objective function Equation ( 17) and penalty function Equation (35).Similarly, the allowable sagbend bending moment activates objective function Equation ( 18) and penalty function Equation (36).Additionally, GA optimization utilizes objective function Equation (25) to minimize the deviation of the criteria at supports and objective function Equation (27) to minimize the cumulative distance between supports and the pipe.Finally, due to the tensioner force optimization, objective function Equation (29) for tensioner force minimization is employed.In total, this yields five objectives for the GA multi-objective optimization process.
Other installation criteria objectives and penalty functions described in Sections 3.2 and 3.3 are inactive (not employed).Table 6 presents five objective function weight factors applied for the GA optimization, calculated in accordance with the corresponding equations provided in Section 3.2.The sum of active objective function weight factors equals one, as per Equation (13).

Test Case 1
Test case 1 involves S-Lay installation of a 32 inch (812.8 mm)-diameter pipeline with a wall thickness of 22.2 mm at a water depth of 50 m.The pipe steel grade is X65 with an SMYS of 450 MPa.To ensure pipeline on-bottom stability, the pipe is coated with a 50-mmthick concrete weight coating.The concrete weight coating does not affect the pipe steel integrity criteria but does increase the overall pipe weight.The allowable overbend strain limit is set to 0.2%, following the concrete coating crushing limiting pipe strain criterion.The allowable sagbend bending moment is calculated according to the DNV-ST-F101 [65] LCC-combined loading criterion, assuming a 30% increase in the static bending moment due to environmental loads (the dynamic bending moment component is assumed to be 30% of the static bending moment).The pipeline parameters and pipelay static analysis criteria for test case 1 are summarized in Table 7.     Figure 7 offers insights into the changes in the best individual tensioner force value throughout the GA optimization runs.Figures 8 and 9 showcase the variations in the best individual objective functions corresponding to the combined loading criteria in the sagbend and overbend during the optimization runs.The evolution of the best individual cumulative distance of supports from the pipe objective function across five GA optimizations is presented in Figure 10.The first GA optimization run yielded the optimal solution (final best individual) configuration, achieving the highest value of the fitness function.Optimized parameters are presented in Table 8.The pipe S-curve geometry and corresponding bending moment, total strain and support reaction distributions for the optimal solution are depicted in Figures 11-14, respectively.The positions of vessel and stinger supports are indicated by red dots and labels.

Test Case 2
Test case 2 involves a distinctive example of S-Lay installation of a pipeline with an outer diameter of 16 inches (406.4 mm) and a wall thickness of 25.4 mm, carried out at a water depth of 400 m.The pipe steel grade is X65 with an SMYS of 450 MPa.The allowable overbend strain limit is set to 0.3%, in accordance with the DNV-ST-F101 guidelines [65].

Test Case 2
Test case 2 involves a distinctive example of S-Lay installation of a pipeline with an outer diameter of 16 inches (406.4 mm) and a wall thickness of 25.4 mm, carried out at a water depth of 400 meters.The pipe steel grade is X65 with an SMYS of 450 MPa.The allowable overbend strain limit is set to 0.3%, in accordance with the DNV-ST-F101 guidelines [65].The allowable sagbend bending moment is calculated based on the DNV-ST-F101 [65] LCC-combined loading criterion, assuming a 30% increase in the static bending moment due to environmental loads (the dynamic bending moment component is assumed to be 30% of the static bending moment).The pipeline parameters and pipelay static analysis criteria for test case 2 are summarized in Table 9.The best individual tensioner force value alterations during GA optimizations are depicted in Figure 16.Figures 17 and 18 showcase the variations in the best individual   The best individual tensioner force value alterations during GA optimizations are depicted in Figure 16.Figures 17 and 18 showcase the variations in the best individual objective functions corresponding to the combined loading criteria in the sagbend and overbend during the optimization runs.Figure 19 presents the evolution of the best individual cumulative distance of supports from the pipe objective function across five GA optimizations for test case 2.

Comparison with Basic Genetic Algorithm
To demonstrate the contribution and efficiency of the developed specialized mutation operators, a comparison with basic integer GA was conducted for test case 1.The basic GA configuration involved setting the probabilities of all specialized mutation operators, except the strong mutation, to zero.The strong mutation operator corresponds to the standard GA mutation operator.The strong mutation operator randomly alters every gene value within the entire range of the corresponding optimization variable.The total mutation probability and the corresponding strong mutation probability were set to 0.65, which falls within the typical range for integer GAs.Five GA optimization runs were executed for test case 1 using the basic GA configuration. Figure 24 compares the mean fitness values across multiple runs for the developed specialized GA runs (as shown in Fig-

Comparison with Basic Genetic Algorithm
To demonstrate the contribution and efficiency of the developed specialized mutation operators, a comparison with basic integer GA was conducted for test case 1.The basic GA configuration involved setting the probabilities of all specialized mutation operators, except the strong mutation, to zero.The strong mutation operator corresponds to the standard GA mutation operator.The strong mutation operator randomly alters every gene value within the entire range of the corresponding optimization variable.The total mutation probability and the corresponding strong mutation probability were set to 0.65, which falls within the typical range for integer GAs.Five GA optimization runs were executed for test case 1 using the basic GA configuration. Figure 24 compares the mean fitness values across multiple runs for the developed specialized GA runs (as shown in Figure 6) with the mean fitness value during basic GA optimization for 200 generations.

Discussion
Two test cases were selected to encompass the variability of two crucial parameters that define subsea pipeline installation scenarios: water depth and pipe OD.The first test case involves a higher pipe OD (32 inches) installed at a moderate water depth of 50 m.The second test case represents a 16-inch pipeline installation at an increased water depth of 400 m.The test pipelay vessel possesses a maximum tensioner capacity of 150 tf, four adjustable supports on the vessel deck and a 50-m rigid stinger with five adjustable supports attached.
Five GA optimization runs were conducted for both test cases over 200 generations.The progression of all GA optimization runs is illustrated in Figures 6 and 15.All five optimization runs for each test case successfully converged to optimal solutions with nearly identical fitness function values for the best individuals.Consequently, the key optimized parameters (tensioner force, maximum sagbend bending moment, maximum overbend strain) also converged to consistent values (Figures 7-9 and 16-18).These outcomes demonstrate the efficacy and validity of the developed specialized GA.The most crucial optimized parameter, tensioner force, maintained a constant value of 25 tf for all optimization runs conducted for test case 1 (Figure 7). Figure 16 demonstrates that an optimal tensioner force of 90 tf is reached for all five optimization runs executed for test case 2.
Test case 1 GA optimization runs shown in Figure 6 demonstrate that the best individual satisfying all constraints and with a fitness value greater than 1 was obtained within the first few generations.This can be attributed to the moderate water depth of 50 m and the substantial available tensioner capacity of 150 tf, well beyond the minimum required tension of 25 tf attained at the end of the optimization process.Figure 15 indicates that a greater number of generations is required to obtain the best individual that satisfies

Discussion
Two test cases were selected to encompass the variability of two crucial parameters that define subsea pipeline installation scenarios: water depth and pipe OD.The first test case involves a higher pipe OD (32 inches) installed at a moderate water depth of 50 m.The second test case represents a 16-inch pipeline installation at an increased water depth of 400 m.The test pipelay vessel possesses a maximum tensioner capacity of 150 tf, four adjustable supports on the vessel deck and a 50-m rigid stinger with five adjustable supports attached.
Five GA optimization runs were conducted for both test cases over 200 generations.The progression of all GA optimization runs is illustrated in Figures 6 and 15.All five optimization runs for each test case successfully converged to optimal solutions with nearly identical fitness function values for the best individuals.Consequently, the key optimized parameters (tensioner force, maximum sagbend bending moment, maximum overbend strain) also converged to consistent values (Figures 7-9 and Figures 16-18).These outcomes demonstrate the efficacy and validity of the developed specialized GA.The most crucial optimized parameter, tensioner force, maintained a constant value of 25 tf for all optimization runs conducted for test case 1 (Figure 7). Figure 16 demonstrates that an optimal tensioner force of 90 tf is reached for all five optimization runs executed for test case 2.
Test case 1 GA optimization runs shown in Figure 6 demonstrate that the best individual satisfying all constraints and with a fitness value greater than 1 was obtained within the first few generations.This can be attributed to the moderate water depth of 50 m and the substantial available tensioner capacity of 150 tf, well beyond the minimum required tension of 25 tf attained at the end of the optimization process.Figure 15 indicates that a greater number of generations is required to obtain the best individual that satisfies all constraints and surpasses a fitness value of 1 in test case 2 compared to test case 1.This is primarily due to the increased water depth, which necessitates a higher optimal tension of 90 tf for that configuration, approaching the maximum tensioner capacity.
An analysis of Figure 6 suggests that all test case 1 optimization runs attained fitness values closely approaching their final values after 125 generations, suggesting that the set number of 200 generations is excessive for this case.The optimization runs for test case 2, as depicted in Figure 15, reveal that some runs could potentially be further enhanced if the optimization process was extended for more generations, while few converged to the maximum fitness values after 100 generations.Test case 1 at a 50-m water depth presents a more straightforward optimization due to the available tensioner capacity, leading to populations saturated with higher-fitness individuals early on.However, the optimization of test case 2 at a higher water depth, requiring a greater tensioner capacity, generates populations with a substantial number of individuals that fail to meet constraints and exhibit poor fitness values.This observation underscores potential benefits of implementing a dynamic GA with adaptable mutation operator probabilities.
The evolution of optimized parameters during optimization (Figures 7-10 for test case 1 and Figures 16-19 for test case 2) demonstrates that all optimization constraints are met relatively promptly and objective functions of optimized parameters exhibit consistent improvement after constraints are satisfied.
The final optimal solutions (Figures 11-14 for test case 1 and Figures 20-23 for test case 2) conclusively demonstrate that the optimization objectives have been achieved and that the best optimal configurations fulfill all specified installation analysis goals.
A comparison with the basic integer GA is conducted for test case 1 and detailed in Section 4.4.Figure 24 clearly illustrates the superiority of the developed specialized GA over the basic GA, which employs only the strong mutation operator.

Conclusions
The crucial and intricate task of a subsea pipeline S-Lay method installation engineering analysis involves identifying the optimal pipelay vessel installation configuration for each distinct pipeline route section.Installation loads in the pipeline are highly sensitive to subtle alterations in the configuration of pipeline supports during pipelay, along with other influential factors like the tensioner force, stinger angle, trim and draft of the pipelay vessel.Consequently, the process of an engineering installation analysis poses a significant challenge, necessitating an automated optimization procedure.
The initial research objective focused on analyzing and systematizing the subsea pipeline installation analysis procedure.By establishing clear procedures, goals and methodology, the conditions that must be met by the optimization process were determined, leading to the formulation of the optimization problem.This optimization problem falls under the category of constrained nonlinear optimization problems involving a combination of continuous and discrete variables.Both the objective and constraint functions cannot be represented analytically but rather derive from the results of static pipelay analysis simulations.Based on these simulation outcomes, the values of individual components of the objective function are calculated and the violations of constraint functions are evaluated.
The second research objective was to develop a comprehensive optimization procedure employing a genetic algorithm that could automatically identify the optimal combination of all critical parameters governing subsea pipeline installation configuration specific to the S-Lay method.A specialized tailored multi-objective GA is designed to be adaptable to various sets of optimization criteria and offshore standards' requirements.This genetic algorithm was specifically tailored to the intricate problem of pipeline installation using the S-Lay method, ensuring compliance with all requirements established during the analysis and formulation of the S-Lay installation analysis methodology.Multi-objective constrained optimization was achieved through the implementation of a penalty function and flexible objective and constraint prioritization via weighting factors.Tailored genetic algorithm operators specially designed for this purpose enable the efficient identification of the optimal combination of all the influential parameters governing the S-Lay installation method.This study represents the first known application of genetic algorithm optimization to the task of finding optimal subsea pipeline S-Lay method configurations.Previous research papers, as mentioned in the introduction, focused on optimizing a limited number of influential laying parameters, whereas the developed GA optimization facilitates simultaneous optimization of all influencing parameters, regardless of their number or type (continuous or discontinuous).
The developed multi-objective genetic algorithm was applied to two representative test cases.These optimized test cases were characterized by adhering to prescribed mixed pipe integrity criteria as per the DNV-ST-F101 standard [65], utilizing a pipelay vessel equipped with a single rigid stinger and a concrete-coated pipe.For the overbend region, a displacement-controlled condition with a defined limiting total strain was implemented, while for the sagbend region, a load-controlled condition with a defined limiting bending moment was employed.
Multiple GA optimization runs yielded optimal solutions that met all the optimization targets.The same minimum optimal tensioner force value was consistently obtained across all optimization runs, despite the stochastic nature of GA optimization.The efficiency of the optimization procedure and the quality of the achieved solutions demonstrate the suitability of the developed genetic algorithm operators and the overall optimization approach for the presented optimization problem and application.
Further research can proceed in several promising directions.One promising direction lies in refining the employed genetic algorithm through a more in-depth examination of the influence of specific parameters and potentially incorporating adaptative parameters (adaptive genetic algorithm) and the development of novel operators.Optimization efficiency could be enhanced by employing a hybrid genetic algorithm that utilizes a deterministic local search method toward the end of the optimization process.Another avenue involves extending the overall optimization process from single-section GA to simultaneous optimization of multiple pipeline route sections, a significantly more complex challenge.

J
. Mar. Sci.Eng.2024, 12, x FOR PEER REVIEW 14 of 38 simulation, making them unavailable in analytical form.A flowchart outlining the genetic algorithm fitness function evaluation process is presented in Figure 5.

Figure 6
Figure 6 depicts the evolution of the best individual fitness function value across five GA optimization runs, alongside the mean fitness function value for all runs.

Figure 6 .
Figure 6.Test case 1: Evolution of the best individual fitness function during GA optimization.

Figure 6 .
Figure 6.Test case 1: Evolution of the best individual fitness function during GA optimization.

Figure 6 .
Figure 6.Test case 1: Evolution of the best individual fitness function during GA optimization.

Figure 7 .
Figure 7. Test case 1: Evolution of the best individual tensioner force during GA optimization.

Figure 7 .
Figure 7. Test case 1: Evolution of the best individual tensioner force during GA optimization.

Figure 8 .
Figure 8. Test case 1: Evolution of the best individual maximum sagbend bending moment objective function (Equation (18)) during GA optimization.

Figure 8 .
Figure 8. Test case 1: Evolution of the best individual maximum sagbend bending moment objective function (Equation (18)) during GA optimization.

Figure 8 .
Figure 8. Test case 1: Evolution of the best individual maximum sagbend bending moment objective function (Equation (18)) during GA optimization.

Figure 9 .
Figure 9. Test case 1: Evolution of the best individual maximum overbend strain objective function (Equation (17)) during GA optimization.

Figure 9 .
Figure 9. Test case 1: Evolution of the best individual maximum overbend strain objective function (Equation (17)) during GA optimization.

Figure 10 .
Figure 10.Test case 1: Evolution of the best individual cumulative distance of supports from pipe objective function (Equation (27)) during GA optimization.The first GA optimization run yielded the optimal solution (final best individual) configuration, achieving the highest value of the fitness function.Optimized parameters are presented in Table 8.The pipe S-curve geometry and corresponding bending moment, total strain and support reaction distributions for the optimal solution are depicted in Figures 11, 12, 13 and 14, respectively.The positions of vessel and stinger supports are indicated by red dots and labels.

Figure 10 .
Figure 10.Test case 1: Evolution of the best individual cumulative distance of supports from pipe objective function (Equation (27)) during GA optimization.

Figure 13 .
Figure 13.Test case 1 optimal solution total strain distribution.

Figure 15 illustrates
Figure 15 illustrates the evolution of the best individual fitness function value across five GA optimization runs, alongside the mean fitness function value for all runs.The best individual tensioner force value alterations during GA optimizations are depicted in Figure16.Figures17 and 18showcase the variations in the best individual objective functions corresponding to the combined loading criteria in the sagbend and overbend during the optimization runs.Figure19presents the evolution of the best individual cumulative distance of supports from the pipe objective function across five GA optimizations for test case 2.

Figure 15 .
Figure 15.Test case 2: Evolution of the best individual fitness function during GA optimization.

Figure 16 .
Figure 16.Test case 2: Evolution of the best individual tensioner force during GA optimization.

Figure 15 . 38 Figure 15 .
Figure 15.Test case 2: Evolution of the best individual fitness function during GA optimization.

Figure 16 .
Figure 16.Test case 2: Evolution of the best individual tensioner force during GA optimization.

Figure 16 .
Figure 16.Test case 2: Evolution of the best individual tensioner force during GA optimization.The third GA optimization run yielded the optimal solution (final best individual) configuration, achieving the highest value of the fitness function.Corresponding optimal solution (final best individual) configuration parameters are presented in 10.The pipe S-curve geometry and corresponding bending moment, total strain and support reaction distributions for the optimal solution are depicted in Figures 20-23, respectively.The

Figure 17 .
Figure 17.Test case 2: Evolution of the best individual maximum sagbend bending moment objective function (Equation (18)) during GA optimization.

Figure 18 .
Figure 18.Test case 2: Evolution of the best individual maximum overbend strain objective function (Equation (17)) during GA optimization.

Figure 18 .
Figure 18.Test case 2: Evolution of the best individual maximum overbend strain objective function (Equation (17)) during GA optimization.

Figure 18 .
Figure 18.Test case 2: Evolution of the best individual maximum overbend strain objective function (Equation (17)) during GA optimization.

Figure 22 .
Figure 22.Test case 2 optimal solution total strain distribution.

Figure 22 .
Figure 22.Test case 2 optimal solution total strain distribution.

Figure 22 .
Figure 22.Test case 2 optimal solution total strain distribution.

Figure 24 .
Figure 24.Test case 1 mean fitness comparison for multiple specialized and basic integer GA runs.

Figure 24 .
Figure 24.Test case 1 mean fitness comparison for multiple specialized and basic integer GA runs.

Table 1 .
Test pipelay vessel and stinger parameters.

Table 6 .
Active objective function weight factors used in GA optimization.

Table 7 .
Test case 1 pipeline parameters and pipelay static analysis criteria.

Table 8 .
Test case 1 optimal solution configuration parameters.

Table 9 .
Test case 2 pipeline parameters and pipelay static analysis criteria.Figure15illustrates the evolution of the best individual fitness function value across five GA optimization runs, alongside the mean fitness function value for all runs.

Table 9 .
Test case 2 pipeline parameters and pipelay static analysis criteria.