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 PDF

Info

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
Application number
CN2011103766381A
Other languages
Chinese (zh)
Other versions
CN102508303B (en
Inventor
李术才
刘斌
聂利超
李树忱
张庆松
林春金
宋杰
刘征宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University
Original Assignee
Shandong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong University filed Critical Shandong University
Priority to CN2011103766381A priority Critical patent/CN102508303B/en
Publication of CN102508303A publication Critical patent/CN102508303A/en
Application granted granted Critical
Publication of CN102508303B publication Critical patent/CN102508303B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment 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

Underground works chromatographic focusing induced polarization forward probe method
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.
CN2011103766381A 2011-11-23 2011-11-23 Advanced detection method for focusing chromatography induced polarization of underground engineering Active CN102508303B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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