Journal of Robotics

Volume 2012 (2012), Article ID 638394, 15 pages

http://dx.doi.org/10.1155/2012/638394

## Reconstruction of Riser Profiles by an Underwater Robot Using Inertial Navigation

^{1}Biomedical Engineering Program (PEB/COPPE), The Alberto Luiz Coimbra Institute for Graduate Studies and Research in Engineering, Universidade Federal do Rio de Janeiro, Avenida Horacio Macedo 2030, Bloco H-338, 21941-914 Rio de Janeiro, RJ, Brazil^{2}Department of Informatics and Automation, Universita Roma Tre, Via della Vasca Navale, 79, I 00146 Roma, Italy^{3}Subsin Engineering, Rua Beneditinos, 16, 12th floor, 20081-050 Rio de Janeiro, RJ, Brazil

Received 10 October 2011; Revised 11 January 2012; Accepted 15 January 2012

Academic Editor: Jorge Manuel Dias

Copyright © 2012 Luciano Luporini Menegaldo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper proposes a kinematic model and an inertial localization system architecture for a riser inspecting robot. The robot scrolls outside the catenary riser, used for underwater petroleum exploration, and is designed to perform several nondestructive tests. It can also be used to reconstruct the riser profile. Here, a realistic simulation model of robot kinematics and its environment is proposed, using different sources of data: oil platform characteristics, riser static configuration, sea currents and waves, vortex-induced vibrations, and instrumentation model. A dynamic finite element model of the riser generates a nominal riser profile. When the robot kinematic model virtually scrolls the simulated riser profile, a robot kinematic pattern is calculated. This pattern feeds error models of a strapdown inertial measurement unit (IMU) and of a depth sensor. A Kalman filter fuses the simulated accelerometers data with simulated external measurements. Along the riser vertical part, the estimated localization error between the simulated nominal and Kalman filter reconstructed robot paths was about 2 m. When the robot model approaches the seabed it assumes a more horizontal trajectory and the localization error increases significantly.

#### 1. Introduction

One of the key elements of deep-water petroleum exploration is the production riser. Risers are the ducts that transport petroleum, water or gases from the exploitation well up to the production platform. Either rigid or flexible types of risers may be used in the oil field. Both types are submitted to a broad spectrum of failure causes [1]: mechanical loads, aging, corrosion, erosion, temperature effects, installation or fabrication nonconformities, and so forth. Therefore, the availability of inspection tools to assess riser integrity status *in situ* is highly desirable. Such procedures are performed mainly by visual inspection with remotely operated vehicles (ROVs) [2] or autonomous underwater vehicles (AUVs) [3]. In some cases, sensors are installed directly on fixed points of the riser surface to measure strain and riser motion [4]. Other types of nondestructive testing (NDT) techniques can be used, such as magnetic, radiographic, or ultrasound methods [5]. In these cases, however, the operational constraints for using human operators are a major problem. A few papers address robotic devices specifically designed for underwater riser inspection. Psarros and his collaborators [6] proposed a robot that moves along the riser by using a mechanism composed of two parts. One part stays attached to the riser body, and the other part moves towards the riser’s side, in a cyclical manner.

A major technical problem in robotic underwater inspection is the navigation and/or localization of the robot in a highly dynamic sea environment. Navigation is especially critical for AUVs and somewhat critical for ROVs. Lee at al. [7] addressed this problem by using several sensors fused by a multirate Extended Kalman Filter (EKF). The sensors set included a strapdown inertial platform, a Doppler velocity log (DVL), magnetic compass, and a depth sensor. However, they had sonar transducers installed in an underwater reference station and in the remote vehicle. Jouffroy and Opderbecke [8] addressed the problem of measuring the horizontal position of a ROV by using a gyro-Doppler together with an ultrashort baseline (USBL) acoustic positioning system. Diffusion-based observers were used to process a trajectory segment, instead of a typical point-by-point localization. He et al. [9] proposed an approach based on an invariance extended Kalman filter (IEKF) to address the problems of using sonar in shallow waters. In the case when the robot is mechanically linked to the inspected structure, the key problem is to localize precisely where it is at every instant of time. Such localization coordinates are associated to NDTs data and the flaws position can be precisely determined.

Recently, our group designed and built a prototype of a robotic device specifically designed to perform nondestructive testing (NDT) in production risers [10]. The robot has neutral floatability and *embraces* the riser by moving along its outside (Figure 1), using a pair of thrusters for propulsion as well as polymeric wheels to guarantee sliding and correct alignment with the riser surface. In Figure 2 it is possible to observe how the robot attaches the riser by opening and closing its motorized arms. This operation is assisted by a human diver. It communicates with the operator’s computer by means of an umbilical cable that transmits power, images and control commands. The dimensions of the robot and additional parameters used along the work are shown in Table 2. This robot will be able to perform several NDT procedures, such as ultrasound, imaging, and mechanical vibration measurements.

This paper proposes a kinematic model of the robot performing a riser profile cast mission, in a realistic simulated environment. Initially, a riser dynamic profile is estimated using a finite element model of the riser subjected to sea and ship motions. The nominal robot kinematic path (including position, velocity and acceleration), as it scrolls by the riser, is *contaminated* with experimental errors, simulated by IMU and depth sensor models. The simulated sensor data is used by a Kalman filter to estimate the original robot path. This path is a good estimate of the actual riser profile, if robot mission time is small, compared to platform motion.

The obtained profile can be used as an imposed displacement data for some structural analysis software based on finite elements techniques, that allows stress to be calculated. In addition, the localization algorithm can be used to associate each NDT measurement with its riser geometrical coordinate. These two aspects are intimately connected, and the localization algorithm can be used either to cast the profile, for fast robot runs, or localize the NDTs.

Reproducing the expected environmental conditions, to test the proposed approach, in a laboratory experiment is essentially impracticable. Field tests, by his turn, should require a expensive positioning system such as a 3D sonar, which does not operate at the required frequency resolution, due to the presence of vortex induced vibrations (VIVs). Therefore, a simulation of the riser application, together with simulated sensors, was used to assess the performance of the localization algorithm.

Actually, a particular environmental and riser configuration scenario is being addressed in this paper. However, the approach is likely to be applicable to similar situations. Other devices that move along subsea pipe systems, such as flowlines, jumpers, and umbilicals, could employ some of the main ideas presented in this paper. No additional localization devices, such as sonar beacons, are needed.

The localization problem formulation used a standard Kalman filter as a sensor fusion algorithm based on a simple kinematic model of a strapdown IMU fused with a depth sensor. More sophisticated sensor fusion algorithms or state-space models for the system (e.g., dynamical models) could also be tested in future implementations, but the problem formulation would be probably quite similar.

#### 2. Riser Simulation Conditions

The particular sea and ship motion conditions selected for running the simulations corresponded to a *severe* condition, relatively similar to that typically found in Campos Basin, in the southeast of South America coast. Actually, they were designed to be worse than the most severe scenario in which a robot is expected to operate (Table 1). Additionally, under milder sea conditions, the localization performance should be expectable to be better than shown here.

Data from a flexible free-hanging riser installed in a PETROBRAS (Brazilian State oil company) Turret Floating Production Storage Offloading (FPSO) oil platform, which is currently in operation in the Campos Basin, was used as the inputs for the riser simulation software FLEXCOM. This is a finite elements software customized for nonlinear static and dynamic analysis of offshore systems, used worldwide by the petroleum industry from the last 20 years, and validated against experimental tests and other finite elements packages [11, 12]. The software allows riser responses to be simulated with several kinds of platform characteristics, sea current profiles, hydrodynamic loads, regular and irregular waves, and so forth. Two situations were studied: static and dynamic. In the static analysis, only the equilibrium configuration of the riser is considered, without motions other than from the robot itself. Three positions of the platform with relation to the wheel head were considered: standard and with a ship offset of 150 m in the directions *near* and *far*. For the dynamic analysis, sea conditions with regular waves, reaching the ship 45° obliquely, were used to generate the riser motion profiles.

To estimate numerically the localization system performance and calculate the associated displacement errors, the arrangement shown in Figure 3 was used. By simulating the kinematic model (Section 4), a *physical* profile of Euler angles, inertial accelerations, and depth was calculated. These simulated variables are then contaminated by noise, using the IMU and depth sensor instrumentation error models (Section 5), providing realistic sensor outputs. The trajectory run by the robot is estimated by a sensor fusion algorithm (Kalman filter), using such noisy sensor data (Section 6). Finally, both *physical* and estimated riser shapes are compared, to estimate the localization error along the run (Sections 7 and 8).

#### 3. Architecture of the Localization System

The proposed localization system, that will be simulated numerically here, is shown in Figure 4. An IMU measures three accelerations, angular velocities, and Euler angles from the robot as it scrolls along the riser. Accelerations are measured in the local reference frame (see definition in Section 4). Using the Euler angles, the accelerations are transformed to the global reference system, using the classic *strapdown* inertial navigation approach [13]. The simulated measurements from the IMU are fused by the Kalman Filter (KF) with the processed simulated external measurements from the depth sensor. Its output is an estimated state vector, that includes robot position, in the global reference frame.

#### 4. Generation of Static and Dynamic Sensor Profiles

A simulated profile of inertial sensor physical excitation was obtained. First, the *nominal* or *static* riser geometry was obtained from the FEM (finite element model) analysis. The global reference system (, , ) was positioned in the turret center, such that axis points to an arbitrary direction (e.g., north), points up, and is orthogonal to both (Figure 5). This reference frame is attached to the FPSO, being thus a slowly moving frame. However, due to the very low frequency of the FPSO movement compared to robot mission time span, the global reference frame was considered as being inertial.

The local reference system (, , ) is localized in the geometric center of the riser section that coincides with the position, relative to robot height, where the inertial sensor is expected to be installed, and moves together with the robot. The center of the local reference system was defined, relatively to the global reference system, using the same nodes of the FEM mesh, which coordinates in the global reference system are , ( is number of nodes in the mesh and the superscript express the reference system). Each vector is found from a unit vector that connects node to node of the FEM mesh. The negative sign for defining in (2) is due to the fact that points down:

In the sequence, a particular was chosen (along with its orthogonal counterpart ), such that it could represent one plausible mission trajectory. In the present version of the model, the robot communicates with the operator in the platform through an umbilical cable, and only a small amount of spin can be allowed, to prevent the cable from curling along the riser. Therefore, is simultaneously orthogonal to a point sufficiently far in the direction and to −. The vector was orthogonal to both and . For each local reference system , a directors cosine matrix (DCM) was defined as follows: where is the unit vector that defines , and similarly and .

Euler angles , , and were defined between the and , and , and axis, respectively. In every step, these angles were derived from the DCM matrix, by using the formulas presented by [14]. The resulting transformation matrix from local to global reference system was posteriorly verified to be the same as DCM, in order to determine if singularities were present. A fixed point expressed in the local reference frame corresponds to the robot body point where the IMU will be possibly installed. This point was arbitrarily chosen within a reasonable distance from the riser section center, where the origin of the local reference frame is positioned. The coordinates of , expressed in the global coordinate frame, were assumed as the sensor displacement profile.

##### 4.1. Sea Current Effect

A linear current profile from the maximum sea current velocity () at sea level to 0 at seabed [15] was considered as being aligned with the riser catenary plane. was assumed as 1.68 m/s, from Campos Basin data: where is the rotation matrix associated with the catenary angle in the plane, given by the following expression:

The robot can move freely along the riser in the direction, but it was constrained in the other directions because it *embraced* the riser on the outside. If no water current was present, the nominal robot velocity propelled by a pair of thrusters should be in the direction of . However, due to the presence of the current, the absolute velocity of the robot was found by considering the sea current velocity component that is projected over the robot’s trajectory, which can change the robot’s progression velocity:
where is the absolute velocity and the dot product. From a preliminary study, was estimated as 1 m/s. Figure 6 illustrates robot velocity in the particular sea current and thrusters conditions adopted in this paper.

The effect of the longitudinal sheer between the robot and the riser due to the transversal current was not taken into account. This current component was expected essentially to increase the normal force that the robot applied to the outer surface of the riser. Because the riser was tightly fitted among the robot’s rigid free wells to avoid longitudinal and torsional slipping, the increase in the shear force of the wheels that could decelerate the robot was considered to be negligible.

Since all the elements of the FEM mesh have approximately the same length, the time steps are no longer uniformly distributed with such variable velocity profile. The resulting variable time array was calculated by:

This nonuniform time array was inconvenient for future calculations of velocity and acceleration profiles, and a new set of was found by spline interpolation using a uniform time array with the same limits of . In the sequence, the local reference systems , , , the DCM matrix and the Euler angles were recalculated using this new set of coordinates, which were sampled by a uniform time array. All profiles were resampled at 5 Hz by spline interpolation, such that the high frequency riser motions could be followed (next subsection). The total mission time was , and a total of points were generated.

##### 4.2. Riser Motion Effect

The dynamic FEM analysis allows finding, along the time, the altered geometry of the riser, with respect to the nominal static profile. The finite elements model had 170 beam elements with flexion, axial, and torsional deformations. Each element had a nominal length of 10 m. Dynamic FEM analysis provided the , , and coordinates of only ten nodes with respect to time at sampling rate. The analysis was run for a 100 second time span, but the first 50 seconds were disregarded to avoid the transient effect of the FEM solution. The second half of the time window was then replicated to reach the total mission time. Thus, the offset, phase, and amplitude parameters were preserved. This adjustment provided a matrix of ten lines (one for each node) and columns, which was resampled by successive spline interpolations. The first interpolation reduced the number of columns to , the second expanded the lines up to the original FEM mesh (171 nodes) and the third resampled the lines again up to . Therefore, three perturbation matrices , , and were obtained, one for each coordinate , , and . In these matrices, each row was a particular riser deviation from the nominal profile and each column was a time step. The three components of a displacement vector were defined for each node :

This vector was decomposed into its normal and tangential parts. Because the robot was free to move along the riser, only the normal component of vector was effectively transmitted to the robot:

Therefore, the new coordinates , , and (in the local reference frame) of the riser path were given by the following sequence of points, ,

##### 4.3. Vortex-Induced Vibration (VIV) Effect

The sea current passing through a circular cylinder produces vortex-shedding in the wake, which causes the structure to vibrate. This complex fluid-structure interaction is called Vortex Induced Vibration (VIV), and it occurs predominantly on the cross-flow direction [16]. Simulating this effect is an arduous numerical problem. In this paper, we used experimental data obtained from a scale-model available at the Open-Source VIV Data Repository of the Center for Ocean Engineering at MIT [17]. Cross-flow displacement data that is available for tests performed in a bare cylinder, 20 mm diameter and 10 m length, which was donated by ExxonMobil, was used in our simulations. Several fluid velocities and both linearly sheared and linear flow conditions may be used. We chose a *strong* condition of regular flow, of approximately 1 m/s, which provided greater displacements compared to the sheared flow for the same nominal velocity.

To adapt the experimental data to the riser that was being analyzed, the displacement was scaled by the test riser diameter and multiplied by the actual riser diameter. The frequency of shedding () in cylinders with cross-flow is given by the following: where is the Strouhal number, is the fluid velocity, and is the cylinder diameter. Keeping fixed, the ratio of between the experiment and the riser was 14.78. This factor was used to scale the time vector, such that the frequency vector of the data spectrum, which was used to simulate riser VIV, was divided by this quantity. The set of ten points along the cylinder where displacements were measured was associated with the closest nodes of the FEM model. A window of data, without transient effects, was replicated ten times until the total mission time was achieved, similar to the previous section, and the same interpolation procedure was applied. Finally, the displacement caused by VIV was rotated to become perpendicular to the catenary plane and added to (10).

Figure 7 shows the RMS profiles of the riser displacements in the , , and directions as a function of the normalized riser length. The figure also shows the RMS of VIV perpendicular to the catenary plane. The number 0 is for the wheel head, and 1 corresponds to the turret.

##### 4.4. Sensor Velocity and Acceleration Profiles

To simulate IMU output, the expected physical acceleration at the sensor installation point must be found. This acceleration was used to feed the sensor model to find realistic sensor signals. The acceleration at point expressed in the global reference frame was given by the following well-known kinematic equation:

The local acceleration and Coriolis terms are zero because the sensor is fixed in the robot body. is the angular acceleration and is the angular velocity. Angular velocity expressed in the global reference frame was calculated by the methods used in [14, 18]:

To solve (12) and (13), the time derivatives of displacement and rotation , , , , and so forth were found by finite differences and considering all variables in the global reference system.

#### 5. Instrumentation Model

A low-cost microstrain 3DMGX-1 IMU was the project choice. This is a compact and integrated device, suitable for a high-depth submarine application, where the electronic case must be as slender as possible, for mechanical structural reasons. It delivers 3D accelerations, angular velocities, and attitude/orientation matrix in a single-serial channel. The error characteristics of each output were modeled as a wide-band noise plus a first order moving bias Markov process [19]:
where: simulated *corrupted* sensor output,: simulated *clean* physical signal,: Scale Factor Error,: Offset Error,: moving bias, 1st order Markov process,: wide-band sensor noise.

The wide-band sensor noise was defined as , (white noise, zero mean). is the wide-band noise variance. Scale factor error, 0.5% by sensor specifications, was considered as follows:
such that is a uniform random variable in and is the *signum* function. This equation introduces a random error limited to 0.5% in the sensor output signal amplitude, when substituted in (14).

The moving bias was found by integrating the following finite difference equation with the Euler method: where: sampling period (s),: sampling frequency (Hz),: bias input noise variance,: random walk variance,: time constant (s).

##### 5.1. Determination of Sensor Parameters

An experiment was performed to collect the error characteristics by keeping the sensor stationary in laboratory while recording signals. Total acquisition time was 4 h 21 min, with a 75 Hz sampling frequency, after the internal temperature was stabilized. The first determined parameter was the time constant . Two techniques were employed: the Allan variance plot [20] and the autocorrelation function [19]. Both techniques gave different results (see Table 2). According to the likeness between the real and the simulated signals, one of the two time-constant techniques was selected. For the attitude sensor, determined by Allan variance was chosen, whereas for acceleration and angular velocity, autocorrelation provided the best results.

Except for angular velocity (AngRate), the experiment time series was low-pass filtered before calculating , to better characterize the low-frequency component of the error. For , the high-frequency part was used. The cut-off frequencies were chosen by trial and error, and the agreement between the simulated and original error signals was observed for each case. Time series, Allan variance plot [21] and fast Fourier transform (FFT) were used to test the agreement between original and simulated error signals. In the case of attitude simulation, the Euler angles were discretized in the interval comprising to , with steps of , corresponding to sensor resolution. For acceleration, the resolution of the sensor was 0.2 mg, and the angular velocity was . However, amplitude discretization was not implemented. Figure 8 shows the results for Pitch angle and Figure 9 shows the results for acceleration in the direction. For other angles, directions and angular velocities, the results were similar.

The water pressure-based depth sensor specified for this robot was the Digiquartz 8CB4000-I Depth Sensor (Paroscientific, Inc., Redmond, WA), provided a 0.01% accuracy and a resolution of m. In the simulations, we considered that the depth measures were corrupted by the addition of a Gaussian noise with = 0.2 meters, according to the estimates made by Jalving [22].

#### 6. Localization Algorithm

The localization problem consisted of estimating a set of states, which included robot position, by reading accelerations as inputs and considering the depth sensor signals as the external measurements. Since a strapdown navigation scheme is used, the localization algorithm assumes that the accelerations measured in the local frame are transformed into the global frame, by using the rotation matrices calculated from Euler angles measured by the IMU. It is also assumed that if a DC-acceleration sensible accelerometer is being used, the vertical acceleration was compensated for gravity. The problem state equations are formulated as [23, 24]:

The state vectors and matrices are given by the following:

is an identity matrix, such that the terms (1,4), (2,5), and (3,6) are equal to ; (4,7), (5,8), and (6,9) assumes the value ; (1,7), (2,8), and (3,9) elements are . is a matrix of zeros except for the terms (1,1), (2,2), and (3,3) that are equal to ; (4,1), (5,2), and (6,3) elements are .

is a matrix of zeros, with the exception of the (i,i) terms equal to 1. In such expressions, is the sample number, is the state vector, the input vector, is the process noise, is the position, is the velocity, is acceleration, is the acceleration error, is sampling period, is the external measurements vector, is the sensor noise, and is the variance of each associated variable. The observability matrix of this system had a full-rank, that is, 9.

An estimate of the state vector is found by a Kalman filter.

Prediction:

Update:
where is the state estimate error covariance matrix, is the sensor noise covariance matrix, is the process noise matrix, and is the Kalman gain. This is a standard, simple and generic implementation of a Kalman filter estimator. More sophisticated approaches could be tested as well, although the results obtained with the standard KF were satisfactory for the intended application, as will be addressed in Section 7. Unlike the kinematic models that are usually used for land robots [25] or those built to operate over flat surfaces [26], the state-space model used here to formulate the KF does not use any *a priori* information of the vehicle characteristics. Nor the dynamical characteristics of the sensors are included in this KF formulation. Thus, in principle, the localization algorithm proposed here could be applied for similar applications with different instrumentation, robot, and environmental characteristics. However, some tuning work on the covariance matrices should likely to be necessary.

##### 6.1. Processing of the External Measurements for the KF

Here, we focus the analysis on the key aspect of this particular localization problem, which is determining the external measurements vector . , and so forth, are the coordinates of the robot position in the global reference frame as measured by an independent absolute sensor. The sensor, that is, in practice, available to perform such measurements is the depth sensor, which can only measure the coordinate. The attempts to use this measurement exclusively in the Update phase of the KF gave nondrifting estimates in the vertical direction, as expectable, but not in the horizontal one.

We proposed the following method to estimate the complete vector using data from , (depth sensor signal), and the Euler angles measured by the IMU. Considering Figure 10, the vector joining to , in the local reference frame centered in , is given by . This vector lies in direction of the local reference frame , with modulus , and coinciding with the instantaneous riser path. In the global reference frame, can be found as:

Thus,

The same difference calculated in (23) can be expressed by a vector of global coordinate differences:

Since and are known, by comparing (23) and (24), it is possible to find as:

Substituting in (23), the KF external measurements vector is found. Alternatively, the displacement could be measured by an odometer. However, in this case, special constructive care should be taken into account to prevent slippage. In any case, this approach assumes a previous observation in the Correction phase of the KF. However, KF theory assumes conditional independence among the observations. Therefore, only *suboptimal* estimation performance is likely to be expected.

#### 7. Simulation Results

##### 7.1. Static Results

In the static case no riser motion is present. The *standard* position was used as the basis to generate the dynamic analysis and is an intermediary *neutral* configuration between the *near* (ship displaces towards wheel head) and the *far* (the opposite direction) positions. Figure 11 shows the robot trajectory and attitude, represented by the successive local reference frames, for the neutral position. In Figure 12, the associated localization errors are shown, as a function of distance from the seabed (using a grid of 1.680 points, or 1.44 Hz). The absolute errors associated with the neutral, far, and near configurations are shown in Figure 13, in a *semilog* scale. To study the impact in the localization error of the grid refinement grade, an additional error curve was generated, for the neutral configuration, with a 7,630 points (5 Hz) grid, the same used for the dynamic analysis.

Figure 12 shows that the greatest source of error occurs in the direction, while in the deepness coordinate , the error is low. In addition, the total error increased as the robot tilted horizontally. This trend was expected due to the relative loss of accuracy of the depth sensor to estimate the parameter when using (21)–(25). When the depth measure difference between the two samples is small, a decrease in signal to Gaussian noise ratio is observed. From an application point of view, the robot should not work too close to the point where the riser periodically touches the seabed, the touch down point (TDP). This safety margin is established as being of 50 m. At this point, the expected positioning error for the static case was approximately 10 m. The error was slightly higher in the far case (Figure 13), where the riser was more horizontal than in the near condition. By increasing the sampling frequency (Figure 12), the error was likely to decrease in most of the path, but when the robot approximates to the TDP, the localization system fails.

##### 7.2. Dynamic Results

By including the wave and VIV effects in the riser nodal displacement profile, the robot trajectory becomes more complex (Figure 14), as expectable. The associated localization errors are shown in Figure 15, using the 5 Hz mesh. Figure 16 shows a zoom in the plane, close to TDP, of both the robot trajectory and the external measurements, calculated with (21)–(25), as well as the Kalman filter position states estimations. The KF estimated the trajectory accurately, up to the point where the processed external measurements lost accuracy, below 15 m distance to seabed. The KF estimates essentially followed the external measurements, in part because the KF Update frequency was the same as the Prediction. By comparing the errors in both static and dynamic cases (Figure 17), analogous profiles were found. It means that the Kalman filter was able to closely track the robot’s trajectory, even in the dynamic case, and the loss of algorithm localization performance was mainly caused, after all, by depth sensor differential signal to noise ratio.

#### 8. Discussion and Conclusions

The simulations show that the robot trajectory using the localization algorithm, in the dynamic case, presented a wavy pattern. The pattern corresponds to the path that is traveled by the robot, and not the nominal riser profile, as expected. If an estimate of the riser catenary shape is desired, the obtained path could be used to fit a smooth profile curve.

In the more vertical part of the riser (above 15 m from seabed), the average estimated error (standard dev.) was 0.76 (0.47) m (Figure 15). This error can be considered sufficiently small for localizing the riser in the oil field and to feed nodal displacement constraints in a structural Finite Elements tension analysis. This accuracy is also satisfactory for localizing flaws detected by the NDTs.

The localization algorithm depends strongly on the external measurements, which are provided by the depth sensor and the IMU, using (21)–(25). However, close to seabed, the simulated robot trajectory tilts towards a more horizontal attitude, and following the riser section which is closer to the TDP. In this part of the trajectory, Gaussian noise from the depth sensor more significantly corrupts the relative depth measurements. This effect increases when a period between two successive measurements becomes shorter. The 5 Hz is the minimum estimated rate necessary to reproduce VIV and the wave motion data. On the other hand, the minimum safety distance between the robot and the seabed coincides with the condition which the localization algorithm loose accuracy. As a possible extreme limit, at 15 m from the seabed, the expected error should be also around 15 m (Figure 15).

In our opinion, the simple and standard KF implementation that was presented in this study is sufficient to provide the localization accuracy required for the proposed application. In the case of rigid risers and other predominantly vertical offshore structures, the observed loss of accuracy is not expected to occur. The kinematic model used to formulate the sensor fusion problem was linear, and therefore the standard KF implementation is appropriate. Other models could be proposed, incorporating robot, riser or sea dynamic features. In such cases, using an extended Kalman filter or a particle filter should be necessary for handling the associated nonlinearities. In any case, the external measurements reliability seems to be key aspect for localization accuracy in this problem.

Using a higher grade low drift IMU instead of a low-cost one could eliminate the need for external measurements in a critical region but at the price of increasing payload, cost and volume. Acoustic positioning systems could be used to localize the robot in the most horizontal parts of the trajectory. However, these systems work at a very low sampling frequency (0.1 to 1 Hz) [8] and, at best, they could only probably cast the external motion envelope of the robot. More practically, improved accuracy could be achieved in the section close to the TDP by using an odometer to estimate in (25). As such, the mechanical setup of an odometer should be accomplished very carefully to prevent slippage. Preinstalled RFID (Radio Frequency Identification) [27] tags or visual marks along the riser, signaling the actual length, are possibilities that could be explored. Switching the depth sensor with other kinds of sensors for external measurements, close to TDP, could be performed by using a fuzzy expert system [28]. In future studies, experimental laboratory and field tests should be performed, using an acoustic localization system to assess the bounds of the navigation errors estimated in this work.

#### Acknowledgments

The authors are gratefully thankful to CAPES (Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior), FINEP (Financiadora de Estudos e Projetos), CNPq (Conselho Nacional de Desenvolvimento Cientifico e Tecnologico), and FAPERJ (Fundacao de Amparo a Pesquisa do Estado do Rio de Janeiro) for financial support.

#### References

- P. O’Brien and and J. Picksley, “State-of-art flexible riser integrity issues,” Tech. Rep. 2-1-4-181/SR01, 2001, UKOOA by MSC, Rev. 04. View at Google Scholar
- G. Chapin, “Inspection and monitoring of girassol hybrid riser towers,” in
*Proceedings of the Offshore Technology Conference (OTC '05)*, Huston, Tex, USA, 2005, OTC paper no. 17696. - K. Hamilton and J. Evans, “Subsea pilotless inspection using an autonomous underwater vehicle (SPINAV): concepts and results,” in
*Proceedings of the OES/IEEE Oceans (OCEANS '05)*, pp. 775–781, Brest, France, June 2005. View at Publisher · View at Google Scholar · View at Scopus - M. Ordonez, M. O. Sonnaillon, D. Murrin, N. Bose, and W. Qiu, “An advanced measurement system for vortex-induced-vibrations characterization in large-scale risers,” in
*Proceedings of the MTS/IEEE Conference Oceans (OCEANS '07)*, pp. 775–781, Vancouver, BC, Canada, 2007. View at Publisher · View at Google Scholar - E. Veith, C. Bucherie, J. L. Lechien, J. L. Larrouse, and B. Rattoni, “Inspection of offshore flexible riser with electromagnetic and radiographic techniques,” in
*Proceedings of the 15th World Conference on Nondestructive Testing*, pp. 404–408, Rome, Italy, 2000. - D. Psarros, V. A. Papadimitriou, P. Chatzakos, V. A. Spais, and K. Hrissagis, “A service robot for Subsea flexible risers: analysis and systematic design,”
*IEEE Robotics and Automation Magazine*, vol. 17, no. 1, pp. 55–63, 2010. View at Publisher · View at Google Scholar · View at Scopus - P.-M. Lee, B.-H. Jun, K. Kim, J. Lee, T. Aoki, and T. Hyakudome, “Simulation of an inertial acoustic navigation system with range aiding for an autonomous underwater vehicle,”
*IEEE Journal of Oceanic Engineering*, vol. 32, pp. 327–345, 2007. View at Google Scholar - J. Jouffroy and J. Opderbecke, “Underwater vehicle navigation using diffusion-based trajectory observers,”
*IEEE Journal of Oceanic Engineering*, vol. 32, pp. 313–326, 2007. View at Google Scholar - C. He, E. Jorge, and L. M. Zurk, “Enhanced kalman filter algorithm using the invariance principle,”
*IEEE Journal of Oceanic Engineering*, vol. 34, no. 4, pp. 575–585, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. F. Santos, L. L. Menegaldo, and M. O. Brito, “An outer device for universal inspectionof risers,” US Patent Application based on Brazilian PI 0705113-1, 2008.
- R. O’Grady, H.-J. Bakkenes, D. Lang, and A. Connaire, “Advancements in response prediction methods for deep water pipe-in-pipe flowline installation,” in
*Proceedings of the Offshore Technology Conference (OTC '08)*, Houston, Tex, USA, 2008, OTC paper no. 19400. - P. J. O'Brien and J. F. McNamara, “Significant characteristics of three-dimensional flexible riser analysis,”
*Engineering Structures*, vol. 54, no. 4, pp. 223–233, 1989. View at Google Scholar · View at Scopus - D. H. Titterton and J. L. Weston,
*Strapdown Inertial Navigation Technology*, Edited by Herts, IEEE, London, UK, 2nd edition, 2004. - B. Siciliano, L. Sciavicco, L. Villani, and G. Oriolo,
*Robotics: Modelling, Planning and Control*, Springer, London, UK, 3rd edition, 2009. - A. L. S. Pinho,
*Tension reduction in TLP platforms rigid risers*, M.S. thesis, Civil Engineering Department, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil, 2001. - C. L. Cunff, F. Biolley, E. Fontaine, S. Etienne, and M. L. Fancchinetti, “Vortex-inducedvibrations of raisers: theoretical, numerical and experimental investigation,”
*Oil and Gas Science and Technology (IFP)*, vol. 5, pp. 59–69, 2002. View at Google Scholar - M. I. of Technology, “Vortex induced vibration data repository,” 2011, http://oe.mit.edu/VIV/.
- H. I. Weber, “Inertial Systems: basic foundations,” in
*Proceedings of V SBEIN 5th Brazilian Symposium on Inertial Engineering, Tutorial A, (in Portuguese)*, Rio de Janeiro,Brazil, 2007. - W. S. Flenniken IV,
*Modeling inertial measurement units and analyzing the effect of their errors in navigation applications*, M.S. thesis, Auburn University, Auburn, Ala, USA, 2005. - D. W. Allan, “Statistics of atomic frequency standards,”
*IEEE Proceedings*, vol. 54, pp. 221–230, 1966. View at Google Scholar - IEEE, “IEEE Recommended Practice for Inertial Sensor Test Equipment, Instrumentation, Data Acquisition, and Analysis,”
*IEEE Std*1554-2005, IEEE Aerospace and Electronic Systems Society, 2005. View at Google Scholar - B. Jalving, “Depth accuracy in seabed mapping with underwater vehicles,” in
*Proceedings of the OES/IEEE Oceans (OCEANS '99)*, pp. 973–978, Seattle, Wash, USA, September 1999. View at Scopus - D. D. S. Santana,
*Terrestrial trajectory estimation using a low cost measurement unity and sensor fusion*, M.S. thesis, Mechatronic Engineering Department, Polytechnic School, University of Sao Paulo, Sao Paulo, Brazil, 2005. - H. Qi and J. B. Moore, “Direct Kalman filtering approach for GPS/INS integration,”
*IEEE Transactions on Aerospace and Electronic Systems*, vol. 38, no. 2, pp. 687–693, 2002. View at Publisher · View at Google Scholar · View at Scopus - T. Sasaki, D. Brščić, and H. Hashimoto, “Human-observation-based extraction of path patterns for mobile robot navigation,”
*IEEE Transactions on Industrial Electronics*, vol. 57, no. 4, pp. 1401–1410, 2010. View at Publisher · View at Google Scholar · View at Scopus - L. L. Menegaldo, G. A. N. Ferreira, M. F. Santos, and R. S. Guerato, “Development and navigation of a mobile robot for floating production storage and offloading ship hull inspection,”
*IEEE Transactions on Industrial Electronics*, vol. 56, no. 9, pp. 3717–3722, 2009. View at Publisher · View at Google Scholar · View at Scopus - S. S. Saad and Z. S. Nakad, “A standalone rfid indoor positioning system using passive
tags,”
*IEEE Transactions on Industrial Electronics*, vol. 58, no. 5, pp. 1961–1970, 2011. View at Google Scholar - S. H. P. Won, F. Golnaraghi, and W. W. Melek, “A fastening tool tracking system using an IMU and a position sensor with Kalman filters and a fuzzy expert system,”
*IEEE Transactions on Industrial Electronics*, vol. 56, no. 5, pp. 1782–1792, 2009. View at Publisher · View at Google Scholar · View at Scopus