Method and device for determining magnetic distribution information based on multi-platform magnetic force data

Abstract

The embodiment of the present application discloses a kind of method and device for determining magnetic distribution information based on multi-platform magnetic force data.The method comprises: obtaining the magnetic force data measured by at least two measuring platforms on the area to be explored;Wherein, the measuring platform includes aerial survey platform, ground survey platform and drilling survey platform;Based on the magnetic force data corresponding to the area to be explored and the preset parameters to be inverted, a target function is constructed;Wherein, the unknown parameters are added to the target function as the parameters to be inverted, and the parameters to be inverted include at least one preset base point difference parameter and the magnetic parameters to be determined;The base point difference parameter is used to adjust the reference gap between different magnetic force data in the target function;Determine the parameter value corresponding to the magnetic parameters to be determined in the target function, and determine the magnetic distribution information generated by the magnetic body contained in the area to be explored based on the parameter value.The technical scheme of the embodiment of the present application can improve the accuracy of magnetic distribution information.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of magnetic exploration technology, and in particular to a method and apparatus for determining magnetic distribution information based on multi-platform magnetic data. Background Technology

[0002] Magnetic exploration is one of the important methods in geophysical exploration, and three-dimensional inversion modeling of magnetic data is a crucial technique in magnetic exploration, enabling the quantitative inversion of the distribution characteristics of subsurface magnetic bodies. Coordinated inversion interpretation of airborne, surface, and borehole data is a developing trend. Because magnetic data measured by various platforms (airborne, surface, and borehole) may have large acquisition time spans, and the baselines of each data point are often different with unknown differences, directly performing coordinated inversion without addressing this issue will lead to difficulties in data fitting and cause the inversion problem to fail to converge.

[0003] In existing technologies, during air-to-ground-well co-inversion, to unify and balance the data, experts typically use their professional knowledge and experience to set baseline differences between different magnetic data from air, ground, and well sources. These baseline differences are then used to achieve co-inversion of the magnetic data from air, ground, and well sources, thereby determining the magnetic distribution information. However, in implementing this invention, it was found that the existing technology has at least the following technical problems: because the baseline differences are determined based on expert experience, the error is relatively large, affecting the inversion effect and resulting in poor accuracy of the determined magnetic distribution information. Summary of the Invention

[0004] This invention provides a method and apparatus for determining magnetic distribution information based on multi-platform magnetic data, in order to improve the accuracy of magnetic distribution information.

[0005] According to one aspect of the present invention, a method for determining magnetic distribution information based on multi-platform magnetic data is provided, comprising:

[0006] Acquire magnetic data from measurements of the area to be explored using at least two measurement platforms; wherein the measurement platforms include an airborne measurement platform, a ground measurement platform, and a borehole measurement platform;

[0007] Based on the magnetic data corresponding to the area to be explored and the preset inversion parameters, an objective function is constructed; wherein, the inversion parameters are added as unknowns to the objective function, and the inversion parameters include at least one preset base point difference parameter and a magnetic parameter to be determined; the base point difference parameter is used in the objective function to adjust the benchmark difference between different magnetic data;

[0008] Determine the parameter value corresponding to the magnetic parameter to be determined in the objective function, and determine the magnetic distribution information generated by the magnetic body contained in the area to be explored based on the parameter value.

[0009] According to another aspect of the present invention, an apparatus for determining magnetic distribution information based on multi-platform magnetic data is provided, the apparatus comprising:

[0010] The data acquisition module is used to acquire magnetic data obtained from measurements of the area to be explored by at least two measurement platforms; wherein, the measurement platforms include an airborne measurement platform, a ground measurement platform, and a borehole measurement platform;

[0011] The function construction module is used to construct an objective function based on the magnetic data corresponding to the area to be explored and the preset inversion parameters; wherein, the inversion parameters are added as unknowns to the objective function, and the inversion parameters include at least one preset base point difference parameter and a magnetic parameter to be determined; the base point difference parameter is used in the objective function to adjust the benchmark difference between different magnetic data;

[0012] The information determination module is used to determine the parameter value corresponding to the magnetic parameter to be determined in the objective function, and to determine the magnetic distribution information generated by the magnetic body contained in the area to be explored based on the parameter value.

[0013] According to another aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0014] At least one processor; and

[0015] A memory communicatively connected to the at least one processor; wherein,

[0016] The memory stores a computer program that can be executed by the at least one processor, which enables the at least one processor to perform the method for determining magnetic distribution information based on multi-platform magnetic data according to any embodiment of the present invention.

[0017] According to another aspect of the present invention, a computer-readable storage medium is provided, the computer-readable storage medium storing computer instructions, the computer instructions being configured to cause a processor to execute and implement the method for determining magnetic distribution information based on multi-platform magnetic data as described in any embodiment of the present invention.

[0018] According to another aspect of the present invention, a computer program product is provided, comprising a computer program that, when executed by a processor, implements the method for determining magnetic distribution information based on multi-platform magnetic data as described in any embodiment of the present invention.

[0019] The technical solution of this invention involves acquiring magnetic data from at least two measurement platforms measuring the area to be explored. These platforms include an aerial measurement platform, a ground measurement platform, and a borehole measurement platform. An objective function is constructed using the magnetic data corresponding to the area to be explored and preset inversion parameters. Since the inversion parameters are added as unknowns to the objective function, these parameters include at least one preset baseline difference parameter and a magnetic parameter to be determined. The baseline difference parameter is used in the objective function to adjust the benchmark difference between different magnetic data. Therefore, there is no need for experts to determine the benchmark difference beforehand; the benchmark difference is directly added as an unknown to the objective function. By determining the parameter value corresponding to the magnetic parameter to be determined in the objective function, the magnetic distribution information generated by the magnetic body contained in the area to be explored is determined based on the parameter value. This technical solution eliminates the need for expert experience to determine the benchmark difference; it directly determines the parameter value corresponding to the magnetic parameter to be determined in the objective function, thereby determining the magnetic distribution information. This avoids errors introduced by the benchmark difference and improves the accuracy of the magnetic distribution information.

[0020] It should be understood that the description in this section is not intended to identify key or essential features of the embodiments of the present invention, nor is it intended to limit the scope of the invention. Other features of the invention will become readily apparent from the following description. Attached Figure Description

[0021] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0022] Figure 1 This is a flowchart of a method for determining magnetic distribution information based on multi-platform magnetic data according to an embodiment of the present invention;

[0023] Figure 2 This is a schematic diagram of a device for determining magnetic distribution information based on multi-platform magnetic data according to an embodiment of the present invention;

[0024] Figure 3 This is a schematic diagram of the structure of an electronic device that implements the method for determining magnetic distribution information based on multi-platform magnetic data according to embodiments of the present invention. Detailed Implementation

[0025] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0026] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "etc." and any variations thereof are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0027] It should be noted that the collection, gathering, updating, analysis, processing, use, transmission, and storage of user personal information involved in this disclosed technical solution all comply with relevant laws and regulations, are used for legitimate purposes, and do not violate public order and good morals. Necessary measures are taken to prevent unauthorized access to user personal information data and to safeguard user personal information security and network security.

[0028] Figure 1 This is a flowchart illustrating a method for determining magnetic distribution information based on multi-platform magnetic data according to an embodiment of the present invention. This embodiment is applicable to determining the distribution of magnetic properties generated by a magnetic body in an area to be explored. The method can be executed by a magnetic distribution information determining device, which can be implemented in hardware and / or software.

[0029] like Figure 1 As shown, the method in this embodiment may specifically include:

[0030] S110. Obtain magnetic data from measurements of the area to be explored by at least two measurement platforms.

[0031] The measurement platforms include airborne measurement platforms, ground measurement platforms, and borehole measurement platforms. The area to be explored is an underground region where the distribution of magnetic properties of magnetic materials needs to be determined. Magnetic data is data reflecting magnetic characteristics obtained by the measurement platforms after being affected by magnetic materials in the area to be explored. For example, airborne magnetic data is data reflecting magnetic characteristics obtained by airborne measurement platforms after being affected by underground magnetic materials in the area to be explored.

[0032] It should be noted that, in order to effectively and comprehensively understand the distribution of underground magnetic materials, magnetic data can be obtained from multiple measurement platforms. For example, magnetic data can be obtained from airborne measurement platforms, ground measurement platforms, and borehole measurement platforms, respectively.

[0033] Specifically, an airborne surveying platform can be an airborne magnetometer deployed on an aircraft; a ground-based surveying platform includes a ground-based magnetometer deployed on the Earth's surface. Magnetic data measured by an airborne surveying platform can be used as airborne magnetic data; large magnetic data acquired by a ground-based surveying platform can be used as ground-based magnetic data; and magnetic data acquired by a borehole surveying platform can be used as borehole magnetic data. Airborne magnetic data, ground-based magnetic data, and borehole magnetic data are all caused by underground magnetic bodies. Airborne magnetic data and ground-based magnetic data are two-dimensional planar data, while borehole magnetic data is one-dimensional data.

[0034] To facilitate the effective and accurate determination of magnetic distribution information using magnetic data, after acquiring magnetic data from at least two measurement platforms measuring the area to be explored, the process also includes: performing data correction processing on the magnetic data to obtain corrected data; performing gridding processing on the corrected data to obtain grid data; and updating the grid data with magnetic data.

[0035] In this embodiment, data quality can be effectively improved and environmental interference can be compensated for by performing data correction on the magnetic force data. Specifically, different data correction processing methods can be used for different magnetic force data, and the data obtained after data correction processing becomes the corrected data.

[0036] It should be noted that, since the magnetic data is discrete and irregularly distributed, it can be converted into data on a regular grid for easier subsequent processing and analysis. Specifically, the calibration data can be preprocessed, including outlier removal and missing value imputation.

[0037] Furthermore, the preprocessed correction data can be gridded to obtain gridded data. Specifically, grid parameters can be determined first. For example, grid parameters include grid spacing and coverage area. Since both airborne magnetic data and ground magnetic data are two-dimensional data and require subsequent collaborative processing, the same grid parameters can be used for gridding of both. Borehole magnetic data is one-dimensional data with dense sampling points. To match the airborne and ground magnetic data, the borehole magnetic data can be thinned so that the point spacing used when dividing the borehole magnetic data matches the grid spacing of the airborne or ground magnetic data. For example, the average point spacing can be a preset multiple of the grid spacing, with the preset multiple ranging from (0.5, 1). The correction data obtained after gridding can then be updated to magnetic data.

[0038] In this embodiment, by calibrating and gridding the magnetic data, the data quality can be effectively improved and environmental interference can be compensated. Furthermore, the magnetic data is converted into data on a regular grid, which facilitates subsequent processing and analysis.

[0039] Optionally, the magnetic data is subjected to data correction processing to obtain corrected data including at least two of the following: normal field correction processing, diurnal variation correction processing, and survey line leveling processing for airborne magnetic data measured by an airborne measurement platform; normal field correction processing and diurnal variation correction processing for ground magnetic data measured by a ground measurement platform; and component conversion processing and normal field correction processing for drilling magnetic data measured by a borehole measurement platform.

[0040] It should be noted that the Earth's magnetic field is a global background field, and its intensity and direction vary in different regions. All magnetic data, including airborne, surface, and well magnetic data, are affected by the Earth's magnetic field, with magnetic anomalies from magnetic bodies being more pronounced. To facilitate subsequent analysis and interpretation, normal field correction methods can be used to process magnetic data from all types of data.

[0041] Furthermore, since both airborne and ground-based magnetic data are collected from the ground, and the Earth's magnetic field is affected by factors such as solar activity, it exhibits diurnal variations—that is, the strength and direction of the Earth's magnetic field change periodically throughout the day. Therefore, diurnal variation correction can be used to eliminate this interference in airborne and ground-based magnetic data, making the magnetic data more accurate and reliable.

[0042] In actual measurements of airborne magnetic data, baseline drift may occur between different survey lines due to factors such as instrument performance and measurement environment, meaning that the overall magnetic field strength of the survey lines may differ. By performing survey line leveling processing on the airborne magnetic data, baseline drift can be eliminated, enabling seamless splicing and comparison of data from different survey lines.

[0043] In this embodiment, the borehole magnetic data measured by the borehole measurement platform is borehole magnetic three-component data, which refers to the magnetic field strength data in two horizontal directions perpendicular to and perpendicular to the borehole axis. To transform the collected borehole magnetic three-component data into a unified coordinate system, the borehole magnetic data can undergo component transformation processing. For example, the three-component data in the local coordinate system can be transformed into the geographic coordinate system for joint analysis with surface magnetic data.

[0044] In this embodiment, by applying matching correction processing to different magnetic force data, the effectiveness of the magnetic force data is improved, and subsequent joint analysis is facilitated.

[0045] S120. Based on the magnetic data corresponding to the area to be explored and the preset inversion parameters, construct an objective function; wherein, the inversion parameters are added as unknowns to the objective function, including at least one preset base point difference parameter and a magnetic parameter to be determined; the base point difference parameter is used in the objective function to adjust the benchmark difference between different magnetic data.

[0046] Among them, the magnetic parameter to be determined can be a parameter used to represent the magnetic intensity. The larger the value of the magnetic parameter to be determined, the greater the magnetic intensity; conversely, the smaller the value, the smaller the magnetic intensity.

[0047] In this embodiment, the area to be explored can be treated as a whole, and the quantitative relationship between the area to be explored and each measuring point in the magnetic data can be determined. For example, the quantitative relationship between airborne magnetic data and the area to be explored is represented by B1, the quantitative relationship between ground magnetic data and the area to be explored is represented by B2, and the quantitative relationship between borehole magnetic data and the area to be explored is represented by B3. The baseline difference parameter may include a first baseline difference parameter η1 for adjusting the baseline difference between the airborne magnetic data and the ground magnetic data, and a second baseline difference parameter η2 for adjusting the baseline difference between the airborne magnetic data and the borehole magnetic data.

[0048] Furthermore, the correspondence between the magnetic data of the area to be explored and the magnetic parameters to be determined is as follows:

[0049]

[0050] Where κ represents the magnetic parameter to be determined, T1 represents airborne magnetic data, T2 represents ground magnetic data, and T3 represents borehole magnetic data. An objective function can be constructed based on the above correspondence.

[0051] Optionally, based on the magnetic data corresponding to the area to be explored and the preset inversion parameters, an objective function is constructed, including: determining the forward modeling operator matrix based on the magnetic data corresponding to the area to be explored; dividing the area to be explored into blocks, and generating a model weighting function matrix based on the location of at least one block, the location of measuring points, and the size of the blocks; and constructing the objective function based on the forward modeling operator matrix, the magnetic data, the model weighting function matrix, and the inversion parameters.

[0052] In this embodiment, the area to be explored is a three-dimensional space. In order to more clearly and meticulously determine the magnetic distribution of the magnetic body corresponding to the area to be explored, the area to be explored can be divided into regular grids into closely arranged cubes. By determining the magnetism of each cube, the magnetic distribution of the underground three-dimensional space can be obtained. This cube can be called a block.

[0053] To obtain the magnetic parameters of each block, a quantitative relationship needs to be established between the magnetic parameters of each block and the magnetic data. This quantitative relationship can be represented by forward operators, and the forward operators corresponding to each block form a forward operator matrix. Optionally, the forward operator for each block can be associated with the relative position between each measurement point corresponding to the magnetic data and the block, the magnetization direction of the block, and the normal field direction of the measurement point. The measurement point can be the endpoint of each grid after the magnetic data is gridded; the magnetization direction of the block and the normal field direction of each measurement point can be the same, both being preset directions.

[0054] In practical implementation, for each measuring point, the forward modeling operator for each block relative to that measuring point is determined; all the forward modeling operators form a forward modeling operator matrix. For example, if the amount of magnetic data is N, and the number of blocks after the exploration area is divided is M, then the size of the forward modeling operator matrix can be N×M. Here, N and M are both positive integers. For magnetic data obtained from different measuring platforms, the corresponding forward modeling operator matrix can be determined. Integrating the magnetic data obtained from different measuring platforms, the forward modeling formula is as follows:

[0055]

[0056] Where T1 represents airborne magnetic data, T2 represents ground magnetic data, T3 represents borehole magnetic data, κ represents the magnetic parameter to be determined, A1 represents the forward operator matrix of T1, A2 represents the forward operator matrix of T2, and A3 represents the forward operator matrix of T3. η1 represents the first base point difference parameter, and η2 represents the second base point difference parameter.

[0057] Furthermore, based on the above forward modeling formula, the data fitting term φ of the objective function... d for:

[0058]

[0059] Furthermore, a model fitting term can be added to the objective function as a regularization term. After adding the model fitting term, the resulting objective function... as follows:

[0060]

[0061] Here, α represents a pre-defined regularization factor, and W represents the model weighting function matrix, which is an M×M diagonal matrix. Each element in the model weighting function matrix corresponds to a block. The model weighting function matrix is ​​crucial for the inversion of airborne, ground, and borehole magnetic survey data. A combined weighting function is proposed, which includes a depth weighting function and a distance weighting function, aiming to counteract the geometric attenuation effect of the magnetic kernel function with increasing depth (for airborne and ground survey points) and with increasing distance (for borehole survey points), respectively.

[0062] In this embodiment, the depth weighting function matrix U m The matrix elements are defined as follows:

[0063]

[0064] Among them, z j Let be the burial depth of the j-th block unit; a small constant z0 is added to avoid singularities, and its value is usually half the block height. The value of j is a positive integer in the range [1, M].

[0065] Furthermore, the distance weighting function matrix V m Defined by the distance between the grid cells and the observation points, the elements of its matrix are defined as follows:

[0066]

[0067] Among them, R ij The distance between the j-th block and the i-th measuring point is given; a small constant R is added. o To avoid singularities, its value is usually taken as half the width of the block unit. The parameter β has a range of 0.5 ≤ β < 1.5, and is usually set to 1.0. The value of j is a positive integer in the range [1, M].

[0068] It should be noted that when the borehole passes through the target body, the aforementioned distance weighting function may lead to inaccurate inversion results at the borehole crossing location, especially since a magnetization intensity distribution that should exist at this location is difficult to reflect in the results. Therefore, for boreholes passing through the target body, the measurement points within them will not participate in the calculation of the weighting function matrix.

[0069] Finally, the model weighting function matrix W combines the depth weighting function matrix and the distance weighting function matrix, and its matrix elements w j by u i and v j Combining and defined as:

[0070] w j =u j v j

[0071] The model weighting function matrix W is the weighting function matrix used for the inversion of airborne, ground, and borehole magnetic survey data.

[0072] It should be noted that incorporating the baseline difference as an unknown into the objective function reduces the complexity of the open-field well inversion method, avoids the influence of baseline difference errors estimated by experts on the inversion results, and enables the inversion to converge successfully.

[0073] It should be noted that existing technologies require experts to determine the baseline difference, which limits the application of collaborative inversion of air-ground-well magnetic data. The objective function provided in this embodiment treats the baseline difference parameter as an unknown variable. By solving the objective function, not only the distribution of underground magnetic parameters can be obtained, but also the baseline difference. This method avoids the problem of pre-setting the baseline difference during collaborative inversion of air-ground-well magnetic data, promotes data fitting, solves the problem of convergence difficulties in inversion, and yields reasonable inversion results. Simultaneously, it reduces the technical difficulty of collaborative inversion of air-ground-well magnetic data and expands the application scope of collaborative inversion of air-ground-well magnetic data.

[0074] In this embodiment, weighting factors can be introduced to balance the differences in data precision. For example, the weighting factors include λ1 and λ2, where λ1 and λ2 represent the degree of precision of ground magnetic data compared to airborne magnetic data, and the degree of precision of borehole magnetic data, respectively. Higher data precision results in a higher weighting factor; conversely, lower data precision results in a lower weighting factor. For example, the value of λ1 is generally between 1 and 0.1; since the measurement precision of borehole data is lower, the value of λ2 is generally between 0.01 and 0.0001.

[0075] Optionally, the objective function includes:

[0076]

[0077] in, Let T1 represent airborne magnetic data, T2 represent ground magnetic data, T3 represent borehole magnetic data, λ1 and λ2 represent pre-set weighting factors, η1 represent the first base point difference parameter, η2 represent the second base point difference parameter, α represent the pre-set regularization factor, W represent the model weighting function matrix, κ represent the magnetic parameters to be determined, A1 represent the forward operator matrix of T1, A2 represent the forward operator matrix of T2, and A3 represent the forward operator matrix of T3.

[0078] S130. Determine the parameter values ​​corresponding to the magnetic parameters to be determined in the objective function, and determine the magnetic distribution information generated by the magnetic bodies contained in the area to be explored based on the parameter values.

[0079] The parameter values ​​are used to reflect the magnitude of the magnetic intensity generated by the magnetic material.

[0080] In this embodiment, for each block, the parameter value of the corresponding magnetic parameter to be determined can be obtained. Based on the position of each block in the area to be explored and the corresponding parameter value, the magnetic distribution information of the magnetic bodies contained in the area to be explored can be determined. This magnetic distribution information reflects the position and range of the magnetic bodies contained in the area to be explored, as well as the variation in the magnetic intensity generated by these magnetic bodies within the area.

[0081] Optionally, determining the parameter values ​​corresponding to the magnetic parameters to be determined in the objective function includes: determining the parameter values ​​corresponding to the magnetic parameters to be determined in the objective function based on at least one optimization algorithm; wherein the optimization algorithm includes at least one of the quasi-Newton method, Bayesian optimization method, and particle swarm optimization algorithm.

[0082] In this embodiment, an optimization algorithm is used to find the parameters to be inverted that minimize the objective function. The values ​​of the undetermined magnetic parameters are then extracted from the identified parameters. Geophysical inversion is achieved using the parameter values ​​corresponding to the identified undetermined magnetic parameters.

[0083] This embodiment provides a variety of optimization algorithms to quickly and effectively determine the parameter values ​​corresponding to the magnetic parameters to be determined in the objective function.

[0084] The technical solution of this invention involves acquiring magnetic data from at least two measurement platforms measuring the area to be explored. These platforms include an aerial measurement platform, a ground measurement platform, and a borehole measurement platform. An objective function is constructed using the magnetic data corresponding to the area to be explored and preset inversion parameters. Since the inversion parameters are added as unknowns to the objective function, these parameters include at least one preset baseline difference parameter and a magnetic parameter to be determined. The baseline difference parameter is used in the objective function to adjust the benchmark difference between different magnetic data. Therefore, there is no need for experts to determine the benchmark difference beforehand; the benchmark difference is directly added as an unknown to the objective function. By determining the parameter value corresponding to the magnetic parameter to be determined in the objective function, the magnetic distribution information generated by the magnetic body contained in the area to be explored is determined based on the parameter value. This technical solution eliminates the need for expert experience to determine the benchmark difference; it directly determines the parameter value corresponding to the magnetic parameter to be determined in the objective function, thereby determining the magnetic distribution information. This avoids errors introduced by the benchmark difference and improves the accuracy of the magnetic distribution information.

[0085] Figure 2 This is a schematic diagram of a device for determining magnetic distribution information based on multi-platform magnetic data according to an embodiment of the present invention. This device is used to execute the magnetic distribution information determination method provided in any of the above embodiments. This device and the magnetic distribution information determination methods of the above embodiments belong to the same inventive concept. Details not described in detail in the embodiments of the magnetic distribution information determination device can be found in the embodiments of the magnetic distribution information determination methods described above. Figure 2 As shown, the device includes:

[0086] Data acquisition module 10 is used to acquire magnetic data obtained from measurements of the area to be explored by at least two measurement platforms; wherein, the measurement platforms include an aerial measurement platform, a ground measurement platform, and a borehole measurement platform;

[0087] The function construction module 11 is used to construct an objective function based on the magnetic data corresponding to the area to be explored and the preset inversion parameters. The inversion parameters are added to the objective function as unknowns. The inversion parameters include at least one preset base point difference parameter and a magnetic parameter to be determined. The base point difference parameter is used in the objective function to adjust the benchmark difference between different magnetic data.

[0088] The information determination module 12 is used to determine the parameter values ​​corresponding to the magnetic parameters to be determined in the objective function, and to determine the magnetic distribution information generated by the magnetic bodies contained in the area to be explored based on the parameter values.

[0089] Based on any optional technical solution in the embodiments of the present invention, optionally, the function construction module 11 includes:

[0090] The matrix determination unit is used to determine the forward operator matrix based on the magnetic data corresponding to the area to be explored;

[0091] Block division unit is used to divide the area to be explored into blocks. Based on the burial depth and height values ​​of at least one block obtained from the division, a model weighting function matrix is ​​generated.

[0092] The function construction unit is used to construct the objective function based on the forward operator matrix, magnetic data, model weighted function matrix, and parameters to be inverted.

[0093] Based on any optional technical solution in the embodiments of the present invention, the objective function may optionally include:

[0094]

[0095] in, Let T1 represent airborne magnetic data, T2 represent ground magnetic data, T3 represent borehole magnetic data, λ1 and λ2 represent pre-set weighting factors, η1 represent the first base point difference parameter, η2 represent the second base point difference parameter, α represent the pre-set regularization factor, W represent the model weighting function matrix, κ represent the magnetic parameters to be determined, A1 represent the forward operator matrix of T1, A2 represent the forward operator matrix of T2, and A3 represent the forward operator matrix of T3.

[0096] Based on any optional technical solution in the embodiments of the present invention, the information determination module 12 may optionally include:

[0097] The parameter value determination unit is used to determine the parameter values ​​corresponding to the magnetic parameters to be determined in the objective function based on at least one optimization algorithm.

[0098] The optimization algorithm includes at least one of the following: quasi-Newton method, Bayesian optimization method, and particle swarm optimization algorithm.

[0099] Optionally, based on any of the optional technical solutions in the embodiments of the present invention, it further includes:

[0100] The correction processing module is used to perform data correction processing on the magnetic data obtained from the measurement of the area to be explored by at least two measurement platforms to obtain corrected data.

[0101] The grid processing module is used to perform grid processing on the correction data to obtain grid data, and then update the grid data with magnetic data.

[0102] Based on any optional technical solution in the embodiments of the present invention, the correction processing module may optionally include at least two of the following units:

[0103] The first processing unit is used to perform normal field correction processing, diurnal variation correction processing, and survey line leveling processing on the airborne magnetic data obtained by the airborne measurement platform.

[0104] The second processing unit is used to perform normal field correction processing and diurnal variation correction processing on the ground magnetic data obtained by the ground measurement platform.

[0105] The third processing unit is used to perform component conversion and normal field correction processing on the drilling magnetic data obtained from the borehole measurement platform.

[0106] The technical solution of this invention involves acquiring magnetic data from at least two measurement platforms measuring the area to be explored. These platforms include an aerial measurement platform, a ground measurement platform, and a borehole measurement platform. An objective function is constructed using the magnetic data corresponding to the area to be explored and preset inversion parameters. Since the inversion parameters are added as unknowns to the objective function, these parameters include at least one preset baseline difference parameter and a magnetic parameter to be determined. The baseline difference parameter is used in the objective function to adjust the benchmark difference between different magnetic data. Therefore, there is no need for experts to determine the benchmark difference beforehand; the benchmark difference is directly added as an unknown to the objective function. By determining the parameter value corresponding to the magnetic parameter to be determined in the objective function, the magnetic distribution information generated by the magnetic body contained in the area to be explored is determined based on the parameter value. This technical solution eliminates the need for expert experience to determine the benchmark difference; it directly determines the parameter value corresponding to the magnetic parameter to be determined in the objective function, thereby determining the magnetic distribution information. This avoids errors introduced by the benchmark difference and improves the accuracy of the magnetic distribution information.

[0107] It is worth noting that in the embodiments of the magnetic distribution information determination device described above, the various units and modules included are only divided according to functional logic, but are not limited to the above division, as long as the corresponding functions can be achieved; in addition, the specific names of each functional unit are only for easy differentiation and are not used to limit the scope of protection of the present invention.

[0108] Figure 3This is a schematic diagram of the structure of an electronic device that implements the method for determining magnetic distribution information based on multi-platform magnetic data according to embodiments of the present invention. The electronic device is intended to represent various forms of digital computers, such as 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 personal digital processors, cellular phones, smartphones, wearable devices (e.g., 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 invention described and / or claimed herein.

[0109] like Figure 3 As shown, the electronic device 20 includes at least one processor 21 and a memory, such as a read-only memory (ROM) 22 or a random access memory (RAM) 23, communicatively connected to the at least one processor 21. The memory stores computer programs executable by the at least one processor. The processor 21 can perform various appropriate actions and processes based on the computer program stored in the ROM 22 or loaded from storage unit 28 into the RAM 23. The RAM 23 can also store various programs and data required for the operation of the electronic device 20. The processor 21, ROM 22, and RAM 23 are interconnected via a bus 24. An input / output (I / O) interface 25 is also connected to the bus 24.

[0110] Multiple components in electronic device 20 are connected to I / O interface 25, including: input unit 26, such as keyboard, mouse, etc.; output unit 27, such as various types of monitors, speakers, etc.; storage unit 28, such as disk, optical disk, etc.; and communication unit 29, such as network card, modem, wireless transceiver, etc. Communication unit 29 allows electronic device 20 to exchange information / data with other devices through computer networks such as the Internet and / or various telecommunications networks.

[0111] Processor 21 can be a variety of general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of processor 21 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special-purpose artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. Processor 21 performs the various methods and processes described above, such as methods for determining magnetic distribution information based on multi-platform magnetic data.

[0112] In some embodiments, the method for determining magnetic distribution information based on multi-platform magnetic data can be implemented as a computer program tangibly contained in a computer-readable storage medium, such as storage unit 28. In some embodiments, part or all of the computer program can be loaded and / or installed on electronic device 20 via ROM 22 and / or communication unit 29. When the computer program is loaded into RAM 23 and executed by processor 21, one or more steps of the method for determining magnetic distribution information based on multi-platform magnetic data described above can be performed. Alternatively, in other embodiments, processor 21 can be configured to perform the method for determining magnetic distribution information based on multi-platform magnetic data by any other suitable means (e.g., by means of firmware).

[0113] Various embodiments of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field-programmable gate arrays (FPGAs), application-specific integrated circuits (ASICs), application-specific standard products (ASSPs), systems-on-a-chip (SoCs), complex programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments may include implementations in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which may be a dedicated or general-purpose programmable processor, capable of receiving data and instructions from a storage system, at least one input device, and at least one output device, and transferring data and instructions to the storage system, the at least one input device, and the at least one output device.

[0114] Computer programs used to implement the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when executed by the processor, the computer programs cause the functions / operations specified in the flowcharts and / or block diagrams to be performed. The computer programs may be executed entirely on a machine, partially on a machine, or as a standalone software package, partially on a machine and partially on a remote machine, or entirely on a remote machine or server.

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

[0116] 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 pointing device (e.g., a mouse or trackball) through which the user provides input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, 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 sound input, voice input, or tactile input).

[0117] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as data servers), or computing systems that include middleware components (e.g., application servers), or computing systems that include frontend components (e.g., user computers with graphical user interfaces or web browsers through which users can interact with implementations of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., communication networks). Examples of communication networks include local area networks (LANs), wide area networks (WANs), blockchain networks, and the Internet.

[0118] A computing system can include clients and servers. Clients and servers are generally located far apart and typically interact through a communication network. The client-server relationship is created by computer programs running on the respective computers and having a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or cloud host, which is a hosting product within the cloud computing service system to address the shortcomings of traditional physical hosts and VPS services, such as high management difficulty and weak business scalability.

[0119] In particular, according to embodiments of the present invention, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, embodiments of the present invention include a computer program product comprising a computer program carried on a non-transitory computer-readable medium, the computer program containing program code for performing the methods shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via communication unit 29, or installed from storage unit 28, or installed from ROM 22. When the computer program is executed by processor 21, it performs the functions defined in the methods of the embodiments of the present invention.

[0120] In implementing the computer program product, computer program code for performing the operations of this invention can be written in one or more programming languages ​​or a combination thereof. Programming languages ​​include object-oriented programming languages ​​such as Java, Smalltalk, and C++, as well as conventional procedural programming languages ​​such as C or similar languages. The program code can be executed entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving remote computers, the remote computer can be connected to the user's computer via any type of network—including a local area network (LAN) or a wide area network (WAN)—or can be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0121] It should be understood that the various forms of processes shown above can be used, with steps reordered, added, or deleted. For example, the steps described in this invention can be executed in parallel, sequentially, or in different orders, as long as the desired result of the technical solution of this invention can be achieved, and this is not limited herein.

[0122] The specific embodiments described above do not constitute a limitation on the scope of protection of this 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 principles of this invention should be included within the scope of protection of this invention.

Claims

1. A method for determining magnetic distribution information based on multi-platform magnetic data, characterized in that, include: Magnetic data from measurements of the area to be explored obtained by multiple measurement platforms; wherein, the measurement platforms include an airborne measurement platform, a ground measurement platform, and a borehole measurement platform; Based on the magnetic data corresponding to the area to be explored and the preset inversion parameters, an objective function is constructed. The inversion parameters are added as unknowns to the objective function, including at least one preset baseline difference parameter and a magnetic parameter to be determined. The baseline difference parameter is used in the objective function to adjust the benchmark difference between different magnetic data. The baseline difference parameter includes a first baseline difference parameter for adjusting the benchmark difference between airborne magnetic data and ground magnetic data, and a second baseline difference parameter for adjusting the benchmark difference between airborne magnetic data and borehole magnetic data. The magnetic parameter to be determined is a parameter representing magnetic intensity. Determine the parameter value corresponding to the magnetic parameter to be determined in the objective function, and determine the magnetic distribution information generated by the magnetic body contained in the area to be explored based on the parameter value; The objective function includes: in, Let be the objective function. Represents airborne magnetic data, Represents ground magnetic data, This represents borehole magnetic data. and These represent the pre-defined weighting factors. This represents the first base point difference parameter. This represents the second base point difference parameter. This represents the pre-defined regularization factor. Let κ represent the model weighting function matrix, and let κ represent the magnetic parameters to be determined. express The forward operator matrix, express The forward operator matrix, express The forward operator matrix.

2. The method according to claim 1, characterized in that, The objective function is constructed based on the magnetic data corresponding to the area to be explored and the preset inversion parameters, including: Based on the magnetic data corresponding to the area to be explored, determine the forward operation matrix; The area to be explored is divided into blocks, and a model weighting function matrix is ​​generated based on the location of at least one block, the location of measuring points, and the size of the block. Based on the forward modeling operator matrix, the magnetic data, the model weighting function matrix, and the parameters to be inverted, an objective function is constructed.

3. The method according to claim 1, characterized in that, Determining the parameter value corresponding to the magnetic parameter to be determined in the objective function includes: Based on at least one optimization algorithm, determine the parameter values ​​corresponding to the magnetic parameters to be determined in the objective function; The optimization algorithm includes at least one of the following: quasi-Newton method, Bayesian optimization method, and particle swarm optimization algorithm.

4. The method according to claim 1, characterized in that, After acquiring the magnetic data obtained from measurements of the area to be explored by multiple measurement platforms, the method further includes: The magnetic force data is subjected to data correction processing to obtain corrected data; The correction data is processed into a grid to obtain grid data, and the grid data is then updated with the magnetic data.

5. The method according to claim 4, characterized in that, The magnetic data is corrected to obtain corrected data including at least two of the following: The airborne magnetic data obtained from the airborne measurement platform are subjected to normal field correction processing, diurnal variation correction processing, and survey line leveling processing. The ground magnetic data obtained from the ground measurement platform are subjected to normal field correction and diurnal variation correction. The drilling magnetic data obtained from the borehole measurement platform are subjected to component conversion processing and normal field correction processing.

6. A device for determining magnetic distribution information based on multi-platform magnetic data, characterized in that, include: The data acquisition module is used to acquire magnetic data obtained from measurements of the area to be explored by multiple measurement platforms; wherein, the measurement platforms include an airborne measurement platform, a ground measurement platform, and a borehole measurement platform; A function construction module is used to construct an objective function based on the magnetic data corresponding to the area to be explored and preset inversion parameters. The inversion parameters are added as unknowns to the objective function, including at least one preset baseline difference parameter and a magnetic parameter to be determined. The baseline difference parameter is used in the objective function to adjust the benchmark difference between different magnetic data. The baseline difference parameter includes a first baseline difference parameter for adjusting the benchmark difference between airborne magnetic data and ground magnetic data, and a second baseline difference parameter for adjusting the benchmark difference between airborne magnetic data and borehole magnetic data. The magnetic parameter to be determined is a parameter representing magnetic intensity. The information determination module is used to determine the parameter value corresponding to the magnetic parameter to be determined in the objective function, and to determine the magnetic distribution information generated by the magnetic body contained in the area to be explored based on the parameter value; The objective function includes: in, Let be the objective function. Represents airborne magnetic data, Represents ground magnetic data, This represents borehole magnetic data. and These represent the pre-defined weighting factors. This represents the first base point difference parameter. This represents the second base point difference parameter. This represents the pre-defined regularization factor. Let κ represent the model weighting function matrix, and let κ represent the magnetic parameters to be determined. express The forward operator matrix, express The forward operator matrix, express The forward operator matrix.

7. An electronic device, characterized in that, The electronic device includes: At least one processor; and A memory communicatively connected to the at least one processor; wherein, The memory stores a computer program that can be executed by the at least one processor, the computer program being executed by the at least one processor to enable the at least one processor to perform the method for determining magnetic distribution information based on multi-platform magnetic data as described in any one of claims 1-5.

8. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions that cause a processor to execute the method for determining magnetic distribution information based on multi-platform magnetic data as described in any one of claims 1-5.

9. A computer program product comprising a computer program that, when executed by a processor, implements the method for determining magnetic distribution information based on multi-platform magnetic data according to any one of claims 1-5.

Images