CN113109876B - Interpretation method and device based on three-dimensional magneto-resistivity method tensor observation mode - Google Patents
Interpretation method and device based on three-dimensional magneto-resistivity method tensor observation mode Download PDFInfo
- Publication number
- CN113109876B CN113109876B CN202110199433.4A CN202110199433A CN113109876B CN 113109876 B CN113109876 B CN 113109876B CN 202110199433 A CN202110199433 A CN 202110199433A CN 113109876 B CN113109876 B CN 113109876B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- abnormal
- tensor
- power supply
- electric dipole
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Abstract
The invention provides a three-dimensional magneto-resistivity method tensor observation mode-based interpretation method and device. Wherein, the method comprises the following steps: carrying out tensor measurement on a target magnetic field under the condition that two electric dipole sources which are perpendicular to each other are respectively powered to obtain an actually measured magnetic field tensor; forward calculation is carried out by utilizing a three-dimensional magneto-resistivity method, and predicted magnetic field tensors of two electric dipole sources in the power supply direction are obtained based on a finite difference numerical simulation method; obtaining total amplitude according to the amplitude absolute values of all magnetic field components in the actually measured magnetic field tensor and the predicted magnetic field tensor, and determining abnormal weight factors of the two electric dipole sources based on the total amplitude; obtaining an abnormal standardized magnetic field according to the ratio of the actually measured magnetic field tensor to the predicted magnetic field tensor; and obtaining a comprehensive abnormal standardized magnetic field according to the abnormal weight factor and the abnormal standardized magnetic field, and further determining the plane position of the abnormal body. By adopting the method disclosed by the invention, the resolution capability of the abnormal body in the complex model can be enhanced, and the interpretation precision and efficiency of the three-dimensional magneto-resistivity method are improved.
Description
Technical Field
The invention relates to the technical field of exploration geophysics, in particular to an interpretation method and device based on a three-dimensional magneto-resistivity method tensor observation mode. In addition, an electronic device and a non-transitory computer readable storage medium are also related.
Background
The magneto-resistivity method is a geophysical electromagnetic prospecting method which measures the magnetic field excited by non-inductive (direct current or low-frequency alternating current) current supplied by manpower between two points on the ground surface and solves engineering and geological problems of underground water flow paths, underground pipelines, prospecting and the like according to the change rule of the magnetic field. However, the geological conditions of actual engineering are often very complex, and the geometric shape of the detected geological abnormal object is mostly irregular, such as a staggered abnormal body. In general, when a geological abnormal body is detected by using a magneto-resistivity method, a unidirectional electric dipole source is adopted for supplying power, current is mainly concentrated in the low-resistance abnormal body parallel to a power supply side, and the current in the low-resistance abnormal body perpendicular to the power supply direction is relatively low, so that the resolving power of the abnormal body perpendicular to the power supply direction is low. Therefore, it is difficult for the conventional magneto-resistivity method power supply mode to realize high-accuracy identification of an abnormal body perpendicular to the power supply direction.
Currently, a common three-dimensional magneto-resistivity interpretation method includes: conjugate gradient least squares, Occam, Gauss-Newton, etc. The common characteristics of the data interpretation methods are as follows: the operation efficiency is low, multiple forward calculations are required in the inversion process, a large amount of memory and a Jacobian matrix of calculation time are required for calculation, and the inversion time is usually calculated in days; the requirement on hardware is high, inversion calculation by using a common notebook computer is almost impossible, and a workstation is required, so that the field work is inconvenient. Therefore, how to realize fast and accurate interpretation of the abnormal plane position by changing the observation mode of the conventional magneto-resistivity method becomes a technical problem to be solved by those skilled in the art.
Disclosure of Invention
Therefore, the invention provides an interpretation method and device based on a three-dimensional magneto-resistivity method tensor observation mode, and aims to solve the problems that in the prior art, the operation efficiency is low when an abnormal body perpendicular to the power supply direction is identified with high precision, and the requirement on hardware is high.
The invention provides an interpretation method based on a three-dimensional magneto-resistivity method tensor observation mode, which comprises the following steps of: tensor measurement is carried out on a target magnetic field under the condition that two electric dipole sources which are perpendicular to each other are respectively powered, and an actually measured magnetic field tensor is obtained;
performing forward calculation by using a three-dimensional magneto-resistivity method, and respectively obtaining earth surface magnetic fields of two electric dipole sources in the power supply direction based on a finite difference numerical simulation method to obtain a predicted magnetic field tensor;
acquiring the total amplitude of each magnetic field component in a measurement area according to the absolute value of the amplitude of each magnetic field component in the actually measured magnetic field tensor and the predicted magnetic field tensor, and determining the abnormal weight factor of the two electric dipole sources in the corresponding power supply direction based on the total amplitude;
obtaining an abnormal standardized magnetic field according to the ratio of the actually measured magnetic field tensor to the predicted magnetic field tensor;
obtaining a comprehensive abnormal standardized magnetic field corresponding to each place of the earth surface according to the abnormal weight factor and the abnormal standardized magnetic field; and determining the plane position of the abnormal body according to the comprehensive abnormal standardized magnetic field.
Further, the expression corresponding to the actually measured magnetic field tensor is as follows:
wherein Hxx、HxyA magnetic field component corresponding to the X, Y power supply direction measured on the X-axis power supply ground surface is provided for the electric dipole source central line; hyx、HyyAnd a magnetic field component corresponding to the X, Y power supply direction measured on the power supply ground surface of the Y-axis for the center line of the electric dipole source.
Further, the expression corresponding to the predicted magnetic field tensor is as follows:
wherein Hxx′、Hxy' a magnetic field component corresponding to X, Y power supply direction predicted by the power supply surface of the X axis of the center line of the electric dipole source; hyx′、Hyy' the predicted X, Y power supply direction corresponding magnetic field component of the Y-axis power supply surface for the center line of the electric dipole source.
Further, the abnormal weight factor includes a first abnormal weight factor ω X corresponding to the power supply in the X-axis direction and a second abnormal weight factor ω Y corresponding to the power supply in the Y-axis direction;
the calculation formula corresponding to the first abnormality weight ω x is:
therein, sigmaΩ|Hxx|、∑Ω|HxyI is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the actually measured magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction; sigmaΩ|Hxx′|、∑Ω|Hxy' I is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the predicted magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction;
the calculation formula corresponding to the second abnormality weight ω y is as follows:
therein, sigmaΩ|Hyy|、∑Ω|HyxI is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the actually measured magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction; sigmaΩ|Hyy′|、∑Ω|Hyx' l is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the predicted magnetic field tensor in the measurement region Ω when power is supplied in the X-axis direction, respectively.
Further, the calculation formula corresponding to the abnormal normalized magnetic field is as follows:
wherein HxsFor abnormal standardized magnetic field, H, of electric dipole source central line when power is supplied on X axisysAn abnormal standardized magnetic field is generated when the electric dipole source central line supplies power on the Y axis; the abnormal standardized magnetic field is used for determining the relative strength of the magnetic field or the relative strength of current distribution in the measurement region under the homogeneous condition;
abnormal normalized magnetic field H at X-axis power supplyxsCorresponding calculation formula and abnormal standardized magnetic field H in power supply of Y axisysThe corresponding calculation formulas are respectively as follows:
wherein HxyAnd HxyRespectively measuring and predicting the magnetic field component H corresponding to the Y power supply direction of the electric dipole source central line on the X-axis power supply ground surfaceyxAnd HyxThe' are respectively the magnetic field components corresponding to the X power supply direction actually measured and predicted by the electric dipole source central line on the Y-axis power supply ground surface.
Further, the comprehensive abnormal normalized magnetic field is a ratio of a sum of products of the abnormal weight factors and the abnormal normalized magnetic field in each power supply direction to a sum of the abnormal weight factors in each power supply direction, and a specific calculation formula is as follows:
wherein, ω X is a first abnormal weight factor corresponding to the power supply in the X-axis direction, and ω Y is a second abnormal weight factor corresponding to the power supply in the Y-axis directionA seed; hxsFor abnormal standardized magnetic field, H, of electric dipole source central line when power is supplied on X axisysAnd (4) an abnormal standardized magnetic field is generated when the central line of the electric dipole source supplies power on the Y axis.
Further, the larger the first abnormality weight ω x is, the larger the current of the measurement region parallel to the power supply direction is, and the smaller the first abnormality weight ω x is, the smaller the current of the measurement region parallel to the power supply direction is; the larger the second abnormality weight ω y, the larger the current of the measurement region parallel to the power supply direction, and the smaller the second abnormality weight ω y, the smaller the current of the measurement region parallel to the power supply direction.
The invention also provides an interpretation device based on the three-dimensional magneto-resistivity method tensor observation mode, which comprises the following components:
the actual measurement magnetic field tensor obtaining unit is used for carrying out tensor measurement on the target magnetic field under the condition that two electric dipole sources which are perpendicular to each other are respectively powered to obtain an actual measurement magnetic field tensor;
the predicted magnetic field tensor obtaining unit is used for performing forward calculation by using a three-dimensional magneto-resistivity method, and respectively obtaining earth surface magnetic fields of the two electric dipole sources in the power supply direction based on a finite difference numerical simulation method to obtain a predicted magnetic field tensor;
the abnormal weight factor determining unit is used for obtaining the total amplitude of each magnetic field component in the measurement area according to the absolute value of the amplitude of each magnetic field component in the actually measured magnetic field tensor and the predicted magnetic field tensor, and determining the abnormal weight factors of the two electric dipole sources in the corresponding power supply directions based on the total amplitude;
an abnormal normalized magnetic field obtaining unit configured to obtain an abnormal normalized magnetic field according to a ratio of the measured magnetic field tensor to the predicted magnetic field tensor;
the abnormal body position determining unit is used for obtaining a comprehensive abnormal standardized magnetic field corresponding to each position of the earth surface according to the abnormal weight factor and the abnormal standardized magnetic field; and determining the plane position of the abnormal body according to the comprehensive abnormal standardized magnetic field.
Further, the expression corresponding to the actually measured magnetic field tensor is as follows:
wherein Hxx、HxyA magnetic field component corresponding to the X, Y power supply direction measured on the X-axis power supply ground surface is provided for the electric dipole source central line; hyx、HyyAnd a magnetic field component corresponding to the X, Y power supply direction measured on the power supply ground surface of the Y-axis for the center line of the electric dipole source.
Further, the expression corresponding to the predicted magnetic field tensor is as follows:
wherein Hxx′、Hyx' respectively predicting X, Y power supply direction corresponding magnetic field components of the power supply surface of the X-axis power supply by the center line of the electric dipole source; hyx′、Hyy' the magnetic field components corresponding to X, Y power supply directions predicted by the power supply surface of the Y-axis of the center line of the electric dipole source are respectively.
Further, the abnormal weight factor includes a first abnormal weight factor ω X corresponding to the power supply in the X-axis direction and a second abnormal weight factor ω Y corresponding to the power supply in the Y-axis direction;
the calculation formula corresponding to the first abnormality weight ω x is:
therein, sigmaΩ|Hxx|、∑Ω|HxyI is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the actually measured magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction; sigmaΩ|Hxx′|、∑Ω|Hxy' I is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the predicted magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction;
the calculation formula corresponding to the second abnormality weight ω y is as follows:
therein, sigmaΩ|Hyy|、∑Ω|HyxI is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the actually measured magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction; sigmaΩ|Hyy′|、∑Ω|Hyx' l is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the predicted magnetic field tensor in the measurement region Ω when power is supplied in the X-axis direction, respectively.
Further, the calculation formula corresponding to the abnormal normalized magnetic field is as follows:
wherein HxsFor abnormal standardized magnetic field, H, of electric dipole source central line when power is supplied on X axisysAn abnormal standardized magnetic field is generated when the electric dipole source central line supplies power on the Y axis; the abnormal standardized magnetic field is used for determining the relative strength of the magnetic field or the relative strength of current distribution in the measurement region under the homogeneous condition;
abnormal normalized magnetic field H at X-axis power supplyxsCorresponding calculation formula and abnormal standardized magnetic field H in power supply of Y axisysThe corresponding calculation formulas are respectively as follows:
wherein HxyAnd HxyThe measured and predicted magnetic field component in the Y power supply direction is the magnetic field component H corresponding to the power supply direction of the electric dipole source central line on the X-axis power supply ground surfaceyxAnd HyxThe' is a magnetic field component corresponding to the X power supply direction which is actually measured and predicted on the power supply ground surface of the Y axis by the central line of the electric dipole source.
Further, the comprehensive abnormal normalized magnetic field is a ratio of a sum of products of the abnormal weight factors and the abnormal normalized magnetic field in each power supply direction to a sum of the abnormal weight factors in each power supply direction, and a specific calculation formula is as follows:
wherein, ω X is a first abnormal weight factor corresponding to the power supply in the X-axis direction, and ω Y is a second abnormal weight factor corresponding to the power supply in the Y-axis direction; hxsFor abnormal standardized magnetic field, H, of electric dipole source central line when power is supplied on X axisysAnd (4) an abnormal standardized magnetic field is generated when the central line of the electric dipole source supplies power on the Y axis.
Further, the larger the first abnormality weight ω x is, the larger the current of the measurement region parallel to the power supply direction is, and the smaller the first abnormality weight ω x is, the smaller the current of the measurement region parallel to the power supply direction is; the larger the second abnormality weight ω y, the larger the current of the measurement region parallel to the power supply direction, and the smaller the second abnormality weight ω y, the smaller the current of the measurement region parallel to the power supply direction.
The invention also provides an electronic device, which comprises a memory, a processor and a computer program stored in the memory and running on the processor, wherein the processor executes the program to realize the interpretation method based on the three-dimensional magneto-resistivity tensor observation mode.
The present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a three-dimensional magneto-resistivity tensor observation pattern based interpretation method as described in any one of the above.
According to the interpretation method based on the three-dimensional magneto-resistivity method tensor observation mode, disclosed by the invention, aiming at a three-dimensional complex geological model, the three-dimensional magneto-resistivity method tensor observation mode is used for measuring magnetic field data of a plurality of power supply directions, and an integrated abnormal standardized magnetic field is further obtained through abnormal weight factors solved through magnetic field amplitude, wherein the integrated abnormal standardized magnetic field represents the integrated strength of magnetic abnormality of each power supply direction, so that an abnormal body in the complex model is easier to determine, the resolution capability of the three-dimensional complex abnormal body is improved, the subsequent solution of the plane position of the three-dimensional complex abnormal body is facilitated, and the interpretation efficiency of the three-dimensional magneto-resistivity method is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings based on the drawings without creative efforts.
Fig. 1 is a schematic flowchart of an interpretation method based on a three-dimensional magneto-resistivity tensor observation mode according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a cross anomaly model provided by an embodiment of the present invention;
FIG. 3a is a schematic diagram of a method H according to an embodiment of the present inventionxx' a contour map of magnetic field strength;
FIG. 3b is a schematic diagram of a process H according to an embodiment of the present inventionxy' a contour map of magnetic field strength;
FIG. 4a is H provided in the embodiment of the present inventionxxA contour map of magnetic field strength;
FIG. 4b is a schematic diagram of a process H according to an embodiment of the present inventionxyA contour map of magnetic field strength;
FIG. 5 is a diagram illustrating the X-axis power supply normalization H according to an embodiment of the present inventionxsA contour map of magnetic field strength;
FIG. 6 is a comprehensive abnormal normalized magnetic field H contour plot provided by an embodiment of the present invention;
fig. 7 is a schematic structural diagram of a weighted anomaly interpretation method apparatus based on a three-dimensional magneto-resistivity tensor observation mode according to an embodiment of the present invention;
fig. 8 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The following describes an embodiment of the method for interpreting the three-dimensional magneto-resistivity tensor observation mode based on the invention in detail. As shown in fig. 1, which is a schematic flow chart of an interpretation method based on a three-dimensional magneto-resistivity tensor observation mode according to an embodiment of the present invention, the specific process includes the following steps:
step 101: and under the condition that two electric dipole sources which are perpendicular to each other are respectively powered, carrying out tensor measurement on the target magnetic field to obtain an actually measured magnetic field tensor.
Specifically, the actually measured magnetic field tensor comprises 4 magnetic field components which are X, Y-direction magnetic fields measured by the central line of the electric dipole source on the X-axis power supply ground surface respectively; the center line of the electric dipole source supplies X, Y directional magnetic fields measured at the earth's surface along the Y-axis.
In the implementation process, two orthogonal electric dipole sources can be used for supplying power respectively according to the actual situation on site, the tensor measurement is carried out on the target magnetic field on the ground, and the actually measured magnetic field tensor is HobsThe actually measured magnetic field tensor is expressed by the following corresponding expression:
wherein Hxx、HxyA magnetic field component corresponding to the X, Y power supply direction measured on the X-axis power supply ground surface is provided for the electric dipole source central line; hyx、HyyA magnetic field component corresponding to the X, Y power supply direction measured on the Y-axis power supply ground surface is provided for the electric dipole source central line; two directional electric dipole sources, two components measured in each direction,for a total of 4 components.
Step 102: forward calculation is carried out by utilizing a three-dimensional magneto-resistivity method, and earth surface magnetic fields of two electric dipole sources in the power supply direction are respectively obtained based on a finite difference numerical simulation method, so that a predicted magnetic field tensor is obtained.
Specifically, the predicted magnetic field tensor also includes 4 magnetic field components and corresponds to the 4 magnetic field components of the measured magnetic field tensor. In the implementation process, a three-dimensional magneto-resistivity method can be used for forward calculation, a measuring area is assumed to be a uniform medium, finite difference numerical simulation methods are adopted to respectively calculate electric dipole source earth surface magnetic fields in two directions, and the tensor of the predicted magnetic field is Hpre. The expression corresponding to the predicted magnetic field tensor is as follows:
wherein Hxx′、Hxy' a magnetic field component corresponding to X, Y power supply direction predicted by the power supply surface of the X axis of the center line of the electric dipole source; hyx′、Hyy' the predicted X, Y power supply direction corresponding magnetic field component of the Y-axis power supply surface for the center line of the electric dipole source.
The forward calculation formula is as follows:
the underground field established in the magneto-resistivity method can be treated as a stable current field, and maxwell's equations under the stable current field can be written as:
j＝σE (5)
in the formulas (3) to (6), E is an electric field intensity vector, H is a magnetic field intensity vector, j is a current density vector, sigma is the conductivity of the underground medium, and mu is the magnetic conductivity; in order to solve the final magnetic field strength H, firstly, the electric field strength E is required to be solved, and a scalar potential V is introduced to decompose the solving process by considering the vector characteristics of the electric field strength E:
i.e. first the scalar field of the potential V is solved and then the vector field H is further solved.
Step 103: and obtaining the total amplitude of each magnetic field component in the measurement area according to the absolute value of the amplitude of each magnetic field component in the actually measured magnetic field tensor and the predicted magnetic field tensor, and determining the abnormal weight factor of the two electric dipole sources in the corresponding power supply direction based on the total amplitude.
The abnormal weight factor is the ratio of the energy sum (modulus of amplitude) of the magnetic field component perpendicular to the power supply direction in the measured magnetic field tensor at each point of the measurement area to the energy sum of the magnetic field component parallel to the power supply direction in the measurement area to the energy sum of the magnetic field component in the corresponding direction in the predicted magnetic field tensor.
Specifically, the abnormality weight ω includes: a first abnormality weight factor ω X corresponding to the power supply in the X-axis direction and a second abnormality weight factor ω Y corresponding to the power supply in the Y-axis direction. The higher the first abnormal weight factor ω x is, the higher the current of the measurement region parallel to the power supply direction is, and the lower the first abnormal weight factor ω x is, the lower the current of the measurement region parallel to the power supply direction is; the larger the second abnormality weight ω y, the larger the current of the measurement region parallel to the power supply direction, and the smaller the second abnormality weight ω y, the smaller the current of the measurement region parallel to the power supply direction.
The calculation formula corresponding to the first abnormality weight ω x is:
therein, sigmaΩ|Hxx|、∑Ω|HxyI is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the actually measured magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction; sigmaΩ|Hxx′|、∑Ω|Hyx' l is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the predicted magnetic field tensor in the measurement region Ω when power is supplied in the X-axis direction, respectively.
The calculation formula corresponding to the second abnormality weight ω y is as follows:
therein, sigmaΩ|Hyy|、∑Ω|HyxI is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the actually measured magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction; sigmaΩ|Hyy′|、∑Ω|Hyx' I is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the predicted magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction; the parameters in ω y are the same as above; larger ω means larger current of the measurement area parallel to the power supply direction, and smaller ω means smaller current of the measurement area parallel to the power supply direction.
Step 104: and obtaining an abnormal standardized magnetic field according to the ratio of the actually measured magnetic field tensor to the predicted magnetic field tensor.
Specifically, the abnormal normalized magnetic field is a ratio of a magnetic field component perpendicular to the power supply direction in the actually measured magnetic field tensor to a magnetic field component in the corresponding direction in the predicted magnetic field tensor, and is a dimensionless numerical value.
In the embodiment of the present invention, an abnormal normalized magnetic field may be obtained according to a ratio of a measured magnetic field to a predicted magnetic field, where a calculation formula corresponding to the abnormal normalized magnetic field is:
wherein HxsFor abnormal standardized magnetic field, H, of electric dipole source central line when power is supplied on X axisysAn abnormal standardized magnetic field is generated when the electric dipole source central line supplies power on the Y axis; the abnormal standardized magnetic field is used for determining the relative strength of the magnetic field or the relative strength of current distribution in the measuring area under the homogeneous condition.
Abnormal normalized magnetic field H at X-axis power supplyxsCorresponding calculation formula and abnormal standardized magnetic field H during power supply of Y axisysThe corresponding calculation formulas are respectively as follows:
wherein HxyAnd HxyRespectively measuring and predicting the magnetic field component H corresponding to the Y power supply direction of the electric dipole source central line on the X-axis power supply ground surfaceyxAnd HyxThe' are respectively the magnetic field components corresponding to the X power supply direction actually measured and predicted by the electric dipole source central line on the Y-axis power supply ground surface. The abnormal standardized magnetic field can visually reflect the relative strength of the magnetic field or the relative strength of current distribution in a measurement area under a homogeneous condition.
Step 105: obtaining a comprehensive abnormal standardized magnetic field corresponding to each place of the earth surface according to the abnormal weight factor and the abnormal standardized magnetic field; and determining the plane position of the abnormal body according to the comprehensive abnormal standardized magnetic field.
Specifically, the comprehensive abnormal standardized magnetic field is the sum of products of the abnormal weight factors of all power supply directions and the abnormal standardized magnetic field, and the sum of the power supply direction weight factors is divided, so that the comprehensive strength of the magnetic abnormality of all the power supply directions is represented, and the resolution capability of the complex model to the abnormal body is improved.
In the embodiment of the present invention, the comprehensive abnormal normalized magnetic field H at each location on the surface can be obtained from the abnormal weight factor and the abnormal normalized magnetic field. The comprehensive abnormal standardized magnetic field is a ratio of the sum of products of the abnormal weight factors and the abnormal standardized magnetic field in each power supply direction to the sum of the abnormal weight factors in each power supply direction, and the specific calculation formula is as follows:
wherein, ω X is a first abnormal weight factor corresponding to the power supply in the X-axis direction, and ω Y is a second abnormal weight factor corresponding to the power supply in the Y-axis direction; hxsFor abnormal standardized magnetic field, H, of electric dipole source central line when power is supplied on X axisysAnd (4) an abnormal standardized magnetic field is generated when the central line of the electric dipole source supplies power on the Y axis. The higher the first abnormal weight factor ω x is, the higher the current of the measurement region parallel to the power supply direction is, and the lower the first abnormal weight factor ω x is, the lower the current of the measurement region parallel to the power supply direction is; the larger the second abnormality weight ω y, the larger the current of the measurement region parallel to the power supply direction, and the smaller the second abnormality weight ω y, the smaller the current of the measurement region parallel to the power supply direction. The comprehensive abnormal standardized magnetic field represents the comprehensive strength of the magnetic abnormality of each power supply direction, so that the resolution capability of an abnormal body in a complex model is improved.
It should be noted that the invention is improved on the basis of the traditional magneto-resistivity method, measures the magnetic field data of two mutually perpendicular power supply directions through a tensor observation mode, and improves the interpretation capability of the underground three-dimensional complex model, and the specific implementation process comprises the following steps: determining the position of a power supply, and measuring magnetic field data of a plurality of power supply directions; forward modeling calculation, namely respectively calculating magnetic fields generated by electric dipole sources in different power supply directions on the earth surface on the assumption that a measurement area is a uniform medium; according to the amplitude absolute values of the magnetic field components in the actually measured magnetic field tensor and the predicted magnetic field tensor, solving an abnormal weight factor and an abnormal standardized magnetic field of the power supply directions of the two electric dipole sources; and obtaining a comprehensive abnormal standardized magnetic field at each part of the earth surface according to the abnormal weight factor and the abnormal standardized magnetic field.
At one endIn a specific embodiment, the designed three-dimensional complex model may be the model shown in fig. 1, the distance between the power supply electrodes AB is 200 meters, the power supply current is 1 ampere, the staggered abnormal body is a combination of two rectangular bodies with the size of 100 × 40 × 10 meters, the two rectangular bodies form a cross shape, the center is located at the origin (0, 0, 0), the top buried depth is 20 meters, the resistivity of the surrounding rock is 100 ohm meters, and the resistivity of the staggered abnormal body is 10 ohm meters. First, the power supply electrode AB is positioned at coordinates (-100, 0, 0), (100, 0, 0), and the predicted magnetic field tensor H at the surface where the cross anomaly model is positioned at Z ═ 0 under the uniform medium condition is calculated, respectivelypreActually measured magnetic field tensor HobsAbnormal normalized magnetic field HxsAs shown in fig. 2, 3a, 3b, 4a and 4b, the total amplitude of each magnetic field component in the measurement region is used to calculate an abnormality weight factor ω x; then, the power supply electrode AB is positioned at coordinates (0, -100, 0) and (0, 100, 0), and the abnormality weight ω y is obtained in the same manner; finally, the comprehensive abnormal normalized magnetic field H at each location on the surface is obtained from the abnormal weight factors and the abnormal normalized magnetic field, and the planar position of the underground abnormal body can be reflected as shown in fig. 5.
By adopting the interpretation method based on the three-dimensional magneto-resistivity method tensor observation mode, disclosed by the embodiment of the invention, aiming at the three-dimensional complex geological model, the three-dimensional magneto-resistivity method tensor observation mode is used for measuring the magnetic field data of a plurality of power supply directions, and the comprehensive abnormal standardized magnetic field is obtained through the abnormal weight factors obtained through the magnetic field amplitude, represents the comprehensive strength of the magnetic abnormality of each power supply direction, so that the abnormal body in the complex model is easier to determine, the resolution capability of the three-dimensional complex abnormal body is improved, the subsequent solution of the plane position of the three-dimensional complex abnormal body is facilitated, and the interpretation efficiency and the accuracy of the three-dimensional magneto-resistivity method are improved.
Corresponding to the interpretation method based on the three-dimensional magneto-resistivity method tensor observation mode, the invention also provides a weighted anomaly interpretation method device based on the three-dimensional magneto-resistivity method tensor observation mode. Since the embodiment of the apparatus is similar to the above method embodiment, the description is simple, and please refer to the description of the above method embodiment, and the embodiment of the apparatus for weighted anomaly interpretation method based on three-dimensional magneto-resistivity tensor observation mode described below is only schematic. Fig. 7 is a schematic structural diagram of a weighted anomaly interpretation method and apparatus based on a three-dimensional magneto-resistivity tensor observation mode according to an embodiment of the present invention. The device for the weighted anomaly interpretation method based on the three-dimensional magneto-resistivity tensor observation mode specifically comprises the following parts:
an actually measured magnetic field tensor obtaining unit 701, configured to perform tensor measurement on the target magnetic field under the condition that the two electric dipole sources perpendicular to each other are powered separately, so as to obtain an actually measured magnetic field tensor;
the predicted magnetic field tensor obtaining unit 702 is configured to perform forward calculation by using a three-dimensional magneto-resistivity method, and obtain surface magnetic fields of the two electric dipole sources in the power supply direction based on a finite difference numerical simulation method, so as to obtain a predicted magnetic field tensor;
an abnormal weight factor determining unit 703, configured to obtain a total amplitude of each magnetic field component in the measurement area according to the actually measured magnetic field tensor and the absolute value of the amplitude of each magnetic field component in the predicted magnetic field tensor, and determine an abnormal weight factor of the two electric dipole sources in the corresponding power supply direction based on the total amplitude;
an abnormal normalized magnetic field obtaining unit 704 configured to obtain an abnormal normalized magnetic field according to a ratio of the measured magnetic field tensor to the predicted magnetic field tensor;
an abnormal body position determining unit 705, configured to obtain a comprehensive abnormal normalized magnetic field corresponding to each place on the earth's surface according to the abnormal weight factor and the abnormal normalized magnetic field; and determining the plane position of the abnormal body according to the comprehensive abnormal standardized magnetic field.
By adopting the weighted anomaly interpretation method device based on the three-dimensional magneto-resistivity method tensor observation mode, which is disclosed by the embodiment of the invention, aiming at the three-dimensional complex geological model, the magnetic field data of a plurality of power supply directions are measured by using the three-dimensional magneto-resistivity method tensor observation mode, and the anomaly weight factor is solved through the magnetic field amplitude, so that the comprehensive anomaly standardized magnetic field is obtained, the comprehensive anomaly standardized magnetic field represents the comprehensive strength of the magnetic anomaly of each power supply direction, the anomaly in the complex model is easier to determine, the resolution capability of the three-dimensional complex anomaly is improved, the subsequent solution of the plane position of the three-dimensional complex anomaly is facilitated, and the interpretation efficiency and the accuracy of the three-dimensional magneto-resistivity method are improved.
Corresponding to the interpretation method based on the three-dimensional magneto-resistivity tensor observation mode, the invention also provides electronic equipment. Since the embodiment of the electronic device is similar to the above method embodiment, the description is simple, and please refer to the description of the above method embodiment, and the electronic device described below is only schematic. Fig. 8 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention. The electronic device may include: a processor (processor)801, a memory (memory)802, and a communication bus 803, wherein the processor 801 and the memory 802 communicate with each other via the communication bus 803. The processor 801 may invoke logic instructions in the memory 802 to perform a method of interpretation based on three-dimensional magneto resistivity tensor observed patterns, the method comprising: carrying out tensor measurement on a target magnetic field under the condition that two electric dipole sources which are perpendicular to each other are respectively powered to obtain an actually measured magnetic field tensor; performing forward calculation by using a three-dimensional magneto-resistivity method, and respectively obtaining earth surface magnetic fields of two electric dipole sources in the power supply direction based on a finite difference numerical simulation method to obtain a predicted magnetic field tensor; obtaining the total amplitude of each magnetic field component in a measuring area according to the absolute value of the amplitude of each magnetic field component in the actually measured magnetic field tensor and the predicted magnetic field tensor, and determining the abnormal weight factor of the two electric dipole sources in the corresponding power supply direction based on the total amplitude; obtaining an abnormal standardized magnetic field according to the ratio of the actually measured magnetic field tensor to the predicted magnetic field tensor; obtaining a comprehensive abnormal standardized magnetic field corresponding to each place of the earth surface according to the abnormal weight factor and the abnormal standardized magnetic field; and determining the plane position of the abnormal body according to the comprehensive abnormal standardized magnetic field.
Furthermore, the logic instructions in the memory 802 may be implemented in software functional units and stored in a computer readable storage medium when sold or used as a stand-alone product. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, an embodiment of the present invention further provides a computer program product, where the computer program product includes a computer program stored on a non-transitory computer-readable storage medium, where the computer program includes program instructions, and when the program instructions are executed by a computer, the computer is capable of executing the three-dimensional magneto-resistivity tensor observation mode-based interpretation method provided by the above-mentioned method embodiments, where the method includes: carrying out tensor measurement on a target magnetic field under the condition that two electric dipole sources which are perpendicular to each other are respectively powered to obtain an actually measured magnetic field tensor; performing forward calculation by using a three-dimensional magneto-resistivity method, and respectively obtaining earth surface magnetic fields of two electric dipole sources in the power supply direction based on a finite difference numerical simulation method to obtain a predicted magnetic field tensor; acquiring the total amplitude of each magnetic field component in a measurement area according to the absolute value of the amplitude of each magnetic field component in the actually measured magnetic field tensor and the predicted magnetic field tensor, and determining the abnormal weight factor of the two electric dipole sources in the corresponding power supply direction based on the total amplitude; obtaining an abnormal standardized magnetic field according to the ratio of the actually measured magnetic field tensor to the predicted magnetic field tensor; obtaining a comprehensive abnormal standardized magnetic field corresponding to each place of the earth surface according to the abnormal weight factor and the abnormal standardized magnetic field; and determining the plane position of the abnormal body according to the comprehensive abnormal standardized magnetic field.
In another aspect, an embodiment of the present invention further provides a non-transitory computer-readable storage medium, on which a computer program is stored, where the computer program is implemented by a processor to perform the three-dimensional magneto-resistivity tensor observation mode-based interpretation method provided by the foregoing embodiments, the method including: carrying out tensor measurement on a target magnetic field under the condition that two electric dipole sources which are perpendicular to each other are respectively powered to obtain an actually measured magnetic field tensor; performing forward calculation by using a three-dimensional magneto-resistivity method, and respectively obtaining earth surface magnetic fields of two electric dipole sources in the power supply direction based on a finite difference numerical simulation method to obtain a predicted magnetic field tensor; acquiring the total amplitude of each magnetic field component in a measurement area according to the absolute value of the amplitude of each magnetic field component in the actually measured magnetic field tensor and the predicted magnetic field tensor, and determining the abnormal weight factor of the two electric dipole sources in the corresponding power supply direction based on the total amplitude; obtaining an abnormal standardized magnetic field according to the ratio of the actually measured magnetic field tensor to the predicted magnetic field tensor; obtaining a comprehensive abnormal standardized magnetic field corresponding to each place of the earth surface according to the abnormal weight factor and the abnormal standardized magnetic field; and determining the plane position of the abnormal body according to the comprehensive abnormal standardized magnetic field.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one position, or may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. An interpretation method based on a three-dimensional magneto-resistivity tensor observation mode is characterized by comprising the following steps:
tensor measurement is carried out on a target magnetic field under the condition that two electric dipole sources which are perpendicular to each other are respectively powered, and an actually measured magnetic field tensor is obtained;
performing forward calculation by using a three-dimensional magneto-resistivity method, and respectively obtaining earth surface magnetic fields of two electric dipole sources in the power supply direction based on a finite difference numerical simulation method to obtain a predicted magnetic field tensor;
acquiring the total amplitude of each magnetic field component in a measurement area according to the absolute value of the amplitude of each magnetic field component in the actually measured magnetic field tensor and the predicted magnetic field tensor, and determining the abnormal weight factor of the two electric dipole sources in the corresponding power supply direction based on the total amplitude;
obtaining an abnormal standardized magnetic field according to the ratio of the actually measured magnetic field tensor to the predicted magnetic field tensor;
obtaining a comprehensive abnormal standardized magnetic field corresponding to each place of the earth surface according to the abnormal weight factor and the abnormal standardized magnetic field; and determining the plane position of the abnormal body according to the comprehensive abnormal standardized magnetic field.
2. The interpretation method based on the three-dimensional magneto-resistivity tensor observation mode as recited in claim 1, wherein the measured magnetic field tensor corresponds to an expression:
wherein Hxx、HxyA magnetic field component corresponding to the X, Y power supply direction actually measured on the X-axis power supply ground surface is provided for the electric dipole source central line; hyx、HyyAnd a magnetic field component corresponding to the X, Y power supply direction measured on the power supply ground surface of the Y-axis for the center line of the electric dipole source.
3. The interpretation method based on the three-dimensional magneto-resistivity tensor observation mode as recited in claim 1, wherein the predicted magnetic field tensor corresponds to an expression:
wherein Hxx′、Hxy' a magnetic field component corresponding to X, Y power supply direction predicted by the power supply surface of the X axis of the center line of the electric dipole source; hyx′、Hyy' the predicted X, Y power supply direction corresponding magnetic field component of the Y-axis power supply surface for the center line of the electric dipole source.
4. The interpretation method based on the three-dimensional magneto-resistivity tensor observation mode as recited in claim 1, wherein the abnormal weight factors include a first abnormal weight factor ω X corresponding to the power supply in the X-axis direction and a second abnormal weight factor ω Y corresponding to the power supply in the Y-axis direction;
the calculation formula corresponding to the first abnormality weight ω x is:
therein, sigmaΩ|Hxx|、∑Ω|HxyI is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the actually measured magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction; sigmaΩ|Hxx′|、∑Ω|Hxy' I is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the predicted magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction;
the calculation formula corresponding to the second abnormality weight ω y is as follows:
therein, sigmaΩ|Hyy|、∑Ω|HyxI is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the actually measured magnetic field tensor in the measurement region omega when power is supplied in the X-axis direction; sigmaΩ|Hyy′|、∑Ω|Hyx' l is the sum of the absolute values of the amplitudes of the X magnetic field component and the Y magnetic field component corresponding to the predicted magnetic field tensor in the measurement region Ω when power is supplied in the X-axis direction, respectively.
5. The interpretation method based on the three-dimensional magneto-resistivity tensor observation mode as recited in claim 1, wherein the abnormal normalized magnetic field corresponds to a calculation formula:
wherein HxsFor abnormal standardized magnetic field, H, of electric dipole source central line when power is supplied on X axisysAn abnormal standardized magnetic field is generated when the electric dipole source central line supplies power on the Y axis; the abnormal standardized magnetic field is used for determining the relative strength of the magnetic field or the relative strength of current distribution in the measurement region under the homogeneous condition;
abnormal normalized magnetic field H at X-axis power supplyxsCorresponding calculation formula and abnormal standardized magnetic field H in power supply of Y axisysThe corresponding calculation formulas are respectively as follows:
wherein HxyAnd HxyRespectively measuring and predicting the magnetic field component H corresponding to the Y power supply direction of the electric dipole source central line on the X-axis power supply ground surfaceyxAnd HyxThe' are respectively the magnetic field components corresponding to the X power supply direction actually measured and predicted by the electric dipole source central line on the Y-axis power supply ground surface.
6. The interpretation method based on the three-dimensional magneto-resistivity tensor observation mode as claimed in claim 1, wherein the comprehensive abnormal normalized magnetic field is a ratio of a sum of products of abnormal weight factors and abnormal normalized magnetic fields in each power supply direction to a sum of abnormal weight factors in each power supply direction, and the specific calculation formula is as follows:
wherein, ω X is a first abnormal weight factor corresponding to the power supply in the X-axis direction, and ω Y is a second abnormal weight factor corresponding to the power supply in the Y-axis direction; hxsFor abnormal standardized magnetic field, H, of electric dipole source central line when power is supplied on X axisysAnd (4) an abnormal standardized magnetic field is generated when the central line of the electric dipole source supplies power on the Y axis.
7. The interpretation method based on the three-dimensional magneto-resistivity tensor observation mode according to claim 4, wherein a larger first abnormal weight factor ω x indicates a larger current of the measurement region parallel to the power supply direction, and a smaller first abnormal weight factor ω x indicates a smaller current of the measurement region parallel to the power supply direction; the larger the second abnormality weight ω y, the larger the current of the measurement region parallel to the power supply direction, and the smaller the second abnormality weight ω y, the smaller the current of the measurement region parallel to the power supply direction.
8. An interpretation apparatus based on three-dimensional magneto-resistivity tensor observation mode, comprising:
the actual measurement magnetic field tensor obtaining unit is used for carrying out tensor measurement on the target magnetic field under the condition that two electric dipole sources which are perpendicular to each other are respectively powered to obtain an actual measurement magnetic field tensor;
the predicted magnetic field tensor obtaining unit is used for performing forward calculation by using a three-dimensional magneto-resistivity method, and obtaining the earth surface magnetic fields of the two electric dipole sources in the power supply direction respectively based on a finite difference numerical simulation method to obtain a predicted magnetic field tensor;
the abnormal weight factor determining unit is used for obtaining the total amplitude of each magnetic field component in the measurement area according to the absolute value of the amplitude of each magnetic field component in the actually measured magnetic field tensor and the predicted magnetic field tensor, and determining the abnormal weight factors of the two electric dipole sources in the corresponding power supply directions based on the total amplitude;
an abnormal normalized magnetic field obtaining unit configured to obtain an abnormal normalized magnetic field according to a ratio of the measured magnetic field tensor to the predicted magnetic field tensor;
the abnormal body position determining unit is used for obtaining a comprehensive abnormal standardized magnetic field corresponding to each position of the earth surface according to the abnormal weight factor and the abnormal standardized magnetic field; and determining the plane position of the abnormal body according to the comprehensive abnormal standardized magnetic field.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program implements the steps of the method for interpreting based on three-dimensional magneto-resistivity tensor observation patterns according to any one of claims 1-7.
10. A non-transitory computer readable storage medium, having stored thereon a computer program, which, when being executed by a processor, carries out the steps of the method for interpreting based on three-dimensional magneto-resistivity tensor observation patterns according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110199433.4A CN113109876B (en) | 2021-02-22 | 2021-02-22 | Interpretation method and device based on three-dimensional magneto-resistivity method tensor observation mode |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110199433.4A CN113109876B (en) | 2021-02-22 | 2021-02-22 | Interpretation method and device based on three-dimensional magneto-resistivity method tensor observation mode |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113109876A CN113109876A (en) | 2021-07-13 |
| CN113109876B true CN113109876B (en) | 2022-06-10 |
Family
ID=76710092
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110199433.4A Active CN113109876B (en) | 2021-02-22 | 2021-02-22 | Interpretation method and device based on three-dimensional magneto-resistivity method tensor observation mode |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113109876B (en) |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4792761A (en) * | 1987-04-06 | 1988-12-20 | King Thomas C | Geophysical prospecting with collimated magnetotelluric fields |
| US7383132B2 (en) * | 2005-03-07 | 2008-06-03 | Exxonmobil Upstream Research Co. | Method for identifying resistivity anomalies in electromagnetic survey data |
| GB2452871B (en) * | 2006-04-06 | 2010-11-03 | Exxonmobil Upstream Res Co | Method for obtaining resistivity from controlled source electromagnetic data |
| US10955579B2 (en) * | 2015-03-20 | 2021-03-23 | Halliburton Energy Services, Inc. | Measurement and control apparatus, systems, and methods |
| CN105158808B (en) * | 2015-07-01 | 2017-10-31 | 长安大学 | A kind of shallow sea transient electromagnetic air-sea detection and its means of interpretation |
| CN105467461B (en) * | 2015-12-01 | 2017-11-17 | 长江地球物理探测（武汉）有限公司 | A kind of method that geological anomalous body is identified using two-dimentional apparent resistivity data |
| CN110346834B (en) * | 2019-07-22 | 2020-11-17 | 中国科学院地球化学研究所 | Forward modeling method and system for three-dimensional frequency domain controllable source electromagnetism |
-
2021
- 2021-02-22 CN CN202110199433.4A patent/CN113109876B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN113109876A (en) | 2021-07-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102798898B (en) | Three-dimensional inversion method for nonlinear conjugate gradient of magnetotelluric field | |
| Persova et al. | Finite element solution to 3-D airborne time-domain electromagnetic problems in complex geological media using non-conforming hexahedral meshes | |
| CN109209338B (en) | Electrical observation system and detection method for detecting abnormities beside well | |
| Farquharson et al. | An algorithm for the three-dimensional inversion of magnetotelluric data | |
| CN109254327B (en) | Exploration method and exploration system of three-dimensional ferromagnetic body | |
| CN111722292B (en) | Three-dimensional resistivity tomography method and system | |
| Lu et al. | 3D finite-volume time-domain modeling of geophysical electromagnetic data on unstructured grids using potentials | |
| CN104656156A (en) | Magnetic reference processing method for acquiring data in three dimensions audio-frequency magnetotelluric sounding | |
| CN110488365B (en) | Multi-polarization magnetotelluric depth measurement method | |
| CN113109876B (en) | Interpretation method and device based on three-dimensional magneto-resistivity method tensor observation mode | |
| CN112503399B (en) | Method and device for determining underground pipeline position and storage medium thereof | |
| CN110119586A (en) | Axial conductivity anisotropy transient electromagnetic three-component three-dimensional FDTD forward modeling method | |
| Slama et al. | Coupling the electromagnetic inverse problem based on genetic algorithms with Moment's method for EMC of circuits | |
| Oh et al. | Interpretation of controlled-source electromagnetic data from iron ores under rough topography | |
| Feng et al. | Contrast between 2D inversion and 3D inversion based on 2D high-density resistivity data | |
| CN115657142A (en) | Method, device, equipment and storage medium for detecting and monitoring water damage of mine stratum | |
| Endo et al. | A multigrid integral equation method for large-scale models with inhomogeneous backgrounds | |
| Lu et al. | 3D large-scale transient electromagnetic modeling using a Shift-and-Invert Krylov subspace method | |
| Ruo et al. | The analysis of CSAMT responses of dyke embedded below conductive overburden | |
| Muttaqien et al. | Two-dimensional inversion modeling of magnetotelluric (MT) synthetic data of a graben structure using SimPEG | |
| CN111597752A (en) | Cross-hole resistivity CT deep learning inversion method and system for balancing inter-hole sensitivity | |
| Turarova et al. | Evaluation of the 3D Topographic Effect of Homogeneous and Inhomogeneous Media on the Results of 2D Inversion of Electrical Resistivity Tomography Data | |
| Chen et al. | Fast and High Precision Numerical Simulation of DC 2.5 D under Multi-source Condition | |
| CN113221411B (en) | Charging potential numerical simulation method, system and terminal for lossy medium with any shape | |
| Imamura et al. | Application of marine controlled-source electromagnetic sounding to submarine massive sulphides explorations |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |