Next Article in Journal
Improved Taillight Detection Model for Intelligent Vehicle Lane-Change Decision-Making Based on YOLOv8
Next Article in Special Issue
Proposal of a Cost-Effective and Adaptive Customized Driver Inattention Detection Model Using Time Series Analysis and Computer Vision
Previous Article in Journal
Simulating Noise, Vibration, and Harshness Advances in Electric Vehicle Powertrains: Strategies and Challenges
Previous Article in Special Issue
Detecting Wheel Slip to Suppress Self-Excited Oscillations in Braking Mode
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Method for Land Vehicle Gravity Anomaly Measurement Combining an Inertial Navigation System, Odometer, and Geo-Information System

Laboratory of Intelligent Control, Rocket Force University of Engineering, Xi’an 710025, China
*
Author to whom correspondence should be addressed.
World Electr. Veh. J. 2024, 15(8), 368; https://doi.org/10.3390/wevj15080368
Submission received: 20 June 2024 / Revised: 5 August 2024 / Accepted: 8 August 2024 / Published: 14 August 2024

Abstract

:
Land vehicle gravity anomaly measurement relies heavily on global navigation satellite system (GNSS). However, when gravity measurement is carried out in special environments such as forests, valleys and tunnels, GNSS observation quality will inevitably decline, which directly affects the accuracy of gravity anomaly measurement. From the point of view of the gravity anomaly measurement principle, obtaining accurate elevation information of the test line is the premise to ensure the accuracy of gravity anomaly measurement. Thus, this paper proposes a strapdown land vehicle dynamic gravity anomaly measurement method combining an odometer and a geo-information system. In this method, strapdown inertial navigation errors are suppressed by observing the velocity of the odometer output. Then, the position information obtained by the combined navigation is entered into the geo-information system to obtain the elevation. The results of a single test line show that the external coincidence accuracy of the proposed method is 1.65 mGal, and the accuracy is comparable to the traditional GNSS assisted land vehicle gravimetry method. In addition, compared with the odometer assisted land vehicle gravimetry method, the external coincidence accuracy is increased by 30%.

1. Introduction

The gravity field is a fundamental physical property of the Earth, reflecting the distribution and variation of materials within it [1]. According to the different physical quantities measured, the methods of measuring the Earth’s gravity field can be divided into absolute gravity measurement and relative gravity measurement. Absolute gravity measurement directly measures the value of gravity at a single point on the ground. Relative gravity measurement measures the gravity difference between two test points and indirectly obtains the gravity value of each test point. Modern absolute gravity measuring instruments are complex precision instruments based on laser technology or atomic interference technology, and the measurement accuracy can reach μGal level [2,3]. However, the adaptability of the absolute gravity measuring instrument to the field environment is not good, and it is relatively complex to realize the gravity measurement of the moving base [4,5,6]. Gravimeters made using the principle of pressure balance or force balance fall into the category of relative gravity measurement [7,8]. At present, the instruments used for relative gravity measurement all adopt the principle of force balance, and the core components are quartz springs or metal wires [9].
To improve the efficiency of gravity measurement, many scholars have studied the gravity measurement methods of the dynamic base carried out by relative gravimeters on underwater robots, ships, hot air balloons, drones, and other carriers [10,11,12,13]. According to different carriers, gravity measurement of the dynamic base can be divided into ship, underwater, airborne, and land vehicle gravity measurements. Among them, the primary use of airborne and land vehicle gravity surveys is the rapid measurement of the ground gravity field. Airborne gravity measurement is highly efficient and can be completed in complex areas such as deserts, glaciers, and forests. However, the airborne gravity measurement method measures the gravity information of the aircraft’s flight path, which often needs to be extended downward [14]. Compared with airborne gravimetry, land vehicle gravimetry can obtain gravity information with a high signal-to-noise ratio and good accuracy without complicated data extension. Land vehicle gravimetry has been widely used in geoid determination, groundwater resource surveying, underground cavity detection, gravity assisted navigation, and other fields [15,16].
For a long time, global satellite navigation systems (GNSS) have played a vital role in dynamic gravity measurement activities [17,18,19,20,21,22,23]. First, observing GNSS position and velocity information can suppress the errors of a strapdown inertial navigation system (SINS) [24]. Second, GNSS positioning information is neededs to solve the carrier acceleration problem [25]. Lastly, GNSS information can be used to calculate the Eotvos correction and normal gravity.
However, vehicles performing gravimetric activities on the ground inevitably encounter complex GNSS observation environments, such as tall urban buildings, tunnels, overpasses, and dense forests. This leads to a decline in the observation quality of GNSS, which directly affects the accuracy of gravity measurement [26]. To eliminate the dependence of land vehicle gravimetry activities on GNSS positioning information, a SINS/velocimeter (VEL) land vehicle gravimetry method was proposed in the literature [27]. This method obtains the carrier acceleration by differentiating the velocity information output by VEL. However, the elevation information of this method is obtained from the integrated navigation, which limits the accuracy of gravity measurement. The literature [28] proposed a land vehicle gravimetry method based on SINS/laser doppler velocimeter (LDV), which uses LDV to output the carrier’s elevation information to further improve gravity measurement accuracy. However, when the carrier travels slowly or encounters severe road turbulence, the optical velocimeter may fail to measure the velocity [29]. The resolution of the measurement results and the experiment’s reliability face severe challenges.
Aiming at the problems of the above methods, this paper proposes a land vehicle gravity anomaly measurement method combining an odometer (OD), geo-information system (GIS), and laser strapdown inertial navigation system. The aim is to free land vehicle gravimetric activities from their heavy reliance on GNSS information, enabling them to maintain the accuracy of gravimetry measurements in extreme environments. The proposed method provides an effective scheme for resource exploration and vehicle gravity assisted navigation.
This paper is organized as follows: Section 2 shows the mathematical model and error equation of strapdown land vehicle gravity anomaly measurement. Section 3 introduces the data processing flow of the SINS/OD/GIS gravity anomaly measurement method. Section 4 describes in detail the comparison test between the gravity anomaly measurement system and the land vehicle gravity measurement. Section 5 compares the accuracy of gravity anomaly measurement of the three methods. See Section 6 for conclusions.

2. Strapdown Land Vehicle Gravity Anomaly Measurement Principle

The reference coordinate systems commonly used in strapdown land vehicle gravimetry are the inertial coordinate system ( i -frame), the Earth coordinate system ( e -frame), the navigation coordinate system ( n -frame), the actual navigation coordinate system ( n -frame), the odometer coordinate system ( m -frame), and the body coordinate system ( b -frame). In this paper, the local geography coordinate system (E-N-U) is chosen as the navigation coordinate system.

2.1. Mathematical Model

The expression of gravity measurement can be obtained from the inertial navigation specific force equation:
g n = v ˙ e n n C b n f b + [ ( 2 ω i e n + ω e n n ) ] v e n n
where g n is the gravity, v ˙ e n n is the acceleration of the vehicle, C b n is the direction cosine matrix from b -frame to n -frame, f b is the specific force output of the accelerometer, ω i e n is the projection of the Earth’s rotational angular velocity in the n -frame, ω e n n is the projection of the angular velocity of the vehicle relative to Earth in the n -frame, v e n n is the velocity of the carrier, and stands for transforming the vector into a skew symmetric matrix.
The detailed expressions of g n , C b n , v e n n , ω i e n and ω e n n are as follows:
g n = [ g E g N g U ] T , C b n = [ c υ c ψ + s υ s ψ s θ s ψ c θ s υ c ψ c ψ s θ c υ s ψ + s υ c ψ s θ c ψ c θ s υ s ψ c υ c ψ s θ s υ c θ s θ c υ c θ ] , v e n n = [ v E v N v U ] T , ω i e n = [ 0 ω i e cos L ω i e sin L ] T , ω e n n = [ v N R M + h v E R N + h v E R N + h tan L ] T .
where g E , g N and g U are eastward gravity, north gravity and upward gravity respectively, c ( · ) and s ( · ) stand for calculating values of sine and cosine, υ is the rolling angle of vehicle, ψ is the yaw angle of the vehicle, θ is the pitch angle of the vehicle, ω i e is the upward component of the Earth’s rotational angular velocity. L is the latitude, R N and R M are the curvature radii of the prime upward and meridian, h is the elevation of the vehicle, and v N , v E and v U are the land vehicle’s north velocity, eastward velocity and upward velocity respectively.
Gravity is considered to be the sum of the normal gravity γ n = [ 0 0 γ ] T ( γ denote the upward component of the normal gravity) and the gravity disturbance δ g n [30]:
g n = γ n + δ g n
The upward component of the normal gravity can be directly calculated using the following equations:
γ = γ 0 ( 1 + h / R M R N ) 2 , γ 0 = 9.780318 ( 1 + 5.3024 × 10 3 sin 2 L 5.9 × 10 6 sin 2 2 L ) m / s 2
So, measuring gravity is equivalent to measuring gravity disturbance. The expression of gravity disturbance can be obtained by combining Equations (1) and (2) [31]:
δ g n = v ˙ e n n C b n f b + [ ( 2 ω i e n + ω e n n ) ] v e n n γ n
In n -frame, gravity disturbance can be decomposed into north component of gravity disturbance, east component of gravity disturbance and gravity anomaly. This paper mainly studies the measurement of gravity anomaly. The other two components are measured in the same way, and the accuracy is closely related to the precision of the gyroscope in SINS.
By expanding Equation (3) in the n -frame, the expression of gravity anomaly is as follows:
δ g U = v ˙ U f U + 2 ω i e v E cos L + v E 2 R N + h + v N 2 R M + h γ
where δ g U is the upward component of the gravity disturbance, v ˙ U is the upward component of the vehicle’s velocity, and f U is the upward component of the specific force.
In geodesy and navigation, the Eotvos correction is denoted as δ a E .
δ a E = 2 ω i e v E cos L + v E 2 R N + h + v N 2 R M + h
Therefore, Equation (4) can be further expressed as:
δ g U = v ˙ U f U + δ a E γ
The relationship of gravity, normal gravity and gravity disturbance is shown in Figure 1. η and ξ are the meridian component and prime unitary component.

2.2. Error Model

Applying the variational operation to Equation (3), the error model of gravity anomaly measurement is obtained:
d δ g n = δ v ˙ e n n f n ϕ C b n δ f b v e n n ( 2 δ ω i e n + δ ω e n n ) + [ ( 2 ω i e n + ω e n n ) ] δ v e n n δ γ n
where d δ g n is the measurement error of gravity anomaly, δ v ˙ e n n is the vehicle’s acceleration measurement error, f n is the projection of the specific force on the n -frame, δ f b is the specific force measurement error, ϕ is the attitude error, δ ω i e n is the Earth rotation angular velocity calculation error, δ ω e n n is the calculation error of the vehicle relative to the earth angular velocity, δ v e n n is the vehicle’s velocity error, and δ γ n is the normal gravity calculation error.
Expand Equation (5) into component form:
[ d δ g E d δ g N d δ g U ] = [ δ v ˙ E δ v ˙ N δ v ˙ U ] [ 0 f U f N f N 0 f E f N f E 0 ] [ ϕ E ϕ N ϕ U ] [ δ f E δ f N δ f U ] [ 0 v U v N v N 0 v E v N v E 0 ] [ δ v N R M + h 2 ω i e δ L cos L + δ v E R N + h 2 ω i e δ L sin L + δ v E tan L R N + h + δ L v E tan L R N + h ] + [ 0 2 ω i e sin L v E tan L R N + h 2 ω i e cos L + v E R N + h 2 ω i e sin L + v E tan L R N + h 0 v N R M + h 2 ω i e cos L v E R N + h v N R M + h 0 ] [ δ v E δ v N δ v U ] [ 0 0 δ γ ]
where d δ g E , d δ g N , and d δ g U are the measurement errors of eastward gravity disturbance, north gravity disturbance and gravity anomaly, respectively, and δ v ˙ E , δ v ˙ N , and δ v ˙ U are the errors of land vehicle’s eastward acceleration, north acceleration, and upward acceleration respectively. f E and f N are the eastward specific force and north specific force, δ f E , δ f N , and δ f U are the eastward specific force measurement error, north specific force measurement error, and upward specific force measurement error, respectively, ϕ E , ϕ N , and ϕ U are the eastward attitude error, north attitude error, and upward attitude error, respectively, δ v E , δ v N , and δ v U are the measurement error of eastward velocity, north velocity and upward velocity, respectively, δ L is the latitude error, and δ γ is the upward normal gravity calculation error.
The upward component of Equation (6) is the measurement error model of gravity anomaly:
d δ g U = δ v ˙ U ( f N ϕ E + f E ϕ N ) δ f U ( 2 ω i e cos L v E R N + h ) δ v E v E δ v E R N + h + v N δ v N R M + h + 2 ω i e v E δ L cos L + v E δ v E R N + h δ γ = δ v ˙ U ( f N ϕ E + f E ϕ N ) δ f U 2 ω i e cos L δ v E + 2 v E δ v E R N + h + 2 ω i e v E δ L cos L δ γ
The influence of latitude on gravity anomaly measurement accuracy is denoted as d δ g ( L ) , which can be obtained from Equation (7):
d δ g ( L ) = 2 ω i e v E δ L cos L
It is assumed that in a land vehicle gravity measurement test, the east velocity of the vehicle is 10 m/s, and the maximum latitude positioning error is 100 m. As seen from Equation (8), the gravity anomaly measurement error caused by the latitude positioning error in this test is no more than 0.1 mGal.
In this paper, the upward gravity gradient of the Earth is taken as 0.3 mGal/m. So, δ γ can be further expressed as:
δ γ = 0.3 δ h
where δ h is the error of elevation.
According to Equation (9), an elevation error of only 3.3 m will cause a gravity anomaly measurement error of 1 mGal.
Similarly, the influences of IMU error, velocity error and acceleration error on gravity anomaly measurement are respectively denoted as d δ g ( SINS ) , d δ g ( v ) , and d δ g ( a c c ) . As can be seen from Equation (7):
d δ g ( SINS ) = ( f N ϕ E + f E ϕ N ) δ f U d δ g ( v ) = 2 ω i e cos L δ v E + 2 v E δ v E R N + h d δ g ( a c c ) = δ v ˙ U
So, Equation (7) can be abbreviated as:
d δ g U = d δ g ( a c c ) + d δ g ( SINS ) + d δ g ( L ) + d δ g ( v ) δ γ
In summary, the measurement accuracy of the land vehicle gravity anomaly is related to the measurement accuracy of carrier acceleration, IMU, latitude positioning, velocity measurement and elevation measurement. The latitude positioning accuracy has little effect on the gravity anomaly measurement accuracy. The accuracy of elevation has great influence on the accuracy of the gravity anomaly measurement.

3. SINS/OD/GIS Gravity Anomaly Measurement Method

As previously discussed, high-precision positioning, attitude error correction, and elevation information have significant positive implications for enhancing the accuracy of land vehicle gravity measurements. Firstly, the latitude, east velocity, and north velocity obtained by integrated navigation are essential information for gravity signal extraction. Secondly, the strapdown land vehicle gravity measurement adopts the mathematical model under the local navigation coordinate system. However, since SINS is firmly connected to the carrier, the output is a specific force of b -frame. Gravity measurement can only be carried out after further calculating the direction cosine matrix C b n to obtain the upward specific force f U . This paper uses an error-state Kalman filter to fuse SINS and OD information. The odometer’s velocity output suppresses SINS errors, improving the accuracy of latitude positioning, combined navigation velocity, and specific force measurement.
In addition, the gravity disturbance information directly extracted by Equation (4) contains high-frequency noise, which a suitable filter must filter before the gravity anomaly can be obtained.
So, gravity anomaly measurement can be regarded as a combination of navigation and gravity extraction.
The specific flow is summarized as follows:
(1)
SINS performs an inertial navigation solution, and the velocity of OD is observed for integrated navigation. The carrier’s specific force, attitude, velocity, and position information are calculated using SINS/OD integrated navigation.
It should be noted that the original OD data can be used for integrated navigation only after being converted to the n -frame using the direction cosine matrices C m b and C b n in turn. The conversion equation for the odometer velocity is as follows:
v O D n = C b n C m b v O D m
where v O D n is the projection of the odometer velocity in the n -frame, and C m b is the direction cosine matrix from m -frame to b -frame.
In real conditions, SINS has an installation error angle, and the calibration coefficient of the odometer is also affected by tire wear and tire pressure changes. Combined with Equation (11), the actual output of the odometer can be obtained after ignoring a small amount:
v ˜ O D n = ( I κ ) C b n ( I η ) ( 1 + δ K O D ) v O D m = v O D n + v O D n κ + C b n ( v O D m ) η + C b n v O D m δ K O D
where v ˜ O D n is the odometer taking into account the actual output of the n -frame after various errors, I is the unit matrix, κ is the misalignment angle between the n -frame and the n -frame, η = [ α θ 0 α ψ ] T is the misalignment angle between the m -frame and the b -frame, α θ and α ψ are SINS installation error angles, and δ K O D is the odometer calibration coefficient error.
SINS/OD integrated navigation selects an 18-dimensional error state vector:
X = [ ϕ n δ v n δ P n ε b b σ ] T
where X is the status vector of the SINS/OD integrated navigation system, ϕ n is the projection of ϕ onto the n -frame, δ v n is the projection of velocity error in n -frame, δ P n is the projection of position error in n -frame, ε b is the gyro zero drift, b is the constant zero bias of the specific force sensor, and σ = [ α θ α ψ δ K O D ] T is the augmented error vector, which can be treated as random constant states, that is,
σ ˙ = 0
Differential equations for navigation error state vectors in Equation (13) are as follows:
ϕ ˙ n = ϕ n ω i n n C b n ε b , δ v ˙ n = f n ϕ n + v n ( 2 δ ω i e n + δ ω e n n ) ( 2 ω i e n + ω e n n ) δ v n + C b n b , δ P ˙ n = M p v δ v n + M p p δ P , M p v = [ 0 1 R M + h 0 sec L R N + h 0 0 0 0 1 ] , M p p = [ 0 0 v N ( R M + h ) 2 v E sec L tan L R N + h 0 v E sec L ( R N + h ) 2 0 0 0 ] .
The differential equation of IMU error state vectors in Equation (13) are as follows:
ε ˙ b = 0 , ˙ b = 0 .
Combining differential equations of error state vectors, we can obtain the 18-dimensional augmented state equations of Kalman filter as follows:
X ˙ = A X + B w
where X ˙ is the state error, A is the state transition matrix, B is the noise-driven matrix, and w ( 0 , q ) is the process noise.
With the OD output velocity as the true value, the SINS output velocity can be expressed as:
v ˜ SINS n + δ v n = v ˜ O D n
where v ˜ SINS n is the actual output velocity of SINS, and v ˜ O D n is the OD indicating the actual output velocity.
Further, by combining the Equations (12) and (15), we can obtain:
δ v n = v ˜ O D n v ˜ SINS n = v O D n κ + C b n ( v O D m ) η + C b n v O D m δ K O D
The difference between the output velocity of SINS and that of the odometer is taken as a measurement, and the measurement equation of the system is as follows:
Z = H X + V Z = [ v ˜ S I N S n v ˜ O D n ] = δ v n = [ δ v E δ v N δ v U ] T H = [ v O D n I 3 × 3 0 3 × 9 B b n v O D ]
where Z is the observation vector, H is the coefficient matrix of the observed vector, and B b n is the equivalent matrix of C b n . The specific expression of B b n is as follows:
B b n = ( C b n ( : , 3 ) C b n ( : , 1 ) C b n ( : , 2 ) )
where C b n ( : , i ) means that all the elements in column i of C b n are taken out to form a new number table.
Combining Equations (14) and (17), we can obtain the 18-dimensional error state equations of Kalman filter as follows:
X ˙ = A X + B w Z = H X + V
(2)
The commercial geographic information software ArcGIS 10.8 was used to extract the elevation information of the survey line. See Figure 2 for the specific process.
First, the longitude and latitude data of SINS/OD integrated navigation are stored. Then, the latitude and longitude data are input into ArcGIS 10.8, and the DEM data of the city where the survey line is located is selected. Second, the elevation information output by ArcGIS is stored in an Excel file. Lastly, the line elevation information stored in Excel is input into the gravity anomaly extraction software.
The accuracy of the elevation information obtained by this method is related to the accuracy of the DEM file.
(3)
The original gravity anomaly is calculated by Equation (4).
Figure 3 illustrates the functions of various sensors in the measurement of gravity anomaly information. OD provides velocity information and is used to calculate the kinematic acceleration of the vehicle. After solving SINS and OD integrated navigation, the high-precision attitude and specific force can be output for upward specific force calculation. The v E , v N , R M , R N , L , h , and ω i e output by SINS/OD integrated navigation are used to calculate the Eotvos correction. The normal gravity is calculated by using the elevation output of ArcGIS 10.8 and L , R M , and R N of the SINS/OD integrated navigation solution.
Further, the original gravity anomaly information containing high-frequency noise is filtered by using an FIR filter to obtain the gravity anomaly information.
(4)
The external coincidence accuracy evaluation is carried out.
This paper uses gravity anomalies calculated by the EGM2190 model as reference values, and REF represents the data.
The calculation formula for the accuracy evaluation of test line data is as follows:
τ = i = 1 N ( Δ g i R E F i ) 2 N
where τ is the external coincidence accuracy of land vehicle gravity measurement, Δ g i is the measured gravity anomalies at point i , R E F i is the gravity values calculated by the EGM2190 model at point i , and N is the measure the total number of sampling points on the test line.
Figure 4 is the data processing flow chart of the SINS/OD/GIS gravity measurement method.

4. Experiments

In March 2024, the research group carried out a land vehicle dynamic gravity measurement test on a highway in Tianjin. This experiment aimed to verify the accuracy of the proposed method and the reliability of the gravity measurement system. The method proposed in this paper was used as the experimental group; the traditional SINS/GNSS method and SINS/OD method were used as the control tests. The specific process of these two methods can be seen in reference [32].
We equipped an SUV with a self-developed dual-axis rotating laser strapdown inertial navigation system, odometer, GNSS positioning equipment, and other necessary equipment. The role of GNSS is to provide a positioning benchmark. The navigation computer, battery, odometer, and strapdown inertial navigation system constitute the land vehicle gravity anomaly measurement system shown in Figure 5. Table 1 shows the detailed parameters of the above sensors.
Figure 6 shows the track of the gravity survey area. The blue line represents the path of the dynamic gravity measurement experiment. In this paper, the relative coordinates are used to represent the latitude of the X-axis and the longitude of the Y-axis. The driver steered the vehicle to run smoothly, the velocity was maintained at 19.6 km/h, and the vehicle did not undergo significant acceleration or deceleration maneuvers. In addition, the GNSS observation signal of the whole gravity sampling section was good.

5. Results and Discussion

5.1. Data Processing

To evaluate whether the positioning accuracy of SINS/OD integrated navigation meets the requirements of gravity measurement, the positioning results of SINS/GNSS integrated navigation are used as reference values to verify the degree of agreement between the positioning results of SINS/OD and the reference values. The results of SINS/OD integrated navigation are shown in Figure 7 and Figure 8.
Figure 7 reflects the consistency of the positioning results of SINS/OD integrated navigation and SINS/GNSS integrated navigation. In addition, as can be seen from Figure 8, the latitude error of SINS/OD integrated navigation is less than 33.3 m. According to Equation (8), the gravity measurement error caused by the latitude positioning error of SINS/OD combined navigation can be ignored.
In fact, the elevation accuracy of GIS is not only related to the accuracy of the DEM file, but also to the horizontal positioning accuracy of SINS/OD integrated navigation. Considering that the horizontal positioning error of SINS/OD integrated navigation will directly affect the elevation accuracy, in this paper, only the data between 1280 s and 1920 s of the test are selected for gravity extraction. During this period, the horizontal positioning error of integrated navigation is small, and the elevation accuracy of GIS output is high.
Figure 9 shows the elevation error of GIS, taking the elevation output of SINS/GNSS integrated navigation as the reference value.
As Figure 9 shows, the maximum error of the GIS output elevation is better than 1.25 m, which meets the requirements of gravity measurement for elevation accuracy.

5.2. Test Results

Using the 100s FIR low-pass filtering result as an example, Figure 10 shows the gravity anomaly results of the three methods in this test.
To improve the efficiency of gravity measurement, the gravity reference point was not established by ground-static gravity measurement. Instead, the gravity anomaly calculated by the EGM2190 model is used as a reference value to evaluate the accuracy of external coincidence.
Table 2 shows the statistical results of external coincidence accuracy. The proposed method has the same accuracy compared to the traditional SINS/GNSS method. Compared with the SINS/OD method, the outer coincidence accuracy of the proposed method is improved from 2.35 mGal to 1.65 mGal, which is a 30% increase in accuracy.

5.3. Discussion

The test results demonstrate that the proposed method achieves higher accuracy in gravity anomaly measurement compared to the SINS/OD method. The primary difference between the two methods lies in the calculation of normal gravity. Specifically, the SINS/OD method computes normal gravity using the elevation, latitude, curvature radius of the prime vertical circle, and curvature radius of the meridian circle derived from integrated navigation. In contrast, the SINS/OD/GIS method utilizes the elevation obtained from the GIS, along with other information provided by the SINS/OD integrated navigation to calculate normal gravity. To further elucidate the relationship between the gravity anomaly measurement performance and the elevation of the proposed method, we analyzed the elevation error of all three methods along the measurement line. The results are presented in Figure 11.
As depicted in Figure 11, the elevation error of the GNSS is very small. In contrast, the elevation error of the SINS/OD integrated navigation system exhibits divergence. The elevation error of the GIS, however, does not diverge over time, with a maximum elevation error of 1.24 m. According to Equation (9), in the SINS/OD/GIS method, the measurement error attributable to elevation is no more than 0.38 mGal. These findings indicate that the precision of gravity measurements can be significantly enhanced by utilizing the high-precision, non-divergent elevation information provided by GIS.

6. Conclusions

It is difficult to maintain the accuracy of gravity anomaly measurement in extreme environments when GNSS is used to assist the land vehicle dynamic gravity anomaly measurement. To solve this problem, a new method combining an odometer and digital map is proposed. From the dynamic test results, the external coincidence accuracy of the proposed method is better than 2 mGal, and the accuracy is comparable to the SINS/GNSS method. The elevation information obtained by GIS does not diverge with time and vehicle mileage, but the accuracy is limited. Further research should continue to focus on how to achieve autonomous and accurate elevation measurement, to ensure the measurement accuracy of gravity anomaly information. This paper can provide technical support for high-precision positioning of electric land vehicles and resource exploration.

Author Contributions

Conceptualization, J.Z. and Z.Z.; Methodology, K.Z.; Software, K.Z.; Writing—original draft preparation, K.Z.; Writing—review and editing, J.Z. and K.Z.; Funding acquisition Z.Z. and J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China (Grant No. 62305393) and the National Key Basic Research Strengthen Foundation of China (Grant No. 2021-JCJQ-JJ-0871).

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to confidentiality.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pistorio, A.; Greco, F.; Currenti, G.; Napoli, R.; Sicali, A.; Del Negro, C.; Fortuna, L. High-Precision Gravity Measurements Using Absolute and Relative Gravimeters at Mount Etna (Sicily, Italy). Ann. Geophys. 2011, 54, 5348. [Google Scholar] [CrossRef]
  2. Wang, H.; Wang, K.; Xu, Y.; Tang, Y.; Wu, B.; Cheng, B.; Wu, L.; Zhou, Y.; Weng, K.; Zhu, D.; et al. A Truck-Borne System Based on Cold Atom Gravimeter for Measuring the Absolute Gravity in the Field. Sensors 2022, 22, 6172. [Google Scholar] [CrossRef]
  3. Wu, B.; Li, D.; Zhou, Y.; Zhu, D.; Zhao, Y.; Qiao, Z.; Cheng, B.; Niu, J.; Guo, X.; Wang, X.; et al. Construction of a Test Field for Relative Gravimeters in a Cave With a Cold Atom Gravimeter. IEEE Sens. J. 2024, 24, 9536–9544. [Google Scholar] [CrossRef]
  4. Zhang, X.; Yan, S.; Li, Q.; Zhu, L.; Yang, J. Mobile gravity surveys in the field based on vehicle-mounted atom interferometer. Chin. J. Sci. Instrum. 2023, 44, 96–103. [Google Scholar] [CrossRef]
  5. Bidel, Y.; Zahzam, N.; Blanchard, C.; Bonnin, A.; Cadoret, M.; Bresson, A.; Rouxel, D.; Lequentrec-Lalancette, M.F. Absolute Marine Gravimetry with Matter-Wave Interferometry. Nat. Commun. 2018, 9, 627. [Google Scholar] [CrossRef]
  6. Bidel, Y.; Zahzam, N.; Bresson, A.; Blanchard, C.; Cadoret, M.; Olesen, A.V.; Forsberg, R. Absolute Airborne Gravimetry with a Cold Atom Sensor. J. Geod. 2020, 94, 20. [Google Scholar] [CrossRef]
  7. Ning, J.S.; Huang, M.T.; Ouyang, Y.Z.; Deng, K.L. Development of Marine and Airborne Gravity Measurement Technologies. Hydrogr. Surv. Charting 2014, 34, 67–72+76. [Google Scholar]
  8. Hu, P.H.; Zhao, M.; Huang, H.; Liu, D.B.; Tang, J.H.; Wei, C. Review on the Development of Airborne/Mairne Gravimetry Instruments. Navig. Position. Timing 2017, 4, 10–19. [Google Scholar] [CrossRef]
  9. Xiu, R.; Guo, G.; Xue, Z.B.; Li, D.M.; Li, H.B. Technical Current Situation and New Application of Marine/aviation Gravimeter. Navig. Control. 2019, 18, 35–43. [Google Scholar]
  10. Jekeli, C. Balloon Gravimetry Using GPS and INS. IEEE Aerosp. Electron. Syst. Mag. 1992, 7, 9–15. [Google Scholar] [CrossRef]
  11. Zhang, Z.; Li, J.; Zhang, K.; Yu, R. Experimental Study on Underwater Moving Gravity Measurement by Using Strapdown Gravimeter Based on AUV Platform. Mar. Geod. 2021, 44, 108–135. [Google Scholar] [CrossRef]
  12. Luo, K.; Cao, J.; Wang, C.; Cai, S.; Yu, R.; Wu, M.; Yang, B.; Xiang, W. First Unmanned Aerial Vehicle Airborne Gravimetry Based on the CH-4 UAV in China. J. Appl. Geophys. 2022, 206, 104835. [Google Scholar] [CrossRef]
  13. Lu, B.; Xu, C.; Li, J.; Zhong, B.; van der Meijde, M. Marine Gravimetry and Its Improvements to Seafloor Topography Estimation in the Southwestern Coastal Area of the Baltic Sea. Remote Sens. 2022, 14, 3921. [Google Scholar] [CrossRef]
  14. Hu, Q.; Xu, X.Y.; Zhao, Y.Q.; Ding, Z.F. Downward Continuation of Airborne Gravity Anomaly based on Lagrange Mean-Value Theorem. J. Geod. Geodyn. 2021, 41, 95–100. [Google Scholar] [CrossRef]
  15. Wang, X.; Liu, J.; Sun, P.; Chen, K.; Li, J. Rapid detection of urban underground cavity based on vehicle gravity measurement platform. J. Jilin Univ. (Earth Sci. Ed.) 2019, 49, 838–845. [Google Scholar] [CrossRef]
  16. Chiang, K.; Lin, C.; Kuo, C. A Feasibility Analysis of Land-Based SINS/GNSS Gravimetry for Groundwater Resource Detection in Taiwan. Sensors 2015, 15, 25039–25054. [Google Scholar] [CrossRef] [PubMed]
  17. Li, X.; Jekeli, C. Ground-Vehicle INS/GPS Vector Gravimetry. Geophysics 2008, 73, I1–I10. [Google Scholar] [CrossRef]
  18. Yu, R.; Cai, S.; Wu, M.; Cao, J.; Zhang, K. An SINS/GNSS Ground Vehicle Gravimetry Test Based on SGA-WZ02. Sensors 2015, 15, 23477–23495. [Google Scholar] [CrossRef] [PubMed]
  19. Li, D.; Zeng, X.; Ma, J.; Wang, W.; Zhang, H. An Airborne Gravimeter Based on the Rate Azimuth Platform INS. Navig. Control. 2014, 13, 7–10. [Google Scholar]
  20. Li, D.; Guo, G.; Xue, Z.; Wang, W. Ground Test Results of a Strap-down Gravimeter Used on the Movable Platform. Navig. Position. Timing 2015, 2, 59–62. [Google Scholar] [CrossRef]
  21. Yu, R.; Qiu, X.; Cao, J.; Cai, S.; Lu, S.; Xu, X.; Wang, L. Improving Land Vehicle Gravimetry Using a New SINS/GNSS/VEL Method. IOP Conf. Ser. Earth Environ. Sci. 2020, 513, 012054. [Google Scholar] [CrossRef]
  22. Yu, R.; Wang, L.; Cao, J.; Cai, S.; Wu, M.; Xu, X. An Improved Method Using SINS/PPP Method for Land Vehicle Gravimetry. J. Phys. Conf. Ser. 2020, 1550, 042052. [Google Scholar] [CrossRef]
  23. Xiao, Y.; Xia, Z.R. Determination of Moving-Base Acceleration in Airborne Gravimetry. Chin. J. Geophys. 2003, 46, 62–67. [Google Scholar]
  24. Hwang, Y.; Jeong, Y.; Kweon, I. Online Misalignment Estimation of Strapdown Navigation for Land Vehicle under Dynamic Condition. Int. J. Automot. Technol. 2021, 22, 1723–1733. [Google Scholar] [CrossRef]
  25. Zhang, K. Research on Methods of Airborne Gravimetry Based on SINS/DGPS. Ph.D. Thesis, National University of Defense Technology, Changsha, China, 2007. [Google Scholar]
  26. Yu, R.; Wu, M.; Cao, J.; Zhang, K.; Cai, S. A New Method of GNSS Fault Data Detection for Strapdown Land Vehicle Gravimetry. In Proceedings of the 2018 IEEE International Conference on Applied System Invention (ICASI), Chiba, Japan, 13–17 April 2018; pp. 299–302. [Google Scholar]
  27. Yu, R.; Wu, M.; Zhang, K.; Cai, S.; Cao, J.; Wang, M.; Wang, L. A New Method for Land Vehicle Gravimetry Using SINS/VEL. Sensors 2017, 17, 766. [Google Scholar] [CrossRef] [PubMed]
  28. Wei, G.; Yang, Z.; Gao, C.; Zhou, J.; Yu, X.; Luo, H.; Deng, B.; Zhou, W.; Cheng, J. Strapdown vehicle autonomous gravimetry method based on two-dimensional laser doppler velocimeter. Infrared Laser Eng. 2023, 52, 339–346. [Google Scholar]
  29. Zhang, Z.; Zhou, Z.; Chen, H. Vehicle Autonomous Combination Positioning and Orientation Technology; National Defense Industry Press: Beijing, China, 2021; ISBN 978-7-118-12233-6. [Google Scholar]
  30. Tie, J.; Cao, J.; Chang, L.; Cai, S.; Wu, M.; Lian, J. A Model of Gravity Vector Measurement Noise for Estimating Accelerometer Bias in Gravity Disturbance Compensation. Sensors 2018, 18, 883. [Google Scholar] [CrossRef]
  31. Yu, R. Research on Key Technologies for Strapdown Ground Vehicle Gravimetry. Ph.D. Thesis, National University of Defense Technology, Changsha, China, 2020. [Google Scholar]
  32. Li, X.; Zhang, Z.; Zhou, Z.; Chang, Z.; Zhao, Z. Vehicle gravity anomaly measurement methods based on SINS/OD/altimeter. J. Beijing Univ. Aeronaut. Astronaut. 2023, 1–15. [Google Scholar] [CrossRef]
Figure 1. Relationship of gravity, normal gravity and gravity disturbance at a certain point in local navigation coordinate system.
Figure 1. Relationship of gravity, normal gravity and gravity disturbance at a certain point in local navigation coordinate system.
Wevj 15 00368 g001
Figure 2. The process of elevation extraction by ArcGIS 10.8.
Figure 2. The process of elevation extraction by ArcGIS 10.8.
Wevj 15 00368 g002
Figure 3. The functions of various sensors in the measurement of gravity anomaly.
Figure 3. The functions of various sensors in the measurement of gravity anomaly.
Wevj 15 00368 g003
Figure 4. Data processing flow of the SINS/OD/GIS method.
Figure 4. Data processing flow of the SINS/OD/GIS method.
Wevj 15 00368 g004
Figure 5. Land vehicle gravity anomaly measurement system.
Figure 5. Land vehicle gravity anomaly measurement system.
Wevj 15 00368 g005
Figure 6. Track of the gravity survey area.
Figure 6. Track of the gravity survey area.
Wevj 15 00368 g006
Figure 7. Position plot of SINS/OD computation compared with the reference data.
Figure 7. Position plot of SINS/OD computation compared with the reference data.
Wevj 15 00368 g007
Figure 8. Latitude error plot of SINS/OD computation compared with the reference data.
Figure 8. Latitude error plot of SINS/OD computation compared with the reference data.
Wevj 15 00368 g008
Figure 9. Elevation error.
Figure 9. Elevation error.
Wevj 15 00368 g009
Figure 10. Comparisons of filtered gravity anomaly.
Figure 10. Comparisons of filtered gravity anomaly.
Wevj 15 00368 g010
Figure 11. Elevation error comparison of the three methods.
Figure 11. Elevation error comparison of the three methods.
Wevj 15 00368 g011
Table 1. Sensor parameter.
Table 1. Sensor parameter.
ItemsParameterPerformance Index
Inertial measurement unitSampling frequency200 Hz
Gyroscope zero bias0.003 (°)/h
Accelerometer zero bias10 μg
OdometerCalibration coefficient0.9986 m/pulse
GPS parameterSampling frequency1 Hz
Positioning accuracy5 m
Table 2. Gravity anomaly measurement external coincidence results (unit: mGal).
Table 2. Gravity anomaly measurement external coincidence results (unit: mGal).
ItemsMaxMinMeanRMSAccuracy
SINS/GNSS2.56−2.231.361.491.49
SINS/OD3.12−5.742.092.352.35
SINS/OD/GIS2.86−3.571.321.651.65
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, K.; Zhao, J.; Zhang, Z. A Method for Land Vehicle Gravity Anomaly Measurement Combining an Inertial Navigation System, Odometer, and Geo-Information System. World Electr. Veh. J. 2024, 15, 368. https://doi.org/10.3390/wevj15080368

AMA Style

Zhang K, Zhao J, Zhang Z. A Method for Land Vehicle Gravity Anomaly Measurement Combining an Inertial Navigation System, Odometer, and Geo-Information System. World Electric Vehicle Journal. 2024; 15(8):368. https://doi.org/10.3390/wevj15080368

Chicago/Turabian Style

Zhang, Kefan, Junyang Zhao, and Zhili Zhang. 2024. "A Method for Land Vehicle Gravity Anomaly Measurement Combining an Inertial Navigation System, Odometer, and Geo-Information System" World Electric Vehicle Journal 15, no. 8: 368. https://doi.org/10.3390/wevj15080368

APA Style

Zhang, K., Zhao, J., & Zhang, Z. (2024). A Method for Land Vehicle Gravity Anomaly Measurement Combining an Inertial Navigation System, Odometer, and Geo-Information System. World Electric Vehicle Journal, 15(8), 368. https://doi.org/10.3390/wevj15080368

Article Metrics

Back to TopTop