CN116804773A - Underwater magnetic target positioning method and device based on high-order deflection magnetic gradient tensor - Google Patents
Underwater magnetic target positioning method and device based on high-order deflection magnetic gradient tensor Download PDFInfo
- Publication number
- CN116804773A CN116804773A CN202311049409.8A CN202311049409A CN116804773A CN 116804773 A CN116804773 A CN 116804773A CN 202311049409 A CN202311049409 A CN 202311049409A CN 116804773 A CN116804773 A CN 116804773A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- point
- axis
- sensor module
- components
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/081—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices the magnetic field is produced by the objects or geological structures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Measuring Magnetic Variables (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于高阶偏导磁梯度张量的水下磁目标定位方法及装置，属于目标方位探测技术领域，该定位方法通过以磁目标所在处为原点O，以任意一点与磁目标的连线作为x轴，建立空间直角坐标系；通过搭载于水下航行器的磁传感器模块对磁目标进行探测，磁传感器模块位于任意一点A处，测量并计算得出A点的磁场三分量和磁梯度值；沿平行于x轴的方向运行水下航行器，使磁传感器模块位于任意一点B处，测量并计算得出B点的磁场三分量和磁梯度值；结合A与B之间的距离和A、B两点的磁场三分量与磁梯度值，计算得出磁目标相对于B点的位置数据。该定位方法使用传感器数量少，占用空间体积小，更容易工程化。
The invention relates to an underwater magnetic target positioning method and device based on a high-order deflection magnetic gradient tensor, which belongs to the technical field of target orientation detection. The positioning method takes the location of the magnetic target as the origin O and connects any point with the magnetic target. The connection line is used as the x-axis to establish a spatial rectangular coordinate system; the magnetic target is detected through the magnetic sensor module mounted on the underwater vehicle. The magnetic sensor module is located at any point A, and the three components of the magnetic field at point A are measured and calculated. and magnetic gradient value; run the underwater vehicle in the direction parallel to the x-axis, so that the magnetic sensor module is located at any point B, measure and calculate the three components of the magnetic field and the magnetic gradient value at point B; combine the relationship between A and B The distance and the three components of the magnetic field and the magnetic gradient value at points A and B are used to calculate the position data of the magnetic target relative to point B. This positioning method uses fewer sensors, takes up less space, and is easier to engineer.
Description
Technical Field
The application belongs to the technical field of target azimuth detection, in particular to a magnetic detection technology, and particularly relates to an underwater magnetic target positioning method and device based on a high-order bias magnetic gradient tensor.
Background
The magnetic detection positioning technology mainly utilizes a magnetic sensor to measure the magnetic field of a magnetic target, thereby realizing the positioning of the target. In practical engineering application, the measured value of the magnetic sensor is the superposition value of the magnetic field of the magnetic target and the geomagnetic field, and the geomagnetic field value is far greater than the magnetic field value of the target, but in a certain range, the gradient value of the geomagnetic field is very small, so that the geomagnetic field interference can be effectively reduced through the magnetic gradient tensor.
In the prior art, nara et al published in 2014 a journal document (Journal of AppliedPhysics,2014,115:17E504-1-17E 504-3) entitled "Moore Penrose generalized inverse of the gradient tensor in Euler's equation forlocating a magnetic dipole" (the inverse of the gradient tensor in the European equation of Mole Penrose popularization) for locating magnetic dipoles, and deduced the linear relationship between the magnetic dipole position information and the magnetic field vector and magnetic gradient tensor matrix at the position, thereby calculating the magnetic dipole position information. The method is not limited to a specific form of a magnetic array structure, and in the application process, the geomagnetic field is firstly removed through the measured value of the magnetic sensor to obtain the three-component magnetic field of the target, but the three-component magnetic field of the target is difficult to directly measure under the geomagnetic field background, so that the method is generally applied to the condition of the known geomagnetic field, otherwise, the positioning error is larger.
In addition, a paper entitled "target positioning improvement method based on magnetic gradient tensor" was published in the domestic in the United states of America et al in 2014 (in the United states of America et al, 2014, system engineering and electronic technology, 36 (7): 1250-1254), a high-order bias magnetic gradient positioning algorithm based on a regular hexahedral structure is provided, and the influence of geomagnetic field estimation errors on positioning accuracy is effectively eliminated through the difference of magnetic field intensity and magnetic gradient tensor between two symmetrical planes of the regular hexahedral structure. However, 8 vector magnetic sensors are needed for the high-order bias magnetic gradient positioning algorithm of the regular hexahedral magnetic array structure, as shown in fig. 1, the 8 vector magnetic sensors 6 are of a regular hexahedral array structure, the structure is complex, the occupied space is large, the application scene of the magnetic gradient positioning algorithm is very limited in engineering application, and particularly in the field of underwater vehicles, the underwater vehicles are generally cylindrical with the diameter smaller than 600mm, and large-volume magnetic array structures are difficult to mount.
Therefore, how to provide a magnetic target positioning method applicable to the magnetic gradient tensor of an underwater vehicle is a technical problem to be solved urgently at present.
Disclosure of Invention
Aiming at the defects existing in the prior art, the application provides an underwater magnetic target positioning method based on a high-order bias magnetic gradient tensor, which is easier to engineer by adopting a plane cross magnetic array structure and can be applied to the field of underwater vehicles and the environment of unknown geomagnetic fields.
The application provides an underwater magnetic target positioning method based on a high-order bias magnetic gradient tensor, which comprises the following steps:
s1, taking the position of a magnetic target as an original point O, and taking a connecting line of any point and the magnetic target as an x-axis to establish a space rectangular coordinate system O-xyz;
s2, detecting a magnetic target through a magnetic sensor module mounted in the underwater vehicle, wherein the magnetic sensor module adopts a plane cross array formed by four vector magnetic sensors, and the four vector magnetic sensors are used as four vertexes of the magnetic sensor module; operating the underwater vehicle, enabling the center point of the magnetic sensor module to be positioned at any point A, and measuring magnetic field three components of four vertexes of the magnetic sensor module in a space rectangular coordinate system O-xyz respectively;
s3, calculating magnetic field three components and magnetic gradient values of a central point A of the magnetic sensor module according to the magnetic field three components of four vertexes of the magnetic sensor module at the point A;
s4, the underwater vehicle moves to any point along the direction parallel to the x axis, the position of the central point of the magnetic sensor module is marked as a point B, and three components of magnetic fields of four vertexes of the magnetic sensor module at the point B in a space rectangular coordinate system O-xyz are measured;
s5, calculating magnetic field three components and magnetic gradient values of a center point B of the magnetic sensor module according to the magnetic field three components of four vertexes of the magnetic sensor module at the point B;
s6, the distance between the point A and the point B is marked as delta x, and the coordinates (a, B and c) of the midpoint M of the connecting line of the point A and the point B are calculated by combining the magnetic field three components and the magnetic gradient values of the point A and the point B, so that the distances of the magnetic target relative to the point B in the directions of the x axis, the y axis and the z axis are respectively a+/-delta x/2, B and c.
In some embodiments, the coordinates (a, B, c) of the midpoint M of the line between the point a and the point B in step S6 are calculated by the formula (1), where the formula (1) is:
（1）；
in the formula (1), B x 、B y 、B z The magnetic field components of the A point in the directions of the x axis, the y axis and the z axis are respectively; b (B) xx 、B xy 、B xz Respectively B x Magnetic gradient values in the x-axis, y-axis and z-axis directions; b (B) yx 、B yy 、B yz Respectively B y Magnetic gradient values in the x-axis, y-axis and z-axis directions; b (B) zx 、B zy 、B zz Respectively B z Magnetic gradient values in the x-axis, y-axis and z-axis directions;
B＇ x 、B＇ y 、B＇ z the magnetic field components of the point B in the directions of the x axis, the y axis and the z axis are respectively; b' xx 、B＇ xy 、B＇ xz Respectively B' x Magnetic gradient values in the x-axis, y-axis and z-axis directions; b' yx 、B＇ yy 、B＇ yz Respectively B' y Magnetic gradient values in the x-axis, y-axis and z-axis directions; b' zx 、B＇ zy 、B＇ zz Respectively B' z Magnetic gradient values in the x-axis, y-axis and z-axis directions.
In some of these embodiments, the magnetic field at midpoint A at step S3 is three-component B x 、B y 、B z Calculated by (2), the magnetic field of the ith vertex of the magnetic sensor module at the point A in the directions of the x axis, the y axis and the z axis is dividedThe amounts are respectively marked as B ix 、B iy 、B iz Wherein i=1, 2, 3, 4, and formula (2) has the expression:
（2），
three components B' of the magnetic field at point B in step S4 x 、B＇ y 、B＇ z The magnetic field components of the ith vertex of the magnetic sensor module at the point B in the directions of the x axis, the y axis and the z axis are respectively marked as B' by the calculation of (3) ix 、B＇ iy 、B＇ iz Wherein i=1, 2, 3, 4, and formula (3) has the expression:
（3）。
in some of these embodiments, the magnetic gradient value B at the point A in step S3 xn 、B yn 、B zn Calculated by the formula (4), the expression of the formula (4) is:
（4），
the magnetic gradient value of the B point, B', in the step S5 xn 、B＇ yn 、B＇ zn Calculated by the formula (5), the expression of the formula (5) is:
（5），
in the formulas (4) and (5), d is a baseline distance of the planar cross array of the magnetic sensor module.
In some of these embodiments, in step S6, when the direction of movement of the underwater vehicle from point a to point B is to be in a direction approaching the magnetic target, the distance of the magnetic target relative to point B in the x-axis direction is a- Δx/2; when the moving direction of the underwater vehicle from the point A to the point B is a direction away from the magnetic target, the distance of the magnetic target relative to the point B in the x-axis direction is a+Deltax/2.
Besides, the application also provides an underwater magnetic target positioning device based on the high-order bias magnetic gradient tensor, which comprises a magnetic sensor module and a processor in communication connection with the magnetic sensor module, wherein the processor is configured to operate the underwater magnetic target positioning method based on the high-order bias magnetic gradient tensor.
In some of these embodiments, the four vector magnetic sensors of the magnetic sensor module are each connected to the processor by wires.
Based on the scheme, the underwater magnetic target positioning method based on the high-order bias magnetic gradient tensor in the embodiment of the application can realize the positioning of the magnetic target by adopting a plane cross magnetic array structure and using 4 magnetic sensors, has the characteristics of small number of the magnetic sensors, small occupied space and easy engineering, and is suitable for being applied to the field of underwater vehicles; in addition, the application eliminates errors caused by geomagnetic fields by utilizing a differential method, and can be effectively applied to the environment without knowing the geomagnetic fields.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application and do not constitute a limitation on the application. In the drawings:
FIG. 1 is a schematic diagram of a prior art regular hexahedral magnetic array structure;
FIG. 2 is a flow chart of a method for locating an underwater magnetic target based on a high-order bias magnetic gradient tensor according to the present application;
FIG. 3 is a schematic diagram of the structure of the underwater magnetic target positioning device based on the high-order bias magnetic gradient tensor of the application;
FIG. 4 is a schematic diagram showing movement of a magnetic sensor module relative to a magnetic target during positioning according to example 1;
FIG. 5 is a comparative positioning error chart of comparative example 1 and example 1;
fig. 6 is a comparative positioning error chart of comparative example 2 and example 1.
In the figure:
1. a first magnetic sensor; 2. a second magnetic sensor; 3. a third magnetic sensor; 4. a fourth magnetic sensor; 5. a magnetic target; 6. a magnetic sensor; 7. a processor; 8. an underwater vehicle.
Detailed Description
The technical solutions in the embodiments will be clearly and completely described below with reference to the drawings in the embodiments of the present application. It will be apparent that the described embodiments are only some, but not all, embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
In the application, a magnetic target refers to an object with magnetism, such as a ship, a submarine, a mine, and the like, which needs to be detected in ocean exploration.
As shown in fig. 2, the method for positioning an underwater magnetic target based on a high-order bias magnetic gradient tensor provided by the application is used for positioning a magnetic target 5 through a magnetic sensor module with a plane cross structure, and comprises the following steps:
s1, taking the position of a magnetic target 5 as an original point O, and taking the connecting line of any point and the magnetic target 5 as an x-axis to establish a space rectangular coordinate system O-xyz;
s2, detecting a magnetic target 5 through a magnetic sensor module mounted in the underwater vehicle 8, wherein the magnetic sensor module adopts a plane cross-shaped array formed by four vector magnetic sensors, the four vector magnetic sensors (specifically, a first magnetic sensor 1, a second magnetic sensor 2, a third magnetic sensor 3 and a fourth magnetic sensor 4 as shown in FIG. 3) are taken as four vertexes of the magnetic sensor module, and the central point of the plane cross-shaped array is taken as the central point of the magnetic sensor module; the center point of the magnetic sensor module is positioned at any point A, and as shown in FIG. 4 (in order to clearly show the measurement process of the magnetic sensor module, the underwater vehicle 8 and the processor 7 are not shown in FIG. 4), three magnetic field components of four vertexes of the magnetic sensor module in a space rectangular coordinate system O-xyz are measured;
s3, calculating magnetic field three components and magnetic gradient values of a central point A of the magnetic sensor module in a space rectangular coordinate system O-xyz according to the magnetic field three components of four vertexes of the magnetic sensor module at the point A obtained in the step S2;
s4, the underwater vehicle 8 moves to any point along the direction parallel to the x axis, the position of the central point of the magnetic sensor module is marked as a point B, and as shown in FIG. 4, three magnetic field components of four vertexes of the magnetic sensor module at the point B in a space rectangular coordinate system O-xyz are measured;
s5, calculating magnetic field three components and magnetic gradient values of a center point B of the magnetic sensor module according to the magnetic field three components of four vertexes of the magnetic sensor module at the point B obtained in the step S4;
s6, the distance between the point A and the point B is recorded as delta x, and the coordinates (a, B and c) of the midpoint M of the connecting line of the point A and the point B are calculated by combining the magnetic field three components and the magnetic gradient values of the point A and the point B, as shown in fig. 4, so that the distances of the magnetic target 5 relative to the point B in the directions of the x axis, the y axis and the z axis are respectively a+/-delta x/2, B and c.
In step S6, coordinates (a, B, c) of a midpoint M of the line between the point a and the point B are calculated by the formula (1), and the formula (1) is expressed as follows:
（1）；
in the formula (1), B x 、B y 、B z The magnetic field components of the A point in the directions of the x axis, the y axis and the z axis are respectively; b (B) xx 、B xy 、B xz Respectively B x Magnetic gradient values in the x-axis, y-axis and z-axis directions; b (B) yx 、B yy 、B yz Respectively B y Magnetic gradient values in the x-axis, y-axis and z-axis directions; b (B) zx 、B zy 、B zz Respectively B z Magnetic gradient values in the x-axis, y-axis and z-axis directions;
B＇ x 、B＇ y 、B＇ z the magnetic field components of the point B in the directions of the x axis, the y axis and the z axis are respectively; b' xx 、B＇ xy 、B＇ xz Respectively B' x Magnetic gradient values in the x-axis, y-axis and z-axis directions; b' yx 、B＇ yy 、B＇ yz Respectively B' y Magnetic gradient values in the x-axis, y-axis and z-axis directions; b' zx 、B＇ zy 、B＇ zz Respectively B' z Magnetic gradient values in the x-axis, y-axis and z-axis directions.
In the above-mentioned exemplary embodiment, the underwater magnetic target positioning method based on the high-order bias magnetic gradient tensor adopts the planar cross magnetic array structure, so that the number of the magnetic sensors is small, the occupied space is small, the engineering is easier, and the method can be applied to the field of underwater vehicles and the environment of unknown geomagnetic fields.
The three components of the magnetic field measured in step S2 and step S4 are the sum of the target magnetic anomaly value and the geomagnetic field value, the target magnetic anomaly value is an inherent property of the magnetic target, the geomagnetic field in a certain region can be considered as a fixed value which does not change with time, and the error caused by the geomagnetic field is eliminated in the step S6 in the formula (1) in a differential manner.
As shown in fig. 2, the magnetic field at the midpoint a in step S3 has three components B x 、B y 、B z Calculated by the method (2), the magnetic field components of the ith vertex of the magnetic sensor module at the point A in the directions of the x axis, the y axis and the z axis are respectively marked as B ix 、B iy 、B iz Wherein i=1, 2, 3, 4, and formula (2) has the expression:
（2），
three components B' of the magnetic field at point B in step S4 x 、B＇ y 、B＇ z The magnetic field components of the ith vertex of the magnetic sensor module at the point B in the directions of the x axis, the y axis and the z axis are respectively marked as B' by the calculation of (3) ix 、B＇ iy 、B＇ iz Wherein i=1, 2, 3, 4, and formula (3) has the expression:
（3）。
as shown in FIG. 2, the magnetic gradient value B at the point A in step S3 xn 、B yn 、B zn Calculated by the formula (4), the expression of the formula (4) is:
（4），
the magnetic gradient value of the B point, B', in the step S5 xn 、B＇ yn 、B＇ zn Calculated by the formula (5), the expression of the formula (5) is:
（5），
in the formulas (4) and (5), d is a baseline distance of the planar cross array of the magnetic sensor module. It will be appreciated that the baseline distance is the distance between two opposing magnetic sensors in a planar cross array, as shown in fig. 4, where the distance between the first magnetic sensor 1 and the third magnetic sensor 3 is d and the distance between the second magnetic sensor 2 and the fourth magnetic sensor 4 is d.
In step S6, when the moving direction of the underwater vehicle 8 from the point a to the point B is a direction to approach the magnetic target 5, the distance of the magnetic target 5 relative to the point B in the x-axis direction is a- Δx/2; when the moving direction of the underwater vehicle 8 from the point a to the point B is a direction away from the magnetic target 5, the magnetic target 5 is a+Δx/2 in the x-axis direction with respect to the point B.
In addition, the application also provides an underwater magnetic target positioning device based on the high-order bias magnetic gradient tensor, which is arranged inside the underwater vehicle 8 and comprises a magnetic sensor module, and further comprises a processor 7 which is in communication connection with the magnetic sensor module, wherein the processor 7 is configured to operate the underwater magnetic target positioning method based on the high-order bias magnetic gradient tensor.
In order to ensure stable transmission between the magnetic sensor module and the processor 7, the four vector magnetic sensors of the magnetic sensor module are connected to the processor 7 through wires, respectively. It should be noted that the processor 7 may be a processor such as a computer, which is capable of running the underwater magnetic target positioning method based on the higher-order bias magnetic gradient tensor of the present application.
The underwater magnetic target positioning method based on the high-order bias magnetic gradient tensor is described below by means of a specific embodiment:
example 1
The method for positioning the underwater magnetic target based on the high-order bias magnetic gradient tensor is used for establishing a simulation parameter model, and positioning the magnetic target 5 in the simulation parameter model, and comprises the following steps:
1) Establishing a space rectangular coordinate system
Taking the position of the magnetic target 5 as an original point O, and taking the connecting line of any point and the magnetic target 5 as an x-axis to establish a space rectangular coordinate system O-xyz;
2) Operating the underwater vehicle 8, three components of the magnetic field at four vertices of the magnetic sensor module at any point A are measured by the magnetic sensor module
The magnetic field components in the x-axis, y-axis and z-axis directions of the first vertex of the magnetic sensor module at the point A (121.5,20,100) measured by the first magnetic sensor 1 are respectively B 1x =37138.604nT、B 1y =-15318.108nT、B 1z -33560.264nT; the magnetic field components of the second vertex in the x-axis, y-axis and z-axis directions measured by the second magnetic sensor 2 are respectively B 2x =37201.123nT、B 2y =-15322.817nT、B 2z -33479.872nT; the magnetic field components of the third vertex in the x-axis, y-axis and z-axis directions measured by the third magnetic sensor 3 are respectively B 3x =37173.631nT、B 3y =-15408.516nT、B 3z -33469.357nT; the magnetic field components of the fourth vertex in the x-axis, y-axis and z-axis directions measured by the fourth magnetic sensor 4 are respectively B 4x =37110.717nT、B 4y =-15402.984nT、B 4z =-33550.586nT；
3) Calculating magnetic field three-component and magnetic gradient value of A point in space rectangular coordinate system O-xyz by processor 7
Calculating the magnetic field three-component B at the point A by the method (2) x =37156.01875nT、B y =-15363.106nT、B z -33515.020nT; calculated by the formula (4), B x Magnetic gradient values in the x-axis, y-axis and z-axis directions are B xx =-35.027nT/m、B xy =90.406nT/m、B xz =-90.907nT/m；B y Magnetic gradient values in the x-axis, y-axis and z-axis directions are B yx =90.408nT/m、B yy =80.167nT/m、B yz =70.714nT/m；B z Magnetic gradient values in the x-axis, y-axis and z-axis directions are B zx =-90.907nT/m、B zy =70.714nT/m、B zz =-45.14nT/m；
4) The underwater vehicle 8 moves to any point along the direction parallel to the x-axis, the center point position of the magnetic sensor module is marked as a point B, and three components of the magnetic fields of four vertexes of the magnetic sensor module at the point B are measured
The underwater vehicle 8 moves to an arbitrary point in a direction parallel to the x-axis toward the magnetic target 5, and the center point position of the magnetic sensor module is denoted as point B, and the movement distance is Δx=1m. The magnetic field components of the first vertex of the magnetic sensor module at point B (120.5,20,100) measured by the first magnetic sensor 1 in the x-axis, y-axis and z-axis directions are B 1x =37173.631nT、B＇ 1y =-15408.516nT、B＇ 1z -33469.357nT; the magnetic field components of the second vertex in the x-axis, y-axis and z-axis directions are measured by the second magnetic sensor 2 as B', respectively 2x =37236.782nT、B＇ 2y =-15413.059nT、B＇ 2z -33387.251nT; the magnetic field components of the third vertex in the x-axis, y-axis and z-axis directions are measured by the third magnetic sensor 3 as B', respectively 3x =37208.547nT、B＇ 3y =-15500.303nT、B＇ 3z -33376.64nT; the magnetic field components of the fourth vertex in the x-axis, y-axis and z-axis directions are measured by the fourth magnetic sensor 4 as B', respectively 4x =37144.998nT、B＇ 4y =-15494.923nT、B＇ 4z =-33459.602nT；
5) Calculating magnetic field three-component and magnetic gradient value of B point in space rectangular coordinate system O-xyz by processor 7
The three components B' of the magnetic field at the point B are calculated by the method (3) x =37190.9895nT、B＇ y =-15454.20025nT、B＇ z -33423.2125nT; calculated by equation (5), B' x Magnetic gradient values in the x-, y-and z-directions are B' xx =-34.916nT/m、B＇ xy =91.784nT/m、B＇ xz =-92.717nT/m；B＇ y Magnetic gradient values in the x-, y-and z-directions are B' yx =91.787nT/m、B＇ yy =81.864nT/m、B＇ yz =72.351nT/m；B＇ z Magnetic gradient values in the x-, y-and z-directions are B' zx =-92.717nT/m、B＇ zy =72.351nT/m、B＇ zz =-46.948nT/m；
6) Calculating the midpoint coordinate between the points AB by combining the distance between the points A and B, and further obtaining the position data of the magnetic target 5 relative to the point B
The coordinates (a, B, c) of the midpoint M of the line between the a point and the B point are calculated by the formula (1) to be (121.1384, 19.9273, 99.8836), and the distance of the magnetic target 5 to the B point in the x-axis direction is a- Δx/2= 120.6384M, the distance of the magnetic target 5 to the B point in the y-axis direction is b= 19.9273M, and the distance of the magnetic target 5 to the B point in the z-axis direction is c= 99.8836M.
Comparative example 1
The magnetic sensor module of the planar cross structure array is established by adopting the linear relation between the magnetic dipole position information, the magnetic field vector at the position and the magnetic gradient tensor matrix, which is proposed by the document Moore Penrose generalized inverse of the gradient tensor in Euler's equation forlocating a magnetic dipole (the magnetic dipole is positioned by the inverse of the gradient tensor in the Euler equation promoted by Mole and Penrose), the magnetic dipole is positioned by the inverse of the gradient tensor in the Euler equation, and the simulation parameter model established in the embodiment 1 is positioned.
Comparative example 2
The simulation parameter model established in example 1 is positioned by adopting a high-order bias magnetic gradient positioning algorithm based on a regular hexahedron structure, which is proposed in paper 'target positioning improvement method based on magnetic gradient tensor'.
Fig. 5 shows a comparison of errors of a plurality of sets of data, respectively measured by the method of comparative example 1 and the method of example 1, in which the distance between the center point of the magnetic sensor module and the magnetic target is within 150m, and as can be seen from fig. 4, the positioning error of example 1 of the present application is much smaller than that of comparative example 1 in the presence of the interference of the geomagnetic field.
Fig. 6 shows a comparison of errors of sets of data, respectively measured with the method of comparative example 2 and the method of example 1, in which the distance between the center point of the magnetic sensor module and the magnetic target is within 150m, and as can be seen from fig. 5, the positioning error of example 1 of the present application is substantially the same as that of comparative example 2.
By way of illustration of various embodiments of the high order bias magnetic gradient tensor-based underwater magnetic target positioning method of the present application, it can be seen that the high order bias magnetic gradient tensor-based underwater magnetic target positioning method embodiments of the present application have at least one or more of the following advantages:
1. the underwater magnetic target positioning method based on the high-order bias magnetic gradient tensor provided by the application adopts four magnetic sensors to form a plane cross array, so that the number of the magnetic sensors is reduced, the occupied space is reduced, the cost is saved, and the engineering application is facilitated;
2. the underwater magnetic target positioning method based on the high-order bias magnetic gradient tensor effectively eliminates the interference of the geomagnetic field by a differential method, and improves the positioning precision and stability.
Finally, it should be noted that: in the present specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, and identical and similar parts between the embodiments are all enough to refer to each other.
The above embodiments are only for illustrating the technical solution of the present application and not for limiting the same; while the application has been described in detail with reference to the preferred embodiments, those skilled in the art will appreciate that: modifications may be made to the specific embodiments of the present application or equivalents may be substituted for part of the technical features thereof; without departing from the spirit of the application, it is intended to cover the scope of the application as claimed.
Claims (7)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311049409.8A CN116804773B (en) | 2023-08-21 | 2023-08-21 | Underwater magnetic target positioning method and device based on high-order deflection magnetic gradient tensor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311049409.8A CN116804773B (en) | 2023-08-21 | 2023-08-21 | Underwater magnetic target positioning method and device based on high-order deflection magnetic gradient tensor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116804773A true CN116804773A (en) | 2023-09-26 |
| CN116804773B CN116804773B (en) | 2023-11-07 |
Family
ID=88080835
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311049409.8A Expired - Fee Related CN116804773B (en) | 2023-08-21 | 2023-08-21 | Underwater magnetic target positioning method and device based on high-order deflection magnetic gradient tensor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116804773B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118566989A (en) * | 2024-06-24 | 2024-08-30 | 清华大学 | Multi-magnetic target positioning method, device, equipment and medium based on magnetic gradient tensor |
| CN119355815A (en) * | 2024-11-06 | 2025-01-24 | 西北工业大学 | A method for locating underwater magnetic targets in marine environments |
| CN119737950A (en) * | 2024-12-18 | 2025-04-01 | 哈尔滨工程大学 | Underwater matching positioning method, program, device and storage medium based on magnetic beacon gradient tensor reference map gridding and two-step attitude measurement |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7561051B1 (en) * | 2005-04-20 | 2009-07-14 | Creare Inc. | Magnet locating apparatus and method of locating a magnet using such apparatus |
| US7688072B1 (en) * | 2007-09-18 | 2010-03-30 | The United States Of America As Represented By The Secretary Of The Navy | Portable magnetic sensing system for real-time, point-by-point detection, localization and classification of magnetic objects |
| CN108931241A (en) * | 2018-03-22 | 2018-12-04 | 南京理工大学 | Magnetic target one-point positioning method based on the total field gradient of magnetic and magnetic tensor modulus gradient |
| US20200326388A1 (en) * | 2019-04-15 | 2020-10-15 | Harbin Institute Of Technology | Method for Analyzing Magnetic Detection Blind Zone |
| CN112050800A (en) * | 2020-08-19 | 2020-12-08 | 哈尔滨工程大学 | A Magnetic Gradient Tensor Localization Method Based on the Symmetrical Planar Array of Japanese-shaped Triaxial Magnetometers |
| CN112611310A (en) * | 2020-12-11 | 2021-04-06 | 哈尔滨工程大学 | Magnetic dipole target distance and direction measuring method |
| CN115728829A (en) * | 2022-11-04 | 2023-03-03 | 中北大学 | Magnetic target positioning method based on magnetic gradient full tensor |
-
2023
- 2023-08-21 CN CN202311049409.8A patent/CN116804773B/en not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7561051B1 (en) * | 2005-04-20 | 2009-07-14 | Creare Inc. | Magnet locating apparatus and method of locating a magnet using such apparatus |
| US7688072B1 (en) * | 2007-09-18 | 2010-03-30 | The United States Of America As Represented By The Secretary Of The Navy | Portable magnetic sensing system for real-time, point-by-point detection, localization and classification of magnetic objects |
| CN108931241A (en) * | 2018-03-22 | 2018-12-04 | 南京理工大学 | Magnetic target one-point positioning method based on the total field gradient of magnetic and magnetic tensor modulus gradient |
| US20200326388A1 (en) * | 2019-04-15 | 2020-10-15 | Harbin Institute Of Technology | Method for Analyzing Magnetic Detection Blind Zone |
| CN112050800A (en) * | 2020-08-19 | 2020-12-08 | 哈尔滨工程大学 | A Magnetic Gradient Tensor Localization Method Based on the Symmetrical Planar Array of Japanese-shaped Triaxial Magnetometers |
| CN112611310A (en) * | 2020-12-11 | 2021-04-06 | 哈尔滨工程大学 | Magnetic dipole target distance and direction measuring method |
| CN115728829A (en) * | 2022-11-04 | 2023-03-03 | 中北大学 | Magnetic target positioning method based on magnetic gradient full tensor |
Non-Patent Citations (2)
| Title |
|---|
| LEI XU ET AL.: "Magnetic Target Linear Location Method Using Two-Point Gradient Full Tensor", IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT, vol. 70, pages 1 - 8, XP011860157, DOI: 10.1109/TIM.2021.3084283 * |
| 刘继昊等: "基于两点磁梯度张量的磁性目标在线定位方法", 地球物理学报, vol. 60, no. 10, pages 3995 - 4003 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118566989A (en) * | 2024-06-24 | 2024-08-30 | 清华大学 | Multi-magnetic target positioning method, device, equipment and medium based on magnetic gradient tensor |
| CN118566989B (en) * | 2024-06-24 | 2024-12-17 | 清华大学 | Multi-magnetic target positioning method, device, equipment and medium based on magnetic gradient tensor |
| CN119355815A (en) * | 2024-11-06 | 2025-01-24 | 西北工业大学 | A method for locating underwater magnetic targets in marine environments |
| CN119737950A (en) * | 2024-12-18 | 2025-04-01 | 哈尔滨工程大学 | Underwater matching positioning method, program, device and storage medium based on magnetic beacon gradient tensor reference map gridding and two-step attitude measurement |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116804773B (en) | 2023-11-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116804773B (en) | Underwater magnetic target positioning method and device based on high-order deflection magnetic gradient tensor | |
| CN101067554B (en) | Device and method for compensating attitude error of geomagnetic sensor | |
| CN114002725B (en) | Lane line auxiliary positioning method, device, electronic device and storage medium | |
| CN114674301B (en) | Electronic compass active compensation method and system under strong magnetic interference environment | |
| CN112050800B (en) | A Magnetic Gradient Tensor Positioning Method Based on Symmetrically Arranged Planar Array of Japanese-shaped Three-axis Magnetometers | |
| CN111412910B (en) | A method and device for positioning the ship's axial frequency magnetic field based on a rotating magnetic dipole | |
| CN108931241A (en) | Magnetic target one-point positioning method based on the total field gradient of magnetic and magnetic tensor modulus gradient | |
| CN109725361B (en) | A Magnetic Target Localization Method Based on Magnetic Gradient Tensor Invariant | |
| CN102313543A (en) | Magnetic azimuth measuring system based on giant magneto-resistance sensor, measurement method and perpendicular compensation method | |
| CN113156355B (en) | A magnetic interference compensation method for a superconducting full tensor magnetic gradient measurement device | |
| CN112051615B (en) | Underwater magnetic anomaly detection system | |
| US20210140795A1 (en) | Method and Device for Creating a Calibration Value for Calibrating an Inertial Measurement Unit for a Vehicle | |
| CN117268372B (en) | INS/GNSS integrated navigation method and system integrating magnetic navigation information | |
| CN115096294B (en) | Multi-parameter underwater magnetic target positioning method | |
| AU2022405422A1 (en) | Tilt compensation apparatus and tilt compensation method therefor | |
| CN109633540B (en) | A real-time positioning system and real-time positioning method of a magnetic source | |
| Chen et al. | A novel calibration method for tri-axial magnetometers based on an expanded error model and a two-step total least square algorithm | |
| CN119124219B (en) | Inertial navigation system installation error correction method, system, and storage medium | |
| Menezes Filho et al. | Accurate and computational-efficient analytical solutions for the extended two-step magnetometer calibration | |
| CN111220932B (en) | Unmanned aerial vehicle magnetic interference calibration method and distributed magnetic anomaly detection system | |
| CN109633541A (en) | A kind of magnetic source positioning device and source localization method | |
| CN112393745B (en) | On-line compensation method for heading error of maglev | |
| KR20170092356A (en) | System for calibrating azimuth of 3-axis magnetic sensor | |
| CN113624219A (en) | Magnetic compass ellipse fitting error compensation method based on OPTICS algorithm | |
| CN113267817A (en) | Underwater magnetic substance positioning method based on magnetic gradient tensor |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20231107 |