CN115236755B - Tensor eigenvalue-based aeromagnetic anomaly boundary detection method and device - Google Patents
Tensor eigenvalue-based aeromagnetic anomaly boundary detection method and device Download PDFInfo
- Publication number
- CN115236755B CN115236755B CN202210874914.5A CN202210874914A CN115236755B CN 115236755 B CN115236755 B CN 115236755B CN 202210874914 A CN202210874914 A CN 202210874914A CN 115236755 B CN115236755 B CN 115236755B
- Authority
- CN
- China
- Prior art keywords
- gradient
- boundary
- data matrix
- aeromagnetic
- tensor data
- 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.)
- Active
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/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- 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/15—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for use during transport, e.g. by a person, vehicle or boat
- G01V3/16—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for use during transport, e.g. by a person, vehicle or boat specially adapted for use from aircraft
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The application discloses a method, a device and a storage medium for detecting a boundary of a magnetic anomaly based on tensor eigenvalue, wherein the method comprises the following steps: acquiring a full magnetic gradient tensor data matrix T; establishing a boundary detection function R and a depth resolution gain factor T according to the full magnetic gradient tensor data matrix T z The method comprises the steps of carrying out a first treatment on the surface of the Respectively calculating the horizontal gradient R of R in the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z The method comprises the steps of carrying out a first treatment on the surface of the According to the horizontal gradient R in the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z Depth resolution gain factor T z Establishing an equilibrium boundary recognition filter MF by using a full magnetic force gradient tensor data matrix T; and according to the equilibrium boundary recognition filter MF and the actually measured aeromagnetic data, the aeromagnetic anomaly boundary detection is realized. The method and the device for detecting the boundary of the aeromagnetic anomaly based on the tensor eigenvalue and the storage medium can improve the calculation stability and the boundary of the false aeromagnetic anomaly geologic body, can clearly give the distribution characteristics of the deeper target geologic body, and have higher resolution and higher precision.
Description
Technical Field
The application relates to the technical field of aeromagnetic measurement, in particular to a method and a device for detecting an aeromagnetic anomaly boundary based on tensor eigenvalues.
Background
The aeromagnetometer (such as an optical pump type, a nuclear rotation type and a fluxgate type) system is installed in an aircraft, and magnetic or magnetic-related ore bodies are searched by observing geomagnetic field parameters (such as total geomagnetic field intensity T or total magnetic field abnormality delta T or gradients thereof) so as to know geological structures, perform magnetic mapping, solve the problems of urban and engineering stability, archaeology and the like.
Aeromagnetic surveying is mainly to study and measure magnetic anomaly fields. The magnetic abnormal field is an additional magnetic field generated by the ferromagnetic geological body in the crust under the action of the geomagnetic field. The aviation magnetic force measurement data are comprehensive reflection of magnetic field information of magnetic bodies with different depths, different forms and different scales on an observation surface. However, due to errors of the measurement data or superposition of magnetic fields, the measurement data is difficult to distinguish, and difficulty is brought to geological interpretation work.
The engineering technology and the research and development technology of the magnetic gradient tensor detection instrument are continuously developed and matured, and the problems are correspondingly developed when the magnetic tensor data are used for analysis and treatment. The magnetic tensor data is the gradient of the vector component of the magnetic field, contains magnetic field information, can reflect vector magnetic moment information of a target body, has the advantages of high precision, high resolution and multiple parameters, can be used for describing the magnetization direction and geometric form of a field source body, and improves the resolution of the target body.
At present, the resolution of a target geologic body can be improved by constructing an equalization boundary recognition filter BS method based on tensor eigenvalues.
Based on this, the inventor of the present application found that the current equalizing boundary recognition filter BS method improves a certain recognition effect, but in the actual aeromagnetic data processing and interpretation process, the method has low precision, poor stability and noise immunity, and can not effectively equalize the amplitude of the deep and shallow anomalies, especially when the positive and negative aeromagnetic anomalies are superimposed and mutually influenced, false boundary results are easy to be generated, and interference and misleading can be generated to later construction interpretation.
The information disclosed in this background section is only for enhancement of understanding of the general background of the application and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person of ordinary skill in the art.
Disclosure of Invention
The application aims to provide a method and a device for detecting aeromagnetic anomaly boundaries based on tensor eigenvalues, which can solve the problems that the existing aeromagnetic data processing results are low in precision, poor in stability and noise resistance and incapable of effectively balancing the amplitudes of deep and shallow anomalies, and particularly false boundary results are easy to generate when positive and negative aeromagnetic anomalies are mutually overlapped and mutually influenced, and interference and misleading can be generated to later construction explanation.
In order to achieve the above object, the present application provides a method for detecting a boundary of an aeromagnetic anomaly based on a tensor eigenvalue, comprising: acquiring a full magnetic gradient tensor data matrix T, wherein the full magnetic gradient tensor data matrix T comprises 9 gradient components in the x, y and z directions respectively under a three-dimensional rectangular coordinate system; establishing a boundary detection function R and a depth resolution gain factor T according to the full magnetic force gradient tensor data matrix T z The method comprises the steps of carrying out a first treatment on the surface of the Respectively calculating the horizontal gradient R of R in the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z The method comprises the steps of carrying out a first treatment on the surface of the According to the x-direction horizontal gradient R x Horizontal gradient R in y direction y And a z-direction vertical gradient R z Depth resolution gain factor T z Establishing an equilibrium boundary recognition filter MF by using a full magnetic force gradient tensor data matrix T; and according to the equilibrium boundary recognition filter MF and the actually measured aeromagnetic data, the aeromagnetic anomaly boundary detection is realized.
In one embodiment of the present application, the horizontal gradient R in the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z Depth resolution gain factor T z Establishing the equilibrium boundary recognition filter MF by the full magnetic force gradient tensor data matrix T includes:
the equalization boundary identifying filter MF is established according to the following formula, which comprises:
wherein k is an adjustment coefficient with a value of 0-1, and max (|R|) is the maximum value of the boundary detection function R.
In one embodiment of the present application, acquiring the full magnetic gradient tensor data matrix T includes:
obtaining actually measured aeromagnetic data; and determining a gradient component of the aeromagnetic data in the x, y and z directions according to the three-dimensional rectangular coordinate system, wherein each gradient component of the aeromagnetic data in the x, y and z directions respectively, and 9 gradient components form full tensor magnetic gradient data.
In one embodiment of the application, the
In one embodiment of the present application, the establishing the boundary detection function R according to the full magnetic force gradient tensor data matrix T includes: calculating three eigenvalues lambda of the full magnetic force gradient tensor data matrix T according to the matrix 1 、λ 2 、λ 3 The method comprises the steps of carrying out a first treatment on the surface of the Calculating according to the full magnetic force gradient tensor data matrix T to obtain a total modulus M of the matrix; according to the eigenvalue lambda of the full magnetic force gradient tensor data matrix 1 、λ 2 、λ 3 And the total modulus M of the full magnetic force gradient tensor data matrix, establishing a boundary detection function R, wherein R=lambda 1 ·λ 2 ·λ 3 ·M。
In one embodiment of the present application, the total module value M of the full magnetic gradient tensor contains information of all 9 tensor elements, and the maximum value thereof corresponds to the boundary of the geological body,
in one embodiment of the present application, the depth resolution gain factor T is established according to the full magnetic force gradient tensor data matrix T z Comprising the following steps: calculating a depth resolution gain factor T according to the following formula z The formula is as follows:
in one embodiment of the present application, the adjustment coefficient k of the balance depth abnormality is calculated according to the following formula:
wherein max|T z I represents T z A maximum value; max|T zz I represents T zz Minimum value.
In order to achieve the above object, the present application further provides an apparatus for detecting an aeromagnetic anomaly boundary based on a tensor eigenvalue, including:
the acquisition module is used for acquiring a full magnetic force gradient tensor data matrix T, wherein the full magnetic force gradient tensor data matrix T comprises 9 gradient components in the x, y and z directions respectively under a three-dimensional rectangular coordinate system;
the boundary detection function establishing module is used for establishing a boundary detection function R and a depth resolution gain factor T according to the full magnetic force gradient tensor data matrix T z ；
Gradient calculation modules for calculating horizontal gradients R in the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z ；
An equalization boundary recognition filter establishing module for establishing a horizontal gradient R according to the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z Depth resolution gain factor T z Establishing an equilibrium boundary recognition filter MF by using a full magnetic force gradient tensor data matrix T;
and the detection module is used for realizing the detection of the aeromagnetic anomaly boundary according to the equalization boundary identification filter MF and the actually measured aeromagnetic data.
To achieve the above object, the present application provides a storage medium storing computer-executable instructions for performing the above-described method for detecting an aeromagnetic anomaly boundary based on tensor feature values.
Compared with the prior art, according to the method for detecting the boundary of the aeromagnetic anomaly based on the tensor eigenvalue, a reasonable method for detecting the boundary of the target geologic body of balanced aeromagnetic data is newly constructed, the boundary of multi-source field objects with different burial depths can be better detected, the boundary identification result is more converged, the interference of magnetization direction and noise on the result is effectively avoided, and the calculation stability and false aeromagnetic geologic body boundary are improved; the effect of the geological bodies with different depths can be balanced by utilizing the ratio function, so that the distribution characteristics of the geological bodies with deeper targets can be clearly given, and the geological bodies with deeper targets have higher resolution and higher precision.
Drawings
FIG. 1 is a flow chart of a method of detecting a boundary of an aeromagnetic anomaly based on tensor eigenvalues according to an embodiment of the present application;
FIG. 2 is a schematic diagram of a device for detecting boundary of aeromagnetic anomaly based on tensor eigenvalues according to an embodiment of the present application.
The main reference numerals illustrate:
the system comprises a 1-acquisition module, a 2-boundary detection function establishment module, a 3-gradient calculation module, a 4-equalization boundary identification filter establishment module and a 5-detection module.
Detailed Description
The following detailed description of embodiments of the application is, therefore, to be taken in conjunction with the accompanying drawings, and it is to be understood that the scope of the application is not limited to the specific embodiments.
Throughout the specification and claims, unless explicitly stated otherwise, the term "comprise" or variations thereof such as "comprises" or "comprising", etc. will be understood to include the stated element or component without excluding other elements or components.
The embodiment of the application provides a method for detecting a boundary of an aeromagnetic anomaly based on a tensor eigenvalue, referring to fig. 1, which is a flow chart of the method for detecting the boundary of the aeromagnetic anomaly based on the tensor eigenvalue, comprising the following steps: step S1-step S5.
Step 1, acquiring a full magnetic gradient tensor data matrix T, wherein the full magnetic gradient tensor data matrix T comprises 9 gradient components in the x, y and z directions respectively under a three-dimensional rectangular coordinate system.
In one implementation, the full magnetic gradient tensor data matrix T may be obtained by actually measuring the full magnetic gradient tensor data, or may be obtained by:
obtaining actually measured aeromagnetic data;
and determining a gradient component of the magnetic field data in the x, y and z directions according to the three-dimensional rectangular coordinate system, wherein 9 gradient components form full tensor magnetic gradient data. Specifically, the aeromagnetic data is a magnetic anomaly signal caused by an underground magnetic geologic body.
Wherein, the liquid crystal display device comprises a liquid crystal display device,
step 2, establishing a boundary detection function R and a depth resolution gain factor T according to the full magnetic gradient tensor data matrix T z 。
In one implementation, establishing the boundary detection function R according to the full magnetic force gradient tensor data matrix T in step 2 may include:
step 201, calculating three eigenvalues λ of the matrix according to the full magnetic force gradient tensor data matrix T 1 、λ 2 、λ 3 Wherein the eigenvalue lambda of the full magnetic force gradient tensor data matrix 1 、λ 2 、λ 3 All correspond to the boundary of the aeromagnetic data target geologic body. But characteristic value lambda 1 、λ 2 、λ 3 The resolution of the deep target geologic body is not high.
Step 202, calculating to obtain a total module value M of the matrix according to the full magnetic force gradient tensor data matrix T.
Specifically, the total modulus M of the full magnetic force gradient tensor comprises all 9The information of the tensor element, the maximum of which corresponds to the boundary of the geologic volume, wherein,its maximum value corresponds to the boundary of the geologic body, but its recognition accuracy is low.
Step 203, according to the eigenvalue lambda of the full magnetic force gradient tensor data matrix 1 、λ 2 、λ 3 And establishing a boundary detection function R by the total modulus M of the full magnetic force gradient tensor data matrix.
R has both tensor eigenvalue and total module value properties, which can improve the recognition accuracy of shallow target geologic bodies, but the resolution of detecting deep target geologic bodies is low, so the capability of detecting the boundary of deep target geologic bodies needs to be further improved.
Wherein r=λ 1 ·λ 2 ·λ 3 ·M。
In one implementation, the depth resolution gain factor T is established in step 2 from the full magnetic force gradient tensor data matrix T z May include:
calculating a depth resolution gain factor T according to the following formula z The formula is as follows:
thus, the vertical detection capability, that is, the capability of detecting the boundary of the deep target geologic body can be improved.
Step 3, respectively calculating the horizontal gradient R of R in the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z Wherein
Step 4, according to the horizontal gradient R in the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z Depth resolution gain factor T z Establishing an equilibrium boundary recognition filter M by using a full magnetic force gradient tensor data matrix TF；
Wherein the expression is:
namely:
wherein k is an adjusting coefficient for balancing the abnormal depth part, and is generally 0-1, and is used for adjusting the amplitude of the boundary recognition result, so that the convergence of the result is improved, and the formula has mathematical significance. The values of K under different geological conditions are different, and the preferred values of K refer to the following steps.
Specifically, the selection of the preferred value of K may be achieved by:
obtaining the maximum value max (R) of the boundary detection function R;
from the full magnetic gradient tensor data matrix T and the depth resolution gain factor T z Calculate max|T z I and max I T zz The value of an adjusting coefficient k of the abnormity of the balance depth part is further determined;
wherein max|T z I represents T z A maximum value; max|T zz I represents T zz Minimum value.
Calculating an improved boundary detection equalization filter MF according to the following formula;
and step 5, detecting the boundary of the aeromagnetic anomaly according to the equalization boundary identification filter MF and the actually measured aeromagnetic data.
And obtaining an accurate aeromagnetic boundary according to the equalization boundary recognition filter MF and the actually measured aeromagnetic data, and further obtaining the position of the geological target body of the aeromagnetic data.
After the equalization boundary MF filter is constructed according to the steps, the position of the aeromagnetic data geological target body can be obtained, so that the boundary, depth, yield, scale, field distribution law, physical properties and the like of a structural body field source can be further accurately deduced, and the method has important significance for the problems of dividing a geodetic structural unit, carrying out structural partition, determining the position of a fracture structural band, distinguishing the distribution of different lithology and stratum, carrying out physical property mapping and the like.
In contrast to existing BS expressions:
compared with the existing BS filter, the method for detecting the boundary of the aeromagnetic data target geologic body is newly constructed, the boundary of the multi-source field objects with different burial depths can be better detected, the boundary identification result is more converged, the interference of magnetization direction and noise on the result is effectively avoided, and the calculation stability and false aeromagnetic data geologic body boundary are improved; the effect of the geological bodies with different depths can be balanced by utilizing the ratio function, so that the distribution characteristics of the geological bodies with deeper targets can be clearly given, and the geological bodies with deeper targets have higher resolution and higher precision.
The embodiment of the application also provides a device for detecting the boundary of the aeromagnetic anomaly based on the tensor eigenvalue, please refer to fig. 2, which is a schematic structural diagram of the device for detecting the boundary of the aeromagnetic anomaly based on the tensor eigenvalue, comprising: the system comprises an acquisition module 1, a boundary detection function establishment module 2, a gradient calculation module 3, an equalization boundary identification filter establishment module 4 and a detection module 5.
The acquisition module 1 is configured to acquire a full magnetic force gradient tensor data matrix T, where the full magnetic force gradient tensor data matrix T includes 9 gradient components in x, y, and z directions of magnetic field components in a three-dimensional rectangular coordinate system.
The boundary detection function establishing module 2 is used for establishing a boundary detection function R and a depth resolution gain factor T according to the full magnetic gradient tensor data matrix T z 。
The gradient calculation module 3 is used for respectively calculating the horizontal gradient R of R in the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z 。
The equalization boundary recognition filter establishment module is used for carrying out the horizontal gradient R according to the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z Depth resolution gain factor T z The full magnetic force gradient tensor data matrix T establishes an equilibrium boundary recognition filter MF.
The detection module 5 is used for realizing the detection of the aeromagnetic anomaly boundary according to the equalization boundary recognition filter MF and the actually measured aeromagnetic data.
The embodiment of the application also provides a storage medium, which stores computer executable instructions, including a program for executing the method for detecting the boundary of the aeromagnetic anomaly based on the tensor eigenvalue, and the computer executable instructions can execute the method in any of the method embodiments.
The storage medium may be any available medium or data storage device that can be accessed by a computer, including, but not limited to, magnetic storage (e.g., floppy disks, hard disks, magnetic tape, magneto-optical disks (MOs), etc.), optical storage (e.g., CD, DVD, BD, HVD, etc.), and semiconductor storage (e.g., ROM, EPROM, EEPROM, nonvolatile storage (NAND FLASH), solid State Disk (SSD)), etc.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing descriptions of specific exemplary embodiments of the present application are presented for purposes of illustration and description. It is not intended to limit the application to the precise form disclosed, and obviously many modifications and variations are possible in light of the above teaching. The exemplary embodiments were chosen and described in order to explain the specific principles of the application and its practical application to thereby enable one skilled in the art to make and utilize the application in various exemplary embodiments and with various modifications as are suited to the particular use contemplated. It is intended that the scope of the application be defined by the claims and their equivalents.
Claims (10)
1. The method for detecting the boundary of the aeromagnetic anomaly based on the tensor eigenvalue is characterized by comprising the following steps:
acquiring a full magnetic gradient tensor data matrix T, wherein the full magnetic gradient tensor data matrix T comprises 9 gradient components in the x, y and z directions respectively under a three-dimensional rectangular coordinate system;
establishing a boundary detection function R and a depth resolution gain factor T according to the full magnetic force gradient tensor data matrix T z ；
Respectively calculating the horizontal gradient R of R in the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z ；
According to the x-direction horizontal gradient R x Horizontal gradient R in y direction y And a z-direction vertical gradient R z Depth resolution gain factor T z Establishing an equilibrium boundary recognition filter MF by using a full magnetic force gradient tensor data matrix T;
and according to the equilibrium boundary recognition filter MF and the actually measured aeromagnetic data, the aeromagnetic anomaly boundary detection is realized.
2. The method for detecting a boundary of an aeromagnetic anomaly according to claim 1, wherein the horizontal gradient R in the x-direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z Depth resolution gain factor T z Establishing the equilibrium boundary recognition filter MF by the full magnetic force gradient tensor data matrix T includes:
the equalization boundary recognition filter MF is established according to the following formula:
wherein k is an adjustment coefficient with a value of 0-1, and max (|R|) is the maximum value of the boundary detection function R.
3. The method of claim 1, wherein obtaining a full magnetic gradient tensor data matrix T comprises:
obtaining actually measured aeromagnetic data;
and determining a gradient component of the aeromagnetic data in the x, y and z directions according to the three-dimensional rectangular coordinate system, wherein each gradient component of the aeromagnetic data in the x, y and z directions respectively, and 9 gradient components form full tensor magnetic gradient data.
4. A method of detecting a boundary of a magnetic anomaly as claimed in claim 1 or claim 3 wherein the
5. A method of boundary detection of aero-magnetic anomalies as claimed in claim 3, wherein said establishing a boundary detection function R from said full magnetic gradient tensor data matrix T comprises:
calculating three eigenvalues lambda of the full magnetic force gradient tensor data matrix T according to the matrix 1 、λ 2 、λ 3 ；
Calculating according to the full magnetic force gradient tensor data matrix T to obtain a total modulus M of the matrix;
according to the eigenvalue lambda of the full magnetic force gradient tensor data matrix 1 、λ 2 、λ 3 And the total modulus M of the full magnetic force gradient tensor data matrix, establishing a boundary detection function R, wherein R=lambda 1 ·λ 2 ·λ 3 ·M。
6. The method for detecting boundary of aeromagnetic anomaly of claim 5, wherein the total modulus M of the full magnetic gradient tensor data matrix comprises all 9 tensor elementsInformation, the maximum of which corresponds to the boundary of the geologic volume,
7. the method for detecting the boundary of the aeromagnetic anomaly according to claim 1 or 5, wherein the depth resolution gain factor T is established according to the full magnetic force gradient tensor data matrix T z Comprising the following steps:
calculating a depth resolution gain factor T according to the following formula z The formula is as follows:
8. the method for detecting a boundary of an aeromagnetic anomaly of claim 5,
calculating an adjustment coefficient k of the abnormality of the equalization depth part according to the following formula:
wherein max|T z I represents T z A maximum value; max|T zz I represents T zz Minimum value.
9. An apparatus for detecting a boundary of a magnetic anomaly based on a tensor eigenvalue, comprising:
the acquisition module is used for acquiring a full magnetic force gradient tensor data matrix T, wherein the full magnetic force gradient tensor data matrix T comprises 9 gradient components in the x, y and z directions respectively under a three-dimensional rectangular coordinate system;
the boundary detection function establishing module is used for establishing a boundary detection function R and a depth resolution gain factor T according to the full magnetic force gradient tensor data matrix T z ；
Gradient calculation modules for calculating horizontal gradients R in the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z ；
An equalization boundary recognition filter establishing module for establishing a horizontal gradient R according to the x direction x Horizontal gradient R in y direction y And a z-direction vertical gradient R z Depth resolution gain factor T z Establishing an equilibrium boundary recognition filter MF by using a full magnetic force gradient tensor data matrix T;
and the detection module is used for realizing the detection of the aeromagnetic anomaly boundary according to the equalization boundary identification filter MF and the actually measured aeromagnetic data.
10. A storage medium having stored thereon computer executable instructions for performing the tensor eigenvalue based aero magnetic anomaly boundary detection method of any one of claims 1-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210874914.5A CN115236755B (en) | 2022-07-25 | 2022-07-25 | Tensor eigenvalue-based aeromagnetic anomaly boundary detection method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210874914.5A CN115236755B (en) | 2022-07-25 | 2022-07-25 | Tensor eigenvalue-based aeromagnetic anomaly boundary detection method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115236755A CN115236755A (en) | 2022-10-25 |
| CN115236755B true CN115236755B (en) | 2023-10-03 |
Family
ID=83676202
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210874914.5A Active CN115236755B (en) | 2022-07-25 | 2022-07-25 | Tensor eigenvalue-based aeromagnetic anomaly boundary detection method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115236755B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116466408B (en) * | 2023-04-13 | 2024-01-23 | 中国自然资源航空物探遥感中心 | Artificial neural network superbedrock identification method based on aeromagnetic data |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111913227A (en) * | 2020-08-11 | 2020-11-10 | 哈尔滨工程大学 | Method and system for calculating frequency characteristics of magnetic anomaly signals |
| CN112633147A (en) * | 2020-12-22 | 2021-04-09 | 电子科技大学 | Magnetic anomaly detection method based on multi-feature fusion and intelligent iterative optimization algorithm |
| CN112749493A (en) * | 2021-01-25 | 2021-05-04 | 中国自然资源航空物探遥感中心 | Geologic body boundary detection method and system based on full magnetic gradient tensor eigenvalue |
| CN113866836A (en) * | 2021-11-23 | 2021-12-31 | 重庆大学 | Multi-target boundary identification method based on normalized magnetic anomaly derivative standard deviation |
| CN113886753A (en) * | 2021-10-09 | 2022-01-04 | 中国自然资源航空物探遥感中心 | Tilt method aeromagnetic boundary detection method and device based on tensor eigenvalue |
| CN113886754A (en) * | 2021-10-09 | 2022-01-04 | 中国自然资源航空物探遥感中心 | Tensor eigenvalue-based Theta Map method aeromagnetic boundary detection method and device |
| CN114488326A (en) * | 2022-02-15 | 2022-05-13 | 中国自然资源航空物探遥感中心 | Method and system for improving capability of detecting aeromagnetic data geologic body boundary |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2689977C (en) * | 2007-06-18 | 2015-11-24 | Commonwealth Scientific And Industrial Research Organisation | Method and apparatus for detection using magnetic gradient tensor |
| RU2664488C1 (en) * | 2015-03-04 | 2018-08-17 | Инститьют Оф Минерал Рисорсиз, Чайниз Акедеми Оф Джиолоджикал Сайенсиз | Method of automatic generation of potential field data structure |
-
2022
- 2022-07-25 CN CN202210874914.5A patent/CN115236755B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111913227A (en) * | 2020-08-11 | 2020-11-10 | 哈尔滨工程大学 | Method and system for calculating frequency characteristics of magnetic anomaly signals |
| CN112633147A (en) * | 2020-12-22 | 2021-04-09 | 电子科技大学 | Magnetic anomaly detection method based on multi-feature fusion and intelligent iterative optimization algorithm |
| CN112749493A (en) * | 2021-01-25 | 2021-05-04 | 中国自然资源航空物探遥感中心 | Geologic body boundary detection method and system based on full magnetic gradient tensor eigenvalue |
| CN113886753A (en) * | 2021-10-09 | 2022-01-04 | 中国自然资源航空物探遥感中心 | Tilt method aeromagnetic boundary detection method and device based on tensor eigenvalue |
| CN113886754A (en) * | 2021-10-09 | 2022-01-04 | 中国自然资源航空物探遥感中心 | Tensor eigenvalue-based Theta Map method aeromagnetic boundary detection method and device |
| CN113866836A (en) * | 2021-11-23 | 2021-12-31 | 重庆大学 | Multi-target boundary identification method based on normalized magnetic anomaly derivative standard deviation |
| CN114488326A (en) * | 2022-02-15 | 2022-05-13 | 中国自然资源航空物探遥感中心 | Method and system for improving capability of detecting aeromagnetic data geologic body boundary |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115236755A (en) | 2022-10-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11460519B2 (en) | Method for making a magnetic gradiometer with high detection accuracy and success rate | |
| CN113886753B (en) | Tilt method aeromagnetic boundary detection method and device based on tensor eigenvalue | |
| CN110007350B (en) | Analysis method for magnetic detection method blind area | |
| Gang et al. | Detection of ferromagnetic target based on mobile magnetic gradient tensor system | |
| CN110333536B (en) | Ranging linear positioning algorithm | |
| CN110133716B (en) | Magnetic anomaly data three-dimensional inversion method based on combined model weighting function | |
| CN115236755B (en) | Tensor eigenvalue-based aeromagnetic anomaly boundary detection method and device | |
| CN112749493A (en) | Geologic body boundary detection method and system based on full magnetic gradient tensor eigenvalue | |
| CN114488326B (en) | Method and system for improving boundary capability of aeromagnetic data geologic body | |
| CN109799540B (en) | Volcanic rock type uranium deposit magnetic susceptibility inversion method based on geological information constraint | |
| CN106842080B (en) | A kind of magnetic field measuring device posture swing interference minimizing technology | |
| Schneider et al. | Inversion of geo-magnetic full-tensor gradiometer data | |
| CN111399066A (en) | Method for processing scalar magnetic anomaly gradient signal based on orthogonal basis function | |
| Yin et al. | Magnetic dipole localization and magnetic moment estimation method based on normalized source strength | |
| NO20120654A1 (en) | Method and device for incorporation of internal gradient and limited diffusion in NMR inversion by well logging | |
| Yin et al. | A closed-form formula for magnetic dipole localization by measurement of its magnetic field vector and magnetic gradient tensor | |
| Yin et al. | Calculation of the magnetic gradient tensor from total magnetic anomaly field based on regularized method in frequency domain | |
| CN113886754B (en) | Tensor eigenvalue-based Theta Map method aeromagnetic boundary detection method and device | |
| Ma et al. | Interpretation of magnetic anomalies by horizontal and vertical derivatives of the analytic signal | |
| Liu et al. | The IDQ curve: A tool for evaluating the direction of remanent magnetization from magnetic anomalies | |
| Li et al. | Application of Helbig integrals to magnetic gradient tensor multi-target detection | |
| CN112668146B (en) | Improved field source position estimation method based on Euler deconvolution method practicability | |
| CN112666615B (en) | Tilt method-based method and device for acquiring accurate aviation magnetic measurement data | |
| Li et al. | Preprocessed method and application of magnetic gradient tensor data | |
| Guo et al. | Three-directional analytic signal analysis and interpretation of 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 |