# A Thrust Allocation Method for Efficient Dynamic Positioning of a Semisubmersible Drilling Rig Based on the Hybrid Optimization Algorithm

Authors: Luman Zhao; Myung-Il Roh
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

A thrust allocation method was proposed based on a hybrid optimization algorithm to efficiently and dynamically position a semisubmersible drilling rig. That is, the thrust allocation was optimized to produce the generalized forces and moment required while at the same time minimizing the total power consumption under the premise that forbidden zones should be taken into account. An optimization problem was mathematically formulated to provide the optimal thrust allocation by introducing the corresponding design variables, objective function, and constraints. A hybrid optimization algorithm consisting of a genetic algorithm and a sequential quadratic programming (SQP) algorithm was selected and used to solve this problem. The proposed method was evaluated by applying it to a thrust allocation problem for a semisubmersible drilling rig. The results indicate that the proposed method can be used as part of a cost-effective strategy for thrust allocation of the rig.


---

## Body

## 1. Introduction

Research Background. A dynamic positioning (DP) system is especially useful to deep-sea working vessel, such as drilling vessel and floating crane. The DP system can maintain the position (with a fixed location or on a predetermined track) and the heading of a working vessel exclusively through the use of active thrusters and propellers. To keep the position means to maintain the desired position within the normal excursion from the desired position and heading in the horizontal plane. For an offshore vessel, which is a semisubmersible drilling rig illustrated in Figure 1 in this case, the DP system generally consists of a position and heading reference system with sensor units, a control system, a power system, and a thruster system [1]. The position and heading reference system, along with the wind sensors and gyrocompasses, can provide the control system with information on the position of the vessel and the magnitude and direction of the environmental forces that affect the position. All of the information obtained is then used to adjust the vessel itself.Figure 1
Dynamic positioning system configured for 8 thrusters of a semisubmersible drilling rig.The DP system must be very responsive to changes in the conditions due to weather, water depth, and so on in order to properly control the operating parameters of the vessel. Despite the importance of the DP system to the vessel, its use in a real working environment still faces many challenges, such as a power failure and a thruster malfunction, due to its complexity. The fuel expenditure of the DP system is a considerable burden since it requires a large amount of energy. This means that the DP system is overactuated, which renders the thruster to have infinite allocation solutions. Therefore, the thrust allocation can be formulated as an optimization process under numerous constraints in order to minimize power consumption.Related Works. In the literature, many methods have been proposed to achieve thrust allocation. Johansen and Fossen [2] provided a recent review of thrust allocation methods, but most of these have focused on proposing optimization algorithms intended for thrust allocation for industries, such as automobile, aerospace [3], and robotics, and not for ocean engineering. The relationship between the power consumption and the resulting thrust can usually be approximated by using a quadratic function, so some researchers have formulated optimization problems for thrust allocation as quadratic programming problems and have solved these with some variant of quadratic programming (QP) algorithms [4].Some of the methods intended for ocean engineering have considered cases with ships or with offshore vessels. Wit [5] proposed a method to achieve the optimal thrust allocation for a ship. The thrust and azimuth direction of the thrusters were selected as design variables for the optimization problem, and the QP algorithm, which is a local optimization algorithm, was used to solve the problem. Similarly, Liang and Cheng [6] and Johansen et al. [7] proposed methods for the ship that include a similar formulation of an optimization problem and a solution through the use of a sequential quadratic programming (SQP) algorithm. Parikshit [8] also proposed a method to achieve the optimal thrust allocation for a drilling rig. The thrust and azimuth direction of thrusters were selected as design variables of an optimization problem, and various global optimization algorithms, such as genetic algorithm (GA) and ITHS (Intelligent Tuned Harmony Search), were used to solve this problem. Zhao et al. [9] also proposed a method that includes a similar formulation of an optimization problem and solution for a drilling rig by using GA.In this manner, most studies in the field of ocean engineering have determined the thrust and the azimuth direction of the thrusters to be design variables for an optimization problem. In fact, the thrust can be expressed as a function of the speed of rotation, the diameter, and the thrust coefficient of a thruster. In order to reflect the specifications of the given thruster and to more easily control the thruster, the speed of rotation can be used as a design variable instead of the thrust. Meanwhile, a local optimization algorithm can be used as the optimization algorithm, such as the SQP algorithm, to derive an accurate optimum. However, this method is very sensitive to a starting point for optimization and sometimes finds a local optimum. On the other hand, a global optimization algorithm, such as the GA, does not need a starting point for the optimization but can only find a rough optimum [10].This study proposes an optimal thrust allocation method for an offshore vessel with the aim of achieving the required generalized forces and moment for dynamic positioning while at the same time minimizing total power. Finding such a method is a particularly important issue to operate an offshore vessel while using minimal energy. The method involves mathematically formulating an optimization problem for thrust allocation and selecting and using a suitable optimization algorithm to solve the problem, not developing a new algorithm. When the optimization problem is formulated, the speed of the rotation and the azimuth direction of thrusters are selected as design variables, the total power is selected as an objective function, and the required forces and moment are selected as constraints. A hybrid optimization algorithm that incorporates global and local optimization algorithms was used for the optimization.Table1 shows a summary of related studies and a comparison of them with this study.Table 1
Summary of related studies and comparison with this study.
Studies
Application
Optimization algorithm
Design variables for formulation
Consideration of thruster interaction
Wit [5]
Ship
Local algorithm (QP)
Thrust and azimuth direction
◯
Liang and Cheng [6], Johansen et al. [7]
Ship
Local algorithm (SQP)
Thrust and azimuth direction
×
Parikshit [8]
Drilling rig
Global algorithm (GA, ITHS, etc.)
Thrust and azimuth direction
◯
Zhao et al. [9]
Drilling rig
Global algorithm (GA)
Thrust and azimuth direction
◯
This study
Drilling rig
Global algorithm (hybrid optimization algorithm)
Rotation speed and azimuth direction
◯


## 2. Thrust Allocation for Dynamic Positioning

A thruster system is important for a drilling rig to maintain its position and heading since it can simultaneously provide transverse and longitudinal thrust. An azimuth thruster consists of an electric podded drive that is fitted to the hull and can rotate 360 degrees along the horizontal angle to provide thrust in all directions. The azimuth thruster has been selected in this study because it can provide thrust in any direction in order to act as a propulsor for dynamic positioning and target optimization of a drilling rig. However, the thruster system is usually overactuated in practice, so a thrust allocation problem can be formulated into a constrained optimization problem for the system.Figure2 depicts a schematic diagram for the thruster arrangement of the drilling rig and the sign conventions for the forces (Fx and Fy) and moment (Mz). In this figure, ∑Fx, ∑Fy, and ∑Mz represent the total forces required in the longitudinal (x) and transverse (y) directions and the total required moment about the vertical (z) direction from the control system, as shown in Figure 1, respectively. Note that the positions of the thrusters are given in a 2D coordinate system (only the horizontal forces and moment are considered). The coordinate system has its origin at the center of gravity (CoG) of the drilling rig.Figure 2
Schematic of thruster arrangement of a drilling rig with eight thrusters (Plan view).The environmental forces and the moment tend to move the drilling rig away from its original position during operation at sea. In order to move the drilling rig back to its original or reference position, the control system in the DP system of the drilling rig should first calculate the total forces (thrusts) and moment required. At this time, thrust allocation should be conducted to determine the thrust and azimuth direction of each thruster so that the required forces (and moment) are generated in each of the longitudinal (∑Fx) and transverse (∑Fy) and vertical (∑Mz) directions. The relationship between the thrust and the azimuth direction of each thruster and the total forces and moment that are required can be stated as in (1) to (3). These equations are the governing equations for the thrust allocation and can be used as equality constraints in an optimization problem:(1)∑Fx=F1cos⁡α1+F2cos⁡α2+⋯+Fncos⁡αn,(2)∑Fy=F1sin⁡α1+F2sin⁡α2+⋯+Fnsin⁡αn,(3)∑Mz=F1cos⁡α1·ly1+F2cos⁡α2·ly2+⋯+Fncos⁡αn·lyn+F1sin⁡α1·lx1+F2sin⁡α2·lx2+⋯+Fnsin⁡αn·lxn.Supposing that the drilling rig is equipped withi azimuth thrusters, Fi is the thrust from ith thruster and αi is the azimuth direction (angle) of ith thruster. lxi and lyi are the distances of ith thruster from the CoG of the drilling rig. Equation (1) indicates that the longitudinal component of the force generated from the thrust and azimuth direction of each thruster should be equal to the total force required in the longitudinal direction. Similarly, (2) indicates that the transverse component of the force from the thrust and azimuth direction of each thruster should be equal to the total force required in the transverse direction. Equation (3) means that the vertical component of the moment generated from the thrust and azimuth direction of each thruster should be equal to the total required moment about the vertical direction.The goal of achieving optimal thrust allocation is to guarantee efficient power use when dynamically positioning the drilling rig in the specific conditions of a given ocean environment. Meanwhile, the economic impact of minimizing power consumption is also taken into account with the ultimate aim of ensuring the stability and capability of the drilling rig during operation.


## 3. Thrust Allocation Method

### 3.1. Optimization Problem for Thrust Allocation

This study mathematically formulated an optimization problem for thrust allocation, and each component of the optimization problem is described below in detail.


#### 3.1.1. Design Variables

Thrust allocation involves determining the thrust (Fi) and azimuth direction (αi) of each thruster to generate the forces in longitudinal and transverse directions and the moment about the vertical direction that are required for the DP system of the drilling rig. Thus, the thrust and azimuth direction of each thruster can be used as design variables in the thrust allocation optimization problem. The objective function that is commonly used for this problem is to minimize the total power of the thrusters due to energy savings. In order to use this power as the objective function, it should be mathematically formulated as a function of the design variables. The total power will be proportional to the sum of the thrusts for each thruster. Some researchers [5–9] used the thrust and the azimuth direction of each thruster as design variables and have formulated the objective function (minimization of total power of thrusters) with the corresponding design variables. However, the rotational speed (ni) and the azimuth direction (αi) of each thruster can be used as design variables to reflect the characteristics of the given thruster, to more easily control the thruster, and to more delicately represent the objective function. Thus, the design variables for the thrust allocation optimization problem in this study are as follows:ni: speed of rotation of each thruster.αi: azimuth direction of each thruster.


#### 3.1.2. Objective Function

In this study, the speed of the rotation (ni) and the azimuth direction (αi) of each of the thrusters were used as design variables, and thus the objective function (minimization of the total power for the thrusters) was mathematically formulated according to the following procedure where in fact the power consumed by the electric propulsion system depends on the speed of rotation and the azimuth direction of the thrusters.The thrustT and torque Q of a thruster can be expressed as functions of the speed of rotation (n):(4)T=ρD4KTn2,(5)Q=ρD5KQn2,where ρ is the density of the sea water, D is the propeller diameter of the thruster, KT and KQ are the thrust and torque coefficients, respectively, and n is the speed of the rotation of the thruster. The expressions for KT and KQ can be obtained from the so-called open water test of the thruster. These nondimensional coefficients for thrust and torque were derived through a regression analysis that can be described according to the following parameters [11]:(6)KT=∑s,t,u,vCs,t,u,vTJsPDtAEAOuZv,KQ=∑s,t,u,vCs,t,u,vQJsPDtAEAOuZv.C
s
,
t
,
u
,
v
T and Cs,t,u,vQ coefficients and s, t, u, and v terms can be found in a report of the Wageningen B-Series propellers [11] and were used in this study. Here, P/D is the pitch diameter ratio, AE/AO is the blade area ratio, and Z is the number of blades of the thruster. The coefficient for the advance J is a nondimensional description of the thruster performance that can be expressed according to the speed of rotation n and the speed of the advance of the thruster defined as VA, which can be represented by the speed near each thruster (V) and the azimuth direction of each thruster (α). Since each of the thrusters has a different speed of rotation in practice, KT and KQ of each thruster can be different:(7)J=VAnD=Vcos⁡αnD.Now, it is time to formulate the total power consumed by the thrusters of the drilling rig. The power consumption (P) can be formulated as(8)P=2πniQi.By combining (5) and (8), the power consumption of the thrusters can be rewritten as(9)P=2πρD5∑i=1NKQini3,where N is the number of thrusters.Therefore, the objective function can be written relating the speed of rotation and the azimuth direction of each thruster. In this study, the total power (power consumption) of the thrusters is minimized by setting the objective function for the thrust allocation optimization problem as(10)MinimizeP=2πρD5∑i=1NKQini3.


#### 3.1.3. Constraints

To obtain a valid thrust allocation, the governing equations shown in (11) to (13) should be satisfied. The equations represent the thrust and azimuth direction of each thruster that generate the forces required in longitudinal (∑Fx) and transverse (∑Fy) directions as well as the moment required about vertical direction (∑Mz) to dynamically position the vessel. If these equations cannot be satisfied, the drilling rig starts to drift. Thus, these equations are used as equality constraints for the optimization problem.The governing equation for the longitudinal force can be stated as(11)∑Fx=ρD2KTin12cos⁡α1+n22cos⁡α2+⋯+ni2cos⁡αi.Similarly, the governing equation for the transverse force can be stated as(12)∑Fy=ρD2KTin12sin⁡α1+n22sin⁡α2+⋯+ni2sin⁡αi.And finally the governing equation for the yawing moment can be stated as(13)∑Mz=ρD2KTin12sin⁡α1lx1+n22sin⁡α2lx2+⋯+ni2sin⁡αilxi+n12cos⁡α1ly1+n22cos⁡α2ly2+⋯+ni2cos⁡αilyi,where ∑Fx, ∑Fy, and ∑Mz are the forces required in the longitudinal and transverse directions and the moment required about the vertical direction, respectively, and are all given from the control system in the dynamic positioning system.It is essential to minimize the thruster interaction of a semisubmersible drilling rig to ensure effective DP operation [5, 12]. In order to find a solution to provide the most effective thrust for operation, additional constraints should be considered by taking into account the physical locations of the thrusters. Thrusters that are close together have an influence on each other, and if one thruster is in the stream of another, the efficiency would drop significantly. To avoid having this happen, certain angles must be prohibited for the azimuth direction. Each thruster is assumed to have a forbidden angle of θf for the closest thruster in order to prevent the direct interaction of the thruster, as shown in Figure 3. With the port-forward pair, the center of the angle of the outboard thruster i is ϕi and for the inboard thruster (i+1) it is ϕi+1. The positive angle is measured clockwise while the negative angle is measured anticlockwise from the horizontal axis. Figure 3 shows the ATRs (Attainable Thrust Regions) for the two port-forward side thrusters i and (i+1). The ATR represents a set of all physically possible surge and sway forces for ith thruster.Figure 3
Interaction of two closely spaced azimuth thrusters.Considering the ATR, the constraints for the azimuth direction of thrustersi and (i+1) can be formulated as (14)θf-αi-ϕi≤0,θf-αi+1-ϕi+1≤0,where θf represents the forbidden angle of the thrusters. αi and αi+1 are the azimuth directions of the thrusters i and (i+1), respectively. ϕi and ϕi+1 are the center of angle for the outboard thruster i and for the inboard thruster (i+1), respectively.


#### 3.1.4. Summary

The optimization problem for thrust allocation can be summarized as follows:(15)Minimizefx=P=2πρD5∑i=1nKQini3,Subject tog1x=∑Fx-ρD2KTin12cos⁡α1+n22cos⁡α2+⋯+ni2cos⁡αi=0,g2x=∑Fy-ρD2KTin12sin⁡α1+n22sin⁡α2+⋯+ni2sin⁡αi=0,g3x=∑Mz-ρD2KTin12sin⁡α1lx1+n22sin⁡α2lx2+⋯+ni2sin⁡αilxi+n12cos⁡α1ly1+n22cos⁡α2ly2+⋯+ni2cos⁡αilyi=0,g4+ix=θf-αi-ϕi≤0,where x=ni,αi,  i=1,…,N(N: number of thrusters).Thus, this problem has one objective function, three equality constraints, andN inequality constraints.


### 3.2. Hybrid Optimization Algorithm for Thrust Allocation

#### 3.2.1. Overview

The optimization algorithms are generally divided into two categories: global and local. Several classes of global optimization algorithms are now available, including the genetic algorithm (GA) [13, 14] and the simulated annealing method. These algorithms are suitable for large-scale problems with many local optima. However, these algorithms require further iteration to obtain an acceptable optimum and not an accurate optimum, in contrast to the local optimization algorithms. Several classes of the local optimization algorithms also exist, including sequential linear programming (SLP) [15] and sequential quadratic programming (SQP) [15] and the method of feasible directions (MFD) [16]. Each of these algorithms can effectively determine an accurate optimum. However, in some cases, these algorithms find the local optimum that is closest to the given starting point.Various attempts have been made to combine global and local optimization algorithms in order to overcome the challenges of using them separately [10, 17–21]. Most of these studies used a combination of a global (e.g., GA) and a local optimization algorithm, which can be referred to as a hybrid optimization algorithm. In this study, the hybrid optimization algorithm was selected and used by incorporating the GA and the SQP algorithm to solve an optimization problem for thrust allocation, as described in Section 3.1. In this case, the SQP algorithm was used to improve the acceptable, global optimum obtained from the GA. That is, at the end of the GA, the SQP algorithm was executed with a starting point that corresponds to the optimum obtained from the GA. Figure 4 shows the general procedure of the hybrid optimization algorithm used in this study.Figure 4
General procedure of the hybrid optimization algorithm used in this study.An experiment on the mathematical optimization problem was performed in order to verify the efficiency, accuracy, and applicability of the hybrid optimization algorithm that was used to optimize the thrust allocation in this study. The selected problem, which is Rastrigin’s problem, one of the benchmark problems, is being widely used to check the efficiency of optimization algorithms [22]. More details about comparative test of the hybrid optimization algorithm can be found in [10].


#### 3.2.2. Optimization Procedure

The following optimization procedure was established to solve the optimization problem that was presented formulated above. First, the initial values were assumed for the design variables. At this time, the values can be randomly generated or can be manually set or extracted from an existing design. Now, these values are transferred to the hybrid optimization algorithm, and the values of an objective function and constraints are then calculated. We then check whether the current values of the design variables are at an optimum or not. If yes, the optimization process finishes and the result will be shown, and if not, the above steps will be repeated until the optimum is found. Figure5 shows this optimization procedure to find an optimal thrust allocation.Figure 5
Optimization procedure for finding optimal thrust allocation.


## 3.1. Optimization Problem for Thrust Allocation

This study mathematically formulated an optimization problem for thrust allocation, and each component of the optimization problem is described below in detail.


### 3.1.1. Design Variables

Thrust allocation involves determining the thrust (Fi) and azimuth direction (αi) of each thruster to generate the forces in longitudinal and transverse directions and the moment about the vertical direction that are required for the DP system of the drilling rig. Thus, the thrust and azimuth direction of each thruster can be used as design variables in the thrust allocation optimization problem. The objective function that is commonly used for this problem is to minimize the total power of the thrusters due to energy savings. In order to use this power as the objective function, it should be mathematically formulated as a function of the design variables. The total power will be proportional to the sum of the thrusts for each thruster. Some researchers [5–9] used the thrust and the azimuth direction of each thruster as design variables and have formulated the objective function (minimization of total power of thrusters) with the corresponding design variables. However, the rotational speed (ni) and the azimuth direction (αi) of each thruster can be used as design variables to reflect the characteristics of the given thruster, to more easily control the thruster, and to more delicately represent the objective function. Thus, the design variables for the thrust allocation optimization problem in this study are as follows:ni: speed of rotation of each thruster.αi: azimuth direction of each thruster.


### 3.1.2. Objective Function

In this study, the speed of the rotation (ni) and the azimuth direction (αi) of each of the thrusters were used as design variables, and thus the objective function (minimization of the total power for the thrusters) was mathematically formulated according to the following procedure where in fact the power consumed by the electric propulsion system depends on the speed of rotation and the azimuth direction of the thrusters.The thrustT and torque Q of a thruster can be expressed as functions of the speed of rotation (n):(4)T=ρD4KTn2,(5)Q=ρD5KQn2,where ρ is the density of the sea water, D is the propeller diameter of the thruster, KT and KQ are the thrust and torque coefficients, respectively, and n is the speed of the rotation of the thruster. The expressions for KT and KQ can be obtained from the so-called open water test of the thruster. These nondimensional coefficients for thrust and torque were derived through a regression analysis that can be described according to the following parameters [11]:(6)KT=∑s,t,u,vCs,t,u,vTJsPDtAEAOuZv,KQ=∑s,t,u,vCs,t,u,vQJsPDtAEAOuZv.C
s
,
t
,
u
,
v
T and Cs,t,u,vQ coefficients and s, t, u, and v terms can be found in a report of the Wageningen B-Series propellers [11] and were used in this study. Here, P/D is the pitch diameter ratio, AE/AO is the blade area ratio, and Z is the number of blades of the thruster. The coefficient for the advance J is a nondimensional description of the thruster performance that can be expressed according to the speed of rotation n and the speed of the advance of the thruster defined as VA, which can be represented by the speed near each thruster (V) and the azimuth direction of each thruster (α). Since each of the thrusters has a different speed of rotation in practice, KT and KQ of each thruster can be different:(7)J=VAnD=Vcos⁡αnD.Now, it is time to formulate the total power consumed by the thrusters of the drilling rig. The power consumption (P) can be formulated as(8)P=2πniQi.By combining (5) and (8), the power consumption of the thrusters can be rewritten as(9)P=2πρD5∑i=1NKQini3,where N is the number of thrusters.Therefore, the objective function can be written relating the speed of rotation and the azimuth direction of each thruster. In this study, the total power (power consumption) of the thrusters is minimized by setting the objective function for the thrust allocation optimization problem as(10)MinimizeP=2πρD5∑i=1NKQini3.


### 3.1.3. Constraints

To obtain a valid thrust allocation, the governing equations shown in (11) to (13) should be satisfied. The equations represent the thrust and azimuth direction of each thruster that generate the forces required in longitudinal (∑Fx) and transverse (∑Fy) directions as well as the moment required about vertical direction (∑Mz) to dynamically position the vessel. If these equations cannot be satisfied, the drilling rig starts to drift. Thus, these equations are used as equality constraints for the optimization problem.The governing equation for the longitudinal force can be stated as(11)∑Fx=ρD2KTin12cos⁡α1+n22cos⁡α2+⋯+ni2cos⁡αi.Similarly, the governing equation for the transverse force can be stated as(12)∑Fy=ρD2KTin12sin⁡α1+n22sin⁡α2+⋯+ni2sin⁡αi.And finally the governing equation for the yawing moment can be stated as(13)∑Mz=ρD2KTin12sin⁡α1lx1+n22sin⁡α2lx2+⋯+ni2sin⁡αilxi+n12cos⁡α1ly1+n22cos⁡α2ly2+⋯+ni2cos⁡αilyi,where ∑Fx, ∑Fy, and ∑Mz are the forces required in the longitudinal and transverse directions and the moment required about the vertical direction, respectively, and are all given from the control system in the dynamic positioning system.It is essential to minimize the thruster interaction of a semisubmersible drilling rig to ensure effective DP operation [5, 12]. In order to find a solution to provide the most effective thrust for operation, additional constraints should be considered by taking into account the physical locations of the thrusters. Thrusters that are close together have an influence on each other, and if one thruster is in the stream of another, the efficiency would drop significantly. To avoid having this happen, certain angles must be prohibited for the azimuth direction. Each thruster is assumed to have a forbidden angle of θf for the closest thruster in order to prevent the direct interaction of the thruster, as shown in Figure 3. With the port-forward pair, the center of the angle of the outboard thruster i is ϕi and for the inboard thruster (i+1) it is ϕi+1. The positive angle is measured clockwise while the negative angle is measured anticlockwise from the horizontal axis. Figure 3 shows the ATRs (Attainable Thrust Regions) for the two port-forward side thrusters i and (i+1). The ATR represents a set of all physically possible surge and sway forces for ith thruster.Figure 3
Interaction of two closely spaced azimuth thrusters.Considering the ATR, the constraints for the azimuth direction of thrustersi and (i+1) can be formulated as (14)θf-αi-ϕi≤0,θf-αi+1-ϕi+1≤0,where θf represents the forbidden angle of the thrusters. αi and αi+1 are the azimuth directions of the thrusters i and (i+1), respectively. ϕi and ϕi+1 are the center of angle for the outboard thruster i and for the inboard thruster (i+1), respectively.


### 3.1.4. Summary

The optimization problem for thrust allocation can be summarized as follows:(15)Minimizefx=P=2πρD5∑i=1nKQini3,Subject tog1x=∑Fx-ρD2KTin12cos⁡α1+n22cos⁡α2+⋯+ni2cos⁡αi=0,g2x=∑Fy-ρD2KTin12sin⁡α1+n22sin⁡α2+⋯+ni2sin⁡αi=0,g3x=∑Mz-ρD2KTin12sin⁡α1lx1+n22sin⁡α2lx2+⋯+ni2sin⁡αilxi+n12cos⁡α1ly1+n22cos⁡α2ly2+⋯+ni2cos⁡αilyi=0,g4+ix=θf-αi-ϕi≤0,where x=ni,αi,  i=1,…,N(N: number of thrusters).Thus, this problem has one objective function, three equality constraints, andN inequality constraints.


## 3.1.1. Design Variables

Thrust allocation involves determining the thrust (Fi) and azimuth direction (αi) of each thruster to generate the forces in longitudinal and transverse directions and the moment about the vertical direction that are required for the DP system of the drilling rig. Thus, the thrust and azimuth direction of each thruster can be used as design variables in the thrust allocation optimization problem. The objective function that is commonly used for this problem is to minimize the total power of the thrusters due to energy savings. In order to use this power as the objective function, it should be mathematically formulated as a function of the design variables. The total power will be proportional to the sum of the thrusts for each thruster. Some researchers [5–9] used the thrust and the azimuth direction of each thruster as design variables and have formulated the objective function (minimization of total power of thrusters) with the corresponding design variables. However, the rotational speed (ni) and the azimuth direction (αi) of each thruster can be used as design variables to reflect the characteristics of the given thruster, to more easily control the thruster, and to more delicately represent the objective function. Thus, the design variables for the thrust allocation optimization problem in this study are as follows:ni: speed of rotation of each thruster.αi: azimuth direction of each thruster.


## 3.1.2. Objective Function

In this study, the speed of the rotation (ni) and the azimuth direction (αi) of each of the thrusters were used as design variables, and thus the objective function (minimization of the total power for the thrusters) was mathematically formulated according to the following procedure where in fact the power consumed by the electric propulsion system depends on the speed of rotation and the azimuth direction of the thrusters.The thrustT and torque Q of a thruster can be expressed as functions of the speed of rotation (n):(4)T=ρD4KTn2,(5)Q=ρD5KQn2,where ρ is the density of the sea water, D is the propeller diameter of the thruster, KT and KQ are the thrust and torque coefficients, respectively, and n is the speed of the rotation of the thruster. The expressions for KT and KQ can be obtained from the so-called open water test of the thruster. These nondimensional coefficients for thrust and torque were derived through a regression analysis that can be described according to the following parameters [11]:(6)KT=∑s,t,u,vCs,t,u,vTJsPDtAEAOuZv,KQ=∑s,t,u,vCs,t,u,vQJsPDtAEAOuZv.C
s
,
t
,
u
,
v
T and Cs,t,u,vQ coefficients and s, t, u, and v terms can be found in a report of the Wageningen B-Series propellers [11] and were used in this study. Here, P/D is the pitch diameter ratio, AE/AO is the blade area ratio, and Z is the number of blades of the thruster. The coefficient for the advance J is a nondimensional description of the thruster performance that can be expressed according to the speed of rotation n and the speed of the advance of the thruster defined as VA, which can be represented by the speed near each thruster (V) and the azimuth direction of each thruster (α). Since each of the thrusters has a different speed of rotation in practice, KT and KQ of each thruster can be different:(7)J=VAnD=Vcos⁡αnD.Now, it is time to formulate the total power consumed by the thrusters of the drilling rig. The power consumption (P) can be formulated as(8)P=2πniQi.By combining (5) and (8), the power consumption of the thrusters can be rewritten as(9)P=2πρD5∑i=1NKQini3,where N is the number of thrusters.Therefore, the objective function can be written relating the speed of rotation and the azimuth direction of each thruster. In this study, the total power (power consumption) of the thrusters is minimized by setting the objective function for the thrust allocation optimization problem as(10)MinimizeP=2πρD5∑i=1NKQini3.


## 3.1.3. Constraints

To obtain a valid thrust allocation, the governing equations shown in (11) to (13) should be satisfied. The equations represent the thrust and azimuth direction of each thruster that generate the forces required in longitudinal (∑Fx) and transverse (∑Fy) directions as well as the moment required about vertical direction (∑Mz) to dynamically position the vessel. If these equations cannot be satisfied, the drilling rig starts to drift. Thus, these equations are used as equality constraints for the optimization problem.The governing equation for the longitudinal force can be stated as(11)∑Fx=ρD2KTin12cos⁡α1+n22cos⁡α2+⋯+ni2cos⁡αi.Similarly, the governing equation for the transverse force can be stated as(12)∑Fy=ρD2KTin12sin⁡α1+n22sin⁡α2+⋯+ni2sin⁡αi.And finally the governing equation for the yawing moment can be stated as(13)∑Mz=ρD2KTin12sin⁡α1lx1+n22sin⁡α2lx2+⋯+ni2sin⁡αilxi+n12cos⁡α1ly1+n22cos⁡α2ly2+⋯+ni2cos⁡αilyi,where ∑Fx, ∑Fy, and ∑Mz are the forces required in the longitudinal and transverse directions and the moment required about the vertical direction, respectively, and are all given from the control system in the dynamic positioning system.It is essential to minimize the thruster interaction of a semisubmersible drilling rig to ensure effective DP operation [5, 12]. In order to find a solution to provide the most effective thrust for operation, additional constraints should be considered by taking into account the physical locations of the thrusters. Thrusters that are close together have an influence on each other, and if one thruster is in the stream of another, the efficiency would drop significantly. To avoid having this happen, certain angles must be prohibited for the azimuth direction. Each thruster is assumed to have a forbidden angle of θf for the closest thruster in order to prevent the direct interaction of the thruster, as shown in Figure 3. With the port-forward pair, the center of the angle of the outboard thruster i is ϕi and for the inboard thruster (i+1) it is ϕi+1. The positive angle is measured clockwise while the negative angle is measured anticlockwise from the horizontal axis. Figure 3 shows the ATRs (Attainable Thrust Regions) for the two port-forward side thrusters i and (i+1). The ATR represents a set of all physically possible surge and sway forces for ith thruster.Figure 3
Interaction of two closely spaced azimuth thrusters.Considering the ATR, the constraints for the azimuth direction of thrustersi and (i+1) can be formulated as (14)θf-αi-ϕi≤0,θf-αi+1-ϕi+1≤0,where θf represents the forbidden angle of the thrusters. αi and αi+1 are the azimuth directions of the thrusters i and (i+1), respectively. ϕi and ϕi+1 are the center of angle for the outboard thruster i and for the inboard thruster (i+1), respectively.


## 3.1.4. Summary

The optimization problem for thrust allocation can be summarized as follows:(15)Minimizefx=P=2πρD5∑i=1nKQini3,Subject tog1x=∑Fx-ρD2KTin12cos⁡α1+n22cos⁡α2+⋯+ni2cos⁡αi=0,g2x=∑Fy-ρD2KTin12sin⁡α1+n22sin⁡α2+⋯+ni2sin⁡αi=0,g3x=∑Mz-ρD2KTin12sin⁡α1lx1+n22sin⁡α2lx2+⋯+ni2sin⁡αilxi+n12cos⁡α1ly1+n22cos⁡α2ly2+⋯+ni2cos⁡αilyi=0,g4+ix=θf-αi-ϕi≤0,where x=ni,αi,  i=1,…,N(N: number of thrusters).Thus, this problem has one objective function, three equality constraints, andN inequality constraints.


## 3.2. Hybrid Optimization Algorithm for Thrust Allocation

### 3.2.1. Overview

The optimization algorithms are generally divided into two categories: global and local. Several classes of global optimization algorithms are now available, including the genetic algorithm (GA) [13, 14] and the simulated annealing method. These algorithms are suitable for large-scale problems with many local optima. However, these algorithms require further iteration to obtain an acceptable optimum and not an accurate optimum, in contrast to the local optimization algorithms. Several classes of the local optimization algorithms also exist, including sequential linear programming (SLP) [15] and sequential quadratic programming (SQP) [15] and the method of feasible directions (MFD) [16]. Each of these algorithms can effectively determine an accurate optimum. However, in some cases, these algorithms find the local optimum that is closest to the given starting point.Various attempts have been made to combine global and local optimization algorithms in order to overcome the challenges of using them separately [10, 17–21]. Most of these studies used a combination of a global (e.g., GA) and a local optimization algorithm, which can be referred to as a hybrid optimization algorithm. In this study, the hybrid optimization algorithm was selected and used by incorporating the GA and the SQP algorithm to solve an optimization problem for thrust allocation, as described in Section 3.1. In this case, the SQP algorithm was used to improve the acceptable, global optimum obtained from the GA. That is, at the end of the GA, the SQP algorithm was executed with a starting point that corresponds to the optimum obtained from the GA. Figure 4 shows the general procedure of the hybrid optimization algorithm used in this study.Figure 4
General procedure of the hybrid optimization algorithm used in this study.An experiment on the mathematical optimization problem was performed in order to verify the efficiency, accuracy, and applicability of the hybrid optimization algorithm that was used to optimize the thrust allocation in this study. The selected problem, which is Rastrigin’s problem, one of the benchmark problems, is being widely used to check the efficiency of optimization algorithms [22]. More details about comparative test of the hybrid optimization algorithm can be found in [10].


### 3.2.2. Optimization Procedure

The following optimization procedure was established to solve the optimization problem that was presented formulated above. First, the initial values were assumed for the design variables. At this time, the values can be randomly generated or can be manually set or extracted from an existing design. Now, these values are transferred to the hybrid optimization algorithm, and the values of an objective function and constraints are then calculated. We then check whether the current values of the design variables are at an optimum or not. If yes, the optimization process finishes and the result will be shown, and if not, the above steps will be repeated until the optimum is found. Figure5 shows this optimization procedure to find an optimal thrust allocation.Figure 5
Optimization procedure for finding optimal thrust allocation.


## 3.2.1. Overview

The optimization algorithms are generally divided into two categories: global and local. Several classes of global optimization algorithms are now available, including the genetic algorithm (GA) [13, 14] and the simulated annealing method. These algorithms are suitable for large-scale problems with many local optima. However, these algorithms require further iteration to obtain an acceptable optimum and not an accurate optimum, in contrast to the local optimization algorithms. Several classes of the local optimization algorithms also exist, including sequential linear programming (SLP) [15] and sequential quadratic programming (SQP) [15] and the method of feasible directions (MFD) [16]. Each of these algorithms can effectively determine an accurate optimum. However, in some cases, these algorithms find the local optimum that is closest to the given starting point.Various attempts have been made to combine global and local optimization algorithms in order to overcome the challenges of using them separately [10, 17–21]. Most of these studies used a combination of a global (e.g., GA) and a local optimization algorithm, which can be referred to as a hybrid optimization algorithm. In this study, the hybrid optimization algorithm was selected and used by incorporating the GA and the SQP algorithm to solve an optimization problem for thrust allocation, as described in Section 3.1. In this case, the SQP algorithm was used to improve the acceptable, global optimum obtained from the GA. That is, at the end of the GA, the SQP algorithm was executed with a starting point that corresponds to the optimum obtained from the GA. Figure 4 shows the general procedure of the hybrid optimization algorithm used in this study.Figure 4
General procedure of the hybrid optimization algorithm used in this study.An experiment on the mathematical optimization problem was performed in order to verify the efficiency, accuracy, and applicability of the hybrid optimization algorithm that was used to optimize the thrust allocation in this study. The selected problem, which is Rastrigin’s problem, one of the benchmark problems, is being widely used to check the efficiency of optimization algorithms [22]. More details about comparative test of the hybrid optimization algorithm can be found in [10].


## 3.2.2. Optimization Procedure

The following optimization procedure was established to solve the optimization problem that was presented formulated above. First, the initial values were assumed for the design variables. At this time, the values can be randomly generated or can be manually set or extracted from an existing design. Now, these values are transferred to the hybrid optimization algorithm, and the values of an objective function and constraints are then calculated. We then check whether the current values of the design variables are at an optimum or not. If yes, the optimization process finishes and the result will be shown, and if not, the above steps will be repeated until the optimum is found. Figure5 shows this optimization procedure to find an optimal thrust allocation.Figure 5
Optimization procedure for finding optimal thrust allocation.


## 4. Application of the Thrust Allocation Method

### 4.1. Optimization Target

The optimization target for this study is a semisubmersible drilling rig equipped with eight azimuth thrusters. Figure5 shows the schematics of this drilling rig with a thruster system and its position. We assumed that each thruster has 4 propeller blades (Z), a propeller diameter (D) of 3.6 m, a pitch diameter ratio (P/D) of 1.1, a blade area ratio (AE/AO) of 0.7, and a maximum thrust of 600 KN. In addition, we assumed that the speed (V) near each thruster was of 1 m/s. The longitudinal and transverse positions of each thruster can be seen in Figure 6.Figure 6
The longitudinal and transverse positions of eight azimuth thrusters (Plan view).


### 4.2. Mathematical Formulation

This drilling rig has eight thrusters. Thus, the design variables in this problem are as follows:(16)x=n1,α1,n2,α2,n3,α3,n4,α4,n5,α5,n6,α6,n7,α7,n8,α8,where ni and αi are the speed of rotation and the azimuth direction of ith thruster, respectively, and thus this problem has 16 design variables.Now, the objective function of this problem can be stated by using (10):(17)Minimizefx=P=2πρD5∑i=18KQini3=2πρD5KQ1n13+KQ2n23+KQ3n33+KQ4n43+KQ5n53+KQ6n63+KQ7n73+KQ8n83,where ρ is the density of sea water, D is the propeller diameter of the thruster, and KT and KQ are the thruster and torque coefficients, respectively. KT and KQ were calculated from (6) by referring to the report on the Wageningen B-Series propellers [11].The constraints for the required forces and moment can be stated using (11) to (13):(18)g1x=∑Fx-ρD2KT1n12cos⁡α1+KT2n22cos⁡α2+KT3n32cos⁡α3+KT4n42cos⁡α4+KT5n52cos⁡α5+KT6n62cos⁡α6+KT7n72cos⁡α7+KT8n82cos⁡α8=0,(19)g2x=∑Fy-ρD2KT1n12sin⁡α1+KT2n22sin⁡α2+KT3n32sin⁡α3+KT4n42sin⁡α4+KT5n52sin⁡α5+KT6n62sin⁡α6+KT7n72sin⁡α7+KT8n82sin⁡α8=0,(20)g3x=∑Mz-ρD2KT1n12sin⁡α1lx1+KT2n22sin⁡α2lx2+KT3n32sin⁡α3lx3+KT4n42sin⁡α4lx4+KT5n52sin⁡α5lx5+KT6n62sin⁡α6lx6+KT7n72sin⁡α7lx7+KT8n82sin⁡α8lx8+KT1n12cos⁡α1ly1+KT2n22cos⁡α2ly2+KT3n32cos⁡α3ly3+KT4n42cos⁡α4ly4+KT5n52cos⁡α5ly5+KT6n62cos⁡α6ly6+KT7n72cos⁡α7ly7+KT8n82cos⁡α8ly8=0.Equations (18) to (20) ensure that the forces generalized by the eight thrusters are equal to the forces demanded from the control system, and the moment should also be satisfied in the same way to keep the drilling rig in the desired position.In order to overcome the energy loss due to the thruster-thruster interaction, the ATR should be considered, excluding forbidden zones, by using (14), as shown in Figure 7. For this, each thruster should satisfy the following equations simultaneously:(21)g4x=25∘-α1-63.44°≤0,g5x=25∘-α2+116.56°≤0,g6x=25∘-α3-116.56°≤0,g7x=25∘-α4+63.44°≤0,g8x=25∘-α5-63.44°≤0,g9x=25∘-α6+116.56°≤0,g10x=25∘-α7-116.56°≤0,g11x=25∘-α8+63.44°≤0.Figure 7
Schematic diagram of thruster interaction layout (Plan view).In addition, there are some limitations on the maximum thrust and azimuth direction for each thruster. Thus, such limitations were also considered as additional constraints in this study. The additional constraints about the maximum thrust and the azimuth direction of each thruster can be stated in the following equations, respectively:(22)g11+i=Ti-600KN≤0,(23)g19+i=αi-180°≤0,where i=1,…,8.Thus, this problem has 3 equality constraints and 24 inequality constraints.


### 4.3. Optimization Result

The problem in (16) to (23) was solved by using the hybrid optimization algorithm described in Section 3.2, and then the problem was also solved by using the GA in order to compare their relative efficiencies. We assumed that the total required forces (∑Fx,∑Fy) and moment (∑Mz) are given from the control system at 14 time steps, as stated in [8] and as shown in Table 2. For example, the values for the first time step are 50 KN, −600 KN, and −64,000 KN·m. For a total of 14 inputs for the total required forces and moment, the optimization was performed by using the GA and the hybrid optimization algorithm, and the optimization results are shown in Table 2. As mentioned earlier, the hybrid optimization algorithm is a combination of the GA and the SQP algorithm. That is, at the end of the GA, the SQP algorithm was executed with a starting point that corresponds to the optimum obtained from the GA. With respect to the parameters for the GA, the number of individuals in the population was 200, the maximum number of generations was 3,000, the crossover probability was 0.8, and the mutation probability was 0.01 [23]. As shown in Table 2, the hybrid optimization algorithm outperforms the GA although the GA yielded similar results for certain time steps. Figure 8 shows a convergence history of the optimization results for the first time step. This figure shows that the objective function is minimized and that the constraints are satisfied at the same time.Table 2
Optimization result for 14 time steps.
Time steps
Input from the control system
Thrusters
GA only  [A]
Hybrid [B](GA + SQP)
Ratio (B/A)
Thruster 1
Thruster 2
Thruster 3
Thruster 4
Thruster 5
Thruster 6
Thruster 7
Thruster 8
Required force
Required moment
RPM
Angle
RPM
Angle
RPM
Angle
RPM
Angle
RPM
Angle
RPM
Angle
RPM
Angle
RPM
Angle
Total power
Total power
(KN)
(KN⋅m)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(kW)
(kW)
Δ
t (s)
∑
F
x
∑
F
y
∑
M
z
n
1
α
1
n
2
α
2
n
3
α
3
n
4
α
4
n
5
α
5
n
6
α
6
n
7
α
7
n
8
α
8
P
P
—
1
50
−600
−64,000
3.10
−93
0.62
−91
0.09
85
0.55
−103
2.2
37
0.18
178
0.22
17.3
1.66
167
9,859
8,250
0.84
2
−500
1,500
−25,000
2.51
115
1.02
171
1.29
152
0.79
93
0.03
−64
1.90
78
1.52
81
2.28
114
12,882
10,824
0.84
3
100
−800
17,000
0.78
−35
1.57
−78
1.02
−2
0.54
−13
2.61
−95
0.26
−15
0.17
−92
0.45
20
7,186
6,580
0.92
4
750
−550
63,000
2026
34
0.61
−177
0.76
170
0.07
32
0.62
−78
1.99
−78
0.03
40
2.64
−48
11,800
9,172
0.78
5
175
560
−32,500
0.08
13
0.92
147
0.05
41
0.73
−128
2.74
74
1.11
−35
0.15
−41
0.29
−146
7,121
7,107
1.00
6
−50
−1,080
3,000
2.17
−116
0.04
129
0.89
−65
0.32
−136
2.45
−85
1.42
−70
0.09
−141
0.04
−4
8,840
7,360
0.83
7
−300
1,560
12,000
0.58
129
2.18
68
2.29
97
0.26
54
1.98
108
1.28
−180
0.07
−65
2.29
112
13,749
9,522
0.70
8
1,050
250
−36,000
0.80
−84
0.52
144
0.64
29
2.40
9
0.1
165
2.78
11
1.32
51
0.13
−143
9,623
7,519
0.78
9
−600
−850
−22,000
2.09
−120
0.38
91
0.34
173
2.42
−102
0.21
17
1.97
−180
0.20
−106
0.8
−91
9,687
7,380
0.76
10
−1,200
250
−11,000
1.27
−165
1.31
169
2.28
−176
0.03
39
2.08
146
0.29
134
1.82
161
0.07
−37
9,206
8,140
0.88
11
−500
1,200
−10,000
0.03
125
0.94
148
2.47
97
0.13
−150
1.89
86
1.10
153
0.20
−9
2.32
138
11,280
10,277
0.91
12
600
550
−7,000
0.36
−10
1.93
26
0.04
−72
2.26
1
2
111
1.38
139
1.4
0
0.24
139
9,575
4,855
0.50
13
50
−930
1,000
1.10
−120
1.18
−92
1.72
−78
0.84
−89
2.07
−107
0.96
−88
1.4
2
0.16
62
6,541
3,949
0.60
14
−700
75
−26,000
1.23
−137
1.04
57
0.18
−12
0.18
−34
0.53
−8
0.06
9
2.27
168
1.86
180
6,243
4,802
0.77
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
Mean
0.79Figure 8
Convergence history of the optimization by the hybrid optimization algorithm for first time step.To evaluate the applicability of the proposed method, the optimization results of this study were compared against those obtained in Parikshit’s study [8]. The same inputs were used for the total forces and moment required from the control system for the comparative test. Parikshit used an ITHS (Intelligent Tuned Harmony Search) algorithm to solve this problem. The ITHS algorithm maintains a proper balance between the diversification and intensification throughout the search process by automatically selecting the proper pitch adjustment strategy based on its harmony memory. Table 3 shows a comparison of the optimization results between the existing study and this study. As shown in this table, the results of this study were slightly better than those of Parikshit’s study in terms of a mean value of the ratio.Table 3
Comparison of the optimization results between the existing study and this study.
Time steps
ITHS [A]
Hybrid [B]
Ratio (B/A)
Total power (kW)
Total power (kW)
Δ
t (s)
P
P
—
1
7,500
8,250
1.10
2
10,850
10,824
0.99
3
7,100
6,580
0.92
4
7,800
9,172
1.17
5
7,900
7,107
0.89
6
5,800
7,360
1.26
7
9,800
9,522
0.97
8
7,000
7,519
1.07
9
7,500
7,379
0.98
10
9,800
9,206
0.83
11
8,900
10,277
1.15
12
6,200
4,855
0.78
13
4,000
3,949
0.98
14
7,600
4,802
0.63
—
—
Mean
0.98


## 4.1. Optimization Target

The optimization target for this study is a semisubmersible drilling rig equipped with eight azimuth thrusters. Figure5 shows the schematics of this drilling rig with a thruster system and its position. We assumed that each thruster has 4 propeller blades (Z), a propeller diameter (D) of 3.6 m, a pitch diameter ratio (P/D) of 1.1, a blade area ratio (AE/AO) of 0.7, and a maximum thrust of 600 KN. In addition, we assumed that the speed (V) near each thruster was of 1 m/s. The longitudinal and transverse positions of each thruster can be seen in Figure 6.Figure 6
The longitudinal and transverse positions of eight azimuth thrusters (Plan view).


## 4.2. Mathematical Formulation

This drilling rig has eight thrusters. Thus, the design variables in this problem are as follows:(16)x=n1,α1,n2,α2,n3,α3,n4,α4,n5,α5,n6,α6,n7,α7,n8,α8,where ni and αi are the speed of rotation and the azimuth direction of ith thruster, respectively, and thus this problem has 16 design variables.Now, the objective function of this problem can be stated by using (10):(17)Minimizefx=P=2πρD5∑i=18KQini3=2πρD5KQ1n13+KQ2n23+KQ3n33+KQ4n43+KQ5n53+KQ6n63+KQ7n73+KQ8n83,where ρ is the density of sea water, D is the propeller diameter of the thruster, and KT and KQ are the thruster and torque coefficients, respectively. KT and KQ were calculated from (6) by referring to the report on the Wageningen B-Series propellers [11].The constraints for the required forces and moment can be stated using (11) to (13):(18)g1x=∑Fx-ρD2KT1n12cos⁡α1+KT2n22cos⁡α2+KT3n32cos⁡α3+KT4n42cos⁡α4+KT5n52cos⁡α5+KT6n62cos⁡α6+KT7n72cos⁡α7+KT8n82cos⁡α8=0,(19)g2x=∑Fy-ρD2KT1n12sin⁡α1+KT2n22sin⁡α2+KT3n32sin⁡α3+KT4n42sin⁡α4+KT5n52sin⁡α5+KT6n62sin⁡α6+KT7n72sin⁡α7+KT8n82sin⁡α8=0,(20)g3x=∑Mz-ρD2KT1n12sin⁡α1lx1+KT2n22sin⁡α2lx2+KT3n32sin⁡α3lx3+KT4n42sin⁡α4lx4+KT5n52sin⁡α5lx5+KT6n62sin⁡α6lx6+KT7n72sin⁡α7lx7+KT8n82sin⁡α8lx8+KT1n12cos⁡α1ly1+KT2n22cos⁡α2ly2+KT3n32cos⁡α3ly3+KT4n42cos⁡α4ly4+KT5n52cos⁡α5ly5+KT6n62cos⁡α6ly6+KT7n72cos⁡α7ly7+KT8n82cos⁡α8ly8=0.Equations (18) to (20) ensure that the forces generalized by the eight thrusters are equal to the forces demanded from the control system, and the moment should also be satisfied in the same way to keep the drilling rig in the desired position.In order to overcome the energy loss due to the thruster-thruster interaction, the ATR should be considered, excluding forbidden zones, by using (14), as shown in Figure 7. For this, each thruster should satisfy the following equations simultaneously:(21)g4x=25∘-α1-63.44°≤0,g5x=25∘-α2+116.56°≤0,g6x=25∘-α3-116.56°≤0,g7x=25∘-α4+63.44°≤0,g8x=25∘-α5-63.44°≤0,g9x=25∘-α6+116.56°≤0,g10x=25∘-α7-116.56°≤0,g11x=25∘-α8+63.44°≤0.Figure 7
Schematic diagram of thruster interaction layout (Plan view).In addition, there are some limitations on the maximum thrust and azimuth direction for each thruster. Thus, such limitations were also considered as additional constraints in this study. The additional constraints about the maximum thrust and the azimuth direction of each thruster can be stated in the following equations, respectively:(22)g11+i=Ti-600KN≤0,(23)g19+i=αi-180°≤0,where i=1,…,8.Thus, this problem has 3 equality constraints and 24 inequality constraints.


## 4.3. Optimization Result

The problem in (16) to (23) was solved by using the hybrid optimization algorithm described in Section 3.2, and then the problem was also solved by using the GA in order to compare their relative efficiencies. We assumed that the total required forces (∑Fx,∑Fy) and moment (∑Mz) are given from the control system at 14 time steps, as stated in [8] and as shown in Table 2. For example, the values for the first time step are 50 KN, −600 KN, and −64,000 KN·m. For a total of 14 inputs for the total required forces and moment, the optimization was performed by using the GA and the hybrid optimization algorithm, and the optimization results are shown in Table 2. As mentioned earlier, the hybrid optimization algorithm is a combination of the GA and the SQP algorithm. That is, at the end of the GA, the SQP algorithm was executed with a starting point that corresponds to the optimum obtained from the GA. With respect to the parameters for the GA, the number of individuals in the population was 200, the maximum number of generations was 3,000, the crossover probability was 0.8, and the mutation probability was 0.01 [23]. As shown in Table 2, the hybrid optimization algorithm outperforms the GA although the GA yielded similar results for certain time steps. Figure 8 shows a convergence history of the optimization results for the first time step. This figure shows that the objective function is minimized and that the constraints are satisfied at the same time.Table 2
Optimization result for 14 time steps.
Time steps
Input from the control system
Thrusters
GA only  [A]
Hybrid [B](GA + SQP)
Ratio (B/A)
Thruster 1
Thruster 2
Thruster 3
Thruster 4
Thruster 5
Thruster 6
Thruster 7
Thruster 8
Required force
Required moment
RPM
Angle
RPM
Angle
RPM
Angle
RPM
Angle
RPM
Angle
RPM
Angle
RPM
Angle
RPM
Angle
Total power
Total power
(KN)
(KN⋅m)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(rpm)
(deg)
(kW)
(kW)
Δ
t (s)
∑
F
x
∑
F
y
∑
M
z
n
1
α
1
n
2
α
2
n
3
α
3
n
4
α
4
n
5
α
5
n
6
α
6
n
7
α
7
n
8
α
8
P
P
—
1
50
−600
−64,000
3.10
−93
0.62
−91
0.09
85
0.55
−103
2.2
37
0.18
178
0.22
17.3
1.66
167
9,859
8,250
0.84
2
−500
1,500
−25,000
2.51
115
1.02
171
1.29
152
0.79
93
0.03
−64
1.90
78
1.52
81
2.28
114
12,882
10,824
0.84
3
100
−800
17,000
0.78
−35
1.57
−78
1.02
−2
0.54
−13
2.61
−95
0.26
−15
0.17
−92
0.45
20
7,186
6,580
0.92
4
750
−550
63,000
2026
34
0.61
−177
0.76
170
0.07
32
0.62
−78
1.99
−78
0.03
40
2.64
−48
11,800
9,172
0.78
5
175
560
−32,500
0.08
13
0.92
147
0.05
41
0.73
−128
2.74
74
1.11
−35
0.15
−41
0.29
−146
7,121
7,107
1.00
6
−50
−1,080
3,000
2.17
−116
0.04
129
0.89
−65
0.32
−136
2.45
−85
1.42
−70
0.09
−141
0.04
−4
8,840
7,360
0.83
7
−300
1,560
12,000
0.58
129
2.18
68
2.29
97
0.26
54
1.98
108
1.28
−180
0.07
−65
2.29
112
13,749
9,522
0.70
8
1,050
250
−36,000
0.80
−84
0.52
144
0.64
29
2.40
9
0.1
165
2.78
11
1.32
51
0.13
−143
9,623
7,519
0.78
9
−600
−850
−22,000
2.09
−120
0.38
91
0.34
173
2.42
−102
0.21
17
1.97
−180
0.20
−106
0.8
−91
9,687
7,380
0.76
10
−1,200
250
−11,000
1.27
−165
1.31
169
2.28
−176
0.03
39
2.08
146
0.29
134
1.82
161
0.07
−37
9,206
8,140
0.88
11
−500
1,200
−10,000
0.03
125
0.94
148
2.47
97
0.13
−150
1.89
86
1.10
153
0.20
−9
2.32
138
11,280
10,277
0.91
12
600
550
−7,000
0.36
−10
1.93
26
0.04
−72
2.26
1
2
111
1.38
139
1.4
0
0.24
139
9,575
4,855
0.50
13
50
−930
1,000
1.10
−120
1.18
−92
1.72
−78
0.84
−89
2.07
−107
0.96
−88
1.4
2
0.16
62
6,541
3,949
0.60
14
−700
75
−26,000
1.23
−137
1.04
57
0.18
−12
0.18
−34
0.53
−8
0.06
9
2.27
168
1.86
180
6,243
4,802
0.77
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
Mean
0.79Figure 8
Convergence history of the optimization by the hybrid optimization algorithm for first time step.To evaluate the applicability of the proposed method, the optimization results of this study were compared against those obtained in Parikshit’s study [8]. The same inputs were used for the total forces and moment required from the control system for the comparative test. Parikshit used an ITHS (Intelligent Tuned Harmony Search) algorithm to solve this problem. The ITHS algorithm maintains a proper balance between the diversification and intensification throughout the search process by automatically selecting the proper pitch adjustment strategy based on its harmony memory. Table 3 shows a comparison of the optimization results between the existing study and this study. As shown in this table, the results of this study were slightly better than those of Parikshit’s study in terms of a mean value of the ratio.Table 3
Comparison of the optimization results between the existing study and this study.
Time steps
ITHS [A]
Hybrid [B]
Ratio (B/A)
Total power (kW)
Total power (kW)
Δ
t (s)
P
P
—
1
7,500
8,250
1.10
2
10,850
10,824
0.99
3
7,100
6,580
0.92
4
7,800
9,172
1.17
5
7,900
7,107
0.89
6
5,800
7,360
1.26
7
9,800
9,522
0.97
8
7,000
7,519
1.07
9
7,500
7,379
0.98
10
9,800
9,206
0.83
11
8,900
10,277
1.15
12
6,200
4,855
0.78
13
4,000
3,949
0.98
14
7,600
4,802
0.63
—
—
Mean
0.98


## 5. Conclusions and Future Works

In this study, a thrust allocation method was proposed for a semisubmersible drilling rig in order to produce the generalized forces and moment that are required to dynamically position the rig while at the same time minimizing the total power consumed. First, a thrust allocation optimization problem was mathematically formulated with the corresponding design variables, objective function, and constraints. In terms of the design variables, the speed of rotation and the azimuth direction of each thruster were selected. As compared with some studies about optimal thruster allocation, the selection of the design variables can not only represent the objection function more precisely but also control the thruster more easily. The objective function aimed to minimize the total power of the thrusters, and, in terms of the constraints, the governing equations for the thrust and azimuth direction of each thruster were used to generate the forces required in longitudinal and transverse directions as well as the moment about the vertical direction to dynamically position the vessel. Some limitations were also used for each thruster, and additional constraints were introduced by considering the energy loss due to the thruster-thruster interaction and the maximum thrust and azimuth direction of each thruster. The hybrid optimization algorithm was then used to solve the formulated problem. Finally, the proposed method was applied to an example to find the optimal thrust allocation for the semisubmersible drilling rig with 8 thrusters. A comparative test was also performed as part of the current study. In this comparative test, the proposed method was observed to produce slightly better results (about 2% in terms of total power) for thrust allocation than existing study, and thus the proposed method could be used to better determine a strategy to allocate the thruster of the drilling rig.In the future, the total forces and moment required from the control system, which was the input of this study, will be estimated by considering the current position and the equations of motion of the drilling rig. That is, a more general method to dynamically position the rig will be further studied. In addition, we will improve the present hybrid optimization method and apply the improved method to the formulated problem in the future.



---

*Source: 183705-2015-10-01.xml*