Next Article in Journal
Deep Learning-Based Geomagnetic Navigation Method Integrated with Dead Reckoning
Previous Article in Journal
Learning Implicit Neural Representation for Satellite Object Mesh Reconstruction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Velocity-Related Magnetic Interference Compensation of Unmanned Underwater Vehicle

College of Intelligence Science and Technology, National University of Defense Technology, Changsha 410073, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(17), 4164; https://doi.org/10.3390/rs15174164
Submission received: 29 June 2023 / Revised: 12 August 2023 / Accepted: 22 August 2023 / Published: 24 August 2023

Abstract

:
Underwater magnetic detection based on unmanned underwater vehicles (UUVs) is an important approach for marine exploration and monitoring. Due to the interference of the carrier magnetic field, the detection accuracy is severely limited. To improve the performance of magnetic detection, various compensation methods have been developed. However, the compensation effectiveness of these methods is greatly diminished when UUVs sail at different velocities. In this paper, we propose a velocity-related interference magnetic compensation method that enhances the conventional Tolles–Lawson model. By introducing a velocity-related term and combining the ridge regression with the linear fitting algorithm, we determine the compensation parameters for various UUV velocities. Field experiments are conducted to verify the effectiveness of the proposed method. The results show that the root mean square error of the total magnetic field after compensation is reduced from 246.90 nT to 2.23 nT. Our study demonstrates that the velocity-related compensation method can significantly improve the accuracy of interference magnetic compensation under different UUV velocities, making it applicable to actual underwater magnetic detection.

Graphical Abstract

1. Introduction

As a non-acoustic detection technology, underwater magnetic detection plays an irreplaceable role in geoscience, especially in oceanography [1,2]. With the rapid development of marine observation and exploration, unmanned underwater vehicles (UUVs) have become attractive underwater carriers for a variety of sensors and navigation equipment [3]. UUVs navigate underwater via remote or automatic control, and they have good performance in terms of adaptability, detection range, and operation efficiency [4]. UUVs equipped with magnetometers for underwater magnetic detection hold enormous promise for marine exploration and monitoring, including earthquake monitoring and research, benthonic geological research, marine mineral resource exploration, underwater salvage and search operations, investigation of undersea oil pipelines, detection of underwater magnetic targets, and autonomous underwater navigation of submersibles [5,6,7,8,9].
However, in addition to sensing the magnetic field of the target, the magnetometer on UUVs inevitably detects interference from the carrier, the local geomagnetic field, and other sources during actual detection [10]. The interference magnetic field generated by carriers can significantly distort underwater magnetic measurements. This field mainly consists of three components: the induced magnetic field, the permanent magnetic field, and the eddy current magnetic field, all of which are generated by ferromagnetic materials that are inevitably present in UUVs. It has been reported that on the surface of an underwater vehicle, the induced magnetic field can be of similar magnitude to the geomagnetic field and continuously changes with the vehicle’s attitude [11]. As stated in reference [12], any measurement error exceeding 30 nT can lead to the failure of underwater magnetic detection. In conclusion, the interference of UUVs can considerably affect the accuracy of underwater magnetic detection. Therefore, compensating for carrier magnetic interference is of critical importance.
Notably, the interference magnetic field of UUVs is also affected by their internal electromechanical system’s working state [13]. For instance, UUVs need to operate at various velocities during magnetic detection to adapt to the complex marine environment and control their heading and path [14]. This is achieved through changes in motor speed and battery current, resulting in fluctuations in the interference of the magnetic field [15]. Therefore, the magnetic field interference varies under different operating velocities of UUVs. Hence, it is necessary to compensate for velocity-related interference on the magnetic field to achieve high-precision measurements.
To reduce magnetic detection distortion, researchers have focused on studying carrier magnetic interference compensation [16,17,18,19]. In 1954, Tolles and Lawson analyzed the material composition and structure of aircraft to represent the interference as the sum of three field strengths: permanent field, induced field, and eddy current field, and proposed the Tolles–Lawson model [20]. Leliak et al. established a method to solve the compensation parameters in the interference magnetic model by representing the interference as 16 compensation functions, with each function composed of a fixed compensation coefficient and a base function [21]. Leach subsequently proposed the least squares method as a means to solve compensation model parameters while analyzing complex collinearity in the process [22]. Meanwhile, Gebre-Egziabher introduced a two-step calibration method that uses the least squares method to estimate intermediate variables and subsequently arrives at a solution for the model using those values [23]. However, Gebre-Egziabher’s calibration method falls short in addressing ill-condition problems [24], since its core algorithm still relies on the least squares method. As computational and mathematical solution methods continue to develop, more algorithms such as recursive least squares (RLS), truncated singular value decomposition, ridge regression, and neural networks have emerged as possible solutions [25,26,27,28]. Although these algorithms can improve computing accuracy, they cannot effectively improve the interference compensation accuracy as they are all based on the conventional Tolles–Lawson model. Currently, scholars studying geomagnetic compensation are increasingly focusing on the electromechanical interference of carriers, which the Tolles–Lawson model does not account for. In reference [29], the authors proposed a compensation method for the magnetic interference caused by the relative motion component and onboard current. In reference [30], the dynamic interference of various airborne equipment is considered, which is useful for reducing interference magnetic field. By studying the sources of these additional magnetic interferences and developing new compensation models, one can effectively improve the precision of magnetic interference compensation. However, available models have not taken the significant interference caused by variations in the velocity of UUV into account, and the compensation for this interference has yet to be addressed.
In this paper, a modified compensation method for UUV magnetic interference that takes velocity into account is proposed. We introduce a velocity-related term to improve the Tolles–Lawson model, which greatly expands the model’s application to UUV sailing at various velocities. We estimate the compensation parameters in two steps. First, we use the ridge estimation algorithm, which has certain advantages in fitting ill-conditioned data, to calculate the initial compensation parameters for several sets of constant velocities. Then, we use linear fitting to determine the functional relation between the initial compensation parameter and UUV velocity. By substituting the functional relation into the initial model, we establish a velocity-related compensation model. Finally, we verify the proposed compensation model’s practicability with a series of field experiments. The contributions made in this paper are as follows:
  • The impact of UUV velocity variation on the accuracy of underwater magnetic detection has been considered. A velocity-related term is introduced to improve the compensation accuracy for carrier interference magnetic field.
  • A two-step approach with ridge regression and linear fitting is employed to determine the compensation parameters, which has good operability. The ridge regression algorithm is applied to address the ill-conditioning of measured date to enhance the reliability of parameter estimation.
  • The proposed compensation method achieves high-precision compensation in field experiments, and it is feasible and effective for underwater magnetic detection.

2. Velocity-Related Compensation Method for UUV Magnetic Interference

2.1. Interference Magnetic Field Model of UUVs

In the measurement system, a triaxial magnetometer is installed on the UUV to capture magnetic signals in real time. However, during practical measurement, these signals often contain interference from ferromagnetic materials present on the UUV. This interference can lead to distortion of underwater magnetic detection. Specifically, the interference typically includes the permanent magnetic field, induced magnetic field, and eddy current magnetic field. The eddy current magnetic field is mainly caused by conductive materials constantly cutting through the geomagnetic field. Since it consisted of high-frequency components, low-pass filtering is commonly employed for reducing its interference during data processing [31].
Moreover, the presence of hard magnetic materials on the UUV can generate a permanent magnetic field over a period of long-term exposure to a strong magnetic field. This field is continuous and typically does not weaken under most circumstances. Its expression is given as:
H p = H p x H p y H p z T ,
where H p x , H p y , H p z denote the projections of the permanent magnetic field on the x, y, and z-axes of the triaxial magnetometer, respectively. Since the carrier is rigidly connected to the underwater magnetic vector measurement system, these projections are treated as constants.
The induced magnetic field is generated by the magnetization of soft magnetic materials due to an external magnetic field. A 3 × 3 matrix with nine parameters is utilized to represent the magnetization and interaction of each axis in the external magnetic field, which can be expressed as:
K = k x x k x y k x z k y x k y y k y z k z x k z y k z z
Furthermore, based on Euler’s theorem [32], the geomagnetic field in the geographic coordinate system can be converted into the measurement field in the magnetometer coordinate system using the rotation matrix R m as shown below:
R m = 1 0 0 0 cos α sin α 0 sin α cos α cos β 0 sin β 0 1 0 sin β 0 cos β cos γ sin γ 0 sin γ cos γ 0 0 0 1 ,
where α , β , γ represent three Euler angles between the geographic coordinate system and the magnetometer coordinate system.
Hence, according to the Tolles–Lawson model [20], the interference magnetic field of UUVs can be expressed as follows:
H e x = K R m H e + H p ,
where H e = H e x H e y H e z T denotes the value of the geomagnetic field in the geographic coordinate system.
The measured magnetic field can be initially modeled as follows:
H m = R m H e + H e x = K + E R m H e + H p ,
where E denotes the identity matrix.

2.2. Velocity-Related Compensation Model

2.2.1. The Mathematical Model

Based on the analysis in Section 2.1, we can derive the expression of the interference compensation model as follows:
H e = R m 1 M k 1 ( H m H p ) M k = ( K + E ) = m 11 m 12 m 13 m 21 m 22 m 23 m 31 m 32 m 33
The internal elements of M k are considered constant, as the magnetization parameters of soft magnetic materials on UUV remain unchanged in all directions.
Apart from the permanent magnetic field and the induced magnetic field, electromechanical devices can introduce extra interference while the UUV is operational. During underwater magnetic detection, changes in velocity are necessary and inevitable. The control of UUV velocity is achieved through changes in motor speed and battery current, which must lead to changes in the magnetic field. Therefore, this velocity-related magnetic field interference cannot be ignored. To confirm this, we secured the UUV with a non-magnetic truss and controlled the motor to run at different speeds during a shore-based experiment. Through spectral analysis, we found that the generated magnetic field interference contains both AC and DC components, where the AC signal can be filtered out by a low-pass filter, while the DC signal reveals a linear relationship with the motor speed. The target signals we need to detect are often low-frequency signals, so it is not sufficient to consider only the permanent magnetic interference in H p .
To compensate for the interference resulting from velocity variations, we introduce a velocity-related term into the compensation model. Accordingly, we establish a linear function to describe the relationship between velocity and the interference it generates, and define the correlation coefficient as follows:
G = [ g x g y g z ] T
The velocity-related term is expressed as:
H p = V G + H p ,
where V is the scalar velocity of UUV, H p is related to the permanent magnetic field.
By substituting H p into Equation (6), we can derive the velocity-related compensation model as follows:
H e = R m 1 M k 1 ( H m H p ) = R m 1 M k 1 ( H m V G H p )

2.2.2. Compensation Principle

The parameter estimation procedure for total magnetic field compensation is mainly divided into two parts. First, we establish a parameter estimation method to solve the values of M k and H p at different velocities of the UUV. Second, we obtain the coefficients related to V in the linear model of H p through fitting.
By reorganizing Equation (9), we obtain the following equation:
H e 2 = H e T H e = ( H m H p ) T ( M k 1 ) T ( R m 1 ) T R m 1 M k 1 ( H m H p ) = ( H m H p ) T C ( H m H p ) C = ( M k 1 ) T M k 1 = c 11 c 12 c 13 c 12 c 22 c 23 c 13 c 23 c 33 ,
where ( R m 1 ) T R m 1 = E for the rotation matrix is the orthogonal matrix, the value of H e 2 can be acquired from the output of the magnetometer.
The above equation can be further expensed as follows:
H e 2 f ( C , H p ) = H P ,
where f ( C , H p ) = c 11 H p x 2 + c 22 H py 2 + c 33 H pz 2 + 2 c 12 H p x H py + 2 c 13 H p x H pz + 2 c 23 H py H pz , H = H m x 2 , H m y 2 , H m z 2 , H m x H m y , H m x H m z , H m y H m z , H m x , H m y , H m z , P = p 1 p 2 p 3 p 4 p 5 p 6 p 7 p 8 p 9 = c 11 c 22 c 33 2 c 12 2 c 13 2 c 23 2 c 11 H m x c 12 H m y c 13 H m z 2 c 22 H m y c 12 H m x c 23 H m z 2 c 33 H m z c 13 H m x c 23 H m y .
H contains only magnetic measured data and P is the parameter vector that needs to be estimated. Since f ( C , H p ) < < H e 2 , it can be ignored during parameter estimation. Thus, the equation can be simplified to a linear form as H e 2 = H P . However, severe ill-conditioning often exists in the measured matrix due to the non-ergodicity of the UUV’s attitude during actual measurement. To address this issue, we use ridge regression [33] to estimate the parameter vector P based on the linear equation:
P ^ = H T H + α E 1 H T H 0 2 ,
where α is the ridge parameter used to suppress the ill-conditioning of the above equation.
As the parameters in P are estimated, the value of f ( C , H p ) can be calculated. Accordingly, iterators can be created to update P in order to reduce the estimation error:
P ^ i = H T H + α E 1 H T [ H 0 2 f i 1 ( C , H p ) ]
As mentioned earlier, compensation parameters C and H p can be determined when the UUV sails at a certain velocity. Using the estimates of H p at various velocities of the UUV as shown in Equation (14), the optimal values of G and H P can be obtained through linear fitting.
H p 1 = V 1 G + H p H p 2 = V 2 G + H p H p 3 = V 3 G + H p
Finally, the compensation for total magnetic field distortion can be achieved through the velocity-related compensation model expressed in Equation (10).
The system flowchart of the velocity-related compensation method is presented in Figure 1.

3. Field Experiment

3.1. Experimental Setup

The field experiment was conducted in an experimental reservoir in Hebei Province, which is undeveloped and far away from artificial construction. The intensity of the local geomagnetic field approximates the IGFR 2012 values monitored with a proton magnetometer. A 200 m caliber UUV with a non-magnetic aluminum alloy shell was used for the experiment, and its main construction is shown in Figure 2. A triaxial fluxgate magnetometer is used to detect the magnetic signal. To minimize stray magnetic fields generated by the motor system installed in the tail and interference from the battery pack, the magnetometer was mounted in a plastic watertight cabin extending from the head of the UUV. Additionally, a Doppler Velocity Log (DVL) and an Inertial Navigation System (INS) were served as the velocity sensor while functioning as the attitude sensor within the UUV, respectively.
The triaxial fluxgate magnetometer used in the experiment can provide highly precision measurements of triaxial static and alternating magnetic fields. Its main performance indicators are listed in Table 1. Additionally, the magnetometer is sealed using poly-ether-ether-ketone (PEEK) with excellent pressure resistance and waterproofness, making it ideal for underwater measurements.
To reduce errors in the measured data, the UUV navigated at various headings, pitches, and roll angles whenever possible. Thus, it was controlled to sail along a typical circular trajectory submerged to a depth of 1.5 m underwater. The diameter of this trajectory was minimized to ensure a nearly constant intensity of the geomagnetic field. The underwater geomagnetic total field strength was measured at the center of the trajectory with a proton magnetometer and was found to be 53,975 nT, which was used as the true value for estimating compensation parameters. Before the field experiment, the magnetometer was calibrated with the method introduced in [34]. After calibration, measurement error was reduced from 43 nT to less than 1.2 nT. Moreover, the UUV sailed along a straight line to ensure the proper operation of the measurement system in the pre-experiment. The entire trajectory of the UUV is illustrated in Figure 3a and the attitude angles measured by the INS are shown in Figure 3b–d.

3.2. Estimation and Compensation

In accordance with the typical cruising velocity of UUVs for marine exploration [35], we set the UUV velocity at around 3 knots in the experiment. Initially, the velocity of UUV was set as 2.1 knots, which corresponded to 20% of its maximum. The actual velocity of UUV was monitored by DVL in real time, which fluctuated by around 2.1 knots as shown I Figure 4c. The measured data of magnetic field components and the calculated total magnetic field are presented in Figure 4a,b. As indicated in the figure, the maximum fluctuation of the total magnetic field calculated based on the true value is nearly 300 nT.
The measured data were subsequently used to calculate compensation parameters C and H p according to Equation (11). For total magnetic field compensation, the exact value of M k is not necessary. To simplify the calculation, we assume that M k is a symmetric matrix, allowing the values of its elements to be calculated easily by matrix decomposition. The resulting solutions are as follows (Table 2):
Finally, these parameters were used to compensate for the interference magnetic field. As demonstrated in Equation (10), the different attitudes of the UUV have no effect on the compensation of the total magnetic field.
Figure 5 compares the total magnetic field before and after compensation. To make the results more intuitive, the true value measured by the proton magnetometer was taken as a reference and presented as a green line. Two primary statistical indicators, maximum fluctuation and root mean square error (RMSE), were employed as evaluation indices of the compensation accuracy in Table 3. After compensation, the fluctuation of the total magnetic field decreased from 315.05 nT to 5.30 nT (1.7% of the value before compensation), while the RMSE dropped from 231.30 nT to 1.61 nT (0.7% of the value before compensation).
In the same way, the compensation parameters when UUV sail at 3.1 knots (30%) and 4.1 knots (40%) can be obtained as follows (Table 4):
The compensation results are presented in Figure 6. After compensation, the maximum fluctuation of the magnetic total field decreased from 298.64 nT to 6.53 nT (2.2%) at 3.1 knots and from 313.45 nT to 9.64 nT (3.1%) at 4.1 knots. The RMSE of the magnetic total field also decreased from 231.16 nT to 2.04 nT (0.9%) at 3.1 knots and from 235.43 nT to 2.84 nT (1.2%) at 4.1 knots (Table 5). These results suggest that the compensation method has reliable performance at a certain velocity.
The estimation results for different velocities reveal that H p exhibits significant variation, while M k remains nearly constant. Minor variations are believed to be caused by the non-uniformity of the magnetic field along the circular trajectory, with M k primarily determined by inherent characteristics of the measurement system based on UUV, such as the magnetization coefficient of the soft magnetic material. Figure 7 provides the fitting results by the linear model based on the H p values at the velocities of 2.1, 3.1, and 4.1 knots.
The G and H p in Equation (8) can be calculated as:
G = 2.30 6.40 24.57 T
H p = [ 341.99 71.27 28.78 ] T
The reliability of the linearity assumption is supported by a goodness-of-fit level that approaches 0.99. Based on the linear relationship between H p and V , it can be concluded that changes in UUV velocity do interfere with the interference magnetic field.

3.3. Verification of Velocity-Related Compensation Method

To assess the effectiveness of the velocity-related compensation method, the UUV was operated to sail at 5.1 knots along the same circular trajectory in the reservoir.
In comparison, Figure 8 presents the velocity-unrelated compensation result. In this case, the parameters estimated from the measured data at 2.1 knots were directly employed to compensate for the interference magnetic field since the velocity-related term was not considered.
After compensation, both the fluctuation and RMSE remained at high levels (Table 6). This indicates that this compensation method is not applicable for compensating interference when the velocity of the UUV varies.
Afterwards, the velocity-related model was employed for compensation. Given that the velocity of the UUV is known and that G and H p were estimated in Section 3.2, H p at 5.1 knots could be calculated using Equation (8). Additionally, M k was determined by averaging the estimates from the three sets of data mentioned above. Using H p and M k , the velocity-related compensation model was applied to achieve interference compensation at 5.1 knots. The resulting plot is presented in Figure 9.
Table 7 shows that after compensation, the fluctuation was reduced to 6.83 nT (2%), and the RMSE was reduced to 2.23 nT (0.9%).
In Figure 10, we present a comparison between the velocity-related compensation results and the velocity-unrelated compensation results. When the sailing velocity of UUV changes, the magnetic interference compensation ability of the initial method greatly diminishes, which highlights the necessity of accounting for velocity-related magnetic disturbances in actual UUV operations. Following the velocity-related compensation, the curve of the calibrated magnetic field closely approximates the reference curve with small fluctuation and RMSE, undoubtedly offering better compensation performance.
In order to further prove the reproducibility of the velocity-related compensation method proposed in this paper, the verification experiment when UUV sails at 3.5 knots is carried out in the same reservoir. The statistical indicators of the total magnetic field before and after compensation are presented in Table 8. It can be found that the compensation parameters estimated based on the velocity-related model are applicable to the measurement data at 3.5 knots. Additionally, Figure 11 shows the statistical indicators at different velocities, providing a validation of the consistency of the results and the effectiveness of the compensation method at varying sailing velocities.
The initial model primarily compensates for the permanent and induced interference magnetic fields, without considering the interference caused by carrier electromechanical equipment. It inevitably leads to inaccuracies in the compensation since the electromechanical interference is unavoidable and exerts a significant influence when the UUV measurement system detects abnormal magnetic field in concerned water areas. In practice, the velocity of UUVs is determined by the operating state of the motor and the battery. As UUV velocity changes, the disturbing magnetic field generated by the motor and battery changes. Hence, the introduction of the velocity-related term in the compensation model becomes necessary. The experiment results prove that the velocity-related compensation model can improve compensation accuracy.

4. Conclusions

In this paper, we present an improved method for compensating carrier interference magnetic fields. First, we establish an initial compensation model based on the Tolles–Lawson model. Next, we estimate the compensation parameters for the UUV at several velocities using the ridge estimation algorithm. We then develop the velocity-related term based on the relationship between the compensation parameters and the velocity and introduce it into the initial magnetic compensation model, achieving the establishment of the velocity-related model. Finally, we verify the compensation method through field experiments. The results show that the initial model based on the Tolles–Lawson model performs well only when the UUV sails at a constant velocity. By contrast, the velocity-related model has good applicability even when the velocity of the UUV varies, with the fluctuation reduced from 334.90 nT to 6.83 nT and the RMSE reduced from 246.90 nT to 2.23 nT after compensation. Additionally, the reproducibility of the compensation method has been proved. These experiment results demonstrate that the velocity-related compensation model can effectively reduce the influence of the variation of UUV velocity on the accuracy of interference magnetic field compensation. Thus, it presents a highly adaptable method for underwater magnetic detection.

Author Contributions

Conceptualization, H.L. and Q.D.; methodology, H.L.; validation, Q.D., J.H. and Z.D.; resources, Q.Z.; writing—original draft preparation, H.L.; writing—review and editing, Q.D.; project administration, M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, W.; Yong, C.; Xiao, S.; Yichao, Z.; Zhihong, Y. A Novel Underwater Target Detection Method Based on Low-Frequency Magnetic Signal Generated by Portable Transmitter. IEEE Sens. J. 2023, 23, 1. [Google Scholar]
  2. Zhao, Y.; Zhang, J.H.; Li, J.H.; Liu, S.; Miao, P.X.; Shi, Y.C.; Zhao, E.M. A brief review of magnetic anomaly detection. Meas. Sci. Technol. 2021, 32, 042002. [Google Scholar] [CrossRef]
  3. Qin, J.; Li, M.; Li, D.; Zhong, J.; Yang, K. A survey on visual navigation and positioning for autonomous UUVs. Remote Sens. 2022, 14, 3794. [Google Scholar] [CrossRef]
  4. Xu, L.; Huang, J.; Zhang, H.; Liao, B. Direction of Arrival Estimation of Acoustic Sources with Unmanned Underwater Vehicle Swarm via Matrix Completion. Remote Sens. 2022, 14, 3790. [Google Scholar] [CrossRef]
  5. Gallimore, E.; Terrill, E.; Pietruszka, A.; Gee, J.; Nager, A.; Hess, R. Magnetic survey and autonomous target reacquisition with a scalar magnetometer on a small AUV. J. Field Robot. 2020, 37, 1246–1266. [Google Scholar] [CrossRef]
  6. Su, W. Magnetic Anomaly Data Detection of Local Marine Geomagnetic Field Model considering Robust Trend Surface Scientific Calculation Algorithm. Sci. Program. 2022, 2022, 4055976. [Google Scholar] [CrossRef]
  7. Gul, K.M.; Kaya, C.; Bektas, A. Design and Control of an Unmanned Underwater Vehicle. In Proceedings of the 2020 4th International Symposium on Multidisciplinary Studies and Innovative Technologies (ISMSIT), Istanbul, Turkey, 22–24 October 2020; pp. 1–9. [Google Scholar]
  8. Walker, C.R.; Stringfield, J.Q.; Wolbrecht, E.T. Measurement of the magnetic signature of a moving surface vessel with multiple magnetometer-equipped AUVs. Ocean. Eng. 2013, 64, 80–87. [Google Scholar] [CrossRef]
  9. Kasaya, T.; Nogi, Y.; Kitada, K. Advanced magnetic survey system and method for detailed magnetic field mapping near the sea bottom using an autonomous underwater vehicle. Explor. Geophys. 2023, 54, 205–216. [Google Scholar] [CrossRef]
  10. Li, S.; Cheng, D.; Wang, Y.; Zhao, J. Compensation Method for the Carrier Magnetic Interference of Underwater Magnetic Vector Measurement System. IEEE Sens. J. 2023, 23, 10694–10705. [Google Scholar] [CrossRef]
  11. Dautta, M.; Faruque, R.B.; Islam, R. A mathematical model of extremely low frequency ocean induced electromagnetic noise. In AIP Conference Proceedings; AIP Publishing LLC.: Melville, NY, USA, 2016; Volume 1754, p. 060010. [Google Scholar]
  12. Zhang, C.; Du, C.; Peng, X.; Han, Q.; Guo, H. An Aeromagnetic Compensation Method for Suppressing the Magnetic Interference Generated by Electric Current with Vector Magnetometer. Sensors 2022, 22, 6151. [Google Scholar] [CrossRef]
  13. Xia, T.; Cui, D.; Chu, Z.; Yu, X. Autonomous Heading Planning and Control Method of Unmanned Underwater Vehicles for Tunnel Detection. J. Mar. Sci. Eng. 2023, 11, 740. [Google Scholar] [CrossRef]
  14. Tao, Q.; Jiapin, C.; Zhenbo, L. The driving control system of micro brushless DC motor used in micro craft. Small Spec. Electr. Mach. 2007, 35, 48–50. [Google Scholar]
  15. Pang, H.; Wan, C.; Mou, S.; Wei, L.; Li, H.; Zhang, Q.; Pan, M.; Yang, D. Integrated calibration of strap-down geomagnetic vector measurement system. IEEE Sens. J. 2022, 22, 10476–10484. [Google Scholar] [CrossRef]
  16. Gao, Q.; Zhao, J.; Cheng, D. A compensation method for the carrier interference of a three-component magnetic measurement system using a Cuckoo search algorithm. Meas. Sci. Technol. 2018, 29, 085003. [Google Scholar] [CrossRef]
  17. Wang, M.; Lin, J.; Yue, L. Compensation for mobile carrier magnetic interference in a SQUID-based full-tensor magnetic gradiometer using the flower pollination algorithm. Meas. Sci. Technol. 2021, 32, 085010. [Google Scholar] [CrossRef]
  18. Wan, C.; Pan, M.; Zhang, Q.; Pang, H.; Zhu, X.; Pan, L.; Sun, X.; Wu, F. Improved component compensation for geomagnetic field vector measurement using Lagrange multiplier method. IET Sci. Meas. Technol. 2017, 11, 1094–1098. [Google Scholar] [CrossRef]
  19. Tolles, W.E. Compensation of Aircraft Magnetic Fields. U.S. Patent 2,692,970, 26 October 1954. [Google Scholar]
  20. Leliak, P. Identification and evaluation of magnetic-field sources of magnetic airborne detector equipped aircraft. IRE Trans. Aerosp. Navig. Electron. 1961, 8, 95–105. [Google Scholar] [CrossRef]
  21. Leach, B.W. Aeromagnetic compensation as a linear regression problem. In Information Linkage between Applied Mathematics and Industry; Academic Press: Cambridge, MA, USA, 1980; pp. 139–161. [Google Scholar]
  22. Gebre-Egziabher, D.; Elkaim, G.H.; David, P.J.; Parkinson, B.W. Calibration of strapdown magnetometers in magnetic field domain. J. Aerosp. Eng. 2006, 19, 87–102. [Google Scholar] [CrossRef]
  23. Gu, Y.; Gui, Q.; Zhang, X.; Han, S. Iterative solution of regularization to ill-conditioned problems in geodesy and geophysics. J. Surv. Mapp. 2020, 2, 59–65. [Google Scholar]
  24. Dou, Z.; Ren, K.; Han, Q. A novel real-time aeromagnetic compensation method based on RLSQ. In Proceedings of the 2014 Tenth International Conference on Intelligent Information Hiding and Multimedia Signal Processing, Kitakyushu, Japan, 27–29 August 2014; pp. 243–246. [Google Scholar]
  25. Su, Z.; Jiao, J.; Zhou, S.; Yu, P.; Zhao, X. Aeromagnetic compensation method based on ridge regression algorithm. Glob. Geol. 2022, 25, 41–48. [Google Scholar]
  26. Gu, B.; Li, Q.; Liu, H. Aeromagnetic compensation based on truncated singular value decomposition with an improved parameter-choice algorithm. In Proceedings of the 2013 6th International Congress on Image and Signal Processing (CISP), Hangzhou, China, 16–18 December 2013; Volume 3, pp. 1545–1551. [Google Scholar]
  27. Mu, Y.; Zhang, X.; Xie, W.; Zheng, Y. Automatic detection of near-surface targets for unmanned aerial vehicle (UAV) magnetic survey. Remote Sens. 2020, 12, 452. [Google Scholar] [CrossRef]
  28. Du, C.; Wang, H.; Wang, H. Extended aeromagnetic compensation modelling including non-manoeuvring interferences. IET Sci. Meas. Technol. 2019, 13, 1033–1039. [Google Scholar] [CrossRef]
  29. Noriega, G.; Marszalkowski, A. Adaptive techniques and other recent developments in aeromagnetic compensation. First Break. 2017, 35, 31–38. [Google Scholar] [CrossRef]
  30. Li, Y.; Han, Q.; Peng, X.; Li, Q.; Tong, X. A Magnetic Interference Compensation Method for Airborne Electronic Equipment without Current Sensors. Remote Sens. 2022, 14, 4151. [Google Scholar] [CrossRef]
  31. Formulae Generales Pro Translatione Quacunque Corporum Rigidorum. Available online: https://math.dartmouth.edu/~euler/docs/originals/E478.pdf (accessed on 2 November 2019).
  32. Ehsanes Saleh, A.K.; Arashi, M.; Tabatabaey, S.M.M. Ridge Regression. In Statistical Inference for Models with Multivariate t-Distributed Errors; Wiley: New York, NY, USA, 2014; pp. 133–170. [Google Scholar]
  33. Li, L.; Liu, W. Compensation of Optical Pump Magnetometer Using the Improved Mayfly Optimization Algorithm. J. Mar. Sci. Eng. 2022, 10, 1982. [Google Scholar] [CrossRef]
  34. Qi, Z.; Hongfeng, P.; Chengbiao, W. Magnetic interference compensation method for geomagnetic field vector measurement. Measurement 2016, 91, 628–633. [Google Scholar]
  35. Thanh, H.; Phung, S.L.; Chapple, P.B.; Bouzerdoum, A.; Ritz, C.H.; Tran, L.C. Deep Gabor Neural Network for Automatic Detection of Mine-Like Objects in Sonar Imagery. IEEE Access 2020, 8, 94126–94139. [Google Scholar] [CrossRef]
Figure 1. System structure diagram of the velocity-related compensation method.
Figure 1. System structure diagram of the velocity-related compensation method.
Remotesensing 15 04164 g001
Figure 2. The UUV magnetic field measurement system: (a) diagram of the UUV; and (b) schematic of the measurement system.
Figure 2. The UUV magnetic field measurement system: (a) diagram of the UUV; and (b) schematic of the measurement system.
Remotesensing 15 04164 g002
Figure 3. UUV trajectory and attitude: (a) the test trajectory; (b) the heading angle when UUV is sailing along the trajectory; (c) the pitch angle; and (d) the roll angle.
Figure 3. UUV trajectory and attitude: (a) the test trajectory; (b) the heading angle when UUV is sailing along the trajectory; (c) the pitch angle; and (d) the roll angle.
Remotesensing 15 04164 g003
Figure 4. The measurement results at 2.1 knots: (a) magnetic-field components; (b) total magnetic field; and (c) actual velocity of UUV.
Figure 4. The measurement results at 2.1 knots: (a) magnetic-field components; (b) total magnetic field; and (c) actual velocity of UUV.
Remotesensing 15 04164 g004
Figure 5. Data comparison before and after compensation at 2.1 knots.
Figure 5. Data comparison before and after compensation at 2.1 knots.
Remotesensing 15 04164 g005
Figure 6. Data comparison before and after compensation: (a) at 3.1 knots and (b) at 4.1 knots.
Figure 6. Data comparison before and after compensation: (a) at 3.1 knots and (b) at 4.1 knots.
Remotesensing 15 04164 g006
Figure 7. The linear fitting of H p and V .
Figure 7. The linear fitting of H p and V .
Remotesensing 15 04164 g007
Figure 8. Data comparison before and after velocity-unrelated compensation at 5.1 knots.
Figure 8. Data comparison before and after velocity-unrelated compensation at 5.1 knots.
Remotesensing 15 04164 g008
Figure 9. Data comparison before and after velocity-related compensation at 5.1 knots.
Figure 9. Data comparison before and after velocity-related compensation at 5.1 knots.
Remotesensing 15 04164 g009
Figure 10. Contrast of velocity-related compensation and velocity-unrelated compensation at 5.1 knots.
Figure 10. Contrast of velocity-related compensation and velocity-unrelated compensation at 5.1 knots.
Remotesensing 15 04164 g010
Figure 11. Indicators comparisons before and after velocity-related compensation at 2.1, 3.1, 3.5, 4.1, 5.1 knots: (a) Fluctuation and (b) RMSE.
Figure 11. Indicators comparisons before and after velocity-related compensation at 2.1, 3.1, 3.5, 4.1, 5.1 knots: (a) Fluctuation and (b) RMSE.
Remotesensing 15 04164 g011
Table 1. Main performance indicators of the triaxial fluxgate magnetometer.
Table 1. Main performance indicators of the triaxial fluxgate magnetometer.
PerformanceParameter Value
Measurement range±100,000 nT
Non-orthogonality error<±0.1°
Linearity<0.0015%
Noise<10 pT RMS√Hz@1Hz
Maximum offset±5 nT at ± 100 μT
Bandwidth3 kHz
Frequency responseDC to 1 kHz
Table 2. Values of compensation parameters when UUV sails at 2.1 knots.
Table 2. Values of compensation parameters when UUV sails at 2.1 knots.
M k H p
1.00626 0.00006 0.00203 0.00006 0.99632 0.00057 0.00203 0.00057 0.99457 347.13 85.06 82.39
Table 3. Measuring accuracy of total magnetic field compensation at 2.1 knots.
Table 3. Measuring accuracy of total magnetic field compensation at 2.1 knots.
Fluctuation (nT)RMSE (nT)
Uncompensation315.05231.66
Compensation5.301.42
Table 4. Values of compensation parameters when UUV sail at 3.1 and 4.1 knots.
Table 4. Values of compensation parameters when UUV sail at 3.1 and 4.1 knots.
Velocity M k H p
3.1 knots 1.00633 0.00007 0.00199 0.00007 0.99670 0.00056 0.00199 0.00056 0.99460 348.98 91.80 106.03
4.1 knots 1.00651 0.00006 0.00211 0.00006 0.99658 0.00056 0.00211 0.00056 0.99457 351.68 97.71 130.93
Table 5. Measuring accuracy of total magnetic field compensation at 3.1 and 4.1 knots.
Table 5. Measuring accuracy of total magnetic field compensation at 3.1 and 4.1 knots.
Fluctuation (nT)RMSE (nT)
3.1 knotsUncompensation298.64231.16
Compensation6.532.04
4.1 knotsUncompensation313.45235.43
Compensation9.642.84
Table 6. Measuring accuracy of velocity-unrelated compensation at 5.1 knots.
Table 6. Measuring accuracy of velocity-unrelated compensation at 5.1 knots.
Fluctuation (nT)RMSE (nT)
Uncompensation334.90246.90
Compensation226.52141.15
Table 7. Measuring accuracy of velocity-related compensation at 5.1 knots.
Table 7. Measuring accuracy of velocity-related compensation at 5.1 knots.
Fluctuation (nT)RMSE (nT)
Uncompensation334.90246.90
Compensation6.832.23
Table 8. Measuring accuracy of velocity-related compensation at 3.5 knots.
Table 8. Measuring accuracy of velocity-related compensation at 3.5 knots.
Fluctuation (nT)RMSE (nT)
Uncompensation310.9234.3
Compensation6.722.16
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

Luo, H.; Pan, M.; Du, Q.; Zhang, Q.; Hu, J.; Ding, Z. Velocity-Related Magnetic Interference Compensation of Unmanned Underwater Vehicle. Remote Sens. 2023, 15, 4164. https://doi.org/10.3390/rs15174164

AMA Style

Luo H, Pan M, Du Q, Zhang Q, Hu J, Ding Z. Velocity-Related Magnetic Interference Compensation of Unmanned Underwater Vehicle. Remote Sensing. 2023; 15(17):4164. https://doi.org/10.3390/rs15174164

Chicago/Turabian Style

Luo, Huihui, Mengchun Pan, Qingfa Du, Qi Zhang, Jiafei Hu, and Zengquan Ding. 2023. "Velocity-Related Magnetic Interference Compensation of Unmanned Underwater Vehicle" Remote Sensing 15, no. 17: 4164. https://doi.org/10.3390/rs15174164

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

Article Metrics

Back to TopTop