A three-dimensional magnetic anomaly inversion method, electronic device and storage medium

Abstract

An embodiment of the present invention discloses a magnetic anomaly three-dimensional inversion method, an electronic device, and a storage medium. The method may include: for a magnetic target located underground, collecting magnetic anomaly observation data of the magnetic target, and obtaining magnetic observation data according to prior magnetic information associated with the magnetic target; acquiring a weight coefficient preset to represent the credibility of the prior magnetic information; respectively establishing data fitting terms for the magnetic anomaly observation data and the magnetic observation data, and combining the established data fitting terms based on the weight coefficient to obtain a total observation data fitting term; and determining a magnetic anomaly three-dimensional inversion result of the distribution of the magnetic target underground based on the total observation data fitting term. The technical solution of the embodiment of the present invention can reduce the difficulty of adding prior magnetic information in magnetic anomaly three-dimensional inversion, and thus improve the accuracy of magnetic anomaly three-dimensional inversion.

Description

Technical Field

[0001] Embodiments of the present invention relate to the technical field of magnetic exploration, and in particular, to a three-dimensional magnetic anomaly inversion method, an electronic device, and a storage medium. Background Art

[0002] The three-dimensional magnetic anomaly inversion method is one of the important means for magnetic exploration interpretation. The addition of prior magnetic information can reduce the non-uniqueness of inversion.

[0003] However, it is difficult to add prior magnetic information in three-dimensional magnetic anomaly inversion, which affects the accuracy of three-dimensional magnetic anomaly inversion and urgently needs to be solved. Summary of the Invention

[0004] Embodiments of the present invention provide a three-dimensional magnetic anomaly inversion method, an electronic device, and a storage medium to reduce the difficulty of adding prior magnetic information in three-dimensional magnetic anomaly inversion and thus improve the accuracy of three-dimensional magnetic anomaly inversion.

[0005] According to one aspect of the present invention, there is provided a three-dimensional magnetic anomaly inversion method, including:

[0006] For a magnetic target located underground, magnetic anomaly observation data of the magnetic target is collected, and magnetic observation data is obtained according to prior magnetic information associated with the magnetic target;

[0007] A weight coefficient for characterizing the credibility of prior magnetic information is obtained;

[0008] Data fitting terms are respectively established for the magnetic anomaly observation data and the magnetic observation data, and each established data fitting term is combined based on the weight coefficient to obtain a total observation data fitting term;

[0009] Based on the total observation data fitting term, a three-dimensional magnetic anomaly inversion result of the underground distribution of the magnetic target is determined.

[0010] According to another aspect of the present invention, there is provided a three-dimensional magnetic anomaly inversion device, including:

[0011] A magnetic observation data obtaining module, configured to collect magnetic anomaly observation data of a magnetic target located underground, and obtain magnetic observation data according to prior magnetic information associated with the magnetic target;

[0012] A weight coefficient obtaining module, configured to obtain a weight coefficient for characterizing the credibility of prior magnetic information;

[0013] A total observation data fitting term obtaining module, configured to respectively establish data fitting terms for the magnetic observation data and the magnetic anomaly observation data, and combine each established data fitting term based on the weight coefficient to obtain a total observation data fitting term;

[0014] The magnetic anomaly three-dimensional inversion result determination module is used to determine the magnetic anomaly three-dimensional inversion result of the underground distribution of the magnetic target based on the total observation data fitting term.

[0015] According to another aspect of the present invention, there is provided an electronic device, including:

[0016] At least one processor; and a memory communicatively connected to the at least one processor; wherein,

[0017] The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor executes any magnetic anomaly three-dimensional inversion method provided in any embodiment of the present invention.

[0018] According to another aspect of the present invention, there is provided a computer-readable storage medium, on which computer instructions are stored, and the computer instructions are used to implement any magnetic anomaly three-dimensional inversion method provided in any embodiment of the present invention when executed by a processor.

[0019] The technical solution of the embodiment of the present invention is as follows: for a magnetic target located underground, magnetic anomaly observation data of the magnetic target is collected, and magnetic observation data is obtained according to the prior magnetic information associated with the magnetic target, so as to use the prior magnetic information to constrain the magnetic anomaly three-dimensional inversion result of the magnetic anomaly observation data, and reduce the non-uniqueness of the magnetic anomaly three-dimensional inversion result; a weight coefficient representing the credibility of the prior magnetic information is obtained in advance, so as to flexibly adjust the constraint of the prior magnetic information with different credibilities on the magnetic anomaly three-dimensional inversion result through the weight coefficient; data fitting terms are established for the magnetic anomaly observation data and the magnetic observation data respectively, and the established data fitting terms are combined based on the weight coefficient to obtain a total observation data fitting term. The magnetic observation data is used as a kind of observation data, and the total observation data fitting term is jointly constructed with the magnetic anomaly observation data through the weight coefficient, which reduces the difficulty of adding prior magnetic information in the magnetic anomaly three-dimensional inversion; based on the total observation data fitting term, the magnetic anomaly three-dimensional inversion result of the underground distribution of the magnetic target is determined. By solving the total observation data fitting term, an accurate magnetic anomaly three-dimensional inversion result can be obtained. The above technical solution, by using the magnetic observation data obtained based on the prior magnetic information as a kind of observation data, using the weight coefficient to represent the credibility of the magnetic observation data, and jointly constructing the total observation data fitting term with the magnetic anomaly observation data, can reduce the difficulty of adding prior magnetic information in the magnetic anomaly three-dimensional inversion, and thus improve the accuracy of the magnetic anomaly three-dimensional inversion. Especially when the credibilities of the prior magnetic information are different, the accuracy of the magnetic anomaly three-dimensional inversion can be significantly improved.

[0020] It should be understood that the content described in this section is not intended to identify the key or important features of the embodiments of the present invention, nor is it used to limit the scope of the present invention. Other features of the present invention will become readily understood through the following description. Description of the Drawings

[0021] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the accompanying drawings required for the description of the embodiments. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention, and those of ordinary skill in the art can also obtain other accompanying drawings based on these drawings without creative efforts.

[0022] Figure 1 is a flowchart of a three-dimensional magnetic anomaly inversion method provided according to an embodiment of the present invention;

[0023] Figure 2 is a flowchart of another three-dimensional magnetic anomaly inversion method provided according to an embodiment of the present invention;

[0024] Figure 3 is a flowchart of yet another three-dimensional magnetic anomaly inversion method provided according to an embodiment of the present invention;

[0025] Figure 4 is a schematic diagram of a magnetic target model in a specific example of yet another three-dimensional magnetic anomaly inversion method provided according to an embodiment of the present invention;

[0026] Figure 5 is a flowchart of a three-dimensional magnetic anomaly inversion method in a specific example of yet another three-dimensional magnetic anomaly inversion method provided according to an embodiment of the present invention;

[0027] Figure 6 is a magnetic anomaly observation data graph in a specific example of yet another three-dimensional magnetic anomaly inversion method provided according to an embodiment of the present invention;

[0028] Figure 7 is a borehole core magnetic information graph in a specific example of yet another three-dimensional magnetic anomaly inversion method provided according to an embodiment of the present invention;

[0029] Figure 8 is a surface rock magnetic information graph in a specific example of yet another three-dimensional magnetic anomaly inversion method provided according to an embodiment of the present invention;

[0030] Figure 9 is a comparison graph of three-dimensional magnetic anomaly inversion results in a specific example of yet another three-dimensional magnetic anomaly inversion method provided according to an embodiment of the present invention;

[0031] Figure 10 is a structural block diagram of a three-dimensional magnetic anomaly inversion device provided according to an embodiment of the present invention;

[0032] Figure 11 It is a schematic structural diagram of an electronic device for implementing the magnetic anomaly three-dimensional inversion method of the embodiment of the present invention. Detailed implementation manners

[0033] In order to enable those skilled in the art to better understand the solution of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without making creative efforts shall fall within the protection scope of the present invention.

[0034] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and do not necessarily need to be used to describe a specific order or sequence. It should be understood that such data can be interchanged under appropriate circumstances so that the embodiments of the present invention described here can be implemented in an order other than those illustrated or described here. The situations of "target", "original", etc. are similar and will not be elaborated here. In addition, the terms "comprising" and "having" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product or device comprising a series of steps or units does not necessarily have to be limited to those steps or units clearly listed, but may include other steps or units not clearly listed or inherent to these process, method, product or device.

[0035] Before introducing the embodiments of the present invention, the relevant content described in the background art will be exemplarily described to better understand the magnetic anomaly three-dimensional inversion method described in the embodiments of the present invention.

[0036] In the field of magnetic exploration technology, the magnetic anomaly three-dimensional inversion method is one of the important means for magnetic exploration interpretation. The addition of prior magnetic information can reduce the non-uniqueness of the magnetic anomaly three-dimensional inversion result and improve the inversion accuracy. There are mainly various methods for adding prior magnetic information in magnetic anomaly three-dimensional inversion. One is the reference model method, that is, a reference model is established using prior magnetic information and added to the objective function. The inversion result of this method is sensitive to the reference model. Therefore, a high credibility requirement is imposed on the added prior magnetic information. The other is the model weighting method, that is, the model weighting matrix is adjusted through prior magnetic information to adjust the distribution of the magnetic anomaly three-dimensional inversion result. The binding force of this method is weak and is suitable for adding prior magnetic information with low credibility.

[0037] However, the variety of prior magnetic information leads to different credibility levels of prior magnetic information. Common prior magnetic information used in 3D magnetic anomaly inversion includes: surface geological information, borehole lithology information, structural information, tectonic dip information, etc. Each type of prior magnetic information has a different credibility level. When adding prior magnetic information with a large difference in credibility during 3D magnetic anomaly inversion, especially when adding prior magnetic information with a large difference in credibility, the addition difficulty is high, thus reducing the accuracy of 3D magnetic anomaly inversion, which urgently needs to be solved. The technical solution of the present invention is proposed based on this background. By using the magnetic observation data obtained based on prior magnetic information as a type of observation data, using a weight coefficient to represent the credibility of prior magnetic information, and jointly constructing a total observation data fitting term with magnetic anomaly observation data, the difficulty of adding prior magnetic information during 3D magnetic anomaly inversion can be reduced, and furthermore, the accuracy of 3D magnetic anomaly inversion can be improved, especially when the credibility levels of the added prior magnetic information are different, the accuracy of 3D magnetic anomaly inversion can be improved.

[0038] Figure 1 It is a flowchart of a 3D magnetic anomaly inversion method provided by an embodiment of the present invention. This embodiment is applicable to the situation of performing 3D inversion on magnetic anomaly observation data in the presence of prior magnetic information, especially applicable to the situation of performing 3D inversion on magnetic anomaly observation data when the credibility levels of the added prior magnetic information are different. This method can be executed by a 3D magnetic anomaly inversion device provided by an embodiment of the present invention. The device can be implemented in software and / or hardware, and can be integrated on an electronic device, which can be various user terminals or servers.

[0039] See Figure 1 , the method of the embodiment of the present invention specifically includes the following steps:

[0040] S110. For a magnetic target located underground, magnetic anomaly observation data of the magnetic target is collected, and magnetic observation data is obtained according to the prior magnetic information associated with the magnetic target.

[0041] Among them, the magnetic target can be understood as a magnetic object located deep underground that can generate a magnetic field, such as mineral resources like magnetite, or substances such as magnetic rocks.

[0042] Magnetic anomaly observation data can be understood as magnetic field data observed on the surface or in the air at the location of the magnetic target. Optionally, a magnetometer carried on an aircraft such as a helicopter or a drone can be used to perform large-scale and high-efficiency magnetic measurements, thereby obtaining magnetic anomaly observation data.

[0043] The prior magnetic information can be understood as the relevant information about the distribution of magnetic targets that has been obtained before the 3D inversion of magnetic anomalies. For example, the magnetic information of surface rocks obtained by directly measuring the magnetic force on the surface and the magnetic information of borehole cores obtained by measuring the magnetic force of the borehole cores drilled at the location of the magnetic target. The 3D inversion results of magnetic anomalies can be constrained by the prior magnetic information to reduce the non-uniqueness of the 3D inversion results of magnetic anomalies.

[0044] The magnetic observation data can be understood as the magnetic observation data for magnetic targets determined based on the prior magnetic information.

[0045] For magnetic targets located underground, magnetic anomaly observation data of the magnetic targets can be collected from the air, and prior magnetic information associated with the magnetic targets can be obtained through various methods, and then magnetic observation data can be obtained.

[0046] S120. Obtain the weight coefficient that characterizes the credibility of the prior magnetic information set in advance.

[0047] Among them, the types of prior information are complex and diverse and their credibility varies. For each type of prior magnetic information, different credibility can be set artificially according to exploration experience, and the credibility can be represented by different weight coefficients.

[0048] S130. Establish data fitting terms for the magnetic anomaly observation data and the magnetic observation data respectively, and combine the established data fitting terms based on the weight coefficient to obtain the total observation data fitting term.

[0049] Among them, different data fitting terms can be established for the magnetic observation data and the magnetic anomaly observation data respectively. Optionally, the credibility of different data fitting terms can be represented based on the weight coefficient. Finally, different weight coefficients are added to the established data fitting terms and then combined to obtain the total observation data fitting term including all observation data.

[0050] S140. Determine the 3D inversion result of the magnetic anomaly of the underground distribution of the magnetic target based on the total observation data fitting term.

[0051] Among them, based on the total observation data fitting term, for example, solving the total observation data fitting term, and according to the solution result, determining the 3D inversion result of the magnetic anomaly of the underground three-dimensional distribution of the magnetic target.

[0052] In the technical solution of the embodiment of the present invention, for a magnetic target located underground, magnetic anomaly observation data of the magnetic target is collected, and magnetic observation data is obtained according to the prior magnetic information associated with the magnetic target, so as to use the prior magnetic information to constrain the three-dimensional magnetic anomaly inversion result of the magnetic anomaly observation data and reduce the non-uniqueness of the three-dimensional magnetic anomaly inversion result; a weight coefficient representing the credibility of the prior magnetic information is obtained, so as to flexibly adjust the constraint of the prior magnetic information with different credibility on the three-dimensional magnetic anomaly inversion result through the weight coefficient; data fitting terms are established for the magnetic anomaly observation data and the magnetic observation data respectively, and based on the weight coefficient, the established data fitting terms are combined to obtain a total observation data fitting term. The magnetic observation data is used as an observation data, and the total observation data fitting term is jointly constructed with the magnetic anomaly observation data through the weight coefficient, which reduces the difficulty of adding prior magnetic information in the three-dimensional magnetic anomaly inversion; based on the total observation data fitting term, the three-dimensional magnetic anomaly inversion result of the underground distribution of the magnetic target is determined. By solving the total observation data fitting term, an accurate three-dimensional magnetic anomaly inversion result can be obtained. In the above technical solution, by using the magnetic observation data obtained based on the prior magnetic information as an observation data, using the weight coefficient to represent the credibility of the magnetic observation data, and jointly constructing the total observation data fitting term with the magnetic anomaly observation data, the difficulty of adding prior magnetic information in the three-dimensional magnetic anomaly inversion can be reduced, and thus the accuracy of the three-dimensional magnetic anomaly inversion can be improved. Especially when the credibility of the prior magnetic information is different, the accuracy of the three-dimensional magnetic anomaly inversion can be significantly improved.

[0053] An optional technical solution for determining the three-dimensional magnetic anomaly inversion result of the underground three-dimensional distribution of the magnetic target based on the total observation data fitting term includes: establishing an objective function based on the total observation data fitting term; and determining the three-dimensional magnetic anomaly inversion result of the underground three-dimensional distribution of the magnetic target based on the objective function.

[0054] Among them, in the process of determining the three-dimensional magnetic anomaly inversion result of the underground three-dimensional distribution of the magnetic target based on the total observation data fitting term, an objective function corresponding to the total observation data fitting term can be established first. Optionally, a pre-determined model fitting term can be used as a regularization term, and the regularization term is added to the total observation data fitting term to establish the objective function, so as to prevent overfitting through the regularization term. Finally, the three-dimensional magnetic anomaly inversion result is determined by solving the objective function.

[0055] In the above technical solution, by transforming the process of three-dimensional magnetic anomaly inversion into a function problem, the three-dimensional magnetic anomaly inversion result can be accurately obtained by solving the objective function.

[0056] On this basis, optionally, based on the objective function, determine the three-dimensional magnetic anomaly inversion result of the underground three-dimensional distribution of the magnetic target, including: solving the objective function; taking the solution that makes the value of the objective function the smallest as the three-dimensional magnetic anomaly inversion result of the underground three-dimensional distribution of the magnetic target.

[0057] Among them, after obtaining the objective function, solution algorithms such as the conjugate gradient algorithm can be used to solve the objective function, and the solution that makes the value of the objective function the smallest is taken as the three-dimensional magnetic anomaly inversion result of the underground three-dimensional distribution of the magnetic target.

[0058] According to the above technical solution, an accurate three-dimensional magnetic anomaly inversion result can be obtained by solving the optimal solution of the objective function.

[0059] Figure 2 It is a flowchart of another three-dimensional magnetic anomaly inversion method provided in an embodiment of the present invention. This embodiment is optimized based on the above technical solutions. In this embodiment, optionally, each data fitting term includes a magnetic anomaly observation data fitting term and a magnetic property observation data fitting term. Data fitting terms are established for the magnetic anomaly observation data and the magnetic property observation data respectively, and each established data fitting term is combined based on the weight coefficient to obtain the total observation data fitting term, including: establishing a magnetic anomaly observation data fitting term based on the magnetic anomaly observation data, and establishing a magnetic property observation data fitting term based on the magnetic property observation data; multiplying the magnetic property observation data fitting term by the weight coefficient and adding it to the magnetic anomaly observation data fitting term to obtain the total observation data fitting term. Among them, the explanations of the same or corresponding terms as those in the above embodiments are not repeated here.

[0060] See Figure 2 , the method of this embodiment may specifically include the following steps:

[0061] S210. For a magnetic target located underground, collect the magnetic anomaly observation data of the magnetic target and obtain the magnetic property observation data according to the prior magnetic property information associated with the magnetic target.

[0062] S220. Obtain the weight coefficient preset to characterize the credibility of the prior magnetic property information.

[0063] S230. Based on the magnetic anomaly observation data, establish a magnetic anomaly observation data fitting term, and based on the magnetic property observation data, establish a magnetic property observation data fitting term.

[0064] Among them, a magnetic anomaly observation data fitting term for characterizing the data fitting result of the magnetic anomaly observation data can be established based on the magnetic anomaly observation data, and a magnetic property observation data fitting term for characterizing the data fitting result of the magnetic property observation data can be established based on the magnetic property observation data.

[0065] S240. After multiplying the fitting term of the magnetic observation data by a weight coefficient, add it to the fitting term of the magnetic anomaly observation data to obtain the fitting term of the total observation data.

[0066] Among them, after obtaining the fitting term of the magnetic observation data, it can be multiplied by the corresponding weight coefficient to characterize the strength of the binding force of the magnetic observation data on the 3D inversion result of the magnetic anomaly through the weight coefficient, and then added to the fitting term of the magnetic anomaly observation data to obtain the fitting term of the total observation data.

[0067] S250. Based on the fitting term of the total observation data, determine the 3D inversion result of the magnetic anomaly of the magnetic target body's underground distribution.

[0068] The technical solution of the embodiment of the present invention constructs the fitting term of the total observation data by using the magnetic observation data as an observation data together with the magnetic anomaly observation data, and flexibly characterizes the credibility of the magnetic observation data through the weight coefficient, thereby reducing the difficulty of adding prior magnetic information and improving the accuracy of the 3D inversion result of the magnetic anomaly.

[0069] An optional technical solution is to establish a fitting term for the magnetic anomaly observation data based on the magnetic anomaly observation data, including: establishing an underground magnet distribution model, and using forward calculation to determine the model magnetic anomaly data emitted by the underground magnet distribution model; establishing a fitting term for the magnetic anomaly observation data based on the difference between the magnetic anomaly observation data and the model magnetic anomaly data.

[0070] Among them, the underground magnet distribution model can be understood as an initial model of the underground magnetic body distribution established according to information such as geological background knowledge. Optionally, the underground magnet distribution model can be a simple geometric body (such as a sphere, a cuboid, etc.), or a complex combined model. Then, using the forward theory of magnetic prospecting, calculate the model magnetic anomaly data generated by the underground magnet distribution model on the ground surface. Finally, establish a fitting term for the magnetic anomaly observation data based on the difference between the magnetic anomaly observation data and the model magnetic anomaly data.

[0071] The above technical solution establishes a fitting term for the magnetic anomaly observation data through the difference between the magnetic anomaly observation data and the model magnetic anomaly data, and can further determine an accurate 3D inversion result of the magnetic anomaly by performing mathematical operations on the fitting term of the magnetic anomaly observation data.

[0072] Figure 3It is a flowchart of another magnetic anomaly three-dimensional inversion method provided in an embodiment of the present invention. This embodiment is optimized based on the above technical solutions. In this embodiment, optionally, the number of prior magnetic information is at least two, and each of the at least two prior magnetic information corresponds to its own credibility; based on the magnetic observation data, a magnetic observation data fitting term is established, including: for each magnetic observation data in the magnetic observation data corresponding to the at least two prior magnetic information, based on the magnetic observation data, a magnetic observation data fitting term is established; after multiplying the magnetic observation data fitting term by a weight coefficient, it is added to the magnetic anomaly observation data fitting term to obtain a total observation data fitting term, including: for each established magnetic observation data fitting term, multiplying the magnetic observation data fitting term by the weight coefficient corresponding to the magnetic observation data fitting term to obtain a weighted observation data fitting term; based on the obtained at least two weighted observation data fitting terms and the summation result of the magnetic anomaly observation data fitting term, a total observation data fitting term is obtained. Among them, the explanations of the same or corresponding terms as those in the above embodiments are not repeated here.

[0073] See Figure 3 , the method of this embodiment may specifically include the following steps:

[0074] S310. For a magnetic target located underground, magnetic anomaly observation data of the magnetic target is collected, and based on the magnetic anomaly observation data, a magnetic anomaly observation data fitting term is established.

[0075] S320. At least two prior magnetic information associated with the magnetic target is obtained, and for each prior magnetic information in the at least two prior magnetic information, S330 - S350 is executed.

[0076] S330. A weight coefficient representing the credibility of the prior magnetic information set in advance is obtained.

[0077] S340. Based on the prior magnetic information, magnetic observation data is obtained, and based on the magnetic observation data, a magnetic observation data fitting term is established.

[0078] S350. The magnetic observation data fitting term is multiplied by the weight coefficient to obtain a weighted observation data fitting term.

[0079] S360. Based on the summation result of all the obtained weighted observation data fitting terms and the magnetic anomaly observation data fitting term, a total observation data fitting term is obtained.

[0080] S370. Based on the total observation data fitting term, a three-dimensional magnetic anomaly inversion result of the underground distribution of the magnetic target is determined.

[0081] Among them, in the case of a relatively large variety of prior magnetic information, for each of at least two prior magnetic information, it is used as magnetic observation data. Optionally, the at least two prior magnetic information may include surface rock magnetic information and borehole core magnetic information. Among them, through the surface rock magnetic information, the types and contents of magnetic minerals contained in the surface rock and their influence on magnetic anomaly data can be understood; through the borehole core magnetic information, deeper and more continuous geological information underground can be understood. Then, based on the magnetic observation data, a magnetic observation data fitting term is established. Finally, the multiple established magnetic observation data fitting terms are respectively multiplied by their corresponding weight coefficients and added to the magnetic anomaly observation data fitting term to obtain the total observation data fitting term.

[0082] In the technical solution of the embodiment of the present invention, in the case where the prior magnetic information is at least two, each prior magnetic information is processed separately, and the credibility of each prior magnetic information is characterized by the weight coefficient corresponding to each prior magnetic information. In the case where the credibility of the prior magnetic information is different, the prior magnetic information can be added to the 3D magnetic anomaly inversion through a simple method, reducing the difficulty of the 3D magnetic anomaly inversion.

[0083] To better understand the above technical solutions, the following provides an exemplary illustration with specific examples. In this specific example, the true magnetic target model is as Figure 4 shown, which is an inclined vein, inclined eastward, with a top burial depth of 0m and a magnetization intensity of 10 A / m. The flowchart of the 3D magnetic anomaly inversion is as Figure 5 shown, and the specific steps are as follows:

[0084] Step 1: Collect magnetic anomaly observation data (i.e., magnetic anomaly data 1) and prior magnetic information (i.e., surface rock magnetic information and borehole core magnetic information)

[0085] a. Based on this magnetic target model, a magnetometer is used to measure the magnetic anomaly data of the magnetic target model at a height of more than 100m above the ground. The measured magnetic anomaly observation data is as Figure 6 shown. The number of measurement points on the measurement plane is 40×30 = 1200 sampling points, and the sampling interval is 30m.

[0086] b. Drill holes above the magnetic target model to obtain borehole core magnetic information. As Figure 4 shown, the simulated 6 drill holes are numbered J1, J2, J3, J4, J5, and J6. The drill hole positions and trajectories are as Figure 4 shown by the black solid lines. Among them, 3 drill holes pass through the magnetic target model, and 3 drill holes are located outside the magnetic target model. The sampling interval in the drill holes is 10m, and the collected borehole core magnetic information is as Figure 7 shown.

[0087] c. Directly collect the magnetic information of surface rocks on the ground above the magnetic target body model. The collected magnetic information of surface rocks is as Figure 8 shown.

[0088] Step 2: Establish data fitting terms

[0089] a. Establish a magnetic anomaly observation data fitting term (i.e., data fitting term 1) φ based on the magnetic anomaly observation data (i.e., magnetic anomaly data 1) d1 d1 ;

[0090] b. Establish magnetic observation data (i.e., magnetic data 1) d2 based on the magnetic information of surface rocks. Establish the first magnetic observation data fitting term (i.e., data fitting term 2) φ based on the magnetic observation data d2 d2 ;

[0091] c. Establish magnetic observation data (i.e., magnetic data 3) d3 based on the magnetic information of the cores of 6 designed boreholes collected. Establish the second magnetic observation data fitting term (i.e., data fitting term 3) φ based on the magnetic observation data d3 d3 ;

[0092] Step 3: Construct the total observation data fitting term (i.e., data fitting term)

[0093] a. Determine the weight coefficient λ2 of the data fitting term φ d2 according to the credibility of the magnetic information of surface rocks. If the credibility of the magnetic information of surface rocks is high, increase the weight coefficient λ2 to achieve strong constraint on the magnetic information of surface rocks; if the credibility of the magnetic information of surface rocks is low, decrease the weight coefficient λ2 to achieve weak constraint on the magnetic information of surface rocks. Multiply the data fitting term φ d2 by λ2 to obtain. Since the magnetic information of surface rocks in this specific example is established based on the theoretical model and the credibility of the magnetic information of surface rocks is high, a relatively large value is selected for the weight coefficient λ2 here to perform strong constraint on the magnetic information of surface rocks.

[0094] b. Determine the weight coefficient λ3 of the third data fitting term φ d3 according to the credibility of the magnetic information of borehole cores. If the credibility of the magnetic information of borehole cores is high, increase the weight coefficient λ3 to achieve strong constraint on the magnetic information of borehole cores; if the credibility of the magnetic information of borehole cores is low, decrease the weight coefficient λ3 to achieve weak constraint on the magnetic information of borehole cores. Since the magnetic information of borehole cores in this specific example is established based on the theoretical model and the credibility of the magnetic information of borehole cores is high, a relatively large value is selected for the weight coefficient λ3 here to perform strong constraint on the magnetic information of borehole cores.

[0095] c. Combine the three data fitting terms through the weight coefficients to obtain the total data fitting term φ d = φd1 +λ2φ d2 +λ3φ d3 ;

[0096] Step 4: Construct the objective function

[0097] a. Select the model fitting term φ m As the regularization term, together with the total data fitting term φ d to form the objective function φ = φ d + αφ m , where α is the regularization factor;

[0098] Step 5: Determine the 3D magnetic anomaly inversion result

[0099] a. Use the conjugate gradient algorithm to solve the optimal solution of the objective function to obtain the 3D magnetic anomaly inversion result. The specific slices of the 3D magnetic anomaly inversion result are as shown in Figure 9 (a), and the slices of the magnetic target body model are as shown in Figure 9 (b). By comparing with the magnetic target body model, it can be seen that the 3D magnetic anomaly inversion result is relatively accurate.

[0100] In the above specific example, by using the magnetic observation data obtained based on the prior magnetic information as a kind of observation data, using the weight coefficient to represent the credibility of the magnetic observation data, and jointly constructing the total observation data fitting term with the magnetic anomaly observation data, the difficulty of adding prior magnetic information in the 3D magnetic anomaly inversion can be reduced, and then the accuracy of the 3D magnetic anomaly inversion can be improved under different credibility levels of the prior magnetic information.

[0101] Figure 10 is the structural block diagram of the 3D magnetic anomaly inversion device provided in the embodiment of the present invention. This device is used to execute the 3D magnetic anomaly inversion method provided in any of the above embodiments. This device and the 3D magnetic anomaly inversion methods in the above embodiments belong to the same inventive concept. For the details not described in detail in the embodiment of the 3D magnetic anomaly inversion device, reference can be made to the embodiments of the above 3D magnetic anomaly inversion methods. Refer to Figure 10 , this device may specifically include: a magnetic observation data obtaining module 410, a weight coefficient obtaining module 420, a total observation data fitting term obtaining module 430, and a 3D magnetic anomaly inversion result determining module 440. Among them,

[0102] The magnetic observation data obtaining module 410 is configured to collect the magnetic anomaly observation data of the magnetic target body located underground and obtain the magnetic observation data according to the prior magnetic information associated with the magnetic target body;

[0103] The weight coefficient obtaining module 420 is configured to obtain the weight coefficient preset to represent the credibility of the prior magnetic information;

[0104] The total observed data fitting term obtaining module 430 is configured to establish data fitting terms for magnetic anomaly observed data and magnetic observed data respectively, and combine the established data fitting terms based on weight coefficients to obtain the total observed data fitting term;

[0105] The magnetic anomaly three-dimensional inversion result determining module 440 is configured to determine the magnetic anomaly three-dimensional inversion result of the three-dimensional underground distribution of the magnetic target body based on the total observed data fitting term.

[0106] In an alternative technical solution, each data fitting term includes a magnetic anomaly observed data fitting term and a magnetic observed data fitting term. The total observed data fitting term obtaining module 430 includes:

[0107] The magnetic observed data fitting term establishing sub-module is configured to establish a magnetic anomaly observed data fitting term based on the magnetic anomaly observed data, and establish a magnetic observed data fitting term based on the magnetic observed data;

[0108] The total observed data fitting term obtaining sub-module is configured to add the magnetic observed data fitting term multiplied by the weight coefficient to the magnetic anomaly observed data fitting term to obtain the total observed data fitting term.

[0109] On this basis, optionally, the magnetic observed data fitting term establishing sub-module includes:

[0110] The model magnetic anomaly data determining unit is configured to establish an underground magnet distribution model and use forward calculation to determine the model magnetic anomaly data emitted by the underground magnet distribution model;

[0111] The magnetic observed data fitting term establishing unit is configured to establish a magnetic anomaly observed data fitting term based on the difference between the magnetic anomaly observed data and the model magnetic anomaly data.

[0112] In another alternative, the number of prior magnetic information is at least two, and each of the at least two prior magnetic information corresponds to its own credibility; the magnetic observed data fitting term establishing sub-module includes:

[0113] The magnetic observed data fitting term establishing unit is configured to establish a magnetic observed data fitting term for each magnetic observed data corresponding to at least two prior magnetic information based on the magnetic observed data;

[0114] The total observed data fitting term obtaining sub-module includes:

[0115] The weighted observed data fitting term obtaining unit is configured to multiply each established magnetic observed data fitting term by the weight coefficient corresponding to the magnetic observed data fitting term to obtain the weighted observed data fitting term;

[0116] The total observed data fitting term obtaining unit is configured to obtain the total observed data fitting term based on the sum result of at least two weighted observed data fitting terms and the magnetic anomaly observed data fitting term obtained.

[0117] On this basis, optionally, the at least two prior magnetic information includes surface rock magnetic information and borehole core magnetic information.

[0118] In another alternative technical solution, the magnetic anomaly three-dimensional inversion result determination module 440 includes:

[0119] The objective function establishment sub-module is configured to establish an objective function based on the total observed data fitting term;

[0120] The magnetic anomaly three-dimensional inversion result determination sub-module is configured to determine the magnetic anomaly three-dimensional inversion result of the three-dimensional distribution of the magnetic target body underground based on the objective function.

[0121] On this basis, optionally, the objective function establishment sub-module includes:

[0122] The regularization term determination unit is configured to use the pre-determined model fitting term as the regularization term;

[0123] The objective function establishment unit is configured to add the regularization term and the total observed data fitting term to establish the objective function.

[0124] In another alternative, the magnetic anomaly three-dimensional inversion result determination sub-module includes:

[0125] The objective function solving unit is configured to solve the objective function;

[0126] The magnetic anomaly three-dimensional inversion result determination unit is configured to use the solution that makes the objective function take the minimum value as the magnetic anomaly three-dimensional inversion result of the three-dimensional distribution of the magnetic target body underground.

[0127] The magnetic anomaly three-dimensional inversion device provided by the embodiments of the present invention includes a magnetic observation data acquisition module, which acquires magnetic anomaly observation data of a magnetic target located underground and magnetic observation data according to prior magnetic information associated with the magnetic target, so as to use the prior magnetic information to constrain the magnetic anomaly three-dimensional inversion result of the magnetic anomaly observation data and reduce the non-uniqueness of the magnetic anomaly three-dimensional inversion result; a weight coefficient acquisition module, which acquires a weight coefficient preset to represent the credibility of the prior magnetic information, so as to flexibly adjust the constraint of the prior magnetic information with different credibility on the magnetic anomaly three-dimensional inversion result through the weight coefficient; a total observation data fitting term generation module, which respectively establishes data fitting terms for the magnetic anomaly observation data and the magnetic observation data, and combines the established data fitting terms based on the weight coefficient to obtain a total observation data fitting term. By taking the magnetic observation data as an observation data and jointly constructing the total observation data fitting term with the magnetic anomaly observation data through the weight coefficient, the difficulty of adding prior magnetic information in the magnetic anomaly three-dimensional inversion is reduced; a magnetic anomaly three-dimensional inversion result determination module, which determines the magnetic anomaly three-dimensional inversion result of the underground distribution of the magnetic target based on the total observation data fitting term. By solving the total observation data fitting term, an accurate magnetic anomaly three-dimensional inversion result can be obtained. The above device, by taking the magnetic observation data obtained based on the prior magnetic information as an observation data, using the weight coefficient to represent the credibility of the magnetic observation data, and jointly constructing the total observation data fitting term with the magnetic anomaly observation data, can reduce the difficulty of adding prior magnetic information in the magnetic anomaly three-dimensional inversion, and thus improve the accuracy of the magnetic anomaly three-dimensional inversion. Especially when the credibility of the prior magnetic information is different, the accuracy of the magnetic anomaly three-dimensional inversion can be significantly improved.

[0128] The magnetic anomaly three-dimensional inversion device provided by the embodiments of the present invention can execute the magnetic anomaly three-dimensional inversion method provided by any embodiment of the present invention, and has corresponding functional modules and beneficial effects for executing the method.

[0129] It should be noted that in the embodiments of the above magnetic anomaly three-dimensional inversion device, the included units and modules are only divided according to the functional logic, but are not limited to the above division, as long as the corresponding functions can be realized; in addition, the specific names of the functional units are only for the convenience of mutual distinction and do not limit the protection scope of the present invention.

[0130] Figure 11FIG. 0 shows a schematic structural diagram of an electronic device 10 that can be used to implement an embodiment of the present invention. The electronic device is intended to represent various forms of digital computers, such as, for example, laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as, for example, personal digital processors, cellular telephones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the present invention described and / or claimed herein.

[0131] As Figure 11 shown, the electronic device 10 includes at least one processor 11, and a memory communicatively connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc. The memory stores a computer program executable by the at least one processor. The processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 into the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0132] Multiple components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a magnetic disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0133] The processor 11 can be various general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various dedicated artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as the three-dimensional magnetic anomaly inversion method.

[0134] In some embodiments, the three-dimensional magnetic anomaly inversion method can be implemented as a computer program tangibly embodied in a computer-readable storage medium, such as storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed onto the electronic device 10 via ROM 12 and / or communication unit 19. When the computer program is loaded into RAM 13 and executed by the processor 11, one or more steps of the three-dimensional magnetic anomaly inversion method described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to execute the three-dimensional magnetic anomaly inversion method by any other suitable means (e.g., by means of firmware).

[0135] The various embodiments of the systems and techniques described above in this document can be implemented in digital electronic circuitry, integrated circuit systems, field programmable gate arrays (FPGA), application specific integrated circuits (ASIC), application specific standard products (ASSP), systems on a chip (SOC), complex programmable logic devices (CPLD), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special-purpose or general-purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit the data and instructions to the storage system, the at least one input device, and the at least one output device.

[0136] The computer programs for implementing the methods of the present invention can be written in any combination of one or more programming languages. These computer programs can be provided to the processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when the computer programs are executed by the processor, the functions / operations specified in the flowcharts and / or block diagrams are implemented. The computer programs can be executed entirely on the machine, partially on the machine, as a stand-alone software package partially on the machine and partially on a remote machine, or entirely on a remote machine or server.

[0137] In the context of the present invention, a computer-readable storage medium can be a tangible medium that can contain or store a computer program for use by or in connection with an instruction execution system, apparatus, or device. The computer-readable storage medium can include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination of the foregoing. Alternatively, the computer-readable storage medium can be a machine-readable signal medium. More specific examples of the machine-readable storage medium would include an electrical connection based on one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0138] In order to provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or a trackball) by which the user can provide input to the electronic device. Other kinds of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0139] The systems and techniques described herein can be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer having a graphical user interface or a web browser through which the user can interact with an implementation of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: local area network (LAN), wide area network (WAN), blockchain network, and the Internet.

[0140] A computing system may include a client and a server. The client and the server are generally far from each other and usually interact via a communication network. The relationship between the client and the server is created by computer programs that run on respective computers and have a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or a cloud host, which is a host product in the cloud computing service system, solving the defects of difficult management and weak business scalability existing in traditional physical hosts and VPS services.

[0141] It should be understood that various forms of the processes shown above can be used, with steps reordered, added, or deleted. For example, the steps recited in the present invention can be executed in parallel, sequentially, or in a different order, as long as the desired results of the technical solution of the present invention can be achieved, and no limitation is made herein.

[0142] The above specific embodiments do not constitute a limitation on the protection scope of the present invention. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims

1. A three-dimensional inversion method for magnetic anomaly, characterized in that: include: For a magnetic target body located underground, magnetic anomaly observation data of the magnetic target body is collected, and magnetic observation data is obtained according to prior magnetic information associated with the magnetic target body; Obtaining a preset weight coefficient representing the credibility of the prior magnetic information; Establishing data fitting items for the magnetic anomaly observation data and the magnetic observation data respectively, and combining the established data fitting items based on the weight coefficient to obtain a total observation data fitting item; Determining a three-dimensional inversion result of magnetic anomaly of the magnetic target body in the underground distribution based on the total observation data fitting item; Each of the data fitting items includes a magnetic anomaly observation data fitting item and a magnetic observation data fitting item. The data fitting items are established for the magnetic anomaly observation data and the magnetic observation data, respectively, and the established data fitting items are combined based on the weight coefficient to obtain a total observation data fitting item, including: Based on the magnetic anomaly observation data, establishing a magnetic anomaly observation data fitting term, and based on the magnetic observation data, establishing a magnetic observation data fitting term; After multiplying the magnetic observation data fitting term by the weight coefficient, the fitting term is added to the magnetic anomaly observation data fitting term to obtain a total observation data fitting term; The number of the prior magnetic information is at least two, and the at least two prior magnetic information respectively correspond to respective credibility; The step of establishing a magnetic observation data fitting term based on the magnetic observation data comprises: For each magnetic observation data in the magnetic observation data respectively corresponding to at least two pieces of the prior magnetic information, a magnetic observation data fitting item is established based on the magnetic observation data; The method of multiplying the magnetic observation data fitting term by the weight coefficient and adding the resultant to the magnetic anomaly observation data fitting term to obtain a total observation data fitting term comprises: For each of the established magnetic observation data fitting items, multiply the magnetic observation data fitting item by the weight coefficient corresponding to the magnetic observation data fitting item to obtain a weighted observation data fitting item; Based on the sum of at least two of the weighted observation data fitting items and the magnetic anomaly observation data fitting item, a total observation data fitting item is obtained.

2. The method according to claim 1, characterized in that The step of establishing a magnetic anomaly observation data fitting term based on the magnetic anomaly observation data comprises: Establishing an underground magnet distribution model, and determining model magnetic anomaly data emitted by the underground magnet distribution model by forward calculation; Based on the difference between the magnetic anomaly observation data and the model magnetic anomaly data, a magnetic anomaly observation data fitting term is established.

3. The method according to claim 1, characterized in that The at least two prior magnetic information include surface rock magnetic information and drill core magnetic information.

4. The method according to claim 1, characterized in that: The method of determining the three-dimensional inversion result of the magnetic anomaly of the magnetic target body in the underground distribution based on the total observation data fitting item includes: Establishing an objective function based on the total observed data fitting term; Based on the objective function, a three-dimensional inversion result of the magnetic anomaly of the magnetic target body in the underground distribution is determined.

5. The method according to claim 4, characterized in that The establishing of the objective function based on the total observed data fitting term comprises: The predetermined model fitting term is used as a regularization term; The regularization term is added to the total observation data fitting term to establish an objective function.

6. The method according to claim 4, characterized in that Determining the three-dimensional inversion result of the magnetic anomaly of the magnetic target body in the underground distribution based on the objective function includes: solving the objective function; The solution that minimizes the value of the objective function is used as the three-dimensional inversion result of the magnetic anomaly of the magnetic target body in the underground distribution.

7. An electronic device, characterized in that: include: at least one processor; as well as a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor performs the three-dimensional inversion method for magnetic anomaly according to any one of claims 1 to 6.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the three-dimensional inversion method for magnetic anomaly as described in any one of claims 1 to 6 when executed.

Images