CN109633764A - A method of the horizontal seepage channel in runoff area is determined using tracer technique - Google Patents
A method of the horizontal seepage channel in runoff area is determined using tracer technique Download PDFInfo
- Publication number
- CN109633764A CN109633764A CN201811549627.7A CN201811549627A CN109633764A CN 109633764 A CN109633764 A CN 109633764A CN 201811549627 A CN201811549627 A CN 201811549627A CN 109633764 A CN109633764 A CN 109633764A
- Authority
- CN
- China
- Prior art keywords
- water
- tracer
- monitoring points
- quality monitoring
- seepage channel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of nuclear radiation, e.g. of natural or induced radioactivity
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of methods for determining the horizontal seepage channel in runoff area using tracer technique, this method is based on the basis of karst area hydrogeologic condition understanding, has hydrogeology peephole first in chosen area, carry out water level unified test, draw ground water field, it determines the hydraulic connection between each hole, and then arranges tracer test, the horizontal view flow velocity in zoning；Then the long term monitoring that representative drilling carries out underground water temperature, conductivity indices is targetedly chosen, the inflection point of ground water movement in hole is judged, determines the covered depth of horizontal seepage channel, calculate the flow velocity of vertical upper underground water in drilling；The last comprehensive analysis vertical difference with flow velocity in horizontal direction determines seepage channel.The present invention provides scientific basis and decision support for the protection and Karst Water Resource management of karst water system.
Description
Technical field
The present invention relates to hydrogeological technical fields, specifically a kind of to determine karst water system water using tracer technique
The method of flat seepage channel.
Background technique
The inhomogeneity of karstic ground water preservation pacifies water head site exploitation, Karst groundwater pollution prevention and control, engineering construction and mining area
Full production causes great harm.Since Karst Fissures develop unevenness, still lack accurate search in rock stratum at different depth at present
There are the technical methods in seepage action of ground water channel.Therefore, judgement precisely determines that the position of karst area seepage channel becomes the country
One of the difficulties of outer karst water science and engineering area research.
Summary of the invention
In view of the above-mentioned problems, the present invention provides a kind of sides for determining the horizontal seepage channel in runoff area using tracer technique
Method, this method scientific in principle, method are simple, easy to operate, and data are easy to obtain, and have stronger practicability.
The technical scheme adopted by the invention to solve the technical problem is that:
A method of the horizontal seepage channel in runoff area is determined using tracer technique, is included the following steps,
The first step investigates the structure of drillings all in survey region, determines hole depth, contains according to hydrogeologic map
Water layer lithology, sealing section depth, borehole coordinate and elevation；
Second step observes the water level in respectively drilling in the first step one by one, draws the ground water field figure in survey region, really
Hydraulic connection between fixed each drilling；
Third step arranges tracer test, calculates horizontal view flow velocity, and according to the horizontal level for determining seepage channel depending on flow velocity
Extending direction；
3.1 determine throwing source point according to the ground water field figure that second step is drawn, and centered on throwing source point, surrounding arranges water level
Observation point, centered on throwing source point, in the downstream of direction of groundwater flow, fan-shaped arrangement water-quality observation point；
3.2 choose tracer；
The Water-quality Monitoring Points determined in 3.3 pairs of steps 3.1 are sampled, and determine the tracer background of each Water-quality Monitoring Points
Value；
3.4 calculate tracer injected volume；
3.5 launch tracer into throwing source point；
The water levels of 3.6 observed stage monitoring points and Water-quality Monitoring Points observes three days once a day, and according to observing
Water level draws water level contour map, then observes the water level contour map drawn, judges to whether there is near Water-quality Monitoring Points
Cone of depression, cone of depression, then be marked the Water-quality Monitoring Points if it exists.
3.7 are sampled monitoring to Water-quality Monitoring Points according to sequence from the near to the distant, record the concentration and appearance of tracer
Time, and according in monitoring data and step 3.6 determine cone of depression situation sampling rate adjusting；
3.8 draw tracer diffusion velocity isochrone and dust trajectory according to the result that step 3.7 monitors, and then calculate water
Look squarely flow velocity；
3.9 calculate the mean flow rate in survey region according to the calculated result of step 3.8；
4th step calculates vertical flow velocity, determines the vertical covered depth of seepage channel.
4.1 choose the Water-quality Monitoring Points that flow velocity is greater than mean flow rate according to the calculated result of step 3.9, and supervise in the water quality
Every 1m monitoring temperature and conductivity in drilling on measuring point；
4.2 draw the temperature profile and conductance profile figure of each Water-quality Monitoring Points according to the monitoring result of step 4.1；
4.3 calculate the groundwater velocity on vertical for a certain Water-quality Monitoring Points；
4.3.1 the temperature profile and conductance profile figure drawn according to step 4.2 determine the temperature of the Water-quality Monitoring Points
The depth that inflection point and conductivity inflection point occur, determines the covered depth of horizontal seepage channel at the Water-quality Monitoring Points；
4.3.2 the vertical flow velocity of the underground water in each inflection point depth that step 4.3.1 is determined is calculated；
4.4 repeat the operation of step 4.3, calculate the underground water vertical flow of each Water-quality Monitoring Points determined in step 4.1
Speed；
5th step draws seepage channel space point according to third step and the 4th step as a result, on the basis of geologic section
Cloth log sheet.
Further, the water level monitoring point of laying and Water-quality Monitoring Points should all include all water-bearing layers in survey region.
Further, the tracer should meet that nontoxic, chemical property is stable, is not easy to be adsorbed and degrade, solubility
High, background value is low and the requirement of high sensitivity.
Further, in step 3.3, same Water-quality Monitoring Points are sampled 2-3 time, and putting down with the tracer concentration that measures
Background value concentration of the mean value as the Water-quality Monitoring Points.
Further, in step 3.5, if throwing source point, there are multiple water-bearing layers, using top layer water-bearing layer as tracer
The dispensing target zone of agent.
Further, in step 3.7, when being detected to the tracer of Water-quality Monitoring Points, with 5 times of background value concentration
Detection standard as tracer.
Further, in step 3.7, the Adjustment principle of sample frequency is,
Increase sample frequency when Water-quality Monitoring Points detect tracer, and the Water-quality Monitoring Points to next group remotely simultaneously
It is monitored；
When the tracer concentration of Water-quality Monitoring Points detection is lower than peak value, sample frequency is reduced；
When the tracer concentration of Water-quality Monitoring Points detection reaches background value, sampling terminates；
There are when cone of depression near Water-quality Monitoring Points, encryption monitoring will be carried out to the Water-quality Monitoring Points.
Further, the principle that monitoring frequency follows in step 4.1 is monthly to test 3 in 2-7 month, 9-December
Secondary, increasing test frequency in January and August is monthly to test 5 times, encryption test 3 days after rain, once a day.
The beneficial effects of the present invention are:
1, scientific in principle, method are simple, easy to operate, and data are easy to obtain, and have stronger practicability.
2, it can work by ongoing geotechnical boring, hydrogeology or engineering geological drilling and carry out water-level observation, be not necessarily to
Special underground water drift net prospecting is put into, the time has been saved, has improved efficiency, had compared with hi-tech economy.
3, for level regards flow velocity, using tracer technique, and using the tracer peak of sample point detection as foundation, just
In understanding, the achievement of acquirement meets that field is practical to have science.
4, for the flow rate detection on vertical, using the observed temperature of different depth, conductivity is foundations in drilling, meter
Result is calculated convenient for explaining and there is no more solutions, the achievement of acquirement meets field reality, has science.
5, on the basis of horizontal flow velocity and accurate vertical flow velocity, the vertical difference with flow velocity in horizontal direction of comprehensive analysis
Also more accurate come the seepage channel that determines, science can provide section for the protection and Karst Water Resource management of karst water system
Learn foundation and decision support.
Detailed description of the invention
Fig. 1 is the ground water field figure in Jinan Spring Region runoff area；
Fig. 2 is tracer diffusion velocity isochrone and dust trajectory；
Fig. 3 is the temperature profile of Water-quality Monitoring Points J2 well；
Fig. 4 is the conductance profile figure of Water-quality Monitoring Points J2 well；
Fig. 5 is ground temperature distribution curve；
Fig. 6 is that there are ground temperature distribution curves when horizontal seepage channel in stratum；
Fig. 7 is Jinan Spring Region runoff area seepage channel spatial distribution log sheet.
Specific embodiment
A kind of method determining the horizontal seepage channel in runoff area using tracer technique includes the following steps.
The first step, drilling investigation.
According to hydrogeologic map, the structure of drillings all in survey region is investigated, determines hole depth, water-bearing layer rock
Property, sealing section depth, borehole coordinate and elevation.
Drilling described herein is including but not limited to drilling and motor-pumped well, all Kong Jun that can be studied hydrogeology
It is contained within.
Second step draws the ground water field figure in survey region.
The water level in respectively drilling in the first step is observed one by one, draws the ground water field figure in survey region, then basis
The ground water field figure of drafting, just, in hole the information such as formation lithology, sealing position determine that the waterpower between each drilling joins to water level
System.
As a kind of specific embodiment, survey region described in the present embodiment is Jinan Spring Region runoff area, and is drawn
The ground water field figure in the region is as shown in Figure 1.
Third step arranges tracer test, calculates horizontal view flow velocity, and according to the horizontal level for determining seepage channel depending on flow velocity
Extending direction.
Ground water field figure in 3.1 survey regions drawn according to second step, which determines, throws source point, is then to throw source point
Center, surrounding arranges water-level observation point, centered on throwing source point, in the downstream of direction of groundwater flow, and fan-shaped arrangement water-quality observation
Point, and the water level monitoring point and Water-quality Monitoring Points laid should all be comprising all water-bearing layers in survey region.
Herein, water level monitoring point and Water-quality Monitoring Points can be same monitoring point, be also possible to different monitorings
Point.
As a kind of specific embodiment, the throwing source point in survey region described in the present embodiment is determined as Xing Jihe, really
Fixed water level monitoring point and Layout of water quality examination points are as shown in Figure 1.
3.2 choose tracers, and the tracer should meet that nontoxic, chemical property is stable, be not easy to be adsorbed and degrade,
Solubility is high, background value is low and the requirement of high sensitivity (being easy to be detected by instrument).
As a kind of specific embodiment, selected tracer is ammonium molybdate in the present embodiment.
The Water-quality Monitoring Points determined in 3.3 pairs of steps 3.1 are sampled, and determine the tracer background of each Water-quality Monitoring Points
Value.
Occurs artificial pollution in order to prevent, same Water-quality Monitoring Points sample the concentration of 2-3 measurement tracer, and flat with it
Background value concentration of the mean value as the Water-quality Monitoring Points.
3.4 calculate tracer injected volume
The calculation formula of tracer injected volume are as follows:
M=a (TQC)b
In formula: M-is injected volume；
A, b-are empirical coefficient, can be obtained by looking into handbook；
T-is migration time of the release position to receiving point；
Q-receiving point karst water flow；
The peak tracer concentration of the preset receiving point of C-.
Release position is to the migration time of receiving point in formula, and there are two types of acquisition modes:
(1) flow velocity of karst water is obtained according to previous existing test, then migration distance is obtained divided by karst water flow velocity
Migration time.
(2) if not doing similar test in the past, water level contour map is drawn according to the water level observed in second step first,
And hydraulic gradient is calculated according to water level contour map, infiltration coefficient is then asked by single-well pumping test, finally uses Darcy formula
Karst water flow velocity is calculated, then migration distance obtains migration time divided by karst water flow velocity.
Water level contour map is drawn herein according to the water level of observation, hydraulic gradient is calculated according to water level contour map, single hole is taken out
Water test is sought infiltration coefficient and is belonged to those skilled in the art by Darcy formula calculating karst water flow velocity known normal
Know, no longer does excessively repeat herein.
As a kind of specific embodiment, the migration time T of release position to receiving point passes through the first side in the present embodiment
Formula obtains.To survey region in this present embodiment, it is previous it is existing test obtained underground karst water flow velocity have 200m/d,
150m/d, 100m/d and 80m/d these types situation, tracer receiving point Baotu Spring flow Q measured value are 160,000 m3/ d, it is default to connect
The peak tracer concentration of sink is 0.015/gm-3, then above-mentioned karst water is calculated according to the calculation formula of tracer injected volume
Under flow velocity, the ammonium molybdate tracer injected volume that source point is thrown in emerging Ji river is as shown in table 1.
The estimation of 1 ammonium molybdate injected volume of table
Interference is formed to prevent rainy season Atmospheric precipitation leaching from infiltrating, the tracer test of this method is selected in no precipitation leaching
Dry season is tested, since dry season view flow velocity is lower, to guarantee that concentrations select injected volume to show for 280kg as this
The injected volume of track test.
Further, in order to avoid injected volume exceeds national standard, injected volume need to be verified, that is, sets flow velocity as most
Big flow velocity 200m/d, when injected volume is 280kg, the peak concentration being calculated is 0.038ppm, is less than national standard
0.07ppm.Therefore the selection of injected volume is reasonable.
3.5 launch tracer to throwing in source point, if throwing source point there are multiple water-bearing layers, using top layer water-bearing layer as showing
The dispensing target zone of track agent.
As a kind of specific embodiment, Xing Jihe described in the present embodiment only includes one, Lower Ordovician limestone and contains
Water layer, therefore Lower Ordovician limestone is the target zone that tracer is launched.
When practical operation, pipeline is goed deep into throwing in source point, and the lower port of pipeline is made to be located at target zone, then pass through pipe
Tracer is launched to target zone in road.
The water levels of 3.6 observed stage monitoring points and Water-quality Monitoring Points observes three days once a day, and according to observing
Water level draws water level contour map, then observes the water level contour map drawn, judges to whether there is near Water-quality Monitoring Points
Cone of depression, cone of depression, then be marked the Water-quality Monitoring Points if it exists.
The main reason for doing so be, the mining of groundwater if it exists near Water-quality Monitoring Points is then likely to form landing
Funnel, in this way the tracer detection time of the Water-quality Monitoring Points and set tracer detection time have difference, general tracer
The time that agent peak concentration occurs can shift to an earlier date, and need to carry out encryption monitoring to the Water-quality Monitoring Points.So after tracer dispensing
Will water level unified test to Water-quality Monitoring Points three times, in order to judge whether there is cone of depression, avoid missing tracer
Peak concentration.
3.7 are sampled monitoring to Water-quality Monitoring Points according to sequence from the near to the distant, record the concentration and appearance of tracer
Time, and according in monitoring data and step 3.6 determine cone of depression situation sampling rate adjusting.That is:
Increase sample frequency when Water-quality Monitoring Points detect tracer, and the Water-quality Monitoring Points to next group remotely simultaneously
It is monitored.
When the tracer concentration of Water-quality Monitoring Points detection is lower than peak value, sample frequency is reduced.
When the tracer concentration of Water-quality Monitoring Points detection reaches background value, sampling terminates.
There are when cone of depression near Water-quality Monitoring Points, encryption monitoring will be carried out to the Water-quality Monitoring Points.I.e. upper one
It, need to be to the Water-quality Monitoring Points (nearby in the presence of landing when criticizing (distance throws the closer a batch of source point) Water-quality Monitoring Points detection tracer
The Water-quality Monitoring Points of funnel) when being sampled monitoring, the Water-quality Monitoring Points (Water-quality Monitoring Points that nearby there is cone of depression)
Sampling frequency is just higher than other Water-quality Monitoring Points that cone of depression is nearby not present.If (nearby there is drop in the Water-quality Monitoring Points
Fall the Water-quality Monitoring Points of funnel) it is the Water-quality Monitoring Points that first is monitored, then when being monitored at the beginning, the water quality monitoring
The sampling frequency of point (Water-quality Monitoring Points that nearby there is cone of depression) is just higher than other water that cone of depression is nearby not present
Quality supervision measuring point.
Herein, the detection standard of tracer is 5 times of background value concentration, i.e., when the tracer concentration of detection reaches background
Tracer is just detected at last at 5 times of value.
As a kind of specific embodiment, the tracer experiment in the present embodiment is tried in dispensing ammonium molybdate on April 18th, 2016
Agent 280kg lasts 95 days altogether, and sampling monitoring sequence is as follows:
(1) analysis is sampled to J3, J5, FJ34, FJ65 quality supervision measuring point first, frequency is once a day.
(2) when J3, J5, FJ34, FJ65 Water-quality Monitoring Points in step (1) have tracer detection, to J3, J5, FJ34,
The monitoring frequencies of FJ65 Water-quality Monitoring Points increases for twice a day, and simultaneously to J1, J2, J4, J6, J7, Q1, Q2, Q3, Q4 and Q5
Water-quality Monitoring Points are sampled, and sampling frequency is once a day.
Wherein J4 Water-quality Monitoring Points are one day due to nearby there is cone of depression, the monitoring frequency of J4 Water-quality Monitoring Points
Twice.
(3) when J1, J2, J6, J7, Q1, Q2, Q3, Q4 and Q5 Water-quality Monitoring Points in step (2) have tracer detection,
The monitoring frequency of J1, J2, J6, J7, Q1, Q2, Q3, Q4 and Q5 Water-quality Monitoring Points is increased as twice a day.
(4) tracer concentration curve is drawn according to the data of step (1), (2) and (3) monitoring, and according to tracer concentration
Curve carries out observation judgement, when the tracer concentration of Water-quality Monitoring Points detection is lower than peak concentration, sampling frequency is reduced, by one
It is changed to once a day twice, and when the tracer concentration of Water-quality Monitoring Points detection reaches background value concentration, sampling monitoring terminates.
Step (4) carries out simultaneously with step (1), (2) and (3).
3.8 calculate horizontal view flow velocity
Tracer diffusion velocity isochrone and dust trajectory are drawn according to the result that step 3.7 monitors, as shown in Fig. 2, so
Horizontal view flow velocity is calculated afterwards.
Level refers to resulting divided by the time of tracer first meeting peak value to monitoring point linear distance from source point is thrown depending on flow velocity
Value.
As a kind of specific embodiment, the horizontal view flow velocity being calculated in the present embodiment is as shown in table 2.
2 tracer test of table regards flow velocity list
3.9 calculate the mean flow rate in survey region according to the calculated result of step 3.8, i.e., calculate in step 3.8
To flow velocity average.
Then flow velocity is greater than there may be seepage channels in the region of mean flow rate.
As a kind of specific embodiment, the mean flow rate in the present embodiment in Jinan Spring Region runoff region is
91.64m/d。
4th step calculates vertical flow velocity, determines the vertical covered depth of seepage channel.
4.1 choose the Water-quality Monitoring Points that flow velocity is greater than mean flow rate according to the calculated result of step 3.9, and supervise in the water quality
Every 1m monitoring temperature and conductivity in drilling on measuring point.
The principle that monitoring frequency follows are as follows:
It is monthly tested 3 times in 2-7 month, 9-December, in the most cold January moon and most hot moon August increase test frequency is every
Moon test 5 times；Test 3 days is encrypted after rain, once a day.
Herein, the borehole test section of casing is located at casing shoe to bottom hole, uncased borehole test section is located at water
Face is to bottom hole.In addition, when instrument readings, 3 minutes are needed to wait for after numerical stability for the accuracy for guaranteeing test result, then
Record need to take when recording data to 2 significant digits.
4.2 draw the temperature profile and conductance profile figure of each Water-quality Monitoring Points according to the test result of step 4.1.
As a kind of specific embodiment, by taking Water-quality Monitoring Points J2 as an example, the temperature curve of the Water-quality Monitoring Points J2 of drafting
Figure and conductance profile figure are as shown in Figure 3 and Figure 4.
4.3 calculate the groundwater velocity on vertical for a certain Water-quality Monitoring Points.
4.3.1 the temperature profile and conductance profile figure drawn according to step 4.2 determine the temperature of the Water-quality Monitoring Points
The depth that inflection point and conductivity inflection point occur, determines the covered depth of horizontal seepage channel at the Water-quality Monitoring Points.
4.3.2 the vertical flow velocity of the underground water in the depth that step 4.3.1 is determined is calculated.
In order to facilitate understanding, now by taking J2 as an example, the calculation method of underground water vertical flow speed is described in detail.
Such as Fig. 3 and Fig. 4 it is found that the depth that J2 temperature inflection point and conductivity inflection point occur jointly be respectively 40m, 60m and
90m.Then there is the seepage channel of 3 depth segments at J2 in the horizontal direction, respectively at 40m, 60m and 90m.
The depth occurred jointly according to temperature inflection point and conductivity inflection point divide computation interval be 25m-51m, 53m-71m and
70m-100m。
According to formulaThe value of f (β, Z/L) can be obtained.
Then, Z/L-f (β, Z/L) typical curve that Bredehoeft and Papadopulos is summarized is looked into, β value can be obtained.β's
Positive and negative positive and negative (underground water, which flows downward, to be positive, and flows up and is negative) depending on Vz.
Then β value is brought into formulaIn, V can be acquiredZ。
In formula: VZ- fluid is in the vertical flow velocity cm/s that depth is at Z；
C0The specific heat capacity of-fluid takes 1cal.g under room temperature-1.℃-1；
ρ0The density of-fluid takes 1g.cm under room temperature-3；
κ-solidliquid mixture pyroconductivity；
T0The temperature of coboundary in-computation interval；
TZThe temperature of lower boundary in-computation interval；
Spacing in L-computation interval between up-and-down boundary；
Distance of the Z-unknown point to computation interval coboundary；
TZThe temperature of-unknown point.
According to above-mentioned calculation method, the vertical flow velocity of underground water that J2 is calculated is as shown in table 3.
The vertical flow velocity calculated result statistical form of 3 J2 underground water of table
Herein, why to calculate underground water it is vertical on flow velocity, be it is vertical upper there is no seepage channel in order to prove,
Belong to runoff area (on vertical there are the region of seepage channel be feeding area).
By upper table analysis it is found that the vertical flow velocity of Water-quality Monitoring Points J2 is -6.67 × 10-7cm/s-3.00×10- 6Between cm/s, there is no the seepage channels on vertical, belong to runoff area.
Determining the theoretical basis of the covered depth of seepage channel by this way is, seeps when drilling through underground water in stratum
When circulation road, influenced by underground water bottom horizontal flow sheet, temperature distribution history will appear inflection point.As shown in figure 5, working as in stratum without ground
When being lauched seepage flow, Temperature Distribution is only related with depth, linearly increasing with depth increase, reflects normal ground temperature distribution.It is local
There are when seepage channel in layer, subsurface flow can disturb the distribution of stratum normal temperature in seepage channel, from the different of temperature curve
It often can be inferred that the information of seepage action of ground water in stratum.As shown in fig. 6, by analysis it is found that 30-50m are to drill through ground
Seepage action of ground water channel in layer, and the temperature of seepage water is lower, forms temperature distribution history spatial abnormal feature recessed down；
The temperature of seepage water is higher at 50-70m, the upward spatial abnormal feature of formation temperature distribution curve, seepage flow between 90-110m
The temperature of water is lower；Form temperature distribution history spatial abnormal feature recessed down.
Herein why simultaneously measuring temperature and conductivity, be in order to ensure reliability, while measuring temperature and conductance
Rate improves the precision and accuracy of this method.When the development of water-bearing layer Karst Fissures, ground water electric conductivity changing sensitivity drop
It is low, therefore temperature can supplement conductivity.
4.4 repeat the operation of step 4.3, calculate the underground water vertical flow of each Water-quality Monitoring Points determined in step 4.1
Speed.
5th step is superimposed horizontal view flow velocity and vertical temperature, conductivity inflection point, draws and seep on the basis of geologic section
Flow channel space is distributed log sheet.
As a kind of specific embodiment, the Jinan Spring Region runoff area seepage channel spatial distribution drawn in the present embodiment
Matter sectional view is as shown in Figure 7.
By the analysis to Fig. 7 it is found that regarding flow velocity between Xing Jihe and springs between 109.19-112.98m/d, hang down
There is the seepage channel between springs between upward 120m to -160m.
Claims (8)
1. a kind of method for determining the horizontal seepage channel in runoff area using tracer technique, it is characterised in that: include the following steps,
The first step investigates the structure of drillings all in survey region, determines hole depth, water-bearing layer according to hydrogeologic map
Lithology, sealing section depth, borehole coordinate and elevation；
Second step observes the water level in respectively drilling in the first step one by one, draws the ground water field figure in survey region, determines each
Hydraulic connection between drilling；
Third step arranges tracer test, calculates horizontal view flow velocity, and according to the horizontal horizontal extension for determining seepage channel depending on flow velocity
Direction；
3.1 determine throwing source point according to the ground water field figure that second step is drawn, and centered on throwing source point, surrounding arranges water-level observation
Point, centered on throwing source point, in the downstream of direction of groundwater flow, fan-shaped arrangement water-quality observation point；
3.2 choose tracer；
The Water-quality Monitoring Points determined in 3.3 pairs of steps 3.1 are sampled, and determine the tracer background value of each Water-quality Monitoring Points；
3.4 calculate tracer injected volume；
3.5 launch tracer into throwing source point；
The water level of 3.6 observed stage monitoring points and Water-quality Monitoring Points observes three days, and according to the water level observed once a day
Water level contour map is drawn, the water level contour map drawn is then observed, is judged near Water-quality Monitoring Points with the presence or absence of landing
Funnel, cone of depression, then be marked the Water-quality Monitoring Points if it exists.
3.7 are sampled monitoring to Water-quality Monitoring Points according to sequence from the near to the distant, record tracer concentration and appearance when
Between, and according to the cone of depression situation sampling rate adjusting determined in monitoring data and step 3.6；
3.8 draw tracer diffusion velocity isochrone and dust trajectory according to the result that step 3.7 monitors, and then calculate horizontal view
Flow velocity；
3.9 according to the mean flow rate in the calculating Structure Calculation survey region of step 3.8；
4th step calculates vertical flow velocity, determines the vertical covered depth of seepage channel.
4.1 choose the Water-quality Monitoring Points that flow velocity is greater than mean flow rate according to the calculated result of step 3.9, and in the Water-quality Monitoring Points
On drilling in every 1m monitoring temperature and conductivity；
4.2 draw the temperature profile and conductance profile figure of each Water-quality Monitoring Points according to the monitoring result of step 4.1；
4.3 calculate the groundwater velocity on vertical for a certain Water-quality Monitoring Points；
4.3.1 the temperature profile and conductance profile figure drawn according to step 4.2 determine the temperature inflection point of the Water-quality Monitoring Points
The depth occurred with conductivity inflection point, determines the covered depth of horizontal seepage channel at the Water-quality Monitoring Points；
4.3.2 the vertical flow velocity of the underground water in each inflection point depth that step 4.3.1 is determined is calculated；
4.4 repeat the operation of step 4.3, calculate the vertical flow velocity of underground water of each Water-quality Monitoring Points determined in step 4.1；
5th step, according to third step and the 4th step as a result, on the basis of geologic section, with drawing seepage channel spatial distribution
Matter sectional view.
2. a kind of method for determining the horizontal seepage channel in runoff area using tracer technique according to claim 1, feature
Be: the water level monitoring point and Water-quality Monitoring Points of laying should all be comprising all water-bearing layers in survey region.
3. a kind of method for determining the horizontal seepage channel in runoff area using tracer technique according to claim 1, feature
Be: the tracer should meet nontoxic, chemical property is stable, be not easy to be adsorbed and degrade, solubility is high, background value is low and
The requirement of high sensitivity.
4. a kind of method for determining the horizontal seepage channel in runoff area using tracer technique according to claim 1, feature
Be: in step 3.3, same Water-quality Monitoring Points are sampled 2-3 times, and the average value of the tracer concentration to measure is as the water
The background value concentration of quality supervision measuring point.
5. a kind of method for determining the horizontal seepage channel in runoff area using tracer technique according to claim 1, feature
Be: in step 3.5, if throwing source point, there are multiple water-bearing layers, the dispensing purpose using top layer water-bearing layer as tracer
Layer.
6. a kind of method for determining the horizontal seepage channel in runoff area using tracer technique according to claim 1, feature
It is: in step 3.7, when being detected to the tracer of Water-quality Monitoring Points, using 5 times of background value concentration as tracer
Detection standard.
7. a kind of method for determining the horizontal seepage channel in runoff area using tracer technique according to claim 1, feature
Be: in step 3.7, the Adjustment principle of sample frequency is,
Increase sample frequency when Water-quality Monitoring Points detect tracer, and the Water-quality Monitoring Points of next group remotely are carried out simultaneously
Monitoring；
When the tracer concentration of Water-quality Monitoring Points detection is lower than peak value, sample frequency is reduced；
When the tracer concentration of Water-quality Monitoring Points detection reaches background value, sampling terminates；
There are when cone of depression near Water-quality Monitoring Points, encryption monitoring will be carried out to the Water-quality Monitoring Points.
8. a kind of method for determining the horizontal seepage channel in runoff area using tracer technique according to claim 1, feature
Be: the principle that monitoring frequency follows in step 4.1 is monthly to test 3 times in 2-7 month, 9-December, in January and August
Increasing test frequency is monthly to test 5 times, encryption test 3 days after rain, once a day.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811549627.7A CN109633764B (en) | 2018-12-18 | 2018-12-18 | Method for determining horizontal seepage channel of runoff area by adopting tracing technology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811549627.7A CN109633764B (en) | 2018-12-18 | 2018-12-18 | Method for determining horizontal seepage channel of runoff area by adopting tracing technology |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109633764A true CN109633764A (en) | 2019-04-16 |
| CN109633764B CN109633764B (en) | 2020-03-31 |
Family
ID=66075042
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811549627.7A Active CN109633764B (en) | 2018-12-18 | 2018-12-18 | Method for determining horizontal seepage channel of runoff area by adopting tracing technology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109633764B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110874976A (en) * | 2019-11-29 | 2020-03-10 | 济南大学 | Method for simulating dynamic state of underground water of karst big spring |
| CN110927816A (en) * | 2019-12-11 | 2020-03-27 | 中国地质科学院岩溶地质研究所 | Detection method of karst underground river system |
| CN111397663A (en) * | 2020-04-20 | 2020-07-10 | 河南省地质矿产勘查开发局第二地质环境调查院 | Method for judging development degree of karst according to temperature of underground water |
| CN111929746A (en) * | 2020-07-28 | 2020-11-13 | 山东科技大学 | Detection method for exploring development degree of collapse column based on tracing technology |
| CN113341479A (en) * | 2021-05-24 | 2021-09-03 | 国电建投内蒙古能源有限公司 | Tracing test method for determining hydraulic connectivity of mining area |
| CN113740225A (en) * | 2021-08-23 | 2021-12-03 | 西安交通大学 | Extraction method of porous medium micro seepage path based on particle tracing |
| CN113899806A (en) * | 2021-09-01 | 2022-01-07 | 中国地质科学院水文地质环境地质研究所 | Novel artificial tracer suitable for putting and tracing analysis method thereof |
| CN114659734A (en) * | 2022-03-02 | 2022-06-24 | 河海大学 | Method for detecting dam leakage by combining high-density electrical method and comprehensive tracing method |
| CN115375203A (en) * | 2022-10-25 | 2022-11-22 | 山东省地质矿产勘查开发局八〇一水文地质工程地质大队(山东省地矿工程勘察院) | Comprehensive analysis system for multi-element water resource optimization configuration |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3678138A (en) * | 1969-04-03 | 1972-07-18 | Imp Metal Ind Kynoch Ltd | Solid propellant charge making in mold having perforated separator means |
| CN1344911A (en) * | 2001-11-09 | 2002-04-17 | 河海大学 | Generalized tracing and diluting method for monoporate penetration flow rate measurement |
| CN103164569A (en) * | 2013-01-31 | 2013-06-19 | 尹立河 | Determination method of different depth ground water circulation volumes of big size basin |
| CN103197039A (en) * | 2013-03-15 | 2013-07-10 | 陈秋阳 | Field underground water dispersion test method in artificial flow field |
| CN104282214A (en) * | 2014-10-29 | 2015-01-14 | 安徽理工大学 | Pipeline flow tracer test comprehensive device of pressure bearing karst aquifer system |
| CN104696010A (en) * | 2015-03-20 | 2015-06-10 | 青岛理工大学 | Comprehensive treatment method of mine water inrush |
| CN108020649A (en) * | 2017-11-22 | 2018-05-11 | 济南大学 | A kind of method of definite karst big spring feed channel and intensity |
| CN108106687A (en) * | 2018-02-08 | 2018-06-01 | 济南大学 | A kind of bedrock underground water drift net containing soft interlayer probes into method and double-capsule water-stop |
| CN108343431A (en) * | 2018-02-08 | 2018-07-31 | 济南大学 | A kind of bedrock fracture ground water regime drift net probes into method |
-
2018
- 2018-12-18 CN CN201811549627.7A patent/CN109633764B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3678138A (en) * | 1969-04-03 | 1972-07-18 | Imp Metal Ind Kynoch Ltd | Solid propellant charge making in mold having perforated separator means |
| CN1344911A (en) * | 2001-11-09 | 2002-04-17 | 河海大学 | Generalized tracing and diluting method for monoporate penetration flow rate measurement |
| CN103164569A (en) * | 2013-01-31 | 2013-06-19 | 尹立河 | Determination method of different depth ground water circulation volumes of big size basin |
| CN103197039A (en) * | 2013-03-15 | 2013-07-10 | 陈秋阳 | Field underground water dispersion test method in artificial flow field |
| CN104282214A (en) * | 2014-10-29 | 2015-01-14 | 安徽理工大学 | Pipeline flow tracer test comprehensive device of pressure bearing karst aquifer system |
| CN104696010A (en) * | 2015-03-20 | 2015-06-10 | 青岛理工大学 | Comprehensive treatment method of mine water inrush |
| CN108020649A (en) * | 2017-11-22 | 2018-05-11 | 济南大学 | A kind of method of definite karst big spring feed channel and intensity |
| CN108106687A (en) * | 2018-02-08 | 2018-06-01 | 济南大学 | A kind of bedrock underground water drift net containing soft interlayer probes into method and double-capsule water-stop |
| CN108343431A (en) * | 2018-02-08 | 2018-07-31 | 济南大学 | A kind of bedrock fracture ground water regime drift net probes into method |
Non-Patent Citations (1)
| Title |
|---|
| 万继伟等: "引汉济渭工程黄三隧洞地下水渗流同位素综合示踪探测分析研究", 《水利水电技术》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110874976B (en) * | 2019-11-29 | 2021-10-22 | 济南大学 | Method for simulating dynamic state of underground water of karst big spring |
| CN110874976A (en) * | 2019-11-29 | 2020-03-10 | 济南大学 | Method for simulating dynamic state of underground water of karst big spring |
| CN110927816A (en) * | 2019-12-11 | 2020-03-27 | 中国地质科学院岩溶地质研究所 | Detection method of karst underground river system |
| CN110927816B (en) * | 2019-12-11 | 2021-08-24 | 中国地质科学院岩溶地质研究所 | Detection method of karst underground river system |
| CN111397663A (en) * | 2020-04-20 | 2020-07-10 | 河南省地质矿产勘查开发局第二地质环境调查院 | Method for judging development degree of karst according to temperature of underground water |
| CN111929746A (en) * | 2020-07-28 | 2020-11-13 | 山东科技大学 | Detection method for exploring development degree of collapse column based on tracing technology |
| CN113341479A (en) * | 2021-05-24 | 2021-09-03 | 国电建投内蒙古能源有限公司 | Tracing test method for determining hydraulic connectivity of mining area |
| CN113341479B (en) * | 2021-05-24 | 2024-05-07 | 国电建投内蒙古能源有限公司 | Trace test method for determining hydraulic connectivity of goaf |
| CN113740225A (en) * | 2021-08-23 | 2021-12-03 | 西安交通大学 | Extraction method of porous medium micro seepage path based on particle tracing |
| CN113899806A (en) * | 2021-09-01 | 2022-01-07 | 中国地质科学院水文地质环境地质研究所 | Novel artificial tracer suitable for putting and tracing analysis method thereof |
| CN114659734A (en) * | 2022-03-02 | 2022-06-24 | 河海大学 | Method for detecting dam leakage by combining high-density electrical method and comprehensive tracing method |
| CN114659734B (en) * | 2022-03-02 | 2022-12-02 | 河海大学 | Method for detecting dam leakage by combining high-density electrical method and comprehensive tracing method |
| CN115375203A (en) * | 2022-10-25 | 2022-11-22 | 山东省地质矿产勘查开发局八〇一水文地质工程地质大队(山东省地矿工程勘察院) | Comprehensive analysis system for multi-element water resource optimization configuration |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109633764B (en) | 2020-03-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109633764A (en) | A method of the horizontal seepage channel in runoff area is determined using tracer technique | |
| WO2019233105A1 (en) | Device and method for measuring flow rate, flow direction, and geological parameter of deep-well cross-hole groundwater | |
| Andersen et al. | Stream-aquifer interactions in the Maules creek catchment, Namoi valley, New South Wales, Australia | |
| CN106437844A (en) | Method for advanced forecast of tunnel water inflow position | |
| CN101476463B (en) | Drill-following natural gamma geosteering method for horizontal well | |
| CN113820078B (en) | Method for detecting dam leakage by combining tracing method and dilution concentration method | |
| CN107289997A (en) | A kind of Karst-fissure water detection system and method | |
| Drogue | Geothermal gradients and ground water circulation in fissured and karstic rocks: the role played by the structure of the permeable network | |
| CN109707373A (en) | A kind of two-way tracing method of horizontal well-straight well based on producing profile testing and inter-well tracer test | |
| CN107703551A (en) | A kind of great burying concrete drain tile line detection method with pressure | |
| CN108241772A (en) | Consider Tunnel Gushing method for determination of amount in the crack artesian aquifer of multiple factors | |
| CN114659734A (en) | Method for detecting dam leakage by combining high-density electrical method and comprehensive tracing method | |
| Feng‐Shan et al. | APPLICATION OF HIGH‐DENSITY RESISTIVITY METHOD IN DETECTING WATER‐BEARING STRUCTURES AT A SEABED GOLD MINE | |
| Guanghui et al. | Hydrogeological characteristics of foot caves in a karst peak-forest plain in South China | |
| CN111912540B (en) | Method for testing and predicting geothermal anomaly tunnel engineering ground temperature in complex and difficult highland mountain areas | |
| CN109143381B (en) | A kind of loess high roadbed aquitard detection method | |
| Zhang et al. | Why mixed groundwater at the outlet of open flowing wells in unconfined-aquifer basins can represent deep groundwater: implications for sampling in long-screen wells. | |
| CN112505792B (en) | Salt lake deep brine zone exploration method based on radon gas indication | |
| CN109187280A (en) | A kind of comprehensive judgement method of injection recovery technique dilation angle in situ | |
| Dapaah‐Siakwan et al. | Estimation of vertical water and heat fluxes in the semi‐confined aquifers in Tokyo metropolitan area, Japan | |
| CN109164018A (en) | The continuous monitor system and monitoring method of injection recovery technique dilation angle in situ | |
| CN111963152B (en) | Ground temperature gradient calculation method based on horizontal hole layered temperature measurement data | |
| CN108343431A (en) | A kind of bedrock fracture ground water regime drift net probes into method | |
| CN114087019A (en) | Method for preventing and controlling water damage area of thick and non-homogeneous sandstone | |
| CN112129681A (en) | Method for acquiring change rule of water content of water-rich soft plastic loess tunnel under seepage effect |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |