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
 workplace
 survey line
 data
 electrodes
 Prior art date
Links
 238000004587 chromatography analysis Methods 0.000 title claims abstract description 12
 238000001514 detection method Methods 0.000 title abstract description 11
 XLYOFNOQVPJJNPUHFFFAOYSAN water Substances data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwJyBoZWlnaHQ9JzMwMCcgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMTAwLjUwMScgeT0nMTcwJyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+SDwvdGV4dD4KPHRleHQgeD0nMTI2LjExNCcgeT0nMTg2JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjI2cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+MjwvdGV4dD4KPHRleHQgeD0nMTM4JyB5PScxNzAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8L3N2Zz4K data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NScgaGVpZ2h0PSc4NScgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMTMuMzEyOCcgeT0nNTMuNTkwOScgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZToyM3B4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPkg8L3RleHQ+Cjx0ZXh0IHg9JzI4LjE1NjknIHk9JzYyLjg2MzYnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MTVweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID4yPC90ZXh0Pgo8dGV4dCB4PSczNS4wNDU1JyB5PSc1My41OTA5JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjIzcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPC9zdmc+Cg== O XLYOFNOQVPJJNPUHFFFAOYSAN 0.000 claims abstract description 64
 238000005259 measurement Methods 0.000 claims abstract description 15
 239000000523 sample Substances 0.000 claims description 17
 239000011159 matrix material Substances 0.000 claims description 12
 230000004807 localization Effects 0.000 claims description 10
 238000003384 imaging method Methods 0.000 claims description 9
 238000000034 method Methods 0.000 claims description 9
 238000010586 diagram Methods 0.000 claims description 7
 238000005553 drilling Methods 0.000 claims description 4
 230000000875 corresponding Effects 0.000 claims description 2
 238000004070 electrodeposition Methods 0.000 claims description 2
 238000010276 construction Methods 0.000 description 6
 238000010291 electrical method Methods 0.000 description 5
 238000005516 engineering process Methods 0.000 description 4
 230000000694 effects Effects 0.000 description 3
 239000011435 rock Substances 0.000 description 3
 230000001105 regulatory Effects 0.000 description 2
 241000893018 Armeria Species 0.000 description 1
 210000002706 Plastids Anatomy 0.000 description 1
 238000009412 basement excavation Methods 0.000 description 1
 230000005540 biological transmission Effects 0.000 description 1
 238000004364 calculation method Methods 0.000 description 1
 230000005684 electric field Effects 0.000 description 1
 238000009114 investigational therapy Methods 0.000 description 1
 239000003129 oil well Substances 0.000 description 1
 201000007094 prostatitis Diseases 0.000 description 1
 230000001360 synchronised Effects 0.000 description 1
Classifications

 Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSSSECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE 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
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 threedimensional 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 twodimensional 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 threedimensional location of an aquifer and solves problem of incapability of location in a traditional focusing electric method during singlepoint 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 waterbearing structure, promptly need realize the threedimensional imaging of waterbearing structure; 2. need predict the water yield in the waterbearing structure, this is to formulate prominent water burst prophylacticotherapeutic 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 nonfocusing type and singlepoint focus type.Wherein, Near the side that three electrode method nonfocusing type forward probe method can't shield the survey line disturbs (like lowresistance 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 nonfocusing 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 BoreTunneling 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 realtime 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
_{1}With central electrode A
_{0}, and link to each other with onthespot main frame with cable; (2) onthespot main frame is through guarded electrode A
_{1}With central electrode A
_{0}The 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) onthespot 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 singlepoint focus type; Single potential electrode just only is set on test surface, does not adopt chromatography to survey, can't obtain the threedimensional 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 abovementioned 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 threedimensional 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
_{0}With 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
_{0}With guarded electrode A
_{1}Infeed electric current in the same way, test survey line receiving electrode M respectively
_{1}With 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
_{1}With 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 threedimensional 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
_{t}Twodimensional 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 48) 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
_{0}With receiving electrode M
_{1}, M
_{2}Place survey line L
_{1}, L
_{2}(receiving electrode M
_{1}, M
_{2}Adjustable positions and with survey line on the electrode position arranged corresponding, be in order to guarantee M
_{1}M
_{2}Between 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 100150m, receiving electrode N places on the base plate;
(2) method that adopts chromatography is image data successively,
1. survey line L
_{1}Electrodes of A
_{0}Infeed positive current I
_{0}, transmitting electrode B infeeds negative currentI
_{0}, guarded electrode A
_{1}Infeed and electrodes of A
_{0}Positive current in the same way keeps electrodes of A
_{0}With potential electrode M
_{1}, M
_{2}Spacing A
_{0}M
_{1}=M
_{1}M
_{2}=nd, wherein d is adjacent two distance between electrodes, n is natural number and 1≤n≤7, moves M
_{1}, M
_{2}The other end from survey line one end to survey line; M
_{1}, M
_{2}Whenever satisfy A
_{0}M
_{1}=M
_{1}M
_{2}=nd moves once, electrode A
_{0}, M
_{1}And M
_{2}Every measurement 1 time is regulated supply current and is made M
_{1}M
_{2}Between electric potential difference U
_{M1M2}=0, and note M
_{1}Electric potential difference U1 between the N
_{M1N}T when partly declining;
2. survey line L
_{2}Repeat sequence L
_{1}Measuring process 1.;
3. electrode A
_{0}Infeed positive current 2I
_{0}, electrode B infeeds negative currentI
_{0}, repeating step 1., 2.;
(3) the 3d inversion imaging of carrying out waterbearing 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
_{2}Electric potential difference data U1
_{M1N}Carry out inverting, obtain the 3 D resistivity image of workplace front geological body, realize the threedimensional localization of water content; Concrete steps are following:
3 D resistivity inversion equation with smooth constraint.
(A
^{T}A+λC
^{T}C)Δ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 threedimensional 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
_{t}With 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
_{t}Twodimensional 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
_{2}Last electrodes of A
_{0}Electric 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 threedimensionalreconstruction of workplace front geological body electrical structure, realized threedimensional localization, solved the spot measurement of traditional focus class methods and a difficult problem that can't threedimensional 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 waterbearing 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 threedimensionalreconstruction of workplace front geological body electrical structure based on smooth constraint; Especially realize the water yield estimation of waterbearing structure, had following characteristic:
1. adopt the focusingdetection 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 threedimensional 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 crosssection of coordinate axis envelope area when partly declining, has realized the water yield estimation of waterbearing 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 3d 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 twodimensional 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 3d inversion location and water yield estimation; Has good forward probe directive property; Realize the poorly threedimensional localization of plastid and the estimation of the water yield such as waterbearing 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 crosssection is an arch, and the high 10m of span 12m water burst occurs near the mileage YK123+672, and possibly there is waterbearing 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
_{0}With receiving electrode M
_{1}, M
_{2}Place 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
_{1}Electrodes of A
_{0}Infeed positive current 0.1A, transmitting electrode B is a negative pole, guarded electrode A
_{1}Infeed positive current 0.1A, keep electrodes of A
_{0}With potential electrode M
_{1}, M
_{2}Spacing A
_{0}M
_{1}=M
_{1}M
_{2}=nd, M
_{1}, M
_{2}Successively from survey line L
_{1}One 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
_{1}And M
_{2}Every measurement 1 time is regulated supply current and is made M
_{1}M
_{2}Between electric potential difference U
_{M1M2}=0, and note M
_{1}Electric potential difference U1 between the N
_{M1N}T when partly declining; Measure survey line L according to above method
_{2}, adjust electrode A then
_{0}Infeed 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
_{2}The electric potential difference data carry out inverting, the inverting objective function is seen formula 3;
Φ＝(ΔdAΔm)
^{T}(ΔdAΔ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 threedimensional localization of water content.The data interpretation threedimensional 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 waterbearing 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
_{t}With 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
_{t}Twodimensional 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
_{0}With 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
_{0}With guarded electrode A
_{1}Infeed electric current in the same way, test survey line receiving electrode M respectively
_{1}With 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
_{1}With 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 threedimensional 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
_{t}Twodimensional 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
_{2}On all establish electrodes of A
_{0}With receiving electrode M
_{1}, M
_{2}, receiving electrode M
_{1}, M
_{2}Adjustable 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
_{1}Electrodes of A
_{0}Infeed positive current I
_{0}, transmitting electrode B infeeds negative currentI
_{0}, guarded electrode A
_{1}Infeed and electrodes of A
_{0}Positive current in the same way keeps electrodes of A
_{0}With potential electrode M
_{1}, M
_{2}Spacing A
_{0}M
_{1}=M
_{1}M
_{2}=nd, wherein d is adjacent two distance between electrodes, n is natural number and 1≤n≤7, moves M
_{1}, M
_{2}The other end from survey line one end to survey line; M
_{1}, M
_{2}Whenever satisfy A
_{0}M
_{1}=M
_{1}M
_{2}=nd moves once, electrode A
_{0}, M
_{1}And M
_{2}Every measurement 1 time, and the adjusting supply current makes M
_{1}M
_{2}Between electric potential difference U
_{M1M2}=0, and note M
_{1}Electric potential difference U1 between the N
_{M1N}T when partly declining;
2. survey line L
_{2}Repeat sequence L
_{1}Measuring process 1.;
3. electrode A
_{0}Infeed positive current 2I
_{0}, electrode B infeeds negative currentI
_{0}, repeating step 1., 2.;
(3) the 3d inversion imaging of carrying out waterbearing structure is located, with the survey line L that records in the step (2)
_{1}, L
_{2}Electric potential difference data U1
_{M1N}Carry out inverting, obtain the 3 D resistivity image of workplace front geological body, realize the threedimensional 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
_{t}With 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
_{t}Twodimensional 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)
_{1}Number 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)
_{0}Span 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
^{T}A+λC
^{T}C)Δ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)  20111123  20111123  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)  20111123  20111123  Advanced detection method for focusing chromatography induced polarization of underground engineering 
Publications (2)
Publication Number  Publication Date 

CN102508303A true CN102508303A (en)  20120620 
CN102508303B CN102508303B (en)  20130605 
Family
ID=46220408
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN2011103766381A CN102508303B (en)  20111123  20111123  Advanced detection method for focusing chromatography induced polarization of underground engineering 
Country Status (1)
Country  Link 

CN (1)  CN102508303B (en) 
Cited By (13)
Publication number  Priority date  Publication date  Assignee  Title 

CN103064121A (en) *  20130115  20130424  山东大学  Underground construction water burst source advance geology forecast threedimensional focus induced polarization equipment 
CN103592697A (en) *  20131031  20140219  山东大学  Combined advanced prediction method based on priori form information constraint for tunnel and other underground constructions 
CN103645296A (en) *  20131203  20140319  中铁十三局集团有限公司  Underground water detection method 
CN103645514A (en) *  20131225  20140319  山东大学  Underground construction advanced detection method and system based on resistivity of multiple same source electrode arrays 
WO2014106308A1 (en) *  20130107  20140710  山东大学  Advanced detection device, system and method using forward threedimensional induced polarization for tbm construction tunnel 
WO2014110689A1 (en) *  20130115  20140724  山东大学  Threedimensional focus induced polarization device for advanced geological forecast of water burst source in underground construction 
CN104088630A (en) *  20140709  20141008  中国石油集团长城钻探工程有限公司  Method for obtaining stratum boundary distance through azimuth resistivity logging instrument while drilling 
US9091779B2 (en)  20130107  20150728  Shandong University  Advanced detector system and method using forward threedimensional induced polarization method for TBM construction tunnel 
CN105044790A (en) *  20150623  20151111  华中科技大学  System and method for determining orientation of anomalous body on the basis of circumferential electrode scanning 
CN106405672A (en) *  20161027  20170215  山东大学  Interference elimination method for threedimensional inducedpolarization advanced prediction in TBM construction tunnel 
CN106772621A (en) *  20170124  20170531  山东大学  A kind of nearly comprehensive resistivity geological extraforecast method 
CN107015282A (en) *  20170605  20170804  山东理工大学  Underground space country rock disaster resistivity method of realtime based on anchor pole 
CN109187850A (en) *  20180703  20190111  浙江农林大学  A kind of laboratory intelligent and safe management system and its working method 
Citations (6)
Publication number  Priority date  Publication date  Assignee  Title 

US3967190A (en) *  19741223  19760629  Zonge Kenneth L  Method using induced polarization for ore discrimination in disseminated earth deposits 
CN1904644A (en) *  20060808  20070131  武汉长盛工程检测技术开发有限公司  Network tunnel real time continuous leading preinforming method and device 
US7388391B2 (en) *  19991124  20080617  Micron Technology, Inc.  Method for evaluating at least one electrical conducting structure of an electronic component 
CN300904644S (en) *  20071213  20090408  房梅宏  English writing grid paper 
CN101672929A (en) *  20090923  20100317  山东大学  Method for exploring water bearing stratum in front of tunnel heading in tunnel and estimating water yield of water bearing stratum 
CN101706585A (en) *  20091028  20100512  中南大学  Method for electrically forecasting danger in underground tunneling engineering 

2011
 20111123 CN CN2011103766381A patent/CN102508303B/en active IP Right Grant
Patent Citations (6)
Publication number  Priority date  Publication date  Assignee  Title 

US3967190A (en) *  19741223  19760629  Zonge Kenneth L  Method using induced polarization for ore discrimination in disseminated earth deposits 
US7388391B2 (en) *  19991124  20080617  Micron Technology, Inc.  Method for evaluating at least one electrical conducting structure of an electronic component 
CN1904644A (en) *  20060808  20070131  武汉长盛工程检测技术开发有限公司  Network tunnel real time continuous leading preinforming method and device 
CN300904644S (en) *  20071213  20090408  房梅宏  English writing grid paper 
CN101672929A (en) *  20090923  20100317  山东大学  Method for exploring water bearing stratum in front of tunnel heading in tunnel and estimating water yield of water bearing stratum 
CN101706585A (en) *  20091028  20100512  中南大学  Method for electrically forecasting danger in underground tunneling engineering 
Cited By (21)
Publication number  Priority date  Publication date  Assignee  Title 

US9091779B2 (en)  20130107  20150728  Shandong University  Advanced detector system and method using forward threedimensional induced polarization method for TBM construction tunnel 
WO2014106308A1 (en) *  20130107  20140710  山东大学  Advanced detection device, system and method using forward threedimensional induced polarization for tbm construction tunnel 
US9256003B2 (en)  20130115  20160209  Shan Dong University  Threedimensional focusing induced polarization equipment for advanced geological prediction of water inrush disaster source in underground engineering 
CN103064121B (en) *  20130115  20140129  山东大学  Underground construction water burst source advance geology forecast threedimensional focus induced polarization equipment 
CN103064121A (en) *  20130115  20130424  山东大学  Underground construction water burst source advance geology forecast threedimensional focus induced polarization equipment 
WO2014110689A1 (en) *  20130115  20140724  山东大学  Threedimensional focus induced polarization device for advanced geological forecast of water burst source in underground construction 
CN103592697A (en) *  20131031  20140219  山东大学  Combined advanced prediction method based on priori form information constraint for tunnel and other underground constructions 
CN103592697B (en) *  20131031  20160817  山东大学  The associating advanced prediction method of the underground engineering priori shape information constraints such as tunnel 
CN103645296A (en) *  20131203  20140319  中铁十三局集团有限公司  Underground water detection method 
CN103645514B (en) *  20131225  20150107  山东大学  Underground construction advanced detection method and system based on resistivity of multiple same source electrode arrays 
CN103645514A (en) *  20131225  20140319  山东大学  Underground construction advanced detection method and system based on resistivity of multiple same source electrode arrays 
CN104088630A (en) *  20140709  20141008  中国石油集团长城钻探工程有限公司  Method for obtaining stratum boundary distance through azimuth resistivity logging instrument while drilling 
CN104088630B (en) *  20140709  20160810  中国石油集团长城钻探工程有限公司  With boring the method that azimuthal resistivity logger obtains stratigraphic boundary distance 
CN105044790B (en) *  20150623  20170725  华中科技大学  A kind of anomalous body orientation scanned based on circumferential electrode judges system and method 
CN105044790A (en) *  20150623  20151111  华中科技大学  System and method for determining orientation of anomalous body on the basis of circumferential electrode scanning 
CN106405672A (en) *  20161027  20170215  山东大学  Interference elimination method for threedimensional inducedpolarization advanced prediction in TBM construction tunnel 
CN106772621A (en) *  20170124  20170531  山东大学  A kind of nearly comprehensive resistivity geological extraforecast method 
CN107015282A (en) *  20170605  20170804  山东理工大学  Underground space country rock disaster resistivity method of realtime based on anchor pole 
CN107015282B (en) *  20170605  20181207  山东理工大学  Underground space country rock disaster resistivity method of realtime based on anchor pole 
CN109187850A (en) *  20180703  20190111  浙江农林大学  A kind of laboratory intelligent and safe management system and its working method 
CN109187850B (en) *  20180703  20210316  浙江农林大学  Intelligent laboratory safety management system and working method thereof 
Also Published As
Publication number  Publication date 

CN102508303B (en)  20130605 
Similar Documents
Publication  Publication Date  Title 

Dahlin et al.  Measuring techniques in induced polarisation imaging  
CN102053270B (en)  Sedimentary formation unitbased seismic facies analysis method  
Li et al.  Detecting and monitoring of water inrush in tunnels and coal mines using direct current resistivity method: a review  
Papadopoulos et al.  Two‐dimensional and three‐dimensional resistivity imaging in archaeological site investigation  
CN101706585B (en)  Method for electrically forecasting danger in underground tunneling engineering  
He et al.  Petroleum electromagnetic prospecting advances and case studies in China  
Tamburini et al.  Retrieving surface deformation by PSInSAR™ technology: A powerful tool in reservoir monitoring  
CN102759745B (en)  Method for forecasting carbonate reservoir based on forward modeling of digital geological outcrop model  
CN102200588B (en)  Method for analyzing waveform similarity body curvature of seismic data  
CN101520517B (en)  Method for accurately evaluating targets containing oil gas in clastic rock basin  
CN103454686B (en)  Small scale sedimentary facies based on strata slicing carries out the method and system of reservoir prediction  
CN101980054B (en)  Method for establishing nearsurface velocity model in highdensity seismic static correction processing  
CN102736124B (en)  Tunnel excavation surrounding rock dynamic refined classification method based on integrated parameters  
CN103675944B (en)  The Comprehensive Geophysics method of a kind of sedimentary basin uranium resources exploration  
CN105221133B (en)  A kind of method and apparatus that content of organic carbon of hydrocarbon source rock is determined based on well logging multiparameter  
CN104407394B (en)  Subway shield section boulder refining detecting method based on resistivity transhole CT (computed tomography)  
CN103995296A (en)  Transient electromagnetic method ground hole detection method and device  
AlFares  Contribution of the geophysical methods in characterizing the water leakage in Afamia B dam, Syria  
US10386531B2 (en)  Geological model analysis incorporating crosswell electromagnetic measurements  
CN102841379B (en)  Method for analyzing prestack time migration and speed based on common scatter point channel set  
CN102879805B (en)  Boreholebased and ground combined seismic wave space exploration method  
CN104360398A (en)  Method for positioning constructed wetland blocked area on basis of twodimensional resistivity imaging technology  
CN104280782B (en)  Onedimensional joint inversion method for timefrequency electromagnetic data and magnetotelluric data  
CN106443765B (en)  Municipal engineering seismic survey integrated imaging method based on multi components observation system  
CN203531883U (en)  Well logging equipment 
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 