Three dimensional magnetic mapping process and apparatus with data filtering using an inverse computation

• US 5,543,714 A
• Filed: 12/06/1993
• Issued: 08/06/1996
• Est. Priority Date: 12/14/1992
• Status: Expired due to Fees

• Current Assignee: Commissariat a L'Energie Atomique (French Alternative Energies and Atomic Energy Commission)
• Original Assignee: Commissariat a L'Energie Atomique (French Alternative Energies and Atomic Energy Commission)
• Inventors: Blanpain, Roland, Flament, Bruno
• Primary Examiner(s): O'Shea, Sandra L.
• Assistant Examiner(s): Patidar, Jay M.
• Application Number:
  US08/161,472
• Time in Patent Office: 974 Days
• Field of Search: 326/331, 326/345, 326/244, 326/245, 326/323, 326/334, 326/357, 326/365, 364/420, 364/724.01, 364/724.02
• US Class Current: 324/345
• CPC Class Codes: G01V 3/40 specially adapted for measuring magnetic field characteristics of the earth

Abstract

In a magnetic cartography process and apparatus, measurements are performed by magnetic field sensors placed in "fish" moved in the area to be mapped. Gradients are determined on the basis of the difference between the measurements supplied by two sensors at the same time, this applying to any pair of sensors, and/or on the basis of the difference between the measurements supplied by the same sensor at two different times, this applying to any sensor, which leads to at least one map of gradients which is broken down into a map of regional field gradients and a map of local field gradients. On the basis of the local field gradients, a determination takes place of the geological field relative to the area, or the gradients of the field and then the local field, by adding the geological field to the regional field determined on the basis of its gradients by an inversion and direct calculation method.

Claims

1. Process for producing a map of a geological field at any point of an area to be mapped, comprising the steps of:

moving two or more magnetic field sensors over the area to be mapped in accordance with predetermined trajectories;

measuring an earth magnetic field of the area to be mapped with the two or more sensors and obtaining earth magnetic field measurements;

preprocessing data including the earth magnetic field measurements and positions of the two or more sensors over time, the preprocessing including

low-pass filtering of the data,

preprocessing of the data yielding preprocessed earth magnetic field measurement data;

obtaining a gradient map from a group of gradients, the group of gradients being determined from a difference between two preprocessed earth magnetic field measurement data corresponding to earth magnetic field measurements obtained by two different ones of the two or more sensors at the same time;

breaking down the obtained gradient map into a local field gradient map of an anomaly field of the area to be mapped, and a regional field gradient map;

choosing a magnetization distribution model of subsoil of the area to be mapped and calculating magnetization distribution and determining, with the magnetization distribution, at least one of the anomaly field and gradients of the anomaly field, in the form of at least one of a field and gradients of the field;

determining a regional field based on the regional field gradient map; and

determining a geological field map by addition of the field corresponding to the anomaly field and the regional field.

2. Process according to claim 1, wherein the preprocessed earth magnetic field measurement data corresponds to a given course trajectory and includes a field measurement and a position of the probe of each of the two or more sensors at a corresponding measurement time.

3. Process according to claim 1, comprising the further step of elimination of geomagnetic noise from the earth magnetic field measurements by subtracting an estimate of geomagnetic noise from the earth magnetic field measurements prior to determining the group of gradients from preprocessed earth magnetic field measurements.

4. Process according to claim 1, wherein a stochastic inversion is used for the calculation of the magnetization distribution.

5. Process according to claim 1, wherein the chosen magnetization distribution is of a dipole type, magnetic dipoles being in the subsoil, and the magnetization distribution extending over an area including the area to be mapped.

6. Process according to claim 5, wherein a curve of variations of the regional field as a function of horizontal space coordinates adopts a shape of a ramp, wherein a maximum depth of a distribution of dipoles is equal to 1/10 of a maximum dimension of a map to be made, wherein a distance between dipoles located at a given depth is at most equal to half the given depth and wherein dipoles are distributed in layers of increasing depths, each layer being located at a depth at the most equal to five times a depth of a preceding layer.

7. Apparatus for producing a geological field map of an area to be mapped, comprising:

a group of magnetic field sensors, each of the sensors having probes, the sensors supplying signals corresponding to an earth magnetic field;

means for simultaneously displacing the sensors over the area to be mapped such that positions of the probes of the sensors as a function of time are known;

a computer, the computer being programmable to

preprocess data including the earth magnetic field measurements and positions of the sensors over time by low-pass filtering of the data to yield preprocessed earth magnetic field measurement data,

obtain a gradient map from a group of gradients, the group of gradients being determined from a difference between two preprocessed earth magnetic field measurement data corresponding to earth magnetic field measurements obtained by two different ones of the sensors at the same time and including gradients of a local field and gradients of a regional field,

obtain, by inversion and direct calculation using the gradients of the local field a local field gradient map of an anomaly field of the area to be mapped, and, by integrating the regional field gradients, a regional field gradient map,

calculate magnetization distribution and determine, with the magnetization distribution, based on a magnetization distribution model of subsoil of the area to be mapped, the anomaly field in the form of a field,

determine a regional field based on the regional field gradient map, and

determine a geological field map, by addition of the field corresponding to the anomaly field and the regional field map, by addition of the gradients of the anomaly field and the gradients of the regional; and

the simultaneous displacement means simultaneously displacing the sensors such that a sufficient clearance between the probes is provided to avoid deterioration of information corresponding to the gradients determined from the difference between the two preprocessed earth magnetic field measurement data corresponding to the earth magnetic field measurements obtained by the two different ones of the sensors at the same time during the production of the gradients.

8. Apparatus according to claim 7, wherein the simultaneous displacing means maintains a spacing between the probes as the probes are moved in a horizontal plane at a determined distance from local magnetic sources, the spacing being at most equal to a smaller of the determined distance and a maximum spacing of the probes wherein the probes are able to have a same electrical conductivity environment for a given accuracy.

9. Apparatus according to claim 7, wherein the sensors are nuclear magnetic resonance sensors.

10. Process for producing a map of a geological field at any point of an area to be mapped, comprising the steps of:

moving one or more magnetic field sensors over the area to be mapped in accordance with predetermined trajectories;

measuring an earth magnetic field of the area to be mapped with the one or more sensors and obtaining earth magnetic field measurements;

preprocessing data including the earth magnetic field measurements and positions of the one or more sensors over time, the preprocessing including

low-pass filtering of the data,

preprocessing of the data yielding preprocessed earth magnetic field measurement data;

obtaining a gradient map from a group of gradients, the group of gradients being determined frown a difference between two preprocessed earth magnetic field measurement data corresponding to earth magnetic field measurements obtained by a single one of the one or more sensors at two different times;

breaking down the obtained gradient map into a local field gradient map of an anomaly field of the area to be mapped, and a regional field gradient map;

choosing a magnetization distribution model of subsoil of the area to be mapped and calculating magnetization distribution and determining, with the magnetization distribution, at least one of the anomaly field and gradients of the anomaly field, in the form of at least one of a field and gradients of the field;

determining a regional field based on the regional field gradient map; and

determining a geological field map by addition of the field corresponding to the anomaly field and the regional field.

11. Process according to claim 10, wherein the preprocessed earth magnetic field measurement data corresponds to a given course trajectory and includes a field measurement and a position of the probe of each of the one or more sensors at a corresponding measurement time.

12. Process according to claim 10, comprising the further step of elimination of geomagnetic noise from the earth magnetic field measurements by subtracting an estimate of geomagnetic noise from the earth magnetic field measurements prior to determining the group of gradients from preprocessed earth magnetic field measurements.

13. Process according to claim 10, wherein a stochastic inversion is used for the calculation of the magnetization distribution.

14. Process according to claim 10, wherein the chosen magnetization distribution is of a dipole type, magnetic dipoles being in the subsoil, and the magnetization distribution extending over an area including the area to be mapped.

15. Process according to claim 14, wherein a curve of variations of the regional field as a function of horizontal space coordinates adopts a shape of a ramp, wherein a maximum depth of a distribution of dipoles is equal to 1/10 of a maximum dimension of a map to be made, wherein a distance between dipoles located at a given depth is at most equal to half the given depth and wherein dipoles are distributed in layers of increasing depths, each layer being located at a depth at the most equal to five times a depth of a preceding layer.

16. Process for producing a map of a geological field at any point of an area to be mapped, comprising the steps of:

moving two or more magnetic field sensors over the area to be mapped in accordance with predetermined trajectories;

measuring an earth magnetic field of the area to be mapped with the two or more sensors and obtaining earth magnetic field measurements;

preprocessing data including the earth magnetic field measurements and positions of the two or more sensors over time, the preprocessing including

low-pass filtering of the data,

preprocessing of the data yielding preprocessed earth magnetic field measurement data;

obtaining a gradient map from a group of gradients, the group of gradients being determined from a difference between two preprocessed earth magnetic field measurement data corresponding to earth magnetic field measurements obtained by two different ones of the two or more sensors at the same time;

breaking down the obtained gradient map into a local field gradient map of an anomaly field of the area to be mapped, and a regional field gradient map;

choosing a magnetization distribution model of subsoil of the area to be mapped and calculating magnetization distribution and determining, with the magnetization distribution, at least one of the anomaly field and gradients of the anomaly field, in the form of at least one of a field and gradients of the field;

determining gradients of the regional field based on the regional field gradient map; and

determining a geological field gradient map by addition of the gradients of the anomaly field and the gradients of the regional field.

17. Process according to claim 16, wherein the preprocessed earth magnetic field measurement data corresponds to a given course trajectory and includes a field measurement and a position of the probe of each of the two or more sensors at a corresponding measurement time.

18. Process according to claim 16, comprising the further step of elimination of geomagnetic noise from the earth magnetic field measurements by subtracting an estimate of geomagnetic noise from the earth magnetic field measurements prior to determining the group of gradients from preprocessed earth magnetic field measurements.

19. Process according to claim 16, wherein a stochastic inversion is used for the calculation of the magnetization distribution.

20. Process according to claim 16, wherein the chosen magnetization distribution is of a dipole type, magnetic dipoles being in the subsoil, and the magnetization distribution extending over an area including the area to be mapped.

21. Process according to claim 20, wherein a curve of variations of the regional field as a function of horizontal space coordinates adopts a shape of a ramp, wherein a maximum depth of a distribution of dipoles is equal to 1/10 of a maximum dimension of a map to be made, wherein a distance between dipoles located at a given depth is at most equal to half the given depth and wherein dipoles are distributed in layers of increasing depths, each layer being located at a depth at the most equal to five times a depth of a preceding layer.

22. Process for producing a map of a geological field at any point of an area to be mapped, comprising the steps of:

moving one or more magnetic field sensors over the area to be mapped in accordance with predetermined trajectories;

measuring an earth magnetic field of the area to be mapped with the one or more sensors and obtaining earth magnetic field measurements;

preprocessing data including the earth magnetic field measurements and positions of the one or more sensors over time, the preprocessing including

low-pass filtering of the data,

preprocessing of the data yielding preprocessed earth magnetic field measurement data;

obtaining a gradient map from a group of gradients, the group of gradients being determined from a difference between two preprocessed earth magnetic field measurement data corresponding to earth magnetic field measurements obtained by a single one of the one or more sensors at two different times;

breaking down the obtained gradient map into a local field gradient map of an anomaly field of the area to be mapped, and a regional field gradient map;

choosing a magnetization distribution model of subsoil of the area to be mapped and calculating magnetization distribution and determining, with the magnetization distribution, at least one of the anomaly field and gradients of the anomaly field, in the form of at least one of a field and gradients of the field;

determining gradients of the regional field based on the regional field gradient map; and

determining a geological field gradient map by addition of the gradients of the anomaly field and the gradients of the regional field.

23. Process according to claim 22, wherein the preprocessed earth magnetic field measurement data corresponds to a given course trajectory and includes a field measurement and a position of the probe of each of the one or more sensors at a corresponding measurement time.

24. Process according to claim 22, comprising the further step of elimination of geomagnetic noise from the earth magnetic field measurements by subtracting an estimate of geomagnetic noise from the earth magnetic field measurements prior to determining the group of gradients from preprocessed earth magnetic field measurements.

25. Process according to claim 22, wherein a stochastic inversion is used for the calculation of the magnetization distribution.

26. Process according to claim 22, wherein the chosen magnetization distribution is of a dipole type, magnetic dipoles being in the subsoil, and the magnetization distribution extending over an area including the area to be mapped.

27. Process according to claim 26, wherein a curve of variations of the regional field as a function of horizontal space coordinates adopts a shape of a ramp, wherein a maximum depth of a distribution of dipoles is equal to 1/10 of a maximum dimension of a map to be made, wherein a distance between dipoles located at a given depth is at most equal to half the given depth and wherein dipoles are distributed in layers of increasing depths, each layer being located at a depth at the most equal to five times a depth of a preceding layer.

28. Apparatus for producing a geological field map of an area to be mapped, comprising:

a group of magnetic field sensors, each of the sensors having probes, the sensors supplying signals corresponding to an earth magnetic field;

means for simultaneously displacing the sensors over the area to be mapped such that positions of the probes of the sensors as a function of time are known;

a computer, the computer being programmable to

preprocess data including the earth magnetic field measurements and positions of the sensors over time by low-pass filtering of the data to yield preprocessed earth magnetic field measurement data,

obtain a gradient map from a group of gradients, the group of gradients being determined from a difference between two preprocessed earth magnetic field measurement data corresponding to earth magnetic field measurements obtained by two different ones of the sensors at the same time and including gradients of a local field and gradients of a regional field,

obtain, by inversion and direct calculation using the gradients of the local field a local field gradient map of an anomaly field of the area to be mapped, and, by integrating the regional field gradients, a regional field gradient map,

calculate magnetization distribution and determine, with the magnetization distribution, based on a magnetization distribution model of subsoil of the area to be mapped, gradients of the anomaly field, in the form of gradients of the field,

determine gradients of the regional field based on the regional field gradient map, and

determine a geological field gradient map, by addition of the gradients of the anomaly field and the gradients of the regional field; and

the simultaneous displacement means simultaneously displacing the sensors such that a sufficient clearance between the probes is provided to avoid deterioration of information corresponding to the gradients determined from the difference between the two preprocessed earth magnetic field measurement data corresponding to the earth magnetic field measurements obtained by the two different ones of the sensors at the same time during the production of the gradients.

29. Apparatus according to claim 28, wherein the simultaneous displacing means maintains a spacing between the probes as the probes are moved in a horizontal plane at a determined distance from local magnetic sources, the spacing being at most equal to a smaller of the determined distance and a maximum spacing of the probes wherein the probes are able to have a same electrical conductivity environment for a given accuracy.

30. Apparatus according to claim 28, wherein the sensors are nuclear magnetic resonance sensors.