A total magnetic field anomaly data processing method and device and related equipment

Abstract

The application discloses a magnetic total field anomaly data processing method and device and related equipment, and the method comprises the following steps: acquiring two-dimensional magnetic total field anomaly data uniformly distributed in a detection area, and acquiring first horizontal distribution data of the modulus of two-dimensional magnetic total field anomaly on each group of unit orthogonal bases based on the two-dimensional magnetic total field anomaly data, wherein the unit orthogonal bases are obtained by decomposing two-dimensional magnetic total field anomaly data generated by a magnetic dipole of a target body at any height; acquiring second horizontal distribution data of the energy of two-dimensional magnetic total field anomaly in a two-dimensional space spanned by each group of unit orthogonal bases based on the first horizontal distribution data; and determining the horizontal coordinate of the target body based on the second horizontal distribution data. The application converts two-dimensional magnetic total field anomaly data into two-dimensional energy distribution, can improve the signal-to-noise ratio of data, effectively improves the accuracy of target detection, can determine the horizontal coordinate of the target, and realizes the positioning of the target in the horizontal direction.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of magnetic exploration, more particularly to a magnetic total field anomaly data processing method and device and related equipment. BACKGROUND

[0002] Magnetic anomaly detection (MAD) is one of the most suitable geophysical techniques for locating and mapping the distribution of ferromagnetic metal objects, and has a wide range of applications in the fields of resource exploration, unexploded ordnance detection, underwater vehicle detection, archaeology, and other fields related to national economic security and cultural construction. The most common way to implement magnetic anomaly detection is to measure one-dimensional or two-dimensional magnetic total field. One-dimensional magnetic total field anomaly data can only be used for target detection, and cannot determine the precise location of the target. Two-dimensional magnetic total field anomaly data can not only achieve target detection, but also can locate the target. The processing method of magnetic total field anomaly data directly affects the detection efficiency and positioning accuracy of the target.

[0003] Magnetic total field anomaly data often contains a variety of noise, making it particularly difficult to detect abnormal signals. In current practical applications, the signal processing method of one-dimensional orthonormal basis function (OBF) decomposition is mainly used to improve the signal-to-noise ratio of one-dimensional data. However, the research on two-dimensional data anomaly signal detection is still in the academic research stage. At the same time, the method of locating the target from two-dimensional data starts from the original data with low signal-to-noise ratio, which is easily affected by noise, thereby reducing the accuracy of target detection. SUMMARY

[0004] Therefore, the present application provides a magnetic total field anomaly data processing method, device and related equipment to improve the accuracy of target detection.

[0005] To achieve the above-mentioned purpose, the first aspect of the present application provides a magnetic total field anomaly data processing method, comprising:

[0006] Obtaining two-dimensional magnetic total field anomaly data, the recording points of which are uniformly distributed in the detection area;

[0007] Based on the two-dimensional magnetic total field anomaly data, obtaining first horizontal distribution data of the modulus of two-dimensional magnetic total field anomaly on each group of unit orthonormal bases, the unit orthonormal bases being obtained by decomposing the two-dimensional magnetic total field anomaly data generated by the magnetic dipole of the target body at any height;

[0008] Based on the first horizontal distribution data, obtaining second horizontal distribution data of the energy of two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each group of unit orthonormal bases;

[0009] Determine a horizontal coordinate of the target object based on the second horizontal distribution data.

[0010] Preferably, the process of obtaining the two-dimensional magnetic total field anomaly data comprises:

[0011] Using the magnetometer, record the magnetic total field data of the magnetic total field sensor in the detection area along the preset scanning path and record the point position, to obtain a plurality of data items;

[0012] Perform a daily variation correction process on the magnetic total field data in each data item to obtain the magnetic total field anomaly data of each data item;

[0013] Based on the magnetic total field anomaly data of each data item and the record point position, obtain the magnetic total field anomaly data of each grid point on the grid using a preset interpolation method, wherein the grid is uniformly distributed in the detection area;

[0014] The two-dimensional magnetic total field anomaly data is composed of the magnetic total field anomaly data of each grid point and the coordinates of each grid point.

[0015] Preferably, the unit orthonormal bases comprise:

[0016]

[0017] Each unit orthonormal base satisfies the following equation:

[0018]

[0019] Wherein, (x, y, z) is the three-dimensional coordinates of the record point position.

[0020] Preferably, the process of obtaining the first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on each group of unit orthonormal bases based on the two-dimensional magnetic total field anomaly data comprises:

[0021] Using a window of a preset size at a preset height to intercept each group of unit orthonormal bases to obtain a plurality of orthonormal base window functions;

[0022] Using each orthonormal base window function, perform inner product calculation on the two-dimensional magnetic total field anomaly data grid by grid to obtain the first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on each group of unit orthonormal bases.

[0023] Preferably, the process of using each unit orthonormal base window function to perform inner product calculation on the two-dimensional magnetic total field anomaly data grid by grid to obtain the first horizontal distribution data of the modulus of the magnetic total field anomaly on each group of unit orthonormal bases comprises:

[0024] The first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on each group of unit orthonormal bases is calculated using the following equation: n(x, y):

[0025]

[0026] wherein I, J are the size of the window in x and y direction respectively, dΔT(x i ,y j ) is the magnetic total field anomaly data of the grid point at (x i ,y j ).

[0027] Preferably, based on the first horizontal distribution data, the process of obtaining the second horizontal distribution data of the energy of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each set of unit orthogonal bases comprises:

[0028] The energy horizontal distribution data E(x, y) of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each set of orthogonal bases is calculated by using the following equation:

[0029]

[0030] wherein α n (x, y) is the first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on the nth set of orthogonal bases;

[0031] The energy horizontal distribution data is normalized to obtain the second horizontal distribution data of the energy of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each set of unit orthogonal bases.

[0032] Preferably, based on the second horizontal distribution data, the process of determining the horizontal coordinate of the target body comprises:

[0033] Based on the second horizontal distribution data, a target sub-region in which a ferromagnetic target body exists in the detection region is determined, and the value of the second horizontal distribution data of the target sub-region is greater than a preset threshold value;

[0034] The horizontal coordinate of the point with the maximum value of the second horizontal distribution data in each target sub-region is determined as the horizontal coordinate of the target body.

[0035] The second aspect of the present application provides a magnetic total field anomaly data processing device, comprising:

[0036] A data acquisition unit is configured to acquire two-dimensional magnetic total field anomaly data, and the record points of the two-dimensional magnetic total field anomaly data are uniformly distributed in a detection region;

[0037] a modulus calculation unit, configured to obtain first horizontal distribution data of modulus of two-dimensional magnetic total field anomaly on each group of unit orthogonal bases based on the two-dimensional magnetic total field anomaly data, the unit orthogonal bases being obtained by decomposing the two-dimensional magnetic total field anomaly data generated by a magnetic dipole of the target body at an arbitrary height;

[0038] an energy calculation unit, configured to obtain second horizontal distribution data of energy of the two-dimensional magnetic total field anomaly in a two-dimensional space spanned by the unit orthogonal bases based on the first horizontal distribution data;

[0039] a coordinate determination unit, configured to determine horizontal coordinates of the target body based on the second horizontal distribution data.

[0040] The third aspect of the present application provides a magnetic total field anomaly data processing device, comprising a memory and a processor;

[0041] The memory is configured to store a program.

[0042] The processor is configured to execute the program to implement each step of the magnetic total field anomaly data processing method.

[0043] The fourth aspect of the present application provides a storage medium having a computer program stored thereon, the computer program being executed by a processor to implement each step of the magnetic total field anomaly data processing method.

[0044] According to the above technical solution, the two-dimensional magnetic total field anomaly data is first obtained, wherein the recording points of the two-dimensional magnetic total field anomaly data are uniformly distributed in the detection area to form the two-dimensional magnetic total field gradient data in a grid form. Then, the first horizontal distribution data of modulus of the two-dimensional magnetic total field anomaly on each group of unit orthogonal bases is obtained based on the two-dimensional magnetic total field anomaly data, wherein the unit orthogonal bases are obtained by decomposing the two-dimensional magnetic total field anomaly data generated by a magnetic dipole of the target body at an arbitrary height. Next, the second horizontal distribution data of energy of the two-dimensional magnetic total field anomaly in a two-dimensional space spanned by the unit orthogonal bases is obtained based on the first horizontal distribution data. It can be understood that, compared with the first horizontal distribution data, the second horizontal distribution data in the two-dimensional space has a higher signal-to-noise ratio, and the area with relatively strong energy, i.e. the area with higher value of the second horizontal distribution data, means that the ferromagnetic target is more likely to exist. Finally, the horizontal coordinates of the target body are determined based on the second horizontal distribution data. The two-dimensional magnetic total field anomaly data is converted into the two-dimensional energy distribution, which can improve the signal-to-noise ratio of the data, effectively improve the accuracy of target detection, and determine the horizontal coordinates of the target to realize the positioning of the target in the horizontal direction. BRIEF DESCRIPTION OF DRAWINGS

[0045] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the accompanying drawings needed to be used in the embodiments or prior art description will be briefly introduced as follows. Obviously, the accompanying drawings in the following description only need to be drawn for the embodiments of the present application, and for those skilled in the art, other drawings can also be obtained without creative labor based on the provided drawings.

[0046] Figure 1 A schematic diagram of the magnetic total field anomaly data processing method disclosed by the embodiments of the present application is shown in the figure.

[0047] Figure 2 A two-dimensional survey line distribution map disclosed by the embodiments of the present application is shown in the figure.

[0048] Figure 3 A two-dimensional magnetic total field gradient surface map obtained after Kriging interpolation provided by the embodiments of the present application is shown in the figure.

[0049] Figure 4 A two-dimensional unit orthogonal base surface map calculated by the first group of unit orthogonal base functions provided by the embodiments of the present application is shown in the figure.

[0050] Figure 5 A two-dimensional unit orthogonal base surface map calculated by the second group of unit orthogonal base functions provided by the embodiments of the present application is shown in the figure.

[0051] Figure 6 A two-dimensional unit orthogonal base surface map calculated by the third group of unit orthogonal base functions provided by the embodiments of the present application is shown in the figure.

[0052] Figure 7 A two-dimensional unit orthogonal base surface map calculated by the fourth group of unit orthogonal base functions provided by the embodiments of the present application is shown in the figure.

[0053] Figure 8 A two-dimensional unit orthogonal base surface map calculated by the fifth group of unit orthogonal base functions provided by the embodiments of the present application is shown in the figure.

[0054] Figure 9 An energy distribution curve provided by the embodiments of the present application is shown in the figure.

[0055] Figure 10 A schematic diagram of the magnetic total field anomaly data processing device disclosed by the embodiments of the present application is shown in the figure.

[0056] Figure 11 A schematic diagram of the magnetic total field anomaly data processing device disclosed by the embodiments of the present application is shown in the figure. DETAILED DESCRIPTION

[0057] With reference to the accompanying drawings, the technical solutions in the embodiments of the present application will be described clearly and completely. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all the other embodiments obtained by a person of ordinary skill in the art without creative work fall within the protection scope of the present application.

[0058] The magnetic total field anomaly data processing method provided by the embodiments of the present application will be described below. Please refer to Figure 1 The magnetic total field anomaly data processing method provided by the embodiments of the present application can include the following steps:

[0059] In step S101, two-dimensional magnetic total field anomaly data is acquired.

[0060] The two-dimensional magnetic total field anomaly data is obtained by applying a daily variation correction and other data processing methods to the magnetic total field data of a magnetic total field sensor arranged in a detection area, and the recording points of the two-dimensional magnetic total field anomaly data are uniformly distributed in the detection area.

[0061] Specifically, after the detection area is determined, a magnetic total field sensor is arranged in the detection area, a magnetometer mounting platform (a flight platform, a water platform or a handheld type, etc.) moves back and forth along several parallel lines, the magnetometer records the magnetic total field data of the magnetic total field sensor and coordinate information in an interval sampling manner, and finally, the two-dimensional magnetic total field anomaly data is calculated after a series of processing such as daily variation correction.

[0062] Figure 2 Exemplarily, the two-dimensional line distribution diagram as shown in may be used to collect the magnetic total field data of the magnetic total field sensor, wherein the lines are along the north-south direction, the distance between each adjacent measuring point on the line is 0.1 m, the distance between each adjacent line is 1 m, and the height of the magnetic total field sensor from the ground is 0.5 m. In order to verify the effectiveness of the present application, five magnetic anomaly bodies are buried in the detection area, and the spatial coordinates are (6, 7, -0.5), (14, 20, -0.5), (15, 10, -0.5), (5, 16, -0.5) and (8, 25, -0.5), and the magnetic moment is 1 Am 2, magnetic moment inclination 5°, declination -10°; background geomagnetic field inclination 30°, declination 0°. Obtain the total magnetic field data of the total magnetic field sensor at each measuring point, and calculate the two-dimensional total magnetic field anomaly data based on the total magnetic field data. Finally, through interpolation processing, the two-dimensional total magnetic field anomaly data of each recording point uniformly distributed in the detection area is obtained. For the total magnetic field data of the total magnetic field sensor at each measuring point, in the actual application of the present application, these data are collected by a magnetometer or other equipment; in the algorithm verification stage, the magnetic total field anomaly data at each measuring point can be directly calculated by simulation calculation method, and the preset signal-to-noise ratio (such as 0 dB) Gaussian white noise is added to simulate the actual observation data.

[0063] Step S102, based on the two-dimensional total magnetic field anomaly data, obtain the first horizontal distribution data of the modulus of the two-dimensional total magnetic field anomaly on each group of unit orthogonal bases.

[0064] Wherein, each group of unit orthogonal bases is obtained by decomposing the two-dimensional total magnetic field anomaly data generated by the magnetic dipole of the target body at any height; the modulus of the two-dimensional total magnetic field anomaly on a group of unit orthogonal bases is defined as the projection of the two-dimensional total magnetic field anomaly on the group of unit orthogonal bases.

[0065] Step S103, based on the first horizontal distribution data, obtain the second horizontal distribution data of the energy of the two-dimensional total magnetic field anomaly in the two-dimensional space spanned by each group of unit orthogonal bases.

[0066] Wherein, the energy is defined as the square of the modulus value. By taking the square of the modulus value of the two-dimensional total magnetic field anomaly, the differences between the two-dimensional total magnetic field anomalies of each recording point in the detection area are amplified, and the signal-to-noise ratio is improved.

[0067] Step S104, based on the second horizontal distribution data, determine the horizontal coordinates of the target body.

[0068] The application first acquires two-dimensional magnetic total field anomaly data, wherein the record points of the two-dimensional magnetic total field anomaly data are uniformly distributed in a detection area to form grid two-dimensional magnetic total field gradient data. Then, based on the two-dimensional magnetic total field anomaly data, first horizontal distribution data of the modulus of two-dimensional magnetic total field anomaly on each group of unit orthogonal bases is acquired, wherein the unit orthogonal bases are obtained by decomposing the two-dimensional magnetic total field anomaly data generated by the magnetic dipole of the target body at any height. Next, based on the first horizontal distribution data, second horizontal distribution data of the energy of two-dimensional magnetic total field anomaly in a two-dimensional space spanned by each group of unit orthogonal bases is acquired. It can be understood that, compared with the first horizontal distribution data, the second horizontal distribution data in the two-dimensional space has a higher signal-to-noise ratio, and the area with relatively strong energy, that is, the area with higher value of the second horizontal distribution data, means that the ferromagnetic target is more likely to exist. Finally, based on the second horizontal distribution data, the horizontal coordinate of the target body is determined. The application converts the two-dimensional magnetic total field anomaly data into two-dimensional energy distribution, which can improve the signal-to-noise ratio of the data, effectively improve the accuracy of target detection, and determine the horizontal coordinate of the target, thereby realizing the positioning of the target in the horizontal direction.

[0069] In some embodiments of the application, the process of acquiring two-dimensional magnetic total field anomaly data in step S101 can include:

[0070] S1, using a magnetometer, recording the magnetic total field data of the magnetic total field sensor in the detection area along the preset scanning path and the record points to obtain a plurality of data items.

[0071] Each data item includes magnetic total field data and record points. It can be understood that the record points are the detection positions when the data is detected.

[0072] S2, performing diurnal correction processing on the magnetic total field data in each data item to obtain the magnetic total field anomaly data of each data item.

[0073] The purpose of diurnal correction processing is to eliminate the influence of geomagnetic field diurnal variation on the observation data (magnetic total field data). Generally, within a range of 100 km 2 , it can be considered that the diurnal variation is the same, so the diurnal variation data of the detection area can be obtained by self-observation or obtained from a nearby geomagnetic station. In addition, the diurnal variation observation results can be plotted into a diurnal variation curve for reference.

[0074] S3, based on the magnetic total field anomaly data and the record points of each data item, using a preset interpolation method to acquire the magnetic total field anomaly data at each grid point of the grid.

[0075] In actual magnetic detection tasks, the spacing between the measuring points is usually not equal to the spacing between the measuring lines, for example, Figure 2In the shown two-dimensional survey line distribution map, the survey point spacing is much smaller than the survey line spacing, thus, the survey line coverage area needs to be gridded, and the obtained grid is uniformly distributed in the detection area, so as to obtain grid subdivision data with the same size in two horizontal directions (i.e., horizontal transverse and longitudinal directions). Commonly used two-dimensional gridding methods include bilinear interpolation, bicubic interpolation, spline interpolation, minimum curvature interpolation, Kriging interpolation and the like. For two-dimensional magnetic total field anomaly data, the Kriging interpolation method can be used. Figure 3 The two-dimensional magnetic total field anomaly surface map obtained after Kriging interpolation provided by the embodiment of the present application has relatively high random noise due to the existence of large detection interference and errors, and thus the magnetic total field anomaly has a relatively low signal-to-noise ratio.

[0076] At this point, the magnetic total field anomaly data at each grid point (recorded point) and the coordinates of each grid point are obtained, and thus the magnetic total field anomaly data at each grid point and the coordinates of each grid point constitute two-dimensional magnetic total field anomaly data.

[0077] In some embodiments of the present application, the unit orthogonal basis mentioned in the above step S102 can be represented by the following five groups of unit orthogonal basis functions:

[0078]

[0079] Each unit orthogonal basis satisfies the following equation:

[0080]

[0081] where (x, y, z) is the three-dimensional coordinates of the recorded point.

[0082] Exemplarily, Figures 4-8 A two-dimensional unit orthogonal basis surface map calculated by each group of unit orthogonal basis functions is provided.

[0083] In some embodiments of the present application, the process of obtaining the first horizontal distribution data of the modulus of two-dimensional magnetic total field anomaly on each group of unit orthogonal bases based on two-dimensional magnetic total field anomaly data in the above step S102 can include:

[0084] S1, each group of unit orthogonal bases is intercepted by using a window of a preset size at a preset height, to obtain a plurality of orthogonal basis window functions.

[0085] S2, the inner product calculation of two-dimensional magnetic total field anomaly data is performed by using each orthogonal basis window function, and the first horizontal distribution data of the modulus of two-dimensional magnetic total field anomaly on each group of unit orthogonal bases is obtained.

[0086] Since the two-dimensional magnetic total field gradient data is discretely distributed in space, a numerical calculation formula for inner product in the discrete case is needed to calculate the first horizontal distribution data.

[0087] Based on this, in some embodiments of the present application, the process of S2 using each unit orthogonal basis window function to perform inner product calculation on the two-dimensional magnetic total field anomaly data to obtain the first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on each group of unit orthogonal bases can include:

[0088] The first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on each group of unit orthogonal bases is calculated by the following equation n (x,y):

[0089]

[0090] where I and J are the sizes of the window in the x direction and the y direction respectively, dΔT(x i ,y j ) is the magnetic total field anomaly data of the grid point at (x i ,y j ).

[0091] In some embodiments of the present application, the process of S103 obtaining the second horizontal distribution data of the energy of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each group of unit orthogonal bases based on the first horizontal distribution data can include:

[0092] S1, the energy horizontal distribution data E(x,y) of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each group of orthogonal bases is calculated by the following equation:

[0093]

[0094] where α n (x,y) is the first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on the nth group of orthogonal bases.

[0095] S2, the energy horizontal distribution data is normalized to obtain the second horizontal distribution data of the energy of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each group of unit orthogonal bases.

[0096] The normalization process can exclude the influence of the magnetic moment size of the target body, facilitating subsequent data processing. For example, Figure 9The second horizontal distribution data (i.e. the energy distribution graph) is provided, and it can be seen from the graph that the second horizontal distribution data can greatly distinguish the energy difference at different grid points. In the simulation experiment of the present application, the signal-to-noise ratio of the E surface is about 17 dB, which is improved by 17 dB compared with the signal-to-noise ratio of 0 dB of the simulation observation data.

[0097] In some embodiments of the present application, the process of determining the horizontal coordinates of the target body based on the second horizontal distribution data in the above step S104 can include:

[0098] S1, determining a target sub-region in which the ferromagnetic target body exists in the detection region based on the second horizontal distribution data.

[0099] Since the energy horizontal distribution data has been normalized in the foregoing step, the energy value in the detection region is between 0 and 1. According to the overall energy distribution, a suitable threshold value is selected, and the region with energy higher than the threshold value is identified as the target sub-region in which the ferromagnetic target exists. As shown in Figure 9 , assuming that the threshold value is 0.5, five target sub-regions are easily obtained.

[0100] S2, determining the horizontal coordinates of the point with the maximum value of the second horizontal distribution data in each target sub-region as the horizontal coordinates of the target body.

[0101] As shown in the second horizontal distribution data Figure 9 , the horizontal coordinates of the energy maximum points of each target sub-region are read as (6, 7.1), (14, 20), (14.9, 10.1), (5, 16.1), and (8, 25), and it can be seen that the maximum horizontal positioning error is about 0.1 m. In actual application, the positioning accuracy can be further improved by encrypting the grid at the cost of reducing the operation efficiency.

[0102] The magnetic total field anomaly data processing device provided in the embodiments of the present application is described below. The magnetic total field anomaly data processing device described below can be correspondingly referred to the magnetic total field anomaly data processing method described above.

[0103] Referring to Figure 10 , the magnetic total field anomaly data processing device provided in the embodiments of the present application can include:

[0104] The data acquisition unit 21 is configured to acquire two-dimensional magnetic total field anomaly data, and the recording points of the two-dimensional magnetic total field anomaly data are uniformly distributed in the detection region.

[0105] The modulus calculation unit 22 is configured to obtain first horizontal distribution data of modulus of the two-dimensional magnetic total field anomaly on each group of unit orthogonal bases based on the two-dimensional magnetic total field anomaly data, the unit orthogonal bases being obtained by decomposing the two-dimensional magnetic total field anomaly data generated by the magnetic dipole of the target body at an arbitrary height;

[0106] The energy calculation unit 23 is configured to obtain second horizontal distribution data of energy of the two-dimensional magnetic total field anomaly in a two-dimensional space spanned by the unit orthogonal bases based on the first horizontal distribution data;

[0107] The coordinate determination unit 24 is configured to determine the horizontal coordinate of the target body based on the second horizontal distribution data.

[0108] In some embodiments of the present application, the data acquisition unit 21 obtains the two-dimensional magnetic total field anomaly data in the following manner:

[0109] The magnetometer records the magnetic total field data of the magnetic total field sensor in the detection area and the recording point along the preset scanning path, and a plurality of data items are obtained;

[0110] The magnetic total field data in each data item is subjected to a daily variation correction process to obtain the magnetic total field anomaly data of each data item;

[0111] The magnetic total field anomaly data of each grid point on the grid is obtained by using a preset interpolation method based on the magnetic total field anomaly data and the recording point of each data item, the grid being uniformly distributed in the detection area;

[0112] The two-dimensional magnetic total field anomaly data is composed of the magnetic total field anomaly data of each grid point and the coordinates of each grid point.

[0113] In some embodiments of the present application, the modulus calculation unit 22 obtains the first horizontal distribution data of modulus of the two-dimensional magnetic total field anomaly on each group of unit orthogonal bases based on the two-dimensional magnetic total field anomaly data in the following manner:

[0114] Each group of unit orthogonal bases is obtained by using a window of a preset size at a preset height, and a plurality of orthogonal base window functions are obtained;

[0115] The two-dimensional magnetic total field anomaly data is subjected to inner product calculation by using each orthogonal base window function and grid point by grid point to obtain the first horizontal distribution data of modulus of the two-dimensional magnetic total field anomaly on each group of unit orthogonal bases.

[0116] The function expression of the unit orthogonal bases is consistent with the 10 groups of unit orthogonal base functions described above.

[0117] In some embodiments of the present application, the modulus calculation unit 22 utilizes each unit orthogonal base window function to perform inner product calculation on the two-dimensional magnetic total field anomaly data grid by grid, and the process of obtaining the first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on each group of unit orthogonal bases can include:

[0118] The first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on each group of unit orthogonal bases is calculated by using the following equation n (x, y):

[0119]

[0120] wherein I and J are the sizes of the window in the x direction and the y direction respectively, dΔT(x i ,y j ) is the magnetic total field anomaly data of the grid point at (x i ,y j ).

[0121] In some embodiments of the present application, the energy calculation unit 23 obtains the second horizontal distribution data of the energy of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each group of unit orthogonal bases based on the first horizontal distribution data, and the process can include:

[0122] The energy horizontal distribution data E(x, y) of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each group of orthogonal bases is calculated by using the following equation:

[0123]

[0124] wherein α n (x, y) is the first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on the nth group of orthogonal bases;

[0125] The energy horizontal distribution data is normalized to obtain the second horizontal distribution data of the energy of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each group of unit orthogonal bases.

[0126] In some embodiments of the present application, the coordinate determination unit 24 determines the horizontal coordinates of the target body based on the second horizontal distribution data, and the process can include:

[0127] Based on the second horizontal distribution data, a target sub-region in which a ferromagnetic target body exists in the detection region is determined, and the value of the second horizontal distribution data of the target sub-region is greater than a preset threshold value;

[0128] The horizontal coordinates of the point with the maximum value of the second horizontal distribution data in each target sub-region are determined as the horizontal coordinates of the target body.

[0129] The magnetic total field anomaly data processing device provided in the embodiments of the present application can be applied to a magnetic total field anomaly data processing device, such as a computer. Figure 11 A hardware structure block diagram of the magnetic total field anomaly data processing device is shown, referring to Figure 11 The hardware structure of the magnetic total field anomaly data processing device can include at least one processor 31, at least one communication interface 32, at least one memory 33 and at least one communication bus 34.

[0130] In the embodiments of the present application, the number of the processor 31, the communication interface 32, the memory 33 and the communication bus 34 is at least one, and the processor 31, the communication interface 32 and the memory 33 complete the communication with each other through the communication bus 34;

[0131] The processor 31 can be a central processing unit CPU, or a specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement the embodiments of the present application, etc.

[0132] The memory 33 can include a high-speed RAM memory, and can also include a non-volatile memory, such as at least one disk memory.

[0133] The memory 33 stores a program, and the processor 31 can call the program stored in the memory 33, and the program is used for:

[0134] Obtaining two-dimensional magnetic total field anomaly data, and the recording points of the two-dimensional magnetic total field anomaly data are uniformly distributed in a detection area;

[0135] Based on the two-dimensional magnetic total field anomaly data, first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on each group of unit orthogonal bases is obtained, and the unit orthogonal bases are obtained by decomposing the two-dimensional magnetic total field anomaly data generated by the magnetic dipole of the target body at any height;

[0136] Based on the first horizontal distribution data, second horizontal distribution data of the energy of the two-dimensional magnetic total field anomaly in a two-dimensional space spanned by each group of unit orthogonal bases is obtained;

[0137] Based on the second horizontal distribution data, the horizontal coordinates of the target body are determined.

[0138] Optionally, the detailed functions and extended functions of the program can refer to the description above.

[0139] The embodiments of the present application also provide a storage medium, which can store a program suitable for the processor to execute, and the program is used for:

[0140] acquire two-dimensional magnetic total field anomaly data, record points of the two-dimensional magnetic total field anomaly data being uniformly distributed in a detection region;

[0141] based on the two-dimensional magnetic total field anomaly data, acquire first horizontal distribution data of a modulus of two-dimensional magnetic total field anomaly on each group of unit orthogonal bases, the unit orthogonal bases being obtained by decomposing two-dimensional magnetic total field anomaly data generated by a magnetic dipole of a target body at an arbitrary height;

[0142] based on the first horizontal distribution data, acquire second horizontal distribution data of energy of two-dimensional magnetic total field anomaly in a two-dimensional space spanned by the unit orthogonal bases;

[0143] based on the second horizontal distribution data, determine a horizontal coordinate of the target body.

[0144] Optionally, the refinement function and the expansion function of the program can refer to the description above.

[0145] In summary:

[0146] The application first acquires two-dimensional magnetic total field anomaly data, record points of the two-dimensional magnetic total field anomaly data being uniformly distributed in a detection region, forming two-dimensional magnetic total field gradient data in a grid. Then, based on the two-dimensional magnetic total field anomaly data, first horizontal distribution data of a modulus of two-dimensional magnetic total field anomaly on each group of unit orthogonal bases is acquired, the unit orthogonal bases being obtained by decomposing two-dimensional magnetic total field anomaly data generated by a magnetic dipole of a target body at an arbitrary height. Next, based on the first horizontal distribution data, second horizontal distribution data of energy of two-dimensional magnetic total field anomaly in a two-dimensional space spanned by the unit orthogonal bases is acquired. It can be understood that, compared with the first horizontal distribution data, the second horizontal distribution data in the two-dimensional space has a higher signal-to-noise ratio, and a region with relatively strong energy, that is, a region with a higher value of the second horizontal distribution data, means that a ferromagnetic target is more likely to exist. Finally, based on the second horizontal distribution data, a horizontal coordinate of the target body is determined. The application converts two-dimensional magnetic total field anomaly data into two-dimensional energy distribution, which can improve the signal-to-noise ratio of the data, effectively improve the accuracy of target detection, and determine the horizontal coordinate of the target, realizing the positioning of the target in the horizontal direction.

[0147] Finally, it should be noted that the terms "first", "second", and the like, herein do not denote any order, quantity, combination, or importance, but rather are used to distinguish one element from another, and are not intended to denote the presence of any such actual relationship or order. Moreover, the terms "include", "have", or any other variant thereof are intended to encompass non-exclusive inclusions, such that processes, methods, articles, or apparatuses that comprise a list of elements are not required to comprise only those elements in the list, but can include other elements not expressly listed, or also include elements inherent in such processes, methods, articles, or apparatuses. Without additional restrictions, an element preceded by "comprises... a" does not exclude the existence of additional identical elements in the process, method, article, or apparatus that comprises the stated element.

[0148] The various embodiments in the specification are described in progressive order with each embodiment building on one or more of the previous embodiments, however the order of the embodiments described is not intended to be construed as a requirement or limitation for these embodiments. Any one or more of the embodiments described with reference to a particular set of one or more other embodiments are optionally employable together with one or more other embodiments and / or in any appropriate combination.

[0149] The above description of disclosed embodiments is intended to enable those skilled in the art to make or use the application. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein can be applied to other embodiments without departing from the spirit or scope of the application. Thus, the present application is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for processing magnetic total field anomaly data, characterized in that, include: Two-dimensional magnetic total field anomaly data is acquired, wherein the recording points of the two-dimensional magnetic total field anomaly data are evenly distributed in the detection area; Based on the two-dimensional magnetic total field anomaly data, the first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on each set of unit orthogonal basis is obtained. The unit orthogonal basis is obtained by decomposing the two-dimensional magnetic total field anomaly data generated by the magnetic dipole of the target body at any height. Based on the first horizontal distribution data, the second horizontal distribution data of the energy of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each set of unit orthogonal bases is obtained; Based on the second horizontal distribution data, determine the horizontal coordinates of the target body; The orthonormal basis includes: ； Each unit orthogonal basis satisfies the following equation: ； in, The three-dimensional coordinates of the recorded point are given.

2. The method according to claim 1, characterized in that, The process of acquiring two-dimensional magnetic total field anomaly data includes: Using a magnetometer, the magnetic field data of the magnetic field sensor in the detection area and the recording points are recorded along a preset scanning path to obtain multiple data items; Diurnal variation correction is applied to the total magnetic field data in each data item to obtain the total magnetic field anomaly data for each data item; Based on the total magnetic field anomaly data and recording points of each data item, the total magnetic field anomaly data at each grid point of the grid is obtained using a preset interpolation method, wherein the grid is uniformly distributed in the detection area; The two-dimensional magnetic field anomaly data is composed of the magnetic field anomaly data at each grid point and the coordinates of each grid point.

3. The method according to claim 1, characterized in that, The process of obtaining the first level distribution data of the modulus of the two-dimensional magnetic total field anomaly on each set of unit orthogonal bases based on the aforementioned two-dimensional magnetic total field anomaly data includes: Each set of unit orthogonal bases is extracted using a window of a preset size at a preset height, resulting in multiple orthogonal base window functions; Using each orthogonal basis window function, the inner product of the two-dimensional magnetic total field anomaly data is calculated point by point on each grid to obtain the first level distribution data of the modulus of the two-dimensional magnetic total field anomaly on each set of unit orthogonal basis.

4. The method according to claim 3, characterized in that, The process of calculating the inner product of the two-dimensional magnetic total field anomaly data point by point using each unit orthogonal basis window function to obtain the first level distribution data of the modulus of the magnetic total field anomaly on each set of unit orthogonal bases includes: The first level distribution data of the modulus of the two-dimensional total magnetic field anomaly on each set of unit orthogonal bases were calculated using the following equations. : ； Where I and J represent the size of the window in the x and y directions, respectively. In order to be in The total magnetic field anomaly data at the grid points.

5. The method according to claim 4, characterized in that, Based on the first horizontal distribution data, the process of obtaining the second horizontal distribution data of the energy of the two-dimensional total magnetic field anomaly in the two-dimensional space spanned by each set of unit orthogonal bases includes: The energy level distribution data of the two-dimensional total magnetic field anomaly in the two-dimensional space spanned by each set of orthogonal bases were calculated using the following equations. : ； in, The first level distribution data of the modulus of the two-dimensional total magnetic field anomaly on the nth orthogonal basis; The energy level distribution data is normalized to obtain the second level distribution data of the energy of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each set of unit orthogonal bases.

6. The method according to any one of claims 1 to 5, characterized in that, The process of determining the horizontal coordinates of the target object based on the second horizontal distribution data includes: Based on the second horizontal distribution data, a target sub-region containing a ferromagnetic target body is determined in the detection area, wherein the value of the second horizontal distribution data of the target sub-region is greater than a preset threshold. The horizontal coordinates of the point with the largest value in the second horizontal distribution data in each target sub-region are determined as the horizontal coordinates of the target body.

7. A magnetic total field anomaly data processing device, characterized in that, include: The data acquisition unit is used to acquire two-dimensional magnetic field anomaly data, wherein the recording points of the two-dimensional magnetic field anomaly data are evenly distributed in the detection area; The modulus calculation unit is used to obtain the first horizontal distribution data of the modulus of the two-dimensional magnetic total field anomaly on each set of unit orthogonal bases based on the two-dimensional magnetic total field anomaly data. The unit orthogonal bases are obtained by decomposing the two-dimensional magnetic total field anomaly data generated by the magnetic dipole of the target body at any height. The energy calculation unit is used to obtain the second horizontal distribution data of the energy of the two-dimensional magnetic total field anomaly in the two-dimensional space spanned by each set of unit orthogonal bases, based on the first horizontal distribution data. The coordinate determination unit is used to determine the horizontal coordinates of the target body based on the second horizontal distribution data; The orthonormal basis includes: ； Each unit orthogonal basis satisfies the following equation: ； in, The three-dimensional coordinates of the recorded point are given.

8. A magnetic total field anomaly data processing device, characterized in that, include: Memory and processor; The memory is used to store programs; The processor is used to execute the program to implement each step of the magnetic total field anomaly data processing method as described in any one of claims 1 to 6.

9. A storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements each step of the magnetic total field anomaly data processing method as described in any one of claims 1 to 6.

Images