Method for acquiring accurate aviation magnetic measurement data based on total horizontal gradient method

Abstract

The invention discloses a method for acquiring accurate aviation magnetic measurement data based on a total horizontal gradient method. The method comprises the following steps: acquiring aviation magnetic measurement data T actually measured in a test area; the horizontal gradient # imgabs0 # of the aviation magnetic measurement data T in the x direction and the horizontal gradient # imgabs1 # of the aviation magnetic measurement data T in the y direction are calculated, and the original total horizontal gradient Hxy is calculated according to the horizontal gradient # imgabs2 # of the aviation magnetic measurement data T in the x direction and the horizontal gradient # imgabs3 # of the aviation magnetic measurement data T in the y direction; calculating a weak signal enhancement function B according to the original total horizontal gradient Hxy; and according to the original total horizontal gradient Hxy and the weak signal enhancement function B, obtaining accurate aviation magnetic survey data. According to the method for obtaining the accurate aeromagnetic survey data based on the total horizontal gradient method provided by the invention, the boundary resolution of a deeper geologic body can be improved, so that the aeromagnetic data processing and conversion quality and the geologic body boundary recognition effect are improved.

Description

Technical Field

[0001] The present invention relates to the technical field of aeromagnetic survey, and in particular to a method for obtaining accurate aeromagnetic survey data based on a total horizontal gradient method. Background Art

[0002] Airborne magnetic survey is to install an airborne magnetometer system (such as optical pumping, nuclear spinning and fluxgate) on an aircraft to observe the geomagnetic field parameters (such as the total geomagnetic field intensity T or the total magnetic field anomaly △T or its gradient) to find magnetic or magnetic-related ore bodies in order to understand geological structure, conduct magnetic mapping, and solve problems such as urban and engineering stability and archaeology.

[0003] Airborne magnetic survey data is a comprehensive reflection of the magnetic field information of magnetic geological bodies of different depths, shapes, and sizes on the observation surface. However, due to errors in the measurement data or the superposition of magnetic fields, the measurement data is difficult to distinguish, which brings difficulties to geological interpretation.

[0004] At present, the resolution of aeromagnetic anomalies can be improved through the total horizontal gradient method NSTD.

[0005] The advantage of the total horizontal gradient method is that it is less affected by magnetization direction and magnetic anomaly components when detecting geological boundaries in aeromagnetic data than the results of vertical derivative processing. However, its disadvantage is that small-scale linear structures can be easily obscured by larger structures and cannot be identified. In other words, the total horizontal gradient method can only determine the boundaries of shallow geological bodies, but has lower resolution for deeper geological bodies and is prone to generating false geological boundaries in aeromagnetic data, which affects its practical application.

[0006] Therefore, in order to improve the quality of aeromagnetic data processing and conversion and the effect of geological body boundary identification, it is necessary to study and improve the ability of the total horizontal gradient method in actually detecting geological body boundaries using aeromagnetic data.

[0007] The information disclosed in this background technology section is only intended to enhance understanding of the overall background of the invention and should not be regarded as an admission or any form of suggestion that the information constitutes the prior art already known to a person skilled in the art. Summary of the Invention

[0008] To solve the above problems, the purpose of the embodiments of the present invention is to provide a method for obtaining accurate aeromagnetic data based on the total horizontal gradient method to improve the boundary resolution of deeper geological bodies, thereby improving the quality of aeromagnetic data processing and conversion and the effect of geological body boundary identification.

[0009] To achieve the above object, the present invention provides a method for obtaining accurate aeromagnetic data based on the total horizontal gradient method, comprising: obtaining the aeromagnetic data T measured in the test area; calculating the horizontal gradient of the aeromagnetic data T in the x direction. The horizontal gradient of the aeromagnetic data T in the y direction According to the horizontal gradient of the aeromagnetic data T in the x direction And the horizontal gradient of the aeromagnetic data T in the y direction Calculate the original total horizontal gradient H xy According to the original total horizontal gradient H xy Calculate the weak signal enhancement function B; according to the original total horizontal gradient H xy and weak signal enhancement function B to obtain accurate aeromagnetic data.

[0010] In one embodiment, after obtaining the measured aeromagnetic data of the test area, the T calculates the horizontal gradient of the aeromagnetic data T in the x direction. The horizontal gradient of the aeromagnetic data T in the y direction Previously, the method further includes: preprocessing the acquired aeromagnetic data T measured in the test area.

[0011] In one embodiment, the preprocessing includes one or more of coordinate conversion, normal field correction, diurnal variation correction, hysteresis correction, and magnetic field level adjustment.

[0012] In one embodiment, the horizontal gradient of the aeromagnetic data T in the x direction is calculated based on the aeromagnetic data T. The horizontal gradient of the aeromagnetic data T in the y direction The method comprises: calling the geological body characteristic distribution data of the test area through an internet connection; evaluating the test area according to the geological body characteristic distribution data, and hierarchically decomposing the test area according to the evaluation result to obtain M hierarchical scale areas, wherein the M hierarchical scale areas have M hierarchical complexity mean identifiers; decomposing the pre-processed aeromagnetic data T with reference to the M hierarchical scale areas to obtain M hierarchical magnetic data t; determining the M horizontal gradients of the M hierarchical magnetic data t in the x-direction according to the corresponding relationship between the two-dimensional Fourier transform and the Hilbert transform in the frequency domain; and M horizontal gradients in the y direction According to the weighted fusion of the M level complexity mean in the x direction, the M horizontal gradients Output horizontal gradient in the x direction According to the weighted fusion of the M level complexity mean in the y direction, the M horizontal gradient Output horizontal gradient in the Y direction

[0013] In one embodiment, the test area is evaluated based on the geological body characteristic distribution data, and the test area is hierarchically decomposed according to the evaluation results to obtain M hierarchical scale areas, including: dividing the test area into multiple grid areas based on a preset grid scale, using the multiple grid areas to divide the geological body characteristic distribution data to obtain multiple grid geophysical characteristic data; performing terrain complexity evaluation based on the multiple grid geological body characteristic data, and outputting multiple regional terrain complexities of the multiple grid areas; presetting a complexity threshold, and mapping and performing multiple rounds of merging of the multiple grid areas based on whether adjacent differences in the terrain complexities of the multiple areas meet the complexity threshold to obtain the M hierarchical scale areas.

[0014] In one embodiment, the total horizontal gradient H xy and weak signal enhancement function B, to obtain accurate aeromagnetic data including:

[0015] Accurate aeromagnetic data is obtained according to the following formula, which includes:

[0016]

[0017] Among them, Max(H xy ) represents the total horizontal gradient H Xy maximum value, T is the measured aeromagnetic data, x and y are the two directions of the spatial coordinates, Δ is the adjustment parameter, and the value range of Δ is 0-1.

[0018] In one embodiment, the method further comprises: Calculate the output first level mean and first level variance; Arrange the M horizontal gradients in the x direction in ascending order To extract the first horizontal extreme value; and so on, perform M horizontal gradients in the y direction After interactively obtaining the average burial depth of the test area, the average burial depth, the first level average, the first level variance, the first level extreme value, the second level average, the second level variance and the second level extreme value are input into a pre-built adjustment parameter perturbation model for perturbation analysis, and real-time adjustment parameters are output.

[0019] In one embodiment, the method also includes: calling multiple groups of sample disturbance variables and multiple sample adjustment parameters online, wherein each group of sample disturbance variables includes sample burial depth, first sample mean, first sample variance, first sample extreme value, second sample mean, second sample variance and second sample extreme value; associating and storing the multiple groups of sample disturbance variables and multiple sample adjustment parameters based on the knowledge graph to complete the construction of the adjustment parameter disturbance model; inputting the average burial depth, first level mean, first level variance, first level extreme value, second level mean, second level variance and second level extreme value into the adjustment parameter disturbance model, performing data similarity evaluation on the multiple groups of sample disturbance variables based on Euclidean distance ratio, and extracting the sample adjustment parameters corresponding to the similarity extreme values as real-time adjustment parameters.

[0020] To achieve the above object, the present invention also provides a device for obtaining accurate aeromagnetic data based on the total horizontal gradient method, comprising: an acquisition module for obtaining the aeromagnetic data T measured in the test area; a gradient calculation module for calculating the horizontal gradient of the aeromagnetic data T in the x direction. The horizontal gradient of the aeromagnetic data T in the y direction The original total horizontal gradient calculation module is used to calculate the horizontal gradient in the x direction according to the aeromagnetic data T And the horizontal gradient of the aeromagnetic data T in the y direction Calculate the original total horizontal gradient H Xy ; Weak signal enhancement function calculation module, used according to the original total horizontal gradient H x y calculates the weak signal enhancement function B; the precise aeromagnetic data acquisition module is used to obtain the original total horizontal gradient H xy and weak signal enhancement function B to obtain accurate aeromagnetic data.

[0021] To achieve the above objectives, the present invention further provides a storage medium storing computer-executable instructions for executing any one of the above-mentioned methods for obtaining accurate aeromagnetic data based on the total horizontal gradient method.

[0022] Compared with the existing technology, the method, device and storage medium for obtaining accurate aeromagnetic data based on the total horizontal gradient method according to the present invention can better detect the boundaries of the target geological body in the aeromagnetic data, make the boundary recognition results more convergent, and solve the low resolution problem of the existing total horizontal gradient method, thereby improving the quality of aeromagnetic data processing and conversion and the geological body boundary recognition effect.

[0023] In addition, it eliminates the interference of geological body boundaries that produces false aeromagnetic data, enhances the signal-to-noise ratio, and improves the boundary position enhancement and extraction capabilities of geological bodies at different burial depths, with higher resolution and accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0025] Figure 1 A flow chart of a method for obtaining accurate aeromagnetic data based on the total horizontal gradient method provided by an embodiment of the present invention is shown;

[0026] Figure 2 A schematic structural diagram of a device for obtaining accurate aeromagnetic data based on the total horizontal gradient method provided by an embodiment of the present invention is shown. DETAILED DESCRIPTION

[0027] The specific embodiments of the present invention are described in detail below with reference to the accompanying drawings, but it should be understood that the protection scope of the present invention is not limited by the specific embodiments.

[0028] In the description of the present invention, it should be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside", "clockwise", "counterclockwise" and the like to indicate orientations or positional relationships based on the orientations or positional relationships shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be understood as limiting the present invention.

[0029] Unless expressly stated otherwise, throughout the specification and claims, the term "comprise" or variations such as "include" or "comprising", etc., will be understood to include the stated elements or components but not to exclude other elements or other components.

[0030] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of the technical features being referred to. Thus, a feature identified as "first" or "second" may explicitly or implicitly include one or more of the features. In the description of the present invention, "plurality" means two or more, unless otherwise specifically defined.

[0031] Flow chart of the method for obtaining accurate aeromagnetic data based on the total horizontal gradient method provided by the embodiment of the present invention, see Figure 1 , including: step S1-step S5.

[0032] Step 1: Obtain the measured aeromagnetic data T of the test area.

[0033] The test area is analyzed for geological characteristics, such as burial depth, scale, and structural complexity, to determine the terrain complexity of the test area. This determines the distribution of terrain complexity, and then optimizes the flight trajectory based on this terrain complexity distribution, generating a dynamic flight trajectory. Raw magnetic data is then collected based on this dynamic flight trajectory. The spacing and altitude of the flight curves vary with terrain complexity. The aeromagnetic data, or T, or aeromagnetic anomaly field, is the additional magnetic field generated by ferromagnetic geological bodies in the Earth's crust under the influence of the Earth's magnetic field.

[0034] In one implementation, after step S1 and before step S2, the method may further include preprocessing the acquired aeromagnetic data T measured in the test area to obtain preprocessed aeromagnetic data T measured in the test area. The preprocessing may include one or more of coordinate conversion, normal field correction, diurnal variation correction, hysteresis correction, and magnetic field level adjustment. Preprocessing can eliminate noise and errors in the aeromagnetic field measurement data caused by various factors.

[0035] In one implementation, the selection order of pre-processing may be, for example, coordinate transformation→normal field correction→diurnal variation correction→lag correction.

[0036] Furthermore, the above methods can be combined. For example, coordinate transformation and diurnal variation correction can be processed simultaneously in two containers. The registered spatial data can then be synchronized with the time-corrected magnetic field data to unify timestamps and location tags. This improves efficiency by parallelizing specific preprocessing steps and aligning and overlaying the results.

[0037] Step 2: Calculate the horizontal gradient of the aeromagnetic data T in the x direction The horizontal gradient of the aeromagnetic data T in the y direction

[0038] In one implementation, step S2 may further include:

[0039] According to the correspondence between the two-dimensional Fourier transform and the Hilbert transform in the frequency domain, the horizontal gradient of the pre-processed aeromagnetic data T in the x direction is determined. and the horizontal gradient of the preprocessed aeromagnetic data T in the y direction.

[0040] Specifically, the original total horizontal gradient H can be calculated according to the following formula: xy ;

[0041]

[0042] Specifically, the two-dimensional Hilbert transform can be calculated directly in frequency, and its Fourier domain transform factor is:

[0043] DH(u,v)=-i·sign(u,v),

[0044]

[0045] Among them, i 2 = -1, u and v are the circular wave numbers in the x and y directions. The components of the Hilbert transform in the x and y directions can be expressed as:

[0046]

[0047] Among them, F[], F -1 [] represent Fourier transform and inverse transform respectively; H x With H y Represent the components of the two-dimensional Hilbert transform in the x and y directions respectively.

[0048] Therefore, the horizontal gradient in the x direction Horizontal gradient in the y direction It can be calculated using the Hilbert transform, and the two have an equivalent relationship:

[0049]

[0050] By leveraging the correspondence between the Hilbert and Fourier transforms, the Hilbert transform is used directly instead of the Fourier transform to calculate the gradients in the x, y, and z directions. Applying this equivalent relationship, i.e., replacing the Fourier transform with the Hilbert transform, has the advantage of not amplifying noise interference in the aeromagnetic data.

[0051] In one implementation, step S2 may further include:

[0052] The geological body characteristic distribution data of the test area is called online; the test area is evaluated according to the geological body characteristic distribution data, and the test area is hierarchically decomposed according to the evaluation result to obtain M hierarchical scale areas, wherein the M hierarchical scale areas have M hierarchical complexity mean identifiers; the pre-processed aeromagnetic data T is decomposed with reference to the M hierarchical scale areas to obtain M hierarchical magnetic data t; according to the corresponding relationship between the two-dimensional Fourier transform and the Hilbert transform in the frequency domain, the horizontal gradient of the M hierarchical magnetic data t in the x direction is determined. and M horizontal gradients in the y direction According to the weighted fusion of the M level complexity mean in the x direction, the M horizontal gradients Output horizontal gradient in the x direction According to the weighted fusion of the M level complexity mean in the y direction, the M horizontal gradient Output horizontal gradient in the Y direction

[0053] In one implementation, the test area is evaluated based on the geological body characteristic distribution data, and the test area is hierarchically decomposed according to the evaluation results to obtain M hierarchical scale areas, including:

[0054] After dividing the test area into multiple grid areas based on a preset grid scale, the multiple grid areas are used to divide the geological body characteristic distribution data to obtain multiple grid geological body characteristic data; terrain complexity evaluation is performed based on the multiple grid geological body characteristic data, and multiple regional terrain complexities of the multiple grid areas are output; a complexity threshold is preset, and based on whether adjacent differences in the terrain complexities of the multiple areas meet the complexity threshold, the multiple grid areas are mapped and merged multiple times to obtain the M hierarchical scale areas.

[0055] It should be understood that in complex geological structure scenarios, the gradient calculation at a single scale may have the problem of excessively high gradient values in strong magnetic anomaly areas, resulting in the neglect of deep weak gradients. Based on this, this embodiment decomposes the magnetic survey data into multiple scales to obtain geological body signals corresponding to different burial depths, and then performs the above horizontal gradient calculation on each scale layer to obtain the horizontal gradient of each scale layer, and further uses weight superposition to perform multi-scale gradient fusion.

[0056] Specifically, in this embodiment, geological body characteristic data covering the test area are called from the geological information platform. The geological body characteristic data include but are not limited to geological structure density information, lithology distribution information, burial depth characteristic information and historical geomagnetic data. The geological body characteristic distribution data provides basic geological background support for subsequent regional evaluation and hierarchical decomposition.

[0057] The test area is divided into multiple grid areas according to a preset grid scale (such as 1 km × 1 km), and then the obtained geological body characteristic distribution data are divided according to the grid area division to obtain multiple grid geophysical characteristic data.

[0058] The terrain complexity of each grid is calculated by inputting the geophysical characteristic data of the multiple grids into a quantitative model. The evaluation indicators of the quantitative model include structural density (based on the distribution density of fault zones or folds), lithologic complexity (reflecting the distribution heterogeneity of different lithologies), burial depth variation coefficient (characterizing the spatial variation degree of the burial depth of the geological body), and magnetic anomaly gradient (the first-order derivative of the magnetic field intensity is calculated in combination with historical geomagnetic data to describe the boundary characteristics of the magnetic body).

[0059] A complexity threshold is set as the criterion for merging adjacent grid regions. If the complexity difference between adjacent grid regions is below the threshold, they are merged into the same hierarchical scale region; otherwise, they remain as independent regions. Multiple rounds of merging algorithms dynamically iterate and adjust the adjacent grid regions, ultimately obtaining M hierarchical scale regions.

[0060] According to the grid area composition of each hierarchical scale area, several corresponding terrain complexities are extracted, the average is calculated, and the average is marked in the corresponding hierarchical scale area.

[0061] Further, referring to the M-level scale regional decomposition preprocessed aeromagnetic data T, M-level magnetic survey data t are obtained; according to the correspondence between the two-dimensional Fourier transform and the Hilbert transform in the frequency domain, the M-level horizontal gradient of the M-level magnetic survey data t in the x direction is determined. and M horizontal gradients in the y direction According to the weighted fusion of the M level complexity mean in the x direction, the M horizontal gradients Output horizontal gradient in the x direction According to the weighted fusion of the M level complexity mean in the y direction, the M horizontal gradient Output horizontal gradient in the Y direction

[0062] In step 1, the horizontal gradients M in the x direction of the M level magnetic measurement data t are determined according to the corresponding relationship between the two-dimensional Fourier transform and the Hilbert transform in the frequency domain. and M horizontal gradients in the y direction In the frequency domain, the corresponding relationship between the two-dimensional Fourier transform and the Hilbert transform can include:

[0063]

[0064] Then, the M horizontal gradients in the x direction are weighted and fused according to the mean value of the complexity of the M levels. Output x-direction horizontal gradient According to the weighted fusion of the M level complexity mean in the y direction, the M horizontal gradient Output y-direction horizontal gradient

[0065] Step 3: According to the horizontal gradient of the aeromagnetic data T in the x direction And the horizontal gradient of the aeromagnetic data T in the y direction Calculate the original total horizontal gradient H xy .

[0066] Among them, the horizontal gradient It can be based on the weighted fusion of the M level complexity mean in the x direction horizontal gradient M The x-direction horizontal gradient output later Horizontal gradient It can be based on the weighted fusion of the M level complexity mean in the y direction. The horizontal gradient in the y direction outputted later

[0067] In this embodiment, the M weights of the M-level scale regions are calculated based on the M-level complexity mean values, which are used for weighted fusion calculation of the x-direction horizontal gradient. And the horizontal gradient in the y direction

[0068] Specifically, the original total horizontal gradient H is calculated according to the following formula: xy ;

[0069]

[0070] Step 4: According to the original total horizontal gradient H xy Calculate the weak signal enhancement function.

[0071] Specifically, step 4 can be calculated using the following formula:

[0072] The weak signal enhancement function B is used to improve the resolution of boundary recognition, and its expression is:

[0073]

[0074] It should be noted that the original total horizontal gradient H xy The weak signal enhancement function can be calculated based on the gradient after weighted fusion according to the mean values of the complexity of the M levels, or can be directly calculated based on the gradient without fusion.

[0075] Step 5: According to the original total horizontal gradient H xy and weak signal enhancement function B to obtain accurate aeromagnetic data.

[0076] Specifically, accurate aeromagnetic data can be obtained according to the following formula, which includes:

[0077]

[0078] Among them, Max(H xy ) represents the total horizontal gradient H xy The maximum value, T is the measured aeromagnetic data, x and y are the two directions of the spatial coordinates, Δ represents the adjustment parameter to prevent the denominator from being 0 and to avoid the "analytic singularity" in the expression. The value of Δ is 0-1.

[0079] In step 5, the weighted fusion gradient values obtained through M hierarchical scale areas are used for calculation, which can make the magnetic survey data more stable in complex geological structure scenarios.

[0080] In one implementation, the adjustment parameters can be obtained by: Calculate the output first level mean and first level variance; Arrange the M horizontal gradients in the x direction in ascending order To extract the first horizontal extreme value; and so on, perform M horizontal gradients in the y direction After interactively obtaining the average burial depth of the test area, the average burial depth, the first level average, the first level variance, the first level extreme value, the second level average, the second level variance and the second level extreme value are input into a pre-built adjustment parameter perturbation model for perturbation analysis, and real-time adjustment parameters are output.

[0081] In one implementation, multiple groups of sample disturbance variables and multiple sample adjustment parameters are called online, wherein each group of sample disturbance variables includes sample burial depth, first sample mean, first sample variance, first sample extreme value, second sample mean, second sample variance and second sample extreme value; the multiple groups of sample disturbance variables and multiple sample adjustment parameters are stored in association with each other based on the knowledge graph to complete the construction of the adjustment parameter disturbance model; the average burial depth, first level mean, first level variance, first level extreme value, second level mean, second level variance and second level extreme value are input into the adjustment parameter disturbance model, and the data similarity of the multiple groups of sample disturbance variables is evaluated based on the Euclidean distance ratio, and the sample adjustment parameters corresponding to the similarity extreme values are extracted as real-time adjustment parameters.

[0082] Exemplarily, part of the knowledge graph is stored as the following three groups of sample perturbation variables and corresponding adjustment parameters Δ.

[0083]

[0084] Therefore, the method for acquiring accurate aeromagnetic data based on the total horizontal gradient method provided in this embodiment can better detect the boundaries of target geological bodies in aeromagnetic data, resulting in more convergent boundary identification results. This solves the low resolution problem of the existing total horizontal gradient method, thereby improving the quality of aeromagnetic data processing and conversion and the effectiveness of geological body boundary identification. Furthermore, it eliminates false geological body boundary interference in aeromagnetic data, enhances the signal-to-noise ratio, and improves the ability to enhance and extract the boundaries of geological bodies at different burial depths, achieving higher resolution and accuracy.

[0085] As a result, the ability to enhance and extract the boundary positions of geological bodies at different burial depths is improved, with higher resolution and accuracy.

[0086] The total horizontal gradient method (THDR) provided in this embodiment is expressed as:

[0087]

[0088] Compared with the existing total horizontal gradient method expression:

[0089]

[0090] The total horizontal gradient method expression provided in this embodiment is based on the fact that the total horizontal gradient method is insufficient for enhancing the boundary of deep-buried field sources. It is proposed to use the maximum value of the total horizontal gradient to enhance H xy B. Normalization processing is performed to balance strong and weak aeromagnetic anomalies by limiting the amplitude of aeromagnetic anomalies with large differences in strength to a specific range, avoiding distortion caused by excessively strong or weak magnetic anomaly signals, thereby improving the stability and noise resistance of the magnetic anomaly signal.

[0091] The total horizontal gradient method expression provided in this embodiment introduces a weak signal enhancement function B and the Hilbert transform to establish a reasonable equalization filter, redefining the total horizontal gradient method. Using this method, the position of geological targets in aeromagnetic data can be determined, thereby accurately inferring the boundaries, depth, occurrence, scale, field distribution patterns, and physical properties of the structural field source. This is of great significance for demarcating tectonic units, conducting structural zoning, determining the location of fault zones, distinguishing the distribution of different lithologies and strata, and conducting physical property mapping.

[0092] An embodiment of the present invention further provides a device for acquiring precise aeromagnetic data based on the total horizontal gradient method, comprising: an acquisition module 1, a gradient calculation module 2, an original total horizontal gradient calculation module 3, a weak signal enhancement function calculation module 4, and a precise aeromagnetic data acquisition module 5.

[0093] The acquisition module 1 is used to obtain the measured aeromagnetic data T of the test area. The gradient calculation module 2 is used to calculate the horizontal gradient of the aeromagnetic data T in the x direction based on the aeromagnetic data T. The horizontal gradient of the aeromagnetic data T in the y direction and the vertical gradient of the aeromagnetic data T in the z direction The original total horizontal gradient calculation module 3 is used to calculate the horizontal gradient in the x direction according to the aeromagnetic data T. And the horizontal gradient of the aeromagnetic data T in the y direction Calculate the original total horizontal gradient H xyThe weak signal enhancement function calculation module 4 is used to calculate the original total horizontal gradient H xy Calculate the weak signal enhancement function B. The precise aeromagnetic data acquisition module 5 is used to obtain the original total horizontal gradient H xy and weak signal enhancement function B to obtain accurate aeromagnetic data.

[0094] An embodiment of the present invention further provides a storage medium storing computer-executable instructions, which includes a program for executing the above-mentioned method for obtaining accurate aeromagnetic data based on the total horizontal gradient method. The computer-executable instructions can execute the method in any of the above-mentioned method embodiments.

[0095] 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 (such as floppy disks, hard disks, magnetic tapes, magneto-optical disks (MO), etc.), optical storage (such as CDs, DVDs, BDs, HVDs, etc.), and semiconductor storage (such as ROMs, EPROMs, EEPROMs, non-volatile memories (NAND FLASH), solid-state drives (SSDs), etc.).

[0096] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any modifications or substitutions that can be easily conceived by a person skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be based on the scope of protection of the claims.

Claims

1. A method for obtaining accurate aeromagnetic data based on the total horizontal gradient method, characterized in that: include: Obtain the measured aeromagnetic data T of the test area; Calculate the horizontal gradient of the aeromagnetic data T in the x direction The horizontal gradient of the aeromagnetic data T in the y direction According to the horizontal gradient of the aeromagnetic data T in the x direction And the horizontal gradient of the aeromagnetic data T in the y direction Calculate the original total horizontal gradient H xy ; According to the original total horizontal gradient H xy Calculate the weak signal enhancement function B; According to the original total horizontal gradient H xy and weak signal enhancement function B to obtain accurate aeromagnetic data.

2. The method according to claim 1, characterized in that After obtaining the measured aeromagnetic data T of the test area, the horizontal gradient of the aeromagnetic data T in the x direction is calculated. The horizontal gradient of the aeromagnetic data T in the y direction Previously, the method further includes: preprocessing the acquired aeromagnetic data T measured in the test area.

3. The method according to claim 2, characterized in that The preprocessing includes one or more of coordinate conversion, normal field correction, diurnal variation correction, hysteresis correction and magnetic field level adjustment.

4. The method according to claim 2, characterized in that Calculating the horizontal gradient of the aeromagnetic data T in the x direction The horizontal gradient of the aeromagnetic data T in the y direction include: Retrieving the geological body characteristic distribution data of the test area through the Internet; The test area is evaluated according to the geological body characteristic distribution data, and the test area is hierarchically decomposed according to the evaluation result to obtain M hierarchical scale areas, wherein the M hierarchical scale areas have M hierarchical complexity mean identifiers; Referring to the M-level scale regional decomposition preprocessed aeromagnetic data T, M-level magnetic survey data t are obtained; According to the correspondence between the two-dimensional Fourier transform and the Hilbert transform in the frequency domain, the horizontal gradient M of the M-level magnetic measurement data t in the x direction is determined. and M horizontal gradients in the y direction According to the weighted fusion of the M level complexity mean in the x direction, the M horizontal gradients Output horizontal gradient in x direction According to the weighted fusion of the M level complexity mean in the y direction, the M horizontal gradient Output horizontal gradient in the Y direction 5. The method according to claim 4, characterized in that The test area is evaluated based on the geological body characteristic distribution data, and the test area is hierarchically decomposed according to the evaluation results to obtain M hierarchical scale areas including: After dividing the test area into a plurality of grid areas based on a preset grid scale, the plurality of grid areas are used to divide the geological body characteristic distribution data to obtain a plurality of grid geological body characteristic data; Performing terrain complexity evaluation based on the plurality of grid geological body feature data, and outputting a plurality of regional terrain complexities of the plurality of grid regions; A complexity threshold is preset, and according to whether adjacent differences in the terrain complexity of the multiple regions meet the complexity threshold, multiple rounds of merging of the multiple grid regions are mapped to obtain the M hierarchical scale regions.

6. The method according to claim 1 or 4, characterized in that According to the original total horizontal gradient H xy and weak signal enhancement function B, to obtain accurate aeromagnetic data including: Accurate aeromagnetic data is obtained according to the following formula, which includes: Among them, Max(H xy ) represents the total horizontal gradient H xy maximum value, T is the measured aeromagnetic data, x and y are the two directions of the spatial coordinates, Δ is the adjustment parameter, and the value range of Δ is 0-1.

7. The method according to claim 6, characterized in that Also includes: According to the M horizontal gradients in the x direction Calculate the output first level mean and first level variance; Arrange M horizontal gradients in the x direction in ascending order To extract the first level extreme value; Similarly, perform M horizontal gradients in the y direction. Data processing, outputting the second level mean, second level variance and second level extreme value; After interactively obtaining the average burial depth of the test area, the average burial depth, first level average, first level variance, first level extreme value, second level average, second level variance and second level extreme value are input into a pre-built adjustment parameter perturbation model for perturbation analysis, and real-time adjustment parameters are output.

8. The method according to claim 7, characterized in that Also includes: Calling multiple groups of sample disturbance variables and multiple sample adjustment parameters online, wherein each group of sample disturbance variables includes a sample burial depth, a first sample mean, a first sample variance, a first sample extreme value, a second sample mean, a second sample variance, and a second sample extreme value; Associatively storing the multiple groups of sample disturbance variables and the multiple sample adjustment parameters based on the knowledge graph, and completing the construction of the adjustment parameter disturbance model; The average burial depth, first level mean, first level variance, first level extreme value, second level mean, second level variance and second level extreme value are input into the adjustment parameter disturbance model, and the data similarity of the multiple groups of sample disturbance variables is evaluated based on the Euclidean distance ratio, and the sample adjustment parameters corresponding to the similarity extreme values are extracted as real-time adjustment parameters.

9. A device for obtaining accurate aeromagnetic data based on the total horizontal gradient method, characterized in that: include: An acquisition module is used to obtain the measured aeromagnetic data T of the test area; Gradient calculation module, used to calculate the horizontal gradient of the aeromagnetic data T in the x direction The horizontal gradient of the aeromagnetic data T in the y direction The original total horizontal gradient calculation module is used to calculate the horizontal gradient in the x direction according to the aeromagnetic data T And the horizontal gradient of the aeromagnetic data T in the y direction Calculate the original total horizontal gradient H xy ; The weak signal enhancement function calculation module is used to calculate the signal intensity of the original total horizontal gradient H xy Calculate the weak signal enhancement function B; Accurate aeromagnetic data acquisition module is used to obtain the original total horizontal gradient H xy and weak signal enhancement function B to obtain accurate aeromagnetic data.

10. A storage medium, characterized in that: The storage medium stores computer-executable instructions, and the computer-executable instructions are used to execute the method for obtaining accurate aeromagnetic data based on the total horizontal gradient method according to any one of claims 1 to 8.

Images