CN113917541A - Method and device for acquiring electrical anisotropy of underground medium based on direct current electrical prospecting - Google Patents
Method and device for acquiring electrical anisotropy of underground medium based on direct current electrical prospecting Download PDFInfo
- Publication number
- CN113917541A CN113917541A CN202111157020.6A CN202111157020A CN113917541A CN 113917541 A CN113917541 A CN 113917541A CN 202111157020 A CN202111157020 A CN 202111157020A CN 113917541 A CN113917541 A CN 113917541A
- Authority
- CN
- China
- Prior art keywords
- background
- apparent resistivity
- resistivity
- potential difference
- normalized
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000004364 calculation method Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims 1
- 208000001491 myopia Diseases 0.000 abstract description 4
- 238000005259 measurement Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 230000005611 electricity Effects 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000010291 electrical method Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
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/02—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current
- G01V3/04—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current using DC
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
- G01N27/041—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a solid body
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
The application relates to a method and a device for acquiring electrical anisotropy of an underground medium based on direct current electrical prospecting. The method comprises the following steps: the method comprises the steps of arranging a transmitting field source at any position of a measuring area, carrying out whole-area observation in the measuring area by adopting mutually orthogonal electrodes to obtain mutually orthogonal potential differences in two directions between different measuring points, calculating apparent resistivity according to the potential differences, calculating normalized apparent resistivity, calculating background apparent resistivity according to a selected background resistivity model, calculating normalized background apparent resistivity, and calculating a resistivity anisotropy coefficient according to actual normalized resistivity and normalized background apparent resistivity. The field source layout is flexible and simple. The influence of the apparent resistivity zero line can be effectively suppressed by adopting the resistivity anisotropy coefficient, the projection and the shape of the underground anisotropic structure body on the ground can be well reflected, tensor resistivity can be avoided to be solved, and the efficiency is high. The error of the near-sighted resistivity attached to the field source can be effectively offset by adopting the normalized apparent resistivity, and the precision is improved.
Description
Technical Field
The application relates to the field of geophysical exploration, in particular to a method, a device, computer equipment and a storage medium for acquiring electrical anisotropy of an underground medium based on direct current electrical prospecting.
Background
The direct current method is a relatively mature geophysical method and is widely applied to the fields of mineral exploration, hydrological environment detection, engineering exploration and the like. The rock on the earth surface layer forms a complex structure and fault due to uneven composition, geological action and the like, has various anisotropic characteristics, and can acquire the electrical anisotropy of the underground medium by observing in the earth surface or in a well by using a direct current electrical methodAnd (5) characterizing. The current major exploration method for dc electrical methods is the P2 invariant apparent resistivity tensor measurement proposed by Bibby (1986). The measuring method comprises emitting from outside of the measuring region by using orthogonal dipole sources AB and CB, and measuring electric field in the measuring region by using orthogonal electrodes
In addition, the current density of the corresponding isotropic background model is calculated
Finally using the expression
Calculating tensor resistivity
In combination with
To characterize the apparent resistivity of the subsurface anisotropic media. Obviously, the method needs to arrange two dipole field sources outside the measuring area, and in addition, tensor resistivity must be solved, and finally, P is obtained through calculation2The method represents the resistivity anisotropy characteristic of the underground medium, so that the field source is complicated to arrange, the processing process is complicated, the working efficiency is low, and particularly when the data volume is large, the processing is quite troublesome.
Disclosure of Invention
Therefore, it is necessary to provide a method, an apparatus, a computer device and a storage medium for acquiring electrical anisotropy of an underground medium based on dc electrical prospecting, which can realize full-area observation, simple field source layout, and simple and reliable data processing.
A method for obtaining electrical anisotropy of a subsurface medium based on direct current electrical prospecting, said method comprising:
after an emission field source is arranged at any position of a measuring area, performing whole-area observation in the measuring area by adopting electrodes which are orthogonal to each other to obtain a first-direction potential difference and a second-direction potential difference between different measuring points; the emission field source is a point source or a dipole source; the first direction and the second direction are orthogonal to each other;
obtaining a first direction apparent resistivity and a second direction apparent resistivity according to the first direction potential difference and the second direction potential difference, and obtaining a normalized apparent resistivity according to the first direction apparent resistivity and the second direction apparent resistivity;
calculating the potential difference of the background between different measuring points in the first direction and the potential difference of the background in the second direction according to a preset resistivity model of the underground medium background;
obtaining a background first direction apparent resistivity and a background second direction apparent resistivity according to the background first direction potential difference and the background second direction potential difference, and obtaining a normalized background apparent resistivity according to the background first direction apparent resistivity and the background second direction apparent resistivity;
and obtaining the resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity.
In one embodiment, the method further comprises the following steps: the emission field source is a point source or a dipole source; when the emission field source is a dipole source, the emission field source can be arranged along any direction.
In one embodiment, the method further comprises the following steps: the mutually orthogonal electrodes may be laid in any direction.
In one embodiment, the method further comprises the following steps: obtaining a first direction apparent resistivity and a second direction apparent resistivity according to the first direction potential difference and the second direction potential difference, and obtaining a background first direction apparent resistivity and a background second direction apparent resistivity according to the background first direction potential difference and the background second direction potential difference, wherein the calculation formula is as follows:
where ρ is apparent resistivity, and may be apparent resistivity ρ in the first direction, respectivelyxOf 1 atApparent resistivity rho in two directionsyBackground first direction apparent resistivity
Background second direction apparent resistivity
Δ U is potential difference corresponding to background apparent resistivity Δ UxPotential difference Δ U in the second directionyBackground first direction potential difference
Background second direction potential difference
K is the device coefficient and I is the current intensity.
In one embodiment, the method further comprises the following steps: the formula for obtaining the normalized apparent resistivity and the normalized background apparent resistivity through calculation is as follows:
where ρ isxyIn order to normalize the apparent resistivity,
normalized to background apparent resistivity.
In one embodiment, the method further comprises the following steps: obtaining a resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity, wherein the resistivity anisotropy coefficient comprises the following steps:
and obtaining the resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity, wherein the resistivity anisotropy coefficient is as follows:
wherein, gxyRepresenting the resistivity anisotropy coefficient of the subsurface medium.
In one embodiment, the method further comprises the following steps: the emission field source can be deployed at the surface or in a well.
An apparatus for obtaining electrical anisotropy of a subsurface medium based on direct current electrical prospecting, said apparatus comprising:
the potential difference observation module is used for carrying out whole-area observation in the measuring area by adopting electrodes which are orthogonal to each other after an emission field source is arranged at any position of the measuring area to obtain potential differences in a first direction and potential differences in a second direction between different measuring points; the emission field source is a point source or a dipole source; the first direction and the second direction are orthogonal to each other;
the normalized apparent resistivity determining module is used for obtaining a first direction apparent resistivity and a second direction apparent resistivity according to the first direction potential difference and the second direction potential difference, and obtaining a normalized apparent resistivity according to the first direction apparent resistivity and the second direction apparent resistivity;
the background potential difference calculation module is used for calculating a background first-direction potential difference and a background second-direction potential difference between different measuring points according to a preset underground medium background resistivity model;
the normalized background apparent resistivity determining module is used for obtaining background first direction apparent resistivity and background second direction apparent resistivity according to the background first direction potential difference and the background second direction potential difference, and obtaining normalized background apparent resistivity according to the background first direction apparent resistivity and the background second direction apparent resistivity;
and the resistivity anisotropy coefficient determination module is used for obtaining the resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity.
A computer device comprising a memory and a processor, the memory storing a computer program, the processor implementing the following steps when executing the computer program:
after an emission field source is arranged at any position of a measuring area, performing whole-area observation in the measuring area by adopting electrodes which are orthogonal to each other to obtain a first-direction potential difference and a second-direction potential difference between different measuring points; the emission field source is a point source or a dipole source; the first direction and the second direction are orthogonal to each other;
obtaining a first direction apparent resistivity and a second direction apparent resistivity according to the first direction potential difference and the second direction potential difference, and obtaining a normalized apparent resistivity according to the first direction apparent resistivity and the second direction apparent resistivity;
calculating the potential difference of the background between different measuring points in the first direction and the potential difference of the background in the second direction according to a preset resistivity model of the underground medium background;
obtaining a background first direction apparent resistivity and a background second direction apparent resistivity according to the background first direction potential difference and the background second direction potential difference, and obtaining a normalized background apparent resistivity according to the background first direction apparent resistivity and the background second direction apparent resistivity;
and obtaining the resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity.
A computer-readable storage medium, on which a computer program is stored which, when executed by a processor, carries out the steps of:
after an emission field source is arranged at any position of a measuring area, performing whole-area observation in the measuring area by adopting electrodes which are orthogonal to each other to obtain a first-direction potential difference and a second-direction potential difference between different measuring points; the emission field source is a point source or a dipole source; the first direction and the second direction are orthogonal to each other;
obtaining a first direction apparent resistivity and a second direction apparent resistivity according to the first direction potential difference and the second direction potential difference, and obtaining a normalized apparent resistivity according to the first direction apparent resistivity and the second direction apparent resistivity;
calculating the potential difference of the background between different measuring points in the first direction and the potential difference of the background in the second direction according to a preset resistivity model of the underground medium background;
obtaining a background first direction apparent resistivity and a background second direction apparent resistivity according to the background first direction potential difference and the background second direction potential difference, and obtaining a normalized background apparent resistivity according to the background first direction apparent resistivity and the background second direction apparent resistivity;
and obtaining the resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity.
According to the method, the device, the computer equipment and the storage medium for acquiring the electrical anisotropy of the underground medium based on the direct current electrical prospecting, after the emission field source is arranged at any position of the measurement area, the electrodes which are orthogonal to each other are adopted to carry out whole-area observation in the measurement area, so that the potential differences in two directions which are orthogonal to each other between different measurement points are obtained, the apparent resistivity is calculated according to the potential differences, then the normalized apparent resistivity is calculated, the background apparent resistivity is calculated according to the selected background resistivity model, then the normalized background apparent resistivity is calculated, and the resistivity anisotropy coefficient is calculated according to the actual normalized resistivity and the normalized background apparent resistivity. The field source arranged by the invention only needs a point source or a dipole source, and the arrangement is flexible and simple. The characteristic of the electrical resistivity anisotropy coefficient representation underground electricity anisotropy is adopted, the influence of an apparent resistivity zero line can be effectively suppressed, and the projection and the shape of the underground anisotropic structure body on the ground can be well reflected. In addition, the complex calculation of solving tensor resistivity can be avoided, and the technical efficiency is improved. By adopting the normalized apparent resistivity, the error of the near-sighted resistivity of the field source can be effectively offset, and the calculation precision is improved.
Drawings
FIG. 1 is a schematic flow chart illustrating a method for obtaining electrical anisotropy of a subsurface medium based on DC electrical prospecting in one embodiment;
FIG. 2 is a schematic diagram of an embodiment of the apparatus;
FIG. 3 is a schematic representation of a model in one embodiment;
FIG. 4 is an apparent resistivity p measured along the x-direction in one embodimentxA plan view;
FIG. 5 is an apparent resistivity ρ measured in the y-direction in one embodimentyA plan view;
FIG. 6 is an example normalized resistivity anisotropy coefficient g for a specific embodimentxyA plan view;
FIG. 7 is a block diagram of an apparatus for obtaining electrical anisotropy of a subsurface medium based on DC electrical prospecting in one embodiment;
FIG. 8 is a diagram illustrating an internal structure of a computer device according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
The method for acquiring the electrical anisotropy of the underground medium based on direct current electrical prospecting can be applied to the following application environments. The terminal executes a method for acquiring the electrical anisotropy of the underground medium based on direct current electrical prospecting, wherein after a transmitting field source is arranged at any position of a measuring area, the whole area observation is carried out in the measuring area by adopting mutually orthogonal electrodes to obtain the potential differences in two mutually orthogonal directions between different measuring points, the apparent resistivity is calculated according to the potential differences, then the normalized apparent resistivity is calculated, the background apparent resistivity is calculated according to a selected background resistivity model, then the normalized background apparent resistivity is calculated, and the resistivity anisotropy coefficient is calculated according to the actual normalized resistivity and the normalized background apparent resistivity. The terminal may be, but is not limited to, various personal computers, notebook computers, tablet computers, and portable devices.
In one embodiment, as shown in FIG. 1, there is provided a method for obtaining electrical anisotropy of a subsurface medium based on direct current electrical prospecting, comprising the steps of:
102, arranging an emission field source at any position of a measuring area, and carrying out whole-area observation in the measuring area by adopting electrodes which are orthogonal to each other to obtain a first-direction potential difference and a second-direction potential difference between different measuring points.
The first direction and the second direction are orthogonal to each other. The emission field source is a point source or a dipole source, and the arrangement is flexible and simple.
And 104, obtaining the apparent resistivity of the first direction and the apparent resistivity of the second direction according to the potential difference of the first direction and the potential difference of the second direction, and obtaining the normalized apparent resistivity according to the apparent resistivity of the first direction and the apparent resistivity of the second direction.
By adopting the normalized apparent resistivity, the error of the near-sighted resistivity of the field source can be effectively offset, and the calculation precision is improved.
And 106, calculating the potential difference of the background in the first direction and the potential difference of the background in the second direction between different measuring points according to a preset resistivity model of the underground medium background.
And 108, obtaining the background apparent resistivity and the background apparent resistivity in the first direction according to the background potential difference in the first direction and the background potential difference in the second direction, and obtaining the normalized background apparent resistivity according to the background apparent resistivity in the first direction and the background apparent resistivity in the second direction.
And step 110, obtaining the resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity.
The characteristic of the electrical resistivity anisotropy coefficient representation underground electricity anisotropy is adopted, the influence of an apparent resistivity zero line can be effectively suppressed, and the projection and the shape of the underground anisotropic structure body on the ground can be well reflected. In addition, the complex calculation of solving tensor resistivity can be avoided, and the technical efficiency is improved.
In the method for acquiring the electrical anisotropy of the underground medium based on the direct current electrical prospecting, after a transmitting field source is arranged at any position of a measuring area, the whole area observation is carried out in the measuring area by adopting electrodes which are orthogonal to each other, so that potential differences in two directions which are orthogonal to each other between different measuring points are obtained, the apparent resistivity is calculated according to the potential differences, then the normalized apparent resistivity is calculated, the background apparent resistivity is calculated according to a selected background resistivity model, then the normalized background apparent resistivity is calculated, and the resistivity anisotropy coefficient is calculated according to the actual normalized resistivity and the normalized background apparent resistivity. The field source arranged by the invention only needs a point source or a dipole source, and the arrangement is flexible and simple. The characteristic of the electrical resistivity anisotropy coefficient representation underground electricity anisotropy is adopted, the influence of an apparent resistivity zero line can be effectively suppressed, and the projection and the shape of the underground anisotropic structure body on the ground can be well reflected. In addition, the complex calculation of solving tensor resistivity can be avoided, and the technical efficiency is improved. By adopting the normalized apparent resistivity, the error of the near-sighted resistivity of the field source can be effectively offset, and the calculation precision is improved.
In one embodiment, the method further comprises the following steps: the emission field source is a point source or a dipole source; when the emission field source is a dipole source, the emission field source can be arranged along any direction. The layout is flexible and simple.
In one embodiment, the method further comprises the following steps: the mutually orthogonal electrodes may be laid in any direction. The layout is flexible and simple.
In one embodiment, the method further comprises the following steps: obtaining the apparent resistivity of the first direction and the apparent resistivity of the second direction according to the potential difference of the first direction and the potential difference of the second direction, and obtaining the apparent resistivity of the first direction and the apparent resistivity of the second direction according to the potential difference of the first direction and the potential difference of the second direction of the background, wherein the calculation formula is as follows:
where ρ is apparent resistivity, and may be apparent resistivity ρ in the first direction, respectivelyxApparent resistivity in the second direction ρyBackground first direction apparent resistivity
Background second direction apparent resistivity
Δ U is potential difference corresponding to background apparent resistivity Δ UxPotential difference Δ U in the second directionyBackground first direction potential difference
Background second direction potential difference
K is the device coefficient and I is the current intensity.
In one embodiment, the method further comprises the following steps: the formula for obtaining the normalized apparent resistivity and the normalized background apparent resistivity through calculation is as follows:
where ρ isxyIn order to normalize the apparent resistivity,
normalized to background apparent resistivity.
In one embodiment, the method further comprises the following steps: obtaining the resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity, wherein the resistivity anisotropy coefficient comprises the following steps:
and obtaining the resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity, wherein the resistivity anisotropy coefficient is as follows:
wherein, gxyRepresenting the resistivity anisotropy coefficient of the subsurface medium.
In one embodiment, the method further comprises the following steps: the emission field source can be deployed at the surface or in the well.
It should be understood that, although the steps in the flowchart of fig. 1 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least a portion of the steps in fig. 1 may include multiple sub-steps or multiple stages that are not necessarily performed at the same time, but may be performed at different times, and the order of performance of the sub-steps or stages is not necessarily sequential, but may be performed in turn or alternately with other steps or at least a portion of the sub-steps or stages of other steps.
Fig. 2 is a schematic diagram of the apparatus of the present invention in an embodiment, wherein M1 and N1 are two electrodes for measuring a potential difference in a first direction, respectively, and M2 and N2 are two electrodes for measuring a potential difference in a second direction, respectively. The distance of M1N1 is equal to the distance of M2N 2.
In one embodiment, a model of an anisotropic anomaly in an isotropic half-space is shown in FIG. 3, wherein the anomaly is a prism whose center point is projected on the ground as the origin of a coordinate system, the prism has a size of 20m × 20m × 10m, a ground burial depth of 10m, and three principal axis resistivities ρ1/ρ2/ρ310/5/10 Ω · m, the euler angle α/β/χ is 30 °/45 °/60 °, and the background resistivity is ρ0100 Ω · m. A point source was used as a field source, located at a position with coordinates (200m, 200m, 0). The receiving electrodes are orthogonal and measure the potentials (differences) in the x-direction and y-direction, respectively. FIG. 4 is apparent resistivity ρ measured along the x-directionxA plan view. FIG. 5 is an apparent resistivity ρ measured in the y-direction for an example of the inventionyA plan view. FIG. 6 is a normalized resistivity anisotropy coefficient gxyA plan view. The black boxes in fig. 4, 5, and 6 indicate the projection position and size of the abnormal body on the ground. As can be seen from FIGS. 4, 5, and 6, the use of the normalized resistivity anisotropy coefficient better reflects the in-situ anisotropyThe position and shape of the table projection.
In one embodiment, as shown in fig. 7, there is provided an apparatus for obtaining electrical anisotropy of a subsurface medium based on direct current electrical prospecting, comprising: potential difference observation module 702, normalized apparent resistivity determination module 704, background potential difference calculation module 706, normalized background apparent resistivity determination module 708, and resistivity anisotropy coefficient determination module 710, wherein:
the potential difference observation module 702 is used for performing whole-area observation in the measurement area by adopting mutually orthogonal electrodes after an emission field source is arranged at any position of the measurement area to obtain a first-direction potential difference and a second-direction potential difference between different measurement points; the first direction and the second direction are mutually orthogonal;
the normalized apparent resistivity determining module 704 is configured to obtain a first direction apparent resistivity and a second direction apparent resistivity according to the first direction potential difference and the second direction potential difference, and obtain a normalized apparent resistivity according to the first direction apparent resistivity and the second direction apparent resistivity;
the background potential difference calculation module 706 is configured to calculate a background first-direction potential difference and a background second-direction potential difference between different measurement points according to a preset underground medium background resistivity model;
a normalized background apparent resistivity determining module 708, configured to obtain a background first-direction apparent resistivity and a background second-direction apparent resistivity according to the background first-direction potential difference and the background second-direction potential difference, and obtain a normalized background apparent resistivity according to the background first-direction apparent resistivity and the background second-direction apparent resistivity;
and the resistivity anisotropy coefficient determination module 710 is used for obtaining the resistivity anisotropy coefficient of the underground medium of the measured area according to the normalized apparent resistivity and the normalized background apparent resistivity.
The resistivity anisotropy coefficient determination module 710 is further configured to obtain a resistivity anisotropy coefficient of the subsurface medium of the measurement area according to the normalized apparent resistivity and the normalized background apparent resistivity, where the resistivity anisotropy coefficient is:
wherein, gxyRepresenting the resistivity anisotropy coefficient of the subsurface medium.
For specific limitations of the apparatus for acquiring electrical anisotropy of subsurface medium based on dc electrical prospecting, reference may be made to the above limitations of the method for acquiring electrical anisotropy of subsurface medium based on dc electrical prospecting, which are not repeated herein. The modules in the device for acquiring the electrical anisotropy of the underground medium based on the direct current electrical prospecting can be wholly or partially realized by software, hardware and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as shown in fig. 8. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to implement a method for obtaining electrical anisotropy of a subsurface medium based on direct current electrical prospecting. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on the shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.
Those skilled in the art will appreciate that the architecture shown in fig. 8 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In an embodiment, a computer device is provided, comprising a memory storing a computer program and a processor implementing the steps of the above method embodiments when executing the computer program.
In an embodiment, a computer-readable storage medium is provided, on which a computer program is stored, which computer program, when being executed by a processor, carries out the steps of the above-mentioned method embodiments.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory, among others. Non-volatile memory can include read-only memory (ROM), Programmable ROM (PROM), Electrically Programmable ROM (EPROM), Electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Link DRAM (SLDRAM), Rambus Direct RAM (RDRAM), direct bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (10)
1. A method for acquiring electrical anisotropy of a subsurface medium based on direct current electrical prospecting, the method comprising:
after an emission field source is arranged at any position of a measuring area, performing whole-area observation in the measuring area by adopting electrodes which are orthogonal to each other to obtain a first-direction potential difference and a second-direction potential difference between different measuring points; the emission field source is a point source or a dipole source; the first direction and the second direction are orthogonal to each other;
obtaining a first direction apparent resistivity and a second direction apparent resistivity according to the first direction potential difference and the second direction potential difference, and obtaining a normalized apparent resistivity according to the first direction apparent resistivity and the second direction apparent resistivity;
calculating the potential difference of the background between different measuring points in the first direction and the potential difference of the background in the second direction according to a preset resistivity model of the underground medium background;
obtaining a background first direction apparent resistivity and a background second direction apparent resistivity according to the background first direction potential difference and the background second direction potential difference, and obtaining a normalized background apparent resistivity according to the background first direction apparent resistivity and the background second direction apparent resistivity;
and obtaining the resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity.
2. The method of claim 1, wherein the transmission field source can be deployed in any direction when it is a dipole source.
3. The method of claim 1, wherein the mutually orthogonal electrodes can be laid in any direction.
4. The method of claim 3, wherein the calculation formula for obtaining the first direction apparent resistivity and the second direction apparent resistivity according to the first direction potential difference and the second direction potential difference, and the calculation formula for obtaining the background first direction apparent resistivity and the background second direction apparent resistivity according to the background first direction potential difference and the background second direction potential difference is as follows:
where ρ is apparent resistivity, and may be apparent resistivity ρ in the first direction, respectivelyxApparent resistivity in the second direction ρyBackground first direction apparent resistivity
Background second direction apparent resistivity
Δ U is potential difference corresponding to background apparent resistivity Δ UxPotential difference Δ U in the second directionyBackground first direction potential difference
Background second direction potential difference
K is the device coefficient and I is the current intensity.
5. The method of claim 4, wherein the formula for calculating the normalized apparent resistivity and the normalized background apparent resistivity is:
where ρ isxyIn order to normalize the apparent resistivity,
normalized to background apparent resistivity.
6. The method of claim 5, wherein obtaining a resistivity anisotropy coefficient for a subsurface medium of the survey area based on the normalized apparent resistivity and the normalized background apparent resistivity comprises:
and obtaining the resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity, wherein the resistivity anisotropy coefficient is as follows:
wherein, gxyRepresenting the resistivity anisotropy coefficient of the subsurface medium.
7. The method of any one of claims 1 to 6, wherein the launch field source is deployable in the earth's surface or in a well.
8. An apparatus for obtaining electrical anisotropy of a subsurface medium based on direct current electrical prospecting, the apparatus comprising:
the potential difference observation module is used for carrying out whole-area observation in the measuring area by adopting electrodes which are orthogonal to each other after an emission field source is arranged at any position of the measuring area to obtain potential differences in a first direction and potential differences in a second direction between different measuring points; the emission field source is a point source or a dipole source; the first direction and the second direction are orthogonal to each other;
the normalized apparent resistivity determining module is used for obtaining a first direction apparent resistivity and a second direction apparent resistivity according to the first direction potential difference and the second direction potential difference, and obtaining a normalized apparent resistivity according to the first direction apparent resistivity and the second direction apparent resistivity;
the background potential difference calculation module is used for calculating a background first-direction potential difference and a background second-direction potential difference between different measuring points according to a preset underground medium background resistivity model;
the normalized background apparent resistivity determining module is used for obtaining background first direction apparent resistivity and background second direction apparent resistivity according to the background first direction potential difference and the background second direction potential difference, and obtaining normalized background apparent resistivity according to the background first direction apparent resistivity and the background second direction apparent resistivity;
and the resistivity anisotropy coefficient determination module is used for obtaining the resistivity anisotropy coefficient of the underground medium of the measuring area according to the normalized apparent resistivity and the normalized background apparent resistivity.
9. A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor implements the steps of the method of any one of claims 1 to 7 when executing the computer program.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111157020.6A CN113917541B (en) | 2021-09-30 | 2021-09-30 | Method and device for acquiring electric anisotropy of underground medium based on direct current electric exploration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111157020.6A CN113917541B (en) | 2021-09-30 | 2021-09-30 | Method and device for acquiring electric anisotropy of underground medium based on direct current electric exploration |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113917541A true CN113917541A (en) | 2022-01-11 |
| CN113917541B CN113917541B (en) | 2024-02-23 |
Family
ID=79237506
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111157020.6A Active CN113917541B (en) | 2021-09-30 | 2021-09-30 | Method and device for acquiring electric anisotropy of underground medium based on direct current electric exploration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113917541B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116088059A (en) * | 2023-04-10 | 2023-05-09 | 山东省煤田地质规划勘察研究院 | Artificial source frequency domain electromagnetic exploration method and system based on double-source error frequency emission |
| CN118392020A (en) * | 2024-05-07 | 2024-07-26 | 中煤江南建设发展集团有限公司 | Direct-insert sediment thickness electronic calculation method, system and equipment |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030197510A1 (en) * | 2002-04-10 | 2003-10-23 | Halliburton Energy Services, Inc. | Galvanic method of measuring electrical anisotropy |
| CN101194262A (en) * | 2005-06-09 | 2008-06-04 | 埃克森美孚上游研究公司 | Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys |
| US20110156709A1 (en) * | 2009-12-30 | 2011-06-30 | Smith International, Inc. | Microresistivity anisotropy logging tool employing a monopole current injection electrode |
| CN107742015A (en) * | 2017-09-30 | 2018-02-27 | 中南大学 | DC lasering electric method Three-dimensional Numerical Simulation Method based on any dipole dipole device |
| CN110187394A (en) * | 2019-06-20 | 2019-08-30 | 甘肃省地震局（中国地震局兰州地震研究所） | Double field source electromagnetic depth methods obtain the anisotropic method and device of formation resistivity |
| CN111551997A (en) * | 2020-03-12 | 2020-08-18 | 上海环联生态科技有限公司 | Investigation system and investigation method for concealed fracture layer |
| CN112709567A (en) * | 2019-10-25 | 2021-04-27 | 斯伦贝谢技术有限公司 | Determination of anisotropic subsurface properties using electromagnetic measurements |
| CN113238285A (en) * | 2021-05-08 | 2021-08-10 | 桂林理工大学 | Resistivity calculation method and system for geophysical charging method exploration and terminal |
-
2021
- 2021-09-30 CN CN202111157020.6A patent/CN113917541B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030197510A1 (en) * | 2002-04-10 | 2003-10-23 | Halliburton Energy Services, Inc. | Galvanic method of measuring electrical anisotropy |
| CN101194262A (en) * | 2005-06-09 | 2008-06-04 | 埃克森美孚上游研究公司 | Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys |
| US20110156709A1 (en) * | 2009-12-30 | 2011-06-30 | Smith International, Inc. | Microresistivity anisotropy logging tool employing a monopole current injection electrode |
| CN107742015A (en) * | 2017-09-30 | 2018-02-27 | 中南大学 | DC lasering electric method Three-dimensional Numerical Simulation Method based on any dipole dipole device |
| CN110187394A (en) * | 2019-06-20 | 2019-08-30 | 甘肃省地震局（中国地震局兰州地震研究所） | Double field source electromagnetic depth methods obtain the anisotropic method and device of formation resistivity |
| CN112709567A (en) * | 2019-10-25 | 2021-04-27 | 斯伦贝谢技术有限公司 | Determination of anisotropic subsurface properties using electromagnetic measurements |
| CN111551997A (en) * | 2020-03-12 | 2020-08-18 | 上海环联生态科技有限公司 | Investigation system and investigation method for concealed fracture layer |
| CN113238285A (en) * | 2021-05-08 | 2021-08-10 | 桂林理工大学 | Resistivity calculation method and system for geophysical charging method exploration and terminal |
Non-Patent Citations (1)
| Title |
|---|
| 徐宏宇: "普通电阻率测井复杂井眼环境校正模板计算正演方法研究", 《中国优秀硕士学位论文全文数据库基础科学辑》, no. 2018 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116088059A (en) * | 2023-04-10 | 2023-05-09 | 山东省煤田地质规划勘察研究院 | Artificial source frequency domain electromagnetic exploration method and system based on double-source error frequency emission |
| CN118392020A (en) * | 2024-05-07 | 2024-07-26 | 中煤江南建设发展集团有限公司 | Direct-insert sediment thickness electronic calculation method, system and equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113917541B (en) | 2024-02-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lei et al. | A numerical study of stress variability in heterogeneous fractured rocks | |
| CN107742015B (en) | DC induced polarization method three-dimensional numerical simulation method based on arbitrary dipole-dipole device | |
| Oldenburg et al. | Three dimensional inversion of multisource time domain electromagnetic data | |
| CN113917541B (en) | Method and device for acquiring electric anisotropy of underground medium based on direct current electric exploration | |
| Weiss | Finite-element analysis for model parameters distributed on a hierarchy of geometric simplices | |
| Xu et al. | A comparison between ES-MDA and restart EnKF for the purpose of the simultaneous identification of a contaminant source and hydraulic conductivity | |
| CN110346835B (en) | Magnetotelluric forward modeling method, forward modeling system, storage medium, and electronic device | |
| Mehanee | Tracing of paleo-shear zones by self-potential data inversion: case studies from the KTB, Rittsteig, and Grossensees graphite-bearing fault planes | |
| Ingeman-Nielsen et al. | Effect of electrode shape on grounding resistances—Part 1: The focus-one protocol | |
| Li et al. | Inverse calculation of in situ stress in rock mass using the surrogate-model accelerated random search algorithm | |
| Zhang et al. | Three-dimensional magnetotelluric modeling using the finite element model reduction algorithm | |
| CN115657142A (en) | Method, device, equipment and storage medium for detecting and monitoring water damage of mine stratum | |
| CN113076678B (en) | Frequency domain two-degree body weight abnormity rapid numerical simulation method and device | |
| CN110766794B (en) | Method and device for acquiring fault geometry and computer equipment | |
| Liu et al. | Performance investigations of auxiliary‐space Maxwell solver preconditioned iterative algorithm for controlled‐source electromagnetic induction problems with electrical anisotropy | |
| Sungkono | Robust interpretation of single and multiple self-potential anomalies via flower pollination algorithm | |
| Gao et al. | Fast 3D modeling of borehole induction measurements in dipping and anisotropic formations using a novel approximation technique | |
| Ai et al. | Stress solution of transversely isotropic surrounding rock for shallow buried elliptic tunnels | |
| CN113673163B (en) | Three-dimensional magnetic abnormal constant rapid forward modeling method, device and computer equipment | |
| Chen et al. | Probabilistic back analysis for geotechnical engineering based on Bayesian and support vector machine | |
| Hu et al. | A new method of 3D direct current resistivity modelling using a long electrode source for forward probing in tunnels | |
| Malecki et al. | A new method for absolute underground positioning based on transient electromagnetics | |
| CN113268702A (en) | Frequency domain magnetic gradient tensor transformation method and device and computer equipment | |
| CN114444363B (en) | Method, device and equipment for calculating response of high-precision receiver at any position | |
| Tang et al. | 2.5-D DC resistivity modeling considering flexibility and accuracy |
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 |