CN110989005A - Weak magnetic anomaly self-adaptive real-time detection method based on scalar magnetometer array - Google Patents
Weak magnetic anomaly self-adaptive real-time detection method based on scalar magnetometer array Download PDFInfo
- Publication number
- CN110989005A CN110989005A CN201911336725.7A CN201911336725A CN110989005A CN 110989005 A CN110989005 A CN 110989005A CN 201911336725 A CN201911336725 A CN 201911336725A CN 110989005 A CN110989005 A CN 110989005A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- noise
- magnetic anomaly
- reconstructed
- detection
- 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
Images
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/40—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for measuring magnetic field characteristics of the earth
Abstract
The invention discloses a weak magnetic anomaly self-adaptive real-time detection method based on a scalar magnetometer array, which is used for solving the technical problem of low detection probability of the existing magnetic anomaly detection method. The technical scheme is that acquired magnetic field signals are decomposed into a plurality of Intrinsic Mode Functions (IMFs) through empirical mode decomposition, and the IMFs are reconstructed into two parts according to the properties of the IMFs: magnetic noise and magnetic anomaly signals. Due to the modal aliasing problem of empirical mode decomposition, the reconstructed magnetic anomaly signal contains a portion of the noise. And updating the probability density function of the noise in real time by using the reconstructed magnetic noise. Based on the updated probability density function, the corresponding geomagnetic entropy is calculated by using a moving window. And when the local magnetic entropy is smaller than the set threshold value, the detection of the magnetic anomaly is realized. The invention improves the probability of magnetic anomaly detection and improves the target detection capability under low signal-to-noise ratio.
Description
Technical Field
The invention relates to a magnetic anomaly detection method, in particular to a weak magnetic anomaly self-adaptive real-time detection method based on a scalar magnetometer array.
Background
The document "Magnetic and analysis detection using error filter, measurement science and technology,2008, Vol19(4), 045205" discloses a Magnetic anomaly detection method based on geomagnetic entropy. The method is a magnetic anomaly detection method constructed based on a mechanism that magnetic noise mode changes are caused by magnetic anomalies generated by a magnetic target. The method takes magnetic field background noise as a stable random process, takes a probability density function as a prior condition, and takes geomagnetic entropy as a criterion for magnetic anomaly detection. In the detection process, the method has low algorithm complexity and simple calculation process, and can meet the requirement of real-time detection. In the method described in the document, when the environmental condition changes or the detection is based on the moving platform, the probability density function of the background magnetic noise as the prior condition is used as the prior condition, the acquisition of the probability density function of the background magnetic noise as the prior condition becomes difficult, and the detection performance of the background magnetic noise is affected due to inaccuracy of the probability density function, so that the detection probability of the method is reduced.
Disclosure of Invention
In order to overcome the defect of low detection probability of the existing magnetic anomaly detection method, the invention provides a weak magnetic anomaly self-adaptive real-time detection method based on a scalar magnetometer array. The method decomposes an acquired magnetic field signal into a plurality of Intrinsic Mode Functions (IMFs) through empirical mode decomposition, and reconstructs the IMFs into two parts according to the properties of the IMFs: magnetic noise and magnetic anomaly signals. Due to the modal aliasing problem of empirical mode decomposition, the reconstructed magnetic anomaly signal contains a portion of the noise. And updating the probability density function of the noise in real time by using the reconstructed magnetic noise. Based on the updated probability density function, the corresponding geomagnetic entropy is calculated by using a moving window. And when the local magnetic entropy is smaller than the set threshold value, the detection of the magnetic anomaly is realized. The invention improves the probability of magnetic anomaly detection and improves the target detection capability under low signal-to-noise ratio.
The technical scheme adopted by the invention for solving the technical problems is as follows: a weak magnetic anomaly self-adaptive real-time detection method based on a scalar magnetometer array is characterized by comprising the following steps:
(a) carrying out self-adaptive decomposition on the acquired magnetic field difference value by using empirical mode decomposition to obtain a plurality of inherent mode functions and a residual error term
In the formula, ci(t) denotes an i-th natural mode function component, and r (t) denotes a residual component. d represents the number of layers of empirical mode decomposition;
(b) according to the property of the inherent mode function, respectively reconstructing magnetic noise and magnetic abnormal signals:
wherein n (t) represents reconstructed magnetic noise, and s (t) represents reconstructed magnetic anomaly signal;
(c) statistically analyzing and updating the probability density function p (-) of the reconstructed magnetic noise N (t);
(d) setting the length L of a sliding window, and calculating the geomagnetic entropy of the reconstructed magnetic anomaly signal S (t) according to the updated noise probability density function:
in the formula, i represents the ith time.
The invention has the beneficial effects that: the method decomposes an acquired magnetic field signal into a plurality of Intrinsic Mode Functions (IMFs) through empirical mode decomposition, and reconstructs the IMFs into two parts according to the properties of the IMFs: magnetic noise and magnetic anomaly signals. Due to the modal aliasing problem of empirical mode decomposition, the reconstructed magnetic anomaly signal contains a portion of the noise. And updating the probability density function of the noise in real time by using the reconstructed magnetic noise. Based on the updated probability density function, the corresponding geomagnetic entropy is calculated by using a moving window. And when the local magnetic entropy is smaller than the set threshold value, the detection of the magnetic anomaly is realized. According to the invention, the acquired magnetic field difference value is subjected to adaptive decomposition by adopting empirical mode decomposition to obtain a plurality of intrinsic mode functions and a residual error item, a magnetic noise part is reconstructed by utilizing the intrinsic mode functions obtained by decomposition, and the approximate probability density function is obtained by carrying out statistical analysis on the constructed magnetic noise, so that the key problem that the magnetic noise probability density function is difficult to obtain in the magnetic anomaly detection process due to the change of external environmental conditions or based on a motion platform is solved, the probability of magnetic anomaly detection is improved, and the target detection capability under the condition of low signal to noise ratio is improved.
The present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
Drawings
FIG. 1 is a flow chart of the weak magnetic anomaly adaptive real-time detection method based on the scalar magnetometer array.
FIG. 2 is a schematic diagram of object detection of a one-dimensional magnetometer array according to the method of the present invention.
FIG. 3 is a magnetic anomaly generated by a distant target.
FIG. 4 is a diagram showing the results of weak magnetic anomaly detection by the method of the present invention.
Detailed Description
Reference is made to fig. 1-4. The weak magnetic anomaly self-adaptive real-time detection method based on the scalar magnetometer array specifically comprises the following steps:
1. and (5) collecting the magnetic field difference value.
Two magnetometers are adopted to form a one-dimensional scalar magnetic sensor array. Sensor spacing in array BxThe magnetic sensor used a CS-L optical pump magnetometer with a sensitivity of 0.6pT at 1 m. The magnetic target moves in a horizontal plane and has a magnetic moment of (0,0,25) A.m2. The target reaches the end position (180,30, -1) m in a uniform motion manner from the starting position (-180,30, -1) m.
2. And (5) reconstructing magnetic noise and magnetic abnormal signals.
Carrying out self-adaptive decomposition on the acquired magnetic field difference value by using empirical mode decomposition, wherein the decomposition layer number d is 8, and obtaining 8 inherent mode functions and 1 residual error term:
in the formula, ci(t) denotes an i-th natural mode function component, and r (t) denotes a residual component. d represents the number of layers of empirical mode decomposition.
According to the property of the natural modal function, the first 4 terms of the natural modal function are used for reconstructing magnetic noise, and the last 4 terms are used for reconstructing magnetic anomaly signals:
wherein n (t) represents reconstructed magnetic noise, and s (t) represents reconstructed magnetic anomaly signal;
3. the probability density function of the magnetic noise is updated.
Statistically analyzing and updating the probability density function p (-) of the reconstructed magnetic noise N (t);
4. and (4) computing the geomagnetic entropy.
Setting the length L of a sliding window, and calculating the geomagnetic entropy of the reconstructed magnetic anomaly signal S (t) according to the updated noise probability density function p (·):
wherein i represents the ith time;
5. and detecting the magnetic anomaly.
And judging whether the magnetic anomaly is detected or not through the threshold value T. When I (S)i) < T, a magnetic anomaly exists, i.e., a magnetic target exists; when I (S)i) T or more, no magnetic anomaly, i.e. no magnetic target.
Claims (1)
1. A weak magnetic anomaly self-adaptive real-time detection method based on a scalar magnetometer array is characterized by comprising the following steps:
(a) carrying out self-adaptive decomposition on the acquired magnetic field difference value by using empirical mode decomposition to obtain a plurality of inherent mode functions and a residual error term
In the formula, ci(t) denotes the ith natural mode function component, and r (t) denotes the residual component; d represents the number of layers of empirical mode decomposition;
(b) according to the property of the inherent mode function, respectively reconstructing magnetic noise and magnetic abnormal signals:
wherein n (t) represents reconstructed magnetic noise, and s (t) represents reconstructed magnetic anomaly signal;
(c) statistically analyzing and updating the probability density function p (-) of the reconstructed magnetic noise N (t);
(d) setting the length L of a sliding window, and calculating the geomagnetic entropy of the reconstructed magnetic anomaly signal S (t) according to the updated noise probability density function:
in the formula, i represents the ith time.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911336725.7A CN110989005B (en) | 2019-12-23 | 2019-12-23 | Weak magnetic anomaly self-adaptive real-time detection method based on scalar magnetometer array |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911336725.7A CN110989005B (en) | 2019-12-23 | 2019-12-23 | Weak magnetic anomaly self-adaptive real-time detection method based on scalar magnetometer array |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110989005A true CN110989005A (en) | 2020-04-10 |
| CN110989005B CN110989005B (en) | 2021-12-28 |
Family
ID=70074259
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911336725.7A Active CN110989005B (en) | 2019-12-23 | 2019-12-23 | Weak magnetic anomaly self-adaptive real-time detection method based on scalar magnetometer array |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110989005B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112651385A (en) * | 2021-01-18 | 2021-04-13 | 西北工业大学青岛研究院 | Magnetic anomaly multi-feature information extraction method based on Hilbert-Huang transform |
| CN112698412A (en) * | 2020-11-27 | 2021-04-23 | 北京自动化控制设备研究所 | Shaft frequency magnetic abnormal signal processing method based on magnetic buoy |
| CN113866834A (en) * | 2021-09-15 | 2021-12-31 | 吉林大学 | Entropy filtering-based field source center position inversion method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105182429A (en) * | 2015-09-29 | 2015-12-23 | 吉林大学 | Marine controlled-source electromagnetic data seawater disturbance noise correction method |
| CN109799535A (en) * | 2019-03-14 | 2019-05-24 | 中船海洋探测技术研究院有限公司 | A kind of filtering method of full tensor magnetic gradient detection and localization data |
| CN109858109A (en) * | 2019-01-14 | 2019-06-07 | 北京工业大学 | A kind of gear signal noise-reduction method combined based on the EMD of correlation with form singular value decomposition |
| CN109871733A (en) * | 2018-09-27 | 2019-06-11 | 南京信息工程大学 | A kind of adaptive sea clutter signal antinoise method |
-
2019
- 2019-12-23 CN CN201911336725.7A patent/CN110989005B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105182429A (en) * | 2015-09-29 | 2015-12-23 | 吉林大学 | Marine controlled-source electromagnetic data seawater disturbance noise correction method |
| CN109871733A (en) * | 2018-09-27 | 2019-06-11 | 南京信息工程大学 | A kind of adaptive sea clutter signal antinoise method |
| CN109858109A (en) * | 2019-01-14 | 2019-06-07 | 北京工业大学 | A kind of gear signal noise-reduction method combined based on the EMD of correlation with form singular value decomposition |
| CN109799535A (en) * | 2019-03-14 | 2019-05-24 | 中船海洋探测技术研究院有限公司 | A kind of filtering method of full tensor magnetic gradient detection and localization data |
Non-Patent Citations (1)
| Title |
|---|
| ARIE SHEINKER 等: "Magnetic anomaly detection using entropy filter", 《MEASUREMENT SCIENCE AND TECHNOLOGY》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112698412A (en) * | 2020-11-27 | 2021-04-23 | 北京自动化控制设备研究所 | Shaft frequency magnetic abnormal signal processing method based on magnetic buoy |
| CN112651385A (en) * | 2021-01-18 | 2021-04-13 | 西北工业大学青岛研究院 | Magnetic anomaly multi-feature information extraction method based on Hilbert-Huang transform |
| CN113866834A (en) * | 2021-09-15 | 2021-12-31 | 吉林大学 | Entropy filtering-based field source center position inversion method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110989005B (en) | 2021-12-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110989005B (en) | Weak magnetic anomaly self-adaptive real-time detection method based on scalar magnetometer array | |
| CN110361778B (en) | Seismic data reconstruction method based on generation countermeasure network | |
| CN110967774B (en) | Magnetic anomaly detection method based on sensor array | |
| CN108983287B (en) | Curvelet transform anti-aliasing seismic data reconstruction method based on convex set projection algorithm | |
| CN114112400A (en) | Mechanical bearing fault diagnosis method based on multi-angle information fusion | |
| CN115200850A (en) | Mechanical equipment anomaly detection method under explicit representation of multi-point sample structure information | |
| CN108549078A (en) | A kind of radar pulse signal is across combining channel and detection method | |
| CN108305265B (en) | Real-time processing method and system for weak and small target image | |
| CN114578356A (en) | Distributed scatterer deformation monitoring method, system and equipment based on deep learning | |
| CN111181634A (en) | Distributed optical fiber vibration signal rapid positioning method | |
| CN112183407B (en) | Tunnel seismic wave data denoising method and system based on time-frequency domain spectral subtraction | |
| CN104679994A (en) | Autonomous underwater vehicle propeller fault detecting method based on wavelet single branch reconstruction | |
| CN106419884B (en) | A kind of rate calculation method and system based on wavelet analysis | |
| CN110287853B (en) | Transient signal denoising method based on wavelet decomposition | |
| CN112033656A (en) | Mechanical system fault detection method based on broadband spectrum processing | |
| CN116595465A (en) | High-dimensional sparse data outlier detection method and system based on self-encoder and data enhancement | |
| CN115828069A (en) | End-to-end magnetic anomaly signal noise reduction method based on deep learning | |
| CN116226770A (en) | Time sequence data anomaly detection method and device | |
| Pan et al. | CIAM: An adaptive 2-in-1 missing data estimation algorithm in wireless sensor networks | |
| CN113466947B (en) | Automatic dead pixel removal method applied to superconducting transient electromagnetic | |
| CN114694014A (en) | SAR image ship target detection method based on multilayer neural network | |
| CN117743836B (en) | Abnormal vibration monitoring method for bearing | |
| CN112504429A (en) | High-precision demodulation algorithm for strong interference DVS | |
| CN114325846B (en) | Magnetic anomaly detection method for suppressing noise by utilizing time coherence | |
| CN117648537B (en) | Atmospheric pollution real-time monitoring method and system based on hyperspectral technology |
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 |