Abstract
In development of a rotating accelerometer gravity gradiometer (RAGG), it is difficult for us to distinguish the measurement signals and error components in the RAGG output without a verified and correctness RAGG analytical model. In addition, many key techniques, such as RAGG analytical model validation, online error compensation, post motion error compensation, are difficult to be verified. RAGG numerical model can provide validation platform for solving all the above problems, which can speed up the development of the RAGG. In this study, based on the principle and configuration of the RAGG, we synthetically consider almost all the error factors, such as circuit gain mismatch, installation errors, accelerometer scale-factor imbalance, and accelerometer second-order error coefficients, construct a parameters adjustable RAGG numerical model. In multi-frequency gravitational gradient simulation experiment, we use the RAGG numerical model simulating the situation that a test mass rotates about the RAGG with time-varying angular velocity to generate multi-frequency gravitational gradient excitations; the experiment results are consistent with the theoretical ones; the RAGG numerical model can recur some phenomenons of a actual RAGG.
You have full access to this open access chapter, Download conference paper PDF
Similar content being viewed by others
Keywords
- Rotating accelerometer gravity gradiometer
- Numerical model
- Center gravitational gradients
- Virtual rotating accelerometer gravity gradiometer
1 Introduction
Airborne gravity gradiometry is an advanced technology for surveying a gravity field; it acquires gravity field information with high efficiency and high spatial resolution. Compared with gravity information, the gravity gradient tensor provides more information on the field source such as orientation, depth, and shape (Tang and Hu 2018; Yan and Ma 2015). There are many different types of gravity gradiometer: rotating accelerometer gravity gradiometers, superconducting gravity gradiometers, cold atomic interferometer gravity gradiometers (Paik 2007; Difrancesco 2007; Moody 2011; Liu 2014; Hao 2013). Rotating accelerometer gravity gradiometer (RAGG) was developed by Ernest Metzger of Bell-Aerospace in the 1980s (Heard 1988). But, over the past decade, rotating accelerometer gravity gradiometers are the only commercial moving base gravity gradiometer in the word; all other types are either in fight testing or in a laboratory setting. Companies that operate commercial rotating accelerometer gravity gradiometer systems are: Bell Geospace (Air-FTG), ARKeX (FTGeX), GEDEX (HD-AGG), and FUGRO (AGG-Falcon) (Rogers 2009; Metzger 1977). In this study, we synthetically consider almost all the unperfect factors, for example accelerometer installation errors, accelerometer scale-factor imbalance, etc., and construct a numerical model. The RAGG numerical model is a virtual RAGG with a comprehensive set of precisely adjustable parameters; based on it, many key techniques, for example RAGG automatic calibration, post error compensation, and self-gradient modeling, can be verified.
2 RAGG Numerical Model
A high-precision numerical model of the RAGG is established, as shown in Fig. 1. In the numerical model, each accelerometer has six mounting error parameters: radial distance, initial phase angle, altitude angle, and misalignment error angles. Among them, the radial distance, initial phase angle, and altitude angle determine the mounting position of the accelerometer; the misalignment error angles determine the orientation deviation between the accelerometer measurement frame and the accelerometer nominal frame of the actual mounting position. The detailed definition of the mounting error parameters can be found in literature (Yu and Cai 2019). Moreover, each accelerometer has nine other output model parameters: zero bias (K j0), linear scale factors (K j1), second-order error coefficients (K j2, K j4, K j5, K j6, K j7, K j8), and current to voltage gain (k jV∕I). We use a test mass to produce gravitational gradients to excite the RAGG. The specific force of the accelerometer A j in the RAGG numerical model is given by:
Where f cmm is the specific force of the RAGG; \({{\dot {\boldsymbol {\omega }} }_{im}}\) is the angular acceleration of the RAGG with respect to the inertial frame; ω im is the angular velocity of the RAGG with respect to the inertial frame; G is the gravitational constant; \({{{\boldsymbol {r}}_{{o_m}{A_j}}}}\) is the position vector of accelerometer A j in the RAGG measurement frame; m represents the weight of the test mass; and A j S is the position vector from accelerometer A j to the test mass. If the test mass is not a point mass, the gravitational acceleration that the RAGG accelerometers undergo produced by the test mass can be calculated using finite element analysis. In addition, the test mass can be in motion with respect to the RAGG; in this case, A j S is time varying. The specific forces of accelerometer A j in the accelerometer nominal frame of the actual mounting position (f jx, f jy, f jz) can be calculated from:
Where Ï„ jx, Ï„ jy, and Ï„ jz are unit vectors of the accelerometer nominal frame of the actual mounting position in the directions of the x-, y-, and z-axes. The specific forces in the accelerometer measurement frame are:
where C is the transformation matrix from the accelerometer nominal frame of the actual mounting position to the accelerometer measurement frame. To make the numerical model approximate the actual RAGG, we add accelerometer noise to the accelerometer model:
The accelerometer noise f jnoise is simulated by a power spectral density model (Jekeli 2006):
where α and b represent the amplitude and low-frequency growth of the red noise, and ω T denotes the amplitude of the white noise.
Figure 1b is the program flow of the RAGG numerical model. Firstly, the RAGG simulation parameters are set up, including test masses parameters, RAGG rotating disk parameters, accelerometer mounting parameters, accelerometer model parameters, RAGG motion parameters, etc. Then substituting the parameters into the formula (1)–(3) calculates the specific force in accelerometer measurement frame at time t. According to the formula (4), calculating the output voltage of the RAGG accelerometer, the RAGG output before demodulation at time t is calculated by: G out(t) = V 1(t) + V 2(t) − V 3(t) − V 4(t). The above process is repeated until time t is equal to the simulation duration time. Finally, the RAGG output data is input to the quadrature amplitude modulation (QAM) demodulator to extract gravitational gradient.
Let Γ xx, Γ xy, Γ xz, Γ yy, and Γ yz represent the five independent gravitational gradient elements at the origin of the RAGG measurement frame. When mass is far enough away from the RAGG, the gravitational acceleration measured by the RAGG accelerometers is a first-order approximation of the gravitational acceleration and gravitational gradient tensor at the center of the rotating disc; in this case, the inline channel measurement and the cross channel measurement of the RAGG approximate Γ xx − Γ xy and Γ xy; otherwise, the inline channel measurement and cross channel measurement of the RAGG are the sum of Γ xx − Γ xy, Γ xy, and high-order gravitational gradient tensor elements. To distinguish between Γ xx − Γ xy, Γ xy and the measurements of the RAGG, Γ xx − Γ xy, Γ xy is called as center gravitational gradients; Γ xx − Γ xy is the inline channel of the center gravitational gradients; Γ xy is the cross channel of the center gravitational gradients. In RAGG analytical model, gravitational acceleration that the RAGG accelerometer undergoes is a first-order taylor approximation of the gravitational acceleration and gravitational gradient tensor at the center of the rotating disc; but in the numerical model, the gravitational accelerations are calculated using Newton’s law of gravitation instead of a linear approximation. Therefore the numerical model are close to the actual RAGG.
3 Multi-frequency Gravitational Gradient Simulation Experiment
In multi-frequency gravitational gradient simulation experiment, a test mass rotates about the RAGG with time-varying angular velocity producing multi-frequency gravitational gradient excitations. Based on the angular velocity of the test mass and its initial coordinate in the RAGG measurement frame, we can obtain the coordinates of the test mass in the RAGG measurement frame at any time. We then calculate the gravitational gradient tensor at the origin of the RAGG measurement frame and then calculate the center gravitational gradients. The RAGG numerical model simulates a perfect RAGG without accelerometer mounting errors, accelerometer scale-factor imbalances, accelerometer second-order error coefficients and accelerometer noise, so we set the accelerometer mounting errors, the accelerometer second-order error coefficients, accelerometer noise parameters to zero. The linear scale factor of the four accelerometers is k j1 = 10 mA/g, the current-to-voltage gain is k jV∕I = 109 ohm, the nominal mounting radius R is 0.1 m, and the rotation frequency of the RAGG disc is 0.25 Hz. Figure 2 shows a point mass of 486 kg with an initial position in the RAGG measurement frame of (1.1, 0.5, 0) and rotating about the RAGG with time-varying angular speed \(\omega (t)=500+400\sin {}(0.0628t)^\circ /\mathrm {{h}}\).
Figure 3 shows the voltage outputs of the four accelerometers in the RAGG numerical model excited by the rotating point mass. Figure 4 shows the output voltage before demodulation of the numerical model. Figures 5 and 6 show the demodulated gravitational gradient comparison among the RAGG model and the center gravitational gradients; we can see that the inline channel and the cross channel of the RAGG numerical model are overlapped with that of the center gravitational gradients and only one curve is displayed; to distinguish the outputs of numerical model and the center gravitational gradients, their differences are calculated and shown in Figs. 7 and 8; the differences are in the order of 10 −1 E, and they are caused by the high-order gravitational gradient tensor.
4 Conclusion
Based on the measurement principle and configuration of the RAGG, we considered the factors of circuit gain mismatch, installation error, accelerometer scale-factor imbalance, and accelerometer second-order error coefficients, then developed a high-precision numerical model. The parameters of the RAGG prototypes are unknowable and uncontrollable, but the parameters of the RAGG numerical model are adjustable and knowable; to some extent, the RAGG numerical model are more suitable for verifying some technique solutions of the RAGG.
References
Difrancesco D (2007) Advances and challenges in the development and deployment of gravity gradiometer systems. In: EGM 2007 international workshop, vol 2007. CSIRO Publishing, p 1. https://doi.org/10.1071/ASEG2007ab034
Hao D (2013) Design of digital balance loop in gravity gradiometer based on ad5791. In: Electronic Measurement & Instruments (ICEMI), 2013 IEEE 11th International Conference on, vol 1. IEEE, pp 396–398. https://doi.org/10.1109/ICEMI.2013.6743105
Heard H (1988) The gravity gradiometer survey system (GGSS). Eos 69(8). https://doi.org/10.1029/88EO00070
Jekeli C (2006) Airborne gradiometry error analysis. Surv Geophys 27(2):257–275. https://doi.org/10.1007/s10712-005-3826-4
Liu H (2014) Design, fabrication and characterization of a micro-machined gravity gradiometer suspension. In: SENSORS, 2014 IEEE, IEEE, pp 1611–1614. https://doi.org/10.1109/ICSENS.2014.6985327
Metzger E (1977) Recent gravity gradiometer developments. In: Guidance and control conference, p 1081. https://doi.org/10.2514/6.1977-1081
Moody M (2011) A superconducting gravity gradiometer for measurements from a moving vehicle. Rev Sci Instrum 82(9):094501. https://doi.org/10.1063/1.3632114
Paik HJ (2007) Geodesy and gravity experiment in earth orbit using a superconducting gravity gradiometer. IEEE Trans Geosci Remote Sens GE-23(4):524–526. https://doi.org/10.1109/TGRS.1985.289444
Rogers MM (2009) An investigation into the feasibility of using a modern gravity gradient instrument for passive aircraft navigation and terrain avoidance. Techreport ADA496707, Air Force Institute of Technology, Wright-Patterson Air Force Base, Ohio, Graduate School of Engineering and Management
Tang J, Hu S (2018) Localization of multiple underwater objects with gravity field and gravity gradient tensor. IEEE Geosci Remote Sens Lett 15(2):247–251. https://doi.org/10.1109/LGRS.2017.2784837
Yan Z, Ma J (2015) Accurate aerial object localization using gravity and gravity gradient anomaly. IEEE Geosci Remote Sens Lett 12(6):1214–1217. https://doi.org/10.1109/LGRS.2015.2388772
Yu M, Cai T (2019) Online error compensation of moving-base rotating accelerometer gravity gradiometer. Rev Sci Instrum 90(7):074501. https://doi.org/10.1063/1.5093078
Acknowledgements
This work is supported by National Key R&D Program of China under Grant No. 2017YFC0601601, 2016YFC0303006 and International Special Projects for Scientific and Technological Cooperation under Grant No. 2014DFR80750.
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Open Access This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license and indicate if changes were made.
The images or other third party material in this chapter are included in the chapter's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the chapter's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
Copyright information
© 2020 The Author(s)
About this paper
Cite this paper
Yu, M., Cai, T. (2020). Numerical Model of Moving-Base Rotating Accelerometer Gravity Gradiometer. In: Freymueller, J.T., Sánchez, L. (eds) 5th Symposium on Terrestrial Gravimetry: Static and Mobile Measurements (TG-SMM 2019). International Association of Geodesy Symposia, vol 153. Springer, Cham. https://doi.org/10.1007/1345_2020_114
Download citation
DOI: https://doi.org/10.1007/1345_2020_114
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-031-25901-2
Online ISBN: 978-3-031-25902-9
eBook Packages: Earth and Environmental ScienceEarth and Environmental Science (R0)