# Simulation Analysis of Magnetic Gradient Full-Tensor Measurement System

Authors: Lei Xu; Ning Zhang; Liqing Fang; Huadong Chen; Pengfei Lin; Chunsheng Lin
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

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 different 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.


---

## Body

## 1. Introduction

The ferromagnetic object magnetized by the geomagnetic field will produce a magnetic field, which will affect the constant distribution of the geomagnetic field. This 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–4]. It is widely used in underwater or underground target detection, aerial magnetic survey, navigation and positioning, mineral exploration, and other fields [5–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]. The 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. This 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–14].It is impossible to directly measure the magnetic gradient full tensor. In actual measurement, the magnetic gradient full-tensor 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]. The former has high sensitivity, good stability, and strong reliability, but the preparation process is complicated and costly. The 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. Therefore, 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.


## 2. Magnetic Gradient Full-Tensor Measurement Elements

The 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. The magnetic gradient full tensor (G) of the magnetic field vector (B⟶) can be shown as(1)G=∂2U∂x2∂2U∂x∂y∂2U∂x∂z∂2U∂y∂x∂2U∂y2∂2U∂y∂z∂2U∂z∂x∂2U∂z∂y∂2U∂z2,=∂Bx∂x∂Bx∂y∂Bx∂z∂By∂x∂By∂y∂By∂z∂Bz∂x∂Bz∂y∂Bz∂z,=BxxBxyBxzByxByyByzBzxBzyBzz.In equation (1), U is the magnetic scalar potential; Bx, By, and Bz are the three components of the magnetic field at any point in space x,y,z; Biji,j=x,y,z is the component of the magnetic gradient full tensor. According to Maxwell’s equations, in a passive ambient static magnetic field, the divergence and curl of the magnetic field are both zero.(2)divB⟶=∇·B⟶=0,(3)rotB⟶=∇×B⟶,=ijk∂∂x∂∂y∂∂zBxByBz=0.According to equations (2) and (3),(4)∂Bx∂x+∂By∂y+∂Bz∂z=0,∂Bx∂y=∂By∂x,∂Bx∂z=∂Bz∂x,∂By∂z=∂Bz∂y.Therefore, in equation (1), only five elements, Bxx, Bxy, Bxz, Byy, and Byz, are required to obtain the magnetic gradient full tensor of the magnetic field.(5)G=BxxBxyBxzBxyByyByzBxzByz−Bxx−Byy.The 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. The difference calculation method is as follows:(6)Bij=∂Bi∂j≈ΔBiΔdj.In equation (6), ΔBi is the difference between the magnetic field components measured by two adjacent vector fluxgates, and Δdj 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.


## 3. Magnetic Gradient Full-Tensor Measurement System

### 3.1. Cross Measurement System

The model of the cross-shaped magnetic gradient full-tensor measurement system is shown in Figure1. The system is composed of four fluxgate magnetometers. The center point is taken as the axis center (O) of the system coordinate axis. The fluxgates are marked as 1, 2, 3, and 4, respectively. The 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].Figure 1
Cross-shaped magnetic gradient full-tensor measurement system.The three-component magnetic fields measured by the four fluxgates are as follows:(7)B1=B1xB1yB1z,B2=B2xB2yB2z,B3=B3xB3yB3z,B4=B4xB4yB4z.In equation (7), Bii=1,2,3,4 is the magnetic field vector of the fluxgate No. i, and Biji=1,2,3,4;j=x,y,z is the component of the fluxgate No. i in the j direction. The expression of the magnetic gradient full tensor at the axis O of the system is as follows:(8)G=1dB1x−B3xB1y−B3yB1z−B3zB2x−B4xB2y−B4yB2z−B4zB1z−B3zB2z−B4zB3x+B4y−B1x+B2y.


### 3.2. Equilateral Triangle Measurement System

The 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. The side length of the equilateral triangle is the baseline distance (d) of the system. The geometric center of the triangle is taken as the axis center (O) of the system coordinate axis. The system model is shown in Figure 2.Figure 2
Equilateral triangle magnetic gradient full-tensor measurement system.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:(9)G=13d3B1x−B3x3B1y−B3y3B1z−B3z2B2x−B1x+B3x2B2y−B1y+B3y2B2z−B1z+B3z3B1z−B3z2B2z−B1z+B3zB1y+B3y−2B2y−3B1x−B3x.


### 3.3. Right Triangle Measurement System

The model of the right triangle magnetic gradient full-tensor measurement system is shown in Figure3. The difference from the equilateral triangle measurement system is that the connection between the centers of the three fluxgate magnetometers is isosceles right triangle. The 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 baseline distance. The full tensor of the magnetic gradient of the axis O is shown as follows:(10)G=1dB1x−B3xB1y−B3yB1z−B3zB2x−B3xB2y−B3yB2z−B3zB1z−B3zB2z−B3zB3x+B3y−B1x+B2y.Figure 3
Right triangle magnetic gradient full-tensor measurement system.


### 3.4. Square Measurement System

Figure4 is a square magnetic gradient full-tensor measurement system model. The center points of the No. 1, No. 2, No. 3, and No. 4 fluxgates are connected to form a square. The center point of the square is taken as the axis center (O) of the coordinate axis of the system. The side length d of the square is the baseline distance, and the magnetic gradient full tensor of the axis O is calculated as(11)G=12dB1x+B4x−B2x+B3xB1y+B4y−B2y+B3yB1z+B4z−B2z+B3zB1x+B2x−B3x+B4xB1y+B2y−B3y+B4yB1z+B2z−B3z+B4zB1z+B4z−B2z+B3zB1z+B2z−B3z+B4zB2x+B3x+B3y+B4y−B1x+B1y+B2y+B4x.Figure 4
Square magnetic gradient full-tensor measurement system.


### 3.5. Regular Tetrahedron Shape Measurement System

The regular tetrahedral magnetic gradient full-tensor measurement system is built by 4 fluxgate magnetometers. The connection at the center of the fluxgate forms a regular tetrahedral structure. The structure model is shown in Figure5. The 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. The magnetic field of the center O is shown as(12)G=1d2B1x−B2x+B3x32B1y−B2y+B3y32B1z−B2z+B3z3B1x+B2x−2B3x3B1y+B2y−2B3y3B1z+B2z−2B3z33B4x−B1x+B2x+B3x63B4y−B1y+B2y+B3y63B4z−B1z+B2z+B3z6.Figure 5
The full-tensor measurement system of the regular tetrahedral magnetic gradient.


### 3.6. Right-Angle Tetrahedron Shape Measurement System

The model of the right-angled tetrahedral magnetic gradient full-tensor measurement system is shown in Figure6. The line connecting the centers of the fluxgate magnetometers No. 1, No. 2, and No. 3 forms the base of an isosceles right triangle, and the length of the right-angle side is the baseline distance (d). The distance between No. 4 and No. 1 fluxgates is also d, and the center of No. 1 fluxgate is the axis center (O) of the coordinate axis of the system. The full tensor of the magnetic gradient system is(13)G=1dB2x−B1xB2y−B1yB2z−B1zB3x−B1xB3y−B1yB3z−B1zB4x−B1xB4y−B1yB4z−B1z.Figure 6
Right-angled tetrahedral magnetic gradient full-tensor measurement system.


### 3.7. Hexahedron Shape Measurement System

As shown in Figure7, the regular hexahedron shape measurement system model is composed of 8 fluxgate magnetometers centered at the apex of the regular hexahedron. The 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. The full-tensor expression of the magnetic gradient at the axis is(14)G=14dB1x+B4x+B5x+B6x−B2x+B3x+B7x+B8xB1y+B4y+B5y+B6y−B2y+B3y+B7y+B8yB1z+B4z+B5z+B6z−B2z+B3z+B7z+B8zB3x+B4x+B6x+B7x−B1x+B2x+B5x+B8xB3y+B4y+B6y+B7y−B1y+B2y+B5y+B8yB3z+B4z+B6z+B7z−B1z+B2z+B5z+B8zB1z+B4z+B5z+B6z−B2z+B3z+B4z+B8zB3z+B4z+B6z+B7z−B1z+B2z+B5z+B8zB2x+B3x+B4x+B8x+B1y+B2y+B5y+B8y−B1x+B4x+B5x+B6x+B3y+B4y+B6y+B7y.Figure 7
The full-tensor measurement system of the hexahedral magnetic gradient.The structure comparison of seven kinds of full-tensor measurement systems is shown in Table1.Table 1
Comparison table of measurement system structure.
Element systemDimensionNumber of magnetometersNumber of independent componentsObservation pointsCross shape246Plane centerEquilateral triangle shape236Plane centerRight triangle shape236Plane centerSquare shape246Plane centerRegular tetrahedron shape349Three-dimensional centerRight-angle tetrahedron shape349Three-dimensional centerHexahedron shape389Three-dimensional centerIt can be seen from Table1 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. The three-dimensional shape measurement system can measure 9 independent components individually, while the two-dimensional plane can measure 6 independent components.


## 3.1. Cross Measurement System

The model of the cross-shaped magnetic gradient full-tensor measurement system is shown in Figure1. The system is composed of four fluxgate magnetometers. The center point is taken as the axis center (O) of the system coordinate axis. The fluxgates are marked as 1, 2, 3, and 4, respectively. The 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].Figure 1
Cross-shaped magnetic gradient full-tensor measurement system.The three-component magnetic fields measured by the four fluxgates are as follows:(7)B1=B1xB1yB1z,B2=B2xB2yB2z,B3=B3xB3yB3z,B4=B4xB4yB4z.In equation (7), Bii=1,2,3,4 is the magnetic field vector of the fluxgate No. i, and Biji=1,2,3,4;j=x,y,z is the component of the fluxgate No. i in the j direction. The expression of the magnetic gradient full tensor at the axis O of the system is as follows:(8)G=1dB1x−B3xB1y−B3yB1z−B3zB2x−B4xB2y−B4yB2z−B4zB1z−B3zB2z−B4zB3x+B4y−B1x+B2y.


## 3.2. Equilateral Triangle Measurement System

The 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. The side length of the equilateral triangle is the baseline distance (d) of the system. The geometric center of the triangle is taken as the axis center (O) of the system coordinate axis. The system model is shown in Figure 2.Figure 2
Equilateral triangle magnetic gradient full-tensor measurement system.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:(9)G=13d3B1x−B3x3B1y−B3y3B1z−B3z2B2x−B1x+B3x2B2y−B1y+B3y2B2z−B1z+B3z3B1z−B3z2B2z−B1z+B3zB1y+B3y−2B2y−3B1x−B3x.


## 3.3. Right Triangle Measurement System

The model of the right triangle magnetic gradient full-tensor measurement system is shown in Figure3. The difference from the equilateral triangle measurement system is that the connection between the centers of the three fluxgate magnetometers is isosceles right triangle. The 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 baseline distance. The full tensor of the magnetic gradient of the axis O is shown as follows:(10)G=1dB1x−B3xB1y−B3yB1z−B3zB2x−B3xB2y−B3yB2z−B3zB1z−B3zB2z−B3zB3x+B3y−B1x+B2y.Figure 3
Right triangle magnetic gradient full-tensor measurement system.


## 3.4. Square Measurement System

Figure4 is a square magnetic gradient full-tensor measurement system model. The center points of the No. 1, No. 2, No. 3, and No. 4 fluxgates are connected to form a square. The center point of the square is taken as the axis center (O) of the coordinate axis of the system. The side length d of the square is the baseline distance, and the magnetic gradient full tensor of the axis O is calculated as(11)G=12dB1x+B4x−B2x+B3xB1y+B4y−B2y+B3yB1z+B4z−B2z+B3zB1x+B2x−B3x+B4xB1y+B2y−B3y+B4yB1z+B2z−B3z+B4zB1z+B4z−B2z+B3zB1z+B2z−B3z+B4zB2x+B3x+B3y+B4y−B1x+B1y+B2y+B4x.Figure 4
Square magnetic gradient full-tensor measurement system.


## 3.5. Regular Tetrahedron Shape Measurement System

The regular tetrahedral magnetic gradient full-tensor measurement system is built by 4 fluxgate magnetometers. The connection at the center of the fluxgate forms a regular tetrahedral structure. The structure model is shown in Figure5. The 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. The magnetic field of the center O is shown as(12)G=1d2B1x−B2x+B3x32B1y−B2y+B3y32B1z−B2z+B3z3B1x+B2x−2B3x3B1y+B2y−2B3y3B1z+B2z−2B3z33B4x−B1x+B2x+B3x63B4y−B1y+B2y+B3y63B4z−B1z+B2z+B3z6.Figure 5
The full-tensor measurement system of the regular tetrahedral magnetic gradient.


## 3.6. Right-Angle Tetrahedron Shape Measurement System

The model of the right-angled tetrahedral magnetic gradient full-tensor measurement system is shown in Figure6. The line connecting the centers of the fluxgate magnetometers No. 1, No. 2, and No. 3 forms the base of an isosceles right triangle, and the length of the right-angle side is the baseline distance (d). The distance between No. 4 and No. 1 fluxgates is also d, and the center of No. 1 fluxgate is the axis center (O) of the coordinate axis of the system. The full tensor of the magnetic gradient system is(13)G=1dB2x−B1xB2y−B1yB2z−B1zB3x−B1xB3y−B1yB3z−B1zB4x−B1xB4y−B1yB4z−B1z.Figure 6
Right-angled tetrahedral magnetic gradient full-tensor measurement system.


## 3.7. Hexahedron Shape Measurement System

As shown in Figure7, the regular hexahedron shape measurement system model is composed of 8 fluxgate magnetometers centered at the apex of the regular hexahedron. The 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. The full-tensor expression of the magnetic gradient at the axis is(14)G=14dB1x+B4x+B5x+B6x−B2x+B3x+B7x+B8xB1y+B4y+B5y+B6y−B2y+B3y+B7y+B8yB1z+B4z+B5z+B6z−B2z+B3z+B7z+B8zB3x+B4x+B6x+B7x−B1x+B2x+B5x+B8xB3y+B4y+B6y+B7y−B1y+B2y+B5y+B8yB3z+B4z+B6z+B7z−B1z+B2z+B5z+B8zB1z+B4z+B5z+B6z−B2z+B3z+B4z+B8zB3z+B4z+B6z+B7z−B1z+B2z+B5z+B8zB2x+B3x+B4x+B8x+B1y+B2y+B5y+B8y−B1x+B4x+B5x+B6x+B3y+B4y+B6y+B7y.Figure 7
The full-tensor measurement system of the hexahedral magnetic gradient.The structure comparison of seven kinds of full-tensor measurement systems is shown in Table1.Table 1
Comparison table of measurement system structure.
Element systemDimensionNumber of magnetometersNumber of independent componentsObservation pointsCross shape246Plane centerEquilateral triangle shape236Plane centerRight triangle shape236Plane centerSquare shape246Plane centerRegular tetrahedron shape349Three-dimensional centerRight-angle tetrahedron shape349Three-dimensional centerHexahedron shape389Three-dimensional centerIt can be seen from Table1 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. The three-dimensional shape measurement system can measure 9 independent components individually, while the two-dimensional plane can measure 6 independent components.


## 4. 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–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 Figure8, the magnetic dipole model is established. In the Cartesian coordinate system, the position of the magnetic target is x0,y0,z0, and the position of the observation point x,y,z is the axis of the magnetic gradient full-tensor measurement system. The projection components at the observation point of the magnetic anomaly field generated are represented by Bx, By, and Bz, which are the three-component magnetic anomaly.Figure 8
Magnetic dipole model.The magnetic anomaly vector field generated by the magnetic target at the observation point is(15)B⟶=μ04πr33m⟶·r⟶r⟶r2−m⟶.In equation (15), B⟶ is the magnetic field vector; μ0=4×10−7H/m is the magnetic permeability of the vacuum; r⟶=x−x0,y−y0,z−z0 is the position vector of the observation point relative to the target point; r is the modulus of the vector r⟶; m⟶=mx,my,mz is the magnetic moment vector of the magnetic dipole. Therefore, the coordinate equation of the magnetic dipole is(16)B⟶=BxByBz,=μ04πr53x2−r23xy3xz3xy3y2−r23yz3xz3yz3z2−r2mxmymz.The observation point in Figure8 is the center point O in Figure 1 to Figure 7. The 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. The 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:(17)r1=x+d2,y,z,r2=x,y+d2,z,r3=x−d2,y,z,r4=x,y−d2,z.In equation (17), r1,r2,r3,r4 are the position coordinates of No. 1 to No. 4, with point O as the reference point. The magnetic field value of each magnetometer can be obtained, and then the full tensor of the magnetic gradient can be gotten.


### 4.2. Simulation Experiment

Assuming that the magnetic dipole is located in a vacuum, the position of the magnetic source is3,4,5. The magnetic declination angle D=15°, the magnetic inclination angle I=30°, the magnetic moment modulus m=200A·m2, mx=mcosIcosd, my=mcosIsind, mz=sinI, 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. The 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. The 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.Figure 9
Comparison of measured values of cross-shaped magnetic gradient full-tensor system and theoretical values.Figure 10
Comparison of the measured values of the magnetic gradient full-tensor system with seven structures and theoretical values.


### 4.3. Simulation Results and Analysis

It can be seen from Figure10 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.(18)REi=fi−yii=1,2,3,….Figure 11
Comparison of measured RE of seven magnetic gradient full-tensor systems.
(a)(b)(c)(d)(e)In equation (18), fi is the measured value of each point, and yi is the theoretical value of each point.It can be seen from Figure11(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. The 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. The 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. The 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.Table 2
Arrangement of magnetic gradient full-tensor measurement system.
TensorsPros and consBxx (nT/m)Square shapeCross shapeHexahedron shapeRegular tetrahedron shapeEquilateral triangle shapeRight triangle shapeRight-angle tetrahedron shapeBxy (nT/m)Cross shapeSquare shapeHexahedron shapeRegular tetrahedron shapeEquilateral triangle shapeRight triangle shapeRight-angle tetrahedron shapeBxz (nT/m)Cross shapeSquare shapeHexahedron shapeRegular tetrahedron shapeEquilateral triangle shapeRight triangle shapeRight-angle tetrahedron shapeByy (nT/m)Cross shapeSquare shapeRight triangle shapeRight-angle tetrahedron shapeEquilateral triangle shapeRegular tetrahedron shapeHexahedron shapeByz (nT/m)Cross shapeRight triangle shapeSquare shapeEquilateral triangle shapeRight triangle shapeRight-angle tetrahedron shapeRegular tetrahedron shapeIn 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. The calculated values are shown in Table 3.(19)RMSE=1n∑i=1nfi−yi2i=1,2,3,....,n.Table 3
RMSE comparison of measured values of seven kinds of magnetic gradient full-tensor system.
TensorStructuresCross shapeEquilateral triangle shapeRight triangle shapeSquare shapeRegular tetrahedron shapeRight-angle tetrahedron shapeRight triangle shapeBxx (nT/m)0.00930.24680.56520.00110.05650.56520.0187Bxy (nT/m)0.00660.12200.41820.01010.10360.41820.0133Bxz (nT/m)0.00800.29230.49770.00860.13660.49770.0161Byy (nT/m)0.00900.12170.39290.01642.91410.39297.8728Byz (nT/m)0.00190.26130.47360.02073.55120.47360.0039In equation (19), fi is the measured value of each point, and yi is the theoretical value of each point.It can be seen from Tables2 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.


## 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–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 Figure8, the magnetic dipole model is established. In the Cartesian coordinate system, the position of the magnetic target is x0,y0,z0, and the position of the observation point x,y,z is the axis of the magnetic gradient full-tensor measurement system. The projection components at the observation point of the magnetic anomaly field generated are represented by Bx, By, and Bz, which are the three-component magnetic anomaly.Figure 8
Magnetic dipole model.The magnetic anomaly vector field generated by the magnetic target at the observation point is(15)B⟶=μ04πr33m⟶·r⟶r⟶r2−m⟶.In equation (15), B⟶ is the magnetic field vector; μ0=4×10−7H/m is the magnetic permeability of the vacuum; r⟶=x−x0,y−y0,z−z0 is the position vector of the observation point relative to the target point; r is the modulus of the vector r⟶; m⟶=mx,my,mz is the magnetic moment vector of the magnetic dipole. Therefore, the coordinate equation of the magnetic dipole is(16)B⟶=BxByBz,=μ04πr53x2−r23xy3xz3xy3y2−r23yz3xz3yz3z2−r2mxmymz.The observation point in Figure8 is the center point O in Figure 1 to Figure 7. The 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. The 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:(17)r1=x+d2,y,z,r2=x,y+d2,z,r3=x−d2,y,z,r4=x,y−d2,z.In equation (17), r1,r2,r3,r4 are the position coordinates of No. 1 to No. 4, with point O as the reference point. The magnetic field value of each magnetometer can be obtained, and then the full tensor of the magnetic gradient can be gotten.


## 4.2. Simulation Experiment

Assuming that the magnetic dipole is located in a vacuum, the position of the magnetic source is3,4,5. The magnetic declination angle D=15°, the magnetic inclination angle I=30°, the magnetic moment modulus m=200A·m2, mx=mcosIcosd, my=mcosIsind, mz=sinI, 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. The 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. The 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.Figure 9
Comparison of measured values of cross-shaped magnetic gradient full-tensor system and theoretical values.Figure 10
Comparison of the measured values of the magnetic gradient full-tensor system with seven structures and theoretical values.


## 4.3. Simulation Results and Analysis

It can be seen from Figure10 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.(18)REi=fi−yii=1,2,3,….Figure 11
Comparison of measured RE of seven magnetic gradient full-tensor systems.
(a)(b)(c)(d)(e)In equation (18), fi is the measured value of each point, and yi is the theoretical value of each point.It can be seen from Figure11(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. The 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. The 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. The 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.Table 2
Arrangement of magnetic gradient full-tensor measurement system.
TensorsPros and consBxx (nT/m)Square shapeCross shapeHexahedron shapeRegular tetrahedron shapeEquilateral triangle shapeRight triangle shapeRight-angle tetrahedron shapeBxy (nT/m)Cross shapeSquare shapeHexahedron shapeRegular tetrahedron shapeEquilateral triangle shapeRight triangle shapeRight-angle tetrahedron shapeBxz (nT/m)Cross shapeSquare shapeHexahedron shapeRegular tetrahedron shapeEquilateral triangle shapeRight triangle shapeRight-angle tetrahedron shapeByy (nT/m)Cross shapeSquare shapeRight triangle shapeRight-angle tetrahedron shapeEquilateral triangle shapeRegular tetrahedron shapeHexahedron shapeByz (nT/m)Cross shapeRight triangle shapeSquare shapeEquilateral triangle shapeRight triangle shapeRight-angle tetrahedron shapeRegular tetrahedron shapeIn 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. The calculated values are shown in Table 3.(19)RMSE=1n∑i=1nfi−yi2i=1,2,3,....,n.Table 3
RMSE comparison of measured values of seven kinds of magnetic gradient full-tensor system.
TensorStructuresCross shapeEquilateral triangle shapeRight triangle shapeSquare shapeRegular tetrahedron shapeRight-angle tetrahedron shapeRight triangle shapeBxx (nT/m)0.00930.24680.56520.00110.05650.56520.0187Bxy (nT/m)0.00660.12200.41820.01010.10360.41820.0133Bxz (nT/m)0.00800.29230.49770.00860.13660.49770.0161Byy (nT/m)0.00900.12170.39290.01642.91410.39297.8728Byz (nT/m)0.00190.26130.47360.02073.55120.47360.0039In equation (19), fi is the measured value of each point, and yi is the theoretical value of each point.It can be seen from Tables2 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.


## 5. Conclusion

In this paper, seven magnetic gradient full-tensor measurement system models are studied in detail, and the full-tensor measurement theory of the tensor measurement array is analyzed. The magnetic dipole model is used to simulate the measurement system. Different measurement systems are quantitatively compared and analyzed. The 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. The three-dimensional magnetic gradient full-tensor measurement system can obtain many independent components and can directly obtain the vertical component. The magnetic gradient full-tensor measurement system should be selected to use according to actual needs. The research results in this paper provide a theoretical reference for the construction and application development of the magnetic gradient full-tensor system.



---

*Source: 6688364-2021-03-20.xml*