CN108873073B - 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
- CN108873073B CN108873073B CN201810343479.7A CN201810343479A CN108873073B CN 108873073 B CN108873073 B CN 108873073B CN 201810343479 A CN201810343479 A CN 201810343479A CN 108873073 B CN108873073 B CN 108873073B
- Authority
- CN
- China
- Prior art keywords
- electrode
- dimensional
- apparent resistivity
- space position
- resistivity
- 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
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000010291 electrical method Methods 0.000 title claims abstract description 17
- 238000003325 tomography Methods 0.000 title claims abstract description 12
- 238000005553 drilling Methods 0.000 claims abstract description 22
- 238000001514 detection method Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 3
- 238000011835 investigation Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
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, the three electrode space positions acquired every time classify 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, specific steps of the invention are as follows:
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, the three electrode space positions acquired every time classify 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, in which:
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 are as follows: set the three electrode apparent resistivity ρ three-dimensional space positions acquired every time as (x0,y0,z0)；
If in three electrodes currently acquiredThen directly give up, carries out the determination of next three electrode space position, ifThen carry out the division of following three types；
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 L in the face yoz0, its direction vector is
OO ",It was O'(0,0,0) straight line along the x axis and straight line L0Intersection 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 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
Are as follows:
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, the three electrode space positions acquired every time classify 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, in which:
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 position detailed process are as follows: set the three electrode apparent resistivity three-dimensional space positions acquired every time as (x0,y0,z0)；If currently adopting
In three electrodes of collectionThen directly give up, carries out the determination of next three electrode space position, ifThen into
The division of the following three types of row:
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 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, NIt is parallel to the straight line L in the face yoz0, its direction vector is OO ",It is
Cross O'(0,0,0) straight line along the x axis and straight line L0Intersection 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 CN108873073A (en) | 2018-11-23 |
| CN108873073B true 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) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112129632B (en) * | 2020-09-24 | 2021-11-09 | 东北大学 | Method for calibrating creep damage and fracture surface of rock by using high-density resistivity |
| CN112666614B (en) * | 2021-01-15 | 2022-09-06 | 西南科技大学 | Debris flow source static reserve calculation method based on electrical prospecting and digital elevation model |
| CN114166892B (en) * | 2021-12-09 | 2024-06-21 | 中国矿业大学 | Loaded rock sample damage self-potential imaging method based on network parallel electrical method |
| CN114994775B (en) * | 2022-08-08 | 2022-11-15 | 山东大学 | Double-line detection device, system and array acquisition method for cross-well induced polarization measurement |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7269515B2 (en) * | 2004-06-15 | 2007-09-11 | Baker Hughes Incorporated | Geosteering in anisotropic formations using multicomponent induction measurements |
| CN104407394B (en) * | 2014-11-26 | 2015-07-22 | 山东大学 | 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 |
| CN106125147B (en) * | 2016-08-18 | 2018-07-10 | 安徽惠洲地质安全研究院股份有限公司 | Orientation electric logging device and its measuring method based on electrical method of network concurrency technology |
-
2018
- 2018-04-17 CN CN201810343479.7A patent/CN108873073B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108873073A (en) | 2018-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108873073B (en) | A kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency | |
| 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 | |
| CN106772644A (en) | mine transient electromagnetic three-component detection method | |
| WO2021243967A1 (en) | Three-dimensional resistivity tomography method and system | |
| CN107742015A (en) | DC lasering electric method Three-dimensional Numerical Simulation Method based on any dipole dipole device | |
| CN103064124B (en) | A kind of ratio approach correcting electromagnetic survey Considering Terrain Effect | |
| Szalai et al. | Geoelectric mapping of near-surface karstic fractures by using null arrays | |
| Cardimona | zywvutsrqponmlihgfed | |
| 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 | |
| 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 | |
| CN112630842B (en) | Dipole differential device and method applied to direct-current resistivity method exploration | |
| 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 | |
| Jiang et al. | Application of three-dimensional electrical resistivity tomography in urban zones by arbitrary electrode distribution survey design | |
| 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 | |
| CN110058318A (en) | A kind of one-dimensional square bearing calibration of the natural electric field based on measuring electrode coordinate | |
| CN114152989B (en) | Method for calculating coordinates of pseudo-section diagram of two-dimensional asymmetric quadrupole electrical method | |
| CN111812728B (en) | Well ground resistivity CT observation system and working method thereof | |
| Monahan | Investigating fault structure using electrical resistivity tomography |
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 |