Simulation Analysis of Magnetic Gradient Full-Tensor Measurement System

The magnetic gradient full-tensor measurement system is diverse, and the magnetometer array structure is complex. Aimed at the problem, seven magnetic gradient full-tensor measurement system models are studied in detail. The full-tensor measurement theories of the tensor measurement arrays are analyzed. Under the same baseline distance, the magnetic dipole model is used to simulate the measurement system. Based on diﬀerent measurement systems, the paper quantitatively compares and analyzes the error of the structure. A more optimized magnetic gradient full-tensor measurement system is suggested. The simulation results show that the measurement accuracy of the planar measurement system is slightly higher than that of the stereo measurement system. Among them, the cross-shaped and square measurement systems have relatively smaller structural errors and higher measurement accuracy.


Introduction
e ferromagnetic object magnetized by the geomagnetic field will produce a magnetic field, which will affect the constant distribution of the geomagnetic field. is phenomenon that causes the abnormality of the geomagnetic field is called magnetic anomaly [1]. Magnetic anomaly detection refers to the technology of obtaining several characteristics of the magnetic target with the process of acquiring and processing the magnetic anomaly information of a magnetic target. Magnetic anomaly detection technology has many advantages such as strong anti-interference ability, rich information, and good concealment performance [2][3][4]. It is widely used in underwater or underground target detection, aerial magnetic survey, navigation and positioning, mineral exploration, and other fields [5][6][7][8][9]. At present, it can be roughly divided into four phases: total magnetic field detection, total magnetic field gradient detection, the magnetic component field detection, and the magnetic gradient full-tensor field detection phase [10]. e phase of magnetic gradient full-tensor detection is the latest development stage of magnetic anomaly detection. It has unique advantages such as being less affected by system attitude changes and less environmental interference. It can get much richer target information and be not affected by the total magnetic background field basically. is technology has gradually developed into a research hotspot in the fields of Earth resource exploration, aerial exploration, and military reconnaissance, and its application prospects are very broad [11][12][13][14].
It is impossible to directly measure the magnetic gradient full tensor. In actual measurement, the magnetic gradient fulltensor measurement system is generally used to measure the components of the magnetic field vector field; and the difference method is used to obtain the full-tensor field of the magnetic gradient approximately. Presently, the magnetic gradient full-tensor measurement system mainly includes two types: superconducting magnetometer and fluxgate magnetometer [15,16]. e former has high sensitivity, good stability, and strong reliability, but the preparation process is complicated and costly. e latter uses several fluxgates to measure the magnetic vector component to obtain the spatial change rate, which is portable, economical, and practical.
At present, the magnetic gradient full-tensor measurement systems constructed by the fluxgate magnetometer can be divided into two types. One is the two-dimensional plane shape including cross, equilateral triangle, right triangle, and square; the other is three-dimensional, which includes regular tetrahedron, right-angle tetrahedron, and regular tetrahedron. Various systems can achieve the measurement purpose under specific environmental conditions. However, under the same measurement environment, the selection of the measurement system is still unclear, and there is very little analysis in the structure of different systems. But the configuration of the fluxgate magnetometer will have a great impact on the measurement of the magnetic gradient full tensor. It even directly affects the feasibility of measurement. erefore, it is necessary to analyze and compare the magnetic gradient full-tensor measurement system to seek a more optimized system. Aiming at this problem, this paper establishes the structural models of various magnetic gradient full-tensor measurement systems in detail. Using magnetic dipole models, simulation experiments are performed to analyze various magnetic gradient full-tensor measurement systems.

Magnetic Gradient Full-Tensor
Measurement Elements e magnetic field is a vector field, and the spatial rate of change in each direction of the component is the tensor of magnetic gradient [17], which includes 9 elements in total. e magnetic gradient full tensor (G) of the magnetic field vector ( B → ) can be shown as erefore, in equation (1) e essence of the magnetic gradient full tensor is the second-order partial differential of the magnetic scalar potential U and the first-order derivative of the magnetic field vector, which cannot be directly measured in actual operation. Generally, the ratio of the magnetic field to the distance is used to approximately obtain each component of the magnetic gradient full tensor. e difference calculation method is as follows: In equation (6), ΔB i is the difference between the magnetic field components measured by two adjacent vector fluxgates, and Δd j is the distance between the two vector fluxgates in the direction of j, which is defined as the baseline of the magnetic gradient full-tensor measurement system distance. 2 Mathematical Problems in Engineering

Cross Measurement
System. e model of the cross-shaped magnetic gradient full-tensor measurement system is shown in Figure 1.
e system is composed of four fluxgate magnetometers. e center point is taken as the axis center (O) of the system coordinate axis. e fluxgates are marked as 1, 2, 3, and 4, respectively. e centers of No. 1 and No. 3 fluxgates are parallel to the x-axis, and their baseline distance is b. No. 2 and No. 4 fluxgates are parallel to the y-axis, and their baseline distance is also b [18]. e three-component magnetic fields measured by the four fluxgates are as follows: In equation (7), B i (i � 1, 2, 3, 4) is the magnetic field vector of the fluxgate No. i, and B ij (i � 1, 2, 3, 4; j � x, y, z) is the component of the fluxgate No. i in the j direction. e expression of the magnetic gradient full tensor at the axis O of the system is as follows:

Equilateral Triangle Measurement
System. e triangular magnetic gradient full-tensor measurement system is composed of three fluxgate magnetometers. No. 1, No. 2, and No. 3 fluxgates are, respectively, located at the apex of the triangle. e side length of the equilateral triangle is the baseline distance (d) of the system. e geometric center of the triangle is taken as the axis center (O) of the system coordinate axis. e system model is shown in Figure 2.
Similar to the cross-shaped measurement system, according to the three components of the magnetic field vector measured by the No. 1, No. 2, and No. 3 fluxgates, the full tensor of the magnetic gradient of the axis can be obtained:

Right Triangle Measurement
System. e model of the right triangle magnetic gradient full-tensor measurement system is shown in Figure 3. e difference from the equilateral triangle measurement system is that the connection between the centers of the three fluxgate magnetometers is isosceles right triangle. e position where the No. 3 fluxgate is located is the axis center (O) of the system coordinate axis, and the right-angle side length d is the

Regular Tetrahedron Shape Measurement
System. e regular tetrahedral magnetic gradient full-tensor measurement system is built by 4 fluxgate magnetometers. e connection at the center of the fluxgate forms a regular tetrahedral structure. e structure model is shown in Figure 5. e centers of the No. 1, No. 2, and No. 3 fluxgates form the bottom surface of a regular tetrahedron, and the center of the bottom triangle is taken as the axis center (O) of the coordinate axis of the system. At the same time, the side length of the regular tetrahedron is the baseline distance (d) of the system. e magnetic field of the center O is shown as

Right-Angle Tetrahedron Shape Measurement System.
e model of the right-angled tetrahedral magnetic gradient full-tensor measurement system is shown in Figure 6 Figure 7, the regular hexahedron shape measurement system model is composed of 8 fluxgate magnetometers centered at the apex of the regular hexahedron. e side length of the regular hexahedron is the baseline distance d, and the center of the regular hexahedron is the coordinate axis center (O) of the system. e full-tensor expression of the magnetic gradient at the axis is

Hexahedron Shape Measurement System. As shown in
e structure comparison of seven kinds of full-tensor measurement systems is shown in Table 1.
It can be seen from Table 1 that the measurement system can be divided into two types: two-dimensional plane and three-dimensional stereo according to the number of dimensions, which also determines the type of observation point. Among all measurement systems, the number of magnetometers used in triangle is the smallest, which is 3; the number of magnetometers used in the regular hexahedron is the largest, which is 8. e threedimensional shape measurement system can measure 9 independent components individually, while the twodimensional plane can measure 6 independent components.

Simulation Experiment of Magnetic Gradient
Full-Tensor Measurement System 4.1. Magnetic Dipole Model. When the distance between the magnetic target point and the observation point is more than 2.5 times the diameter or length of the magnetic target point, the magnetic field generated by the magnetic target point can be equivalent to a magnetic dipole [19][20][21]. In the process of measuring the magnetic gradient full tensor, the magnetic gradient full-tensor measurement system is often far away from the magnetic target. According to actual experience, a magnetic dipole model can be used to replace the magnetic target for modeling. As shown in Figure 8, the magnetic dipole model is established. In the Cartesian coordinate system, the position of the magnetic target is (x 0 , y 0 , z 0 ), and the position of the observation point (x, y, z) is the axis of the magnetic gradient full-tensor measurement system. e projection components at the observation point of the magnetic anomaly field generated are represented by B x , B y , and B z , which are the three-component magnetic anomaly.
e magnetic anomaly vector field generated by the magnetic target at the observation point is    Mathematical Problems in Engineering In equation (15) e observation point in Figure 8 is the center point O in Figure 1 to Figure 7. e magnetic field value of point O can be obtained according to the coordinate equation (16). In order to obtain the full tensor of the magnetic gradient at point O, it is necessary to determine the location based on the structural parameters of the measurement system. e magnetic field value of each point can be determined by using the magnetic dipole equation several times. Take Figure 1 as an example to calculate the coordinates of No. 1 to No. 4 fluxgates, as shown in the following equation: In equation (17), r 1 , r 2 , r 3 , r 4 are the position coordinates of No. 1 to No. 4, with point O as the reference point. e magnetic field value of each magnetometer can be obtained, and then the full tensor of the magnetic gradient can be gotten.

Simulation Experiment.
Assuming that the magnetic dipole is located in a vacuum, the position of the magnetic source is (3,4,5). e magnetic declination angle D � 15°, the magnetic inclination angle I � 30°, the magnetic moment modulus m � 200A · m 2 , m x � m cos I cos d, m y � m cos I sin d, m z � sin I, and the grid area is 20m × 60m. During the sampling process, the attitude of the tensor measurement system is unchanged, and it moves unidirectionally along the x-axis. e sampling interval is set to 1 m; the baseline distance is set to 0.3 m. Taking the cross-shaped magnetic gradient full-tensor measurement system as an example, the comparison between distribution value and theoretical value of tensor is shown in Figure 9. e full-tensor distribution value of the seven structures is shown in Figure 10. In order to highlight the impact of the structural error in measurement system on the tensor measurement, the output values of the three-axis fluxgate are ideal during the simulation.

Simulation Results and Analysis.
It can be seen from Figure 10 that each measurement system can measure the full tensor within the allowable error range, but the error between the tensor value and theoretical value obtained from different structures is different. Among them, the errors of cross, equilateral triangle, right triangle, and square are relatively small. In order to describe the measurement error of tensor more clearly, equation (18) is used to calculate the relative error (RE) of each point in the simulation. In Figure 11, the measurement errors of each point of the measurement system are compared and shown.
In equation (18), f(i) is the measured value of each point, and y(i) is the theoretical value of each point.
It can be seen from Figure 11(a) that the Bxx measured by the square magnetic gradient full-tensor measurement system has the highest accuracy. Within the range of ±30 nT/ m, the maximum measurement error does not exceed 0.02 nT/m. e accuracy of right-angle tetrahedron is the worst, exceeding 2 nT/m. As can be seen from Figure 11(b), the Bxy measured by the cross-shaped magnetic gradient full-tensor measurement system has the highest accuracy. Within the range of ±30 nT/m, the maximum measurement error is not about 0.03 nT/m. e measurement accuracy of the right-angled tetrahedron is the worst, reaching 1.7 nT/m. It can be seen in Figure 11(c) that the Bxz measured by the cross-shaped magnetic gradient full-tensor measurement system has the highest accuracy. In the range of ±30 nT/m, the maximum measurement error does not exceed 0.02 nT/ m. e measurement accuracy of the right-angled tetrahedron is the worst, exceeding 2 nT/m. It can be seen in Figure 11(d) that Byy measured by the cross-shaped magnetic gradient full-tensor measurement system has the highest accuracy. Within the range of ±30 nT/m, the maximum measurement error does not exceed 0.02 nT/m, and the measurement accuracy of the regular hexahedron is the worst, reaching 30 nT/m. It can be seen in Figure 11(e) that the Byz accuracy measured by the cross-shaped magnetic gradient full-tensor measurement system is the highest. Within the range of ±30 nT/m, the maximum measurement error does not exceed 0.01 nT/m, and the measurement accuracy of the regular tetrahedron is the worst, reaching 15 nT/m. From Figure 11, it is basically possible to analyze the pros and cons of the seven measurement systems for each component, and the order from the best to the worst is shown in Table 2.
In order to further test the overall measurement accuracy of each system, equation (19) is used to calculate the root mean square error (RMSE) [22] to quantitatively analyze the overall deviation of the measurement system. e calculated values are shown in Table 3.      Equilateral triangle shape Regular tetrahedron shape Right-angle tetrahedron shape Right triangle shape (e) Figure 11: Comparison of measured RE of seven magnetic gradient full-tensor systems.
In equation (19), f i is the measured value of each point, and y i is the theoretical value of each point.
It can be seen from Tables 2 and 3 that, under the same measurement conditions, the cross-shaped and square measurement systems have relatively high accuracy, while the right-angled triangle and right-angled tetrahedral measurement systems have relatively great errors. Among them, the overall measurement accuracy of the cross shape is the highest in the analyzed measurement accuracy.

Conclusion
In this paper, seven magnetic gradient full-tensor measurement system models are studied in detail, and the fulltensor measurement theory of the tensor measurement array is analyzed. e magnetic dipole model is used to simulate the measurement system. Different measurement systems are quantitatively compared and analyzed. e analysis shows that, under the same measurement conditions, the cross-shaped and square measurement systems have smaller structural errors and higher measurement accuracy. In addition, the cross-shaped measurement system has the best overall performance. Of course, in engineering practice, the planar magnetic gradient full-tensor measurement system is easy to install, and the center point of the structure is easier to find. e three-dimensional magnetic gradient full-tensor measurement system can obtain many independent components and can directly obtain the vertical component. e magnetic gradient full-tensor measurement system should be selected to use according to actual needs. e research results in this paper provide a theoretical reference for the construction and application development of the magnetic gradient full-tensor system.

Data Availability
e data that support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.