Performance of an Advanced Intelligent Control Strategy in a Dynamic Positioning (DP) System Applied to a Semisubmersible Drilling Platform a Dynamic Positioning

: Oil drilling and extraction platforms are currently being used in many offshore areas around the world. Whilst those operating in shallow seas are secured to the seabed, for deeper water operations, Dynamic Positioning (DP) is essential for the platforms to maintain their position within a safe zone. Operating DP requires intelligent and reliable control systems. Nearly all DP accidents have been caused by a combination of technical and human failures; however, according to the International Marine Contractors Association (IMCA) DP Incidents Analysis, DP control and thruster system failures have been the leading causes of incidents over the last ten years. This paper will investigate potential operational improvements for DP system accuracy by adding a Predictive Neural Network (PNN) control algorithm in the thruster allocation along with a nonlinear Proportional Integral derivative (PID) motion control system. A DP system’s performance on a drilling platform in oil and gas deep-water ﬁelds and subject to real weather conditions is simulated with these advanced control methods. The techniques are developed for enhancing the safety and reliability of DP operations to improve the positioning accuracy, which may allow faster response to a critical situation during DP drilling operations. The semisubmersible drilling platform’s simulation results using the PNN strategy show improved control of the platform’s positioning.


Introduction
Historically, dynamic positioning (DP) was developed as a result of increasing oil demands and the need to extract it from deeper water. The offshore oil and gas industry has become more dependent on technology development to achieve growth, with low risk, over the last four decades [1]. The DP system's basic principle is a computer-controlled system to calculate the necessary thrust to counteract the environmental disturbance forces to automatically maintain the platform's required position and heading by using its thrusters' forces. It uses information provided by reference sensors to determine the deviation between the actual position and the required position, and the forces from the wind, waves, and current, which act upon the vessel, and it then calculates the forces that the thrusters must produce to make the deviation as small as possible [2]. A DP configuration is a machine system running under a qualified operator's supervision to achieve the prime purpose of controlling surge and sway to maintain the drilling platform's DP time-domain simulators have been rapidly developed during the last decade [4,17,18], and the results are inherently accurate having been compared with experimental analysis. Studies have included the use of PID control on semisubmersibles using a wave filter to account for the thruster forbidden zones [19] and Kalman filtering for wave motion [20].
The purpose of this paper is to develop an intelligent, reliable DP control system by integrating Predictive Neural Network (PNN) thruster allocation control into an existing basic control system to perform full dynamic positioning, enhance the accuracy of the control and minimize the risks of losing its position. In drilling operations, the DP system needs a quick response to apply precise amounts of force in certain directions to maintain position and heading in harsh weather, and keep updating the environmental forces and direction data in real-time. In this study, the traditional nonlinear PID motion control methods and the Moore-Penrose pseudo-inverse thruster allocation algorithm have been applied as a first stage for comparison with the PNN strategy. The study has utilized the thrust force and the thrusters' azimuth angle as design variables for an optimization problem. The thrusters' forbidden zone and interaction were neglected and the model is implemented without a wave filter in order to investigate the PNN algorithm's capability. The main aim is to develop a framework for improvement of accuracy and reliability of DP drilling activities by gaining real-time optimal thruster allocation in terms of directions and forces while maintaining safe drilling operations. This will form the basis for developing a generic intelligent DP control framework operated by any vessel that uses the DP system. It has been concluded that adding artificial intelligence methods to the thruster allocation control is more reliable and faster than the DP operator response during drilling operations.
The paper is organized as follows: Section 2 introduces the semisubmersible platform model, describing the mathematical model that has been used. Section 3 describes the design of PNN control and demonstrates numerical simulation. Section 4 presents the simulation results in the time-domain and discussion of those results. Finally, the study is concluded in Section 5.

Semisubmersible Platform Mathematical Modeling
This study investigates a method to improve the real-time DP control accuracy and reliability while taking into account semisubmersible platform motion, environmental loads and thruster allocation using PID and PNN controls. This combination is intended to develop a framework for an intelligent DP control strategy for application to semisubmersible drilling vessels and other vessels that use DP systems. The most common method of modeling a vessel's motion is the vectorial notation from References [17,18], which has become the standard for marine control systems.

Semisubmersible Drilling Platform Model
A model of the hull of an existing semisubmersible drilling platform with a mass of 52,476 tonnes, has been built by using the Det Norske Veritas (DNV) Sesam (GeniE) software and is depicted in Figure 1. Hydrodynamic response of the platform model was generated by using DNV Seasam (HydroD) and Wave Analysis at Massachusetts Institute of Technology (WAMIT) industry-standard numerical tools [21][22][23], and in accordance with DNV standards in order to accurately to simulate the platform behavior and provides the data needed for the time-domain simulation.

Mathematical Modeling
The mathematical model describing a drilling platform's dynamics is separated into low frequency (LF) and wave frequency (WF) models. The LF motions are caused by second-order wave loads (drift forces, including viscous effects and slowly varying forces), current loads, wind loads, and the thrusters' forces. The first-order wave loads cause the WF motions. In the modeling of marine dynamic control systems, it is considered sufficient to consider the horizontal-plane dynamics [24], which has been assumed in this paper. Typically, linear damping and wave drift loads make up a significant part of the LF forces, which again is considered in this model. The viscous forces appear as a huge effect at the wave splash zone on the semisubmersible platform columns, especially in storm conditions [25]. The model is designed to compensate for WF motions and LF motions with a maneuvering model to describe the relation between the control activity, the motion-induced and seakeeping models to describe the motion due to the wave loads. In general, the first-order wave loads are more significant than the second-order wave loads. Both can be determined by employing quadratic transfer functions [26,27].
The study of the dynamic equations of marine system's motion can be split into two aspects: The kinematic equations of motion, which relate to the geometrical aspects, and the kinetic equations of motion, which correspond to the motion analysis caused by the forces [18]. The operational drilling mode requires a specific mathematical model and parameters as the primary physical properties will depend on how the platform is operated. The vessel motions in the horizontal plane are in the surge, sway, and yaw directions; these are considered in the DP control used in this study. The reference frames used in the paper, platform kinematics, and thrusters' locations [1][2][3][4][5][6][7][8] are illustrated in Figure 2. The North-East-Down (NED) coordinate system is defined with respect to an earth-fixed , reference frame. In addition, the reference-parallel , frame is earth-fixed rotated to the required heading angle , and the origin is translated to the desired , position coordinates. The body-fixed , frame is fixed to the body of the platform rotated to the platform heading and the origin platform , position coordinate. By using the generalized position , the velocity is given with respect to a body-fixed frame, and ( ) is the transformation matrix in a vectorial representation [18]. These six differential equations are lumped together into a three Degree of Freedom (DOF) reference frame, [28] as:

Mathematical Modeling
The mathematical model describing a drilling platform's dynamics is separated into low frequency (LF) and wave frequency (WF) models. The LF motions are caused by second-order wave loads (drift forces, including viscous effects and slowly varying forces), current loads, wind loads, and the thrusters' forces. The first-order wave loads cause the WF motions. In the modeling of marine dynamic control systems, it is considered sufficient to consider the horizontal-plane dynamics [24], which has been assumed in this paper. Typically, linear damping and wave drift loads make up a significant part of the LF forces, which again is considered in this model. The viscous forces appear as a huge effect at the wave splash zone on the semisubmersible platform columns, especially in storm conditions [25]. The model is designed to compensate for WF motions and LF motions with a maneuvering model to describe the relation between the control activity, the motion-induced and seakeeping models to describe the motion due to the wave loads. In general, the first-order wave loads are more significant than the second-order wave loads. Both can be determined by employing quadratic transfer functions [26,27].
The study of the dynamic equations of marine system's motion can be split into two aspects: The kinematic equations of motion, which relate to the geometrical aspects, and the kinetic equations of motion, which correspond to the motion analysis caused by the forces [18]. The operational drilling mode requires a specific mathematical model and parameters as the primary physical properties will depend on how the platform is operated. The vessel motions in the horizontal plane are in the surge, sway, and yaw directions; these are considered in the DP control used in this study. The reference frames used in the paper, platform kinematics, and thrusters' locations [1][2][3][4][5][6][7][8] are illustrated in Figure 2. The North-East-Down (NED) coordinate system is defined with respect to an earth-fixed X E , Y E reference frame. In addition, the reference-parallel X R , Y R frame is earth-fixed rotated to the required heading angle ψ d , and the origin is translated to the desired x d , y d position coordinates. The body-fixed X, Y frame is fixed to the body of the platform rotated to the platform heading ψ and the origin platform x, y position coordinate. By using the generalized position η, the velocity V is given with respect to a body-fixed frame, and R(ψ) is the transformation matrix in a vectorial representation [18]. These six differential equations are lumped together into a three Degree of Freedom (DOF) reference frame, [28] as: The rigid-body and hydrodynamic equations of motion are in general described by three differential equations, one for each degree of freedom. These three differential equations are lumped together into a vectorial equation of motion for the nonlinear unified seakeeping and maneuvering model [18]: where is the platform inertia matrix, including the added mass, and ( ) is the skew-symmetric Coriolis and centripetal matrix. The effects of sea current on the platform are divided into two parts: the potential part is formulated as ( ) which includes the Munk moments, and the viscous part [18]. The damping vector ( ) is divided into linear and nonlinear components ( ) = + ( , ). ( ) defines the restoring vector, is the wind load vector, is the first-order wave loads, defines the second-order wave loads, and represents the thruster forces. The platform inertia mass matrix , including the added mass, the linear damping , and nonlinear damping ( , ), are defined as: where m is the platform mass, is the moment of inertia about the z-axis, and , , , , and are the zero-frequency added mass in the surge, sway and yaw directions; hence, is symmetrical and positive definite. defines the strictly positive damping caused by linear wave drift and laminar skin friction damping, The rigid-body and hydrodynamic equations of motion are in general described by three differential equations, one for each degree of freedom. These three differential equations are lumped together into a vectorial equation of motion for the nonlinear unified seakeeping and maneuvering model [18]: where M is the platform inertia matrix, including the added mass, and C RB (ν)ν is the skew-symmetric Coriolis and centripetal matrix. The effects of sea current on the platform are divided into two parts: the potential part is formulated as C A (V r )V r which includes the Munk moments, and the viscous part [18]. The damping vector D(V r )V r is divided into linear and nonlinear components D(V r )V r = D L V + D NL (V r , γ r ). G(η) defines the restoring vector, τ wind is the wind load vector, τ wave1 is the first-order wave loads, τ wave2 defines the second-order wave loads, and τ thr represents the thruster forces. The platform inertia mass matrix M, including the added mass, the linear damping D L V, and nonlinear damping D NL (V r , γ r ), are defined as: where m is the platform mass, I Z is the moment of inertia about the z-axis, and X .
r are the zero-frequency added mass in the surge, sway and yaw directions; hence, M is symmetrical and positive definite. D L ν defines the strictly positive damping caused by linear wave drift and laminar skin friction damping, where X u , Y v , Y r , N v , and N r are the hydrodynamic potential damping which can be calculated by the DNV Sesam (HydroD) and WAMIT software tools. The nonlinear damping D NL (V r , γ r ) needs C cx (γ r ), C cy (γ r ) and C cψ (γ r ), which are non-dimensional current coefficients in the horizontal plane, as illustrated in Figure 3, and can be found by model tests using the DNV Sesam Simulation of Marine Operations (SIMO) software for the platform with some defined location of the origin. In Equation [6], ρ w is the water density, L pp is the platform length between perpendicular, and D is the platform drilling draft. where , , , , and N are the hydrodynamic potential damping which can be calculated by the DNV Sesam (HydroD) and WAMIT software tools. The nonlinear damping ( , ) needs ( ), ( ) and ( ), which are non-dimensional current coefficients in the horizontal plane, as illustrated in Figure 3, and can be found by model tests using the DNV Sesam Simulation of Marine Operations (SIMO) software for the platform with some defined location of the origin. In Equation [6], is the water density, is the platform length between perpendicular, and D is the platform drilling draft. Furthermore, the sea current loads typically included in the LF model define the relative velocity vector and drag angle according to: = 2( , ).
The horizontal sea current formulation in surge and sway are defined as and , noting that, in yaw, r is small, almost equal to zero, where and are the sea current speed and direction, respectively, as follows: Generally, the motions are nonlinear, but linear approximations about specific points can be calculated. To linearize the motion equations, there is a time-domain model for the linear zero speed state-space model, which has been proposed by the Cummins equation in seakeeping [29]: where ( − ) is the retardation function, V(τ) is a unit impulse, and µ (t) is an impulse response function. A state-space formulation has been developed for the potential damping in the equation [30]. Consequently, µ (t) is represented by a linear state-space model with as the state vector: Furthermore, the sea current loads typically included in the LF model define the relative velocity vector V r and drag angle γ r according to: γ r = atan2(u r , v r ).
The horizontal sea current formulation in surge and sway are defined as u c and v c , noting that, in yaw, r is small, almost equal to zero, where V c and β c are the sea current speed and direction, respectively, as follows: Generally, the motions are nonlinear, but linear approximations about specific points can be calculated. To linearize the motion equations, there is a time-domain model for the linear zero speed state-space model, which has been proposed by the Cummins equation in seakeeping [29]: where K(t − τ) is the retardation function, V(τ) is a unit impulse, and µ (t) is an impulse response function. A state-space formulation has been developed for the potential damping in the equation [30]. Consequently, µ (t) is represented by a linear state-space model with X as the state vector: .
where A r , B r , C r , and D r are constant matrices, and µ is calculated using the state-space model equations by solving the linear seakeeping motion equations in the time-domain, which includes fluid memory effects. Then, assuming linear force superposition, wave loads can be added for different speeds and sea states. The formulation of the wind velocities u w and v w are defined according to the following equations, where the total wind speed V w , and relative wind angle γ w may be simplified: The wind loads τ wind in the surge, sway, and yaw directions are defined as follows: where ρ a is the air density, L oa is the platform length overall, and A wx , A wy are the lateral and longitudinal areas of the platform freeboard projected on the xz-plane and yz-plane. C wx (γ w ), C wy (γ w ), and C wψ (γ w ) are the non-dimensional wind coefficients in the horizontal plane, which can be calculated by using DNV Sesam (SIMO), as shown in Figure 4, or found by employing semi-empirical formulae as presented in Reference [31]. where , , , and are constant matrices, and is calculated using the state-space model equations by solving the linear seakeeping motion equations in the time-domain, which includes fluid memory effects. Then, assuming linear force superposition, wave loads can be added for different speeds and sea states. The formulation of the wind velocities and are defined according to the following equations, where the total wind speed , and relative wind angle may be simplified: The wind loads in the surge, sway, and yaw directions are defined as follows: where is the air density, is the platform length overall, and , are the lateral and longitudinal areas of the platform freeboard projected on the xz-plane and yzplane.  By assuming linear force superposition, and are obtained by using the Quadratic transfer function to get the Response Amplitude Operator (RAO) forces [32], as illustrated in Figures 5 and 6. The Joint North Sea Wave Project (JONSWAP) spectra were applied for the platform wave parameters. By assuming linear force superposition, τ wave1 and τ wave2 are obtained by using the Quadratic transfer function to get the Response Amplitude Operator (RAO) forces [32], as illustrated in Figures 5 and 6. The Joint North Sea Wave Project (JONSWAP) spectra were applied for the platform wave parameters.
where τ thr defines the generalized forces generated by the thruster system, and the thrust configuration T(α) with the azimuth angles α, the control forces obtained by the thruster system, f = Ku with the magnitudes of the force produced by each thruster vector u, and K is a diagonal force coefficient matrix. T(α) can be expressed as, where the angle α i is the angle of the i-th actuator i = 1 . . . 8, determining the force direction produced in the platform body-fixed coordinate system, l xi , l yi are the locations of the thrusters on the platform by using the extended thrust vector τ c [18] defined according to where u e is the extended control input vector, and T e is the extended thrust configuration, which is constant while T i (α i ) depends on α i . u e can be solved directly from Equation (23) using a least-squares optimization method, then the azimuth control can be derived from u e by mapping the pairs using [33]: where , , , and are constant matrices, and is calculated using the state-space model equations by solving the linear seakeeping motion equations in the time-domain, which includes fluid memory effects. Then, assuming linear force superposition, wave loads can be added for different speeds and sea states. The formulation of the wind velocities and are defined according to the following equations, where the total wind speed , and relative wind angle may be simplified: The wind loads in the surge, sway, and yaw directions are defined as follows: where is the air density, is the platform length overall, and , are the lateral and longitudinal areas of the platform freeboard projected on the xz-plane and yzplane.      The azimuth angles must be calculated with the control inputs subject to amplitude and rate saturation. The solution to this problem by using Lagrange multipliers is defined as the generalized inverse: In the case where W = I, an equal weighting is given to each control force where W is a positive definite matrix. Thus, Equation (27) reduces to the Moore-Penrose pseudo-inverse.
Since f e = T * e τ c , the control input can be computed as: The locations of the thrusters, as shown in Figure 2, are presented in Table 1. The semisubmersible platform in the simulation was commanded to maintain the position and heading in the horizontal-plane. The nonlinear PID controller is modeled, with an error signal η e computed by using the reference system inputs for the semisubmersible platform. The nonlinear PID motion control receives the error signals and provides the force required to the platform. The position and speed of actual values are obtained from the platform dynamic equations. The PID control has been designed modeled using marine systems simulator (MSS) Simulink toolbox [34].
where K p , K i , and K d are the non-negative controller gains for the proportional, time integral, and time derivative of the error signal, respectively. The controller's gains used for drilling operations were obtained through tuning using the Simulink optimization toolbox as presented in Table 2.

PNN Model
PNN control is an advanced control method in which a Neural Network (NN) model of the nonlinear DP control has been used to predict future platform performance utilizing the Newton-Raphson optimization algorithm to minimize the cost function [35]. The computationally efficient derivation of PNN will minimize the cost function with a gradientbased iterative solution in real-time. The PNN provides the transformation between the required forces as input and the thruster commands as output in three DOF. The PNN determines the thruster allocation control input to optimize the platform performance over a specified time horizon. Figure 7 and Equation (30) illustrate the NN procedure in which the sum of inputs is multiplied by the weights, added to the bias terms, and then passed through an activation function F(u).
where N is defined as the number of inputs, j is the input number, W j the weights of the input j, y j is the input j and θ is the bias.
J. Mar. Sci. Eng. 2021, 9, x FOR PEER REVIEW 10 of 25 mersible platform. The nonlinear PID motion control receives the error signals and provides the force required to the platform. The position and speed of actual values are obtained from the platform dynamic equations. The PID control has been designed modeled using marine systems simulator (MSS) Simulink toolbox [34].
where , , and are the non-negative controller gains for the proportional, time integral, and time derivative of the error signal, respectively. The controller's gains used for drilling operations were obtained through tuning using the Simulink optimization toolbox as presented in Table 2.

PNN Model
PNN control is an advanced control method in which a Neural Network (NN) model of the nonlinear DP control has been used to predict future platform performance utilizing the Newton-Raphson optimization algorithm to minimize the cost function [35]. The computationally efficient derivation of PNN will minimize the cost function with a gradientbased iterative solution in real-time. The PNN provides the transformation between the required forces as input and the thruster commands as output in three DOF. The PNN determines the thruster allocation control input to optimize the platform performance over a specified time horizon. Figure 7 and Equation (30) illustrate the NN procedure in which the sum of inputs is multiplied by the weights, added to the bias terms, and then passed through an activation function F(u).
where N is defined as the number of inputs, j is the input number, the weights of the input j, is the input j and is the bias. For PNN multiple-layer feedforward networks with time-delayed line (TDL) there is one hidden layer of tan-sigmoid transfer function (tan-sigmoid neurons), including input weight (IW) and input bias . This is followed by an output layer of a linear transfer function (linear neurons) used for the fitting problem, including layer weight (LW) and For PNN multiple-layer feedforward networks with time-delayed line (TDL) there is one hidden layer of tan-sigmoid transfer function (tan-sigmoid neurons), including input weight (IW) and input bias θ 1 . This is followed by an output layer of a linear transfer function (linear neurons) used for the fitting problem, including layer weight (LW) and input bias θ 2 . Multiple layers of neurons with nonlinear transfer functions allow the network to learn nonlinear relationships between PID control input and thruster force output, as demonstrated in Figure 8. input bias . Multiple layers of neurons with nonlinear transfer functions allow the network to learn nonlinear relationships between PID control input and thruster force output, as demonstrated in Figure 8. This approach yielded a PNN structure of with 3 input neurons, 20 hidden layers, and 8 output neurons of . A numerical cost function optimization program used the predictions to specify the control signal that would decrease the next performance criterion over the specified horizon: where , , and describe the horizons over which the signal error and the control inputs are estimated. The term defines the tentative control signal that minimizes the cost function of . describes the required response from the PID motion control, and defines the neural network model response. The value represents the influence that the sum of the squares of the control augmentations has on the performance index. The first phase of the PNN is to train the NN model to define the platform thruster allocation′s forward dynamics. The prediction error between the thruster system outputs and the NN model outputs is used for training the signal. The PNN involves the NN model and an optimization block to calculate the tentative control signal's values that decrease the performance criterion, then integrate the optimal control signal to the thruster allocation. This was simulated by using the neural network predictive control toolbox in the Simulink software.
The advanced intelligent DP motion control system is designed to control the LF and WF motions of a drilling platform, which is based on five interconnected systems: (1) setposition to the platform model so that the reference system estimates the position, course and distance travelled that in turn is needed to control the vessel's course, (2) the reference system which determines the position and heading errors to maintain the set-position, (3) the PID control system which determines the necessary motion control forces and moments produced by the platform to satisfy a specific control objective, (4) the PNN allocation control system which determines the thruster action required to keep the platform in the safe drilling zone through optimizing the PID motion control signal by training the NN model using Bayesian regularization to predict the optimum performance, and (5) the thruster system that produces the forces needed to maintain the position, and thruster This approach yielded a PNN structure of τ PID with 3 input neurons, 20 hidden layers, and 8 output neurons of u PNN . A numerical cost function optimization program used the predictions to specify the control signal that would decrease the next performance criterion over the specified horizon: where N 1 , N 2 , and N u describe the horizons over which the signal error and the control inputs are estimated. The term u defines the tentative control signal that minimizes the cost function of u PNN . τ PID describes the required response from the PID motion control, and τ PNN defines the neural network model response. The ρ value represents the influence that the sum of the squares of the control augmentations has on the performance index. The first phase of the PNN is to train the NN model to define the platform thruster allocation s forward dynamics. The prediction error between the thruster system outputs and the NN model outputs is used for training the signal. The PNN involves the NN model and an optimization block to calculate the tentative control signal's values that decrease the performance criterion, then integrate the optimal control signal to the thruster allocation. This was simulated by using the neural network predictive control toolbox in the Simulink software.
The advanced intelligent DP motion control system is designed to control the LF and WF motions of a drilling platform, which is based on five interconnected systems: (1) setposition to the platform model so that the reference system estimates the position, course and distance travelled that in turn is needed to control the vessel's course, (2) the reference system which determines the position and heading errors to maintain the set-position, (3) the PID control system which determines the necessary motion control forces and moments produced by the platform to satisfy a specific control objective, (4) the PNN allocation control system which determines the thruster action required to keep the platform in the safe drilling zone through optimizing the PID motion control signal by training the NN model using Bayesian regularization to predict the optimum performance, and (5) the thruster system that produces the forces needed to maintain the position, and thruster allocation in order to calculate the optimized control forces for each thruster. The overall simulation also includes three environment modules, in which waves, current, and wind are generated, as shown in Figure 9. allocation in order to calculate the optimized control forces for each thruster. The overall simulation also includes three environment modules, in which waves, current, and wind are generated, as shown in Figure 9. The simulation study presented in this paper investigates the potential operational impacts of using the advanced intelligent control systems through three significant steps: First, the modeling of a hull of an existing semisubmersible drilling platform to calculate the hydrodynamic data and environmental coefficient data, using DNV Sesam (GeniE and HydroD), SIMO, and WAMIT software. The second step in the modeling is building a time-domain simulation with nonlinear PID control, as in primary DP systems, following Fossen′s (Unified Model) method in which maneuvering and seakeeping ship motion theories are modeled in three degrees of freedom (DoF). This has been completed by using MATLAB/Simulink software. Thirdly, the DP control has been enhanced by adding the PNN control method, and this has been developed using the Simulink Model Predictive Control Toolbox.

Results and Discussion
This study used a model that simulates the environmental loads affecting the platform, namely wind, current, and first-order and second-order waves. Platform response data is taken from the platform model's simulation with integrated hydrodynamic and hydrostatic matrices (Inertia mass, Coriolis and centripetal, damping and restoring forces). The Nonlinear PID control collects the position and heading difference data and then sends the motion control unit's input to the thruster allocation to configure the required thruster force using the Pseudo-inverse (Moore-Penrose) matrix. The platform initial position and heading are , = 0 m and = 0 degrees in the NED reference frame. The PNN control determines the best action to send to the thruster system to provide the platform′s required force. The system is evaluated by simulating DP drilling operations in oil and gas deep-water fields in real weather conditions to drilling operation limits, as shown in Table 3. The simulation study presented in this paper investigates the potential operational impacts of using the advanced intelligent control systems through three significant steps: First, the modeling of a hull of an existing semisubmersible drilling platform to calculate the hydrodynamic data and environmental coefficient data, using DNV Sesam (GeniE and HydroD), SIMO, and WAMIT software. The second step in the modeling is building a time-domain simulation with nonlinear PID control, as in primary DP systems, following Fossen s (Unified Model) method in which maneuvering and seakeeping ship motion theories are modeled in three degrees of freedom (DoF). This has been completed by using MATLAB/Simulink software. Thirdly, the DP control has been enhanced by adding the PNN control method, and this has been developed using the Simulink Model Predictive Control Toolbox.

Results and Discussion
This study used a model that simulates the environmental loads affecting the platform, namely wind, current, and first-order and second-order waves. Platform response data is taken from the platform model's simulation with integrated hydrodynamic and hydrostatic matrices (Inertia mass, Coriolis and centripetal, damping and restoring forces). The Nonlinear PID control collects the position and heading difference data and then sends the motion control unit's input to the thruster allocation to configure the required thruster force using the Pseudo-inverse (Moore-Penrose) matrix. The platform initial position and heading are x, y = 0 m and ψ = 0 degrees in the NED reference frame. The PNN control determines the best action to send to the thruster system to provide the platform s required force. The system is evaluated by simulating DP drilling operations in oil and gas deep-water fields in real weather conditions to drilling operation limits, as shown in Table 3. The advanced intelligent DP control framework developed in this paper is based on two weather condition scenarios. The use of the first weather scenarios investigates the characteristics of the PID motion control system based on comparison with other published results. The second weather case investigates the improvements achieved by adding PNN thruster allocation to the conventional PID motion control in DP drilling operations.

Conventional PID Model
Comparisons with DP time-domain simulators for motion control and thruster allocation systems have been made to demonstrate the similarity of responses between the developed framework in this study and the published literature. References [19,20] both simulated semisubmersible platforms at 45 degree environmental angle of attack (quarter sea). Figure 10 presents simulation results obtained in this study for the position and heading under the same environmental conditions and thruster power as case 3 in Reference [19], at the required position and heading x, y = 10 m and ψ = 0 degrees within a 1000 s time-frame.  The advanced intelligent DP control framework developed in this paper is based on two weather condition scenarios. The use of the first weather scenarios investigates the characteristics of the PID motion control system based on comparison with other published results. The second weather case investigates the improvements achieved by adding PNN thruster allocation to the conventional PID motion control in DP drilling operations.

Conventional PID Model
Comparisons with DP time-domain simulators for motion control and thruster allocation systems have been made to demonstrate the similarity of responses between the developed framework in this study and the published literature. References [19,20] both simulated semisubmersible platforms at 45 degree environmental angle of attack (quarter sea). Figure 10 presents simulation results obtained in this study for the position and heading under the same environmental conditions and thruster power as case 3 in Reference [19], at the required position and heading , = 10 m and = 0 degrees within a 1000 s time-frame. The studies in Reference [19] presented experimental results for a smaller semisubmersible model but using similar environmental load tools. Table 4 shows a comparison of the key control performance parameters for the results in Figure 10 and [19]. The response from the platform in this study, when compared with the response shown in Reference [19] demonstrates similar behavior during both transient and steady state periods with regard to the Y position and heading setting. In the case of the response for X position, the steady state result is similar, as expected; however, the transient response in this study shows improvements when compared with the result from Reference [19]. The studies in Reference [19] presented experimental results for a smaller semisubmersible model but using similar environmental load tools. Table 4 shows a comparison of the key control performance parameters for the results in Figure 10 and [19]. The response from the platform in this study, when compared with the response shown in Reference [19] demonstrates similar behavior during both transient and steady state periods with regard to the Y position and heading setting. In the case of the response for X position, the steady state result is similar, as expected; however, the transient response in this study shows improvements when compared with the result from Reference [19]. The improvements in performance may be attributed to the higher thruster allocation from the more aggressive PID control algorithm used in this study. Figure 11 shows the thruster force for thrusters 3 and 6 of the platform shown in Figure 2. Both this platform and the one modeled in Reference [19] have thrusters with a maximum rating of 720 kN; however, higher thruster capability is used in this study with 88.2% and 87.3% for thrusters 3 and 6 here versus 55.5% and 51.1% for thrusters 2 and 7 (the equivalent locations) in Reference [19]. The improvements in performance may be attributed to the higher thruster allocation from the more aggressive PID control algorithm used in this study. Figure 11 shows the thruster force for thrusters 3 and 6 of the platform shown in Figure 2. Both this platform and the one modeled in Reference [19] have thrusters with a maximum rating of 720 kN; however, higher thruster capability is used in this study with 88.2% and 87.3% for thrusters 3 and 6 here versus 55.5% and 51.1% for thrusters 2 and 7 (the equivalent locations) in Reference [19]. Further comparison was made of this study′s results with those of Reference [20], which reported the performance of a smaller semisubmersible model at a required position and heading = −4 m, = −3 m and = 0 degrees within a 2000 s time-frame, as shown in Figure 12. Further comparison was made of this study s results with those of Reference [20], which reported the performance of a smaller semisubmersible model at a required position and heading x = −4 m, y = −3 m and ψ = 0 degrees within a 2000 s time-frame, as shown in Figure 12. The model in Reference [20] used PID motion control with wind feedforward and a Kalman filter for wave motion instead of the conventional PID model used in this study. Table 5 shows the key control performance parameters, comparing this study with reference [20]. It can be seen that there is similar transient response in terms of position and heading displacements. It can be seen that there are differences in terms of settling time and overshoot between the two sets of data, but the model in this study reaches the peak more quickly. Since the thruster forces in Reference [20] were mainly in the range of 300 to 400 kN, the thrusters in this paper were constrained to have maximum and minimum thrust forces by using a least-square optimization method, as shown in Figure 13. The model in Reference [20] used PID motion control with wind feedforward and a Kalman filter for wave motion instead of the conventional PID model used in this study. Table 5 shows the key control performance parameters, comparing this study with reference [20]. It can be seen that there is similar transient response in terms of position and heading displacements. It can be seen that there are differences in terms of settling time and overshoot between the two sets of data, but the model in this study reaches the peak more quickly. Since the thruster forces in Reference [20] were mainly in the range of 300 to 400 kN, the thrusters in this paper were constrained to have maximum and minimum thrust forces by using a least-square optimization method, as shown in Figure 13.
As a result, the thruster capability in the two cases was very similar, with 89.2% and 88.8% for thrusters 3 and 6 compared to 86.1% and 88.1% for thrusters 2 and 7 in Reference [20]. Figures 10-13 demonstrate the position and heading keeping capabilities of the PID motion control with the pseudo-inverse thruster allocation algorithm and provide a good basis for examining the PNN thruster allocation capability.

Settling time (s)
375 470 500 400 380 450 It can be seen that there are differences in terms of settling time and overshoot between the two sets of data, but the model in this study reaches the peak more quickly. Since the thruster forces in Reference [20] were mainly in the range of 300 to 400 kN, the thrusters in this paper were constrained to have maximum and minimum thrust forces by using a least-square optimization method, as shown in Figure 13. As a result, the thruster capability in the two cases was very similar, with 89.2% and 88.8% for thrusters 3 and 6 compared to 86.1% and 88.1% for thrusters 2 and 7 in Reference [20]. Figures 10-13 demonstrate the position and heading keeping capabilities of the PID motion control with the pseudo-inverse thruster allocation algorithm and provide a good basis for examining the PNN thruster allocation capability.

Intelligent PNN Model
Recently, [36] studied a neural network approach to control the allocation of a shipshaped vessel and different types of thrusters using DP. The study proposed a new static NN strategy for the thruster allocation control using PID motion control with a wave filter. It took a similar approach to the proposed scheme for DP using PNN allocation but without prediction optimization. Research studies in this sector have been on the drill ships, meaning that more research is needed on semisubmersible platforms. In this paper, the DP simulation results show the semisubmersible platform position and heading using the developed intelligent algorithm based on a combination of PID and PNN controls.
The framework developed in this study has been employed in the simulation of the semisubmersible model under particular environmental conditions for drilling operations to investigate whether the platform could maintain its position in the drilling zone. Safety zones have been applied to assess the safety of DP drilling operation position and heading [37]. Drilling safety zones are defined as the Red region viewed as unacceptable, leading to an incident (24 m diameter), Yellow region viewed as caution and may lead to an incident (18 m diameter), and the Green region viewed as satisfactory and safe for drilling operations (12 m diameter). The semisubmersible platform in the simulation was commanded to maintain the position and heading at , = 0 m and = 10 degrees, and it was modeled as full scale in the DP simulation. The figures presented here show the semisubmersible platform position, heading performance map and the corresponding thruster force and angle efficiency. Environmental angles of attack from the north (head sea) and east (beam sea) were tested using PID and PNN controls. The north environmental angle of attack (head sea) condition is presented in Figures 14-18. vspace-12pt

Intelligent PNN Model
Recently, [36] studied a neural network approach to control the allocation of a shipshaped vessel and different types of thrusters using DP. The study proposed a new static NN strategy for the thruster allocation control using PID motion control with a wave filter. It took a similar approach to the proposed scheme for DP using PNN allocation but without prediction optimization. Research studies in this sector have been on the drill ships, meaning that more research is needed on semisubmersible platforms. In this paper, the DP simulation results show the semisubmersible platform position and heading using the developed intelligent algorithm based on a combination of PID and PNN controls.
The framework developed in this study has been employed in the simulation of the semisubmersible model under particular environmental conditions for drilling operations to investigate whether the platform could maintain its position in the drilling zone. Safety zones have been applied to assess the safety of DP drilling operation position and heading [37]. Drilling safety zones are defined as the Red region viewed as unacceptable, leading to an incident (24 m diameter), Yellow region viewed as caution and may lead to an incident (18 m diameter), and the Green region viewed as satisfactory and safe for drilling operations (12 m diameter). The semisubmersible platform in the simulation was commanded to maintain the position and heading at x, y = 0 m and ψ = 10 degrees, and it was modeled as full scale in the DP simulation. The figures presented here show the semisubmersible platform position, heading performance map and the corresponding thruster force and angle efficiency. Environmental angles of attack from the north (head sea) and east (beam sea) were tested using PID and PNN controls. The north environmental angle of attack (head sea) condition is presented in Figures 14-18.           It is clear that the PID control could not maintain the platform position within the limit of the green, safe zone, as shown in Figure 14. In addition, the heading performance is improved in terms of settling time when using the PNN as compared to the PID only control, as presented in Figure 16. These findings show that the PNN algorithm is able to stabilize the force and direction of the thrusters effectively. The results clearly show that PID control requires the thrusters to operate at or close to their maximum 700 kN rating and with large thrust direction deviations that lead to high wear and tear, as illustrated in Figures 17 and 18. By using the PNN algorithm for thruster allocation, the wear on these thrusters is significantly reduced. The impact of using the computationally efficient PNN control to minimize the cost function in the development of an intelligent DP system (PNN algorithm) can be seen in the simulation results. The dynamic response of the PID controller is not fast enough to control the nonlinear system, and it requires the use of the PNN where the gradient can be computed rapidly by the Newton-Raphson optimization algorithm. The intelligent DP system improved the platform position and heading performance by distributing the thruster maximum forces with smaller angles than in the PID case.
The east environmental angle of attack (beam sea) condition results are presented in Figures 19-23, showing comparisons of the semisubmersible platform position, heading performance map, corresponding thruster force and angle efficiency using PID and PNN controls.

case.
The east environmental angle of attack (beam sea) condition results are presented in Figures 19-23, showing comparisons of the semisubmersible platform position, heading performance map, corresponding thruster force and angle efficiency using PID and PNN controls.       From the DP simulation using PID control only with head and beam sea angles of attack, as shown in Figures 14 and 19, it is clear that the DP control has difficulty in maintaining the platform within the limit of the green, safe zone. Meanwhile, the DP simulations using PID and PNN controls with head and beam sea angles of attack are shown in From the DP simulation using PID control only with head and beam sea angles of attack, as shown in Figures 14 and 19, it is clear that the DP control has difficulty in maintaining the platform within the limit of the green, safe zone. Meanwhile, the DP simula- From the DP simulation using PID control only with head and beam sea angles of attack, as shown in Figures 14 and 19, it is clear that the DP control has difficulty in maintaining the platform within the limit of the green, safe zone. Meanwhile, the DP simulations using PID and PNN controls with head and beam sea angles of attack are shown in Figures 15 and 20. The heading performance shows a faster response and less overshoot, shorter settling time to the set heading using PNN control in Figures 16 and 21, but they are almost the same peak and rise time in the head seas. One of the advantages of the intelligent DP framework developed in this study is the rapid response of the PNN thruster allocation in allowing the DP system to easily adjust to a step change and minimize the position error. The commands provided by the PNN are sensitive to variations in the generalized force command. These simulation results show improved position accuracy when using the PNN thruster allocation but less impact the platform's heading because the head sea effect on the platform heading is the most significant. It was clear that the PID plus PNN maintained the platform in the drilling zone, giving it a higher margin of safety when compared to using PID alone. It improves the position and heading accuracy to keep the platform in the green drilling zone as shown in Table 6. The PID thruster force and angle efficiency in Figures 17, 18, 22 and 23 show high wear and tear are likely to develop with maximum forces in every direction. In addition, the corresponding PNN thruster forces and angle efficiency illustrate that the forces and directions are more stable, and less wear and tear would result in terms of forces with directions in excess of ±40 degrees minimizing the thrust loss effect. In addition, the average and the standard deviation of the thruster force and angle, as shown in Table 7, clearly demonstrates that the PID plus PNN has larger average force, but lower average angle. In addition, the standard deviation for both force and angle is lower than for PID alone.
It is clear that PID control alone cannot handle the high sea condition typically experienced by vessels operating in sea state seven or more. However, the intelligent DP system based on the PNN control can increase the accuracy and maintain safe drilling operation. The results further show that the semisubmersible platform could maintain the position and heading in the safe drilling zone under different environmental angles of attack by using PID control for motion and PNN control for organizing the thruster forces. In addition, the PID control would potentially cause significant wear and tear in terms of the thruster forces and directions, while the PNN control is able to minimize this through appropriate thruster allocation. Although drilling is not expected to continue in such extreme conditions, one of the work targets is extending the weather window in which drilling can occur. This study offers a timely input to the increasing demand for reliability in intelligent DP control technology. The system developed in the study would also allow faster and more precise response action to avoid experiencing a critical situation during DP drilling operations. The study also offers a framework for improving environmental pollution protection by minimizing the risk of drilling platform accidents instead of the continued use of the conventional control strategy that has a significant failure record.

Conclusions
In this study, a hull of an existing semisubmersible drilling platform has been modeled in accordance with certified DNV Classification Society standards using DNV Sesam (GeniE and HydroD), SIMO and WAMIT software. The study developed a framework for a generic intelligent DP control strategy for drilling vessels. The results reported in this paper confirmed improvements in the DP thruster allocation control accuracy by adding a PNN to the existing basic DP system, enhancing safety and zone keeping of DP operations in harsh weather. Additionally, the study shows that by integrating artificial intelligence into the thruster allocation control, the system becomes safer and more reliable than the conventional dependence on the DP operator to handle all drilling operations themselves. It also demonstrates that an intelligent DP control system instead of a conventional control strategy, can offer environmental pollution protection by minimizing errors that could cause a drilling platform to ground leading to potential oil spillage from a collapsed drilling well. The forces and directions predicted by the PNN allocation control are more stable even for a small angle deviation. They offer less wear and tear in terms of the forces, minimizing the thrust loss effect. It is expected that further research will be needed to fully develop thruster allocation control performance, using advanced motion control based on PNN with sliding mode and backstepping controls. Using new prediction methods in the thruster allocation control, such as nonlinear automatic regression with recessive exogenous (NARX), requires further investigation.