A method and system for identifying the spatial position of a concealed magnetic rock mass based on synchronous extraction W transformation, an electronic device, and a storage medium

Abstract

The application discloses a kind of concealed magnetic rock mass spatial position identification method, system based on synchronous extraction W transform, electronic equipment and storage medium, it is related to the automatic interpretation technical field of aerial magnetic data, obtains grid aerial magnetic data;Setting W transform related parameters, W transform is carried out to single channel data to obtain empty wave coefficient, according to this, calculate instantaneous empty wave number, and utilize Dirac function to construct synchronous extraction operator, obtain synchronous extraction W transform coefficient matrix by matrix operation.Complete whole work area aerial magnetic data processing and integration coefficient, slice along wave number direction, in combination with the position and amplitude variation of energy group in slice, identify the spatial position and buried depth of concealed magnetic rock mass.The present application combines W transform and synchronous extraction algorithm, and gives consideration to spatial and wave number resolution, and the algorithm has reversibility, and can completely reconstruct the original signal.The present application is stable in operation, high in recognition accuracy, can efficiently analyze aerial magnetic anomaly, provides effective technical means for concealed magnetic rock mass geological exploration and data interpretation.

Description

Technical Field

[0001] This invention relates to the field of automatic interpretation technology of airborne magnetic data, and more specifically to a method, system, electronic device and storage medium for spatial location identification of concealed magnetic rock masses based on synchronous extraction of W-transform. Background Technology

[0002] Currently, aeromagnetic anomalies refer to the differences between Earth's magnetic field data obtained through airborne magnetic measurements and the theoretical normal magnetic field (reference field), reflecting the magnetic characteristics of underground rocks or ore bodies. Since magnetic field data is a superposition of potential fields from geological bodies of different depths and scales, high wavenumber information from deep magnetic bodies is largely attenuated, primarily retaining low wavenumber information. In contrast, shallow magnetic bodies retain not only low wavenumber information but also abundant high wavenumber information. Therefore, aeromagnetic data can be analyzed in the wavenumber domain to obtain information about underground magnetic anomalies.

[0003] Traditional airwave analysis methods have significant limitations. They suffer from insufficient spatial resolution in the low wavenumber range, making it difficult to accurately analyze the low wavenumber characteristics of aeromagnetic data and effectively distinguish concealed magnetic rock masses at different depths and ranges. Existing conventional data processing methods are inefficient in identifying and separating aeromagnetic anomalies, have weak anti-interference capabilities, are prone to losing original signal information during processing, and cannot simultaneously consider spatial and wavenumber resolution. Ultimately, this leads to significant deviations in the location and depth determination of concealed rock masses, failing to meet the application requirements of high-precision geological exploration.

[0004] Therefore, how to provide an identification method to improve the efficiency of identifying and separating aeromagnetic anomalies is a problem that urgently needs to be solved by those skilled in the art. Summary of the Invention

[0005] In view of this, the present invention provides a method, system, electronic device and storage medium for spatial location identification of concealed magnetic rock masses based on synchronous extraction of W transformation, so as to solve the problems existing in the background art.

[0006] To achieve the above objectives, the present invention adopts the following technical solution: A method for spatial location identification of concealed magnetic rock masses based on synchronous extraction of W-transform includes: Acquire gridded aeromagnetic data; Set the basic parameters of the W transform, perform the W transform on one of the data in the gridded aeromagnetic data, and obtain the W transform air wave coefficient; The instantaneous space wave number is calculated based on the W-transform space wave coefficient, and the synchronization extraction operator is defined using the Dirac function. The distribution of the synchronization extraction operator value is calculated based on the instantaneous space wave number. Perform a matrix point-to-point multiplication operation between the W-transform spatial wave coefficients and the synchronous extraction operator to obtain the synchronous extraction W-transform spatial wave coefficient matrix. Each data point in the gridded aeromagnetic data within the work area is processed to obtain the synchronously extracted W-transform air wave coefficients for each aeromagnetic data point, and the obtained synchronously extracted W-transform coefficients are stored in the same matrix. By slicing the matrix along the wavenumber direction, aeromagnetic data information under different wavenumber distributions is obtained. Based on the changes in the position and amplitude of energy clusters on different wavenumber slices, the location and burial depth of concealed magnetic rock masses are identified.

[0007] Optionally, the acquisition of gridded aeromagnetic data includes the following steps: acquiring raw aeromagnetic data, performing pole conversion calculations on the raw aeromagnetic data to obtain aeromagnetic pole conversion data; and regularly gridding the aeromagnetic pole conversion data to obtain gridded aeromagnetic data in text-readable grd format.

[0008] Optionally, the basic parameters of the W transform include window function parameters, wavenumber range parameters, and calculation control parameters.

[0009] Optionally, the expression for the W-transform spatial wave coefficient is:

[0010] in, x For spatial location, k For space wavenumber, This is single-track aeromagnetic data. ≠0 represents the dominant wavenumber of the waveform. K is a scalar factor, and is the initial width of the Gaussian window function.

[0011] Optionally, the expression for the synchronous extraction operator is:

[0012] In the formula, For the Dirac function, It represents the instantaneous wave number.

[0013] A spatial location identification system for concealed magnetic rock masses based on synchronous extraction of W-transform, comprising: The acquisition module acquires gridded aeromagnetic data; The parameter setting module allows you to set the basic parameters for the W transform. The W-transform module performs a W-transform on one of the data streams of the gridded aeromagnetic data to obtain the W-transformed air wave coefficients. The instantaneous space wavenumber calculation module calculates the instantaneous space wavenumber based on the W-transform space wave coefficient. The synchronous extraction operator calculation module defines the synchronous extraction operator using the Dirac function and calculates the distribution of synchronous extraction operator values ​​based on the instantaneous wavenumber. The matrix operation module performs a matrix point-to-point multiplication operation on the W-transform spatial wave coefficients and the synchronous extraction operator to obtain the synchronous extraction W-transform spatial wave coefficient matrix. The aeromagnetic data processing module processes each data point in the gridded aeromagnetic data within the work area to obtain the synchronously extracted W-transform air wave coefficients for each aeromagnetic data point, and stores the obtained synchronously extracted W-transform coefficients in the same matrix. The aeromagnetic anomaly identification module slices the matrix along the wavenumber direction to obtain aeromagnetic data information under different wavenumber distributions. Based on the changes in the position and amplitude of energy clusters on different wavenumber slices, it identifies the location and burial depth of concealed magnetic rock masses.

[0014] An electronic device includes a processor and a memory, the memory storing a computer program, the processor executing the computer program stored in the memory to implement the steps of a method for identifying the spatial location of concealed magnetic rock masses based on synchronous extraction W-transformation.

[0015] A storage medium storing a computer program, which, when executed by a processor, implements the steps of a method for identifying the spatial location of concealed magnetic rock masses based on synchronous extraction of W-transform.

[0016] As can be seen from the above technical solution, compared with the prior art, this invention discloses a method, system, electronic device, and storage medium for spatial location identification of concealed magnetic rock masses based on synchronous extraction of W-transform, which acquires gridded aeromagnetic data; sets basic parameters of W-transform, including window function parameters, wavenumber range parameters, and calculation control parameters; and performs W-transform on one channel of the gridded aeromagnetic data to obtain the W-transformed airwave coefficient. The invention involves: differentiating the W-transform spatial wavenumbers with respect to space to obtain the instantaneous spatial wavenumbers; defining a synchronization extraction operator using the Dirac function and calculating the distribution of the operator values ​​based on the instantaneous spatial wavenumbers; performing matrix point-to-point multiplication of the W-transform spatial wavenumbers with the synchronization extraction operator to obtain the synchronized extraction W-transform spatial wavenumber matrix; processing each track of the gridded aeromagnetic data within the work area to obtain the synchronized extraction W-transform spatial wavenumbers for each track, and storing these coefficients in the same matrix; slicing the matrix along the wavenumber direction to obtain aeromagnetic data information under different wavenumber distributions; and identifying the location and burial depth of concealed magnetic rock masses based on the position and amplitude variations of energy clusters on different wavenumber slices. This invention converts aeromagnetic data from the spatial-amplitude domain to the spatial-wavenumber domain, effectively extracting the wavenumber characteristics of aeromagnetic anomalies during this process. The range of energy clusters on the wavenumber slices can be used to delineate the extent of magnetic rock masses. Furthermore, the distribution of aeromagnetic anomaly energy clusters along the wavenumber can effectively determine the relative burial depth between magnetic rock masses. This method addresses the issue of insufficient resolution in the original data. It also exhibits good stability and anti-interference capabilities, and can perform inverse transformations on the synchronously extracted and transformed results, demonstrating reversibility. The prediction results can be provided to geological interpreters as an objective basis for interpreting aeromagnetic anomalies. Attached Figure Description

[0017] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0018] Figure 1 A flowchart of the method provided by the present invention; Figure 2 This is a schematic diagram of the aeromagnetic polarization data after polarization provided by the present invention; Figure 3a The present invention provides a wavenumber of 0.10m. -1 Below are the output results based on different frequency slices of SEWT aeromagnetic anomalies; Figure 3b The present invention provides a wavenumber of 0.50m -1 Below are the output results based on different frequency slices of SEWT aeromagnetic anomalies; Figure 3c The present invention provides a wavenumber of 0.80m -1 Below are the output results based on different frequency slices of SEWT aeromagnetic anomalies; Figure 3d The present invention provides a wavenumber of 1.00m. -1Below are the output results based on different frequency slices of SEWT aeromagnetic anomalies; Figure 3e The present invention provides a wavenumber of 1.40m. -1 Below are the output results based on different frequency slices of SEWT aeromagnetic anomalies; Figure 3f The present invention provides a wavenumber of 2.00m. -1 The following are the output results based on different frequency slices of SEWT aeromagnetic anomalies. Detailed Implementation

[0019] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0020] This invention combines the W-transform with the Synchronous Extraction Transform (SET) to introduce the Synchronous Extraction W-transform (SEWT). Based on the SEWT method, aeromagnetic data is processed, transforming it from the spatial amplitude domain to the spatial wavenumber domain. This yields wavenumber slices of aeromagnetic data under different wavenumber distributions in the wavenumber domain. Based on the position and amplitude variations of energy clusters on these wavenumber slices, the location and burial depth of concealed magnetic rock masses are inferred, thus completing the spatial location identification of concealed magnetic rock masses.

[0021] First, the aeromagnetic data of the target work area is cropped into a rectangular grid, i.e., regular gridding, and then the data is saved as a text-readable grd file (grid file). Next, the W-transform is used to process the aeromagnetic data of a specific channel in the data matrix to obtain the airwave coefficients. Then, the spatial derivative of the W-transform airwave coefficients is calculated to obtain the instantaneous airwave number distribution of the aeromagnetic data for that channel. Subsequently, a synchronization extraction operator is calculated based on the instantaneous airwave number, and the airwave coefficients at the instantaneous airwave number positions are extracted using the synchronization extraction operator. This removes divergent airwave coefficients, retaining only the airwave coefficients at the instantaneous airwave number positions, thus completing the directional extraction of the W-transform airwave coefficients and obtaining a two-dimensional synchronized extracted airwave coefficient matrix for that channel. An airwave spectrum is then generated based on the coefficient matrix. After performing synchronized extraction and W-transform processing on all data in the work area using the same method, all synchronized compression transformation values ​​are finally stored in the same three-dimensional matrix. The coefficient matrix is ​​plotted as an airwave spectrum, and then sliced ​​along the wavenumber direction of the three-dimensional matrix for extraction and observation, thereby obtaining aeromagnetic data information at different wavenumbers. Based on the changes in the position and amplitude of energy clusters on slices with different wavenumbers, the location and burial depth of the concealed magnetic rock mass can be inferred.

[0022] This invention discloses a method for spatial location identification of concealed magnetic rock masses based on synchronous extraction of W-transform, such as... Figure 1 As shown, it includes: S1. Acquire gridded aeromagnetic data; Gridded aeromagnetic data is typically processed raw aeromagnetic data into a gridded form. It is in text-readable GRD format and is a matrix of data composed of multiple rows and columns. Therefore, a W transform can be performed on any data point in the gridded aeromagnetic data, such as a row or column, to obtain the W transform coefficients corresponding to that row or column.

[0023] The process of obtaining gridded aeromagnetic data includes the following steps: obtaining raw aeromagnetic data, performing pole conversion calculations on the raw aeromagnetic data to obtain aeromagnetic pole conversion data; and regularly gridding the aeromagnetic pole conversion data to obtain gridded aeromagnetic data in text-readable grd format.

[0024] Specifically, raw aeromagnetic data can be obtained by acquiring publicly available data through open data platforms, or by using drones for high-precision exploration of small areas. Typically, data acquisition is performed on specific geographical areas requiring aeromagnetic anomaly identification. Polarization calculation and aeromagnetic polarization data gridding are both standard techniques in this field.

[0025] S2. Set the basic parameters of W transformation, and perform W transformation on one of the data of the gridded aeromagnetic data to obtain the W transformation air wave coefficient; The basic parameters of the W-transform include window function parameters, wavenumber range parameters, and calculation control parameters. The window function parameters include window type and window width parameters, used to control the width of the Gaussian window to achieve wavenumber-adaptive spatial wave analysis. In this invention, the window type is a Gaussian window. The wavenumber range parameters include the signal sampling wavenumber, and parameters for suppressing noise or enhancing specific wavebands through wavenumber weighting. The calculation control parameters include the sliding step size Δx.

[0026] The Gaussian window is a Gaussian window function with wavenumber weights. This window function uses a non-stationary wavenumber weight definition. The new Gaussian window function aligns the centroid of the spatial spectrum energy distribution with the dominant wavenumber of the waveform. As the absolute value of the difference between the instantaneous spatial wavenumber and the dominant wavenumber increases, the wavenumber weights gradually decrease. Therefore, the spatial spectrum energy coefficients are concentrated near the dominant wavenumber and exhibit symmetry at the dominant wavenumber of the signal. The spatial spectrum obtained through the W transform has high spatial resolution in both the high and low wave components. The expression for the Gaussian window function used is:

[0027] in, The dominant wavenumber of the waveform. , representing the deviation between the actual instantaneous wavenumber and the current analyzed wavenumber; K is a scalar factor that controls the width of the window; x For spatial location, k For space wavenumber, It is a spatial location variable.

[0028] Furthermore, the value of the scalar factor K should be between 0.1 and 1. The value of the scalar factor K is usually determined based on the characteristics of the data being processed. Based on simulations and testing with actual data, it is recommended that the value of the scalar factor K be between 0.1 and 1 in aeromagnetic data processing. A scalar factor K greater than 1 will lead to insufficient spatial resolution, making it impossible to accurately delineate the anomaly range. A scalar factor K less than 0.1 will lead to insufficient wavenumber resolution, making it impossible to determine the relative burial depth of the anomaly. Parameter selection is limited by the Heisenberg uncertainty principle: the smaller the scalar factor K, the better the spatial resolution and the worse the wavenumber resolution, and vice versa. Therefore, in application, the value of the scalar factor K should be maintained to ensure a relative balance between wavenumber resolution and spatial resolution in the analysis results, but the preference for wavenumber or spatial resolution can be adjusted according to actual needs.

[0029] As a preferred implementation, in establishing the Gaussian window function, the standard deviation σ of the Gaussian window function should be ≤ 4Δs, where Δs is the sampling interval.

[0030] Specifically, regarding the determination of the standard deviation of the Gaussian window function, the standard deviation σ should be ≤ 4Δs, where Δs is the sampling interval. A standard deviation σ greater than 4Δs will cause the air wave coefficients in the low-frequency part to diverge, which is not conducive to the identification of the aeromagnetic anomaly range.

[0031] The formula for the W transformation is:

[0032] in, ≠0 represents the dominant wavenumber of the waveform. K is a scalar factor, representing the initial width of the Gaussian window function.

[0033] The dominant wavenumber is approximated by a weighted average of instantaneous wavenumbers in the W variation. Its expression is:

[0034] in, L (s) is a low-pass filter defined by a Gaussian function; r This is a scaling factor, the size of which is inversely proportional to the standard deviation of the low-pass filter, and its specific value can be obtained through trial and error. It is the instantaneous amplitude envelope, used to estimate the dominant wavenumber. The weighting function; The instantaneous wavenumber of the input signal.

[0035] S3. Calculate the instantaneous space wave number based on the W-transform space wave coefficient, and define the synchronization extraction operator using the Dirac function. Calculate the distribution of the synchronization extraction operator value based on the instantaneous space wave number. By differentiating the W-transform spatial coefficients with respect to space, the signal can be obtained. instantaneous wavenumber Its expression is:

[0036] Specifically, a synchronization extraction operator is calculated based on the instantaneous space number, and then the space wave coefficients at the instantaneous space number location are extracted using the synchronization extraction operator. This process removes divergent space wave coefficients, retaining only the space wave coefficients at the instantaneous space number location, thus completing the directional extraction of W-transform space wave coefficients. The expression for this is:

[0037] S4. Perform matrix point-to-point multiplication on the W-transform spatial wave coefficients and the synchronous extraction operator to obtain the synchronous extraction W-transform spatial wave coefficient matrix. Specifically, the synchronous extraction W transform can be abbreviated as SEWT, and its expression is:

[0038] The algorithm for synchronous extraction is reversible, and the original signal can be reconstructed through inverse transformation.

[0039] Because the synchronous extraction process involves a matrix point-to-point multiplication of the W-transform spatial wave coefficients with the synchronous extraction operator, the algorithm is reversible. The original signal can be reconstructed through the inverse transform, and the original signal can be restored from the SEWT spatial wave representation, ensuring no information loss during analysis. This ensures a strict correspondence between the spatial wave analysis results and the physical characteristics of the original signal, and verifies the accuracy and reliability of SEWT through reconstruction errors. The reconstructed expression is:

[0040] in, A This represents the amplitude coefficient of the signal.

[0041] S5. Process each track of gridded aeromagnetic data within the work area to obtain the synchronous extraction W-transform airwave coefficients for each track of aeromagnetic data, and store the obtained synchronous extraction W-transform coefficients in the same matrix; specifically, after processing each track of gridded aeromagnetic data within the work area, store the synchronous extraction W-transform results of all data in the entire work area in the same three-dimensional matrix.

[0042] S6. Slice the matrix along the wavenumber direction to obtain aeromagnetic data information under different wavenumber distributions. Based on the changes in the position and amplitude of energy clusters on different wavenumber slices, identify the location and burial depth of the concealed magnetic rock mass.

[0043] By slicing the matrix along the wavenumber direction, aeromagnetic data under different wavenumber distributions are obtained. Based on the changes in the position and amplitude of energy clusters on different wavenumber slices, the location and burial depth of concealed magnetic rock masses are inferred. According to the order in which energy clusters appear on the wavenumber slices, if energy clusters are only present on low-wavenumber slices, the magnetic body at that location is relatively deep; if energy clusters are only present on high-wavenumber slices, the magnetic body at that location is relatively shallow; if energy clusters are present on both low-wavenumber and high-wavenumber slices, the magnetic body at that location is shallow and has a certain burial depth.

[0044] The following case study further illustrates the specific process of the method of the present invention.

[0045] Figure 2 For the aeromagnetic data after pole conversion, the corresponding data was gridded and then subjected to synchronous W-transform extraction. A Gaussian window function with wavenumber weights was used, with specific parameters set as follows: scalar factor K = 0.4, and the standard deviation of the Gaussian low-pass filter was 5Δt. The final results are as follows. Figures 3a-3f As shown. Where 3a- Figure 3f These correspond to wavenumbers of 0.10m respectively. -1 Wave number 0.50m -1 Wave number 0.80m -1 Wave number 1.00m -1 Wave number 1.40m -1 Wave number 2.0m -1 The output results of the slice below show that, as can be seen from the figure, at a low wavenumber of 0.10m... -1 A large area of ​​energy distribution exists on the right side of the spatiotemporal spectrum, which corresponds to the distribution of aeromagnetic positive anomalies. When the wavenumber rises to 0.80m... -1 At that time, most of the energy distribution on the right side disappeared, leaving only the energy within the black circle. The two energy distributions in the red circle and the energy distribution in the lower left corner were well preserved. Based on the rate of energy disappearance, it can be determined that the magnetic body causing the magnetic anomaly in the right area is relatively buried deep. By further increasing the wavenumber, the relative burial depth of other anomalies can be determined based on the rate of energy disappearance.

[0046] A spatial location identification system for concealed magnetic rock masses based on synchronous extraction of W-transform, comprising: The acquisition module acquires gridded aeromagnetic data; The parameter setting module allows you to set the basic parameters for the W transform. The W-transform module performs a W-transform on one of the data streams of the gridded aeromagnetic data to obtain the W-transformed air wave coefficients. The instantaneous wave number calculation module calculates the instantaneous wave number based on the W-transform wave coefficients. The synchronous extraction operator calculation module defines the synchronous extraction operator using the Dirac function and calculates the distribution of synchronous extraction operator values ​​based on the instantaneous wavenumber. The matrix operation module performs a matrix point-to-point multiplication operation on the W-transform spatial wave coefficients and the synchronous extraction operator to obtain the synchronous extraction W-transform spatial wave coefficient matrix. The aeromagnetic data processing module processes each data point in the gridded aeromagnetic data within the work area to obtain the synchronously extracted W-transform air wave coefficients for each aeromagnetic data point, and stores the obtained synchronously extracted W-transform coefficients in the same matrix. The aeromagnetic anomaly identification module slices the matrix along the wavenumber direction to obtain aeromagnetic data information under different wavenumber distributions. Based on the changes in the position and amplitude of energy clusters on different wavenumber slices, it identifies the location and burial depth of concealed magnetic rock masses.

[0047] Specifically, the acquisition module includes a polarization unit and a gridding unit. The polarization unit is used to acquire the raw aeromagnetic data and perform polarization calculations on the raw aeromagnetic data to obtain aeromagnetic polarization data. The gridding unit is used to regularly grid the aeromagnetic polarization data to obtain gridded aeromagnetic data in text-readable grd format.

[0048] The parameter setting module includes a function unit and a parameter unit. The function unit is used to establish a Gaussian window function with wavenumber weights. The parameter unit is used to set the scalar factor K in the W transform to a value between 0.1 and 1, and in establishing the Gaussian window function, to set the standard deviation σ of the Gaussian window function to σ≤4Δs.

[0049] An electronic device includes a processor and a memory, the memory storing a computer program, the processor executing the computer program stored in the memory to implement the steps of a method for identifying the spatial location of concealed magnetic rock masses based on synchronous extraction W-transformation.

[0050] A storage medium storing a computer program, which, when executed by a processor, implements the steps of a method for identifying the spatial location of concealed magnetic rock masses based on synchronous extraction of W-transform.

[0051] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the apparatus disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the description is relatively simple; relevant parts can be referred to the method section.

[0052] The above description of the disclosed embodiments enables those skilled in the art to make or use the invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the invention. Therefore, the invention is not 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 spatial location identification of concealed magnetic rock masses based on synchronous extraction of W-transform, characterized in that, include: Acquire gridded aeromagnetic data; Set the basic parameters of the W transform, perform the W transform on one of the data in the gridded aeromagnetic data, and obtain the W transform air wave coefficient; The instantaneous space wave number is calculated based on the W-transform space wave coefficient, and the synchronization extraction operator is defined using the Dirac function. The distribution of the synchronization extraction operator value is calculated based on the instantaneous space wave number. Perform a matrix point-to-point multiplication operation between the W-transform spatial wave coefficients and the synchronous extraction operator to obtain the synchronous extraction W-transform spatial wave coefficient matrix. Each data point in the gridded aeromagnetic data within the work area is processed to obtain the synchronously extracted W-transform air wave coefficients for each aeromagnetic data point, and the obtained synchronously extracted W-transform coefficients are stored in the same matrix. By slicing the matrix along the wavenumber direction, aeromagnetic data information under different wavenumber distributions is obtained. Based on the changes in the position and amplitude of energy clusters on different wavenumber slices, the location and burial depth of concealed magnetic rock masses are identified.

2. The method for spatial location identification of concealed magnetic rock masses based on synchronous extraction of W-transform as described in claim 1, characterized in that, The process of obtaining gridded aeromagnetic data includes the following steps: obtaining raw aeromagnetic data, performing pole conversion calculations on the raw aeromagnetic data to obtain aeromagnetic pole conversion data; and regularly gridding the aeromagnetic pole conversion data to obtain gridded aeromagnetic data in text-readable grd format.

3. The method for spatial location identification of concealed magnetic rock masses based on synchronous extraction of W-transform as described in claim 1, characterized in that, The basic parameters of the W transform include window function parameters, wavenumber range parameters, and calculation control parameters.

4. The method for spatial location identification of concealed magnetic rock masses based on synchronous extraction of W-transform as described in claim 1, characterized in that, The expression for the W-transform spatial wave coefficient is: in, x For spatial location, k For space wavenumber, This is single-track aeromagnetic data. ≠0 represents the dominant wavenumber of the waveform. , is a scalar factor, and is the initial width of the Gaussian window function.

5. The method for spatial location identification of concealed magnetic rock masses based on synchronous extraction of W-transform as described in claim 4, characterized in that, The expression for the synchronous extraction operator is: In the formula, For the Dirac function, It represents the instantaneous wave number.

6. A spatial location identification system for concealed magnetic rock masses based on synchronous extraction of W-transform, characterized in that, The method for spatial location identification of concealed magnetic rock masses based on synchronous extraction of W-transform as described in any one of claims 1-5 includes: The acquisition module acquires gridded aeromagnetic data; The parameter setting module allows you to set the basic parameters for the W transform. The W-transform module performs a W-transform on one of the data streams of the gridded aeromagnetic data to obtain the W-transformed air wave coefficients. The instantaneous space wavenumber calculation module calculates the instantaneous space wavenumber based on the W-transform space wave coefficient. The synchronous extraction operator calculation module defines the synchronous extraction operator using the Dirac function and calculates the distribution of synchronous extraction operator values ​​based on the instantaneous wavenumber. The matrix operation module performs a matrix point-to-point multiplication operation on the W-transform spatial wave coefficients and the synchronous extraction operator to obtain the synchronous extraction W-transform spatial wave coefficient matrix. The aeromagnetic data processing module processes each data point in the gridded aeromagnetic data within the work area to obtain the synchronously extracted W-transform airwave coefficients for each aeromagnetic data point, and stores the obtained synchronously extracted W-transform coefficients in the same matrix. The aeromagnetic anomaly identification module slices the matrix along the wavenumber direction to obtain aeromagnetic data information under different wavenumber distributions. Based on the changes in the position and amplitude of energy clusters on different wavenumber slices, it identifies the location and burial depth of concealed magnetic rock bodies.

7. An electronic device comprising a processor and a memory, the memory storing a computer program, characterized in that, The processor executes the computer program stored in the memory to implement the steps of the method for identifying the spatial location of concealed magnetic rock masses based on synchronous extraction W-transformation as described in any one of claims 1 to 5.

8. A storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the steps of the method for identifying the spatial location of concealed magnetic rock masses based on synchronous extraction W-transformation as described in any one of claims 1 to 5.

Images