Automation, Control and Intelligent Systems
Volume 4, Issue 3, June 2016, Pages: 59-65

A Deterministic Approach to Measurement of Noise Attenuation in Oil-Rig Drill Ship Positioning Systems

E. C. Obinabo1, T. C. Nwaoha2, *, F. I. Ashiedu3, C. O. Izelu3

1Department of Electrical and Electronic Engineering, Ambrose Alli University, Ekpoma, Edo State, Nigeria

2Department of Marine Engineering, Federal University of Petroleum Resources, Effurun, Delta State, Nigeria

3Department of Mechanical Engineering, Federal University of Petroleum Resources, Effurun, Delta State, Nigeria

Email address:

(T. C. Nwaoha)

*Corresponding author

To cite this article:

E. C. Obinabo, T. C. Nwaoha, F. I. Ashiedu, C. O. Izelu. A Deterministic Approach to Measurement of Noise Attenuation in Oil-Rig Drill Ship Positioning Systems. Automation, Control and Intelligent Systems. Vol. 4, No. 3, 2016, pp. 59-65. doi: 10.11648/j.acis.20160403.12

Received: May 30, 2016; Accepted: June 12, 2016; Published: June 29, 2016


Abstract: This paper presents a qualitative evaluation of wave-induced motions in an oil-rig drill ship positioning system which incorporates a priori knowledge of noise contamination in the measured data. The noise contamination β defined in the function of the known form (P(X, β)) and X takes the specific values z, which from Cramer-Rao bound, gives the smallest possible variance with which the estimate of β can be determined. A conceptual model of the problem based on the maximum likelihood techniques in terms of joint probability distribution functions enhanced convergence of the iteration process. A filter was postulated to define the error covariance matrix which yielded unbiased estimates of the measured data.

Keywords: Wave-Induced Motions, Oil-Rig Drill Ship, Maximum Likelihood, Error Covariance Matrix, Noise Attenuation, Filter, Minimum Variance Estimate, Error Covariance Matrix, Iteration Process


1. Introduction

Noise attenuation in oil-rig drill ship positioning systems has been addressed severally by recourse to filter design [1] [2] and control system modelling assuming all the noise processes are independent, and the filter designed to give a minimum variance estimate. This reduces the problem to one of optimal control whereby the best estimate using the measured values of the input and output of the system are required. Because measurements invariably contain errors [3] [4] [5]; the approach to the problem utilizes concepts of probability and statistics in which state estimation is addressed because noise is usually known to be correlated with the measured data [6] [7].

In an exploration method by acoustic means, [8] employed digital filtering and signal processing to establish a guaranteed data interpretation of the liquid level position in noisy well bores. Signal filtration in this application provides additional processing techniques under operator control to obtain accurate results in oil wells with shallow liquid levels and noisy well bores. This is enhanced by inclusion of an analogue filter [8] [9] to reduce the interference from signals generated by the reflections at the tubing collars, and to ensure that the liquid level signal could be detected in the majority of the cases. Here, the fundamental objective is markedly enhanced by increasing the signal-to-noise ratio in the recorded data. One method of doing this is to filter measured values to reject undesirably high frequency components [8] [9]. This is achieved by passing the measurements through simple low-pass filters before connecting to the input. The cut-off frequency of the filter would then be selected in conjunction with the sampling and conversion rates in order to satisfy a given sampling requirement while preserving the signal components of interest.

In digital signal processing [8] [9] [10] [11] [12], the analogue-to-digital converter (ADC) gives high resolution of the data based on continuous integration and oversampling in combination with carefully designed low-pass filters (fig. 1) which eliminate gaps in the sampling process. The negative feedback from the output (fig. 1) along the  digital-to-analogue converter (DAC) and the sum amplifier at the input are performed at a high sample rate, transferring the quantization noise into the stop band of the digital low-pass filter, for decimation. In the process, the sampled analogue signal is fed to the sum amplifier along with the output of the  digital-to-analogue converter (DAC). The integrated difference signal is fed to the strobed comparator whose output samples the difference signal at the sampling frequency many times that of the analogue signal frequency. The output of the comparator provides the digital input for the  digital-to-analogue converter (DAC). Thus, the system functions on a negative feedback loop which minimizes the difference signal by tracking the input. The integrator is continuously fed with the differential signal and there are no gaps in the analogue input signal as introduced by sample and hold devices. The digital information representing the analogue input voltage was coded in the polarities of the pulse train appearing at the output of the comparator which can be retrieved as a parallel binary data word applying a digital filter operator.

Fig. 1. A Continuous Sampling Signal Modulator [8].

The Kalman filter has received considerable attention in the existing literature [2] [13] [14] [15] [16] and has been applied in wave-induced motions of ships [17]. Given a signal model that consists of a linear dynamic system driven by stochastic white noise processes, the Kalman filter [2] [18] exploits a state space model for optimal filtration of noisy measurements. A precursor to the Kalman filter was the Weiner filter which was derived independently by Weiner and Kolmogorov [19] and which gives a method of optimally attenuating noise in process measurements. However, the Weiner filter is limited to time-invariant problems involving stationary noise sequence. the purpose of this study was to develop a statistical quality control scheme for measurements of position and heading of a ship in wave-induced motions, and to show that a filter can be postulated to define a covariance matrix to yield unbiased measurements.

2. Formulation of the Problem

Generally, estimation of process variables contaminated with noise is formulated on the basis of maximum likelihood [17] [2] using statistical information in terms of joint probability distribution functions. When the additive, zero-mean white Gaussian measurement noise is defined in terms of mean values and variances, which will be appropriate for many practical applications, the least-squares solution is formulated as a deterministic problem [20] [21] [22] with appropriate weighting that leads to the maximum likelihood estimate. the ship dynamics considered in this study were assumed to be completely controllable and observable. The autopilot maintains the ship on a set or desired heading while being subjected to disturbances such as wind, waves and current. The main elements of the autopilot system are shown in fig. 2. The actual heading is measured by a gyro-compass (or magnetic compass in a smaller vessel), and compared with the desired heading, dialed into the autopilot by the ship’s master. The autopilot, or controller, computes the demanded rudder angle and sends a control signal to the steering gear. The actual rudder angle is monitored by a rudder angle sensor and compared with the demanded rudder angle, to form a control loop.

Fig. 2. Ship Autopilot Control System [23].

The rudder provides a control moment on the hull to drive the actual heading towards the desired heading while the wind, waves and current produce noise moments that may help or hinder this action. The block diagram of the system is shown in fig. 3.

Fig. 3. Block Diagram of Ship Autopilot Control System [23].

For the purpose of this study, the ship dynamics were represented by a linear discrete-time process using the difference equation [24] [1].

, (1)

where  is an r-dimensional random vector process with the following statistical properties

(2)

(3)

(4)

The analytical approach which often has some computational advantages is to derive a deterministic equation whose solutions are the desired mean and covariance matrix of . Taking the expectation of both sides of (1) yields

, (5)

This is a deterministic linear difference equation which yields as its solution the expected value of . To derive the covariance of , we subtract (5) from (1) to obtain  and then postmultiply both sides of the resulting equation by  and then take the expectations of both sides.

3. Maximum Likelihood Estimate

Estimation of the process considered in this study was based on the assumption that some or all of the parameters may be unknown even though the structure of the differential equation characterizing the system as well as the initial and boundary conditions may be available. This reduced the problem to one of optimal control whereby the best estimate using the measured values of the input and output the system are required. Because measurements invariably contain errors, solution of the problem should utilize concepts of probability and statistics in which the problem should also address state estimation because noise is usually known to be correlated with the measured data. This must be estimated at the same time as the parameters during which filtration is mandatory in the processing of the data. Computation of the optimal estimates should consequently rely on convergence of the iteration employed, which is accomplished through sequential filtration of the estimate.

We considered a time-series model of the form [25]

(6)

The partial fraction expansion required for obtaining a time solution to (6) is defined for distinct roots using the method of residues, as

(7)

where  are the poles of the function  with multiplicities .

The inversion integral was then written in the form

, (8)

so that the function  with contaminated noise becomes

(9)

where ,  are the parameters of the system, and n is the order of the system

Using the shift operator  defined by , equation (6) becomes

, (10)

where  is the observed output signal,  is the applied input signal, N is the number of samples and  is the noise sequence. The polynomial operators A, B, and C are defined as follows

(11)

If the continuous model in (4) is discretized we shall obtain a discrete time model of the form

(12)

where the parameter set a, b is related to the parameter set ,  via the relation

and

(13)

where T is the sampling interval in seconds. The infinite noise of the discrete observations due to aliasing that may result from the above discretization process is negligible. This is justified if  and  are independent of all k and s. this is a reasonable assumption as long as the identification is performed for data acquired from experiments, where  is a priori known sequence. In some practical situations, this assumption is often violated when operating records are used because in such a case the input may depend on the output through feedback.

The canonical model in (10) was made equivalent to the discrete-time model of (12) by satisfying the conditions

(i)

(ii)             (14)

The statistical method of maximum likelihood was employed to optimize the probability of obtaining the expected result. Consequently, a loss function  was defined as follows

(15)

which will be minimized with respect to the system parameter set . The residues were defined by

(16)

so that the values of the parameter set a’ and b’ that make  in (15) minimum will be the estimates of the parameters of the system. The approach considered is one of finding the coefficients of the prediction model

(17)

so that the mean square prediction error

(18)

is as small as possible. By doing so the assumption of Gaussian distribution of the noise sequence  may be relaxed. In the model given in (12), the residues were obtained as

(19)

For system order greater than n,  as given in (19) becomes computationally difficult and highly nonlinear, whereas in the model of (12), the residues given by (16) may be computed quite easily for any given system order. Since the two models are equivalent via the transformations defined in (8), we are free to use any of them. Thus the form of residues given in (16) is recommended and used in this study for parameter estimation; application of a filter for signal noise attenuation becomes evident.

Now, if there exists an unknown parameter  defined in the function of the known form  and  takes the specific values , Cramer-Rao bound gives the smallest possible variance with which an estimate of  can be determined. Now, let  be an unbiased estimate of the given function . Since it is unbiased, we write

(20)

i.e.,

(21)

since  was drawn from the process which generated , it has the same probability density function

(22)

Re-writing (12)

(23)

Differentiating with respect to  gives

(24)

Since  is a probability function

(25)

and its differential with respect to  is zero, that is,

(26)

then

(27)

Substituting (21) from (24) gives

(28)

We now consider the log-likelihood function as follows:

(29)

(30)

Substituting in (25) gives

(31)

Now from, the Schwartz inequality

(32)

where

(33)

and

(34)

Thus for any unbiased estimate the lower bound on the variance of the estimate is given by Cramer-Rao bound as follows:

(35)

If  is an unbiased estimate of , i.e.,

Therefore

(36)

(37)

4. Application

In an oil-rig drill ship positioning system which incorporates a priori knowledge of noise contamination in the measured data, it is required to obtain an estimate of the mean value of the position of the ship based on independent measurements  with the time between successive data acquisition to be distributed with a probability function of the form

(38)

which describes the probability that the parameters  and  caused the measurements  to occur.

Here, we obtained an estimate of the mean position of the ship and expressed it as a product of the individual functions:

(39)

For covariance, we define the log-likelihood function

(40)

which we maximized with respect to  by equating to zero the partial derivatives of  as follows:

(41)

and

(42)

Maximizing with respect to  and  gives

(43)

(44)

Letting  gives the maximum likelihood of the estimates as follows

(45)

(46)

From the Cramer-Rao bound

(47)

we obtained, for

(48)

so that

(49)

showing that the estimate is consistent with  as .

The wave-induced motions of a ship considered in this paper was concerned essentially with measurements which are correlated with coloured noise assumed to be additive to the output.

Fig. 4. Block Diagram of Ship Dynamics System.

The modelling of (19) was enhanced using the discrete-time stochastic sequence as follows

(50)

Multiplying through (50) by :

(51)

i.e.,

(52)

where  is the correlated noise sequence (53)

In matrix form we write (50) as follows

(54)

Using least squares:

(55)

(56)

Taking the expected values of both sides of (50) gives:

(57)

Since  depends on  then  and e are correlated.

Thus

(58)

Therefore the least square estimates of  are biased.

From (55),

(59)

An unbiased estimate of  can only be obtained if the noise sequence  is reduced to a white noise sequence. In order to reduce the noise sequence to a white noise sequence, we considered (51) and postulated a filter F [19] [17] such that

(60)

We then obtained:

(61)

Now, we chose  so that (61) becomes:

(62)

thereby reducing the additive noise  to a white noise sequence . The least square estimate  is now unbiased. Fu and Fy are filtered inputs and outputs.

(63)

5. Conclusions

A stochastic time-invariant model has been developed for linear filtering and prediction of state in wave-induced motions of a ship. The model provides a processing technique under operator control for eliminating the interference from signals generated by the wave-induced motion. The result presented has shown that the wave-induced motions of an oil-rig drill ship positioning system which incorporates a priori knowledge of noise contamination in the measured data can be estimated using the maximum likelihood techniques in terms of joint probability distribution functions. The ship dynamics itself is represented by a stable, linear, parameter-dependent state model.


References

  1. Grimble M. J. (1981). Adaptive Kalman Filter for Control Systems with Unknown Disturbance IEE Proc. 128, pt. D, 6, 263–267.
  2. Kalman, R. E. (1960). A new approach to linear filtering and prediction problems, Trans. ASME J. Basic Eng., 82, 35–45.
  3. Ekejiuba, A. O. M. (1986). Locating By-passed Hydrocarbon Accumulation in Old Reservoirs, Nigerian Association of Petroleum Explorationists, Fourth Annual Conference.
  4. Ahonsi, B. and Gdula, J. (1993). A Case History of By-passed Oil in the Niger Delta, Society of Petroleum Engineers, Nigeria Council, 35–46.
  5. Weber, K. J. (1971). Sedimentological Aspects of Oil Fields in the Niger Delta, Goelogie en Mijnbouw Sc., 559–576.
  6. Dalley, R. M., Gevers, E. C. A. Stamfli, G. M. Davies, D. J. Gastaldi, C. N. Ruijtenberg, P. R. and Vermeer, G. I. O. (1989). DAZIL (Dip Azimuth, illumination), First Break, 7 (3), 86–100.
  7. Esho, K. M. (1993). Prediction of Remaining Oil in the Imo River D2.0 Sand Using 3D Seismic Data, Society of Petroleum Engineers, Nigeria Council, 23–28.
  8. Van Wagoner, J. C., Mitchum, R. M. K. Campion, M. and Rahmanian, V. D. (1990). Siliciclastic Sequence Stratigraphy in Well Logs, Core and Outcrops: Concept for High-Resolution Correlation of Time and Facies, AAPC Methods in Exploration Series, No. 7.
  9. Weber, K. J. (1986). Hydrocarbon Distribution Patterns in Nigeria Growth Fault Structures Controlled by Structural Style and Stratigraphy, Journal of Petroleum Science and Engineering, 1, 91–104.
  10. Benato, G., D’Andrea, V., Cattadori, C., and Riboldi, S. (2015). Improvement of the GERDA Ge detectors energy resolution by an optimized digital signal processing, Physics Procedia, 61, 673-682.
  11. Lépy, M. C., Cissé, O. I., and Pierre, S. (2014). Comparison of digital signal processing modules in gamma-ray spectrometry" Applied Radiation and Isotopes, 87, 402-406.
  12. Seyfried, D., Brueckner, S., and Schoebel, J. (2014). Comparison of antenna dispersion and digital signal processing effects in ultrawideband ground penetrating radar systems, Journal of Applied Geophysics, 101, 20-26.
  13. Hsiao C. H. and Wang, W. J. (2000). State Analysis and Parameter Estimation of Bilinear Systems Via Haar Wavelets, IEEE Trans. on Circuits and Systems 1: Fundamental Theory and Applications, 47 (2), 246–250.
  14. Hong, L., Girsang, I. P., Dhupia, J. S. (2016). Identification and control of stick–slip vibrations using Kalman estimator in oil-well drill strings, Journal of Petroleum Science and Engineering, 140, 119-127.
  15. Tan, C., Dai, W., Yeung, H., and Dong, F. (2015). A Kalman estimation based oil–water two-phase flow measurement with CRCC, International Journal of Multiphase Flow, 72, 306-317.
  16. Liu, M., Lai, J., Li, Z., and Liu, J. (2016). An adaptive cubature Kalman filter algorithm for inertial and land-based navigation system, Aerospace Science and Technology, 51, 52-60.
  17. Athans, M. (1971). The role and use of the stochastic linear quadratic Gaussian problem in control system design, IEEE Trans. Automatic Control AC-16, 529–552.
  18. Foster, M. R. and Saffman, P. G. (1970). The Drag of a Body Moving Transversely in a Confined Stratified Fluid, Journal of Fluid Mech., 43 (2), 407 – 418.
  19. Obinabo, E. C. and Ojieabu, C. E. (2010). Measurement Noise Filtrationand State Estimationofa Discrete-time Stochastic Process, International Journal of Soft Computing, 5 (2), .
  20. Bacher, R. A., Gray, W. and Murray–Smith, D. J. (1981). Time Domain System Identification Applied to Non – invasive Estimation Cardiopulmonary Quantities ibid, 128 (2), 56–64.
  21. Ho Y. C. (1962). On the Stochastic approximation method and optimal filtration theory. Journal Math anal Applications, 6, 152–154.
  22. Sage A. P. (1967). Least Squares Curve Fitting and Discrete Optimum Fitting, IEEE Trans. EDUC. E-10, 1, 29–38.
  23. Obinabo, E. C. (2015). the Development of Oil-rig Drill Ship positioning systems, Bowens Nigeria Limited, Warri, Nigeria, 425.
  24. Ljung I. (1979). Asymptotic behaviour of the extended Kalman filter as a parameter Estimation for linear Systems, IEEE trans, ac-24, 36-50.
  25. Kochhar, A. K. (1978). Comparison of Stochastic Identification Techniques for Dynamic Modelling of Plastic Extrusion Processes, Proc. Instn. Mech. Engrs., 192 (28), 299–309.

Article Tools
  Abstract
  PDF(592K)
Follow on us
ADDRESS
Science Publishing Group
548 FASHION AVENUE
NEW YORK, NY 10018
U.S.A.
Tel: (001)347-688-8931