MAGNETIC ANOMALY SENSING SYSTEM AND METHODS FOR MANEUVERABLE SENSING PLATFORMS

• US 20020171427A1
• Filed: 04/24/2001
• Published: 11/21/2002
• Est. Priority Date: 04/24/2001
• Status: Active Grant

• Current Assignee: the united states of america as represented by the secretary of the navy
• Original Assignee: The United States of America As Represented By The Secretary of Agriculture
• Inventors: Price, Brian L., Wiegert, Roy F.
• Granted Patent: US 6,476,610 B1
• Time in Patent Office: Days
• US Class Current: 324/345
• CPC Class Codes: B63G 7/06 of electromagnetic type
• CPC Class Codes: G01V 3/15 specially adapted for use during transport, e.g. by a person, vehicle or boat

Abstract

A system and method for sensing magnetic anomalies uses a gradiometer having at least one pair of triaxial magnetometer-accelerometer (TMA) sensors. Each TMA sensor has X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes that are parallel to one another and to the X,Y,Z magnetic sensing and acceleration axes of all other TMA sensors. Each TMA sensor outputs components (B_x,B_y,B_z) of local magnetic fields and components (A_x,A_y,A_z) of local gravitational acceleration fields. The components (B_x,B_y,B_z) and (A_x,A_y,A_z) output from each TMA sensor are processed to generate motion-compensated components (B_cx,B_cy,B_cz) of local magnetic fields. A difference is generated between the motion-compensated components (B_cx,B_cy,B_cz) for each pair of TMA sensors thereby generating differential vector field components (ΔB_x,ΔB_y,ΔB_z). For improved accuracy, the differential vector field components (ΔB_x,ΔB_y,ΔB_z) are adjusted using the local gravitational acceleration field components and motion-compensated local magnetic field components in order to compensate for gradiometer motion. Gradient components are generated using the differential vector field components (ΔB_x,ΔB_y,ΔB_z). In general, for each of the magnetic sensing axes, the gradient components G are defined by (ΔB_x/Δ_j,ΔB_y/Δ_j,ΔB_z/Δ_j), wherein Δ_j is a distance between a pair of TMA sensors relative to a j-th one of the X,Y,Z magnetic sensing axes. A scalar-quantity gradient contraction defined as $C^2=\sum _{i,j}{\left(G_{\mathrm{ij}}\right)}^2$

is generated for each pair of TMA sensors. The gradient contraction C² is a robust, rotationally-invariant quantity that changes monotonically with proximity to a magnetic target.

Claims

1. A magnetic anomaly sensing system, comprising:

a support that is electrically non-conductive and non-magnetic;

at least one pair of triaxial magnetometer-accelerometer (TMA) sensors coupled to said support and separated by a known distance, each of said TMA sensors having X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes that are parallel to one another and parallel to said X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes of all others of said TMA sensors, wherein each of said TMA sensors outputs X,Y,Z components (B_x,B_y,B_z) of local magnetic fields and X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields;

means for processing said X,Y,Z components (B_x,B_y,B_z) and (A_x,A_y,A_z) for each of said TMA sensors to generate motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields for each of said TMA sensors;

means for generating a difference between said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) for each said pair of TMA sensors to generate differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z);

means for generating gradient components G_ij using said differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) where i={x,y,z} and j={x,y,z} and wherein, for each of said X,Y,Z magnetic sensing axes, said gradient components G_ij are defined by (Δ

B_x/Δ

_j,Δ

B_y/Δ

_j,Δ

B_z/Δ

_j) wherein Δ

_j is a distance between said pair of TMA sensors relative to a j-th one of said X,Y,Z magnetic sensing axes; and

means for generating a scalar-quantity gradient contraction

$C^2=\underset{i,j}{\sum <br><br>}<br><br>\text{\hspace{1em}<br><br>}<br><br>{\left(G_{\mathrm{ij}}\right)}^2$

 

for each said pair of TMA sensors wherein said gradient contraction C² changes monotonically with proximity to a magnetic target.

2. A magnetic anomaly sensing system as in claim 1 wherein each said pair of TMA sensors is positioned on said support such that only one said distance A_j is non-zero.

3. A magnetic anomaly sensing system as in claim 1 further comprising:

means for generating a scalar-quantity total magnetic field B_T=[(B_cx)²+(B_cy)²+(B_cz)²]^1/2 for each of said TMA sensors; and

means for generating a sum Σ

B_T of said total magnetic field for each of said pair of TMA sensors wherein changes in said sum Σ

B_T can be correlated with said gradient contraction C².

4. A magnetic anomaly sensing system as in claim 1 wherein said at least one pair of TMA sensors comprises four TMA sensors, each of said four TMA sensors being positioned at a vertex of a planar quadrilateral, said system further comprising means for determining range and relative bearing to said magnetic target using values of said gradient contraction C² for at least two pair of said four TMA sensors positioned on opposite sides of said planar quadrilateral.

5. A magnetic anomaly sensing system as in claim 1 wherein said at least one pair of TMA sensors comprises eight TMA sensors, each of said eight TMA sensors being positioned at a vertex of a rectangular parallelepiped, said system further comprising:

means for selecting a first plane of said rectangular parallelepiped defined by four TMA sensors of said eight TMA sensors that is closest to said magnetic target as defined by values of said gradient contraction C² determined using said four TMA sensors defining said first plane; and

means for determining range and relative bearing to said magnetic target using values of said gradient contraction C² determined using said four TMA sensors defining said first plane and values of said gradient contraction C² determined using four TMA sensors of said eight TMA sensors defining a second plane, said second plane opposing said first plane and being parallel thereto.

6. A magnetic anomaly sensing system as in claim 1 further comprising means for adjusting said differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) using said X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields and said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields in order to compensate for motion of said support.

7. A method of sensing magnetic anomalies, comprising the steps of:

providing a non-magnetic and electrically non-conductive support with at least one pair of triaxial magnetometer-accelerometer (TMA) sensors coupled thereto and separated by a known distance, each of said TMA sensors having X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes that are parallel to one another and parallel to said X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes of all others of said TMA sensors, wherein each of said TMA sensors outputs X,Y,Z components (B_x,B_y,B_z) of local magnetic fields and X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields;

generating motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields for each of said TMA sensors using said X,Y,Z components (B_x,B_y,B_z) and (A_x,A_y,A_z) for each of said TMA sensors;

generating a difference between said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) for each said pair of TMA sensors to generate differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z);

generating gradient components G_ij using said differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) where i={x,y,z} and j={x,y,z} and wherein, for each of said X,Y,Z magnetic sensing axes, said gradient components G_ij are defined by (Δ

B_x/Δ

_j,Δ

B_y/Δ

_j,Δ

B_z/Δ

_j), wherein Δ

_j is a distance between said pair of TMA sensors relative to a j-th one of said X,Y,Z magnetic sensing axes; and

generating a scalar-quantity gradient contraction

$C^2=\underset{i,j}{\sum <br><br>}<br><br>\text{\hspace{1em}<br><br>}<br><br>{\left(G_{\mathrm{ij}}\right)}^2$

 

for each said pair of TMA sensors wherein said gradient contraction C² changes monotonically with proximity to a magnetic target.

8. A method according to claim 7 further comprising the step of positioning each said pair of TMA sensors on said support such that only one said distance Δ

_j is non-zero for each said pair of TMA sensors.

9. A method according to claim 7 further comprising the steps of:

generating a scalar-quantity total magnetic field B_T=[(B_cx)²+(B_cy)²+(B_cz)²]^1/2 for each of said TMA sensors;

generating a sum Σ

B_T of said total magnetic field for each of said pair of TMA sensors; and

correlating changes in said sum Σ

B_T with said gradient contraction C².

10. A method according to claim 7 wherein said at least one pair of TMA sensors comprises four TMA sensors, said method further comprising the steps of:

positioning each of said four TMA sensors at a vertex of a planar quadrilateral; and

determining range and relative bearing to said magnetic target using values of said gradient contraction C² for at least two pair of said four TMA sensors positioned on opposite sides of said planar quadrilateral.

11. A method according to claim 7 wherein said at least one pair of TMA sensors comprises eight TMA sensors, said method further comprising the steps of:

positioning each of said eight TMA sensors at a vertex of a rectangular parallelepiped;

selecting a first plane of said rectangular parallelepiped defined by four TMA sensors of said eight TMA sensors that is closest to said magnetic target as defined by values of said gradient contraction C² determined using said four TMA sensors defining said first plane; and

determining range and relative bearing to said magnetic target using values of said gradient contraction C² determined using said four TMA sensors defining said first plane and values of said gradient contraction C² determined using four TMA sensors of said eight TMA sensors defining a second plane, wherein said second plane opposes said first plane and is parallel thereto.

12. A method according to claim 7 further comprising the step of adjusting said differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) using said X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields and said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields in order to compensate for motion of said support.

14. A magnetic anomaly sensing system as in claim 13 further comprising means for adjusting said differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) using said X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields and said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields in order to compensate for motion of said support.

16. A method according to claim 15 further comprising the step of adjusting said differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) using said X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields and said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields in order to compensate for motion of said support.

18. A magnetic anomaly sensing system as in claim 17 further comprising means for adjusting said differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) using said X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields and said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields in order to compensate for motion of said support.

20. A method according to claim 19 further comprising the step of adjusting said differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) using said X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields and said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields in order to compensate for motion of said support.

13. A magnetic anomaly sensing system, comprising:

a support that is electrically non-conductive and non-magnetic;

a pair of triaxial magnetometer-accelerometer (TMA) sensors coupled to said support, each of said TMA sensors having X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes that are parallel to one another and parallel to said X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes of another of said TMA sensors, said pair of TMA sensors being positioned such that a known distance A separates said TMA sensors along only one coordinate axis of said X,Y,Z magnetic sensing axes and said X,Y,Z acceleration sensing axes, wherein each of said TMA sensors outputs X,Y,Z components (B_x,B_y,B_z) of local magnetic fields and X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields;

means for processing said X,Y,Z components (B_x,B_y,B_z) and (A_x,A_y,A_z) for each of said TMA sensors to generate motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields for each of said TMA sensors;

means for generating a difference between said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) for said pair of TMA sensors to generate differential vector field components (Δ

B^X,Δ

B_y,Δ

B_z);

means for generating gradient components G_i using said differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) where i={x,y,z} and wherein, for each of said X,Y,Z magnetic sensing axes, said gradient components G_i are defined by (Δ

B_x/Δ

,Δ

B_y/Δ

,Δ

B_z/Δ

); and

means for generating a scalar-quantity gradient contraction

$C^2=\underset{i}{\sum <br><br>}<br><br>\text{\hspace{1em}<br><br>}<br><br>{\left(G_i\right)}^2$

 

for said pair of TMA sensors wherein said gradient contraction C² changes monotonically with proximity to a magnetic target.

15. A method for sensing magnetic anomalies, comprising the steps of:

a) providing a non-magnetic and electrically non-conductive support with a pair of triaxial magnetometer-accelerometer (TMA) sensors coupled thereto, each of said TMA sensors having X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes that are parallel to one another and parallel to said X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes of another of said TMA sensors;

b) positioning said pair of TMA sensors on said non-magnetic support such that a known distance A separates said TMA sensors along only one coordinate axis of said X,Y,Z magnetic sensing axes and said X,Y,Z acceleration sensing axes;

c) moving said non-magnetic support with said pair of TMA sensors coupled thereto to a measurement location wherein each of said TMA sensors outputs X,Y,Z components (B_x,B_y,B_z) of local magnetic fields and X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields;

d) generating motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields for each of said TMA sensors using said X,Y,Z components (B_x,B_y,B_z) and (A_x,A_y,A_z) for each of said TMA sensors;

e) generating a difference between said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) for said pair of TMA sensors to generate differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z);

f) generating gradient components G_i using said differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) where i={x,y,z} and wherein, for each of said X,Y,Z magnetic sensing axes, said gradient components G_i are defined by (Δ

B_x/Δ

,Δ

B_y/Δ

,Δ

B_z/Δ

);

g) generating a scalar-quantity gradient contraction

$C^2=\underset{i}{\sum <br><br>}<br><br>\text{\hspace{1em}<br><br>}<br><br>{\left(G_i\right)}^2$

 

for said pair of TMA sensors; and

h) repeating steps c)-g) for a second measurement location wherein said gradient contraction C² changes monotonically with proximity to a magnetic target.

17. A magnetic anomaly sensing system, comprising:

a support that is electrically non-conductive and non-magnetic;

four triaxial magnetometer-accelerometer (TMA) sensors coupled to said support, each of said four TMA sensors being positioned at a vertex of a planar square of side length S and having X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes that are parallel to one another and parallel to said X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes of all others of said TMA sensors, wherein each of said TMA sensors outputs X,Y,Z components (B_x,B_y,B_z) of local magnetic fields and X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields;

means for processing said X,Y,Z components (B_x,B_y,B_z) and (A_x,A_y,A_z) for each of said TMA sensors to generate motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields for each of said TMA sensors;

means for generating a difference between said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) for pairs of said TMA sensors positioned on opposite sides of said planar square to generate differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) for each of said pairs;

means for generating gradient components G_i for each of said pairs using said differential vector field components (Δ

B_y,Δ

B_y,Δ

B_z) where i={X,Y,Z} and wherein, for each of said X,Y,Z magnetic sensing axes, said gradient components G_i are defined by (Δ

B_x/S,Δ

B_y/S,Δ

B_z/S);

means for generating a scalar-quantity gradient contraction

$C^2=\underset{i}{\sum <br><br>}<br><br>\text{\hspace{1em}<br><br>}<br><br>{\left(G_i\right)}^2$

 

for each of said pairs; and

means for determining range and relative bearing to a magnetic target using values of said gradient contraction C² for said pairs positioned on opposite sides of said planar square.

19. A method for sensing magnetic anomalies, comprising the steps of:

a) providing a non-magnetic and electrically non-conductive support with four triaxial magnetometer-accelerometer (TMA) sensors coupled thereto, each of said TMA sensors being positioned at a vertex of a planar square of side length S and having X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes that are parallel to one another and parallel to said X,Y,Z magnetic sensing axes and X,Y,Z acceleration sensing axes of all others of said TMA sensors;

b) moving said non-magnetic support with said four TMA sensors coupled thereto to a measurement location wherein each of said TMA sensors outputs X,Y,Z components (B_x,B_y,B_z) of local magnetic fields and X,Y,Z components (A_x,A_y,A_z) of local gravitational acceleration fields;

c) generating motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) of local magnetic fields for each of said TMA sensors using said X,Y,Z components (B_x,B_y,B_z) and (A_x,A_y,A_z) for each of said TMA sensors;

d) generating a difference between said motion-compensated X,Y,Z components (B_cx,B_cy,B_cz) for pairs of TMA sensors positioned on opposite sides of said planar square to generate differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) for each of said pairs;

e) generating gradient components G_i for each of said pairs using said differential vector field components (Δ

B_x,Δ

B_y,Δ

B_z) where i={x,y,z} and wherein, for each of said X,Y,Z magnetic sensing axes, said gradient components G_i are defined by (Δ

B_x/Δ

,Δ

B_y/Δ

,Δ

B_z/Δ

);

f) generating a scalar-quantity gradient contraction

$C^2=\underset{i}{\sum <br><br>}<br><br>\text{\hspace{1em}<br><br>}<br><br>{\left(G_i\right)}^2$

 

for each of said pairs of TMA sensors; and

g) repeating steps b)-f) for a second measurement location wherein said gradient contraction C² changes monotonically with proximity to a magnetic target.