Open AccessArticle

A Combined Gravity Compensation Method for INS Using the Simplified Gravity Model and Gravity Database

Xiao Zhou

^1,2,*,

Gongliu Yang

^1,2,

Jing Wang

³ and

Zeyang Wen

^1,2

School of Instrument Science and Opto-Electronics Engineering, Beihang University, Beijing 100191, China

Science and Technology on Inertial Laboratory, Beihang University, Beijing 100191, China

School of Mechatronices and Information Engineering, China Mining University, Beijing 100083, China

Author to whom correspondence should be addressed.

Sensors 2018, 18(5), 1552; https://doi.org/10.3390/s18051552

Submission received: 21 March 2018 / Revised: 29 April 2018 / Accepted: 8 May 2018 / Published: 14 May 2018

(This article belongs to the Section Physical Sensors)

Download

Browse Figures

Figure 1
Description of gravity disturbance vector (GDV). "> Figure 2
The deflection of the vertical (DOV). "> Figure 3
The north position error caused by gravity vertical deflection <math display="inline"> <semantics> <mi>ζ</mi> </semantics> </math>. "> Figure 4
The common logarithm of <math display="inline"> <semantics> <mrow> <mi>cov</mi> <msub> <mrow> <mo stretchy="false">(</mo> <msub> <mi>C</mi> <mrow> <mi>n</mi> <mi>m</mi> </mrow> </msub> <mo>,</mo> <mi>V</mi> <mo stretchy="false">)</mo> </mrow> <mi>P</mi> </msub> </mrow> </semantics> </math> when <math display="inline"> <semantics> <mrow> <msub> <mi>N</mi> <mrow> <mi>max</mi> </mrow> </msub> <mo>=</mo> <mn>70</mn> </mrow> </semantics> </math>. "> Figure 5
The common logarithm of <math display="inline"> <semantics> <mrow> <mrow> <mo>|</mo> <mrow> <mi>cov</mi> <msub> <mrow> <mo stretchy="false">(</mo> <msub> <mi>S</mi> <mrow> <mi>n</mi> <mi>m</mi> </mrow> </msub> <mo>,</mo> <mi>V</mi> <mo stretchy="false">)</mo> </mrow> <mi>P</mi> </msub> </mrow> <mo>|</mo> </mrow> </mrow> </semantics> </math> when <math display="inline"> <semantics> <mrow> <msub> <mi>N</mi> <mrow> <mi>max</mi> </mrow> </msub> <mo>=</mo> <mn>70</mn> </mrow> </semantics> </math>. "> Figure 6
The flow chart of the combined gravity compensation method in the inertial navigation system (INS). "> Figure 7
Field test device. "> Figure 8
Test profiles on Google map. (a) Test 1; (b) Test 2. "> Figure 9
The gravity anomaly and DOVs of two tests. (a) South-north DOV in test 1; (b) East-west DOV in test 1; (c) Gravity anomaly in test 1; (d) South-north DOV in test 2; (e) East-west DOV in test 2; (f) Gravity anomaly in test 2. "> Figure 10
Position errors of two tests. (a) North position error of test 1; (b) East position error of test 1; (c) Position error of test 1; (d) North position error of test 2; (e) East position error of test 2; (f) Position error of test 2. ">

Versions Notes

Abstract

In recent decades, gravity compensation has become an important way to reduce the position error of an inertial navigation system (INS), especially for a high-precision INS, because of the extensive application of high precision inertial sensors (accelerometers and gyros). This paper first deducts the INS’s solution error considering gravity disturbance and simulates the results. Meanwhile, this paper proposes a combined gravity compensation method using a simplified gravity model and gravity database. This new combined method consists of two steps all together. Step 1 subtracts the normal gravity using a simplified gravity model. Step 2 first obtains the gravity disturbance on the trajectory of the carrier with the help of ELM training based on the measured gravity data (provided by Institute of Geodesy and Geophysics; Chinese Academy of sciences), and then compensates it into the error equations of the INS, considering the gravity disturbance, to further improve the navigation accuracy. The effectiveness and feasibility of this new gravity compensation method for the INS are verified through vehicle tests in two different regions; one is in flat terrain with mild gravity variation and the other is in complex terrain with fierce gravity variation. During 2 h vehicle tests, the positioning accuracy of two tests can improve by 20% and 38% respectively, after the gravity is compensated by the proposed method.

Keywords:

error modelling; gravity model; extreme learning machine (ELM); gravity compensation; high precision free-INS

1. Introduction

Achieving accurate navigation is crucial for the modern motion carrier. Nowadays the most common way to navigate for most modern vehicles and aircrafts is utilizing the global positioning system (GPS). The significant advantage of GPS is its real-time high precision positioning with the help of satellites’ signals. However, GPS cannot receive the satellites’ signals in some circumstances due to the physical blockages, such as in a tunnel or underwater. The inertial navigation system (INS) utilizes Newton’s law to realize autonomous navigation worldwide, and the initialization errors propagate throughout the trajectory with the increase in time. Although the long-term navigation accuracy of free-INS cannot compare with that of GPS, it is still essential during times of loss of GPS. Because INS can work in all weather conditions, has high anti-jamming capacity, and radiates no signal to the outside, INS has been used for many military and civil applications [1]. In recent years, with the emergence of cold atom interferometry, developing high precision INS has become possible. It is predicted that the accelerometer and gyro of this INS can have an accuracy of 10⁻⁸ m/s² and 10⁻⁶ deg/h, respectively, which is at least 100 times better than the current INS. The significant improvement of INS, especially on inertial sensors, leaves the gravity compensation as the most important part effecting the positioning accuracy of INS, particularly for rough topological areas [2].

Normally, there are mainly three viable means to realize gravity compensation for INS [3]. Firstly there is the conventional method that uses gravitational gradiometers to gain the gravity disturbance on the trajectory [4,5]. This case is hard to apply to engineering applications because of the high cost of devices. The second method is using gravity models (like WGS84, DQM2000, EGM2008, etc.) to calculate the gravity disturbance and compensate it into the error equations of INS [6,7]. This method is easy to operate but the complexity of these models leads to a long time needed for a solution, which will hinder the real-time navigation of INS. The third method is to obtain the gravity disturbance with the help of the interpolation method based on a gravity database and compensate it into error equations of INS incorporated with gravity disturbance, which has a relatively ideal compensation result when gravity data is available [8,9,10,11,12]. Considering the characteristics of the second and the third method, in this paper, a combined gravity compensation method for INS using the simplified gravity model and gravity data is proposed to restrain the error propagation of INS. At last, two vehicle tests were applied to prove the effectiveness and precision of the proposed gravity compensation method.

The rest of this paper is organized as follows: Section 2 states the error analysis of INS solution, considering gravity disturbance. Section 3 presents the principle of the simplified gravity model. Section 4 introduces the theory and frame of the gravity data-based compensation method using ELM. Section 5 exhibits the framework of this combined gravity compensation method in INS. Section 6 performs vehicle tests in two different regions. Section 7 concludes the whole paper.

2. Error Analysis of INS Solution Caused by Gravity Disturbance

2.1. Definition of Gravity Disturbance Vector

As Figure 1 shows, the gravity disturbance vector (GDV) is the vector difference between the normal gravity and the actual gravity on the same point in space. GDV is composed of two parts: the orthogonal component is denoted as gravity anomaly and the tangential component is denoted as vertical deflection [13].

In Figure 1, g_p and γ_p are defined as the gravity vector and the normal gravity vector at the point P, respectively. The difference between g_p and γ_p is defined as gravity disturbance vector δg:

δ g = g_{P} - γ_{P}

(1)

As Figure 2 shows, the difference in direction is denoted as the deflections of the vertical (DOVs) [14,15]. The DOVs are composed of two components, ζ and η are represented as the north–south component and the east–west component, respectively.

δ g = [Δ g_{E} Δ g_{N} Δ g_{U}]

(2)

{\begin{cases} δ g_{N} = - γ_{0} ζ \\ δ g_{E} = - γ_{0} η \end{cases}

(3)

where, Δg_N, Δg_E, and Δg_U are presented as the north, east, and vertical components of gravity disturbance, respectively. γ₀ is the value of normal gravity.

2.2. INS Error Equations Incorporated with Gravity Disturbance

The following INS error equations are considered with gravity disturbance [16,17,18,19]:

\begin{array}{l} δ {\dot{V}}^{n} = - ϕ^{n} \times f^{n} + C_{b}^{n} (δ K_{A} + δ A) f^{b} + δ V^{n} \times (2 ω_{i e}^{n} + ω_{e n}^{n}) \\ + V \times (2 δ ω_{i e}^{n} + δ ω_{e n}^{n}) + \nabla^{n} + δ g^{n} \end{array}

(4)

{\begin{cases} δ \dot{L} = \frac{δ V_{N}}{R_{M} + h} \\ δ \dot{λ} = \frac{δ V_{E}}{R_{N} + h} \sec L + δ L \frac{V_{E}}{R_{N} + h} \tan L \sec L \\ δ \dot{h} = δ V_{U} \end{cases}

(5)

\dot{ϕ} = ϕ \times ω_{i n}^{n} + δ ω_{i n}^{n} - C_{b}^{n} (δ K_{G} + δ G) ω_{i b}^{b} + ε^{n}

(6)

in the above formulas,

δ V^{n}

is defined as the velocity error in navigation frame,

ϕ^{n}

is defined as attitude error of the motion carrier,

f^{n}

is defined as the specific force in navigation frame,

C_{b}^{n}

is a direction cosine matrix which can be used to transform the body acceleration vector into the navigation frame,

δ K_{A}

and

δ A

are defined as scale coefficient error and installation angle error of accelerometer respectively,

δ G

and

δ K_{G}

are defined as installation angle error of gyros and scale coefficient error respectively,

\nabla^{n}

is defined as the bias of accelerometers in navigation frame,

δ g^{n}

is defined as the gravity disturbance in navigation frame and

ε^{n}

is the drift of gyros in navigation frame,

δ ω_{i n}^{n}

can be calculated from the following equation:

δ ω_{i n}^{n} = δ ω_{i e}^{n} + δ ω_{e n}^{n}

(7)

where

ω_{i e}^{n}

and

ω_{e n}^{n}

are presented as earth’s rotation rate and the navigation frame’s rotation with respect to earth respectively; both are expressed in the navigation frame. They can be calculated from the following formula:

δ ω_{e n}^{n} = [\begin{matrix} - \frac{δ V_{N}}{R_{M} + h} \\ \frac{δ V_{E}}{R_{N} + h} \\ \frac{δ V_{E} \tan L}{R_{N} + h} + δ L \frac{V_{E} \sec^{2} L}{R_{N} + h} \end{matrix}]

(8)

δ ω_{i e}^{n} = {[0 - δ L ω_{i e} \sin L δ L ω_{i e} \cos L]}^{T}

(9)

where

L

λ

, and

h

are denoted as the latitude, longitude, and altitude of the motion carrier respectively,

R_{M}

R_{N}

are denoted as meridian radius and prime vertical radius, respectively.

From Equations (4)–(9), we could draw the conclusion that INS velocity errors

δ V^{n}

are mainly caused by the accelerometer’s output error

(δ K_{A} + δ A) f^{b} + \nabla^{n}

, velocity error

δ V^{n}

, and gravity disturbance

δ g^{n}

. And the velocity error

δ V^{n}

is the main error source that causes the position errors of INS.

Based on the above analysis, we can see that gravity disturbance

δ g^{n}

first leads to INS velocity errors in Equation (4), then further affects the position and attitude accuracy of INS through Equations (5) and (6). With the improvement of the accuracy of inertial sensors, INS solution errors due to gravity disturbance are becoming more and more significant and thus cannot be neglected; some measures must be taken to compensate for it.

Here we calculate the north position error caused by three typical values of north–south gravity vertical deflection

ζ

as an example to better represent the influence of gravity disturbance on INS. And the simulation results are shown in Table 1 and Figure 3.

Table 1 and Figure 3 show that the north position error of INS due to gravity vertical deflection propagates according to the Schuler cycle (about 84.4 min), and the amplitude of the position error is in direct proportion to the value of gravity vertical deflection. Thus, the influence of gravity disturbance in a high-precision INS cannot be neglected; it must be compensated to realize a more accurate navigation result.

3. The Principle of the Simplified Gravity Model

For a simplified reference, the Earth can be assumed to take the shape of an ellipsoid of revolution with a uniform mass distribution. There are some Earth models that can describe the Earth’s characteristics, such as WGS84, DQM2000, EGM2008, etc., which are normally used to describe the Earth’s gravitational field and are expressed with high degree spherical harmonic polynomials. Normally, in most inertial navigators, the WGS84 reference ellipsoid is used to calculate the normal gravity and the order of it only expands to 2 degrees. However, with the development of the gravity models, WGS84 cannot perfectly describe the gravitational field to meet the demands of high accurate inertial navigators. To further restrain the errors caused by the gravity model, the more accurate gravity model (EGM2008) is used to replace the WGS84 model to calculate the normal gravity. Although the EGM2008 model can well describe the characteristics of the Earth, it is hard to apply this gravity model with a high number of degrees to the solution of the INS because of the time and space complexity [7]. Thus, here a simplified gravity model is proposed to solve this problem.

3.1. Expression of Spherical Harmonic Gravity Model (SHM)

The gravitational potential of the earth to external point is a harmonic function, it can be represented by the infinite order number of a spherical harmonic function [20,21]:

V_{(ρ, θ, λ)} = \frac{f M}{ρ} [1 - \sum_{n = 2}^{\infty} {(\frac{a_{e}}{ρ})}^{n} \sum_{m = 0}^{n} ({\bar{C}}_{n m} \cos (m λ) + {\bar{S}}_{n m} \sin (m λ) {\bar{P}}_{n m} (\cos θ))],

(10)

where

(ρ, θ, λ)

is the spherical coordinate with the center of the earth as the origin of the coordinate;

f

is denoted as the Newton’s gravitational constant;

a_{e}

is the major radius of reference ellipsoid;

M

is the mass of the Earth;

{\bar{P}}_{n m} (\cos θ)

is associated Legendre polynomials,

{\bar{C}}_{n m}, {\bar{S}}_{n m}

are denoted as constant coefficients of the spherical harmonic of degree n and order m. Therefore, if these coefficient values are known, the value of the gravitational potentiation can be obtained. In order to get the value of the gravitational potentiation, the number of coefficients is limited. Thus, Equation (10) can be rewritten as follows:

V_{(ρ, θ, λ)} = \frac{f M}{ρ} [1 - \sum_{n = 2}^{N} {(\frac{a_{e}}{ρ})}^{n} \sum_{m = 0}^{n} ({\bar{C}}_{n m} \cos (m λ) + {\bar{S}}_{n m} \sin (m λ) {\bar{P}}_{n m} (\cos θ))],

(11)

where N is the highest order of the known coefficients.

3.2. The Selection of the Degree for the Simplified Gravity Model

The covariance between the gravity potential and harmonic coefficients considers the degree of their correlation. The effect of certain coefficients on the gravity potential is reflected by a high value of the covariance. Therefore, we could theoretically use it to determine a suitable degree of SHM. In the gravity potential theory, Equation (11) can be changed into the following equation by multiplying both sides with a specific

P_{n m} (\cos θ) \cos (m λ)

P_{n m} (\cos θ) \sin (m λ)

and integrates over the unit sphere [22,23]:

\begin{matrix} C_{n m} \\ S_{n m} \end{matrix}} = \frac{ρ}{4 π f M} {(\frac{ρ}{a})}^{n} \iint_{σ} V P_{n m} (\cos θ) {\begin{matrix} \cos (m λ) \\ \sin (m λ) \end{matrix}} d σ,

(12)

where

d σ = \sin θ d θ d λ

is defined as the surface element of the unit sphere.

Defining

L^{C}

and

L^{S}

as a linear functional, so the Equation (12) can be simplified as follows:

\begin{matrix} C_{n m} \\ S_{n m} \end{matrix}} = {\begin{matrix} L_{n m}^{C} \\ L_{n m}^{S} \end{matrix}} V .

(13)

The linear functional is defined as the continuous analogue to the usual concept of a linear function in n-dimensional vector space. Therefore, the covariance between

C_{n m}, S_{n m}

and T satisfy the following condition:

cov ({\begin{matrix} C_{n m} \\ S_{n m} \end{matrix}}, V) = {\begin{matrix} L_{n m}^{C} \\ L_{n m}^{S} \end{matrix}} K (P),

(14)

where

K (P)

is the covariance of the gravity potential on point P, described as

\begin{array}{l} K (P) = \sum_{n = 2}^{\infty} \frac{K_{n}}{2 n + 1} {(\frac{R}{ρ})}^{2 (n + 1)} \\ \times \sum_{m = 0}^{n} P_{n m} {(\cos θ)}^{2} \cos (m λ), \end{array}

(15)

where R is the Earth’s average radius,

ρ

is the geocentric radius distance of P, and

K_{n}

is the degree variance of coefficients:

K_{n} = {(\frac{f M}{R})}^{2} \sum_{m = 2}^{n} (C_{n m}^{2} + S_{n m}^{2}) .

(16)

Consequently, the covariance between the gravity potential and the harmonic coefficients can be described as

cov ({\begin{matrix} C_{n m} \\ S_{n m} \end{matrix}}, V) = \frac{K_{n}}{2 n + 1} \frac{R}{f M} {(\frac{R}{ρ})}^{n + 1} P_{n m} (\cos θ) {\begin{matrix} \cos (m λ) \\ \sin (m λ) \end{matrix}} .

(17)

The maximum number of degrees of EGM2008 is 2160 degrees, however, the high order terms have only little influence on the results of calculation. Thus, to better illustrate the number of orders of EGM2008 affecting on the results of calculation, here we set the maximum order (

N_{\max} = 70

) to calculate the absolute magnitudes of

cov {(C_{n m}, V)}_{P}

and

cov {(S_{n m}, V)}_{P}

. Since the covariance of the points around the world are similar, the point

P (θ = 45^{\circ}, λ {= 45}^{\circ})

on the Earth’s surface is used for discussion. When

N_{\max} = 70

, the absolute magnitudes of

cov {(C_{n m}, V)}_{P}

and

cov {(S_{n m}, V)}_{P}

range from

10^{- 3}

10^{- 15}

and

10^{- 3}

10^{- 13}

, respectively. To better visualize the results, Figure 4 and Figure 5 are plotted with common logarithms of

| cov {(C_{n m}, V)}_{P} |

and

| cov {(S_{n m}, V)}_{P} |

respectively, shown as follows:

In Figure 4 and Figure 5, horizontal ordinate (m) and longitudinal coordinate (n) represent the number of degrees of the EGM2008. In these two figures, the darker place means the greater value of logarithms of

| cov {(C_{n m}, V)}_{P} |

and

| cov {(S_{n m}, V)}_{P} |

. From these two figures, we could see that the values of the covariance are mainly focused below 12 degrees, which indicates the value of gravity mainly focuses on the SHM’s lower degree. Thus, in this paper, we set the maximum order of the simplified SHM as

N_{\max} = 12

to calculate the normal gravity.

4. The Principle and Procedure of the Data-Based Gravity Disturbance Compensation Method for INS Using ELM

4.1. The Principle of ELM

In recent years, ELM has been becoming a popular research method to solve many engineering problems. ELM has its own distinguishing features, which are depicted as follows [24]:

(1): Only the predefined network structure needs to be modulated;
(2): ELM has the ability to do fast learning;
(3): High generalization performance can be achieved through ELM;
(4): A wide selection range of activation functions can be used in ELM.

Due to its many advantages, in this paper, an estimation method based on ELM is used to estimate the gravity disturbance on the trajectory of the motion carrier, then compensate it into the INS solution equations to improve the navigation accuracy.

The principle of the ELM is described as follows [25]:

Given that we have a training set

{(x_{i}, t_{i})}_{i = 1}^{N}

with N distinct examples, where

x_{i} = {[x_{i 1}, x_{i 2}, x_{i 3}, \dots, x_{i n}]}^{T}

has n inputs and

t_{i} = {[t_{i 1}, t_{i 2}, t_{i 3}, \dots, t_{i m}]}^{T}

has m outputs. Here we define each

x_{i}

has longitude and latitude two values, described as

x_{i} = [λ_{i} L_{i}]

and the output

t_{i}

is the gravity disturbance on geoid, described as

t_{i} = [Δ g_{E i} Δ g_{N i} Δ g_{U i}]

. Here

l

is denoted as the number of hidden neurons, ω is denoted as the

l \times n

input weight matrix where

ω_{j} = {[ω_{j 1}, ω_{j 2}, ω_{j 3}, \dots, ω_{j n}]}^{T}

, b is the

l \times 1

biases vector and

β

is the

l \times m

output weight matrix where

β_{j} = {[β_{j 1}, β_{j 2}, β_{j 3}, \dots, β_{j m}]}^{T}

. Generally, the ELM network function is shown as

t_{i} = \sum_{j = 1}^{l} β_{j} g (ω_{j} \cdot x_{i} + b_{j}), i = 1, 2, \dots, N;

(18)

where

j \in {1, 2, \dots, l}, ω_{j} \cdot x_{i}

is defined as the inner product of

ω_{j}

and

x_{i}

, the activation function g(x) chooses sigmoid function in the proposed method, which is formulated as

g (x) = \frac{1}{1 + e x p [- (ω \cdot x + b)]}

(19)

Equation (18) can be simplified as the following simple form:

H β = T

(20)

where

H = {[\begin{matrix} g (ω_{1} \cdot x_{1} + b_{1}) & \dots & g (ω_{l} \cdot x_{1} + b_{l}) \\ ⋮ & ⋱ & ⋮ \\ g (ω_{1} \cdot x_{N} + b_{1}) & \dots & g (ω_{l} \cdot x_{N} + b_{l}) \end{matrix}]}_{N \times l}

(21)

H is denoted as the hidden layer output matrix of the network. T is the output matrix and

T = {[t_{1}, t_{2}, \dots, t_{N}]}^{T}

According to ELM theory, the input weights and hidden biases can be selected at random, and the output weights can be obtained by calculating the following formula:

β = H^{+} T

(22)

where H⁺ is the Moore-Penrose (MP) generalized inverse of matrix H.

4.2. The Procedure of the Data-Based Gravity Disturbance Compensation Method in INS Using ELM

The process of the data-based gravity disturbance compensation method for INS using ELM is described as follows:

Obtain the motion carrier’s position. Get the position information (λ, L) of the motion carrier through the INS calculation.
Choose the gravity data base. Find the suitable gravity data base taking the position obtained by step 1 as the center. Normally, 5′ × 5′ gravity grid data base is chosen as the training database.
ELM training. Set the motion carrier’s position information (λ, L) as the inputs of the network and acquire the gravity disturbance on geoid ( $Δ g_{E 0}, Δ g_{N 0}, Δ g_{U 0}$ ) through training the gravity data base acquired by step 2.
Upward continuation. Calculate the gravity disturbance with upward continuation to the point where INS is. The height of INS is acquired through altimeter. In the geographic engineering application, the most practical upward continuation method is free air correction. The computational equation is described as follows [26]:

$Δ g = Δ g_{0} - 0.3086 H$

(23)

where $Δ g$ is the gravity disturbance where the motion carrier is and $Δ g_{0}$ is the gravity disturbance on the geoid. The units of $Δ g$ and $Δ g_{0}$ are milligal. H is the height between the geoid and the motion carrier. The unit of H is meter.

5. The Framework of the Combined Gravity Compensation Method for INS

The flow chart of the combined gravity compensation method for INS is illustrated in Figure 6:

The process of the combined gravity compensation method is described as follows:

Obtain the motion carrier’s position. Get the position information of the motion carrier through the INS calculation.
Calculate the gravity and the gravity disturbance. Based on the position information from step 1, the gravity g and the gravity disturbance $δ g$ are obtained respectively through the simplified gravity model and ELM training.
Compensate the gravity and the gravity disturbance from step 2 into INS equations, considering the gravity disturbance to restrain the error propagation.

6. Experiment

To prove the effectiveness and feasibility of the combined gravity compensation method in INS, two vehicle tests were applied in the city of Xi’an in China and denoted as test 1 and test 2. Test 1 and test 2 were applied in the areas with flat terrain and complex terrain respectively. On the test vehicle, an altimeter, a GPS receiver, a set of batteries, a high-precision INS, and a data processing computer (PC) were carried. The high-precision INS and the batteries were fixed on the overplate at the trunk. The computer was used to collect the data of GPS and INS in real time, then the data was resolved for navigation. The test vehicle is shown in Figure 7. The accuracy of the inertial sensors and GPS are listed in Table 2.

The travel profiles, gravity anomalies and DOVs of the two tests are illustrated in Figure 8 and Figure 9. In Figure 9, the gravity anomalies and DOVs of the two tests were obtained from the estimation based on the measured gravity data using ELM training. The gravity data is provided by the Institute of Geodesy and Geophysics; Chinese Academy of sciences.

From Figure 9, we can see that the values of DOVs and gravity anomaly in the area of test 1 are both larger than that of test 2. To better prove the effectiveness and accuracy of the combined gravity compensation method, the maximum position errors were used to prove this. Here we regarded the position result of GPS as the true value, and respectively compared it with the position results compensated with three different gravity compensation methods: 1 only with the reference ellipse (shown as blue colour); 2 with the reference ellipse and gravity disturbance (shown as black colour); 3 with the proposed compensation method (shown as red colour).The position error results are shown in Figure 10, in which (a), (b), (c) are the position results of test 1; (d), (e), (f) are the position results of test 2.

According to Figure 10, the maximum position errors of INS compensated with three different methods in two tests are listed in Table 3.

From Figure 10 and Table 3, we could draw the conclusion that the proposed gravity compensation method performed well in two tests, especially in test 2 with a wider variation range of gravity disturbance. During the 2 h vehicle tests, compared with the solution results compensated with the reference ellipse only, the maximum value of position errors can reduce by 20% and 38% respectively, in test 1 and test 2, with the proposed gravity compensation method. Thus, the effectiveness of this combined gravity compensation method for INS is proved.

7. Conclusions

This paper proposed a combined gravity compensation method for high-precision INS. The method utilizes the simplified gravity model to subtract the normal gravity from the navigation equations, and meanwhile compensates the gravity disturbance through the ELM training based on measured gravity data. The effectiveness and accuracy of the combined gravity compensation method were verified by vehicle tests, which were applied with a high-precision INS, and the solution results show that the maximum position errors can reduce by 20% and 38%, respectively, in flat terrain and complex terrain. It should be noted that in this work the tests were carried on an experimental vehicle running on the surface of the Earth with relatively low speed. Therefore, in the follow-up study, a flight test with a relatively high travel speed should applied to further test the applicability of this gravity compensation method.

Author Contributions

In this paper, X.Z. proposed the combined gravity compensation method for high precision INS and wrote this paper. G.Y. analyzed the feasibility of the proposed method and modified the English. J.W. provided the source code of gravity field model. Z.W. analyzed the data.

Funding

The project is supported by the National Natural Science Foundation of China (61340044) and the Fundamental Research Funds for the Central Universities (YWF-10-01-B30).

Acknowledgments

The writers would like to thank the Institute of Geodesy and Geophysics, Chinese Academy of Sciences, for providing the gravity data. Many thanks are also given to Hua Chai at the Institute of Geodesy and Geophysics, Chinese Academy of Sciences, for his valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

Titterton, D.H.; Weston, J.L. Strapdown Inertial Navigation Technology, 2nd ed.; Institution of Engineering and Technology: Stevenage, UK, 2004; pp. 25–30. ISBN 9780863413582. [Google Scholar]
Zhou, X.; Yang, G.L.; Cai, Q.Z.; Wang, J. A novel gravity compensation method for high precision free-INS based on Extreme Learning Machine. Sensors 2016, 16, 2019. [Google Scholar] [CrossRef] [PubMed]
Kwon, J.H. Gravity Compensation Methods for Precision INS. In Proceedings of the 60th Annual Meeting of the Institute of Navigation, Dayton, OH, USA, 7–9 June 2014; pp. 483–490. [Google Scholar]
Jekeli, C. Precision free-inertial navigation with gravity compensation by an onboard gradiometer. J. Guid. Control Dyn. 2006, 29, 704–713. [Google Scholar] [CrossRef]
Mandour, I.M.; EI-Dakiky, M.M. Inertial Navigation System Synthesis Approach and Gravity-Induced Error Sensitivity. IEEE Trans. Aerosp. Electron. Syst. 1988, 24, 40–50. [Google Scholar] [CrossRef]
Siouris, G.M. Gravity modeling in aerospace applications. Aerosp. Sci. Technol. 2009, 13, 301–315. [Google Scholar] [CrossRef]
Wang, J.; Yang, G.; Li, J.; Zhou, X. An Online Gravity Modeling Method Applied for High Precision Free-INS. Sensors 2016, 16, 1541. [Google Scholar] [CrossRef] [PubMed]
Arora, N.; Russell, R.P. Fast efficient and adaptive interpolation of the geopotential. In Proceedings of the AAS/AIAA Astrodynamics Specialist Conference, Girdwood, AK, USA, 31 July–4 August 2011. [Google Scholar]
Sarzeaud, O.M.; Rouxel, D. Optimal Interpolation of Gravity Maps Using a Modified Neural Network. Int. Assoc. Math. Geosci. 2009, 41, 379–395. [Google Scholar] [CrossRef]
Zhao, Y.F.; Ren, X.H. CNC Thermal Compensation Based on Mind Evolutionary Algorithm Optimized BP Neural Network. World J. Eng. Technol. 2016, 4, 38–44. [Google Scholar] [CrossRef]
Chen, Y.S.; Yan, S. Mind-Evolution-Based Machine Learning Framework and the Implementation of Optimization. In Proceedings of the IEEE International Conference on Intelligent Engineering Systems, Vienna, Austria, 17–19 September 1998; pp. 355–359. [Google Scholar]
Hang, L. Methods of Statistical Learning; Tsinghua University Press: Beijing, China, 2012; ISBN 978-7-302-27595-4. [Google Scholar]
Wellenhof, B.H.; Moritz, H. Physical Geodesy, 2nd ed.; Springer: Vienna, Austria, 2006; ISBN 9783211335444. [Google Scholar]
Lu, Z.L. The Theory and Method of the Earth Gravity Field; Liberation Army Press: Beijing, China, 1996; ISBN 7-5065-3033-3. [Google Scholar]
Hofmann-Wellenhof, B.; Moritz, H. Physical Geodesy; Springer: Berlin, Germany, 2006; ISBN 3211335447. [Google Scholar]
Jordan, S.K. Effects of geodetic uncertainties on a damped inertial navigation system. IEEE Trans. Aerosp. Electron. Syst. 1973, 741–752. [Google Scholar] [CrossRef]
Vanderwerf, K. Schuler pumping of inertial velocity errors due to gravity anomalies along a popular North Pacific airway. In Proceedings of the IEEE Position Location and Navigation Symposium, Atlanta, GA, USA, 22–25 April 1996; pp. 642–648. [Google Scholar]
Levine, S.A.; Gelb, A. Effect of reflections of the vertical on the performance of a terrestrial inertial navigation system. J. Spacecr. Rocket. 1969, 6, 978–984. [Google Scholar] [CrossRef]
Leonard, J.M.; Nievinski, F.G.; Born, G.H. Gravity error compensation using second-order Gauss-Markov processes. J. Spacecr. Rocket. 2013, 50, 217–229. [Google Scholar] [CrossRef]
Hirt, C.; Marti, U.; Bürki, B. Assessment of EGM2008 in Europe using accurate astrogeodetic vertical deflections and omission error estimates from SRTM/DTM2006.0 residual terrain model data. J. Geophys. Res. 2010, 115. [Google Scholar] [CrossRef]
Fu, G.Y.; Zhu, Y.Q.; Gao, S.H.; Liang, W.; Jin, H.; Yang, G.; Zhou, X.; Guo, S.; Xu, Y.; Du, W. Discrepancies between free air gravity anomalies from EGM2008 and the ones from dense gravity/GPS observations at west Sichuan Basin. Chin. J. Geophys. 2013, 56, 3761–3769. [Google Scholar] [CrossRef]
Abeyratne, P.G.V.; Featherstone, W.E.; Tantrigoda, D.A. Assessment of EGM2008 over Sri Lanka, an area where “fill-in” data were used in EGM2008. Newtons Bull. 2009, 4, 284–316. [Google Scholar]
Dawod, G.; Mohamed, H.; Ismail, S. Evaluation and adaptation of the EGM2008 geopotential model along the Northern Nile Valley, Egypt: Case study. J. Surv. Eng. 2009, 136, 36–40. [Google Scholar] [CrossRef]
Cho, D.; Ham, J.; Oh, J.; Park, J.; Kim, S.; Lee, N.K.; Lee, B. Detection of Stress Levels from Biosignals Measured in Virtual Reality Environments Using a Kernel-Based Extreme Learning Machine. Sensors 2017, 17, 2435. [Google Scholar] [CrossRef] [PubMed]
Khosravi, V.; Ardejani, F.D.; Yousefi, S. Monitoring soil lead and zinc contents via combination of spectroscopy with extreme learning machine and other data mining methods. Geoderma 2018, 318, 29–41. [Google Scholar] [CrossRef]
Jekeli, C.; Lee, J.K.; Kwon, J.H. Modeling errors in upward continuation for INS gravity compensation. J. Geod. 2007, 81, 297–309. [Google Scholar] [CrossRef]

Figure 1. Description of gravity disturbance vector (GDV).

Figure 2. The deflection of the vertical (DOV).

Figure 3. The north position error caused by gravity vertical deflection

ζ

Figure 3. The north position error caused by gravity vertical deflection

ζ

Figure 4. The common logarithm of

cov {(C_{n m}, V)}_{P}

when

N_{\max} = 70

Figure 4. The common logarithm of

cov {(C_{n m}, V)}_{P}

when

N_{\max} = 70

Figure 5. The common logarithm of

| cov {(S_{n m}, V)}_{P} |

when

N_{\max} = 70

Figure 5. The common logarithm of

| cov {(S_{n m}, V)}_{P} |

when

N_{\max} = 70

Figure 6. The flow chart of the combined gravity compensation method in the inertial navigation system (INS).

Figure 7. Field test device.

Figure 8. Test profiles on Google map. (a) Test 1; (b) Test 2.

Figure 9. The gravity anomaly and DOVs of two tests. (a) South-north DOV in test 1; (b) East-west DOV in test 1; (c) Gravity anomaly in test 1; (d) South-north DOV in test 2; (e) East-west DOV in test 2; (f) Gravity anomaly in test 2.

Figure 10. Position errors of two tests. (a) North position error of test 1; (b) East position error of test 1; (c) Position error of test 1; (d) North position error of test 2; (e) East position error of test 2; (f) Position error of test 2.

Table 1. The simulation results of north position error caused by typical gravity vertical deflection.

	$ζ = 1 \sec$	$ζ = 3 \sec$	$ζ = 5 \sec$
$Δ g_{N}$ (m Gal *)	5	14	24
North position error (m)	62	185	309

* 1 m Gal = 1 × 10⁻⁵ m/s².

Table 2. Accuracy of the sensors for the test.

Sensors Types	Characteristics	Magnitude (1 σ)
Gyroscope	Constant Bias	0.003°/h
Accelerometer	Constant Bias	10 μg
GPS velocity	Horizontal error	0.03 m/s
GPS velocity	Height error	0.05 m/s
GPS position	Horizontal error	2 m
GPS position	Height error	5 m

Table 3. The maximum value of position error compared with GPS result (Unit: meter).

	The Reference Ellipse Only	The Reference Ellipse with DOVs	With the Proposed Gravity Compensation	Position Improvement (Compared with the Reference Ellipse Only)
Test 1	1050	1012	837	213(20%)
Test 2	1120	876	689	431(38%)

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Zhou, X.; Yang, G.; Wang, J.; Wen, Z. A Combined Gravity Compensation Method for INS Using the Simplified Gravity Model and Gravity Database. Sensors 2018, 18, 1552. https://doi.org/10.3390/s18051552

AMA Style

Zhou X, Yang G, Wang J, Wen Z. A Combined Gravity Compensation Method for INS Using the Simplified Gravity Model and Gravity Database. Sensors. 2018; 18(5):1552. https://doi.org/10.3390/s18051552

Chicago/Turabian Style

Zhou, Xiao, Gongliu Yang, Jing Wang, and Zeyang Wen. 2018. "A Combined Gravity Compensation Method for INS Using the Simplified Gravity Model and Gravity Database" Sensors 18, no. 5: 1552. https://doi.org/10.3390/s18051552

APA Style

Zhou, X., Yang, G., Wang, J., & Wen, Z. (2018). A Combined Gravity Compensation Method for INS Using the Simplified Gravity Model and Gravity Database. Sensors, 18(5), 1552. https://doi.org/10.3390/s18051552

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

A Combined Gravity Compensation Method for INS Using the Simplified Gravity Model and Gravity Database

Abstract

1. Introduction

2. Error Analysis of INS Solution Caused by Gravity Disturbance

2.1. Definition of Gravity Disturbance Vector

2.2. INS Error Equations Incorporated with Gravity Disturbance

3. The Principle of the Simplified Gravity Model

3.1. Expression of Spherical Harmonic Gravity Model (SHM)

3.2. The Selection of the Degree for the Simplified Gravity Model

4. The Principle and Procedure of the Data-Based Gravity Disturbance Compensation Method for INS Using ELM

4.1. The Principle of ELM

4.2. The Procedure of the Data-Based Gravity Disturbance Compensation Method in INS Using ELM

5. The Framework of the Combined Gravity Compensation Method for INS

6. Experiment

7. Conclusions

Author Contributions

Funding

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI