An underground rock mass magnetic target positioning method, system, device, medium and product based on aeromagnetic three-component data

Abstract

The application discloses a kind of underground rock mass magnetic target positioning method, system, equipment, medium and product based on aeromagnetic three-component data, it is related to geophysical prospecting field, the method includes obtaining the aeromagnetic three-component data of multiple measuring points in measurement area;According to the aeromagnetic three-component data of multiple measuring points in the measurement area, the magnetization direction of the magnetic target in the measurement area is determined;The magnetization direction includes total magnetization dip angle and total magnetization declination;According to the aeromagnetic three-component data, the position and magnetic susceptibility of the magnetic target in the measurement area are determined using step-by-step weighting inversion method.The present application can effectively utilize the rich geological information in aeromagnetic three-component data, combined with step-by-step weighting inversion method to improve the stability and resolution of inversion, can quickly and accurately obtain the position and direction information of underground rock mass.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of geophysical exploration, and particularly relates to a method, system, device, medium and product for locating a magnetic target of a subsurface rock mass based on three-component aeromagnetic data. BACKGROUND

[0002] Traditional aeromagnetic survey mainly focuses on the measurement of total field intensity, and then the north, east and vertical components of the magnetic field are obtained through mathematical processing to infer the spatial distribution of underground magnetic anomaly bodies. However, this method has error accumulation in data conversion, especially in low-latitude areas where the geomagnetic pole problem is prominent. SUMMARY

[0003] The purpose of the present application is to provide a method, system, device, medium and product for locating a magnetic target of a subsurface rock mass based on three-component aeromagnetic data, so as to improve the detection resolution and positioning accuracy of the subsurface rock mass.

[0004] To achieve the above purpose, the present application provides the following solutions:

[0005] In a first aspect, the present application provides a method for locating a magnetic target of a subsurface rock mass based on three-component aeromagnetic data, comprising:

[0006] obtaining three-component aeromagnetic data of a plurality of measurement points in a measurement area;

[0007] determining a magnetization direction of a magnetic target in the measurement area according to the three-component aeromagnetic data of the plurality of measurement points in the measurement area; the magnetization direction includes a total magnetization inclination and a total magnetization declination;

[0008] determining the position and magnetic susceptibility of the magnetic target in the measurement area by using a step-by-step weighted inversion method according to the three-component aeromagnetic data.

[0009] Optionally, the determination of the magnetization direction of the magnetic target in the measurement area according to the three-component aeromagnetic data of the measurement area specifically comprises:

[0010] calculating the magnetization direction of each measurement point in the measurement area according to the three-component aeromagnetic data of the plurality of measurement points 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 in the measurement area, the measurement area is divided into multiple independent regions by using a boundary recognition method; one independent region contains one magnetic target;

[0014] calculating the magnetization direction of each measuring point in each of the independent regions;

[0015] taking the magnetization direction of the measuring point with the least magnetization direction change as the magnetization direction of the magnetic target in the independent region.

[0016] Optionally, the magnetization direction of the measuring point is calculated according to the aeromagnetic three-component data of the measuring point in the measuring region, and specifically includes:

[0017] integrating the aeromagnetic three-component data of the measuring point by using the Helbig method to obtain the magnetic moment component of the measuring point;

[0018] calculating the magnetization direction of the measuring point according to the magnetic moment component of the measuring point.

[0019] Optionally, the magnetic moment component of the measuring point is obtained by integrating the aeromagnetic three-component data of the measuring point by using the Helbig method, and specifically includes:

[0020] determining the magnetic moment component of the measuring point by using the formula ; wherein, m x is the x-direction magnetic moment component of the measuring point; m y is the y-direction magnetic moment component of the measuring point; m z is the z-direction magnetic moment component of the measuring point; x and y are coordinate values of the measuring point in a right-handed coordinate system; B x is the x-direction component in the aeromagnetic three-component data of the measuring point; B z is the z-direction component in the aeromagnetic three-component data of the measuring point.

[0021] Optionally, the magnetization direction of the measuring point is calculated according to the magnetic moment component of the measuring point, and specifically includes:

[0022] calculating the total magnetization inclination of the measuring point by using the formula ; wherein, I m is the total magnetization inclination; m x is the x-direction magnetic moment component of the measuring point; m y is the y-direction magnetic moment component of the measuring point; m z is the z-direction magnetic moment component of the measuring point.

[0023] calculating the total magnetization declination of the measuring point by using the formula ; wherein, D m is the total magnetization declination.

[0024] Optionally, the magnetization direction of the magnetic target in the measuring region is determined according to the aeromagnetic three-component data of the plurality of measuring points in the measuring region, and the method further includes:

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

[0026] In a second aspect, the application provides a system for locating a magnetic target in a subsurface rock mass based on aeromagnetic three-component data, which is used to implement the method for locating a magnetic target in a subsurface rock mass based on aeromagnetic three-component data described above. The system for locating a magnetic target in a subsurface rock mass based on aeromagnetic three-component data comprises:

[0027] a data acquisition module configured to acquire aeromagnetic three-component data of a plurality of measurement points in a measurement region;

[0028] a magnetization direction determination module configured to determine a magnetization direction of a magnetic target in the measurement region based on the aeromagnetic three-component data of the plurality of measurement points in the measurement region, wherein the magnetization direction comprises a total magnetization inclination and a total magnetization declination;

[0029] a position determination module configured to determine a position and a magnetic susceptibility of the magnetic target in the measurement region based on the aeromagnetic three-component data by using a step-by-step weighted inversion method.

[0030] In a third aspect, the application provides a computer device, which comprises 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 method for locating a magnetic target in a subsurface rock mass based on aeromagnetic three-component data according to any one of the above embodiments.

[0031] In a fourth aspect, the application provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program is executable by a processor to implement the method for locating a magnetic target in a subsurface rock mass based on aeromagnetic three-component data according to any one of the above embodiments.

[0032] In a fifth aspect, the application provides a computer program product comprising a computer program, wherein the computer program is executable by a processor to implement the method for locating a magnetic target in a subsurface rock mass based on aeromagnetic three-component data according to any one of the above embodiments.

[0033] According to the embodiments provided in the application, the following technical effects are achieved:

[0034] The application provides a method, system, device, medium and product for locating a magnetic target of an underground rock mass based on airborne magnetic three-component data, acquires airborne magnetic three-component data of multiple measuring points in a measuring area, determines a magnetization direction of the magnetic target in the measuring area according to the airborne magnetic three-component data of the multiple measuring points in the measuring area, wherein the magnetization direction comprises a total magnetization inclination and a total magnetization declination, and determines the position and magnetic susceptibility of the magnetic target in the measuring area by using a step-by-step weighted inversion method according to the airborne magnetic three-component data. The application can effectively utilize the rich geological information in the airborne magnetic three-component data, improve the stability and resolution of inversion by combining the step-by-step weighted inversion method, and quickly and accurately obtain the position and direction information of the underground rock mass. BRIEF DESCRIPTION OF DRAWINGS

[0035] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed in the embodiments will be briefly introduced as follows. Obviously, the drawings in the following description only constitute some embodiments of the present application, and other drawings can be obtained by those skilled in the art without any creative effort on the basis of these drawings.

[0036] Figure 1 A flowchart of a method for locating a magnetic target of an underground rock mass based on airborne magnetic three-component data is provided for an embodiment of the present application.

[0037] Figure 2 A step-by-step weighted inversion diagram is provided.

[0038] Figure 3 A structural diagram of a computer device is provided for an embodiment of the present application. DETAILED DESCRIPTION

[0039] The technical solutions in the embodiments of the present application will be described in detail below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments only constitute some embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without any creative effort fall within the scope of protection of the present application.

[0040] In order to make the above-mentioned purposes, features and advantages of the present application more obvious and easy to understand, the present application will be further described in detail below with reference to the drawings and specific embodiments.

[0041] With the development of airborne magnetic three-component observation technology, it is possible to directly use airborne magnetic three-component data for inversion calculation, which helps to improve the detection resolution and positioning accuracy.

[0042] The application provides a method for locating a magnetic target of an underground rock mass based on three-component aeromagnetic data, which can effectively utilize the rich geological information in the three-component aeromagnetic data and improve the accuracy and efficiency of inversion results.

[0043] In one exemplary embodiment, as shown in Figure 1 a method for locating a magnetic target of an underground rock mass based on three-component aeromagnetic data is provided, comprising the following steps:

[0044] S1: Obtain three-component aeromagnetic data of multiple measurement points in a measurement area. The three-component aeromagnetic data includes a north component X, an east component Y, and a vertical component Z. The collected data is subjected to noise removal and data normalization processing.

[0045] S2: Determine the magnetization direction of the magnetic target in the measurement area according to the three-component aeromagnetic data of the multiple measurement points in the measurement area; the magnetization direction includes a total magnetization inclination and a total magnetization declination.

[0046] In actual application, the measured three-component aeromagnetic data is used to calculate the magnetization direction of the rock mass (including the total magnetization inclination and the total magnetization declination).

[0047] As an optional implementation, S2 further comprises, before the step of determining the magnetization direction of the magnetic target in the measurement area according to the three-component aeromagnetic data of the multiple measurement points in the measurement area:

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

[0049] S3: Determine the position and magnetic susceptibility of the magnetic target in the measurement area according to the three-component aeromagnetic data by using a step-by-step weighted inversion method.

[0050] In actual application, the step-by-step weighted inversion method specifically comprises the following steps:

[0051] Fast inversion preliminary screening: an automatic inversion method (such as Euler deconvolution) is used in combination with a structure index to estimate the depth of the field source, thereby providing an initial model for fine inversion, and the position (x, y, z) of the underground rock mass can be obtained.

[0052] Physical property inversion: the underground is divided into prismatic units, the three-component aeromagnetic data is used as input, and a step-by-step weighted method is used for 3D magnetic susceptibility inversion, thereby obtaining the position and physical property information (x, y, z, m) of the underground rock mass, wherein m is the magnetic susceptibility.

[0053] Result fusion: the rock mass direction and position information calculated by the three-component aeromagnetic data in combination with geological data are used for comprehensive interpretation of the flight data acquisition area.

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

[0055] S21: calculating the magnetization direction of each measuring point in the measuring area according to the aeromagnetic three-component data of the measuring points in the measuring area.

[0056] S22: when there is only one magnetic target in the measuring area, the magnetization direction of each measuring point is calculated for multiple times.

[0057] S23: the magnetization direction of the measuring point with the least change of magnetization direction is taken as the magnetization direction of the magnetic target.

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

[0059] S25: the magnetization direction of each measuring point in each independent area is calculated for multiple times.

[0060] S26: the magnetization direction of the measuring point with the least change of magnetization direction is taken as the magnetization direction of the magnetic target in the independent area.

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

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

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

[0064] the magnetic moment component of the measuring point is determined by using the formula ; wherein, m x is the x-direction magnetic moment component of the measuring point; m y is the y-direction magnetic moment component of the measuring point; m z is the z-direction magnetic moment component of the measuring point; x and y are the coordinate values of the measuring point in the right-hand coordinate system; B x is the x-direction component in the aeromagnetic three-component data of the measuring point; B z is the z-direction component in the aeromagnetic three-component data of the measuring point.

[0065] S212: calculating the magnetization direction of the measuring point according to the magnetic moment component of the measuring point.

[0066] As an optional S212, specifically comprises:

[0067] the total magnetization inclination of the measuring point is calculated by using the formula ; wherein, I m is the total magnetization inclination; m xis the x-direction magnetic moment component of the measurement point; m y is the y-direction magnetic moment component of the measurement point; m z is the z-direction magnetic moment component of the measurement point.

[0068] The total magnetization inclination of the measurement point is calculated by using the formula m is the total magnetization inclination.

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

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

[0071]

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

[0073]

[0074] wherein, is the total magnetic moment modulus.

[0075] From the formula (1)-(3), it can be seen that infinite integral along the z plane is needed, in order to facilitate the application of the aeromagnetic gridding data, a sliding window is used to calculate the two-dimensional trapezoidal quadrature formula, and the Simpson quadrature formula has higher precision than the trapezoidal quadrature formula under the same calculation amount, therefore, the Simpson quadrature formula is used to calculate the total magnetization direction in the embodiment.

[0076] Since the magnetization direction calculated by the Helbig method at the horizontal position of the magnetic source does not change with the change of the window size, the magnetization direction of each point in the measurement area can be calculated by using the Helbig method through the selection of multiple sliding windows with different sizes, when there is only one magnetic target in the measurement area, the horizontal position of the measurement point with the smallest change in the magnetization direction obtained by multiple calculations is the horizontal position of the magnetic target, and the magnetic inclination and the magnetic declination of the point are the total magnetization direction of the magnetic target, when there are multiple target bodies in the measurement area, the boundary identification method is used for partitioning, and then the magnetization direction of each independent area containing only one target body is estimated.

[0077] A weighted matrix is constructed based on single or multiple components, and a step-by-step weighted inversion is carried out through the weight value to position the underground rock mass. ​

[0078] In view of the fact that multi-component data combination can improve the inversion result, the step-by-step weighted inversion method is applied to the inversion of aeromagnetic three-component data. First, the depth-weighted inversion is performed, then the weight matrix is constructed based on multiple components, and finally the weighted components are used for the next step of inversion calculation. Through weight guidance, the place with large weight participates in the inversion. The information in different directions of each component of aeromagnetic data can be effectively utilized to quickly obtain the inversion result.

[0079] Figure 2 The flow of the step-by-step weighted inversion method is given. Figure 2 The method involves two different weight matrices: one (W2) is constructed from a single component and a depth weighting function W1, and the other (W3,..., W i+1 ) is constructed from multiple components. The overall construction process of the weight matrix is divided into three parts.

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

[0081]

[0082] where the element a ij in the sensitivity matrix A represents the influence of the unit magnetization of cell j (j = 1,..., N) at data position i (i = 1,..., M); β is a constant, and the larger the value, the stronger the weight. Usually, the value of β is close to 1. In the construction process of the weight matrix, only the fitting function is used. 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 shows the method of constructing the weight matrix based on multiple components. After obtaining W2, another component B2 is selected to obtain the inversion result m. The absolute value of the component of m is taken, and m is normalized. Then, a small enough value (10 -6 or less) is assigned to the elements equal to 0, to obtain W3. By analogy, the weight matrix W i based on W i and B i+1 is obtained.

[0086] In the third part, the last component B α in the component combination is selected for inversion. Based on B α and W i+1The objective function, i.e. formula (9), is established, and a nonlinear conjugate gradient algorithm is applied to solve it. Finally, which component is involved in the inversion can be determined according to the actual geological conditions. Since the vertical component alone can be inverted to obtain a result basically consistent with the geological interpretation, the vertical component B z As B α .

[0087]

[0088] Compared with related art, the present application has the following advantages:

[0089] Reducing errors: directly using aeromagnetic three-component data reduces errors in the data conversion process.

[0090] Improving stability and resolution: the step-by-step weighted inversion method improves the stability and resolution of the inversion.

[0091] Quick and accurate information acquisition: the position and direction information of the underground rock mass can be quickly and accurately obtained, providing a new technical means for mineral resource exploration.

[0092] Based on the same inventive concept, the embodiments of the present application also provide an aeromagnetic three-component data-based underground rock mass magnetic target positioning system for implementing the above-mentioned aeromagnetic three-component data-based underground rock mass magnetic target positioning method. The problem-solving implementation scheme provided by the system is similar to the implementation scheme described in the above method, so the specific limitations in the following aeromagnetic three-component data-based underground rock mass magnetic target positioning system embodiments can be referred to the limitations of the aeromagnetic three-component data-based underground rock mass magnetic target positioning method in the above text, which will not be repeated here.

[0093] In one exemplary embodiment, an aeromagnetic three-component data-based underground rock mass magnetic target positioning system is provided, comprising:

[0094] A data acquisition module is configured to acquire aeromagnetic three-component data of multiple measurement points in a measurement region.

[0095] A magnetization direction determination module is configured to determine the magnetization direction of a magnetic target in the measurement region according to the aeromagnetic three-component data of the multiple measurement points in the measurement region; the magnetization direction includes a total magnetization inclination and a total magnetization declination.

[0096] A position determination module is configured to determine the position and magnetic susceptibility of the magnetic target in the measurement region by using a step-by-step weighted inversion method according to the aeromagnetic three-component data.

[0097] In an exemplary embodiment, a computer device is provided, comprising a memory and a processor, the memory storing a computer program, and the processor implementing the above-mentioned method for locating a magnetic target of a subsurface rock mass based on three-component aeromagnetic data when executing the computer program.

[0098] In an exemplary embodiment, a computer readable storage medium is provided, storing a computer program, and the computer program implementing the above-mentioned method for locating a magnetic target of a subsurface rock mass based on three-component aeromagnetic data when executed by a processor.

[0099] In an exemplary embodiment, a computer program product is provided, comprising a computer program, and the computer program implementing the above-mentioned method for locating a magnetic target of a subsurface rock mass based on three-component aeromagnetic data when executed by a processor.

[0100] In an exemplary embodiment, a computer device is provided, which can be a server or a terminal, and its internal structure diagram can be as shown in Figure 3 The computer device comprises a processor, a memory, an input / output interface (I / O) and a communication interface. The processor, the memory and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device comprises a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operating system and the computer program in the non-volatile storage medium to run. The input / output interface of the computer device is used to exchange information between the processor and external devices. The communication interface of the computer device is used to communicate with external terminals through network connection. The computer program is executed by the processor to implement a method for locating a magnetic target of a subsurface rock mass based on three-component aeromagnetic data.

[0101] Those skilled in the art can understand that Figure 3 The structure shown in the above-mentioned figure is only a block diagram of part of the structure related to the scheme of the present application, and does not constitute a limitation on the computer device to which the scheme of the present application is applied. The specific computer device can comprise more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[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 for analysis, stored data, displayed data, etc.) involved in the present application are all information and data authorized by the user or authorized by all parties, and the collection, use and processing of related data need to comply with relevant regulations.

[0103] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer readable storage medium, and when the computer program is executed, the processes of the above-mentioned embodiments of the methods can be included. Any reference to memory, databases or other media used in the embodiments provided in the present 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 storage, high-density embedded non-volatile memory, resistive memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. As an illustration but not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc.

[0104] The database involved in the embodiments provided in the present application can include at least one of a relational database and a non-relational database. The non-relational database can include a distributed database based on a blockchain, etc., without being limited thereto. The processor involved in the embodiments provided in the present application can be a general-purpose processor, a central processing unit, a graphics processing unit, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., without being limited thereto.

[0105] The technical features of the above embodiments can be combined arbitrarily. In order to make the description concise, all possible combinations of the technical features in the above embodiments are not described, but as long as the combinations of the technical features do not exist, they should be considered as the scope of the present application.

[0106] The principles and implementation modes of the present application are described by using specific examples in this paper, and the above-mentioned embodiments are only used to help understand the method and its core idea of the present application. Meanwhile, for those skilled in the art, according to the idea of the present application, the specific implementation mode and application range will be changed. In conclusion, the content of the present application should not be understood as a limitation.

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 magnetic susceptibility of the magnetic target within the measurement area are determined using a stepwise weighted inversion method. A weighted matrix is ​​constructed based on one or more components, and a step-by-step weighted inversion is carried out guided by the weights to locate the underground rock mass; The 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. Through weight guidance, areas with larger weights participate in the inversion. The specific process of the stepwise weighted inversion method is as follows: Rapid inversion initial screening: An automatic inversion method is used, combined with the structural index to estimate the source depth, to provide an initial model for fine inversion and obtain the location (x, y, z) of the underground rock mass; 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. Results fusion: Combining geological data with magnetization direction, location information, and magnetic susceptibility calculated based on aeromagnetic three-component data, a comprehensive interpretation of the measurement area is provided.

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; wherein, The x-direction magnetic moment component of the measurement point; The y-direction magnetic moment component of the measurement point; 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. The x-direction component is the aeromagnetic three-component data of the measurement point; 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; wherein, This is the total magnetization tilt angle; The x-direction magnetic moment component of the measurement point; The y-direction magnetic moment component of the measurement point; The z-direction magnetic moment component of the measurement point; Using formula Calculate the total magnetization deflection at the measurement point; where, 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