CN108873073A - A kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency - Google Patents
A kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency Download PDFInfo
- Publication number
- CN108873073A CN108873073A CN201810343479.7A CN201810343479A CN108873073A CN 108873073 A CN108873073 A CN 108873073A CN 201810343479 A CN201810343479 A CN 201810343479A CN 108873073 A CN108873073 A CN 108873073A
- Authority
- CN
- China
- Prior art keywords
- electrode
- dimensional
- apparent resistivity
- resistivity
- coordinate
- 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
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
Abstract
Across the hole resistivity tomography method of three-dimensional based on electrical method of network concurrency that the invention discloses a kind of, first then drilling carries out electrode Hydrographic General Line；And then multiple three electrodes apparent resistivity is carried out to survey line top electrode using electrical method of network concurrency and is acquired, then classify to the position of three different electrodes, according to different types of three electrode space position determine its corresponding apparent resistivity locating for three-dimensional space position calculation method, and then obtain three-dimensional space position locating for different types of apparent resistivity, ultimately produce the apparent resistivity distribution map that three-dimensional bodies are surrounded in multiple drillings, it can reflect the spatial distribution characteristic of subsurface anomaly body, according to distribution map convenient for providing data supporting when subsequent underground space development.
Description
Technical field
It is specifically a kind of based on electrical method of network concurrency the present invention relates to a kind of underground dc resistivity imaging detection method
Across the hole resistivity tomography method of three-dimensional.
Background technique
3 D resistivity tomography detection method is one of dc resistivity, be detection electrode is put into it is multiple
Detection method in drilling (three or three or more) in (or in well) a kind of hole of acquisition electric current and electric potential signal.This method pair
Drill enclosing region internal abnormality body detection accuracy and resolution ratio be substantially better than customary DC resistivity method, mineral products, the energy,
It is used widely in the fields such as hydrogeological disaster exploration.With the fast development of Urban Underground Space Excavation, using three-dimensional across
Hole Resistivity Tomography method is fine to the carry out of the geological disasters such as boulder, the solution cavity encountered in the processs of construction such as subway
It detects, has important practical significance to urban underground space security development.
Three-dimensional across hole Resistivity Tomography method mainly uses finite difference or limited when carrying out into figure at present
The three-dimensional-electroded method method of member.For it is multiple drilling (three or three or more) encirclement three-dimensional bodies apparent resistivity at
Drawing method is there is not yet research.
Summary of the invention
Across the hole electricity of three-dimensional that in view of the above existing problems in the prior art, the present invention provides a kind of based on electrical method of network concurrency
Resistance rate chromatography imaging method is scattered in using the apparent resistivity that electrical method of network concurrency surrounds three-dimensional bodies for multiple drillings
Figure, and then according to the spatial distribution characteristic for reflecting subsurface anomaly body at figure, convenient for providing data when subsequent underground space development
Support.
To achieve the goals above, the technical solution adopted by the present invention is that：It is a kind of based on electrical method of network concurrency it is three-dimensional across
Hole resistivity tomography method, the specific steps are that：
A, at least three vertical drillings are drilled on ground, and the subpoint of all vertical drillings on the ground is same
On round circumference；
B, electrode survey line is laid in each vertical drilling；
C, it carries out multiple three electrodes apparent resistivity to survey line top electrode with electrical method of network concurrency to acquire, three electrode is by one
A current electrode and two measuring electrode compositions；
D, classify according to the three electrode space positions acquired every time and determine its type；
E, according to different types of three electrode space position determine its corresponding apparent resistivity locating for three-dimensional space position
Calculation method, and then obtain three-dimensional space position locating for different types of apparent resistivity；
F, step D to E is repeated, is stopped after obtaining the apparent resistivity three-dimensional space position of all three electrodes of acquisition；
G, by the apparent resistivity three-dimensional space position of all three electrodes obtained, underground three is obtained after handling using interpolation method
Tie up apparent resistivity spatial distribution map；
H, it is finally completed according to apparent resistivity space distribution situation between detection architectonic drilling.
Compared with prior art, the present invention first acquires the apparent resistivity and its seat of multiple three electrodes using electrical method of network concurrency
Mark, then classifies to the position of three different electrodes, determines that it is corresponding according to different types of three electrode space position
The calculation method of three-dimensional space position locating for apparent resistivity, and then obtain three-dimensional space meta position locating for different types of apparent resistivity
It sets, ultimately produces multiple drillings and surround the apparent resistivity distribution map of three-dimensional bodies, and then reflect that underground is different according to distribution map
The spatial distribution characteristic of normal body, convenient for providing data supporting when subsequent underground space development.
Detailed description of the invention
Fig. 1 is the three-dimensional simulation schematic diagram of multiple drilling surveys line in the present invention；
Fig. 2 is the three-dimensional coordinate figure established using the midpoint Fig. 1 O' as origin in the present invention；
Fig. 3 is the schematic diagram of the three-dimensional electrical method model of Fig. 1；
Fig. 4 is the y-coordinate of A, M, N and the identical schematic diagram of z coordinate in the present invention.
In figure：ρ is apparent resistivity.
Specific embodiment
The invention will be further described below.
As shown, it is of the invention the specific steps are：
A, at least three vertical drillings are drilled on ground, and the subpoint of all vertical drillings on the ground is same
On round circumference；
B, electrode survey line is laid in each vertical drilling；
C, multiple three electrodes apparent resistivity is carried out to each vertical drilling with electrical method of network concurrency to acquire, three electrode by
One current electrode and two measuring electrode compositions；
D, classify according to the three electrode space positions acquired every time and determine its type；
E, according to different types of three electrode space position determine its corresponding apparent resistivity locating for three-dimensional space position
Calculation method, and then obtain three-dimensional space position locating for different types of apparent resistivity；
F, step D to E is repeated, is stopped after obtaining the apparent resistivity three-dimensional space position of all three electrodes of acquisition；
G, by the apparent resistivity three-dimensional space position of all three electrodes obtained, underground three is obtained after handling using interpolation method
Tie up apparent resistivity spatial distribution map；
H, it is finally completed according to apparent resistivity space distribution situation between detection architectonic drilling.
Further, if the radius of vertical drilling is R, the bottom surface center of circle is O' point, and the length of electrode survey line is Xi, each electrode
Three-dimensional coordinate be (x, y z)；Wherein [0, Xi] x ∈, y ∈ [- R, R], z ∈ [- R, R].
Further, if it is respectively M and N that current electrode, which is A and two measuring electrode, in three electrodes, coordinate is respectively A
(x1, y1, z1), M (x2, y2, z2), N (x3, y3, z3)；It sets up an officeIt is the midpoint of measuring electrode M and N；
If investigation depth h=h0OA, depth coefficient h0Value range 0.2~1, wherein：
Further, the type of the three electrode space positions acquired every time is determining and obtains locating for different types of apparent resistivity
The detailed process of three-dimensional space position is：If the three electrode apparent resistivity ρ three-dimensional space positions acquired every time are (x0,y0,z0)；
If in three electrodes currently acquiredThe division of following three types is then carried out, ifThen directly give up, carries out
The determination of next three electrode space position：
If a, x1=x2=x3, then judge:
1. if (x1- h) ∈ [0, Xi], then
2. if (x1+ h) ∈ [0, Xi], then
3. if (x1- h) and (x1+ h) in section [0, Xi], then coordinateWithIt is apparent resistivity ρ three-dimensional space position (x0,y0,z0) coordinate；
If b, the y-coordinate of A, M, N are identical with z coordinate, i.e. y1=y2=y3And z1=z2=z3, A, M, N are in same at this time
On electrode survey line；Cross the midpoint M, NIt is parallel to the straight line L0 in the face yoz, its direction vector is
OO ",It was O'(0,0,0) intersection point of straight line along the x axis and straight line L0 (OO ")；
So as to obtain three-dimensional space position：
If c, electrode A, M, N spatial coordinate location do not meet the condition of a and b, apparent resistance is determined using following process
Three-dimensional space position (x locating for rate ρ0,y0,z0)；
If the direction vector of straight line L is (a, b, c), it is known that straight line L is perpendicular to face AMN.
AM vectorial coordinate is (x1-x2,y1-y2,z1-z2), the vectorial coordinate of AN is (x1-x3,y1-y3,z1-z3), MN vector
Coordinate is (x2-x3,y2-y3,z2-z3)
1. if x1≠x3, then
B=1
2. if x1=x3, then
B=1
According to
?：
Or：
In the result obtained, if x0∈ [0, Xi], then point (x0,y0,z0) in exploration volume, determine that it is apparent resistance
The three-dimensional space position of rate ρ, otherwise gives up.
Claims (4)
1. a kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency, which is characterized in that its specific steps
For：
A, at least three vertical drillings are drilled on ground, and the subpoint of all vertical drillings on the ground is in the same circle
On circumference；
B, electrode survey line is laid in each vertical drilling；
C, it carries out multiple three electrodes apparent resistivity to survey line top electrode with electrical method of network concurrency to acquire, three electrode is supplied by one
Electrode and two measuring electrode compositions；
D, classify according to the three electrode space positions acquired every time and determine its type；
E, according to different types of three electrode space position determine its corresponding apparent resistivity locating for three-dimensional space position calculating
Method, and then obtain three-dimensional space position locating for different types of apparent resistivity；
F, step D to E is repeated, is stopped after obtaining the apparent resistivity three-dimensional space position of all three electrodes of acquisition；
G, by the apparent resistivity three-dimensional space position of all three electrodes obtained, subsurface three-dimensional view is obtained after handling using interpolation method
Resistivity spatial distribution map；
H, it is finally completed according to apparent resistivity space distribution situation between detection architectonic drilling.
2. a kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency according to claim 1,
It is characterized in that, if the radius of vertical drilling is R, the bottom surface center of circle is O' point, and the length of electrode survey line is Xi, the three-dimensional of each electrode
Coordinate is (x, y, z)；Wherein [0, Xi] x ∈, y ∈ [- R, R], z ∈ [- R, R].
3. a kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency according to claim 2,
It is characterized in that, if it is respectively M and N that current electrode, which is A and two measuring electrode, in three electrodes, coordinate is respectively A (x1, y1,
z1), M (x2, y2, z2), N (x3, y3, z3)；It sets up an officeIt is the midpoint of measuring electrode M and N；
If investigation depth h=h0OA, depth coefficient h0Value range 0.2~1, wherein：
4. a kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency according to claim 3,
It is characterized in that, the type of the three electrode space positions acquired every time is determining and obtains three-dimensional space locating for different types of apparent resistivity
Between the detailed process of position be：If the three electrode apparent resistivity three-dimensional space positions acquired every time are (x0,y0,z0)；If currently adopting
In three electrodes of collectionThe division of following three types is then carried out, ifThen directly give up, carries out next
The determination of three electrode space positions：
If a, x1=x2=x3, then judge:
1. if (x1- h) ∈ [0, Xi], then
2. if (x1+ h) ∈ [0, Xi], then
3. if (x1- h) and (x1+ h) in section [0, Xi], then coordinateWith
It is three-dimensional space position (x0,y0,z0) coordinate；
If b, the y-coordinate of A, M, N are identical with z coordinate, i.e. y1=y2=y3And z1=z2=z3, cross the midpoint M, N
It is parallel to the straight line L0 in the face yoz, its direction vector is OO ",It was O'(0,0,0) along the x axis
Straight line and straight line L0 (OO ") intersection point；
So as to obtain three-dimensional space position：
If c, electrode A, M, N spatial coordinate location do not meet the condition of a and b, apparent resistivity institute is determined using following process
Three-dimensional space position (the x at place0,y0,z0)；
If the direction vector of straight line L is (a, b, c), it is known that straight line L is perpendicular to face AMN.
AM vectorial coordinate is (x1-x2,y1-y2,z1-z2), the vectorial coordinate of AN is (x1-x3,y1-y3,z1-z3), MN vectorial coordinate
For (x2-x3,y2-y3,z2-z3)
1. if x1≠x3, then
B=1
2. if x1=x3, then
B=1
According to
?：
Or：
In the result obtained, if x0∈ [0, Xi], then point (x0,y0,z0) in exploration volume, determine that it is the three of apparent resistivity
Dimension space position, otherwise gives up.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810343479.7A CN108873073B (en) | 2018-04-17 | 2018-04-17 | A kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810343479.7A CN108873073B (en) | 2018-04-17 | 2018-04-17 | A kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108873073A true CN108873073A (en) | 2018-11-23 |
| CN108873073B CN108873073B (en) | 2019-06-28 |
Family
ID=64326938
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810343479.7A Active CN108873073B (en) | 2018-04-17 | 2018-04-17 | A kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108873073B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112129632A (en) * | 2020-09-24 | 2020-12-25 | 东北大学 | Method for calibrating creep damage and fracture surface of rock by using high-density resistivity |
| CN112666614A (en) * | 2021-01-15 | 2021-04-16 | 西南科技大学 | Debris flow source static reserve calculation method based on electrical prospecting and digital elevation model |
| CN114994775A (en) * | 2022-08-08 | 2022-09-02 | 山东大学 | Double-line detection device, system and array acquisition method for cross-well induced polarization measurement |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070210799A1 (en) * | 2004-06-15 | 2007-09-13 | Baker Hughes Incorporated | Geosteering In Earth Formations Using Multicomponent Induction Measurements |
| CN104407394A (en) * | 2014-11-26 | 2015-03-11 | 山东大学 | Subway shield section boulder refining detecting method based on resistivity trans-hole CT (computed tomography) |
| CN104678443A (en) * | 2015-01-22 | 2015-06-03 | 深圳市市政设计研究院有限公司 | Lonestone detector and control method of boulder detector |
| CN106125147A (en) * | 2016-08-18 | 2016-11-16 | 安徽惠洲地质安全研究院股份有限公司 | Orientation based on electrical method of network concurrency technology electric logging device and measuring method thereof |
-
2018
- 2018-04-17 CN CN201810343479.7A patent/CN108873073B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070210799A1 (en) * | 2004-06-15 | 2007-09-13 | Baker Hughes Incorporated | Geosteering In Earth Formations Using Multicomponent Induction Measurements |
| CN104407394A (en) * | 2014-11-26 | 2015-03-11 | 山东大学 | Subway shield section boulder refining detecting method based on resistivity trans-hole CT (computed tomography) |
| CN104678443A (en) * | 2015-01-22 | 2015-06-03 | 深圳市市政设计研究院有限公司 | Lonestone detector and control method of boulder detector |
| CN106125147A (en) * | 2016-08-18 | 2016-11-16 | 安徽惠洲地质安全研究院股份有限公司 | Orientation based on electrical method of network concurrency technology electric logging device and measuring method thereof |
Non-Patent Citations (1)
| Title |
|---|
| 刘征宇 等: "用于地铁盾构区间孤石探测的三维电阻率跨孔CT方法研究及物理模型试验", 《岩石力学与工程学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112129632A (en) * | 2020-09-24 | 2020-12-25 | 东北大学 | Method for calibrating creep damage and fracture surface of rock by using high-density resistivity |
| CN112129632B (en) * | 2020-09-24 | 2021-11-09 | 东北大学 | Method for calibrating creep damage and fracture surface of rock by using high-density resistivity |
| CN112666614A (en) * | 2021-01-15 | 2021-04-16 | 西南科技大学 | Debris flow source static reserve calculation method based on electrical prospecting and digital elevation model |
| CN114994775A (en) * | 2022-08-08 | 2022-09-02 | 山东大学 | Double-line detection device, system and array acquisition method for cross-well induced polarization measurement |
| CN114994775B (en) * | 2022-08-08 | 2022-11-15 | 山东大学 | Double-line detection device, system and array acquisition method for cross-well induced polarization measurement |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108873073B (en) | 2019-06-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Okpoli | Sensitivity and resolution capacity of electrode configurations | |
| Chambers et al. | River terrace sand and gravel deposit reserve estimation using three-dimensional electrical resistivity tomography for bedrock surface detection | |
| CN108873073B (en) | A kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency | |
| WO2021243967A1 (en) | Three-dimensional resistivity tomography method and system | |
| CN103064124B (en) | A kind of ratio approach correcting electromagnetic survey Considering Terrain Effect | |
| Cardimona | zywvutsrqponmlihgfed | |
| Szalai et al. | Geoelectric mapping of near-surface karstic fractures by using null arrays | |
| AU2012363900B2 (en) | Focused array laterolog tool | |
| CN109100808A (en) | A kind of horizontal magnetic polarization field detection method of multi-thread source transient electromagnetic | |
| CN104330832A (en) | Transient electromagnetic quick three-dimensional human-machine interaction inversion method | |
| CN110208867A (en) | A kind of three-dimensional electrical prospecting method based on composite profile | |
| CN103278855A (en) | Method for eliminating influence of roadways and terrains on apparent resistivity in direct-current exploration | |
| CN107065019A (en) | Applied to road disaster and the 3 D electromagnetic imaging device and application method that collapse detection | |
| CN109343131B (en) | Lateral power supply ground well direct current method detection method and device | |
| CN105842738B (en) | A kind of buried target body quantitatively judges method and apparatus extremely | |
| CN110989001A (en) | Pile foundation three-dimensional drilling CT detection method applying HDD technology | |
| CN112630842B (en) | Dipole differential device and method applied to direct-current resistivity method exploration | |
| Kirsch et al. | Geoelectrical methods | |
| CN113238285A (en) | Resistivity calculation method and system for geophysical charging method exploration and terminal | |
| CN112882124A (en) | Method for delineating beneficial parts of three-dimensional uranium mineralization | |
| CN205958772U (en) | Position electric logging device based on parallel electric method technique of network | |
| Zhang et al. | A new method of anti-collision while drilling based on radial gradient measurement | |
| Domenzain et al. | 3D DC inversion, visualization, and processing of dense time-lapse data in fine domains applied to remediation monitoring | |
| KR20150021623A (en) | Geological resource monitoring method using electrical resistivity | |
| CN114152989B (en) | Method for calculating coordinates of pseudo-section diagram of two-dimensional asymmetric quadrupole electrical method |
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 |