CN108982317A - A kind of determination method of high roadbed soil body large scale infiltration coefficient - Google Patents
A kind of determination method of high roadbed soil body large scale infiltration coefficient Download PDFInfo
- Publication number
- CN108982317A CN108982317A CN201810613729.4A CN201810613729A CN108982317A CN 108982317 A CN108982317 A CN 108982317A CN 201810613729 A CN201810613729 A CN 201810613729A CN 108982317 A CN108982317 A CN 108982317A
- Authority
- CN
- China
- Prior art keywords
- coefficient
- infiltration
- large scale
- scale
- soil body
- 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
- 230000008595 infiltration Effects 0.000 title claims abstract description 167
- 238000001764 infiltration Methods 0.000 title claims abstract description 167
- 239000002689 soil Substances 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000004088 simulation Methods 0.000 claims abstract description 63
- 238000012360 testing method Methods 0.000 claims abstract description 61
- 230000035699 permeability Effects 0.000 claims abstract description 51
- 125000002619 bicyclic group Chemical group 0.000 claims abstract description 43
- 238000009826 distribution Methods 0.000 claims abstract description 39
- 238000005520 cutting process Methods 0.000 claims abstract description 27
- 230000000694 effects Effects 0.000 claims abstract description 14
- 230000003204 osmotic effect Effects 0.000 claims abstract description 10
- 238000009738 saturating Methods 0.000 claims description 3
- 230000003321 amplification Effects 0.000 claims 1
- 238000003199 nucleic acid amplification method Methods 0.000 claims 1
- 239000003673 groundwater Substances 0.000 abstract description 5
- 238000011065 in-situ storage Methods 0.000 abstract description 3
- 238000011160 research Methods 0.000 abstract description 3
- 238000007619 statistical method Methods 0.000 abstract description 3
- 238000005422 blasting Methods 0.000 description 15
- 238000005070 sampling Methods 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000001186 cumulative effect Effects 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000013138 pruning Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of determination methods of high roadbed soil body large scale infiltration coefficient, comprising: obtains and fills soil body osmotic coefficient on cutting ring scale in the changing rule in space, anisotropy and heterogeneity；The bicyclic infiltration test in field is carried out at high roadbed soil body different location, obtains the permeability coefficient on cutting ring scale；In conjunction with scale effect principle, analogue simulation is carried out to the bicyclic infiltration test in field；The bicyclic infiltration test analogue simulation in field is amplified in large scale and is simulated, in conjunction with the original distribution for filling soil body osmotic coefficient on cutting ring scale, obtains the distribution of large scale infiltration coefficient.Required large scale infiltration coefficient when the Groundwater movement rule following the present invention relates to the strong high roadbed soil body needs research prediction of heterogeneity, i.e. by largely the methods of test, statistical analysis, condition simulation, scale effect and means, the large scale infiltration coefficient for the high roadbed soil body that large-scale in-situ test of having analyzed and researched is difficult to carry out and heterogeneity is high.
Description
Technical field
The present invention relates to soil body osmotic coefficient technical fields, more particularly relate to a kind of high roadbed soil body large scale infiltration system
Several determination methods.
Background technique
Loess is widely distributed in China, and since loess sheet has the characteristic easily corroded, hills gully becomes its mainly
Looks type, even if there is also a large amount of large-scale coombes in Loess Tableland, and the scale of coombe is also constantly expanding.In recent years
Come, with the development of economy and society, land resource day is becoming tight.In order to meet needs, Canal in Loess Area starts to occur extensive
The engineering for filling out ditch level land, produce a large amount of high roadbed loess filling body.The change gesture of landform shape in great geotechnical engineering
The characteristics of motion of underground water must be influenced, numerical simulation is to predict the effective way of ground water movement rule under future condition.?
In numerical simulation, the accuracy of infiltration coefficient is to influence the key factor of calculated result.
Largely studies have shown that infiltration coefficient can change with the difference of Research scale, that is, there is certain ruler
Spend effect.Schulze-Makuch etc. demonstrates infiltration system by the data of 39 kinds of geologic medium infiltration coefficients and test scale
Number increases with observing result and the rules that increase and tend towards stability, and before tending towards stability infiltration coefficient and observing result right
It is in a linear relationship on number coordinate value.The determination method of infiltration coefficient mainly has permeameter method, bicyclic infiltration method, bailing test at present
Method, water injection test method etc..However, bailing test and water injection test are often difficult to carry out, and the number in region since condition limits
Value simulation is built upon on the basis of large scale subdivision grid, and indoor permeability test is carried out on the scale of very little, volume
Only 120cm3, calculated result is difficult to represent the infiltration coefficient in large scale.In addition, landfill area's loess heterogeneity is stronger,
The infiltration test of field minority can not react the distribution of large area landfill area's infiltration coefficient comprehensively.So how to determine high roadbed
The infiltration coefficient of soil body large scale becomes the critical issue for studying such groundwater simulation.
In conclusion infiltration coefficient in the prior art determines method, high roadbed soil body large scale can not be determined by existing
The problem of infiltration coefficient.
Summary of the invention
The embodiment of the present invention provides a kind of determination method of high roadbed soil body large scale infiltration coefficient, to solve existing skill
There are problems that not determining the infiltration coefficient of high roadbed soil body large scale in art.
The embodiment of the present invention provides a kind of determination method of high roadbed soil body large scale infiltration coefficient, comprising: in high roadbed
Soil body different level position and vertical position are sampled, and with permeameter method is vertical to soil sample and horizontal infiltration coefficient into
Row measurement obtains and fills soil body osmotic coefficient on cutting ring scale in the changing rule in space, anisotropy and heterogeneity；
The bicyclic infiltration test in field is carried out at high roadbed soil body different location, obtains the infiltration coefficient on cutting ring scale
Value；
The bicyclic infiltration test in field is imitated in conjunction with scale effect principle according to the permeability coefficient on cutting ring scale
True simulation；
The bicyclic infiltration test analogue simulation in field is amplified in large scale and is simulated, it is native in conjunction with being filled on cutting ring scale
Body infiltration coefficient obtains the distribution of large scale infiltration coefficient in the changing rule in space, anisotropy and heterogeneity.
The embodiment of the present invention provides a kind of determination method of high roadbed soil body large scale infiltration coefficient, further includes:
Bicyclic infiltration test simulation model is established, the top layer of simulation model is arranged using square volume mesh for uniformly subdivision
To determine head, head value is 0m, will be also configured as determining head within the scope of simulation model bottom inner ring and outer rings, head value is equal to imitative
Distance of the true mode top layer to bottom；
In conjunction with scale effect principle, chooses certain infiltration coefficient typical value, by isotropic medium, calculate different subdivisions
Different random is distributed resulting permeability coefficient under scale, and compares with the infiltration coefficient measured value of the bicyclic infiltration test in field,
Choose suitable subdivision size；
On determining subdivision size basis, large scale equivalent coefficient of permeability calculated result is analyzed to different infiltration coefficient generations
Tabular value quantity and anisotropic sensibility.
Further, the described bicyclic infiltration test analogue simulation in field is amplified in large scale is simulated, coupling collar
Soil body osmotic coefficient is filled on knife scale in the changing rule in space, anisotropy and heterogeneity, obtains large scale infiltration system
Several distributions；Include:
On determining subdivision size basis, choose certain typical value quantity, by isotropic medium by simulation context
Gradually amplify, the infiltration coefficient under random distribution different under each simulation context is calculated, by permeability coefficient with simulation context
Simulation context value of increase when tending towards stability be determined as the simulation context of large scale infiltration coefficient；
Large scale equivalent coefficient of permeability calculated result is analyzed to the sensibility of the quantity of infiltration coefficient typical value, it is suitable to determine
Typical value number；
Large scale equivalent coefficient of permeability calculated result is analyzed to anisotropic sensibility, determines tested out cutting ring scale
Large scale permeability coefficient under infiltration coefficient large sample；
According to sample and overall relationship, the distribution of large scale infiltration coefficient is determined.
In the embodiment of the present invention, a kind of determination method of high roadbed soil body large scale infiltration coefficient is provided, with the prior art
It compares, its advantages are as follows:
When needing to study prediction future Groundwater movement rule the present invention relates to a kind of high roadbed soil body that heterogeneity is strong
Required large scale infiltration coefficient, by largely testing, statisticalling analyze, the methods of condition simulation, scale effect and means,
The large scale infiltration coefficient for the high roadbed soil body that large-scale in-situ test of having analyzed and researched is difficult to carry out and heterogeneity is high, big ruler
Degree infiltration coefficient is the important parameter of forecast analysis Groundwater movement rule, and the medium-and-large-sized field in-situ test of the high roadbed soil body is past
It is past to be difficult to carry out, therefore the invention proposes a kind of high roadbed soil bodys based on the bicyclic infiltration test of results of permeability tests and field
The determination method of large scale infiltration coefficient, since the high roadbed soil body has the strong feature of heterogeneity, by largely sampling point
Analysis, obtains the infiltration coefficient of cutting ring scale, statisticallys analyze to cutting ring size infiltration coefficient, in conjunction with scale effect principle, to field
Bicyclic infiltration test carries out analogue simulation, on the basis of gradually analyzing each factor for influencing analog result, will simulate model
It encloses and is expanded within the scope of large scale, according to sample and overall relationship, determine the distribution of large scale infiltration coefficient.
Detailed description of the invention
Fig. 1 is a kind of determination method flow diagram of high roadbed soil body large scale infiltration coefficient provided in an embodiment of the present invention；
Fig. 2 is sample position figure provided in an embodiment of the present invention；
Fig. 3 is the native same position different directions infiltration coefficient comparison diagram of landfill provided in an embodiment of the present invention；
Fig. 4 is the landfill soil same position same direction infiltration coefficient comparison diagram that inventive embodiments provide；
Fig. 5 is the different sampling number vertical permeability coefficient box traction substations that inventive embodiments provide；
Fig. 6 is the different sampling number horizontal-hole blasting box traction substations that inventive embodiments provide；
Fig. 7 is the different sampling high perpendicular infiltration coefficient box traction substations that inventive embodiments provide；
Fig. 8 is different sampling height level's infiltration coefficient box traction substations that inventive embodiments provide；
Fig. 9 is the vertical permeability coefficient histogram that inventive embodiments provide；
Figure 10 is the horizontal-hole blasting histogram that inventive embodiments provide；
Figure 11 is the cumulative frequency figure for the infiltration coefficient that inventive embodiments provide；
Infiltration coefficient random distribution figure within the scope of the bicyclic test simulation that Figure 12 provides for inventive embodiments；
Figure 13 is that the bicyclic test model that inventive embodiments provide determines head boundary profile figure；
Figure 14 is that the model that inventive embodiments provide calculates and the bicyclic infiltration test comparative result figure in field；
Figure 15 is bicyclic infiltration test calculated result comparison diagram under the anisotropy various combination that inventive embodiments provide；
Figure 16 is the bicyclic infiltration test calculated result comparison diagram of different number typical value that inventive embodiments provide；
Figure 17 is the different simulation context large scale infiltration coefficient calculated result comparison diagrams that inventive embodiments provide；
Figure 18 is large scale infiltration coefficient calculated result box traction substation under the different typical value quantity that inventive embodiments provide；
Figure 19 is large scale infiltration coefficient calculated result comparison diagram under the anisotropy various combination that inventive embodiments provide；
Figure 20 is the large scale infiltration coefficient calculated result distribution map that inventive embodiments provide.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Referring to Fig. 1, a kind of determination method of high roadbed soil body large scale infiltration coefficient provided in an embodiment of the present invention, the party
Method includes:
Step 1, it is sampled in high roadbed soil body different level position and vertical position, and with permeameter method to soil
Sample is vertical and horizontal infiltration coefficient is measured, obtain landfill soil body osmotic coefficient on cutting ring scale space changing rule,
Anisotropy and heterogeneity.
It should be noted that in order to analyze the changing rule in infiltration coefficient space on cutting ring scale, anisotropy, heterogeneous
Property, the high roadbed soil body is largely sampled, same position sampling number should be 2-3, and sample position is in entire space average
Distribution.
Step 2, the bicyclic infiltration test in field is carried out at high roadbed soil body different location, obtains the infiltration on cutting ring scale
Coefficient value.
It should be noted that can be chosen in space a small number of due to the bicyclic permeability test in field since workload is bigger
Representative point tested.
Step 3, according to the permeability coefficient on cutting ring scale, in conjunction with scale effect principle, to the bicyclic infiltration test in field
Carry out analogue simulation.
It should be noted that the bicyclic infiltration test analogue simulation in field, is square grid by simulation uniformly subdivision, it is square
The infiltration coefficient of volume mesh is chosen according to the infiltration coefficient of cutting ring scale, the distribution of infiltration coefficient in different grids, according to a large amount of
Cutting ring scale infiltration coefficient statistical result determine.According to the principle of pit permeability test, set fixed for the top layer of model
Head, head value are 0m；It is also configured as determining head within the scope of model bottom inner ring and outer rings, head value is equal to model top layer on earth
The distance of layer, with guarantee bottom to top layer the average current gradient for 1.The size of square and its representative infiltration coefficient are big
It is small very big on analog result influence.It in conjunction with the scale effect of infiltration coefficient, is verified, is come with the bicyclic permeability test result in field
Determine the suitable subdivision size of square.
Step 4, the bicyclic infiltration test analogue simulation in field is amplified in large scale and is simulated, in conjunction on cutting ring scale
Soil body osmotic coefficient is filled in the changing rule in space, anisotropy and heterogeneity, obtains the distribution of large scale infiltration coefficient.
It should be noted that the bicyclic infiltration test analogue simulation in field to be amplified to the simulation of large scale, top layer and bottom
It is respectively set to determine head boundary, head difference is equal to the vertical path length that water is flowed through from bottom to top layer.
Further, the determination method of a kind of high roadbed soil body large scale infiltration coefficient provided in an embodiment of the present invention, also
Include:
(1) bicyclic infiltration test simulation model is established, uniformly subdivision is set the top layer of simulation model using square volume mesh
It is set to and determines head, head value is 0m, will be also configured as determining head within the scope of simulation model bottom inner ring and outer rings, head value is equal to
Distance of the simulation model top layer to bottom.
(2) scale effect principle is combined, certain infiltration coefficient typical value is chosen, by isotropic medium, calculates different
Different random is distributed resulting permeability coefficient, and the infiltration coefficient measured value pair with the bicyclic infiltration test in field under subdivision scale
Than choosing suitable subdivision size.
(3) on determining subdivision size basis, analysis large scale equivalent coefficient of permeability calculated result is to different infiltrations
Number typical value quantity and anisotropic sensibility.
It should be noted that subdivision grid chooses certain infiltration coefficient Distribution value indoors on the basis of permeability test
Scheme carries out analogue simulation to bicyclic permeability test with condition simulation method in conjunction with scale effect, and with the bicyclic permeability test in field
As a result it is verified, determines the subdivision size of model, and the sensitivity of the different schemes to permeability coefficient selected by subdivision grid
Property is analyzed.
Based on the above method, the bicyclic infiltration test analogue simulation in field is amplified in large scale and is simulated,
In conjunction with soil body osmotic coefficient is filled on cutting ring scale in the changing rule in space, anisotropy and heterogeneity, large scale is obtained
The distribution of infiltration coefficient；It specifically includes:
(1) it on determining subdivision size basis, chooses certain typical value quantity, model will be simulated by isotropic medium
It encloses and gradually amplifies, calculate the infiltration coefficient under random distribution different under each simulation context, by permeability coefficient with simulation model
Simulation context value when the increase enclosed tends towards stability is determined as the simulation context of large scale infiltration coefficient.
(2) analysis large scale equivalent coefficient of permeability calculated result determines the sensibility of the quantity of infiltration coefficient typical value
Suitable typical value number.
(3) analysis large scale equivalent coefficient of permeability calculated result determines tested out cutting ring to anisotropic sensibility
Large scale permeability coefficient under scale infiltration coefficient large sample.
(4) according to sample and overall relationship, the distribution of large scale infiltration coefficient is determined.
It should be noted that subdivision grid chooses certain infiltration coefficient Distribution value on determining subdivision size basis
Simulation context is gradually expanded in scheme, analogue simulation is carried out to the infiltration of large scale landfill soil, until equivalent infiltration in simulation context
Coefficient tends towards stability, and determines simulation context corresponding to large scale infiltration coefficient；On the basis of determining subdivision scale, to subdivision
The sensibility of the different schemes of permeability coefficient selected by grid is analyzed, and combines cutting ring scale infiltration coefficient in space
Changing rule, anisotropy and heterogeneity, the final distribution for determining large scale infiltration coefficient.
Embodiment:
By taking a cheuch of Canal in Loess Area as an example, the loess for pruning top is filled out in ditch, becomes level land.One layer of process of banketing
One layer of progress, every layer of depth of fill are about 1m, and the depth of fill is about 90m, produce the high roadbed soil body therefrom.Specifically
Steps are as follows:
(1) the high roadbed soil body is sampled with work progress, the embankment soil body is sampled, the distribution of sample point is shown in
Fig. 2.Each sample point samples 2-3 soil sample, and soil sample is taken back interior, and it is real to carry out varying head infiltration with T-55 type permeameter
It tests, 133 groups of permeability coefficients is obtained, wherein 66 groups of vertical hydraulic conductivity, 67 groups of horizontal-hole blasting.
(2) for statistical analysis to the infiltration coefficient for filling native cutting ring scale.
In order to analyze the anisotropy and heterogeneity of landfill soil, hang down to what two samples of the same sample position measured
The infiltration coefficient of straight infiltration coefficient and horizontal-hole blasting and two same directions is compared, and sees Fig. 3 and Fig. 4 respectively,
It obtains, the heterogeneity for filling soil is very big, has very strong randomness, the infiltration tested on individual points by varying head
Coefficient is difficult to extrapolate the infiltration coefficient of large scale.
In order to further analyze infiltration coefficient with the changing rule of construction speed and landfill height, different samplings are drawn respectively
The box traction substation of vertical permeability coefficient and horizontal-hole blasting, is shown in Fig. 5~8 in number, different sampling altitude ranges.It obtains, vertically
Infiltration coefficient and horizontal-hole blasting are with sampling number and fill the variation of height without apparent changing rule, therefore can be with
United analysis is carried out as a whole.
The vertical permeability coefficient and horizontal-hole blasting of all landfill earth samples are statisticallyd analyze respectively, are shown in Table 1
With Fig. 9~10.As can be seen that the distribution of vertical permeability coefficient and horizontal-hole blasting is almost the same from table and figure.In order into
The analysis of one step, the cumulative frequency figure of infiltration coefficient is depicted with logarithmic axis, sees Figure 11.As can be seen from the figure vertical permeability coefficient
Frequency distribution with horizontal-hole blasting essentially coincides.
1 infiltration coefficient statistical form of table (unit: m/d)
| Statistical items | Vertical permeability coefficient | Horizontal-hole blasting |
| It is average | 0.0317 | 0.0322 |
| Standard error | 0.0042 | 0.0045 |
| Median | 0.02 | 0.0243 |
| Standard deviation | 0.0338 | 0.0368 |
| Variance | 0.0011 | 0.0014 |
| Kurtosis | 0.70 | 3.73 |
| The degree of bias | 1.25 | 1.84 |
| Region | 0.129 | 0.1716 |
| Minimum value | 0.00031 | 0.0003 |
| Maximum value | 0.1293 | 0.1719 |
| Observe number | 66 | 67 |
| Confidence level (95.0%) | 0.0083 | 0.009 |
(3) the bicyclic infiltration test in field has also been carried out in landfill area.Seepage ring high 25cm, outer annular diameter 50cm, inner ring are straight
Diameter 25cm.Obtain the infiltration coefficient result that the landfill soil area's infiltration test of 5 groups of field measures be respectively 0.030,0.030,0.041,
0.048、0.051m/d。
(4) the bicyclic infiltration test analogue simulation in field
Intermediate point is chosen in every 10% frequency range, as typical value, with the method for condition simulation, using Visual
Modflow carries out analogue simulation to the bicyclic infiltration test in field.In order to not influence boundary is arrived, planar range is selected as 3~4.5m, hangs down
It is selected as 0.4~0.45m to range, model is split into square.Each square has one permeability coefficient.Infiltration system
Number uses independent equiprobability random distribution in the distribution of each square, and the infiltration coefficient of a square is in different realizations
It can be different in (random distribution).
The distance between two soil samples due to the same sample position are within 15cm, the test result of infiltration coefficient
Still it differs greatly, so the side length of the square of model facetization is selected as 7.5,7.0,6.5 and 6.0cm, tetra- kinds of sizes and is carried out respectively
Simulation calculates.Meanwhile the value of infiltration coefficient is also calculated according to corresponding size according to scale effect principle, calculated result
It is shown in Table 2.
2 difference of table measures size permeability coefficient value (unit: m/d)
The infiltration coefficient of corresponding size is distinguished into assignment to model in the way of the distribution of isotropism equiprobability random independent
In each square, as shown in figure 12.It sets the top layer of model to determine head, head value is 0m；Model bottom inner ring and outer
It is also configured as determining head within the scope of ring, head value is equal to model top layer to the distance of bottom, as shown in figure 13.The initial water of model
Head is set as 0m.
It calculates by simulation it can be concluded that determining the water that head boundary enters model from bottom inner ring in the unit time.With wild
The Computing Principle of exosmosis coefficient is identical, and the ratio of the water and inner ring area is the equivalent infiltration that soil layer is simulated within the scope of inner ring
Saturating coefficient, calculated result are shown in Figure 14.By comparing with the bicyclic infiltration test in field, the subdivision side length for obtaining square is 7.0cm
When, result of the analog result closest to pit permeability test.
(5) sensitivity analysis is simulated in bicyclic infiltration test
In order to analyze infiltration coefficient to anisotropic sensibility, on the basis of the model of above-mentioned 7.0cm subdivision size, equally
Above-mentioned 10 typical values are chosen, and keep vertical permeability coefficient constant, it is randomly ordered by being carried out to horizontal-hole blasting, it obtains
To different combinations.Under the identical random distribution of comparative analysis, different vertical hydraulic conductivities and horizontal-hole blasting combination are calculated
Permeability coefficient, see Figure 15.It obtains, the simulation bicyclic infiltration test in field can count it according to isotropic medium
It calculates.
For the model base for analyzing infiltration coefficient to the sensibility of typical value number, equally in above-mentioned 7.0cm subdivision size
On plinth, 5,10 and 20 typical values are chosen respectively with identical method, and on the basis of identical random distribution, compare
The calculated result for analyzing infiltration coefficient, is shown in Figure 16.When obtaining simulating bicyclic infiltration test, selected typical value quantity is to result
It influences little.
(6) calculating of large scale infiltration coefficient
It, most can be outside reaction field when to obtain square size of mesh opening subdivision be 7cm by simulating the bicyclic infiltration test in field
Actual conditions.Above-mentioned 10 typical values are selected first, and mould is carried out to the square that simulation context is respectively 0.7,2.1 and 3.5m
Quasi-, i.e., plane difference subdivision is 100,900 and 2500 grids, and vertical subdivision respectively is 10,30 and 50 layers, top layer and bottom point
It is not set as determining head boundary, head difference is equal to the vertical path length that water is flowed through from bottom to top layer.For each subdivision
Small square presses isotropic medium, carries out equiprobability independent random tax to each square with 10 selected typical values
Value, and it is for statistical analysis to infiltration coefficient calculated under different probability distribution, see Figure 17.It obtains, simulation context 2.1m
The infiltration coefficient of large scale can be represented substantially with analog result when 3.5m.Later research is in the simulation context of 3.5m
It carries out.
In order to analyze influence of the typical value quantity to analog result, with same method, 5,10,20 are selected respectively
Typical value is 3.5m with simulation context, and subdivision is simulated having a size of the model of 7cm, and is compared and analyzed to result, sees figure
18.It obtains, the calculated result of infiltration coefficient has a certain distance under different number typical value, but gap is little, and relative difference exists
9% or so.Because the bigger closer actual conditions of typical value quantity, but the size of calculation amount is also considered simultaneously, so after
Analysis select 20 typical values.
In order to analyze influence of the anisotropy to large scale infiltration coefficient, on the basis of above-mentioned model, vertical infiltration system
Several and horizontal permeation equally chooses above-mentioned 20 typical values, and keeps vertical permeability coefficient constant, and horizontal-hole blasting is arranged at random
Column obtain different combinations, analyze the distribution of results situation of infiltration coefficient under various combination, the result is shown in Figure 19.It obtains, vertical infiltration
There are a certain distance for the result of resulting equivalent coefficient of permeability under saturating coefficient and horizontal-hole blasting various combination, but also differ
Less, relative difference is 10% or so.The average value of the equivalent coefficient of permeability obtained under various combinations and distribution is 0.0221m/
d。
The calculating of the above large scale infiltration coefficient is on the basis of certain pedotheque analyzes result, by statistics
Analysis chooses what the average value in different range was carried out as typical value.However, between certain amount sample and Large-scale Stratum
Belong to sample and overall relationship, the mean value of sample can not be fully equivalent to overall mean value.According to central-limit theorem, when
It (is generally required more than or equal to 30) when sample size is sufficiently large, mean valueSampling distribution approximation Normal Distribution, can be with
The confidence interval of population mean is obtained with following formula approximation:
In formulaFor the average value of sample, the critical probability or confidence level that α is confidence level when being (1- α), Zα/2It is significant
Critical value under horizontal α.According to vertical permeability coefficient and horizontal-hole blasting statistically analyze as a result, looking into standardized normal distribution
Table obtains the section of population mean under different confidence levels.The infiltration coefficient knot that varying head permeability test in large sample room is obtained
The mean value of fruit is approximate with field large scale infiltration coefficient to regard proportional relation as, it is possible thereby to calculate big ruler under different cumulative frequencies
The distribution for spending infiltration coefficient, is as a result shown in Figure 20.
Disclosed above is only several specific embodiments of the invention, and those skilled in the art can carry out the present invention
Various modification and variations without departing from the spirit and scope of the present invention, if these modifications and changes of the present invention belongs to the present invention
Within the scope of claim and its equivalent technologies, then the present invention is also intended to include these modifications and variations.
Claims (3)
1. a kind of determination method of high roadbed soil body large scale infiltration coefficient characterized by comprising
It is sampled in high roadbed soil body different level position and vertical position, and vertical to soil sample and horizontal with permeameter method
Be measured to infiltration coefficient, obtain cutting ring scale on landfill soil body osmotic coefficient the changing rule in space, anisotropy and
Heterogeneity；
The bicyclic infiltration test in field is carried out at high roadbed soil body different location, obtains the permeability coefficient on cutting ring scale；
According to the permeability coefficient on cutting ring scale, in conjunction with scale effect principle, emulation mould is carried out to the bicyclic infiltration test in field
It is quasi-；
The bicyclic infiltration test analogue simulation in field is amplified in large scale and is simulated, is seeped in conjunction with the soil body is filled on cutting ring scale
Saturating coefficient obtains the distribution of large scale infiltration coefficient in the changing rule in space, anisotropy and heterogeneity.
2. the determination method of high roadbed soil body large scale infiltration coefficient as described in claim 1, which is characterized in that further include:
Bicyclic infiltration test simulation model is established, uniformly subdivision is set as fixed using square volume mesh, by the top layer of simulation model
Head, head value are 0m, will be also configured as determining head within the scope of simulation model bottom inner ring and outer rings, head value is equal to emulation mould
Distance of the type top layer to bottom；
In conjunction with scale effect principle, chooses certain infiltration coefficient typical value, by isotropic medium, calculate different subdivision scales
Lower different random is distributed resulting permeability coefficient, and compares with the infiltration coefficient measured value of the bicyclic infiltration test in field, chooses
Suitable subdivision size；
On determining subdivision size basis, large scale equivalent coefficient of permeability calculated result is analyzed to different infiltration coefficient typical values
Quantity and anisotropic sensibility.
3. the determination method of high roadbed soil body large scale infiltration coefficient as claimed in claim 2, which is characterized in that it is described will be wild
Outer bicyclic infiltration test analogue simulation, which is amplified in large scale, is simulated, and is existed in conjunction with soil body osmotic coefficient is filled on cutting ring scale
Changing rule, anisotropy and the heterogeneity in space, obtain the distribution of large scale infiltration coefficient；Include:
On determining subdivision size basis, choose certain typical value quantity, by isotropic medium by simulation context gradually
Amplification, calculates the infiltration coefficient under random distribution different under each simulation context, by permeability coefficient with the increasing of simulation context
Simulation context value when tending towards stability greatly is determined as the simulation context of large scale infiltration coefficient；
Large scale equivalent coefficient of permeability calculated result is analyzed to the sensibility of the quantity of infiltration coefficient typical value, determines suitable generation
Tabular value number；
Large scale equivalent coefficient of permeability calculated result is analyzed to anisotropic sensibility, determines tested out cutting ring scale infiltration
Large scale permeability coefficient under coefficient large sample；
According to sample and overall relationship, the distribution of large scale infiltration coefficient is determined.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810613729.4A CN108982317B (en) | 2018-06-14 | 2018-06-14 | Method for determining large-scale permeability coefficient of high fill soil body |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810613729.4A CN108982317B (en) | 2018-06-14 | 2018-06-14 | Method for determining large-scale permeability coefficient of high fill soil body |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108982317A true CN108982317A (en) | 2018-12-11 |
| CN108982317B CN108982317B (en) | 2021-01-29 |
Family
ID=64540459
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810613729.4A Expired - Fee Related CN108982317B (en) | 2018-06-14 | 2018-06-14 | Method for determining large-scale permeability coefficient of high fill soil body |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108982317B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118392755A (en) * | 2024-06-25 | 2024-07-26 | 中国电建集团西北勘测设计研究院有限公司 | Method for obtaining equivalent permeability coefficient of heterogeneous mixed soil stratum |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105631078A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Numerical value simulation method for adaptive medium of natural fractured reservoir |
| CN106202746A (en) * | 2016-07-14 | 2016-12-07 | 南京大学 | The Yeh multi-level finite element modeling method of simulation Water in Porous Medium stream Darcy velocity |
| CN106884442A (en) * | 2017-04-25 | 2017-06-23 | 西安理工大学 | A kind of implementation of cheuch shape high fill foundation multi- scenarios method monitoring system |
-
2018
- 2018-06-14 CN CN201810613729.4A patent/CN108982317B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105631078A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Numerical value simulation method for adaptive medium of natural fractured reservoir |
| CN106202746A (en) * | 2016-07-14 | 2016-12-07 | 南京大学 | The Yeh multi-level finite element modeling method of simulation Water in Porous Medium stream Darcy velocity |
| CN106884442A (en) * | 2017-04-25 | 2017-06-23 | 西安理工大学 | A kind of implementation of cheuch shape high fill foundation multi- scenarios method monitoring system |
Non-Patent Citations (1)
| Title |
|---|
| 吴健华等: "高填方土体大尺度渗透系数的确定", 《中国地质学会2015学术年会论文摘要汇编（下册）》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118392755A (en) * | 2024-06-25 | 2024-07-26 | 中国电建集团西北勘测设计研究院有限公司 | Method for obtaining equivalent permeability coefficient of heterogeneous mixed soil stratum |
| CN118392755B (en) * | 2024-06-25 | 2024-10-11 | 中国电建集团西北勘测设计研究院有限公司 | Method for obtaining equivalent permeability coefficient of heterogeneous mixed soil stratum |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108982317B (en) | 2021-01-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Abbaspour et al. | Uncertainty in estimation of soil hydraulic parameters by inverse modeling: example lysimeter experiments | |
| Zhang et al. | Rainfall-runoff modelling in a catchment with a complex groundwater flow system: application of the Representative Elementary Watershed (REW) approach | |
| CN106770928B (en) | Ecology fragility mining area vegetation constrains the measuring method of lower WATER LEVEL CHANGES threshold | |
| Augusto Filho et al. | Landslide analysis of unsaturated soil slopes based on rainfall and matric suction data | |
| Caiqiong et al. | Application of HYDRUS-1D model to provide antecedent soil water contents for analysis of runoff and soil erosion from a slope on the Loess Plateau | |
| CN104331639B (en) | A kind of soil moisture content indirect gain and fast appraisement method | |
| KR101078297B1 (en) | System for infinite slope stability analysis considering saturation depth ratio of rainfall | |
| CN109061105B (en) | A kind of calculation method of the critical underground water buried depth of the soil salinization | |
| CN108876918A (en) | A kind of high roadbed loess body sedimentation and deformation prediction technique | |
| CN108961402A (en) | Space-time precision calibration method of the multi-satellite remote sensing precipitation inverting in large scale complexity basin | |
| CN108427741B (en) | DEM relative error evaluation method based on large number of high-precision control points | |
| CN108303512B (en) | Method for in-situ testing soil-water characteristic curve | |
| Malota et al. | Modelling mid-span water table depth and drainage discharge dynamics using DRAINMOD 6.1 in a sugarcane field in Pongola, South Africa | |
| Wang et al. | The methods and program implementation for river longitudinal profile analysis—RiverProAnalysis, a set of open-source functions based on the Matlab platform | |
| CN108982317A (en) | A kind of determination method of high roadbed soil body large scale infiltration coefficient | |
| CN102435540A (en) | Measurement system and method for permeability coefficient of shallow-layer sand soil | |
| CN117094167A (en) | Landslide whole-process surface three-dimensional deformation Beidou radar surface scene data simulation method | |
| Apaydin et al. | Determination of the drainage basin characteristics using vector GIS | |
| Tombul | Mapping field surface soil moisture for hydrological modeling | |
| Warrach-Sagi et al. | Streamflow data assimilation for soil moisture analysis | |
| Zhang et al. | The effect of texture and irrigation on the soil moisture vertical-temporal variability in an urban artificial landscape: A case study of Olympic Forest Park in Beijing | |
| Scott et al. | Illinois soil moisture under sod experiment | |
| Holberg | Downward model development of the soil moisture accounting loss method in HEC-HMS: Revelations concerning the soil profile | |
| Jabro et al. | Evaluation and comparison of five simulation models for estimating water drainage fluxes under corn | |
| CN111308045A (en) | Method for determining soil texture |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210129 |