A magnetic moment calculation method and device based on scalar magnetic anomaly and related equipment

Abstract

The application discloses a magnetic moment calculation method and device based on scalar magnetic anomaly and related equipment, and the method comprises the following steps: acquiring two-dimensional scalar magnetic anomaly data, and determining the three-dimensional coordinates of a target body by using the two-dimensional scalar magnetic anomaly data, wherein the record points of the two-dimensional scalar magnetic anomaly data are uniformly distributed in a detection area; combining the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data to obtain the modulus of the two-dimensional scalar magnetic anomaly on each group of unit orthogonal bases, wherein the unit orthogonal bases are obtained by decomposing the two-dimensional scalar magnetic anomaly data generated by the magnetic dipole of the target body at an arbitrary height; substituting the modulus of the two-dimensional scalar magnetic anomaly on each group of unit orthogonal bases into a preset overdetermined equation set, obtaining the minimum variance solution of the overdetermined equation set, and obtaining the three components of the vector magnetic moment of the target body; then, the size, inclination and declination information of the magnetic moment can be obtained by appropriately transforming the three components of the vector magnetic moment.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of magnetic exploration, and more particularly to a magnetic moment calculation method and device based on scalar magnetic anomaly and related equipment. BACKGROUND

[0002] Magnetic anomaly detection (MAD) is one of the most suitable geophysical techniques for locating and mapping the distribution of ferromagnetic metal objects, and has a wide range of applications in the fields of resource exploration, unexploded bomb detection, underwater submarine detection, archaeology, and other fields related to national economic security and cultural construction. Due to the low noise level of scalar magnetic sensors, and the fact that they are not sensitive to mechanical noise such as rotation and vibration, scalar magnetic anomaly detection is currently the mainstream method of magnetic anomaly detection, and the calculation of ferromagnetic target parameters based on two-dimensional scalar magnetic anomaly data is also an important research direction in magnetic anomaly detection.

[0003] In some existing methods for determining the magnetic target parameters based on single magnetic sensor data, initial parameters such as the size of the target body magnetic moment, the inclination angle, the declination angle, and the position are assumed, and then a nonlinear optimization theory is applied to iteratively solve the parameters to fit the observed data. In the iteration process, the various parameters to be solved affect each other, which is not conducive to solving the magnetic moment information of the target. SUMMARY

[0004] Therefore, the present application provides a magnetic moment calculation method and device based on scalar magnetic anomaly and related equipment to realize the calculation of the magnetic moment of the target body.

[0005] To achieve the above object, the first aspect of the present application provides a magnetic moment calculation method based on scalar magnetic anomaly, comprising:

[0006] Obtaining two-dimensional scalar magnetic anomaly data, and determining the three-dimensional coordinates of the target body using the two-dimensional scalar magnetic anomaly data, wherein the recording points of the two-dimensional scalar magnetic anomaly data are uniformly distributed in the detection area;

[0007] Combining the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data, obtaining the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases, wherein the unit orthogonal bases are obtained by decomposing the two-dimensional scalar magnetic anomaly data generated by the magnetic dipole of the target body at any height;

[0008] Substituting the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases into a predetermined overdetermined equation set, and obtaining the minimum variance solution of the overdetermined equation set to obtain the three components of the vector magnetic moment of the target body.

[0009] Preferably, the process of obtaining the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases in combination with the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data comprises:

[0010] Substituting the three-dimensional coordinates of the target body into the expression of each set of unit orthogonal bases of the two-dimensional scalar magnetic anomaly, the two-dimensional orthogonal base expression of the target body is obtained.

[0011] Based on the preset size window and the two-dimensional orthogonal base expression, the inner product calculation is performed on the two-dimensional scalar magnetic anomaly data at each grid point to obtain the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases.

[0012] Preferably, the two-dimensional orthogonal base expression comprises:

[0013]

[0014] Each unit orthogonal base satisfies the following equation:

[0015]

[0016] Wherein, (x0, y0, z0) is the three-dimensional coordinates of the target body, and (x, y, z) is the three-dimensional coordinates of the record point of the two-dimensional scalar magnetic anomaly data.

[0017] Preferably, the process of obtaining the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases based on the preset size window and the two-dimensional orthogonal base expression, and performing the inner product calculation on the two-dimensional scalar magnetic anomaly data at each grid point comprises:

[0018] The modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases is calculated by the following equation: n (x0, y0):

[0019]

[0020] Wherein, I and J are the size of the window in x direction and y direction respectively, and ΔT(x i ,y j ) is the two-dimensional scalar magnetic anomaly data of the grid point at (x i ,y j ).

[0021] Preferably, the process of obtaining the three components of the vector magnetic moment of the target body by substituting the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases into the preset overdetermined equation set and obtaining the minimum variance solution of the overdetermined equation set comprises:

[0022] The three components (m x ,m y ,mz

[0023]

[0024] wherein, α i is the modulus of the two-dimensional scalar magnetic anomaly on the i-th set of unit orthogonal bases, dz is the distance from the target body to the observation plane, I and A are the geomagnetic inclination and geomagnetic declination of the place where the detection matrix is located, respectively.

[0025] Preferably, the process of obtaining the two-dimensional scalar magnetic anomaly data comprises:

[0026] Using the magnetometer, the scalar total magnetic field data of the magnetic total field sensor in the detection area is recorded along the grid points of the preset grid, and a plurality of data items are obtained, each data item comprising the scalar total magnetic field data and the coordinates of the grid point corresponding to the scalar total magnetic field data, the grid being uniformly distributed in the detection area.

[0027] The scalar total magnetic field data in each data item is subjected to a daily variation correction process to obtain the scalar magnetic anomaly data of each data item.

[0028] The two-dimensional scalar magnetic anomaly data is composed of the scalar magnetic anomaly data at each grid point and the coordinates of each grid point.

[0029] Preferably, the process of determining the three-dimensional coordinates of the target body using the two-dimensional scalar magnetic anomaly data comprises:

[0030] The three-dimensional coordinates of the target body are determined using the two-dimensional scalar magnetic anomaly data by the Euler deconvolution method or a preset target three-dimensional positioning method based on scalar magnetic anomaly.

[0031] The second aspect of the present application provides a scalar magnetic anomaly-based magnetic moment calculation device, comprising:

[0032] A coordinate calculation unit is configured to obtain two-dimensional scalar magnetic anomaly data and determine the three-dimensional coordinates of the target body using the two-dimensional scalar magnetic anomaly data, the recording points of the two-dimensional scalar magnetic anomaly data being uniformly distributed in the detection area.

[0033] A modulus calculation unit is configured to obtain the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases in combination with the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data, the unit orthogonal bases being obtained by decomposing the two-dimensional scalar magnetic anomaly data generated by the magnetic dipole of the target body at any height.

[0034] A magnetic moment calculation unit is configured to substitute the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases into a preset overdetermined equation set and obtain the minimum variance solution of the overdetermined equation set to obtain the three-component vector magnetic moment of the target body. ​

[0035] The third aspect of the present application provides a magnetic moment calculation device based on scalar magnetic anomaly, comprising a memory and a processor;

[0036] The memory is configured to store a program;

[0037] The processor is configured to execute the program to implement each step of the magnetic moment calculation method based on scalar magnetic anomaly.

[0038] The fourth aspect of the present application provides a storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements each step of the magnetic moment calculation method based on scalar magnetic anomaly.

[0039] According to the above technical solution, the two-dimensional scalar magnetic anomaly data is first obtained, and the three-dimensional coordinates of the target body are determined by using the two-dimensional scalar magnetic anomaly data, wherein the recording points of the two-dimensional scalar magnetic anomaly data are uniformly distributed in the detection area. Then, the modulus of the two-dimensional scalar magnetic anomaly on each group of unit orthogonal bases is obtained in combination with the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data, wherein the unit orthogonal bases are obtained by decomposing the two-dimensional scalar magnetic anomaly data generated by the magnetic dipole of the target body at any height. Finally, the modulus of the two-dimensional scalar magnetic anomaly on each group of unit orthogonal bases is substituted into a predetermined overdetermined equation set, and the minimum variance solution of the overdetermined equation set is obtained to obtain the three components of the vector magnetic moment of the target body. The size, inclination and declination information of the magnetic moment can be obtained by appropriately transforming the three components of the vector magnetic moment, and the magnetic moment calculation of the target body is realized. BRIEF DESCRIPTION OF DRAWINGS

[0040] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed in the embodiments or prior art description will be briefly introduced. Obviously, the drawings in the following description are only embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor on the basis of the provided drawings.

[0041] Figure 1 The schematic diagram of the magnetic moment calculation method based on scalar magnetic anomaly disclosed in the embodiments of the present application is shown;

[0042] Figure 2 The magnetic anomaly detection two-dimensional survey line distribution diagram disclosed in the embodiments of the present application is shown;

[0043] Figure 3 The two-dimensional magnetic total field gradient surface diagram provided by the embodiments of the present application is shown;

[0044] Figure 4The figure of the two-dimensional unit orthogonal base surface calculated by the first group of unit orthogonal base functions provided in the embodiment of the present application is shown.

[0045] Figure 5 The figure of the two-dimensional unit orthogonal base surface calculated by the second group of unit orthogonal base functions provided in the embodiment of the present application is shown.

[0046] Figure 6 The figure of the two-dimensional unit orthogonal base surface calculated by the third group of unit orthogonal base functions provided in the embodiment of the present application is shown.

[0047] Figure 7 The figure of the two-dimensional unit orthogonal base surface calculated by the fourth group of unit orthogonal base functions provided in the embodiment of the present application is shown.

[0048] Figure 8 The figure of the two-dimensional unit orthogonal base surface calculated by the fifth group of unit orthogonal base functions provided in the embodiment of the present application is shown.

[0049] Figure 9 The schematic diagram of the magnetic moment calculation device based on scalar magnetic anomaly disclosed in the embodiment of the present application is shown.

[0050] Figure 10 The schematic diagram of the magnetic moment calculation device based on scalar magnetic anomaly disclosed in the embodiment of the present application is shown. DETAILED DESCRIPTION

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

[0052] The magnetic moment calculation method based on scalar magnetic anomaly provided in the embodiment of the present application will be introduced below. Please refer to Figure 1 The magnetic moment calculation method based on scalar magnetic anomaly provided in the embodiment of the present application can include the following steps:

[0053] In step S101, the two-dimensional scalar magnetic anomaly data is acquired, and the three-dimensional coordinates of the target body are determined by using the two-dimensional scalar magnetic anomaly data.

[0054] The two-dimensional total magnetic field anomaly data is obtained by applying the data processing method such as the daily variation correction to the total magnetic field data of the total magnetic field sensor arranged in the detection area, and the recording points of the two-dimensional total magnetic field anomaly data are uniformly distributed in the detection area.

[0055] Specifically, after the detection area is determined, a magnetic total field sensor is arranged in the detection area, and a magnetometer mounting platform (a flight platform, a water platform or a handheld device, etc.) moves along several parallel lines, and the magnetometer records the magnetic total field data of the magnetic total field sensor and coordinate information in an interval sampling manner, and finally, after a series of processing such as diurnal variation correction, two-dimensional magnetic total field anomaly data are calculated.

[0056] Exemplarily, the two-dimensional line distribution diagram as shown in Figure 2 is used to collect the magnetic total field data of the magnetic total field sensor, wherein the lines are along the north-south direction, in order to avoid the error introduced by the two-dimensional data gridding, the interval between adjacent measuring points on the line is 0.1 m, and the interval between adjacent lines is also 0.1 m (only the approximate position relationship of the lines is shown, and all the lines are not shown), and the height of the magnetic total field sensor from the ground is 1 m. In order to verify the effectiveness of the present application, a magnetic anomaly body is buried underground in the detection area, with a spatial coordinate of (10, 10, 0.5), a magnetic moment of 2 Am 2 , a magnetic moment inclination of 40° and a magnetic moment declination of 0°, a background geomagnetic field inclination of 30° and a background geomagnetic field declination of 0°. The scalar magnetic total field data of the magnetic total field sensor at each measuring point are obtained, and two-dimensional scalar magnetic anomaly data are calculated based on the scalar magnetic total field data. For the scalar magnetic total field data of the magnetic total field sensor at each measuring point, in the actual application of the present application, the data are collected by a magnetometer or the like; in the algorithm verification stage, the scalar magnetic anomaly data at each measuring point can be directly calculated by simulation calculation, and a preset (such as a mean value of 0 and a standard deviation of 5 nT) Gaussian white noise is added to simulate the actual observation data.

[0057] In the determination of the three-dimensional coordinates of the target body by using the two-dimensional scalar magnetic anomaly data, a mature Euler deconvolution method or other three-dimensional positioning methods without considering the magnetic moment information of the target body can be used.

[0058] In step S102, the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases is obtained in combination with the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data.

[0059] Each set of unit orthogonal bases is obtained by decomposing the two-dimensional scalar magnetic anomaly data generated by the magnetic dipole of the target body at any height; and the modulus of the two-dimensional scalar magnetic anomaly on a set of unit orthogonal bases is defined as the projection of the two-dimensional scalar magnetic anomaly on the set of unit orthogonal bases.

[0060] In step S103, the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases is substituted into a preset overdetermined equation set, and a minimum variance solution of the overdetermined equation set is obtained to obtain the three components of the vector magnetic moment of the target body.

[0061] Wherein, the over-determined equations describe the mathematical relationship between the modulus of the two-dimensional scalar magnetic anomaly data on each unit orthogonal basis and the three components of the vector magnetic moment of the target body.

[0062] The application firstly acquires two-dimensional scalar magnetic anomaly data, and determines the three-dimensional coordinates of the target body by using the two-dimensional scalar magnetic anomaly data, wherein the record points of the two-dimensional scalar magnetic anomaly data are uniformly distributed in the detection area. Then, the modulus of the two-dimensional scalar magnetic anomaly on each group of unit orthogonal bases is acquired in combination with the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data, wherein the unit orthogonal bases are obtained by decomposing the two-dimensional scalar magnetic anomaly data generated by the magnetic dipole of the target body at any height. Finally, the modulus of the two-dimensional scalar magnetic anomaly on each group of unit orthogonal bases is substituted into the preset over-determined equations, and the minimum variance solution of the over-determined equations is acquired to obtain the three components of the vector magnetic moment of the target body. The magnetic moment size, inclination and declination information can be acquired by properly transforming the three components of the vector magnetic moment to realize the magnetic moment calculation of the target body.

[0063] In some embodiments of the application, the process of acquiring the two-dimensional scalar magnetic anomaly data in the above step S101 can include:

[0064] S1, the scalar total magnetic field data of the magnetic total field sensor in the detection area is recorded by using the magnetometer along the grid points of the preset grid to obtain a plurality of data items.

[0065] Wherein, each data item includes the scalar total magnetic field data and the coordinates of the grid point corresponding to the scalar total magnetic field data, and the grid is uniformly distributed in the detection area.

[0066] S2, the scalar total magnetic field data in each data item is subjected to a daily variation correction process to obtain the scalar magnetic anomaly data of each data item.

[0067] Wherein, the purpose of the daily variation correction process is to eliminate the influence of the geomagnetic field daily variation on the observation data (the magnetic total field data). Generally, within the range of 100km 2 , it can be considered that the daily variation is the same, therefore, the daily variation data of the detection area can be observed by itself or obtained from the nearby geomagnetic station, in addition, the daily variation observation results can be plotted into a daily variation curve for reference.

[0068] The two-dimensional scalar magnetic anomaly data is composed of the scalar magnetic anomaly data on each grid point and the coordinates of each grid point. Figure 3 The two-dimensional scalar magnetic anomaly surface graph provided by the embodiments of the application is illustrated.

[0069] In some embodiments of the application, the process of determining the three-dimensional coordinates of the target body by using the two-dimensional scalar magnetic anomaly data in the above step S101 can include:

[0070] The three-dimensional coordinates of the target body are determined by using the two-dimensional scalar magnetic anomaly data, through Euler deconvolution method or a preset target three-dimensional positioning method based on scalar magnetic anomaly.

[0071] The preset target three-dimensional positioning method based on scalar magnetic anomaly is an invention patent application filed on the same day as the present application, and specific reference can be made to the corresponding patent document.

[0072] It can be understood that neither of the two methods uses the magnetic moment information of the target body, so the three-dimensional coordinates can be excluded from the subsequent influence of the vector magnetic moment information.

[0073] In some embodiments of the present application, the process of obtaining the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal basis, in the above step S102 in combination with the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data, comprises:

[0074] S1, the three-dimensional coordinates of the target body are substituted into the expression of each set of unit orthogonal basis of the two-dimensional scalar magnetic anomaly, to obtain the two-dimensional orthogonal basis expression of the target body.

[0075] S2, based on the preset size window and the two-dimensional orthogonal basis expression, the inner product calculation of the two-dimensional scalar magnetic anomaly data is carried out grid by grid, to obtain the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal basis.

[0076] By decomposing the expression of the two-dimensional scalar magnetic gradient generated by the magnetic dipole of the target body at any observation height, a set of unit orthogonal basis function expressions corresponding to the height can be obtained. Based on this, in some embodiments of the present application, the unit orthogonal basis mentioned in the above step S102 can be expressed as the following 5 sets of two-dimensional orthogonal basis expressions:

[0077]

[0078] Each unit orthogonal basis satisfies the following equation:

[0079]

[0080] Wherein, (x0, y0, z0) is the three-dimensional coordinates of the target body, and (x, y, z) is the three-dimensional coordinates of the record point of the two-dimensional scalar magnetic anomaly data.

[0081] Exemplarily, Figures 4 to 8 The two-dimensional unit orthogonal basis surface diagram calculated by each set of unit orthogonal basis functions is provided.

[0082] In some embodiments of the present application, the process of calculating the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases based on the preset size of the window and the two-dimensional orthogonal basis expression, and performing inner product calculation on the two-dimensional scalar magnetic anomaly data point by point, in the step S2, can include:

[0083] The modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases is calculated by using the following equation n (x0,y0):

[0084]

[0085] wherein I and J are the size of the window in the x direction and the y direction respectively, and ΔT(x i ,y j ) is the two-dimensional scalar magnetic anomaly data of the grid point at (x i ,y j ).

[0086] Exemplarily, the size of the window can be set to 5, and then, since the width of the grid is 0.1 m, I and J are both 50.

[0087] In some embodiments of the present application, the process of substituting the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases into the preset overdetermined equation set, and obtaining the minimum variance solution of the overdetermined equation set to obtain the three components of the vector magnetic moment of the target body, in the step S103, can include:

[0088] The three components of the vector magnetic moment of the target body (m x ,m y ,m z ) are calculated by using the following overdetermined equation set:

[0089]

[0090] wherein α i is the modulus of the two-dimensional scalar magnetic anomaly on the i-th set of unit orthogonal bases, dz is the distance from the target body to the observation plane, and I and A are the geomagnetic inclination and the geomagnetic declination of the place where the detection matrix is located respectively. It can be understood that the observation plane is the plane where the scalar magnetic total field data is collected, that is, the plane where the grid is located.

[0091] The three components of the vector magnetic moment of the target body can be obtained by solving the minimum variance solution of the equation set. In the embodiments of the present application, the value of dz is 1.48 m, and the calculation result is (m x = 1.4711 Am 2 , m y = -0.0035 Am 2 , m z = 1.2716 Am 2), after conversion, the magnetic moment size is 1.94 Am, the magnetic moment inclination angle is 40.84°, and the magnetic moment declination angle is -0.14°, which has a small error with the actual value. 2

[0092] The scalar magnetic anomaly-based magnetic moment calculation device provided in the embodiments of the present application is described below. The scalar magnetic anomaly-based magnetic moment calculation device described below can be correspondingly referred to the scalar magnetic anomaly-based magnetic moment calculation method described above.

[0093] Please refer to Figure 9 The scalar magnetic anomaly-based magnetic moment calculation device provided in the embodiments of the present application can include:

[0094] The coordinate calculation unit 21 is configured to acquire two-dimensional scalar magnetic anomaly data, and determine the three-dimensional coordinates of the target body by using the two-dimensional scalar magnetic anomaly data, wherein the recording points of the two-dimensional scalar magnetic anomaly data are uniformly distributed in the detection area.

[0095] The modulus calculation unit 22 is configured to acquire the modulus of the two-dimensional scalar magnetic anomaly on each group of unit orthogonal bases by combining the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data, wherein the unit orthogonal bases are obtained by decomposing the two-dimensional scalar magnetic anomaly data generated by the magnetic dipole of the target body at any height.

[0096] The magnetic moment calculation unit 23 is configured to substitute the modulus of the two-dimensional scalar magnetic anomaly on each group of unit orthogonal bases into a preset overdetermined equation set, and acquire the minimum variance solution of the overdetermined equation set to obtain the three components of the vector magnetic moment of the target body.

[0097] In some embodiments of the present application, the process of acquiring the two-dimensional scalar magnetic anomaly data by the coordinate calculation unit 21 can include:

[0098] The scalar total field data of the magnetic total field sensor in the detection area is recorded by using the magnetometer along the grid points of the preset grid to obtain a plurality of data items, each data item including the scalar total field data and the coordinates of the grid point corresponding to the scalar total field data, and the grid points are uniformly distributed in the detection area.

[0099] The scalar magnetic anomaly data of each data item is obtained by performing the diurnal correction processing on the scalar total field data in each data item.

[0100] The two-dimensional scalar magnetic anomaly data is composed of the scalar magnetic anomaly data on each grid point and the coordinates of each grid point.

[0101] In some embodiments of the present application, the process of determining the three-dimensional coordinates of the target body by the coordinate calculation unit 21 using the two-dimensional scalar magnetic anomaly data can include:

[0102] ​The three-dimensional coordinates of the target body are determined by using the two-dimensional scalar magnetic anomaly data and an Euler deconvolution method or a preset target three-dimensional positioning method based on scalar magnetic anomaly.

[0103] In some embodiments of the present application, the modulus calculation unit 22 obtains the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases by combining the three-dimensional coordinates of the target body with the two-dimensional scalar magnetic anomaly data, and the process includes:

[0104] The three-dimensional coordinates of the target body are substituted into the expression of each set of unit orthogonal bases of the two-dimensional scalar magnetic anomaly to obtain the two-dimensional orthogonal base expression of the target body.

[0105] Based on the preset size window and the two-dimensional orthogonal base expression, the inner product calculation of the two-dimensional scalar magnetic anomaly data is performed grid point by grid point to obtain the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases.

[0106] In some embodiments of the present application, the two-dimensional orthogonal base expression can include:

[0107]

[0108] Each unit orthogonal base satisfies the following equation:

[0109]

[0110] where (x0, y0, z0) is the three-dimensional coordinates of the target body, and (x, y, z) is the three-dimensional coordinates of the record point of the two-dimensional scalar magnetic anomaly data.

[0111] In some embodiments of the present application, the modulus calculation unit 22 obtains the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases by performing the inner product calculation of the two-dimensional scalar magnetic anomaly data grid point by grid point based on the preset size window and the two-dimensional orthogonal base expression, and the process can include:

[0112] The modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases is calculated by using the following equation: n (x0, y0):

[0113]

[0114] where I and J are the size of the window in the x direction and the y direction, respectively, and ΔT(x i ,y j ) is the two-dimensional scalar magnetic anomaly data of the grid point at (x i ,y j ).

[0115] In some embodiments of the present application, the magnetic moment calculation unit 23 substitutes the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases into a preset overdetermined equation set, and obtains a minimum variance solution of the overdetermined equation set to obtain a process of obtaining three components of the vector magnetic moment of the target body, which can include:

[0116] The three components of the vector magnetic moment of the target body (m x ,m y ,m z ) are calculated by using the following overdetermined equation set:

[0117]

[0118] Wherein, α i is the modulus of the two-dimensional scalar magnetic anomaly on the i-th set of unit orthogonal bases, dz is the distance from the target body to the observation plane, I and A are the geomagnetic inclination and geomagnetic declination of the place where the detection matrix is located, respectively.

[0119] The magnetic moment calculation device based on scalar magnetic anomaly provided in the embodiments of the present application can be applied to a magnetic moment calculation device based on scalar magnetic anomaly, such as a computer and the like. Optionally, Figure 10 A hardware structure block diagram of the magnetic moment calculation device based on scalar magnetic anomaly is shown, referring to Figure 10 The hardware structure of the magnetic moment calculation device based on scalar magnetic anomaly can include at least one processor 31, at least one communication interface 32, at least one memory 33 and at least one communication bus 34.

[0120] In the embodiments of the present application, the number of the processor 31, the communication interface 32, the memory 33 and the communication bus 34 is at least one, and the processor 31, the communication interface 32 and the memory 33 complete the communication among each other through the communication bus 34;

[0121] The processor 31 can be a central processing unit CPU, or a specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement the embodiments of the present application, etc.

[0122] The memory 33 can include a high-speed RAM memory, and can also include a non-volatile memory, such as at least one disk memory, etc.

[0123] The memory 33 stores a program, and the processor 31 can call the program stored in the memory 33, and the program is used for:

[0124] Obtain two-dimensional scalar magnetic anomaly data, and determine three-dimensional coordinates of the target body by using the two-dimensional scalar magnetic anomaly data, wherein the recording points of the two-dimensional scalar magnetic anomaly data are uniformly distributed in a detection area;

[0125] Obtain the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases by combining the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data, wherein the unit orthogonal bases are obtained by decomposing the two-dimensional scalar magnetic anomaly data generated by the magnetic dipole of the target body at an arbitrary height;

[0126] Substitute the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases into a preset overdetermined equation set, obtain the minimum variance solution of the overdetermined equation set, and obtain the three components of the vector magnetic moment of the target body.

[0127] Optionally, the refinement function and the expansion function of the program can refer to the description above.

[0128] The embodiment of the application further provides a storage medium which can store a program suitable for being executed by a processor, and the program is used for:

[0129] Obtain two-dimensional scalar magnetic anomaly data, and determine three-dimensional coordinates of the target body by using the two-dimensional scalar magnetic anomaly data, wherein the recording points of the two-dimensional scalar magnetic anomaly data are uniformly distributed in a detection area;

[0130] Obtain the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases by combining the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data, wherein the unit orthogonal bases are obtained by decomposing the two-dimensional scalar magnetic anomaly data generated by the magnetic dipole of the target body at an arbitrary height;

[0131] Substitute the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases into a preset overdetermined equation set, obtain the minimum variance solution of the overdetermined equation set, and obtain the three components of the vector magnetic moment of the target body.

[0132] Optionally, the refinement function and the expansion function of the program can refer to the description above.

[0133] In summary:

[0134] The application firstly acquires two-dimensional scalar magnetic anomaly data, and determines three-dimensional coordinates of a target body by using the two-dimensional scalar magnetic anomaly data, wherein the recording points of the two-dimensional scalar magnetic anomaly data are uniformly distributed in a detection area. Since the three-dimensional coordinates of the target body are determined without using the magnetic moment information of the target body, the three-dimensional coordinates can be excluded from the influence on subsequent vector magnetic moment information. Then, the modulus of two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases is acquired by combining the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data, wherein the unit orthogonal bases are obtained by decomposing the two-dimensional scalar magnetic anomaly data generated by a magnetic dipole of the target body at an arbitrary height. Finally, the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases is substituted into a preset overdetermined equation set, and a minimum variance solution of the overdetermined equation set is acquired to obtain three components of the vector magnetic moment of the target body. Subsequently, the magnetic moment size, inclination and declination information can be acquired by appropriately transforming the three components of the vector magnetic moment to realize the magnetic moment calculation of the target body.

[0135] Finally, it should be noted that the terms such as first and second, etc. are merely used to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "contain" or any other variants thereof are intended to cover non-exclusive inclusion, so that the process, method, article or equipment including a series of elements not only includes those elements, but also includes other elements not explicitly listed or inherent to such process, method, article or equipment. Without more limitations, the element defined by the statement "including a" does not exclude the presence of another identical element in the process, method, article or equipment including the element.

[0136] The various embodiments in the specification are described in a progressive manner, each of which focuses on the difference from other embodiments, and the various embodiments can be combined as needed, and the same and similar parts refer to each other.

[0137] The above description of the disclosed embodiments enables a person skilled in the art to implement or use the application. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the application. Therefore, the application will not be limited to these embodiments shown herein, but will conform to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for calculating a magnetic moment based on a scalar magnetic anomaly, characterized in that, The method comprises the following steps: acquiring two-dimensional scalar magnetic anomaly data, and determining three-dimensional coordinates of a target body by using the two-dimensional scalar magnetic anomaly data, wherein the recording points of the two-dimensional scalar magnetic anomaly data are uniformly distributed in a detection area; combining the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data to acquire a modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases, and the process specifically comprises the following steps: substituting the three-dimensional coordinates of the target body into an expression of each set of unit orthogonal bases of the two-dimensional scalar magnetic anomaly to obtain a two-dimensional orthogonal base expression of the target body; performing inner product calculation on the two-dimensional scalar magnetic anomaly data at each grid point based on a preset window size and the two-dimensional orthogonal base expression to obtain the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases; the unit orthogonal bases are obtained by decomposing two-dimensional scalar magnetic anomaly data generated by a magnetic dipole of the target body at an arbitrary height; substituting the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases into a preset overdetermined equation set, and acquiring a minimum variance solution of the overdetermined equation set to obtain three components of a vector magnetic moment of the target body, and the process specifically comprises the following steps: The vector magnetic moment three components of the target body are calculated by using the following overdetermined equation set : ； wherein, is the modulus of the two-dimensional scalar magnetic anomaly on the i-th set of orthonormal basis, dz is the distance from the target body to the observation plane, I and A are the geomagnetic inclination and geomagnetic declination at the location of the detection matrix, respectively.

2. The method of claim 1, wherein, the two-dimensional orthogonal base expression comprises: ； each unit orthogonal base satisfies the following equation: ； wherein, is a three-dimensional coordinate of the target body, is a three-dimensional coordinate of the recording point of the two-dimensional scalar magnetic anomaly data.

3. The method of claim 2, wherein, the process of performing inner product calculation on the two-dimensional scalar magnetic anomaly data at each grid point based on a preset window size and the two-dimensional orthogonal base expression to obtain the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases comprises: The modulus of the two-dimensional scalar magnetic anomaly in each set of orthonormal basis is calculated by the following equation : ； wherein I, J are the size of the window in x and y direction respectively, is the two-dimensional scalar magnetic anomaly data at the grid point of .

4. The method of claim 1, wherein, the process of acquiring the two-dimensional scalar magnetic anomaly data comprises: using a magnetometer to record scalar total magnetic field data of a magnetic total field sensor in the detection area along grid points of a preset grid to obtain a plurality of data items, each data item comprising the scalar total magnetic field data and coordinates of the grid point corresponding to the scalar total magnetic field data, and the grid points are uniformly distributed in the detection area; performing daily variation correction processing on the scalar total magnetic field data in each data item to obtain scalar magnetic anomaly data of each data item; the two-dimensional scalar magnetic anomaly data is composed of the scalar magnetic anomaly data at each grid point and the coordinates of each grid point.

5. The method of claim 1, wherein, the process of determining the three-dimensional coordinates of the target body by using the two-dimensional scalar magnetic anomaly data comprises: determining the three-dimensional coordinates of the target body by using the two-dimensional scalar magnetic anomaly data through an Euler deconvolution method or a preset target three-dimensional positioning method based on scalar magnetic anomaly.

6. A magnetic moment calculation device based on scalar magnetic anomaly, characterized by, The method comprises the following steps: a coordinate calculation unit is configured to acquire two-dimensional scalar magnetic anomaly data, and determine three-dimensional coordinates of a target body by using the two-dimensional scalar magnetic anomaly data, wherein the recording points of the two-dimensional scalar magnetic anomaly data are uniformly distributed in a detection area; a modulus calculation unit is configured to combine the three-dimensional coordinates of the target body and the two-dimensional scalar magnetic anomaly data to acquire a modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases, and the process specifically comprises the following steps: substituting the three-dimensional coordinates of the target body into an expression of each set of unit orthogonal bases of the two-dimensional scalar magnetic anomaly to obtain a two-dimensional orthogonal base expression of the target body; performing inner product calculation on the two-dimensional scalar magnetic anomaly data at each grid point based on a preset window size and the two-dimensional orthogonal base expression to obtain the modulus of the two-dimensional scalar magnetic anomaly on each set of unit orthogonal bases; the unit orthogonal bases are obtained by decomposing two-dimensional scalar magnetic anomaly data generated by a magnetic dipole of the target body at an arbitrary height; The unit orthogonal bases are obtained by decomposing two-dimensional scalar magnetic anomaly data generated by magnetic dipoles of the target body at any height; The magnetic moment calculation unit is used for substituting the modulus of the two-dimensional scalar magnetic anomaly on each group of unit orthogonal bases into a preset overdetermined equation set, and obtaining a minimum variance solution of the overdetermined equation set to obtain three components of the vector magnetic moment of the target body, and the process specifically includes: The vector magnetic moment three components of the target body are calculated by using the following overdetermined equation set : ； wherein, is the modulus of the two-dimensional scalar magnetic anomaly on the i-th set of orthonormal basis, dz is the distance from the target body to the observation plane, I and A are the geomagnetic inclination and geomagnetic declination at the location of the detection matrix, respectively.

7. A magnetic moment calculation device based on scalar magnetic anomaly, characterized by, The magnetic moment calculation method based on the scalar magnetic anomaly comprises the following steps: A memory and a processor; The memory is used for storing a program; The processor is used for executing the program to realize each step of the magnetic moment calculation method based on the scalar magnetic anomaly according to any one of claims 1-5.

8. A storage medium having stored thereon a computer program, characterized in that The computer program is executed by the processor to realize each step of the magnetic moment calculation method based on the scalar magnetic anomaly according to any one of claims 1-5.

Images