CN102508303A - Advanced detection method for focusing chromatography induced polarization of underground engineering - Google Patents
Advanced detection method for focusing chromatography induced polarization of underground engineering Download PDFInfo
- Publication number
- CN102508303A CN102508303A CN2011103766381A CN201110376638A CN102508303A CN 102508303 A CN102508303 A CN 102508303A CN 2011103766381 A CN2011103766381 A CN 2011103766381A CN 201110376638 A CN201110376638 A CN 201110376638A CN 102508303 A CN102508303 A CN 102508303A
- Authority
- CN
- China
- Prior art keywords
- electrode
- survey line
- workplace
- electrodes
- data
- 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
Images
Classifications
-
- 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
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
The invention relates to an advanced detection method for focusing chromatography induced polarization of underground engineering. The advanced detection method comprises the following steps of: arranging a circle of shield electrodes A1 on a contour line of an operating surface of the underground engineering by utilizing a focusing operating manner and arranging a plurality of parallel measuring lines on parts with different heights from a base plate to a top plate firstly, wherein the measuring lines are provided with power supply electrodes and receiving electrodes, and the base plate is also provided with power supply electrodes and receiving electrodes; acquiring data layer by layer by utilizing a chromatography method; carrying out inversion on measuring line potential difference data obtained through measurement to obtain a three-dimensional resistivity image of a geologic body in the front of the operating surface; drawing a relation coordinate graph of difference data and horizontal distances of half time by using the half time obtained through measurement, calculating an envelope area St of half time difference and abscissa axis, and drawing a two-dimensional profile map of the envelope area St to realize estimation of the quantity of water in the front of the operating surface. The advanced detection method for the focusing chromatography induced polarization of the underground engineering, disclosed by the invention, has excellent advanced detection directivity, realizes three-dimensional location of an aquifer and solves problem of incapability of location in a traditional focusing electric method during single-point measurement.
Description
Technical field
The present invention relates to unfavorable geology geology forward probe method in the underground engineering construction processes such as a kind of tunnel or tunnel.
Background technology
At present, China has become in the world underground works such as tunnel and has built largestly, builds the highest country of difficulty.The geologic hazard that the underground engineering construction phase takes place frequently becomes the huge challenge that underground engineering construction faces, and gushing water, prominent mud disaster are all occuping the prostatitis aspect frequency and the death toll.It is significant in the underground engineering construction advanced prediction to be implemented in gushing water water burst disaster source.For moisture architectonic forward probe, main task is divided into two: 1. need verify the characteristics such as position, scale or even form of water-bearing structure, promptly need realize the three-dimensional imaging of water-bearing structure; 2. need predict the water yield in the water-bearing structure, this is to formulate prominent water burst prophylactico-therapeutic measures and the key factor of estimating tunnel safety the constructing tunnel phase.
Induced polarization method is to survey the resistivity prospecting method of geological condition according to the effect of induced polarization of rock, and present induced polarization method forecasting procedure can be divided into three electrode method non-focusing type and single-point focus type.Wherein, Near the side that three electrode method non-focusing type forward probe method can't shield the survey line disturbs (like low-resistance water content, hardware etc.) unusually; Than the useful information that is difficult to go out development end the place ahead under the complex environment from the background interference extracting data; Cause forecast precision to reduce easily even wrong report, this problem not can solve always, has seriously restricted the development of three electrode method non-focusing type electrical method class forecasting technique.Ruan Baiyao etc. have proposed a kind of tunnel focus dc resistivity method forward probe method, and BEAM (the Bore-Tunneling Electrical Ahead Monitoring) technology of Germany's research and development is unique a kind of focusing class induced polarization methods that the moisture situation in the place ahead, tunnel is surveyed that are used for.
Patent CN1904644A discloses a kind of tunnel real-time continuous advanced prediction method and device of networking; Step is following: (1) according to the site operation concrete condition, adopts that the tunnel is outer, the metalwork at tunnel inwall or workplace scene arranges infinite distance potential electrode B, guarded electrode A respectively
1With central electrode A
0, and link to each other with on-the-spot main frame with cable; (2) on-the-spot main frame is through guarded electrode A
1With central electrode A
0The loop that forms with infinite distance potential electrode B receives the electric signal of central electrode simultaneously with the pulse current of time domain or the different voltage same polarities of frequency field mode of operation synchronized transmissions, treatedly converts digital signal into A/D, is input to host cpu; (3) on-the-spot host cpu is calculated the apparent resistivity of detection of a target body and is looked the induced polarization rate; (4) according to the apparent resistivity that calculates in the step (3) with look the induced polarization rate and compare with the numerical range of looking the induced polarization rate with the apparent resistivity of the plurality of target body that presets, with apparent resistivity with look the scope area that the induced polarization rate falls into jointly and come recognition objective body object and visualization display and warning.
But still there is following problem in present focusing class electrical method: the 1. present existing thinking that type electrical survey (-ing) method has all been used for reference oil well logging that focuses on; Belong to the single-point focus type; Single potential electrode just only is set on test surface, does not adopt chromatography to survey, can't obtain the three-dimensional information of geologic body; 2. utilize the curve of each measurement result and tunnel mileage to infer the moisture situation in development end the place ahead, detection range is less, focuses on a type electrical survey (-ing) technology at present and does not get into the practical stage as yet.3. present electrical method class forecasting procedure all can not realize the estimation to the water content water yield, has become technical barrier of domestic and international Tunnel Engineering circle based on the estimation of the water yield of geophysical exploration technology.
Summary of the invention
The objective of the invention is in order to overcome the deficiency of above-mentioned prior art; A kind of underground works chromatographic focusing induced polarization forward probe method is provided, and this method is estimated gushing water, prominent mud disaster source water yield to the three-dimensional imaging of workplace front geological condition chromatography; And do not need continuity to survey, efficient is high.
The technical scheme that the present invention adopts is:
Underground works chromatographic focusing induced polarization forward probe method: at first adopt the working method that focuses on, on underground works workplace outline line, arrange a circle guarded electrode A
1,, establish electrodes of A on the survey line arranging many parallel surveys line to the differing heights at top apart from base plate
0With receiving electrode M
1, M
2, also be provided with transmitting electrode B and receiving electrode N on the base plate; The method that adopts chromatography then is image data successively, with electrodes of A
0With guarded electrode A
1Infeed electric current in the same way, test survey line receiving electrode M respectively
1With the electric potential difference of base plate receiving electrode N with when partly declining, and change the survey line receiving electrode M that size of current test different electric flows down
1With the electric potential difference of base plate receiving electrode N with when partly declining; Survey line electric potential difference data with recording are carried out inverting, obtain the 3 D resistivity image of workplace front geological body, realize the three-dimensional localization of water content; With record partly decline the time draw difference data and horizontal range when partly declining concern coordinate diagram, the difference when calculating partly declines and the envelope area S of abscissa axis
t, draw out envelope area S
tTwo-dimensional cross section, realize the estimation of workplace the place ahead water yield.
Described underground works chromatographic focusing induced polarization forward probe method preferably includes following steps:
(1) adopts the working method that focuses on, at first on (tunnel etc.) underground works workplace outline line, arrange a circle equally spaced (preferred 4-8) guarded electrode A
1, apart from base plate 1/3H, two parallel survey line L of 2/3H height placement
1, L
2, every survey line is arranged electrode number N
1=N
2=15, adjacent two distance between electrodes are d, and wherein H is the height of vault apart from base plate; Electrodes of A
0With receiving electrode M
1, M
2Place survey line L
1, L
2(receiving electrode M
1, M
2Adjustable positions and with survey line on the electrode position arranged corresponding, be in order to guarantee M
1M
2Between electric potential difference etc. zero, guarantee focusing effect, measure the electric potential difference between M1 and the N simultaneously), another transmitting electrode B utmost point places on the base plate of workplace rear apart from workplace 100-150m, receiving electrode N places on the base plate;
(2) method that adopts chromatography is image data successively,
1. survey line L
1Electrodes of A
0Infeed positive current I
0, transmitting electrode B infeeds negative current-I
0, guarded electrode A
1Infeed and electrodes of A
0Positive current in the same way keeps electrodes of A
0With potential electrode M
1, M
2Spacing A
0M
1=M
1M
2=nd, wherein d is adjacent two distance between electrodes, n is natural number and 1≤n≤7, moves M
1, M
2The other end from survey line one end to survey line; M
1, M
2Whenever satisfy A
0M
1=M
1M
2=nd moves once, electrode A
0, M
1And M
2Every measurement 1 time is regulated supply current and is made M
1M
2Between electric potential difference U
M1M2=0, and note M
1Electric potential difference U1 between the N
M1NT when partly declining;
2. survey line L
2Repeat sequence L
1Measuring process 1.;
3. electrode A
0Infeed positive current 2I
0, electrode B infeeds negative current-I
0, repeating step 1., 2.;
(3) the 3-d inversion imaging of carrying out water-bearing structure is located, and adopts the processing such as least square inverting iteration based on smooth constraint, with the survey line L that records in the step (2)
1, L
2Electric potential difference data U1
M1NCarry out inverting, obtain the 3 D resistivity image of workplace front geological body, realize the three-dimensional localization of water content; Concrete steps are following:
3 D resistivity inversion equation with smooth constraint.
(A
TA+λC
TC)Δm=A
TΔd (1)
m
(k+1)=m
(k)+Δm (2)
In the formula, Δ d is the difference of observation data, and A is a partial derivative matrix, and Δ m is the model parameter correction of grid, and C is the smoothness matrix, and λ is Lagrangian constant, has determined the weight of smooth constraint, m
(k)Be the k time model parameter vector in the iteration, m
(k+1)Be the model parameter vector in the inverting iteration next time.
When carrying out inversion imaging (prior art), at first setting model parameter initial value is just being drilled calculating, and the observation data of the theory of computation compares judgement with data measured, if variance less than setting value, output model parameter then; Otherwise, calculate partial derivative matrix A, the smoothness Matrix C is utilized equation (2) can obtain the model parameter increment in each inverting iteration, thereby is tried to achieve the model parameter m in the next iteration
(k+1), until the variance of the theoretical observation data of just drilling and measurement data less than setting value (precision according to inverting that the user requires is different, and setting value is different, and inversion accuracy is high, and then setting value is less, otherwise setting value is bigger), output model parameter.Thereby inverting obtains the 3 D resistivity image of workplace front geological body, realizes the three-dimensional localization of water content;
Data were carried out the water yield estimation of workplace the place ahead when (4) employing partly declined; With record in the step (2) partly decline the time T data; When decline in two and half of size two primary currents when power supply identical point position data do poor; And draw the coordinate diagram that concerns of difference data and horizontal range when partly declining respectively, calculate difference and the envelope area S of abscissa axis when partly declining
t, and envelope area S
tWith the linear positive correlation of big or small V of the water yield, be the response of the water yield on the occasion of part, by the size on the occasion of the area estimation water yield of part, nearly again one to go on foot horizontal direction with workplace be that the vertical direction of X axle, workplace is the Y axle, draws out envelope area S
tTwo-dimensional cross section, realize the estimation of workplace the place ahead water yield.
The present invention focuses on induced polarization method, adopts the guarded electrode A that is arranged in the workplace outline line
1, make survey line L
1, L
2Last electrodes of A
0Electric current receive electric field extruding and flow into to workplace the place ahead, strengthened the directive property of electric current and improved investigation depth, solved the problem of classic method directive property difference; The chromatography induced polarization method; Adopt a plurality of parallel surveys line to obtain lot of data; Many survey line data are carried out inverting; Realized the three-dimensionalreconstruction of workplace front geological body electrical structure, realized three-dimensional localization, solved the spot measurement of traditional focus class methods and a difficult problem that can't three-dimensional localization underground engineering gushing water, prominent mud disaster source; Gushing water, the water yield estimation of prominent mud disaster source, the difference and the positive correlation of coordinate axis envelope area and the water-bearing structure water yield when partly declining, thus realization is to the estimation of underground engineering gushing water, the mud disaster source water yield of dashing forward.
The invention has the beneficial effects as follows:
The present invention proposes the advanced prediction method in a kind of underground works gushing water, prominent mud disaster source based on the focused resistivity imaging; Through at workplace outline line arranging shielding electrode; In shielding area, arrange two parallel surveys line, adopt least square inversion method, realized the three-dimensionalreconstruction of workplace front geological body electrical structure based on smooth constraint; Especially realize the water yield estimation of water-bearing structure, had following characteristic:
1. adopt the focusing-detection pattern, had good forward probe directive property, solved the difficult problem of previous methods side serious interference;
2. the method that chromatographic focusing is surveyed has been proposed; Can obtain the mass data of tunnel front geological information, and distribute, realize the three-dimensional localization of water content through the 3 D resistivity that Inversion Calculation obtains workplace the place ahead; Solved and focused on a type electrical method spot measurement, the difficult problem that can't locate in the past;
3. the difference data in the time of can obtaining induced polarization and partly decline has been drawn difference and the two dimensional cross-section of coordinate axis envelope area when partly declining, has realized the water yield estimation of water-bearing structure, has solved the crucial difficult problem of water content water yield prediction in the underground engineering construction;
4. need not follow the excavation continuous probe of workplace, improve efficient, practice thrift the advanced prediction cost, improve operating efficiency.
Description of drawings
Fig. 1 is a detection method process flow diagram of the present invention.
Fig. 2 a is the planimetric map of measurement mechanism of the present invention,
Fig. 2 b is the elevation drawing of measurement mechanism of the present invention,
Fig. 2 c is an observed pattern synoptic diagram of the present invention.
Fig. 3 is an inversion method flow process of the present invention.
Fig. 4 is a focused resistivity 3-d inversion image of the present invention.
Fig. 5 is the difference data figure of the present invention when partly declining.
Fig. 6 is difference and the two-dimensional cross section of coordinate axis envelope area of the present invention when partly declining.
Embodiment
Below in conjunction with accompanying drawing and embodiment the present invention is further specified.
Among detection method of the present invention such as Fig. 1; Flow process comprises: focus on working method and lay, chromatography image data, flow processs such as 3-d inversion location and water yield estimation; Has good forward probe directive property; Realize the poorly three-dimensional localization of plastid and the estimation of the water yield such as water-bearing structure at last, solved and focused on a type electrical method spot measurement, the difficult problem that can't locate in the past.
Embodiment 1:
With certain tunnel is example, and tunnel cross-section is an arch, and the high 10m of span 12m water burst occurs near the mileage YK123+672, and possibly there is water-bearing structure in the place ahead.In order to verify the geological condition in the place ahead,, adopt the chromatographic focusing induced polarization method to carry out forward probe for construction provides foundation and guidance.
(1) concrete working method is shown in Fig. 2 a, b, c.On the tunnel tunnel face outline line, arrange equally spaced 8 the guarded electrode A of a circle
1, apart from tunnel floor 2m, two parallel survey line L of 6m height placement
1, L
2, every survey line is arranged electrode number N
1=N
2=15; Electrodes of A
0With receiving electrode M
1, M
2Place survey line L
1, L
2, another transmitting electrode B utmost point places on the base plate of 150m, and receiving electrode N places on the base plate of 120m;
The idiographic flow of image data is: survey line L
1Electrodes of A
0Infeed positive current 0.1A, transmitting electrode B is a negative pole, guarded electrode A
1Infeed positive current 0.1A, keep electrodes of A
0With potential electrode M
1, M
2Spacing A
0M
1=M
1M
2=nd, M
1, M
2Successively from survey line L
1One end moves to the other end of survey line, and wherein d is adjacent two distance between electrodes, and n is natural number and 1≤n≤7; Electrode A
0, M
1And M
2Every measurement 1 time is regulated supply current and is made M
1M
2Between electric potential difference U
M1M2=0, and note M
1Electric potential difference U1 between the N
M1NT when partly declining; Measure survey line L according to above method
2, adjust electrode A then
0Infeed positive current 0.2A, repeat above surveying work.
(2) the electric potential difference data of gathering are carried out inversion imaging and explain, as shown in Figure 3, adopt processing such as least square inverting iteration, with survey line L based on smooth constraint
1, L
2The electric potential difference data carry out inverting, the inverting objective function is seen formula 3;
Φ=(Δd-AΔm)
T(Δd-AΔm)+λ(CΔm)
T(CΔm) (3)
In the formula, Φ is the inverting objective function, and Δ d is the difference of observation data, and A is a partial derivative matrix, and Δ m is the model parameter correction of grid, and C is the smoothness matrix, and λ is Lagrangian constant, has determined the weight of smooth constraint.
Formula (3) is carried out can obtaining equation (2), when carrying out inversion imaging after the variation; Setting model parameter initial value at first, this initial value resistivity all is set to 1000 Ω .m, is just drilling calculating; The observation data of the theory of computation; Compare judgement with data measured, if variance is less than setting value (this instance relative value 3%), then output model parameter; Otherwise, calculate partial derivative matrix A, the susceptibility Matrix C is utilized equation (2) can obtain the model parameter increment in each inverting iteration, thereby is tried to achieve the model parameter m in the next iteration
(k+1), until the variance of theoretical observation data of just drilling and measurement data less than setting value, the output model parameter.Thereby inverting obtains the 3 D resistivity image of workplace front geological body, realizes the three-dimensional localization of water content.The data interpretation three-dimensional imaging result of this detection, as shown in Figure 4, there are many water guide cracks in YK123+672~YK123+652, this section country rock, YK123+652~YK123+642, this section country rock watery is stronger.
(3) carry out the estimation of the water-bearing structure water yield; Data were carried out the water yield estimation of workplace the place ahead when employing partly declined; With record partly decline the time T data, two and half of identical point position data when declining are done poorly when 0.1A, the power supply of 0.2A two primary currents, and draw the coordinate diagram that concerns of difference data and horizontal range when partly declining respectively; As shown in Figure 5, the difference when calculating partly declines and the envelope area S of abscissa axis
t, and envelope area S
tWith the linear positive correlation of big or small V of the water yield, be the response of the water yield on the occasion of part, by the size of estimating the water yield on the occasion of the area of part; A nearly again step is that the vertical direction of X axle, workplace is the Y axle with the horizontal direction of workplace, draws out envelope area S
tTwo-dimensional cross section, realize the estimation of workplace the place ahead water yield, as shown in Figure 6, the water yield and Fig. 4 are comparatively identical.
Claims (5)
1. underground works chromatographic focusing induced polarization forward probe method is characterized in that, may further comprise the steps: at first adopt the working method that focuses on, on underground works workplace outline line, arrange a circle guarded electrode A
1,, establish electrodes of A on the survey line arranging many parallel surveys line to the differing heights at top apart from base plate
0With receiving electrode M
1, M
2, also be provided with transmitting electrode B and receiving electrode N on the base plate; The method that adopts chromatography then is image data successively, with electrodes of A
0With guarded electrode A
1Infeed electric current in the same way, test survey line receiving electrode M respectively
1With the electric potential difference of base plate receiving electrode N with when partly declining, and change the survey line receiving electrode M that size of current test different electric flows down
1With the electric potential difference of base plate receiving electrode N with when partly declining; Survey line electric potential difference data with recording are carried out inverting, obtain the 3 D resistivity image of workplace front geological body, realize the three-dimensional localization of water content; With record partly decline the time draw difference data and horizontal range when partly declining concern coordinate diagram, the difference when calculating partly declines and the envelope area S of abscissa axis
t, draw out envelope area S
tTwo-dimensional cross section, realize the estimation of workplace the place ahead water yield.
2. underground works chromatographic focusing induced polarization forward probe method according to claim 1 is characterized in that, may further comprise the steps:
(1) adopts the working method that focuses on, at first on underground works workplace outline line, arrange the equally spaced guarded electrode A of a circle
1, apart from base plate 1/3H, two parallel survey line L of 2/3H height placement
1, L
2, every survey line arranges that electrode number is respectively N
1, N
2, and N
1=N
2, adjacent two distance between electrodes are d, wherein H is the height of vault apart from base plate; Survey line L
1, L
2On all establish electrodes of A
0With receiving electrode M
1, M
2, receiving electrode M
1, M
2Adjustable positions and with survey line on the electrode position arranged corresponding, also be provided with transmitting electrode B and receiving electrode N on the base plate;
(2) method that adopts chromatography is image data successively
1. survey line L
1Electrodes of A
0Infeed positive current I
0, transmitting electrode B infeeds negative current-I
0, guarded electrode A
1Infeed and electrodes of A
0Positive current in the same way keeps electrodes of A
0With potential electrode M
1, M
2Spacing A
0M
1=M
1M
2=nd, wherein d is adjacent two distance between electrodes, n is natural number and 1≤n≤7, moves M
1, M
2The other end from survey line one end to survey line; M
1, M
2Whenever satisfy A
0M
1=M
1M
2=nd moves once, electrode A
0, M
1And M
2Every measurement 1 time, and the adjusting supply current makes M
1M
2Between electric potential difference U
M1M2=0, and note M
1Electric potential difference U1 between the N
M1NT when partly declining;
2. survey line L
2Repeat sequence L
1Measuring process 1.;
3. electrode A
0Infeed positive current 2I
0, electrode B infeeds negative current-I
0, repeating step 1., 2.;
(3) the 3-d inversion imaging of carrying out water-bearing structure is located, with the survey line L that records in the step (2)
1, L
2Electric potential difference data U1
M1NCarry out inverting, obtain the 3 D resistivity image of workplace front geological body, realize the three-dimensional localization of water content;
Data were carried out the water yield estimation of workplace the place ahead when (4) employing partly declined; With record in the step (2) partly decline the time T data; When decline in two and half of size two primary currents when power supply identical point position data do poor; And draw the coordinate diagram that concerns of difference data and horizontal range when partly declining respectively, calculate difference and the envelope area S of abscissa axis when partly declining
t, and envelope area S
tWith the linear positive correlation of big or small V of the water yield, be the response of the water yield on the occasion of part, by the size on the occasion of the area estimation water yield of part, nearly again one to go on foot horizontal direction with workplace be that the vertical direction of X axle, workplace is the Y axle, draws out envelope area S
tTwo-dimensional cross section, realize the estimation of workplace the place ahead water yield.
3. underground works chromatographic focusing induced polarization forward probe method according to claim 2 is characterized in that, guarded electrode A in the step (1)
1Number be 4~8, N
1=N
2=15.
4. underground works chromatographic focusing induced polarization forward probe method according to claim 2 is characterized in that, I in the step (2)
0Span 0.1A~0.5A.
5. underground works chromatographic focusing induced polarization forward probe method according to claim 2 is characterized in that, step (3) is specially with having the 3 D resistivity inversion equation of smooth constraint:
(A
TA+λC
TC)Δm=A
TΔd (1)
m
(k+1)=m
(k)+Δm (2)
In the formula, Δ d is the difference of observation data, and A is a partial derivative matrix, and Δ m is the model parameter correction of grid, and C is the smoothness matrix, and λ is Lagrangian constant, has determined the weight of smooth constraint, m
(k)Be the k time model parameter vector in the iteration, m
(k+1)Be the model parameter vector in the inverting iteration next time;
At first setting model parameter initial value is just being drilled calculating, and the observation data of the theory of computation compares judgement with data measured, if variance less than setting value, output model parameter then; Otherwise, calculate partial derivative matrix A, the smoothness Matrix C is utilized equation (2) to obtain the model parameter increment in each inverting iteration, thereby is tried to achieve the model parameter m in the next iteration
(k+1), less than setting value, the output model parameter obtains the 3 D resistivity image of workplace front geological body with the model parameter inverting that obtains until the variance of theoretical observation data of just drilling and measurement data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103766381A CN102508303B (en) | 2011-11-23 | 2011-11-23 | Advanced detection method for focusing chromatography induced polarization of underground engineering |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103766381A CN102508303B (en) | 2011-11-23 | 2011-11-23 | Advanced detection method for focusing chromatography induced polarization of underground engineering |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102508303A true CN102508303A (en) | 2012-06-20 |
CN102508303B CN102508303B (en) | 2013-06-05 |
Family
ID=46220408
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011103766381A Active CN102508303B (en) | 2011-11-23 | 2011-11-23 | Advanced detection method for focusing chromatography induced polarization of underground engineering |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102508303B (en) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064121A (en) * | 2013-01-15 | 2013-04-24 | 山东大学 | Underground construction water burst source advance geology forecast three-dimensional focus induced polarization equipment |
CN103592697A (en) * | 2013-10-31 | 2014-02-19 | 山东大学 | Combined advanced prediction method based on priori form information constraint for tunnel and other underground constructions |
CN103645514A (en) * | 2013-12-25 | 2014-03-19 | 山东大学 | Underground engineering advanced detection method and system for resistivity of multi-same-source electrode array |
CN103645296A (en) * | 2013-12-03 | 2014-03-19 | 中铁十三局集团有限公司 | Underground water detection method |
WO2014106308A1 (en) * | 2013-01-07 | 2014-07-10 | 山东大学 | Advanced detection device, system and method using forward three-dimensional induced polarization for tbm construction tunnel |
WO2014110689A1 (en) * | 2013-01-15 | 2014-07-24 | 山东大学 | Three-dimensional focus induced polarization device for advanced geological forecast of water burst source in underground construction |
CN104088630A (en) * | 2014-07-09 | 2014-10-08 | 中国石油集团长城钻探工程有限公司 | Method for obtaining stratum boundary distance through azimuth resistivity logging instrument while drilling |
US9091779B2 (en) | 2013-01-07 | 2015-07-28 | Shandong University | Advanced detector system and method using forward three-dimensional induced polarization method for TBM construction tunnel |
CN105044790A (en) * | 2015-06-23 | 2015-11-11 | 华中科技大学 | System and method for determining orientation of anomalous body on the basis of circumferential electrode scanning |
CN106405672A (en) * | 2016-10-27 | 2017-02-15 | 山东大学 | Interference elimination method for three-dimensional induced-polarization advanced prediction in TBM construction tunnel |
CN106772621A (en) * | 2017-01-24 | 2017-05-31 | 山东大学 | A kind of nearly comprehensive resistivity geological extra-forecast method |
CN107015282A (en) * | 2017-06-05 | 2017-08-04 | 山东理工大学 | Underground space country rock disaster resistivity method of real-time based on anchor pole |
CN109187850A (en) * | 2018-07-03 | 2019-01-11 | 浙江农林大学 | A kind of laboratory intelligent and safe management system and its working method |
CN110133726A (en) * | 2019-04-25 | 2019-08-16 | 中铁二院工程集团有限责任公司 | Railway tunnel airborne electromagnetic prospecting arrangement of measuring-line method |
CN110889249A (en) * | 2019-11-07 | 2020-03-17 | 中国建筑第四工程局有限公司 | Resistivity karst cave identification method based on population evolution algorithm |
CN111239819A (en) * | 2020-02-12 | 2020-06-05 | 西安理工大学 | Direct envelope inversion method with polarity based on seismic channel attribute analysis |
CN113391359A (en) * | 2021-07-26 | 2021-09-14 | 中铁工程装备集团有限公司 | Real-time detection electric advanced geological prediction system and method |
CN115508897A (en) * | 2022-09-16 | 2022-12-23 | 山东大学 | Method and system for rapidly calculating half decay time by induced polarization |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3967190A (en) * | 1974-12-23 | 1976-06-29 | Zonge Kenneth L | Method using induced polarization for ore discrimination in disseminated earth deposits |
CN1904644A (en) * | 2006-08-08 | 2007-01-31 | 武汉长盛工程检测技术开发有限公司 | Network tunnel real time continuous leading preinforming method and device |
US7388391B2 (en) * | 1999-11-24 | 2008-06-17 | Micron Technology, Inc. | Method for evaluating at least one electrical conducting structure of an electronic component |
CN300904644S (en) * | 2007-12-13 | 2009-04-08 | 房梅宏 | English writing grid paper |
CN101672929A (en) * | 2009-09-23 | 2010-03-17 | 山东大学 | Method for exploring water bearing stratum in front of tunnel heading in tunnel and estimating water yield of water bearing stratum |
CN101706585A (en) * | 2009-10-28 | 2010-05-12 | 中南大学 | Method for electrically forecasting danger in underground tunneling engineering |
-
2011
- 2011-11-23 CN CN2011103766381A patent/CN102508303B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3967190A (en) * | 1974-12-23 | 1976-06-29 | Zonge Kenneth L | Method using induced polarization for ore discrimination in disseminated earth deposits |
US7388391B2 (en) * | 1999-11-24 | 2008-06-17 | Micron Technology, Inc. | Method for evaluating at least one electrical conducting structure of an electronic component |
CN1904644A (en) * | 2006-08-08 | 2007-01-31 | 武汉长盛工程检测技术开发有限公司 | Network tunnel real time continuous leading preinforming method and device |
CN300904644S (en) * | 2007-12-13 | 2009-04-08 | 房梅宏 | English writing grid paper |
CN101672929A (en) * | 2009-09-23 | 2010-03-17 | 山东大学 | Method for exploring water bearing stratum in front of tunnel heading in tunnel and estimating water yield of water bearing stratum |
CN101706585A (en) * | 2009-10-28 | 2010-05-12 | 中南大学 | Method for electrically forecasting danger in underground tunneling engineering |
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014106308A1 (en) * | 2013-01-07 | 2014-07-10 | 山东大学 | Advanced detection device, system and method using forward three-dimensional induced polarization for tbm construction tunnel |
US9091779B2 (en) | 2013-01-07 | 2015-07-28 | Shandong University | Advanced detector system and method using forward three-dimensional induced polarization method for TBM construction tunnel |
WO2014110689A1 (en) * | 2013-01-15 | 2014-07-24 | 山东大学 | Three-dimensional focus induced polarization device for advanced geological forecast of water burst source in underground construction |
CN103064121B (en) * | 2013-01-15 | 2014-01-29 | 山东大学 | Underground construction water burst source advance geology forecast three-dimensional focus induced polarization equipment |
US9256003B2 (en) | 2013-01-15 | 2016-02-09 | Shan Dong University | Three-dimensional focusing induced polarization equipment for advanced geological prediction of water inrush disaster source in underground engineering |
CN103064121A (en) * | 2013-01-15 | 2013-04-24 | 山东大学 | Underground construction water burst source advance geology forecast three-dimensional focus induced polarization equipment |
CN103592697A (en) * | 2013-10-31 | 2014-02-19 | 山东大学 | Combined advanced prediction method based on priori form information constraint for tunnel and other underground constructions |
CN103592697B (en) * | 2013-10-31 | 2016-08-17 | 山东大学 | The associating advanced prediction method of the underground engineering priori shape information constraints such as tunnel |
CN103645296A (en) * | 2013-12-03 | 2014-03-19 | 中铁十三局集团有限公司 | Underground water detection method |
CN103645514A (en) * | 2013-12-25 | 2014-03-19 | 山东大学 | Underground engineering advanced detection method and system for resistivity of multi-same-source electrode array |
CN103645514B (en) * | 2013-12-25 | 2015-01-07 | 山东大学 | Underground engineering advanced detection method and system for resistivity of multi-same-source electrode array |
CN104088630A (en) * | 2014-07-09 | 2014-10-08 | 中国石油集团长城钻探工程有限公司 | Method for obtaining stratum boundary distance through azimuth resistivity logging instrument while drilling |
CN104088630B (en) * | 2014-07-09 | 2016-08-10 | 中国石油集团长城钻探工程有限公司 | With boring the method that azimuthal resistivity logger obtains stratigraphic boundary distance |
CN105044790B (en) * | 2015-06-23 | 2017-07-25 | 华中科技大学 | A kind of anomalous body orientation scanned based on circumferential electrode judges system and method |
CN105044790A (en) * | 2015-06-23 | 2015-11-11 | 华中科技大学 | System and method for determining orientation of anomalous body on the basis of circumferential electrode scanning |
CN106405672A (en) * | 2016-10-27 | 2017-02-15 | 山东大学 | Interference elimination method for three-dimensional induced-polarization advanced prediction in TBM construction tunnel |
CN106772621A (en) * | 2017-01-24 | 2017-05-31 | 山东大学 | A kind of nearly comprehensive resistivity geological extra-forecast method |
CN107015282A (en) * | 2017-06-05 | 2017-08-04 | 山东理工大学 | Underground space country rock disaster resistivity method of real-time based on anchor pole |
CN107015282B (en) * | 2017-06-05 | 2018-12-07 | 山东理工大学 | Underground space country rock disaster resistivity method of real-time based on anchor pole |
CN109187850A (en) * | 2018-07-03 | 2019-01-11 | 浙江农林大学 | A kind of laboratory intelligent and safe management system and its working method |
CN109187850B (en) * | 2018-07-03 | 2021-03-16 | 浙江农林大学 | Intelligent laboratory safety management system and working method thereof |
CN110133726A (en) * | 2019-04-25 | 2019-08-16 | 中铁二院工程集团有限责任公司 | Railway tunnel airborne electromagnetic prospecting arrangement of measuring-line method |
CN110133726B (en) * | 2019-04-25 | 2024-03-22 | 中铁二院工程集团有限责任公司 | Method for arranging exploration survey lines of railway tunnel aviation electromagnetic method |
CN110889249B (en) * | 2019-11-07 | 2023-07-25 | 中国建筑第四工程局有限公司 | Method for identifying resistivity karst cave based on population evolution algorithm |
CN110889249A (en) * | 2019-11-07 | 2020-03-17 | 中国建筑第四工程局有限公司 | Resistivity karst cave identification method based on population evolution algorithm |
CN111239819A (en) * | 2020-02-12 | 2020-06-05 | 西安理工大学 | Direct envelope inversion method with polarity based on seismic channel attribute analysis |
CN111239819B (en) * | 2020-02-12 | 2022-05-03 | 西安理工大学 | Direct envelope inversion method with polarity based on seismic channel attribute analysis |
CN113391359B (en) * | 2021-07-26 | 2022-03-11 | 中铁工程装备集团有限公司 | Real-time detection electric advanced geological prediction system and method |
CN113391359A (en) * | 2021-07-26 | 2021-09-14 | 中铁工程装备集团有限公司 | Real-time detection electric advanced geological prediction system and method |
CN115508897A (en) * | 2022-09-16 | 2022-12-23 | 山东大学 | Method and system for rapidly calculating half decay time by induced polarization |
Also Published As
Publication number | Publication date |
---|---|
CN102508303B (en) | 2013-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102508303B (en) | Advanced detection method for focusing chromatography induced polarization of underground engineering | |
CN103645514B (en) | Underground engineering advanced detection method and system for resistivity of multi-same-source electrode array | |
CN102508308B (en) | Transient electromagnetic multipoint array-type detection method of tunnel | |
CN102419455B (en) | Interwell parallel resistivity CT (computed tomography) testing method | |
CN102156301B (en) | Advanced-prediction observation system while drilling | |
Ungureanu et al. | Use of electric resistivity tomography (ERT) for detecting underground voids on highly anthropized urban construction sites | |
CN107742015A (en) | DC lasering electric method Three-dimensional Numerical Simulation Method based on any dipole dipole device | |
CN102866417A (en) | Device and method for seismic cross hole computed tomography (CT) detection and tomography of underground cave | |
CN109031428A (en) | A kind of hilly pond embankment seepage detection method | |
CN106772678A (en) | A kind of wellhole many reference amounts method for surveying of rock deformation destructive characteristics | |
CN102768369A (en) | Roadway drivage drilling induced polarization advance water probing forecasting method, device and probe | |
CN202837558U (en) | Underground karst cave earthquake cross-hole CT (computer tomography) detection and tomographic imaging device | |
CN110187398A (en) | The multiple electrode system detection method of objective body between a kind of searching well | |
CN103967476A (en) | Drilling geophysical prospecting advanced detecting device and method | |
CN110221340A (en) | A kind of set of tunneling construction method for forecasting advanced geology | |
CN103454680A (en) | Method for calculating vertical coverage times of Walk-away VSP observing system | |
CN113156518B (en) | Real-time advanced detection method for vector resistivity of water-containing disaster body | |
CN103955000A (en) | Three-dimensional electrical sounding method | |
CN104614770A (en) | Exploration device and method of resource gobs | |
CN102182437B (en) | Method for determining and eliminating hydraulic fracture stress boundary of coal mine underground drilling | |
CN112255693A (en) | Goaf filling space detection method under coal mine fully mechanized caving coal mining process condition | |
CN101793973A (en) | While-drilling electric method | |
CN104502951A (en) | Rayleigh wave railway substructure monitoring cavity three-dimensional positioning method | |
CN106772637B (en) | A kind of forward probe method of the underground coal mine based on slot wave | |
CN112485829A (en) | Method for calculating dynamic water content by resistivity method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |