Magnetic resistivity method forward modeling method and device, electronic equipment and storage medium

Abstract

The invention provides a magnetoresistivity method forward modeling simulation method and device, electronic equipment and a storage medium, and belongs to the technical field of geophysical exploration, and the method comprises the steps: obtaining the modeling related data of a to-be-measured geological region and the actual measurement value of the resistivity of each position of the to-be-measured geological region; establishing a three-dimensional structure model of the geological area to be measured based on the modeling related data; assigning the resistivity of the three-dimensional structure model based on the measured value of the resistivity to obtain a three-dimensional resistivity model of the geological area to be measured; and determining three-dimensional space magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model. According to the magnetoresistivity method forward modeling method provided by the invention, the three-dimensional structure model of the geological area to be measured is established, the three-dimensional resistivity model is established on the basis of the three-dimensional structure model, and the three-dimensional space magnetic field distribution of the geological area to be measured is determined according to the three-dimensional resistivity model. Through a magnetoresistivity forward modeling method, calculation of three-dimensional space magnetic field distribution of a geological area to be measured is realized.

Description

Technical Field

[0001] The present invention relates to the technical field of geophysical exploration, and particularly to a forward simulation method, device, electronic device and storage medium for magneto-resistivity method. Background Art

[0002] The field of geophysical exploration includes various sub-fields such as mineral resource exploration, geological disaster prediction, oil and gas exploration, metal ore exploration, dam leakage detection, hydrogeological survey, underground pipeline detection, etc. Each sub-field usually needs to calculate the three-dimensional spatial magnetic field distribution of a geological area to provide technical support.

[0003] That is, each sub-field needs to consider how to obtain the three-dimensional spatial magnetic field distribution of a geological area in actual production. Summary of the Invention

[0004] The present invention provides a forward simulation method, device, electronic device and storage medium for magneto-resistivity method, which is used to solve the technical problem of how to obtain the three-dimensional spatial magnetic field distribution of a geological area, and realizes the calculation of the three-dimensional spatial magnetic field distribution of a geological area through the forward simulation method of magneto-resistivity.

[0005] The present invention provides a forward simulation method for magneto-resistivity method, including the following steps: Obtain the modeling-related data of the geological area to be measured and the measured values of the resistivity at each position of the geological area to be measured; Based on the modeling-related data, establish a three-dimensional structure model of the geological area to be measured; Based on the measured values of the resistivity, assign the resistivity to the three-dimensional structure model to obtain a three-dimensional resistivity model of the geological area to be measured; Determine the three-dimensional spatial magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model.

[0006] According to the forward simulation method for magneto-resistivity method provided by the present invention, the step of assigning the resistivity to the three-dimensional structure model based on the measured values of the resistivity includes: Divide the three-dimensional structure model into multiple polyhedron units; Based on the measured values of the resistivity, assign the resistivity to each of the polyhedron units.

[0007] According to the forward simulation method for magneto-resistivity method provided by the present invention, the step of assigning the resistivity to each of the polyhedron units based on the measured values of the resistivity includes: Determine the target area in the three-dimensional structure model; Divide each polyhedron unit in each target area into multiple polyhedron sub-units; Assign resistivity values to each of the polyhedral sub-units and each of the polyhedral units outside the target area based on the measured values of the resistivity.

[0008] According to a forward simulation method of magneto-resistivity method provided by the present invention, the dividing each polyhedral unit within each target area into a plurality of polyhedral sub-units includes: Dividing each polyhedral unit within the target area into a plurality of the polyhedral sub-units by an adaptive mesh refinement technique.

[0009] According to a forward simulation method of magneto-resistivity method provided by the present invention, the determining the three-dimensional spatial magnetic field distribution of the geological area to be measured based on the three-dimensional resistivity model includes: Establish an integral expression between the magnetic field intensity and the resistivity, and transform the integral expression into a discrete equation set by the finite volume method; Substitute the resistivity values of the three-dimensional resistivity model into the discrete equation set, and solve the discrete equation set to obtain the three-dimensional spatial magnetic field distribution.

[0010] According to a forward simulation method of magneto-resistivity method provided by the present invention, the solving the discrete equation set includes: Solve the discrete equation set by the preconditioned stabilized bi-conjugate gradient method.

[0011] The present invention also provides a forward simulation device for magneto-resistivity method, including the following modules: An acquisition module, configured to acquire the modeling-related data of the geological area to be measured and the measured values of the resistivity at each position of the geological area to be measured; A modeling module, configured to establish a three-dimensional structure model of the geological area to be measured based on the modeling-related data; An assignment module, configured to assign resistivity values to the three-dimensional structure model based on the measured values of the resistivity to obtain a three-dimensional resistivity model of the geological area to be measured; A determination module, configured to determine the three-dimensional spatial magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model.

[0012] The present invention also provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, where when the processor executes the program, it implements the forward simulation method of magneto-resistivity method as described in any one of the above.

[0013] The present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, it implements the forward simulation method of magneto-resistivity method as described in any one of the above.

[0014] The present invention also provides a computer program product, including a computer program, which when executed by a processor, implements the forward simulation method of magneto-resistivity method as described in any one of the above.

[0015] The forward simulation method, device, electronic device and storage medium of magneto-resistivity method provided by the present invention The forward simulation method of magneto-resistivity method provided by the present invention establishes a three-dimensional structure model of a geological area to be measured and a three-dimensional resistivity model based on the three-dimensional structure model, determines the three-dimensional spatial magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model, and realizes the calculation of the three-dimensional spatial magnetic field distribution of the geological area to be measured through the forward simulation method of magneto-resistivity. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following will briefly introduce the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings in the following description are some embodiments of the present invention. For those of ordinary skill in the art, other drawings can also be obtained based on these drawings without creative efforts.

[0018] Figure 1 is one of the flow schematic diagrams of the forward simulation method of magneto-resistivity method provided by the present invention.

[0019] Figure 2 is a schematic diagram of the three-dimensional structure model provided by the present invention.

[0020] Figure 3 is a schematic diagram of the three-dimensional structure model after coarse grid division provided by the present invention.

[0021] Figure 4 is a schematic diagram of the three-dimensional structure model after fine grid division provided by the present invention.

[0022] Figure 5 is another schematic diagram of the three-dimensional structure model after fine grid division provided by the present invention.

[0023] Figure 6 is a schematic diagram of the three-dimensional resistivity model provided by the present invention.

[0024] Figure 7 is a schematic diagram of the three-dimensional spatial magnetic field distribution provided by the present invention.

[0025] Figure 8 is a schematic diagram of the structure of the forward simulation device of magneto-resistivity method provided by the present invention.

[0026] Figure 9 is a schematic diagram of the structure of the electronic device provided by the present invention.

[0027] Description of the reference numerals in the drawings: 10 - formation; 20 - dam body; 30 - water body; 40 - leakage channel; 50 - area corresponding to the power supply electrode; 60 - project area. Detailed implementation manners

[0028] In order to make the objectives, technical solutions and advantages of the present invention clearer, the technical solutions in the present invention will be clearly and completely described below with reference to the drawings in the present invention. Apparently, the described embodiments are some but not all of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments in the present invention without making creative efforts shall fall within the protection scope of the present invention.

[0029] It should be noted that in the description of the present invention, the terms "comprise", "include" or any other variants thereof are intended to cover non - exclusive inclusion, so that a process, method, article or device including a series of elements not only includes those elements but also includes other elements not explicitly listed, or further includes elements inherent to such a process, method, article or device. Without further limitation, an element defined by the statement "comprising one..." does not exclude the existence of another identical element in the process, method, article or device including the said element. The orientation or positional relationship indicated by terms such as "upper", "lower", etc. is based on the orientation or positional relationship shown in the drawings, and is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus cannot be construed as a limitation to the present invention. Unless otherwise clearly specified and defined, the terms "mounted", "connected" and "coupled" should be understood in a broad sense. For example, it may be a fixed connection, a detachable connection or an integral connection; it may be a mechanical connection or an electrical connection; it may be directly connected or indirectly connected through an intermediate medium, and may be the internal communication of two elements. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0030] The terms "first", "second", etc. in the present invention are used to distinguish similar objects and are not 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 can be implemented in an order other than those illustrated or described herein, and the objects distinguished by "first", "second", etc. are generally of the same type, and do not limit the number of objects. For example, the first object may be one or multiple. In addition, "and / or" indicates at least one of the connected objects, and the character " / " generally indicates an "or" relationship between the associated objects before and after.

[0031] The following will combine with Figures 1 - 8 to describe the forward simulation method, device, electronic device, and storage medium of the magneto-resistivity method provided by the present invention, aiming to calculate the three-dimensional space magnetic field distribution of the geological area to be measured through the magneto-resistivity forward simulation method.

[0032] Figure 1 is one of the flow schematic diagrams of the forward simulation method of the magneto-resistivity method provided by the present invention. As Figure 1 shown, it includes but is not limited to the following steps: Step S1, obtain the modeling-related data of the geological area to be measured and the measured values of the resistivity at each position of the geological area to be measured.

[0033] In some embodiments, the geological area to be measured can be various geological areas, such as the area where the dam is located. As Figure 2 shown, the area where the dam is located includes components such as air, formation 10, dam body 20, water body 30, and leakage channel 40.

[0034] The modeling-related data can include point cloud data and Building Information Modeling (BIM) data. The point cloud data is used to accurately reflect the actual shape and structure of the geological area to be measured, and the BIM data includes various component information of the geological area to be measured. The point cloud data can be obtained by scanning the geological area to be measured with a three-dimensional scanning device, and the BIM data can be obtained from engineering drawing materials and data obtained from geological exploration.

[0035] In some embodiments, the measured values of the resistivity at each position of the geological area to be measured can be obtained through geological exploration borehole experiments.

[0036] Step S2, establish a three-dimensional structure model of the geological area to be measured based on the modeling-related data.

[0037] In some embodiments, the three-dimensional structure model of the geological area to be measured can be established through three-dimensional geological modeling software and engineering structure BIM modeling methods. Exemplarily, the three-dimensional geological modeling software can be Revit software. Revit is a BIM-based software with the following characteristics: using parametric components, users can adjust the attributes of the components, such as materials and dimensions, by modifying parameters, thus quickly generating new designs; most components in Revit are related. When one component changes, the related components will be automatically updated, reducing the tediousness of manual adjustment; Revit supports multi-person collaboration, and design files can be managed on a central server, and updated in real time to ensure the consistency of the design.

[0038] Among them, in Revit, the "Insert" tab can be selected, and then the "Point Cloud" command can be clicked to import the point cloud data obtained in step S1 into the project.

[0039] Among them, in step S2, a three-dimensional structural model of the engineering area 60 (including the dam body 20, the water body 30, and the geological body) needs to be established first, and then the three-dimensional structural model of the engineering area is enlarged to obtain the three-dimensional structural model of the geological area to be measured.

[0040] Exemplarily, the three-dimensional structural model of the area where the dam is located established in step S2 can be as Figure 2 shown, the file format is STL. In the figure, the light green part is the water body 30, the blue and red parts are the dam body 20, the gray part is the stratum 10, and the yellow part is the leakage channel 40.

[0041] Step S3: Assign the resistivity values to the three-dimensional structural model based on the measured resistivity values to obtain the three-dimensional resistivity model of the geological area to be measured.

[0042] It can be understood that each position in the geological area to be measured has a corresponding position in the three-dimensional resistivity model. Step S3 specifically assigns the measured resistivity values of each position in the geological area to be measured to the corresponding positions in the three-dimensional structural model.

[0043] Step S4: Determine the three-dimensional spatial magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model.

[0044] In some embodiments, the three-dimensional spatial magnetic field distribution of the geological area to be measured can be determined according to the three-dimensional resistivity model by the finite difference method (such as the staggered grid finite difference method) and the finite element method.

[0045] Thus, steps S1 - S4 actually constitute a magneto-resistivity forward modeling method, that is, the three-dimensional spatial magnetic field distribution of the geological area to be measured can be calculated by the magneto-resistivity forward modeling method. The magneto-resistivity method inputs direct current or low-frequency alternating current between two points in the geological area to be measured through the power supply electrodes, measures the magnetic field excited by the underground conductors on the ground surface, and solves engineering geological problems such as dam leakage detection, hydrogeological survey, underground pipeline detection, and metal ore exploration based on the variation law of the magnetic field.

[0046] As can be seen from the above, the present invention establishes a three-dimensional structural model of the geological area to be measured and a three-dimensional resistivity model of the geological area to be measured based on the three-dimensional structural model, and determines the three-dimensional spatial magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model, which is equivalent to calculating the three-dimensional spatial magnetic field distribution of the geological area to be measured by the magneto-resistivity forward modeling method, realizing the calculation of the three-dimensional spatial magnetic field distribution of the geological area to be measured, and can provide technical support for the emergency rescue of dam leakage risks.

[0047] In step S3, assigning resistivity values to the three-dimensional structure model based on the measured resistivity values may include: dividing the three-dimensional structure model of the geological area to be measured into multiple polyhedron units; assigning resistivity values to each polyhedron unit based on the measured resistivity values.

[0048] Among them, the three-dimensional structure model of the geological area to be measured can be meshed and evenly divided into multiple polyhedron units; for example, the polyhedron units can include hexahedron units and tetrahedron units, and most of the area of the three-dimensional structure model is divided into multiple hexahedron units, and a small part of the transition area of the three-dimensional structure model is divided into multiple tetrahedron units, resulting in a three-dimensional structure model as Figure 3 shown.

[0049] It can be understood that each position in the geological area to be measured corresponds to a polyhedron unit in the three-dimensional structure model, so that the measured resistivity values of each position can be assigned to the corresponding polyhedron unit.

[0050] In this way, a three-dimensional resistivity model can be established by assigning resistivity values to the polyhedron units, and it can also ensure that the resistivity distribution of the three-dimensional resistivity model is consistent with the geological area to be measured, thereby improving the calculation accuracy of the three-dimensional space magnetic field distribution.

[0051] It can be understood that dividing the three-dimensional structure model of the geological area to be measured into multiple polyhedron units is equivalent to performing a coarse mesh division on the three-dimensional structure model. Each obtained polyhedron unit is relatively large and the total quantity is small. After assignment, the resistivity data volume of the three-dimensional resistivity model is small, and the calculation result accuracy is low.

[0052] Generally, when calculating by the magneto-resistivity forward simulation method, the finer the mesh division of the three-dimensional structure model and the larger the resistivity data volume of the three-dimensional resistivity model after assignment, the higher the calculation result accuracy, but the corresponding calculation efficiency is lower. In order to improve the accuracy of the calculation results, in some embodiments, each polyhedron unit of the three-dimensional structure model can be further divided into multiple polyhedron sub-units. For example, each hexahedron unit can be further divided into multiple hexahedron sub-units, and each tetrahedron unit can be further divided into multiple tetrahedron sub-units, but this will result in an excessive resistivity data volume and reduce the calculation efficiency.

[0053] Therefore, in some embodiments, assigning resistivity values to each polyhedron unit based on the measured resistivity values may further include: determining the target area in the three-dimensional structure model; dividing each polyhedron unit in each target area into multiple polyhedron sub-units; assigning resistivity values to each polyhedron sub-unit and each polyhedron unit outside the target area based on the measured resistivity values.

[0054] Among them, the target area is the area of key concern, such as the area 50 corresponding to the power supply electrode and the engineering area 60 (including the water body 30, the dam body 20 and the geological body). The target area can be determined in the way of [Xmin,Xmax;Ymin,Ymax;Zmin,Zmax;β], where Xmin and Xmax are the two boundary coordinates of the target area on the X-axis, Ymin and Ymax are the two boundary coordinates of the target area on the Y-axis, Zmin and Zmax are the two boundary coordinates of the target area on the Z-axis, and β is the multiple of grid refinement based on the coarse grid.

[0055] After determining the target area, each polyhedron unit in each target area is divided into multiple polyhedron sub-units, and the polyhedron units outside the target area are not subdivided. Exemplarily, a three-dimensional structure model of the area where the dam is located as shown in Figure 4 and Figure 5 can be obtained.

[0056] It can be understood that each position in the geological area to be measured corresponds to a polyhedron unit or a polyhedron sub-unit in the three-dimensional structure model. Thus, the measured value of the resistivity at each position can be assigned to the corresponding polyhedron unit or polyhedron sub-unit, and the resistivity values of each polyhedron unit and polyhedron sub-unit are constants. Exemplarily, after conductivity assignment, a three-dimensional resistivity model as shown in Figure 6 can be obtained, where different colors in the figure represent different resistivities.

[0057] Since the present invention can further refine the grid only for the areas of key concern, it can not only improve the calculation accuracy of the areas of key concern, but also avoid excessive reduction of the calculation efficiency.

[0058] In some embodiments, in step S3, dividing each polyhedron unit in each target area into multiple polyhedron sub-units may further include: dividing each polyhedron unit in the target area into multiple polyhedron sub-units by using the adaptive mesh refinement (AMR) technology.

[0059] The AMR technology can use a fine grid in the target area while maintaining a coarser grid in other areas, which can effectively reduce the amount of calculation and memory requirements, thereby improving the calculation efficiency while maintaining the calculation accuracy; and since the physical properties of the geological body in the engineering area may vary significantly over a large range, and the AMR technology allows local refinement of the grid in key areas, it can flexibly adapt to undulating terrain, complex structures, and multi-scale detection of deep and shallow structures; and since the AMR technology uses structured grids or quasi-structured grids, it has better orthogonality in the key refinement areas, making the gradient calculation more stable and accurate, and effectively reducing error accumulation.

[0060] As mentioned above, the three-dimensional spatial magnetic field distribution of the geological area to be measured can be determined according to the three-dimensional resistivity model by the finite difference method (such as the staggered grid finite difference method) and the finite element method. The finite difference method has high calculation efficiency, but its accuracy is relatively low, and it is only applicable to the electromagnetic field simulation of flat terrain and simple geometric shapes. The finite element method uses unstructured grids to better simulate complex terrain, and its accuracy is high. However, the simulation requires a large amount of unknown solutions to be calculated and stored, and the calculation efficiency is low. The forward simulation is the core technical part of the inversion interpretation imaging and the basis for the three-dimensional fine inversion of magneto-resistivity method data. The forward simulation methods based on the finite difference method or the finite element method are difficult to balance the calculation efficiency and calculation accuracy, and it is difficult to meet the requirements of three-dimensional fine simulation of complex structures at the engineering scale and fast inversion imaging.

[0061] In order to balance the calculation efficiency and calculation accuracy of the three-dimensional magneto-resistivity method for fine simulation of complex structures with undulating terrain, in some embodiments, step S4 may include: establishing an integral expression between the magnetic field intensity and the resistivity, and converting the integral expression into a discrete equation set by the finite volume method; substituting the resistivity value of the three-dimensional resistivity model into the discrete equation set and solving the discrete equation set to obtain the three-dimensional spatial magnetic field distribution.

[0062] Among them, the integral expression between the magnetic field intensity and the resistivity may be: ; In the formula, r is an arbitrary point on the three-dimensional resistivity model, is the magnetic field intensity at point r, is the magnetic permeability, is the position of the current density element, JS represents the power supply current density (also known as the primary current density), is the reciprocal of the resistivity (i.e., the conductivity), is the Hamiltonian operator, is the electric field scalar potential, is the derivative with respect to for derivation.

[0063] Among them, the finite volume method is based on the conservation equation in integral form rather than the differential equation, and focuses on constructing discrete equations from a physical perspective. Each discrete equation is an expression of the conservation of a certain physical quantity on a finite-sized volume. The derivation process has clear physical concepts, the coefficients of the discrete equations have certain physical meanings, and it can ensure that the discrete equations have conservation characteristics. By converting the integral expression into a discrete equation set by the finite volume method, the calculation efficiency and calculation accuracy of the three-dimensional magneto-resistivity method for fine simulation of complex structures with undulating terrain can be balanced.

[0064] In some embodiments, in step S4, discrete equations can be solved by software such as COMSOL using methods such as the Biconjugate Gradient Stabilized Method (BICGSTAB). In some embodiments, in step S4, solving the discrete equations can also include: solving the discrete equations by the Preconditioned BICGSTAB (PBICGSTAB).

[0065] PBICGSTAB is an iterative algorithm for solving large sparse linear equations. It adds a preprocessing step on the basis of the BICGSTAB algorithm, which can reduce the condition number of the coefficient matrix, thereby accelerating the convergence speed. Preprocessing usually uses incomplete LU decomposition or other preprocessing techniques, such as incomplete diagonal Cholesky decomposition (DIC) and incomplete LU decomposition (DILU). These preprocessing techniques can effectively improve the condition number of the matrix, making the iterative process more efficient and more suitable for solving asymmetric and sparse matrices, especially suitable for complex three-dimensional electromagnetic field solving scenarios, and can improve convergence stability.

[0066] Finally, after successively establishing a three-dimensional structure model of the geological area to be measured, performing coarse grid division on the three-dimensional structure model, performing fine grid division through the AMR technique, establishing a three-dimensional resistivity model, realizing magnetotelluric forward simulation by the finite volume method, and solving the discrete equations by the PBICGSTAB method, the three-dimensional spatial magnetic field distribution of the geological area to be measured as shown in Figure 7 can be obtained. Different colors in the figure represent different magnetic field intensities. Considering that the positions of the power supply electrodes corresponding in the three-dimensional structure model are different, the obtained three-dimensional spatial magnetic field distributions are also different. In practice, it is necessary to change the positions of the power supply electrodes corresponding in the three-dimensional structure model multiple times to obtain a variety of different three-dimensional spatial magnetic field distributions to reveal the magnetic field variation law of the geological area to be measured.

[0067] Furthermore, the present invention can also develop a visualization module for the entire simulation calculation process, supporting multi-information visualization interactive analysis of BIM models, grid meshing, calculation results, and other three-dimensional data of scientific calculations.

[0068] Figure 8 is one of the schematic diagrams of the magnetotelluric forward simulation device provided by the present invention, as shown in Figure 7 and includes but is not limited to the following modules: An acquisition module for acquiring the modeling-related data of the geological area to be measured and the measured values of the resistivity at each position of the geological area to be measured; A modeling module for establishing a three-dimensional structure model of the geological area to be measured based on the modeling-related data; An assignment module, configured to assign resistivity values to the three-dimensional structure model based on the measured resistivity values, so as to obtain a three-dimensional resistivity model of the geological area to be measured; A determination module, configured to determine the three-dimensional spatial magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model.

[0069] In some embodiments, the assignment module may specifically be configured to: Divide the three-dimensional structure model into multiple polyhedron units; Assign resistivity values to each polyhedron unit based on the measured resistivity values.

[0070] In some embodiments, the assignment module may also specifically be configured to: Determine the target area in the three-dimensional structure model; Divide each polyhedron unit within each target area into multiple polyhedron sub-units; Assign resistivity values to each polyhedron sub-unit within each target area and each polyhedron unit outside the target area based on the measured resistivity values.

[0071] In some embodiments, the assignment module may also specifically be configured to: Divide each polyhedron unit within the target area into multiple polyhedron sub-units by using an adaptive grid encryption technique.

[0072] In some embodiments, the determination module may specifically be configured to: Establish an integral expression between the magnetic field intensity and the resistivity, and convert the integral expression into a discrete equation set by using the finite volume method; Substitute the resistivity values of the three-dimensional resistivity model into the discrete equation set, and solve the discrete equation set to obtain the three-dimensional spatial magnetic field distribution.

[0073] In some embodiments, the determination module may also specifically be configured to: Solve the discrete equation set by using the preconditioned stabilized bi-conjugate gradient method.

[0074] It should be noted that when the forward simulation device of the magneto-resistivity method provided by the present invention is specifically operating, it can execute the magneto-resistivity method forward simulation method described in any of the above embodiments, and details thereof are not elaborated in this embodiment.

[0075] Figure 9 is a schematic structural diagram of an electronic device provided by the present invention, as Figure 8As shown in the figure, the electronic device may include: a processor, a communication interface, a memory, and a communication bus. Among them, the processor, the communication interface, and the memory complete communication with each other through the communication bus. The processor can call the logical instructions in the memory to execute the forward modeling simulation method of magneto-resistivity method, and this method includes: obtaining the modeling-related data of the geological area to be measured and the measured values of the resistivity at various positions in the geological area to be measured; establishing a three-dimensional structure model of the geological area to be measured based on the modeling-related data; assigning resistivity values to the three-dimensional structure model based on the measured values of the resistivity to obtain a three-dimensional resistivity model of the geological area to be measured; determining the three-dimensional spatial magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model.

[0076] In addition, when the logical instructions in the above-mentioned memory are implemented in the form of software functional units and sold or used as an independent product, they can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention. And the aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

[0077] On the other hand, the present invention also provides a computer program product. The computer program product includes a computer program stored on a non-transitory computer-readable storage medium. The computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the forward modeling simulation method of magneto-resistivity method provided in the above-mentioned various embodiments. This method includes: obtaining the modeling-related data of the geological area to be measured and the measured values of the resistivity at various positions in the geological area to be measured; establishing a three-dimensional structure model of the geological area to be measured based on the modeling-related data; assigning resistivity values to the three-dimensional structure model based on the measured values of the resistivity to obtain a three-dimensional resistivity model of the geological area to be measured; determining the three-dimensional spatial magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model.

[0078] In another aspect, the present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, it implements the forward simulation method of magnetoresistivity method provided in the above embodiments. The method includes: obtaining the modeling-related data of the geological area to be measured and the measured values of the resistivity at each position of the geological area to be measured; establishing a three-dimensional structure model of the geological area to be measured based on the modeling-related data; assigning resistivity values to the three-dimensional structure model based on the measured values of the resistivity to obtain a three-dimensional resistivity model of the geological area to be measured; and determining the three-dimensional spatial magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model.

[0079] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art can understand and implement it without creative effort.

[0080] Through the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on such an understanding, the above technical solutions, in essence, or the part that contributes to the prior art, can be embodied in the form of a software product. The computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0081] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements for some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A forward modeling method of magnetoresistivity method, characterized in that: include: Acquire modeling related data of the geological area to be measured and measured values ​​of resistivity at various locations in the geological area to be measured; Establishing a three-dimensional structural model of the geological area to be measured based on the modeling related data; Assigning a value to the resistivity of the three-dimensional structural model based on the measured value of the resistivity to obtain a three-dimensional resistivity model of the geological area to be measured; The three-dimensional magnetic field distribution of the geological area to be measured is determined according to the three-dimensional resistivity model.

2. The magnetoresistivity forward modeling method according to claim 1, characterized in that: The assigning a value to the resistivity of the three-dimensional structural model based on the measured value of the resistivity includes: Dividing the three-dimensional structural model into a plurality of polyhedral units; The resistivity of each of the polyhedral units is assigned a value based on the measured value of the resistivity.

3. The magnetoresistivity forward modeling method according to claim 2, characterized in that: The assigning of the resistivity of each of the polyhedral units based on the measured value of the resistivity includes: determining a target area in the three-dimensional structural model; Dividing each of the polyhedral units in each of the target areas into a plurality of polyhedral sub-units; The resistivity of each of the polyhedral sub-units and each of the polyhedral units outside the target area is assigned a value based on the measured value of the resistivity.

4. The magnetoresistivity forward modeling method according to claim 3, characterized in that: The step of dividing each of the polyhedral units in each of the target areas into a plurality of polyhedral sub-units comprises: Each of the polyhedral units in the target area is divided into a plurality of polyhedral sub-units by using an adaptive mesh encryption technology.

5. The magnetoresistivity forward modeling method according to claim 1, characterized in that: Determining the three-dimensional magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model includes: Establishing an integral expression between the magnetic field intensity and the resistivity, and converting the integral expression into a discrete equation group by a finite volume method; The resistivity value of the three-dimensional resistivity model is substituted into the discrete equation group, and the discrete equation group is solved to obtain the three-dimensional spatial magnetic field distribution.

6. The magnetoresistivity forward modeling method according to claim 5, characterized in that: The step of solving the discrete equations comprises: The discrete system of equations is solved by a preconditioned stable biconjugate gradient method.

7. A forward modeling device for magnetoresistivity method, characterized in that: include: An acquisition module, used to acquire modeling-related data of the geological area to be measured and measured values ​​of resistivity at various locations in the geological area to be measured; A modeling module, used for establishing a three-dimensional structural model of the geological area to be measured based on the modeling-related data; An assignment module, used for assigning the resistivity of the three-dimensional structural model based on the measured value of the resistivity to obtain a three-dimensional resistivity model of the geological area to be measured; The determination module is used to determine the three-dimensional magnetic field distribution of the geological area to be measured according to the three-dimensional resistivity model.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the magnetoresistivity forward modeling method according to any one of claims 1 to 6 is implemented.

9. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the magnetoresistivity forward modeling method according to any one of claims 1 to 6 is implemented.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the magnetoresistivity forward modeling method according to any one of claims 1 to 6 is implemented.

Images