Polarization method and device for complex magnetic anomaly data containing residual magnetism

Abstract

The invention belongs to the field of geophysical magnetic prospecting, and particularly discloses a remanence-containing complex magnetic anomaly data polarization method and device, and the method comprises the steps: obtaining the geophysical magnetic survey data of an underground magnetic body; converting the geophysical magnetic measurement data into regularized magnetic source intensity by using a frequency domain processing method, and determining the horizontal position of the center of a magnetic source according to the maximum value of the regularized magnetic source intensity; performing integral calculation according to the geophysical magnetic measurement data to obtain an integral result of magnetization characteristics, and obtaining a 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; and according to the total magnetization direction, carrying out magnetization direction variable polarity processing on the geophysical magnetic measurement data to obtain a polarity field of complex magnetic anomaly data containing residual magnetism. According to the invention, the precision of the total magnetization direction can be improved, and the pole conversion problem containing residual magnetism abnormity under the condition of multiple field sources is solved.

Description

Technical Field

[0001] The present application belongs to the field of geophysical magnetic exploration, and more specifically, to a method and device for digitizing complex magnetic anomalies containing residual magnetism. Background Art

[0002] Magnetic prospecting is an important exploration method in the fields of mineral exploration and geological survey. Underground magnetic bodies generally have residual magnetism, and its direction is generally different from the contemporary geomagnetic field. In the case of strong residual magnetism from multiple sources, the processing and interpretation of magnetic anomaly data will be affected. In actual work, magnetic anomaly data are usually polarized to eliminate the influence of oblique magnetization. However, if the magnetic anomaly is generated by multiple sources, the traditional polarization method cannot achieve good results, which in turn affects the subsequent processing and interpretation.

[0003] At present, in the treatment of complex magnetic anomalies containing 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 change the magnetization direction in the case of multiple magnetic sources. Summary of the invention

[0004] In view of the defects of the prior art, the purpose of the present application is to provide a method and device for polarization of complex magnetic anomaly data containing residual magnetism, aiming to solve the problem of large error in determining the total magnetization direction in the case of multiple magnetic sources, and the inability to change the magnetization direction in the case of multiple magnetic sources.

[0005] To achieve the above objectives, in a first aspect, the present application provides a method for digitizing complex magnetic anomalies containing residual magnetism, comprising: Obtain geophysical magnetic survey data of underground magnetic bodies; The geophysical magnetic survey data is converted into regularized magnetic source intensity by using a frequency domain processing method, and the horizontal position of the magnetic source center is determined according to the maximum value of the regularized magnetic source intensity; Obtaining the total magnetization direction of the entire domain according to the integration result; performing integration calculation according to the geophysical magnetic survey data to obtain the integration result of the magnetization characteristics, and obtaining the total magnetization direction of the entire domain according to the integration result; Obtaining the total magnetization direction of the magnetic body according to the horizontal position of the center of the magnetic source and the total magnetization direction of the entire domain; The geophysical magnetic survey data is subjected to polarization processing with variable magnetization direction according to the total magnetization direction, so as to obtain a polarization field of complex magnetic anomaly data containing residual magnetism.

[0006] Optionally, the method for determining the regularized magnetic source intensity includes: A magnetic gradient tensor matrix is ​​calculated based on the geophysical magnetic survey data; the magnetic gradient tensor matrix is ​​a symmetric matrix; Determine the eigenvalue matrix and the eigenvectors corresponding to the eigenvector of the symmetric matrix, and arrange the multiple eigenvalues ​​of the eigenvalue matrix in a monotonically decreasing order; The regularized magnetic source strength is determined according to a plurality of eigenvalues.

[0007] Optionally, the method for acquiring the magnetic gradient tensor matrix includes: Performing Fourier transformation on the geophysical magnetic survey data to determine the transverse wave number and the longitudinal wave number corresponding to the x-direction and the y-direction of the geophysical magnetic survey data, and determining the geomagnetic inclination and the geomagnetic declination of the geophysical magnetic survey data; Determine a first conversion factor in the frequency domain direction according to the transverse wave number, the longitudinal wave number, the geomagnetic inclination and the geomagnetic declination; Calculate the three components of magnetic anomaly in the x-direction, y-direction and z-direction respectively according to the geophysical magnetic survey data, the first conversion factor in the frequency domain direction, the transverse wave number and the longitudinal wave number; Partial derivatives are calculated based on the three components of magnetic anomalies in the x-direction, the y-direction and the z-direction, and a magnetic gradient tensor matrix is ​​constructed based on the partial derivative calculation results.

[0008] Optionally, the integration result includes an x-integral corresponding to the x direction, a y-integral corresponding to the y direction, and a z-integral corresponding to the z direction; The total magnetization direction of the entire domain is obtained according to the integration result, including: Perform an arcsine calculation based on the x integral, y integral, and z integral to obtain a total magnetization inclination; The judgment is made based on the x-integral and the y-integral, and the total magnetization deflection angle is calculated based on the judgment result combined with the x-integral and the y-integral.

[0009] Optionally, performing polarization processing of changing the magnetization direction on the geophysical magnetic survey data according to the total magnetization direction to obtain a polarization field of complex magnetic anomaly data containing remanent magnetization includes: Determine the inclination average value according to the total magnetization inclination, and determine the deflection average value according to the total magnetization deflection, and respectively determine the inclination difference and deflection difference; Determine a second conversion factor in the frequency domain direction according to the transverse wave number, the longitudinal wave number, the average inclination angle, and the average deflection angle; Performing two-dimensional Fourier transform and inverse transform according to the geophysical magnetic survey data, the transverse wave number, the longitudinal wave number, the first conversion factor and the second conversion factor to obtain a mean value of the magnetic field data; According to the mean value of the magnetic field data, partial derivatives are calculated for the average value of the inclination angle and the average value of the deflection angle respectively to obtain the partial derivative of the inclination angle and the partial derivative of the deflection angle; According to the mean value of the magnetic field data, the inclination difference, the inclination partial derivative, the deflection difference and the deflection partial derivative, the variable magnetization direction polarization processing is performed to obtain the polarization field of complex magnetic anomaly data containing remanent magnetization.

[0010] Optionally, calculating the total magnetization deflection angle 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, an inverse tangent calculation is performed based on the ratio of the x integral to the y integral to obtain a total magnetization deflection angle; If the x integral is less than 0 and the y integral is greater than or equal to 0, perform an inverse tangent calculation based on the ratio of the y integral to the x integral, add π to the result, and obtain the total magnetization deflection angle; If the x integral is less than 0 and the y integral is less than 0, perform an inverse tangent calculation based on the ratio of the y integral to the x integral, subtract π from the calculation result, and obtain the total magnetization deflection angle.

[0011] Optionally, obtaining the total magnetization direction of the magnetic body according to the horizontal position of the center of the magnetic source and the total magnetization direction of the entire domain includes: According to the horizontal position of the center of the magnetic source, the total magnetization inclination and the total magnetization deflection corresponding to the global total magnetization direction are selected, and the total magnetization direction is obtained according to the total magnetization inclination and the total magnetization deflection.

[0012] In a second aspect, the present application provides a complex magnetic anomaly data pole device containing residual magnetism, comprising: A data acquisition module, used to acquire geophysical magnetic survey data of underground magnetic bodies; A magnetic source center determination module, used to convert the geophysical magnetic survey data into 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; A global total magnetization direction module is used to obtain the global total magnetization direction according to the integration result, perform integral calculation according to the geophysical magnetic survey data, obtain the integral result of the magnetization characteristics, and obtain the global total magnetization direction according to the integration result; A total magnetization direction module, used to obtain the total magnetization direction of the magnetic body according to the horizontal position of the center of the magnetic source and the total magnetization direction of the whole domain; The polarization module is used to perform polarization processing on the geophysical magnetic survey data by changing the magnetization direction according to the total magnetization direction, so as to obtain the polarization field of complex magnetic anomaly data containing remanent magnetization.

[0013] In a third aspect, the present application provides an electronic device comprising: at least one memory for storing programs; and at least one processor for executing the programs 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 possible implementation of the first aspect.

[0014] In a fourth aspect, the present application provides a computer-readable storage medium, which stores a computer program. When the computer program runs on a processor, the processor executes the method described in the first aspect or any possible implementation of the first aspect.

[0015] In a fifth aspect, the present application provides a computer program product. When the computer program product runs on a processor, the processor executes the method described in the first aspect or any possible implementation of the first aspect.

[0016] It can be understood that the beneficial effects of the second to fifth aspects mentioned above can be found in the relevant description of the first aspect mentioned above, and will not be repeated here.

[0017] In general, the above technical solutions conceived by this application have the following beneficial effects compared with the prior art: (1) By considering the total magnetization direction of multiple field sources, the present application can more accurately reflect the magnetic field characteristics under complex geological environments. 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 the polarization processing, the accuracy of the total magnetization direction calculated by the Helbig integral is improved, while the error caused by residual magnetism is reduced. Therefore, the present application can meet the needs under complex geological conditions and solve the polarization problem of residual magnetism anomalies under multiple magnetic sources.

[0018] (2) This application uses the Helbig integral to estimate the total magnetization direction of multiple field sources, and uses the variable magnetization direction polarization technology to obtain the polarization field of non-isolated field sources, providing data support for the processing and interpretation of magnetic anomaly data under residual magnetization conditions. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] Figure 1 This is one of the flow charts of the method for digitizing complex magnetic anomalies containing residual magnetism provided in the embodiment of the present application; Figure 2 A schematic diagram of an underground magnetic body model according to an embodiment of the present application; Figure 3 The geophysical magnetic survey data of the embodiment of the present application and the polarization field under theoretical conditions; Figure 4 It is the normalized source strength NSS (Normalized Source Strength) and the horizontal position of the magnetic source center calculated according to the geophysical magnetic survey data in the embodiment of the present application; Figure 5 The total magnetization inclination and total magnetization deflection calculated based on the geophysical magnetic survey data in the embodiment of the present application; Figure 6The geomagnetic direction-oriented pole field of the experimental model in the embodiment of the present application and the variable inclination pole field considering the residual magnetism; Figure 7 This is the second flow chart of the method for digitizing complex magnetic anomalies containing residual magnetism provided in the embodiment of the present application; Figure 8 It is a structural schematic diagram of a complex magnetic anomaly data pole device containing residual magnetism provided in an embodiment of the present application; Fig. 9 It is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0020] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0021] The term "and / or" in this article is a description of the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. The symbol " / " in this article indicates that the associated objects are in an or relationship, for example, A / B means A or B.

[0022] The terms "first" and "second" in the specification and claims herein are used to distinguish different objects rather than to describe a specific order of the objects. For example, a first response message and a second response message are used to distinguish different response messages rather than to describe a specific order of the response messages.

[0023] In the embodiments of the present application, words such as "exemplary" or "for example" are used to indicate examples, illustrations or descriptions. Any embodiment or design described as "exemplary" or "for example" in the embodiments of the present application should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of words such as "exemplary" or "for example" is intended to present related concepts in a specific way.

[0024] In the description of the embodiments of the present application, unless otherwise specified, "multiple" means two or more than two. For example, multiple processing units refer to two or more processing units, etc.; multiple elements refer to two or more elements, etc.

[0025] The embodiments of the present application are described below in conjunction with the drawings in the embodiments of the present application.

[0026] Reference Figure 1 The present application provides a method for digitizing complex magnetic anomalies containing residual magnetism, including: S101. Obtain geophysical magnetic data of underground magnetic bodies; S102. Converting the geophysical magnetic survey data into regularized magnetic source intensity 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; S103. Obtaining the total magnetization direction of the entire domain according to the integral result. Performing integral calculation according to the geophysical magnetic survey data to obtain the integral result of the magnetization characteristics, and obtaining the total magnetization direction of the entire domain according to the integral result; S104. Obtaining the total magnetization direction of the magnetic body according to the horizontal position of the center of the magnetic source and the total magnetization direction of the entire domain; S105. Performing polarization processing of changing the magnetization direction on the geophysical magnetic survey data according to the total magnetization direction, so as to obtain a polarization field of complex magnetic anomaly data containing remanent magnetization.

[0027] Specifically, in this embodiment S101, it is first necessary to obtain relevant magnetic survey data of underground magnetic bodies through geophysical detection technology.

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

[0029] It should be noted that Easting in this figure and subsequent figures means easting, Northing means northing, and Depth means depth.

[0030] The x direction in the embodiments of the present application refers to the east-west direction (east-west direction). The y direction generally refers to the north-south direction (south-north direction). The z direction refers to the depth direction, which refers to the direction perpendicular to the ground.

[0031] Reference Figure 3 , Figure 3 (a) is the magnetic anomaly generated by the model, which is equivalent to the geophysical magnetic survey data in actual work. Figure 3 (b) in the figure represents the polarization field in the theoretical case, which is the expected output.

[0032] Then, in step S102, the geophysical magnetic survey data is converted into regularized magnetic source intensity using a frequency domain processing method. NSS .

[0033] In this step, the acquired geophysical magnetic survey data is converted by frequency domain processing to obtain the regularized magnetic source intensity. This process involves Fourier transforming the magnetic survey data to convert 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 results. Finally, the center position of the magnetic source is further determined based on the regularized magnetic source intensity. In this embodiment, NSS The position of the maximum value is taken as the center of the magnetic source.

[0034] Furthermore, an integral calculation, namely Helbig integral, is performed according to the geophysical magnetic survey data in step S103. In this step, the geophysical magnetic survey data is analyzed using Helbig integral to obtain the global total magnetization direction.

[0035] It should be noted that Helbig integration is a mathematical method widely used in magnetic field analysis. By integrating the acquired magnetic measurement data, it is possible to identify magnetic field effects from various sources.

[0036] Through S104, the total magnetization direction of the magnetic body is obtained according to the horizontal position of the center of the magnetic source and the total magnetization direction of the whole domain. According to the horizontal position of the center of the magnetic source determined in the previous step and the total magnetization direction of the whole domain obtained by the Helbig integral. According to the position of the magnetic source and its magnetization characteristics, a comprehensive description of the magnetization direction is formed. The total magnetization direction model can reflect the overall magnetization state.

[0037] Finally, in S105, the geophysical magnetic survey data is subjected to polarization processing of changing the magnetization direction by using the total magnetization direction obtained previously. This processing process mainly involves adjusting the measured magnetic field data to a form consistent with the total magnetization direction, thereby reducing interference caused by residual magnetic anomalies.

[0038] Optionally, the method for determining the regularized magnetic source intensity includes: A magnetic gradient tensor matrix is ​​calculated based on the geophysical magnetic survey data; the magnetic gradient tensor matrix is ​​a symmetric matrix; Determine the eigenvalue matrix and the eigenvectors corresponding to the eigenvector of the symmetric matrix, and arrange the multiple eigenvalues ​​of the eigenvalue matrix in a monotonically decreasing order; The regularized magnetic source strength is determined according to a plurality of eigenvalues.

[0039] Furthermore, the method for obtaining the magnetic gradient tensor matrix includes: Performing Fourier transformation on the geophysical magnetic survey data to determine the transverse wave number and the longitudinal wave number corresponding to the x-direction and the y-direction of the geophysical magnetic survey data, and determining the geomagnetic inclination and the geomagnetic declination of the geophysical magnetic survey data; Determine a first conversion factor in the frequency domain direction according to the transverse wave number, the longitudinal wave number, the geomagnetic inclination and the geomagnetic declination; Calculate the three components of magnetic anomaly in the x-direction, y-direction and z-direction respectively according to the geophysical magnetic survey data, the first conversion factor in the frequency domain direction, the transverse wave number and the longitudinal wave number; Partial derivatives are calculated based on the three components of magnetic anomalies in the x-direction, the y-direction and the z-direction, and a magnetic gradient tensor matrix is ​​constructed based on the partial derivative calculation results.

[0040] Specifically, this embodiment calculates the regularized magnetic source intensity in the frequency domain NSS The specific process is as follows: According to the geophysical magnetic survey data, the regularized magnetic source intensity is calculated in the frequency domain NSS , as the first set of data, such as Figure 4 As shown, the data size is a real matrix of 41×41, and the calculation process is as follows: Geophysical magnetic survey data Calculate the magnetic gradient tensor

[0041] According to the potential field theory, we have , so this gradient tensor is a symmetric matrix containing only six independent components, so It can be expressed as:

[0042]

[0043] in, Symmetrical array The characteristic value of is the eigenvalue The corresponding eigenvectors arrange the eigenvalues ​​of the magnetic gradient tensor in monotonically decreasing order: , then the regularized magnetic source intensity NSS It is expressed as:

[0044] Last Selection NSS The position of the maximum value is taken as the horizontal position of the center of the magnetic source, such as Figure 3 The position of the "white cross" in the picture.

[0045] Optionally, the integration result includes an x-integral corresponding to the x direction, a y-integral corresponding to the y direction, and a z-integral corresponding to the z direction; The total magnetization direction of the entire domain is obtained according to the integration result, including: Perform an arcsine calculation based on the x integral, y integral, and z integral to obtain a total magnetization inclination; The judgment is made based on the x-integral and the y-integral, and the total magnetization deflection angle is calculated based on the judgment result combined with the x-integral and the y-integral.

[0046] Further, the total magnetization deflection angle is calculated according to the judgment result in combination with the x integral and the y integral, including: If the x integral is greater than or equal to 0, an inverse tangent calculation is performed based on the ratio of the x integral to the y integral to obtain a total magnetization deflection angle; If the x integral is less than 0 and the y integral is greater than or equal to 0, perform an inverse tangent calculation based on the ratio of the y integral to the x integral, add π to the result, and obtain the total magnetization deflection angle; If the x integral is less than 0 and the y integral is less than 0, perform an inverse tangent calculation based on the ratio of the y integral to the x integral, subtract π from the calculation result, and obtain the total magnetization deflection angle.

[0047] Specifically, this embodiment uses geophysical magnetic survey data To calculate the global magnetization, the steps are as follows: First, the calculation formula for Helbig integral is as follows:

[0048] in, The geophysical magnetic data are The three components of the calculated magnetic anomaly are then used to obtain the total magnetization direction of the entire domain using the following formula: , , in, and Represent the total magnetization inclination and total magnetization deflection in the total magnetization direction of the whole domain, such as Figure 5 As shown, Figure 5 (a) and Figure 5 In (b), the horizontal coordinates are all east and the vertical coordinates are all north.

[0049] Optionally, obtaining the total magnetization direction of the magnetic body according to the horizontal position of the center of the magnetic source and the total magnetization direction of the entire domain includes: According to the horizontal position of the center of the magnetic source, the total magnetization inclination and the total magnetization deflection corresponding to the global total magnetization direction are selected, and the total magnetization direction is obtained according to the total magnetization inclination and the total magnetization deflection.

[0050] According to the center position of the magnetic source obtained in the above embodiment, the total magnetization inclination and total magnetization deflection corresponding to the total magnetization direction of the entire domain are selected as the total magnetization direction of the magnetic source. The estimated error is shown in Table 1.

[0051] Table 1 Total magnetization direction estimation error

[0052] Optionally, performing polarization processing of changing the magnetization direction on the geophysical magnetic survey data according to the total magnetization direction to obtain a polarization field of complex magnetic anomaly data containing remanent magnetization includes: Determine the inclination average value according to the total magnetization inclination, and determine the deflection average value according to the total magnetization deflection, and respectively determine the inclination difference and deflection difference; Determine a second conversion factor in the frequency domain direction according to the transverse wave number, the longitudinal wave number, the average inclination angle, and the average deflection angle; Performing two-dimensional Fourier transform and inverse transform according to the geophysical magnetic survey data, the transverse wave number, the longitudinal wave number, the first conversion factor and the second conversion factor to obtain a mean value of the magnetic field data; According to the mean value of the magnetic field data, partial derivatives are calculated for the average value of the inclination angle and the average value of the deflection angle respectively to obtain the partial derivative of the inclination angle and the partial derivative of the deflection angle; According to the mean value of the magnetic field data, the inclination difference, the inclination partial derivative, the deflection difference and the deflection partial derivative, the variable magnetization direction polarization processing is performed to obtain the polarization field of complex magnetic anomaly data containing remanent magnetization.

[0053] Specifically, according to the total magnetization direction of the magnetic source, the geophysical magnetic survey data Carry out polarization treatment by changing the magnetization direction to obtain the polarization field , the calculation formula is as follows:

[0054] in,

[0055]

[0056]

[0057] represents geophysical magnetic survey data, denote the two-dimensional Fourier transform and its inverse transform, respectively. is an imaginary unit, After Fourier transformation, and The wave number of the direction, , is the frequency domain direction conversion factor, namely the first conversion factor and the second conversion factor, are the geomagnetic inclination and declination, respectively, which are Subscript The direction indicated, are the average values ​​of the total magnetization inclination and total magnetization deflection respectively, that is, Subscript The direction indicated, and Represent the difference between the total magnetization inclination and the total magnetization deflection and their average values, respectively.

[0058] Reference Figure 6 , the polarization field obtained when the remanence is not considered is as follows Figure 6 As shown in (a), the high value of magnetic anomaly does not correspond to the center of magnetic source. Considering the variable magnetization direction of remanent magnetization, the polar field is as follows: Figure 6 As shown in (b), the oblique magnetization is effectively eliminated, and the center of the magnetic source corresponds to the high-value area of ​​the magnetic anomaly.

[0059] It should be noted that the above calculation formula requires the geomagnetic direction and the total magnetization direction as input. In actual work, the geomagnetic direction can be obtained by querying the International Geomagnetic Reference Field (IGRF). The total magnetization direction input is to select a total magnetization direction within the set total magnetization direction range. Reference Figure 7 , Figure 7 The present application embodiment provides a complete flow chart of a method for digitizing a complex magnetic anomaly containing residual magnetism, including the following steps: Input geophysical magnetic survey data; Calculate the Helbig integral to obtain the total magnetization direction of the entire domain; Calculate the regularized magnetic source strength NSS Get the horizontal position of the center of the magnetic source; Construct a model of the total magnetization direction; To change the magnetization direction and polarity of geophysical magnetic survey data; Input geophysical magnetic survey data.

[0060] Reference Figure 8 The present application provides a complex magnetic anomaly data pole device containing residual magnetism, comprising: The data acquisition module 810 is used to acquire geophysical magnetic survey data of underground magnetic bodies; A magnetic source center determination module 820 is used to convert the geophysical magnetic survey data into regularized magnetic source intensity 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; A global total magnetization direction module 830 is used to obtain the global total magnetization direction according to the integration result, perform integration calculation according to the geophysical magnetic survey data, obtain the integration result of the magnetization characteristics, and obtain the global total magnetization direction according to the integration result; A total magnetization direction module 840 is used to obtain the total magnetization direction of the magnetic body according to the horizontal position of the center of the magnetic source and the total magnetization direction of the entire domain; The polarization module 850 is used to perform polarization processing on the geophysical magnetic survey data by changing the magnetization direction according to the total magnetization direction, so as to obtain a polarization field of complex magnetic anomaly data containing remanent magnetization.

[0061] Optionally, the method for determining the regularized magnetic source intensity includes: A magnetic gradient tensor matrix is ​​calculated based on the geophysical magnetic survey data; the magnetic gradient tensor matrix is ​​a symmetric matrix; Determine the eigenvalue matrix and the eigenvectors corresponding to the eigenvector of the symmetric matrix, and arrange the multiple eigenvalues ​​of the eigenvalue matrix in a monotonically decreasing order; The regularized magnetic source strength is determined according to a plurality of eigenvalues.

[0062] Optionally, the method for acquiring the magnetic gradient tensor matrix includes: Performing Fourier transformation on the geophysical magnetic survey data to determine the transverse wave number and the longitudinal wave number corresponding to the x-direction and the y-direction of the geophysical magnetic survey data, and determining the geomagnetic inclination and the geomagnetic declination of the geophysical magnetic survey data; Determine a first conversion factor in the frequency domain direction according to the transverse wave number, the longitudinal wave number, the geomagnetic inclination and the geomagnetic declination; Calculate the three components of magnetic anomaly in the x-direction, y-direction and z-direction respectively according to the geophysical magnetic survey data, the first conversion factor in the frequency domain direction, the transverse wave number and the longitudinal wave number; Partial derivatives are calculated based on the three components of magnetic anomalies in the x-direction, the y-direction and the z-direction, and a magnetic gradient tensor matrix is ​​constructed based on the partial derivative calculation results.

[0063] Optionally, the integration result includes an x-integral corresponding to the x direction, a y-integral corresponding to the y direction, and a z-integral corresponding to the z direction; The total magnetization direction of the entire domain is obtained according to the integration result, including: Perform an arcsine calculation based on the x integral, y integral, and z integral to obtain a total magnetization inclination; The judgment is made based on the x-integral and the y-integral, and the total magnetization deflection angle is calculated based on the judgment result combined with the x-integral and the y-integral.

[0064] Optionally, performing polarization processing of changing the magnetization direction on the geophysical magnetic survey data according to the total magnetization direction to obtain a polarization field of complex magnetic anomaly data containing remanent magnetization includes: Determine the inclination average value according to the total magnetization inclination, and determine the deflection average value according to the total magnetization deflection, and respectively determine the inclination difference and deflection difference; Determine a second conversion factor in the frequency domain direction according to the transverse wave number, the longitudinal wave number, the average inclination angle, and the average deflection angle; Performing two-dimensional Fourier transform and inverse transform according to the geophysical magnetic survey data, the transverse wave number, the longitudinal wave number, the first conversion factor and the second conversion factor to obtain a mean value of the magnetic field data; According to the mean value of the magnetic field data, partial derivatives are calculated for the average value of the inclination angle and the average value of the deflection angle respectively to obtain the partial derivative of the inclination angle and the partial derivative of the deflection angle; According to the mean value of the magnetic field data, the inclination difference, the inclination partial derivative, the deflection difference and the deflection partial derivative, the variable magnetization direction polarization processing is performed to obtain the polarization field of complex magnetic anomaly data containing remanent magnetization.

[0065] Optionally, calculating the total magnetization deflection angle 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, an inverse tangent calculation is performed based on the ratio of the x integral to the y integral to obtain a total magnetization deflection angle; If the x integral is less than 0 and the y integral is greater than or equal to 0, perform an inverse tangent calculation based on the ratio of the y integral to the x integral, add π to the result, and obtain the total magnetization deflection angle; If the x integral is less than 0 and the y integral is less than 0, perform an inverse tangent calculation based on the ratio of the y integral to the x integral, subtract π from the calculation result, and obtain the total magnetization deflection angle.

[0066] Optionally, obtaining the total magnetization direction of the magnetic body according to the horizontal position of the center of the magnetic source and the total magnetization direction of the entire domain includes: According to the horizontal position of the center of the magnetic source, the total magnetization inclination and the total magnetization deflection corresponding to the global total magnetization direction are selected, and the total magnetization direction is obtained according to the total magnetization inclination and the total magnetization deflection.

[0067] It can be understood that the detailed functional implementation of each of the above-mentioned units / modules can be found in the introduction of the aforementioned method embodiment, and will not be repeated here.

[0068] It should be understood that the above-mentioned device is used to execute the method in the above-mentioned embodiment. The implementation principle and technical effect of the corresponding program module in the device are similar to those described in the above-mentioned method. The working process of the device can refer to the corresponding process in the above-mentioned method, which will not be repeated here.

[0069] Reference Fig. 9Based on the method in the above embodiment, the 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, wherein the processor 910, the communication interface 920, and the memory 930 communicate with each other through the communication bus 940. The processor 910 can call the logic instructions in the memory 930 to execute the method in the above embodiment.

[0070] In addition, the logic instructions in the above-mentioned memory 930 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product, which is stored in a storage medium and includes a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present application.

[0071] Based on the method in the above embodiment, an embodiment of the present application provides a computer-readable storage medium, which stores a computer program. When the computer program runs on a processor, the processor executes the method in the above embodiment.

[0072] Based on the method in the above embodiment, an embodiment of the present application provides a computer program product. When the computer program product runs on a processor, the processor executes the method in the above embodiment.

[0073] It is understandable that the processor in the embodiment of the present application may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field programmable gate arrays (FPGA) 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.

[0074] The method steps in the embodiments of the present application can be implemented by hardware or by a processor executing software instructions. The software instructions can be composed of corresponding software modules, and the software modules can be stored in random access memory (RAM), flash memory, read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), registers, hard disks, mobile hard disks, CD-ROMs, or any other form of storage medium 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 can also be a component of the processor. The processor and the storage medium can be located in an ASIC.

[0075] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware or any combination thereof. When implemented by 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 process or function described in the embodiment of the present application is generated in whole or in part. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions may be stored in a computer-readable storage medium or transmitted through the computer-readable storage medium. The computer instructions may be transmitted from a website site, computer, server or data center to another website site, computer, server or data center by wired (e.g., coaxial cable, optical fiber, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium may be any available medium that a computer can access or a data storage device such as a server or data center that includes one or more available media integrated. The available medium may be a magnetic medium (e.g., a floppy disk, a hard disk, a tape), an optical medium (e.g., a DVD), or a semiconductor medium (e.g., a solid-state drive (SSD)), etc.

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

[0077] It will be easily understood by those skilled in the art that the above description is only a preferred embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present application shall be included in the scope of protection of the present application.

Claims

1. A method for digitizing complex magnetic anomalies containing residual magnetism, characterized in that: include: Obtain geophysical magnetic survey data of underground magnetic bodies; The geophysical magnetic survey data is converted into regularized magnetic source intensity by using a frequency domain processing method, and the horizontal position of the magnetic source center is determined according to the maximum value of the regularized magnetic source intensity; Performing an integral calculation based on the geophysical magnetic survey data to obtain an integral result of the magnetization characteristics, and obtaining a global total magnetization direction based on the integral result; Obtaining the total magnetization direction of the magnetic body according to the horizontal position of the center of the magnetic source and the total magnetization direction of the entire domain; The geophysical magnetic survey data is subjected to polarization processing with variable magnetization direction according to the total magnetization direction, so as to obtain a polarization field of complex magnetic anomaly data containing residual magnetism.

2. The method for digitizing complex magnetic anomalies containing residual magnetism according to claim 1, characterized in that: The method for determining the regularized magnetic source intensity includes: A magnetic gradient tensor matrix is ​​calculated based on the geophysical magnetic survey data; the magnetic gradient tensor matrix is ​​a symmetric matrix; Determine the eigenvalue matrix and the eigenvectors corresponding to the eigenvector of the symmetric matrix, and arrange the multiple eigenvalues ​​of the eigenvalue matrix in a monotonically decreasing order; The regularized magnetic source strength is determined according to a plurality of eigenvalues.

3. The method for digitizing complex magnetic anomalies containing residual magnetism according to claim 2, characterized in that: The method for obtaining the magnetic gradient tensor matrix includes: Performing Fourier transformation on the geophysical magnetic survey data to determine the transverse wave number and the longitudinal wave number corresponding to the x-direction and the y-direction of the geophysical magnetic survey data, and determining the geomagnetic inclination and the geomagnetic declination of the geophysical magnetic survey data; Determine a first conversion factor in the frequency domain direction according to the transverse wave number, the longitudinal wave number, the geomagnetic inclination and the geomagnetic declination; Calculate the three components of magnetic anomaly in the x-direction, y-direction and z-direction respectively according to the geophysical magnetic survey data, the first conversion factor in the frequency domain direction, the transverse wave number and the longitudinal wave number; Partial derivatives are calculated based on the three components of magnetic anomalies in the x-direction, the y-direction and the z-direction, and a magnetic gradient tensor matrix is ​​constructed based on the partial derivative calculation results.

4. The method for digitizing complex magnetic anomalies containing residual magnetism according to claim 1, characterized in that: The integration result includes an x-integral corresponding to the x direction, a y-integral corresponding to the y direction, and a z-integral corresponding to the z direction; The total magnetization direction of the entire domain is obtained according to the integration result, including: Perform an arcsine calculation based on the x integral, y integral, and z integral to obtain a total magnetization inclination; The judgment is made based on the x-integral and the y-integral, and the total magnetization deflection angle is calculated based on the judgment result combined with the x-integral and the y-integral.

5. The method for digitizing complex magnetic anomalies containing residual magnetism according to claim 4, characterized in that: The method of performing polarization processing of changing the magnetization direction on the geophysical magnetic survey data according to the total magnetization direction to obtain a polarization field of complex magnetic anomaly data containing remanent magnetization includes: Determine the inclination average value according to the total magnetization inclination, and determine the deflection average value according to the total magnetization deflection, and respectively determine the inclination difference and deflection difference; Determine a second conversion factor in the frequency domain direction according to the transverse wave number, the longitudinal wave number, the average inclination angle, and the average deflection angle; Performing two-dimensional Fourier transform and inverse transform according to the geophysical magnetic survey data, the transverse wave number, the longitudinal wave number, the first conversion factor and the second conversion factor to obtain a mean value of the magnetic field data; According to the mean value of the magnetic field data, partial derivatives are calculated for the average value of the inclination angle and the average value of the deflection angle respectively to obtain the partial derivative of the inclination angle and the partial derivative of the deflection angle; According to the mean value of the magnetic field data, the inclination difference, the inclination partial derivative, the deflection difference and the deflection partial derivative, the variable magnetization direction polarization processing is performed to obtain the polarization field of complex magnetic anomaly data containing residual magnetism.

6. The method for digitizing complex magnetic anomalies containing residual magnetism according to claim 4, characterized in that: The method of calculating the total magnetization deflection angle according to the judgment result combined with the x integral and the y integral includes: If the x integral is greater than or equal to 0, an inverse tangent calculation is performed based on the ratio of the x integral to the y integral to obtain a total magnetization deflection angle; If the x integral is less than 0 and the y integral is greater than or equal to 0, perform an inverse tangent calculation based on the ratio of the y integral to the x integral, add π to the result, and obtain the total magnetization deflection angle; If the x integral is less than 0 and the y integral is less than 0, perform an inverse tangent calculation based on the ratio of the y integral to the x integral, subtract π from the calculation result, and obtain the total magnetization deflection angle.

7. The method for digitizing complex magnetic anomalies containing residual magnetism according to claim 1, characterized in that: The method of obtaining the total magnetization direction of the magnetic body according to the horizontal position of the center of the magnetic source and the total magnetization direction of the entire domain includes: According to the horizontal position of the center of the magnetic source, the total magnetization inclination and the total magnetization deflection corresponding to the global total magnetization direction are selected, and the total magnetization direction is obtained according to the total magnetization inclination and the total magnetization deflection.

8. A complex magnetic anomaly data polarization device containing residual magnetism, characterized in that: include: A data acquisition module, used to acquire geophysical magnetic survey data of underground magnetic bodies; A magnetic source center determination module, used to convert the geophysical magnetic survey data into 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; A global total magnetization direction module is used to obtain the global total magnetization direction according to the integration result, perform integral calculation according to the geophysical magnetic survey data, obtain the integral result of the magnetization characteristics, and obtain the global total magnetization direction according to the integration result; A total magnetization direction module, used to obtain the total magnetization direction of the magnetic body according to the horizontal position of the center of the magnetic source and the total magnetization direction of the entire domain; The polarization module is used to perform polarization processing on the geophysical magnetic survey data by changing the magnetization direction according to the total magnetization direction, so as to obtain the polarization field of complex magnetic anomaly data containing remanent magnetization.

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

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

Images