A method and device for reducing the pole of complex magnetic anomaly data containing remanent magnetism

Abstract

This application belongs to the field of geophysical magnetic exploration, and specifically discloses a method and device for reducing the pole of complex magnetic anomaly data containing remanent magnetism. The method includes: obtaining geophysical magnetic measurement data of underground magnetic bodies; using the frequency domain processing method to convert the geophysical magnetic measurement data into regularized magnetic source intensity, and determining the horizontal position of the magnetic source center according to the maximum value of the regularized magnetic source intensity; performing integral calculation according to the geophysical magnetic measurement data to obtain the integral result of the magnetization characteristics, and obtaining the global total magnetization direction according to the integral result; obtaining the total magnetization direction of the magnetic body according to the horizontal position of the magnetic source center and the global total magnetization direction; performing pole reduction processing with variable magnetization direction on the geophysical magnetic measurement data according to the total magnetization direction to obtain the pole-reduced field of the complex magnetic anomaly data containing remanent magnetism. Through this application, the accuracy of the total magnetization direction can be improved, and the problem of pole reduction of the remanent magnetism anomaly in the case of multiple field sources can be solved.

Description

Technical Field

[0001] This application belongs to the field of geophysical magnetic prospecting, and more specifically, relates to a method and device for digitizing the poles of complex magnetic anomaly data with remanent magnetism. Background Art

[0002] Magnetic prospecting is an important exploration method in the fields of mineral exploration, geological survey, etc. Underground magnetic bodies generally have remanent magnetism, and their directions are generally different from the contemporary geomagnetic field. In the case of strong remanent magnetism of multiple field sources, the processing and interpretation of magnetic anomaly data will be affected. In actual work, the magnetic anomaly data is usually subjected to pole digitization to eliminate the influence of inclined magnetization. However, if the magnetic anomalies are generated by multiple field sources, the traditional pole digitization means cannot achieve good results, thus affecting subsequent processing and interpretation.

[0003] Currently, in the pole digitization processing of complex magnetic anomalies with remanent magnetism, there are still problems such as large errors in determining the total magnetization direction in the case of multiple magnetic sources, and the inability to perform pole digitization with variable magnetization directions in the case of multiple magnetic sources. Summary of the Invention

[0004] Aiming at the defects of the prior art, the purpose of this application is to provide a method and device for digitizing the poles of complex magnetic anomaly data with remanent magnetism, aiming to solve the problems of large errors in determining the total magnetization direction in the case of multiple magnetic sources and the inability to perform pole digitization with variable magnetization directions in the case of multiple magnetic sources.

[0005] To achieve the above purpose, in the first aspect, this application provides a method for digitizing the poles of complex magnetic anomaly data with remanent magnetism, including:

[0006] Obtain the geophysical magnetic measurement data of the underground magnetic body;

[0007] Use the frequency domain processing method to convert the geophysical magnetic measurement data into regularized magnetic source intensity, and determine the horizontal position of the magnetic source center according to the maximum value of the regularized magnetic source intensity;

[0008] Perform integral calculation according to the geophysical magnetic measurement data to obtain the integral result of the magnetization characteristics, and obtain the global total magnetization direction according to the integral result;

[0009] Obtain the total magnetization direction of the magnetic body according to the horizontal position of the magnetic source center and the global total magnetization direction;

[0010] Perform pole digitization with variable magnetization directions on the geophysical magnetic measurement data according to the total magnetization direction to obtain the pole digitized field of the complex magnetic anomaly data with remanent magnetism.

[0011] Optionally, the method for determining the regularized magnetic source intensity includes:

[0012] Calculate the magnetic gradient tensor matrix based on the geophysical magnetic measurement data; the magnetic gradient tensor matrix is a symmetric matrix;

[0013] Determine the eigenvalue matrix of the symmetric matrix and the eigenvectors corresponding to the eigenvectors, and arrange the multiple eigenvalues of the eigenvalue matrix in a monotonically decreasing order;

[0014] Determine the regularized magnetic source intensity according to multiple eigenvalues.

[0015] Optionally, the method for obtaining the magnetic gradient tensor matrix includes:

[0016] Perform Fourier transform on the geophysical magnetic measurement data, determine the transverse wave numbers and longitudinal wave numbers corresponding to the x-direction and y-direction of the geophysical magnetic measurement data, and determine the geomagnetic inclination and geomagnetic declination of the geophysical magnetic measurement data;

[0017] Determine the first conversion factor in the frequency domain direction according to the transverse wave number, longitudinal wave number, geomagnetic inclination and geomagnetic declination;

[0018] Calculate the three magnetic anomaly components in the x-direction, y-direction and z-direction respectively according to the geophysical magnetic measurement data, the first conversion factor in the frequency domain direction, the transverse wave number and the longitudinal wave number;

[0019] Perform partial derivative calculation according to the three magnetic anomaly components in the x-direction, y-direction and z-direction, and construct a magnetic gradient tensor matrix according to the partial derivative calculation results.

[0020] Optionally, the integration result includes the x-integration corresponding to the x-direction, the y-integration corresponding to the y-direction, and the z-integration corresponding to the z-direction;

[0021] Obtain the global total magnetization direction according to the integration result, including:

[0022] Perform arcsine calculation according to the x-integration, y-integration and z-integration to obtain the total magnetization inclination;

[0023] Judge according to the x-integration and y-integration, and calculate the total magnetization declination according to the judgment result in combination with the x-integration and y-integration.

[0024] Optionally, the method for performing variable magnetization direction reduction on the geophysical magnetic measurement data according to the total magnetization direction to obtain the reduction field of the complex magnetic anomaly data containing remanence includes:

[0025] Determine the average inclination according to the total magnetization inclination, and determine the average declination according to the total magnetization declination, and respectively determine the inclination difference and declination difference;

[0026] Determine the second conversion factor in the frequency domain direction according to the horizontal wave number, vertical wave number, average dip angle, and average declination angle;

[0027] Perform two-dimensional Fourier transform and inverse transform on the geophysical magnetic survey data, horizontal wave number, vertical wave number, first conversion factor, and second conversion factor to obtain the average magnetic field data;

[0028] Calculate the partial derivatives of the average dip angle and average declination angle respectively based on the average magnetic field data to obtain the dip angle partial derivative and declination angle partial derivative;

[0029] Perform reduction to the pole processing with variable magnetization directions based on the average magnetic field data, dip angle difference, dip angle partial derivative, declination angle difference, and declination angle partial derivative to obtain the reduced-to-pole field of the complex magnetic anomaly data containing remanence.

[0030] Optionally, the calculating the total magnetization declination by combining the x-integral and y-integral according to the judgment result includes:

[0031] If the x-integral is greater than or equal to 0, perform arctangent calculation based on the ratio of the x-integral and y-integral to obtain the total magnetization declination;

[0032] If the x-integral is less than 0 and the y-integral is greater than or equal to 0, perform arctangent calculation based on the ratio of the y-integral and x-integral, and add π to the calculation result to obtain the total magnetization declination;

[0033] If the x-integral is less than 0 and the y-integral is less than 0, perform arctangent calculation based on the ratio of the y-integral and x-integral, and subtract π from the calculation result to obtain the total magnetization declination.

[0034] Optionally, the obtaining the total magnetization direction of the magnetic body according to the horizontal position of the magnetic source center and the global total magnetization direction includes:

[0035] According to the horizontal position of the magnetic source center, select the total magnetization dip angle and total magnetization declination corresponding to the global total magnetization direction, and obtain the total magnetization direction according to the total magnetization dip angle and total magnetization declination.

[0036] In a second aspect, the present application provides a device for reducing to the pole of complex magnetic anomaly data containing remanence, including:

[0037] A data acquisition module for acquiring geophysical magnetic survey data of underground magnetic bodies;

[0038] A magnetic source center determination module for converting the geophysical magnetic survey data into a regularized magnetic source intensity by using the frequency domain processing method, and determining the horizontal position of the magnetic source center according to the maximum value of the regularized magnetic source intensity;

[0039] A global total magnetization direction module for obtaining the global total magnetization direction according to the integration result, performing integration calculation on the geophysical magnetic measurement data to obtain the integration result of the magnetization characteristics, and obtaining the global total magnetization direction according to the integration result;

[0040] A total magnetization direction module for obtaining the total magnetization direction of the magnetic body according to the horizontal position of the magnetic source center and the global total magnetization direction;

[0041] A pole conversion module for performing pole conversion processing on the geophysical magnetic measurement data according to the total magnetization direction to obtain the pole conversion field of the complex magnetic anomaly data containing remanence.

[0042] In a third aspect, the present application provides an electronic device, including: at least one memory for storing a program; at least one processor for executing the program stored in the memory. When the program stored in the memory is executed, the processor is used to execute the method described in the first aspect or any one of the possible implementation manners of the first aspect.

[0043] In a fourth aspect, the present application provides a computer-readable storage medium storing a computer program. When the computer program runs on a processor, the processor is caused to execute the method described in the first aspect or any one of the possible implementation manners of the first aspect.

[0044] In a fifth aspect, the present application provides a computer program product. When the computer program product runs on a processor, the processor is caused to execute the method described in the first aspect or any one of the possible implementation manners of the first aspect.

[0045] It can be understood that the beneficial effects of the above second aspect to fifth aspect can refer to the relevant descriptions in the above first aspect and will not be repeated here.

[0046] Generally speaking, compared with the prior art, the above technical solutions conceived by the present application have the following beneficial effects:

[0047] (1) By considering the total magnetization directions of multiple field sources, the present application can more accurately reflect the magnetic field characteristics in a complex geological environment. By comprehensively analyzing the regularized magnetic source intensity and the Helbig integral, the accuracy of calculating the total magnetization direction is effectively improved. By establishing a total magnetization direction model and considering the change of the magnetization direction during pole conversion processing, the accuracy of calculating the total magnetization direction by the Helbig integral is improved, and the error caused by remanence is reduced. Therefore, it meets the requirements under complex geological conditions and solves the pole conversion problem of the remanence anomaly in the case of multiple magnetic sources.

[0048] (2) This application uses the Helbig integral to estimate the total magnetization direction of multiple field sources, and uses the variable magnetization direction reduction-to-pole technique to obtain the reduction-to-pole field of non-isolated field sources, providing data support for the processing and interpretation of magnetic anomaly data in the case of remanence. Description of the Drawings

[0049] Figure 1 is one of the schematic flowcharts of the reduction-to-pole method for complex magnetic anomaly data with remanence provided by an embodiment of this application;

[0050] Figure 2 is a schematic diagram of the underground magnetic body model in an embodiment of this application;

[0051] Figure 3 are the geophysical magnetic measurement data in an embodiment of this application and the reduction-to-pole field in the theoretical case;

[0052] Figure 4 are the normalized source strength NSS (Normalized Source Strength) calculated according to the geophysical magnetic measurement data and the horizontal position of the magnetic source center in an embodiment of this application;

[0053] Figure 5 are the global total magnetization inclination angle and total magnetization declination angle calculated according to the geophysical magnetic measurement data in an embodiment of this application;

[0054] Figure 6 are the geomagnetic direction reduction-to-pole field of the experimental model and the variable inclination reduction-to-pole field considering remanence in an embodiment of this application;

[0055] Figure 7 is the second schematic flowchart of the reduction-to-pole method for complex magnetic anomaly data with remanence provided by an embodiment of this application;

[0056] Figure 8 is the structural schematic diagram of the reduction-to-pole device for complex magnetic anomaly data with remanence provided by an embodiment of this application;

[0057] Figure 9 is the structural schematic diagram of the electronic device provided by an embodiment of this application. Detailed Embodiments

[0058] In order to make the objectives, technical solutions and advantages of this application clearer, the following further elaborates on this application in combination with the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain this application and are not used to limit this application.

[0059] In this text, the term "and / or" describes the relationship between associated objects, indicating three possible relationships. For example, A and / or B can represent three cases: A exists alone, A and B exist simultaneously, and B exists alone. In this text, the symbol " / " indicates that the associated objects are in an "or" relationship. For example, A / B means A or B.

[0060] In the description of the specification and claims of this application, terms such as "first" and "second" are used to distinguish different objects, rather than to describe a specific order of the objects. For example, the first response message and the second response message are used to distinguish different response messages, rather than to describe a specific order of the response messages.

[0061] In the embodiments of this application, words such as "exemplary" or "for example" are used to indicate examples, illustrations, or explanations. Any embodiment or design solution described as "exemplary" or "for example" in the embodiments of this application should not be construed as being more preferred or having more advantages than other embodiments or design solutions. Rather, the use of words such as "exemplary" or "for example" is intended to present relevant concepts in a specific manner.

[0062] In the description of the embodiments of this application, unless otherwise specified, the meaning of "a plurality of" refers to two or more. For example, a plurality of processing units refers to two or more processing units, etc.; a plurality of elements refers to two or more elements, etc.

[0063] The embodiments of this application will be described below with reference to the accompanying drawings in the embodiments of this application.

[0064] Refer to Figure 1 , this application provides a method for reducing the pole of complex magnetic anomaly data with remanence, including:

[0065] S101. Obtain the geophysical magnetic measurement data of the underground magnetic body;

[0066] S102. Use the frequency domain processing method to convert the geophysical magnetic measurement data into a regularized magnetic source intensity, and determine the horizontal position of the magnetic source center according to the maximum value of the regularized magnetic source intensity;

[0067] S103. Integrate the geophysical magnetic measurement data to obtain the integration result of the magnetization characteristics, and obtain the total magnetization direction of the whole domain according to the integration result;

[0068] S104. Obtain the total magnetization direction of the magnetic body according to the horizontal position of the magnetic source center and the total magnetization direction of the whole domain;

[0069] S105. Perform pole reduction processing on the geophysical magnetic measurement data in the variable magnetization direction according to the total magnetization direction to obtain the pole reduction field of the complex magnetic anomaly data with remanence.

[0070] Specifically, in this embodiment S101, relevant magnetic measurement data of the underground magnetic body needs to be obtained through geophysical exploration technology first.

[0071] The embodiment of the present application adopts the analytical forward modeling method to simulate the geophysical magnetic measurement data generated by the underground magnetic body, such as Figure 2 shown Figure 2 is the underground magnetic body model. The analytical forward formula of a regular hexahedron is used for forward modeling. The central burial depths of the four underground magnetic bodies used in the simulation experiment are all 300 meters, the shape is a cube, the length, width and height are all 200 meters, the total magnetization intensity is 1 A / m, the geomagnetic inclination is 45°, the geomagnetic declination is 0°, the total magnetization inclinations of the four magnetic bodies are 30°, 60°, -30° and -60° respectively, and the total magnetization declinations of the four magnetic bodies are 30°, 60°, -30° and -60° respectively. The ground measurement point grid size is 41×41, the north-south spacing of the grid is 50 meters, and the entire measurement area range is 2000m×2000m.

[0072] It should be noted that in this drawing and subsequent drawings, Easting represents the eastward direction, Northing represents the northward direction, and Depth represents the depth.

[0073] In the embodiment of the present application, the x direction: refers to the east-west direction (east - west direction). The y direction: usually refers to the north-south direction (south - north direction). The z direction: refers to the depth direction, which is the direction perpendicular to the ground.

[0074] Refer to Figure 3 , Figure 3 in which (a) is the magnetic anomaly generated by the model, which is equivalent to the geophysical magnetic measurement data in actual work, Figure 3 in which (b) represents the reduced-to-pole field in the theoretical case, that is, the expected output.

[0075] Then through S102, the geophysical magnetic measurement data is converted into regularized magnetic source intensity by using the frequency domain processing method NSS .

[0076] In this step, the frequency domain processing method is used to convert the obtained geophysical magnetic measurement data to obtain the regularized magnetic source intensity. This process involves performing a Fourier transform on the magnetic measurement data to transform it into the frequency domain. During the processing, by selecting appropriate regularization parameters, the noise and irregularities in the data are smoothed out to ensure the stability and reliability of the result. Finally, based on the regularized magnetic source intensity, the central position of the magnetic source is further determined. In this embodiment, NSS the position of the maximum value is selected as the magnetic source center position.

[0077] Further, through S103, integral calculation is performed based on the geophysical magnetic measurement data, namely the Helbig integral. This step uses the Helbig integral to analyze the geophysical magnetic measurement data to obtain the total magnetization direction of the entire region.

[0078] It should be noted that the Helbig integral is a mathematical method widely used in magnetic field analysis. By performing integral calculation on the obtained magnetic measurement data, magnetic field effects from multiple sources can be identified.

[0079] Through S104, the total magnetization direction of the magnetic body is obtained based on the horizontal position of the magnetic source center and the total magnetization direction of the entire region. Based on the horizontal position of the magnetic source center determined in the previous step and the total magnetization direction of the entire region obtained from the Helbig integral. According to the magnetic source position and its magnetization characteristics, a comprehensive description of the magnetization direction is formed. The total magnetization direction model can reflect the overall magnetization state.

[0080] Finally, through S105, the pole transformation of the geophysical magnetic measurement data with variable magnetization direction is implemented by adopting the previously obtained total magnetization direction. This processing mainly involves adjusting the measured magnetic field data to a form consistent with the total magnetization direction, thereby reducing the interference caused by remanence anomalies.

[0081] Optionally, the method for determining the regularized magnetic source strength includes:

[0082] Calculating a magnetic gradient tensor matrix based on the geophysical magnetic measurement data; the magnetic gradient tensor matrix is a symmetric matrix;

[0083] Determining the eigenvalue matrix of the symmetric matrix and the eigenvectors corresponding to the eigenvectors, and arranging the multiple eigenvalues of the eigenvalue matrix in a monotonically decreasing order;

[0084] Determining the regularized magnetic source strength based on the multiple eigenvalues.

[0085] Further, the method for obtaining the magnetic gradient tensor matrix includes:

[0086] Performing Fourier transform on the geophysical magnetic measurement data, determining the transverse wave numbers and longitudinal wave numbers corresponding to the x - direction and y - direction of the geophysical magnetic measurement data, and determining the geomagnetic inclination and geomagnetic declination of the geophysical magnetic measurement data;

[0087] Determining a first conversion factor in the frequency - domain direction based on the transverse wave number, longitudinal wave number, geomagnetic inclination, and geomagnetic declination;

[0088] Calculating the three magnetic anomaly components in the x - direction, y - direction, and z - direction respectively based on the geophysical magnetic measurement data, the first conversion factor in the frequency - domain direction, the transverse wave number, and the longitudinal wave number;

[0089] Calculate the partial derivatives based on the three magnetic anomaly components in the x, y, and z directions, and construct a magnetic gradient tensor matrix according to the partial derivative calculation results.

[0090] Specifically, this embodiment calculates the regularized magnetic source intensity in the frequency domain NSS The specific process is as follows:

[0091] Calculate the regularized magnetic source intensity in the frequency domain based on the geophysical magnetic survey data NSS , as the first group of data, as Figure 4 shown, the data scale is a 41×41 real matrix, and the calculation process is as follows:

[0092] First, calculate the magnetic gradient tensor from the geophysical magnetic survey data

[0093]

[0094] According to the potential field theory, there is , so this gradient tensor is a symmetric matrix and only contains six independent components. Therefore can be expressed as:

[0095]

[0096]

[0097] Among them, is the eigenvalue of the symmetric matrix , is the eigenvector corresponding to the eigenvalue . Arrange the eigenvalues of the magnetic gradient tensor in a monotonically decreasing order: , then the regularized magnetic source intensity NSS is expressed as:

[0098]

[0099] Finally, select NSS the position of the maximum value as the horizontal position of the magnetic source center, as Figure 3 the position of the "white cross" in

[0100] Optionally, the integration results include the x-integration corresponding to the x direction, the y-integration corresponding to the y direction, and the z-integration corresponding to the z direction;

[0101] Obtain the global total magnetization direction according to the integration results, including:

[0102] Perform an arcsine calculation based on the x-integration, y-integration, and z-integration to obtain the total magnetization dip angle;

[0103] Judgment is made according to the x-integral and the y-integral, and the total magnetization declination is calculated by combining the x-integral and the y-integral according to the judgment result.

[0104] Furthermore, the calculating the total magnetization declination by combining the x-integral and the y-integral according to the judgment result includes:

[0105] If the x-integral is greater than or equal to 0, arctangent calculation is performed according to the ratio of the x-integral and the y-integral to obtain the total magnetization declination;

[0106] If the x-integral is less than 0 and the y-integral is greater than or equal to 0, arctangent calculation is performed according to the ratio of the y-integral and the x-integral, and π is added to the calculation result to obtain the total magnetization declination;

[0107] If the x-integral is less than 0 and the y-integral is less than 0, arctangent calculation is performed according to the ratio of the y-integral and the x-integral, and π is subtracted from the calculation result to obtain the total magnetization declination.

[0108] Specifically, in this embodiment, according to the geophysical magnetic survey data calculate the global magnetization, and the steps are as follows:

[0109] First, the calculation formula of the Helbig integral is as follows:

[0110]

[0111] Among them, are the three components of the magnetic anomaly calculated from the geophysical magnetic survey data The global total magnetization direction is obtained by the following formula:

[0112] ,

[0113] ,

[0114] Among them, and represent the total magnetization inclination angle and the total magnetization declination angle in the global total magnetization direction respectively. As Figure 5 shown, Figure 5 The abscissa of (a) in Figure 5 and the abscissa of (b) in

[0115] are both eastward, and the ordinates are both northward.

[0116] Optionally, the obtaining the total magnetization direction of the magnetic body according to the horizontal position of the magnetic source center and the global total magnetization direction includes:

[0117] Select the total magnetization inclination angle and the total magnetization declination angle corresponding to the global total magnetization direction according to the horizontal position of the magnetic source center, and obtain the total magnetization direction according to the total magnetization inclination angle and the total magnetization declination angle.According to the magnetic source center position obtained in the above embodiments, select the total magnetization dip angle and the total magnetization declination corresponding to the global total magnetization direction as the total magnetization direction of the magnetic source. The estimation errors are shown in Table 1.

[0118] Table 1 Estimation Errors of Total Magnetization Direction

[0119]

[0120] Optionally, the processing of polarizing the geophysical magnetic measurement data in a variable magnetization direction according to the total magnetization direction to obtain a polarized field of complex magnetic anomalies containing remanent magnetization includes:

[0121] Determine the average dip angle according to the total magnetization dip angle, and determine the average declination according to the total magnetization declination, and respectively determine the dip angle difference and the declination difference;

[0122] Determine the second conversion factor in the frequency domain direction according to the transverse wave number, the longitudinal wave number, the average dip angle, and the average declination;

[0123] Perform two-dimensional Fourier transform and inverse transform on the geophysical magnetic measurement data, the transverse wave number, the longitudinal wave number, the first conversion factor, and the second conversion factor to obtain the mean magnetic field data;

[0124] Calculate the partial derivatives of the average dip angle and the average declination respectively according to the mean magnetic field data to obtain the dip angle partial derivative and the declination partial derivative;

[0125] Perform variable magnetization direction polarizing processing on the mean magnetic field data, the dip angle difference, the dip angle partial derivative, the declination difference, and the declination partial derivative to obtain a polarized field of complex magnetic anomalies containing remanent magnetization.

[0126] Specifically, according to the total magnetization direction of the magnetic source, the geophysical magnetic measurement data is processed by variable magnetization direction polarizing to obtain a polarized field , and the calculation formula is as follows:

[0127]

[0128] Among them,

[0129]

[0130]

[0131]

[0132] represents the geophysical magnetic measurement data, respectively represent two-dimensional Fourier transform and its inverse transform, is the imaginary unit, After Fourier transform and the wavenumbers in the 、 are direction conversion factors in the frequency domain, i.e., the first conversion factor and the second conversion factor. are the geomagnetic inclination and declination respectively, i.e., is the subscript of the direction indicated. are the average values of the obtained total magnetization inclination and total magnetization declination respectively, i.e., is the subscript of the direction indicated. and represent the differences between the total magnetization inclination and total magnetization declination and their average values respectively.

[0133] Referring to Figure 6 , when the remanent magnetization is not considered, the reduced-to-pole field obtained is as shown in (a) of Figure 6 . The high value of the magnetic anomaly does not correspond to the magnetic source center. The reduced-to-pole field with variable magnetization direction considering the remanent magnetization is as shown in (b) of Figure 6 . The inclined magnetization is effectively eliminated, and the magnetic source center corresponds to the high value area of the magnetic anomaly.

[0134] It should be noted that the above calculation formula takes the geomagnetic direction and the total magnetization direction as inputs. The geomagnetic direction can be obtained by querying the International Geomagnetic Reference Field (IGRF) in actual work, and the input of the total magnetization direction is to select a total magnetization direction within the set range of the total magnetization direction.

[0135] Referring to Figure 7 , Figure 7 The embodiment of the present application provides a complete flow schematic diagram of a method for reducing complex magnetic anomaly data with remanent magnetization, including the following steps:

[0136] Input geophysical magnetic measurement data;

[0137] Calculate the Helbig integral to obtain the global total magnetization direction;

[0138] Calculate the regularized magnetic source strength NSS Obtain the horizontal position of the magnetic source center;

[0139] Construct a total magnetization direction model;

[0140] Perform reduced-to-pole with variable magnetization direction on the geophysical magnetic measurement data;

[0141] Input geophysical magnetic measurement data.

[0142] Referring to Figure 8 , the present application provides a device for reducing complex magnetic anomaly data with remanent magnetization, including:

[0143] A data acquisition module 810, configured to acquire geophysical magnetic measurement data of underground magnetic bodies;

[0144] A magnetic source center determination module 820, configured to convert the geophysical magnetic measurement data into a regularized magnetic source intensity by using a frequency domain processing method, and determine the horizontal position of the magnetic source center according to the maximum value of the regularized magnetic source intensity;

[0145] A global total magnetization direction module 830, configured to perform integral calculation according to the geophysical magnetic measurement data to obtain an integral result of magnetization characteristics, and obtain the global total magnetization direction according to the integral result;

[0146] A total magnetization direction module 840, configured to obtain the total magnetization direction of the magnetic body according to the horizontal position of the magnetic source center and the global total magnetization direction;

[0147] A polarizing module 850, configured to perform polarizing processing of variable magnetization direction on the geophysical magnetic measurement data according to the total magnetization direction to obtain a polarized field of complex magnetic anomalies containing remanent magnetism.

[0148] Optionally, the method for determining the regularized magnetic source intensity includes:

[0149] Calculating a magnetic gradient tensor matrix according to the geophysical magnetic measurement data; the magnetic gradient tensor matrix is a symmetric matrix;

[0150] Determining an eigenvalue matrix of the symmetric matrix and eigenvectors corresponding to the eigenvectors, and arranging multiple eigenvalues of the eigenvalue matrix in a monotonically decreasing order;

[0151] Determining the regularized magnetic source intensity according to multiple eigenvalues.

[0152] Optionally, the method for obtaining the magnetic gradient tensor matrix includes:

[0153] Performing Fourier transform on the geophysical magnetic measurement data, determining transverse wave numbers and longitudinal wave numbers corresponding to the x direction and y direction of the geophysical magnetic measurement data, and determining the geomagnetic inclination and geomagnetic declination of the geophysical magnetic measurement data;

[0154] Determining a first conversion factor in the frequency domain direction according to the transverse wave number, longitudinal wave number, geomagnetic inclination, and geomagnetic declination;

[0155] Calculating magnetic anomaly three components in the x direction, y direction, and z direction respectively according to the geophysical magnetic measurement data, the first conversion factor in the frequency domain direction, the transverse wave number, and the longitudinal wave number;

[0156] Perform partial derivative calculations based on the three magnetic anomaly components in the x, y, and z directions, and construct a magnetic gradient tensor matrix according to the partial derivative calculation results.

[0157] Optionally, the integration results include the x-integration corresponding to the x direction, the y-integration corresponding to the y direction, and the z-integration corresponding to the z direction;

[0158] Obtain the global total magnetization direction based on the integration results, including:

[0159] Perform arcsine calculations based on the x-integration, y-integration, and z-integration to obtain the total magnetization inclination;

[0160] Judge based on the x-integration and y-integration, and calculate the total magnetization declination according to the judgment result in combination with the x-integration and y-integration.

[0161] Optionally, the variable magnetization direction reduction processing of the geophysical magnetic measurement data according to the total magnetization direction to obtain the reduction field of the complex magnetic anomaly data containing remanent magnetization includes:

[0162] Determine the average inclination angle according to the total magnetization inclination angle, and determine the average declination angle according to the total magnetization declination angle, and respectively determine the inclination angle difference and declination angle difference;

[0163] Determine the second conversion factor in the frequency domain direction according to the transverse wave number, longitudinal wave number, average inclination angle, and average declination angle;

[0164] Perform two-dimensional Fourier transform and inverse transform on the geophysical magnetic measurement data, transverse wave number, longitudinal wave number, first conversion factor, and second conversion factor to obtain the mean magnetic field data;

[0165] Perform partial derivative calculations on the average inclination angle and average declination angle respectively according to the mean magnetic field data to obtain the inclination angle partial derivative and declination angle partial derivative;

[0166] Perform variable magnetization direction reduction processing according to the mean magnetic field data, inclination angle difference, inclination angle partial derivative, declination angle difference, and declination angle partial derivative to obtain the reduction field of the complex magnetic anomaly data containing remanent magnetization.

[0167] Optionally, the calculation of the total magnetization declination according to the judgment result in combination with the x-integration and y-integration includes:

[0168] If the x-integration is greater than or equal to 0, perform arctangent calculation according to the ratio of the x-integration and y-integration to obtain the total magnetization declination;

[0169] If the x-integration is less than 0 and the y-integration is greater than or equal to 0, perform arctangent calculation according to the ratio of the y-integration and x-integration, and add π to the calculation result to obtain the total magnetization declination;

[0170] If the x-integral is less than 0 and the y-integral is less than 0, perform an arctangent calculation based on the ratio of the y-integral to the x-integral, and subtract π from the calculation result to obtain the total magnetization declination.

[0171] Optionally, obtaining the total magnetization direction of the magnetic body according to the horizontal position of the magnetic source center and the global total magnetization direction includes:

[0172] According to the horizontal position of the magnetic source center, select the total magnetization inclination angle and the total magnetization declination corresponding to the global total magnetization direction, and obtain the total magnetization direction according to the total magnetization inclination angle and the total magnetization declination.

[0173] It can be understood that the detailed function implementation of each of the above units / modules can be referred to the introduction in the foregoing method embodiments, and will not be elaborated here.

[0174] It should be understood that the above device is used to execute the method in the above embodiments. For the corresponding program modules in the device, the implementation principle and technical effects are similar to the description in the above method. The working process of the device can refer to the corresponding process in the above method, and will not be elaborated here.

[0175] Refer to Figure 9 , based on the method in the above embodiments, an embodiment of the present application provides an electronic device, which may include: a processor (Processor) 910, a communication interface (Communications Interface) 920, a memory (Memory) 930, and a communication bus 940. Among them, the processor 910, the communication interface 920, and the memory 930 complete mutual communication through the communication bus 940. The processor 910 can call the logical instructions in the memory 930 to execute the method in the above embodiments.

[0176] In addition, when the logical instructions in the above memory 930 are implemented in the form of a software functional unit and sold or used as an independent product, they can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present application, in essence, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present application.

[0177] Based on the method in the above embodiments, an embodiment of the present application provides a computer-readable storage medium. The computer-readable storage medium stores a computer program, and when the computer program runs on a processor, it causes the processor to execute the method in the above embodiments.

[0178] Based on the method in the above embodiments, an embodiment of the present application provides a computer program product. When the computer program product runs on a processor, it causes the processor to execute the method in the above embodiments.

[0179] It can be understood that the processor in the embodiments of the present application may be a central processing unit (CPU), or may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field programmable gate arrays (FPGAs), or other programmable logic devices, transistor logic devices, hardware components, or any combination thereof. The general-purpose processor may be a microprocessor or any conventional processor.

[0180] The method steps in the embodiments of the present application may be implemented in a hardware manner or by a processor executing software instructions. The software instructions may be composed of corresponding software modules, and the software modules may be stored in a random access memory (RAM), flash memory, read-only memory (ROM), programmable ROM (PROM), erasable PROM (EPROM), electrically erasable PROM (EEPROM), registers, hard disks, removable hard disks, CD-ROMs, or any other form of storage medium well known in the art. An exemplary storage medium is coupled to the processor so that the processor can read information from the storage medium and write information to the storage medium. Of course, the storage medium may also be a component of the processor. The processor and the storage medium may be located in an ASIC.

[0181] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented using software, it can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the processes or functions described in the embodiments of the present application are generated in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions can be stored in a computer-readable storage medium or transmitted through the computer-readable storage medium. The computer instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center in a wired manner (such as coaxial cable, optical fiber, digital subscriber line (DSL)) or wirelessly (such as infrared, wireless, microwave, etc.). The computer-readable storage medium can be any available medium that can be accessed by a computer or a data storage device such as a server or a data center that includes one or more integrated available media. The available medium can be a magnetic medium (such as a floppy disk, a hard disk, a magnetic tape), an optical medium (such as a DVD), or a semiconductor medium (such as a solid state disk (SSD)), etc.

[0182] It can be understood that the various numerical numbers involved in the embodiments of the present application are only for the convenience of description and are not used to limit the scope of the embodiments of the present application.

[0183] Those skilled in the art can easily understand that the above are only the preferred embodiments of the present application and are not intended to limit the present application. Any modifications, equivalent replacements, and improvements made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method for reducing the pole of complex magnetic anomaly data containing remanent magnetism, characterized in that, Including: Obtaining geophysical magnetic measurement data of underground magnetic bodies; Using the frequency domain processing method to convert the geophysical magnetic measurement data into regularized magnetic source intensity, and determining the horizontal position of the magnetic source center according to the maximum value of the regularized magnetic source intensity; Performing integral calculation on the geophysical magnetic measurement data to obtain an integral result of magnetization characteristics, and obtaining the global total magnetization direction according to the integral result; Obtaining the total magnetization direction of the magnetic body according to the horizontal position of the magnetic source center and the global total magnetization direction; Performing variable magnetization direction reduction processing on the geophysical magnetic measurement data according to the total magnetization direction to obtain a reduction field of complex magnetic anomaly data containing remanence.

2. The method for reducing the pole of complex magnetic anomaly data with remanent magnetism according to claim 1, wherein, The method for determining the regularized magnetic source intensity includes: Calculating a magnetic gradient tensor matrix according to the geophysical magnetic measurement data; the magnetic gradient tensor matrix is a symmetric matrix; Determining the eigenvalue matrix of the symmetric matrix and the eigenvectors corresponding to the eigenvectors, and arranging the multiple eigenvalues of the eigenvalue matrix in a monotonically decreasing order; Determining the regularized magnetic source intensity according to multiple eigenvalues.

3. The method for reducing the pole of complex magnetic anomaly data with remanence according to claim 2, characterized in that, The method for obtaining the magnetic gradient tensor matrix includes: Performing Fourier transform on the geophysical magnetic measurement data, determining the transverse wave number and longitudinal wave number corresponding to the x direction and y direction of the geophysical magnetic measurement data, and determining the geomagnetic inclination and geomagnetic declination of the geophysical magnetic measurement data; Determining a first conversion factor in the frequency domain direction according to the transverse wave number, longitudinal wave number, geomagnetic inclination, and geomagnetic declination; Calculating the three magnetic anomaly components in the x direction, y direction, and z direction respectively according to the geophysical magnetic measurement data, the first conversion factor in the frequency domain direction, the transverse wave number, and the longitudinal wave number; Performing partial derivative calculation according to the three magnetic anomaly components in the x direction, y direction, and z direction, and constructing a magnetic gradient tensor matrix according to the partial derivative calculation result.

4. The method for reducing the pole of complex magnetic anomaly data with remanent magnetism according to claim 3, wherein The integral result includes the x integral corresponding to the x direction, the y integral corresponding to the y direction, and the z integral corresponding to the z direction; Obtaining the global total magnetization direction according to the integral result includes: Performing arcsine calculation according to the x integral, y integral, and z integral to obtain the total magnetization dip angle; Judging according to the x integral and y integral, and calculating the total magnetization declination according to the judgment result in combination with the x integral and y integral.

5. The method for reducing the pole of complex magnetic anomaly data with remanence according to claim 4, characterized in that Performing variable magnetization direction reduction processing on the geophysical magnetic measurement data according to the total magnetization direction to obtain a reduction field of complex magnetic anomaly data containing remanence includes: Determining the average dip angle according to the total magnetization dip angle, and determining the average declination angle according to the total magnetization declination angle, and respectively determining the dip angle difference and declination angle difference; Determining a second conversion factor in the frequency domain direction according to the transverse wave number, longitudinal wave number, average dip angle, and average declination angle; Performing two-dimensional Fourier transform and inverse transform on the geophysical magnetic measurement data, the transverse wave number, longitudinal wave number, first conversion factor, and second conversion factor to obtain the magnetic field data mean value; Performing partial derivative calculation on the average dip angle and average declination angle respectively according to the magnetic field data mean value to obtain the dip angle partial derivative and declination angle partial derivative; Perform variable magnetization direction reduction on the basis of the mean value of the magnetic field data, the dip difference, the dip partial derivative, the declination difference, and the declination partial derivative to obtain the reduced field of complex magnetic anomalies with remanent magnetization.

6. The method for reducing the pole of complex magnetic anomaly data with remanent magnetism according to claim 4, characterized in that The calculating the total magnetization declination according to the judgment result in combination with the x-integral and the y-integral includes: If the x-integral is greater than or equal to 0, perform arctangent calculation according to the ratio of the x-integral to the y-integral to obtain the total magnetization declination; If the x-integral is less than 0 and the y-integral is greater than or equal to 0, perform arctangent calculation according to the ratio of the y-integral to the x-integral, and add π to the calculation result to obtain the total magnetization declination; If the x-integral is less than 0 and the y-integral is less than 0, perform arctangent calculation according to the ratio of the y-integral to the x-integral, and subtract π from the calculation result to obtain the total magnetization declination.

7. The method for reducing the pole of complex magnetic anomaly data with remanent magnetism according to claim 1, characterized in that The obtaining the total magnetization direction of the magnetic body according to the horizontal position of the magnetic source center and the global total magnetization direction includes: According to the horizontal position of the magnetic source center, select the total magnetization inclination and the total magnetization declination corresponding to the global total magnetization direction, and obtain the total magnetization direction according to the total magnetization inclination and the total magnetization declination.

8. A reduction to the pole device for complex magnetic anomaly data containing remanent magnetism, characterized in that, Including: A data acquisition module for acquiring geophysical magnetic measurement data of underground magnetic bodies; A magnetic source center determination module for converting the geophysical magnetic measurement data into a regularized magnetic source intensity by using a frequency domain processing method, and determining the horizontal position of the magnetic source center according to the maximum value of the regularized magnetic source intensity; A global total magnetization direction module for performing integral calculation according to the geophysical magnetic measurement data to obtain an integral result of magnetization characteristics, and obtaining the global total magnetization direction according to the integral result; A total magnetization direction module for obtaining the total magnetization direction of the magnetic body according to the horizontal position of the magnetic source center and the global total magnetization direction; A reduction module for performing variable magnetization direction reduction on the geophysical magnetic measurement data according to the total magnetization direction to obtain the reduced field of complex magnetic anomalies with remanent magnetization.

9. An electronic device, characterized in that, Including: At least one memory for storing a computer program; At least one processor for executing the program stored in the memory, and when the program stored in the memory is executed, the processor is used to execute the method according to any one of claims 1-7.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program runs on the processor, the processor is caused to execute the method according to any one of claims 1-7.

Images