Underground rock mass magnetic target positioning method, system and equipment based on aeromagnetic three-component data, medium and product

Abstract

The invention discloses an underground rock mass magnetic target positioning method, system and device based on aeromagnetic three-component data, a medium and a product, and relates to the field of geophysical exploration, and the method comprises the steps: obtaining the aeromagnetic three-component data of a plurality of measurement points in a measurement region; determining the magnetization direction of a magnetic target in the measurement area according to the aeromagnetic three-component data of the plurality of measurement points in the measurement area; the magnetization direction comprises a total magnetization inclination angle and a total magnetization deflection angle; and according to the aeromagnetic three-component data, determining the position and the magnetic susceptibility of the magnetic target in the measurement area by using a step-by-step weighted inversion method. According to the method, rich geological information in the aeromagnetic three-component data can be effectively utilized, the inversion stability and resolution are improved in combination with the step-by-step weighted inversion method, and the position and direction information of the underground rock mass can be rapidly and accurately obtained.

Description

Technical Field

[0001] This application relates to the field of geophysical exploration, and in particular to a method, system, equipment, medium and product for locating magnetic targets in underground rock masses based on aeromagnetic three-component data. Background Technology

[0002] Traditional aeromagnetic measurements primarily focus on measuring the total field strength, followed by mathematical processing to obtain the northward, eastward, and vertical components of the magnetic field to infer the spatial distribution of underground magnetic anomalies. However, this method suffers from accumulated errors during data conversion, particularly prominent in low-latitude regions where magnetic pole reorientation is a significant issue. Summary of the Invention

[0003] The purpose of this application is to provide a method, system, equipment, medium, and product for locating magnetic targets in underground rock masses based on aeromagnetic three-component data, so as to improve the detection resolution and positioning accuracy of underground rock masses.

[0004] To achieve the above objectives, this application provides the following solution:

[0005] Firstly, this application provides a method for locating magnetic targets in underground rock masses based on aeromagnetic three-component data, including:

[0006] Acquire aeromagnetic three-component data from multiple measurement points within the measurement area;

[0007] Based on the aeromagnetic three-component data from multiple measurement points in the measurement area, the magnetization direction of the magnetic target in the measurement area is determined; the magnetization direction includes the total magnetization tilt angle and the total magnetization deflection angle.

[0008] Based on the aeromagnetic three-component data, the position and magnetization of the magnetic target within the measurement area are determined using a stepwise weighted inversion method.

[0009] Optionally, the magnetization direction of the magnetic target within the measurement area is determined based on the aeromagnetic three-component data of the measurement area, specifically including:

[0010] Based on the aeromagnetic three-component data of multiple measurement points in the measurement area, calculate the magnetization direction of each measurement point in the measurement area;

[0011] When there is only one magnetic target in the measurement area, the magnetization direction of each measurement point is calculated multiple times.

[0012] The magnetization direction of the measurement point with the smallest change in magnetization direction is taken as the magnetization direction of the magnetic target;

[0013] When there are multiple magnetic targets within the measurement area, the measurement area is divided into multiple independent regions using a boundary recognition method; each independent region contains one magnetic target.

[0014] The magnetization direction of each measurement point within each independent region is calculated multiple times;

[0015] The magnetization direction of the measurement point with the smallest change in magnetization direction is taken as the magnetization direction of the magnetic target within the independent region.

[0016] Optionally, the magnetization direction of the measurement point is calculated based on the aeromagnetic three-component data of the measurement point within the measurement area, specifically including:

[0017] The magnetic moment component of the measurement point is obtained by integrating the aeromagnetic three-component data of the measurement point using the Helbig method.

[0018] The magnetization direction of the measurement point is calculated based on the magnetic moment component of the measurement point.

[0019] Optionally, the aeromagnetic three-component data of the measurement point are integrated using the Helbig method to obtain the magnetic moment component of the measurement point, specifically including:

[0020] Using formula Determine the magnetic moment components at the measurement point; where m x m represents the x-direction magnetic moment component of the measurement point. y m represents the y-direction magnetic moment component of the measurement point. z B represents the z-direction magnetic moment component of the measurement point; x and y are the coordinates of the measurement point in the right-hand coordinate system; x B is the x-direction component of the aeromagnetic three-component data at the measurement point; z The z-direction component is the aeromagnetic three-component data of the measurement point.

[0021] Optionally, the magnetization direction of the measurement point is calculated based on the magnetic moment component of the measurement point, specifically including:

[0022] Using formula Calculate the total magnetization tilt angle at the measurement point; where, I m The total magnetization tilt angle; m x m represents the x-direction magnetic moment component of the measurement point. y m represents the y-direction magnetic moment component of the measurement point. z The z-direction magnetic moment component of the measurement point;

[0023] Using formula Calculate the total magnetization deflection at the measurement point; where, D m This represents the total magnetization deflection.

[0024] Optionally, the magnetization direction of the magnetic target within the measurement area is determined based on the aeromagnetic three-component data from multiple measurement points in the measurement area, further including:

[0025] The aeromagnetic three-component data are subjected to geomagnetic field correction, cutting line leveling, and gridding to obtain the processed aeromagnetic three-component data.

[0026] Secondly, this application provides a magnetic target positioning system for underground rock masses based on aeromagnetic three-component data. This system is used to implement the aforementioned magnetic target positioning method for underground rock masses based on aeromagnetic three-component data. The system includes:

[0027] The data acquisition module is used to acquire aeromagnetic three-component data from multiple measurement points within the measurement area;

[0028] The magnetization direction determination module is used to determine the magnetization direction of a magnetic target within the measurement area based on the aeromagnetic three-component data of multiple measurement points within the measurement area; the magnetization direction includes the total magnetization tilt angle and the total magnetization deflection angle;

[0029] The location determination module is used to determine the location and magnetic susceptibility of magnetic targets within the measurement area based on the aeromagnetic three-component data and using a stepwise weighted inversion method.

[0030] Thirdly, this application provides a computer device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the underground rock mass magnetic target positioning method based on aeromagnetic three-component data as described above.

[0031] Fourthly, this application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for locating underground rock magnetic targets based on aeromagnetic three-component data as described above.

[0032] Fifthly, this application provides a computer program product, including a computer program that, when executed by a processor, implements the method for locating underground rock magnetic targets based on aeromagnetic three-component data as described above.

[0033] According to the specific embodiments provided in this application, this application has the following technical effects:

[0034] This application provides a method, system, equipment, medium, and product for locating magnetic targets in underground rock masses based on aeromagnetic three-component data. The method involves acquiring aeromagnetic three-component data from multiple measurement points within a measurement area; determining the magnetization direction of magnetic targets within the measurement area based on the aeromagnetic three-component data; the magnetization direction including the total magnetization dip angle and the total magnetization deflection angle; and determining the location and magnetic susceptibility of magnetic targets within the measurement area using a stepwise weighted inversion method based on the aeromagnetic three-component data. This application effectively utilizes the rich geological information in the aeromagnetic three-component data and improves the stability and resolution of the inversion by combining it with a stepwise weighted inversion method, enabling rapid and accurate acquisition of the location and orientation information of underground rock masses. Attached Figure Description

[0035] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0036] Figure 1 A flowchart illustrating a method for locating underground rock magnetic targets based on aeromagnetic three-component data, provided in an embodiment of this application;

[0037] Figure 2 This is a schematic diagram of the stepwise weighted inversion.

[0038] Figure 3 This is a schematic diagram of the structure of a computer device provided in an embodiment of this application. Detailed Implementation

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

[0040] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, the application will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0041] With the advancement of aeromagnetic three-component observation technology, it has become possible to directly use aeromagnetic three-component data for inversion calculations, which helps to improve detection resolution and positioning accuracy.

[0042] This application provides a method for locating magnetic targets in underground rock masses based on aeromagnetic three-component data. This method can effectively utilize the rich geological information in the aeromagnetic three-component data to improve the accuracy and efficiency of the inversion results.

[0043] In one exemplary embodiment, such as Figure 1 As shown, a method for locating magnetic targets in underground rock masses based on aeromagnetic three-component data is provided, including the following steps:

[0044] S1: Acquire aeromagnetic three-component data from multiple measurement points within the measurement area. The aeromagnetic three-component data includes the northward component (X), the eastward component (Y), and the vertical component (Z). Noise removal and data normalization are performed on the acquired data.

[0045] S2: Based on the aeromagnetic three-component data of multiple measurement points in the measurement area, determine the magnetization direction of the magnetic target in the measurement area; the magnetization direction includes the total magnetization tilt angle and the total magnetization deflection angle.

[0046] In practical applications, measured aeromagnetic three-component data are used to calculate the magnetization direction of the rock mass (including the total magnetization dip angle and the total magnetization deflection angle).

[0047] As an optional implementation, S2 further includes:

[0048] The aeromagnetic three-component data are subjected to geomagnetic field correction, cutting line leveling, and gridding to obtain processed aeromagnetic three-component data, eliminating noise and interference.

[0049] S3: Based on the aeromagnetic three-component data, the position and magnetic susceptibility of the magnetic target within the measurement area are determined using a stepwise weighted inversion method.

[0050] In practical applications, the specific process of the stepwise weighted inversion method is as follows:

[0051] Rapid inversion initial screening: Using automatic inversion methods (such as the Euler deconvolution method) and combining tectonic indices to estimate the source depth, an initial model is provided for fine inversion, which can obtain the location (x, y, z) of the underground rock mass.

[0052] Physical property inversion: The subsurface is divided into prismatic units. Using aeromagnetic three-component data as input, 3D magnetic susceptibility inversion is performed using a stepwise weighting method to obtain the location and physical property information (x, y, z, m) of the subsurface rock mass, where m is the magnetic susceptibility.

[0053] Results fusion: By combining geological data and rock mass orientation and location information calculated from aeromagnetic three-component data, a comprehensive interpretation of the flight data acquisition area is performed.

[0054] As an optional implementation, S2 specifically includes:

[0055] S21: Calculate the magnetization direction of each measurement point in the measurement area based on the aeromagnetic three-component data of multiple measurement points in the measurement area.

[0056] S22: When there is only one magnetic target in the measurement area, calculate the magnetization direction of each measurement point multiple times.

[0057] S23: The magnetization direction of the measurement point with the smallest change in magnetization direction is taken as the magnetization direction of the magnetic target.

[0058] S24: When there are multiple magnetic targets in the measurement area, the measurement area is divided into multiple independent areas using a boundary recognition method; each independent area contains one magnetic target.

[0059] S25: Calculate the magnetization direction of each measurement point in each of the independent regions multiple times.

[0060] S26: The magnetization direction of the measurement point with the smallest change in magnetization direction is taken as the magnetization direction of the magnetic target in the independent region.

[0061] As an optional implementation, S21 specifically includes:

[0062] S211: Integrate the aeromagnetic three-component data of the measurement point using the Helbig method to obtain the magnetic moment component of the measurement point.

[0063] As an optional implementation, S211 specifically includes:

[0064] Using formula Determine the magnetic moment components at the measurement point; where m x m represents the x-direction magnetic moment component of the measurement point. y m represents the y-direction magnetic moment component of the measurement point. z B represents the z-direction magnetic moment component of the measurement point; x and y are the coordinates of the measurement point in the right-hand coordinate system; x B is the x-direction component of the aeromagnetic three-component data at the measurement point; z The z-direction component is the aeromagnetic three-component data of the measurement point.

[0065] S212: Calculate the magnetization direction of the measurement point based on the magnetic moment component of the measurement point.

[0066] As an optional configuration, S212 specifically includes:

[0067] Using formula Calculate the total magnetization tilt angle at the measurement point; where, I m The total magnetization tilt angle; m xm represents the x-direction magnetic moment component of the measurement point. y m represents the y-direction magnetic moment component of the measurement point. z Let be the z-direction magnetic moment component of the measurement point.

[0068] Using formula Calculate the total magnetization deflection at the measurement point; where, D m This represents the total magnetization deflection.

[0069] In this embodiment, the aeromagnetic three-component data are first integrated using the Helbig method to obtain the magnetic moment component (m). x m y m z Then, the magnetic inclination angle I is obtained according to formulas (5) and (6). m And magnetic declination D m .

[0070] The integral relationship between Helbig's magnetic moment and the magnetic field components is as follows:

[0071]

[0072] After calculating the magnetic moments in the three directions, the total magnetization direction can be further calculated.

[0073]

[0074] in, This is the total magnetic moment modulus.

[0075] As can be seen from formulas (1)-(3), infinite integration is required along the z-plane. In order to facilitate the application of aeromagnetic gridded data, a sliding window is used to calculate the two-dimensional trapezoidal quadrature formula. Under the same amount of computation, the Simpson quadrature formula is more accurate than the trapezoidal quadrature formula. Therefore, in this embodiment, the Simpson quadrature formula is used to calculate the total magnetization direction.

[0076] Since the magnetization direction calculated by the Helbig method at the horizontal position of the magnetic source does not change with the window size, the magnetization direction at each point within the measurement area can be calculated using the Helbig method by selecting multiple sliding windows of different sizes. When the measurement area contains only one magnetic target, the horizontal position of the measurement point with the smallest change in magnetization direction obtained from multiple calculations is the horizontal position of the magnetic target, and the magnetic inclination and magnetic declination of this point are the total magnetization direction of the magnetic target. When the measurement area contains multiple targets, the area is first divided into regions using boundary identification methods, and then the magnetization direction is estimated for each independent region containing only one target.

[0077] A weighted matrix is ​​constructed based on one or more components, and the weights are used to guide the stepwise weighted inversion to locate the underground rock mass.

[0078] Given that combining multi-component data can improve inversion results, a stepwise weighted inversion method is applied to the inversion of aeromagnetic three-component data. First, a deep weighted inversion is performed, then a weight matrix is ​​constructed based on multiple components, and finally, the weighted components are used for the next inversion calculation. Guided by weights, areas with larger weights participate in the inversion. This effectively utilizes information from different directions in each aeromagnetic component data, quickly obtaining inversion results.

[0079] Figure 2 The flowchart of the stepwise weighted inversion method is given. Figure 2 The process involves two different weighting matrices: one (W2) is constructed from a single component and a depth weighting function W1, and the other (W3, ..., W...). i+1 It is constructed from multiple components. The overall process of constructing the weighted matrix consists of three parts.

[0080] In the first part, a weighting matrix based on a single component was obtained. The data is B1, and its depth weighting function is W1. To avoid introducing other information, a depth weighting function based on the sensitivity matrix is ​​used. The expression for W1 is defined as:

[0081]

[0082] Among them, the element a in the sensitivity matrix A ij This represents the influence of unit magnetization at data position i (i = 1, ..., M) for cell j (j = 1, ..., N); β is a constant, with larger values ​​indicating stronger weighting. Typically, β is close to 1. Only the fitting function is used in constructing the weighted matrix. The objective function is:

[0083]

[0084] The conjugate gradient algorithm is used to solve this problem. The obtained inversion result is m.

[0085] The second part demonstrates the method for constructing a weight matrix based on multiple components. After obtaining W2, another component B2 is selected to obtain the inversion result m. ​​The absolute values ​​of the components of m are taken, and m is normalized. Then, sufficiently small values ​​(10) are assigned to elements equal to 0. -6 (or smaller), resulting in W3. By analogy, we obtain W3 based on... i and B i The weighted matrix W i+1 .

[0086] In the third part, select the last component B in the component combination. α Perform inversion. Based on B... α and W i+1The objective function, i.e., formula (9), is established and solved using the nonlinear conjugate gradient algorithm. The component ultimately used in the inversion depends on the actual geological conditions. Since the vertical component alone can yield results that basically conform to geological interpretation, the vertical component B is selected first. z As B α .

[0087]

[0088] Compared with related technologies, this application has the following advantages:

[0089] Reduced errors: Directly utilizing aeromagnetic three-component data reduces errors during the data conversion process.

[0090] Improved stability and resolution: The stepwise weighted inversion method improves the stability and resolution of the inversion.

[0091] Rapid and accurate information acquisition: It can quickly and accurately obtain the location and orientation information of underground rock masses, providing new technical means for mineral resource exploration.

[0092] Based on the same inventive concept, this application also provides a system for locating underground rock magnetic targets based on aeromagnetic three-component data, which is used to implement the above-mentioned method for locating underground rock magnetic targets based on aeromagnetic three-component data. The solution provided by this system is similar to the solution described in the above method. Therefore, the specific limitations in the embodiments of the system for locating underground rock magnetic targets based on aeromagnetic three-component data provided below can be found in the limitations of the method for locating underground rock magnetic targets based on aeromagnetic three-component data described above, and will not be repeated here.

[0093] In one exemplary embodiment, a magnetic target localization system for underground rock masses based on aeromagnetic three-component data is provided, comprising:

[0094] The data acquisition module is used to acquire aeromagnetic three-component data from multiple measurement points within the measurement area.

[0095] The magnetization direction determination module is used to determine the magnetization direction of a magnetic target within the measurement area based on the aeromagnetic three-component data of multiple measurement points in the measurement area; the magnetization direction includes the total magnetization tilt angle and the total magnetization deflection angle.

[0096] The location determination module is used to determine the location and magnetic susceptibility of magnetic targets within the measurement area based on the aeromagnetic three-component data and using a stepwise weighted inversion method.

[0097] In one exemplary embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the above-described method for locating underground rock magnetic targets based on aeromagnetic three-component data.

[0098] In one exemplary embodiment, a computer-readable storage medium is provided storing a computer program that, when executed by a processor, implements the above-described method for locating underground rock magnetic targets based on aeromagnetic three-component data.

[0099] In one exemplary embodiment, a computer program product is provided, including a computer program that, when executed by a processor, implements the above-described method for locating underground rock magnetic targets based on aeromagnetic three-component data.

[0100] In one exemplary embodiment, a computer device is provided, which may be a server or a terminal, and its internal structure diagram may be as follows. Figure 3 As shown, this computer device includes a processor, memory, input / output (I / O) interfaces, and a communication interface. The processor, memory, and I / O interfaces are connected via a system bus, and the communication interface is also connected to the system bus via the I / O interfaces. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system, computer programs, and databases. The internal memory provides the environment for the operating system and computer programs stored in the non-volatile storage media to run. The I / O interfaces are used for exchanging information between the processor and external devices. The communication interface is used for communicating with external terminals via a network connection. When the computer program is executed by the processor, it implements a method for locating underground rock magnetic targets based on aeromagnetic three-component data.

[0101] Those skilled in the art will understand that Figure 3 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0102] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of the relevant data must comply with relevant regulations.

[0103] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0104] The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0105] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0106] This document uses specific examples to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the methods and core ideas of this application. Furthermore, those skilled in the art will recognize that, based on the ideas of this application, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A method for locating magnetic targets in underground rock masses based on aeromagnetic three-component data, characterized in that, include: Acquire aeromagnetic three-component data from multiple measurement points within the measurement area; Based on the aeromagnetic three-component data from multiple measurement points in the measurement area, the magnetization direction of the magnetic target in the measurement area is determined; the magnetization direction includes the total magnetization tilt angle and the total magnetization deflection angle. Based on the aeromagnetic three-component data, the position and magnetization of the magnetic target within the measurement area are determined using a stepwise weighted inversion method.

2. The method for locating underground rock magnetic targets based on aeromagnetic three-component data according to claim 1, characterized in that, Based on the aeromagnetic three-component data of the measurement area, the magnetization direction of the magnetic target within the measurement area is determined, specifically including: Based on the aeromagnetic three-component data of multiple measurement points in the measurement area, calculate the magnetization direction of each measurement point in the measurement area; When there is only one magnetic target in the measurement area, the magnetization direction of each measurement point is calculated multiple times. The magnetization direction of the measurement point with the smallest change in magnetization direction is taken as the magnetization direction of the magnetic target; When there are multiple magnetic targets within the measurement area, the measurement area is divided into multiple independent regions using a boundary recognition method; each independent region contains one magnetic target. The magnetization direction of each measurement point within each independent region is calculated multiple times; The magnetization direction of the measurement point with the smallest change in magnetization direction is taken as the magnetization direction of the magnetic target within the independent region.

3. The method for locating underground rock magnetic targets based on aeromagnetic three-component data according to claim 2, characterized in that, Based on the aeromagnetic three-component data of the measurement points within the measurement area, the magnetization direction of the measurement points is calculated, specifically including: The magnetic moment component of the measurement point is obtained by integrating the aeromagnetic three-component data of the measurement point using the Helbig method. The magnetization direction of the measurement point is calculated based on the magnetic moment component of the measurement point.

4. The method for locating underground rock mass magnetic targets based on aeromagnetic three-component data according to claim 3, characterized in that, The magnetic moment components of the measurement point are obtained by integrating the three aeromagnetic components using the Helbig method, specifically including: Using formula Determine the magnetic moment components at the measurement point; where m x m represents the x-direction magnetic moment component of the measurement point. y The y-direction magnetic moment component of the measurement point; m z B represents the z-direction magnetic moment component of the measurement point; x and y are the coordinates of the measurement point in the right-hand coordinate system; x B is the x-direction component of the aeromagnetic three-component data at the measurement point; z The z-direction component is the aeromagnetic three-component data of the measurement point.

5. The method for locating underground rock mass magnetic targets based on aeromagnetic three-component data according to claim 3, characterized in that, Calculating the magnetization direction of the measurement point based on its magnetic moment components includes: Using formula Calculate the total magnetization tilt angle at the measurement point; where, I m The total magnetization tilt angle; m x m represents the x-direction magnetic moment component of the measurement point. y The y-direction magnetic moment component of the measurement point; m z The z-direction magnetic moment component of the measurement point; Using formula Calculate the total magnetization deflection at the measurement point; where, D m This represents the total magnetization deflection.

6. The method for locating underground rock magnetic targets based on aeromagnetic three-component data according to claim 1, characterized in that, Based on the aeromagnetic three-component data from multiple measurement points in the measurement area, the magnetization direction of the magnetic target within the measurement area is determined, which also includes: The aeromagnetic three-component data are subjected to geomagnetic field correction, cutting line leveling, and gridding to obtain the processed aeromagnetic three-component data.

7. A magnetic target positioning system for underground rock masses based on aeromagnetic three-component data, characterized in that, The underground rock mass magnetic target positioning system based on aeromagnetic three-component data is used to implement the underground rock mass magnetic target positioning method based on aeromagnetic three-component data as described in any one of claims 1-6. The underground rock mass magnetic target positioning system based on aeromagnetic three-component data includes: The data acquisition module is used to acquire aeromagnetic three-component data from multiple measurement points within the measurement area; The magnetization direction determination module is used to determine the magnetization direction of a magnetic target within the measurement area based on the aeromagnetic three-component data of multiple measurement points within the measurement area; the magnetization direction includes the total magnetization tilt angle and the total magnetization deflection angle; The location determination module is used to determine the location and magnetic susceptibility of magnetic targets within the measurement area based on the aeromagnetic three-component data and using a stepwise weighted inversion method.

8. A computer device, comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that the processor executes the computer program to implement the method for locating underground rock magnetic targets based on aeromagnetic three-component data as described in any one of claims 1-6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the method for locating underground rock magnetic targets based on aeromagnetic three-component data as described in any one of claims 1-6.

10. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the method for locating underground rock magnetic targets based on aeromagnetic three-component data as described in any one of claims 1-6.

Images