Additional mass method numeral template is asked for the method for rockfill density
Technical field
The present invention relates to additional mass method numeral template and ask for the method for rockfill density, specifically on rockfill, set up set of data samples and set up plane coordinate system according to the additional mass method test of one group of basis measuring point and the correlative value of hole survey method test, the set of data samples of setting up is carried out triangle gridding and contour tracing obtains digital template, and utilize the digital template of setting up to calculate the method for rockfill density in the difference position of same rockfill, belong to the exploration geophysics technical field.
Background technology
The compacted dry density of rockfill dam project is an important indicator of Engineering Quality Control.Equal and less than the density measurement method of the thick enrockment soil of 300mm for maximum particle diameter, Ministry of Water Resources's " earthwork test rule " (SL237-1999), Ministry of Geology and Mineral Resources's " earthwork test rule " (DT-82) proposed douche and fill with micromicrofarad, and the size of testing pits (diameter and the degree of depth) is defined as 3~5 times of medium maximum particle diameter." rolling earth-rock dam job specifications " be middle regulation (DL/T5129-2001): the stockpile density test, should adopt douche, the diameter of testing pits is not less than 2~3 times of maximum particle diameter, and maximum is no more than 2m, the degree of depth of testing pits is rolling layers thickness, also can be aided with Surface Wave Method, sedimentation.Obviously the size of testing pits of (DL/T5129-2001) regulation and stipulation wants little many than what " earthwork test rule " stipulated.Collateral security measuring accuracy angle, application are that large particle diameter is tested pits greatly; From the engineering construction angle, test pits unsuitable excessive; Testing requirement and construction requirement have certain contradiction.Therefore, many engineering technical personnel attempt to seek a kind of accurately and fast method and address this problem for a long time.It is reported, since about 20 century 70s, just proposed more than 10 and planted enrockment soil Density Detection method, facts have proved, for super-size enrockment soil, method more likely or dynamic parameter method.Namely measure the method for density by the kinetic parameter (elastic modulus, foundation stiffness, elastic wave velocity, Poisson ratio) of measuring enrockment soil, because theoretical research proof kinetic parameter and Media density have close relationship.The additional mass rule is to measure the simple and direct effective ways of foundation stiffness and mass of vibration.
Chinese patent ZL96117338.6 discloses a kind of additional mass method of rockfill density measurement, it belongs to the method for asking in early days rockfill density is the calibration Y-factor method Y, the method need to be in experimental stage by a large amount of measuring points, additional mass method is carried out calibration with hole survey method set up the calibration matrix of coefficients, need all will carry out calibration to various different lithology dam material, different gradation, different construction parameter.In the stage of filling, then ask for the density value of measuring point in conjunction with inverse in the additional mass method test parameter substitution density calculation formula by the calibration matrix of coefficients of setting up.The shortcoming of the method is: 1. need a large amount of contrast measuring points could set up complete calibration matrix of coefficients; 2. because the complicacy of dam material and the diversity of construction parameter, the method of this calibration tends to cause between the calibration coefficient to produce conflict, namely the calibration coefficient of certain dam material and construction parameter be applied on the contrary calculate in the test result of another kind of dam material and construction parameter just suitable or all improper; 3. the inverse density value of asking for additional mass method test can only pass through point-to-point corresponding relation.
Summary of the invention
Utilize the calibration Y-factor method Y to calculate the problem that rockfill density exists for additional mass method, the object of the present invention is to provide a kind of additional mass method to utilize digital template to calculate the method for rockfill density, for achieving the above object, technical solution of the present invention is:
A kind of additional mass method utilization numeral template is asked for the method for rockfill density, on rockfill, set up plane coordinate system according to the additional mass method test of one group of basis measuring point and the correlative value of hole survey method test, the set of data samples of setting up is carried out triangle gridding and contour tracing obtains digital template, and utilize the digital template of setting up to calculate rockfill density in the difference position of same rockfill, its concrete steps are as follows:
One sets up additional mass method numeral template
1. on rockfill, choose arbitrarily one group of basis measuring point, and on each measuring point, carry out respectively:
A) ginseng that is comprised of additional mass and rockfill at the measuring point system of shaking utilizes the additional mass method test to obtain joining to shake the mass of vibration of body
Dynamic stiffness of foundation K
0
B) adopt survey method in hole to dig pit, take a sample, weigh, measure the wet density that volume obtains cheating survey
Dry density
And water percentage η
0
C) calculate the vibratory volume of joining the body that shakes
Join the self-vibration angular frequency of the body that shakes
2. with the shake vibratory volume of body of the ginseng of each measuring point
With water percentage η
0For parameter forms set of data samples (K
0, ω
0 -2,
η
0), again with sample set (K
0, ω
0 -2,
η
0) middle ω
0 -2Multiply by experience enlargement factor n and obtain sample set (K
0, n * ω
0 -2,
η
0);
3. revise sample set (K
0, n * ω
0 -2,
η
0) precision, obtain new sample set (K
1, ω
1 -2,
η
1);
4. with dynamic stiffness of foundation K
1For ordinate, with the shake self-vibration angular frequency of body of ginseng
1 -2Be horizontal ordinate, set up plane coordinate system, with sample set (K
1, ω
1 -2,
η
1) be drawn in the coordinate system;
5. the sample set in the coordinate system is carried out triangle gridding and contour tracing, obtain the vibratory volume of basic measuring point
With water percentage η
1Isogram;
6. vibratory volume isogram and water percentage isogram are carried out the gridding covering, realize additional mass method numeral template;
Two utilize additional mass method numeral template to ask for rockfill density, and calculation procedure is as follows:
1. draw dynamic stiffness of foundation K in the difference position of same rockfill by the data analysis of additional mass method collection in worksite
2, the shake mass of vibration of body of ginseng
At first, calculate the self-vibration angular frequency of joining the body that shakes
2 -2, utilize K
2And n * ω
2 -2Find vibratory volume on the corresponding template at digital template
With water percentage η
2Then, according to formula
The wet density that obtains asking on the template
At last, according to formula ρ
2 do=ρ
2 is wet-ρ
2 is wet* η
2The dry density ρ that obtains asking on the template
2 do
Describedly choose arbitrarily one group of basis measuring point number and be no less than 10.
Described experience enlargement factor n=40.
Described mesh spacing is take 0.1 as unit.
Owing to having adopted above technical scheme, the present invention relates to additional mass method numeral template and ask for the method for rockfill density and overcome the calibration Y-factor method Y and can only pass through point-to-point corresponding relation, realized aspectant corresponding relation.The method is the same some contrast experiment who surveys method and additional mass method according to certain hole, set up the relation of limited several different lithology dam material, different gradation and additional mass method test parameter, obtain the additional mass method test parameter regularity of distribution, set up the digital template of several dams material additional mass method test parameter few in number, then progressively expand and revise this numeral template, directly ask for again the density value of measuring point by this template.The advantage of the method is: 1. do not need a large amount of contrast measuring points, only need limited several different lithology dam material, different gradation opening relationships, digital template can constantly be revised; 2. the data of setting up model are open-ended, realize the mobilism of data, and correction result is constantly convergence; 3. the method achievement has portability or applicability, achievement can be generalized in the Other Engineering to use.
Positive effect of the present invention is:
(1) invented the new method that additional mass method is asked for rockfill density, the method has realized aspectant corresponding relation, improved the precision of bulk density;
(2) realized that surveying contrast points with a small amount of hole can calculate all rockfill density;
(3) have portability, the result of calculation of a certain dam material can be used for the density calculation of other dam material, reduce the quantity of hole survey method, save the engineering funds.
Description of drawings
Accompanying drawing 1 sample set precision correction figure
The triangulation diagram that accompanying drawing 2 sample sets are set up
Accompanying drawing 3 volumetric water content isograms
Accompanying drawing 4 digital templates
Accompanying drawing 5 triangle interpolations calculate schematic diagram
Accompanying drawing 6 digital templates 1
Accompanying drawing 7 digital templates 2
Linear regression graph between accompanying drawing 8 digital templates
Embodiment
The method of additional mass method utilization numeral template of the present invention being asked for rockfill density below in conjunction with drawings and Examples further describes.
A kind of additional mass method utilization numeral template is asked for the method for rockfill density, the method is to set up data set and set up plane coordinate system according to the additional mass method test of one group of basis measuring point and the correlative value of hole survey method test on rockfill, the set of data samples of setting up is carried out triangle gridding and contour tracing obtains digital template, and utilize the digital template of setting up to calculate rockfill density in the difference position of same rockfill, its concrete steps are as follows:
One sets up additional mass method numeral template
1. on rockfill, choose arbitrarily one group of basis measuring point, and on each measuring point, carry out respectively additional mass method test and the test of hole survey method;
2. utilize additional mass method to record the mass of vibration of joining the body that shakes
Dynamic stiffness of foundation K
0
3. take to cheat the survey method in this measuring point point position, dig pit, take a sample, weigh, measure the hole hygrometric density that obtains this some position of volume
Dry density
And water percentage η
0
4. calculate the vibratory volume of joining the body that shakes
Join the self-vibration angular frequency of the body that shakes
5. with the shake vibratory volume of body of the ginseng of each measuring point
With water percentage η
0For parameter forms set of data samples (K
0, ω
0 -2,
η
0), again with sample set (K
0, ω
0 -2,
η
0) middle ω
0 -2Multiply by experience enlargement factor n and obtain sample set (K
0, n * ω
0 -2,
η
0);
6. revise sample set (K
0, n * ω
0 -2,
η
0) precision, Euclidean distance between sample is pressed less than the sample of threshold value H
Revise, obtain new sample set (K
1, ω
1 -2,
η
1);
7. with dynamic stiffness of foundation K
1For ordinate, with the shake self-vibration angular frequency of body of ginseng
1 -2Be horizontal ordinate, set up plane coordinate system, with sample set (K
1, ω
1 -2,
η
1) be drawn in the coordinate system.
8. set up sample set (K
1, ω
1 -2,
η
1) the triangulation network;
9. to sample set (K
1, ω
1 -2,
η
1) respectively to the shake vibratory volume of body of ginseng
With water percentage η
1Carry out contour tracing, the formation base measuring point is joined the vibratory volume of the body that shakes
η with water percentage
1Isogram;
10. volume isogram and water percentage isogram are carried out the gridding covering, realize additional mass method numeral template.
Two utilize additional mass method numeral template to ask for rockfill density, and calculation procedure is as follows:
1. draw dynamic stiffness of foundation K in the difference position of same rockfill by the data analysis of additional mass method collection in worksite
2, the shake mass of vibration of body of ginseng
2. calculate the self-vibration angular frequency of joining the body that shakes
3. in the digitizing template, find out a position p (K
2, n * ω
2 -2) be arranged in certain triangle T of the triangulation network;
4. choose a certain summit V in the triangle T wantonly
E0Cross some p do straight line and with limit V
E1V
E2Give mutually a p
t(K
t, ω
t -2), utilize formula (1) (2) calculation level p
tVibratory volume
And water percentage η
t
Wherein, K
E1, ω
-2 E1,
Be triangle summit V
E1Dynamic stiffness of foundation, self-vibration angular frequency, vibratory volume;
K
E2, ω
-2 E2,
Be triangle summit V
E2Dynamic stiffness of foundation, self-vibration angular frequency, vibratory volume.The coordinate of recycling formula (3) (4) calculation level p (
η
p), thereby obtain vibratory volume
With water percentage η
2
Wherein, K
E1, ω
-2 E1,
Be triangle summit V
E1Dynamic stiffness of foundation, self-vibration angular frequency, vibratory volume;
K
E2, ω
-2 E2,
Be triangle summit V
E2Dynamic stiffness of foundation, self-vibration angular frequency, vibratory volume.
5. calculate wet density
Dry density
Describedly choose arbitrarily one group of basis measuring point number and be no less than 10.
Described experience enlargement factor n=40.
Described mesh spacing is take 0.1 as unit.
The numeral template is asked for density method and other methods of asking for density of additional mass method have a unique distinction, and digital template has portability, and the density calculation of the rockfill of available and other characteristics reduces the quantity of hole survey method, saving engineering funds.The portability of numeral template may be embodied in: exist linear approximate relationship between the digital template of different qualities rockfill, its calculation procedure is as follows:
1. respectively characteristic 1 and the rockfill of characteristic 2 are set up digital template;
2. the rockfill in characteristic 1 utilizes the additional mass method test to obtain mass of vibration M
0With dynamic stiffness K;
3. in the digital template of characteristic 1 rockfill, use mass of vibration M
0Calculate ground rockfill point position dry density ρ with dynamic stiffness K
1In the digital template of the rockfill of characteristic 2, use mass of vibration M
0Calculate ground rockfill point position dry density ρ with dynamic stiffness K
2
4. other different some positions utilize the additional mass method test to obtain mass of vibration on the rockfill of characteristic 1 and characteristic 2
With dynamic stiffness K
1, obtain sample set (
K
1);
5. with sample set (
K
1) in sample utilize the digital template of characteristic 1 rockfill to calculate ground rockfill point position dry density set (ρ
1), and utilize the digital template of characteristic 2 rockfills to calculate ground rockfill point position dry density set (ρ
2);
6. set up dry density set (ρ by linear regression analysis
1) gather (ρ with dry density
2) linear approximate relationship, obtain equation of linear regression ρ
1≈ a ρ
2+ b;
7. realize being no less than linear approximate relationship between the digital template of two specific character stockpiles by same mode.
The present invention be multiply by enlargement factor n when setting up set of data samples, be because:
1. ω in sample
-2Differ very little, when setting up the triangulation network, will cause like this triangulation network to be set up unusual, thereby it is unusual to cause isoline to generate.Therefore, to ω
-2Multiply by an enlargement factor n, can be correct set up the triangulation network and contour tracing.In addition, search and dynamic stiffness K and self-vibration angular frequency at digital template
-2Corresponding vibratory volume V
0During with water percentage η, the self-vibration angular frequency
-2Also multiply by an enlargement factor n, just can in digital template, find smoothly corresponding vibratory volume V
0With water percentage η, result of calculation is not produced any Accuracy.
2. in addition, since ω
-2Can not be too little, therefore can not be too large, otherwise it is unusual also can to cause the triangulation network to be set up, thus it is unusual to cause isoline to generate.So enlargement factor n can not get too large, through experiment showed, that n gets 40 and can accomplish correct the attractive in appearance of the triangulation network, contour tracing and isogram of setting up.
Why to carry out the precision correction of sample, be because:
Consider in the situation that in the two-dimensional coordinate system Euclidean distance close, sample K, M
0Corresponding volume V
0Differing with water percentage η should be very little, and V
0, the trend consistent data of the variation of η and isoline graded is just reasonable.And usually dig pit and additional mass method can exist inevitable error and accidental error, cause V
0Thereby the inconsistent situation of trend of the variation of η and isoline graded causes V
0, the value of η differs larger.For this type of situation, need to carry out accuracy correction, the present invention has taked the method for method of mean error in the sample close quarters to carry out the accuracy correction of sample.
The performing step of method of mean error is as follows in the sample close quarters:
1. set up take K as ordinate ω
-2Two-dimensional coordinate system for horizontal ordinate.2. calculate the ultimate range between all samples, and two sample number of recording distance maximum.2. respectively according to the sample around predetermined two samples of zone radius search.3. with all samples in the zone radius
Average calculating, average error, simultaneously, with coordinate
Also average calculating, set up new sample (K
1, ω
1 -2, V
1, η
1) and replace all interior samples of sample areas.4. the sample labeling that mark is new was for calculating.5. finding out does not have calculated sample, and continues to calculate until all samples have all calculated according to above step.
Specific embodiment
Below in conjunction with embodiment the present invention is further described in detail:
Embodiment one
Select a measuring point to carry out the additional mass method collecting work in test site, draw rockfill mass of vibration M by acquisition software
0=510 (kg) and dynamic stiffness of foundation K=127.6 (MN/m).Dig pit at measuring point, weigh up the rockfill mass M of digging with platform scale
Heap, waterproof plastics are deployed troops on garrison duty in the hole, pour water, and calculate the volume of pouring water, calculate rockfill wet density ρ according to formula ρ=M/V
Wet=2.260 * 10
3Kg.m
-3, the rockfill that digs out is dried, the quality after obtaining drying, the difference of the quality after the front quality of oven dry and the oven dry is not mass M moisture in the rockfill
Water, according to formula η=M
Water/ M
Heap=1.5 (%) calculate water percentage, thereby draw rockfill dry density ρ
Do=ρ
Wet* (1-η)=2.23 * 10
3Kg.m
-3Get vibratory volume V by formula ρ=M/V
0=M
0/ ρ
Wet=225.66372 (m
3).
Similarly, material carries out the contrast work of other 24 measuring points on the dam of identical characteristics, and each secondary data of surveying sees Table 1.
The glutinous Zha Du of table 1 power station rockfill density hole survey method and additional mass method contrast table
Set up (K, ω
-2, V
0, sample set η) sees Table 2.
Table 2 is set up set of data samples
K |
ω
-2 |
V
0 |
η |
127.6 |
4.00 |
225.664 |
1.50 |
82 |
4.26 |
156.152 |
2.10 |
123.3 |
4.26 |
260.317 |
1.60 |
88.6 |
3.60 |
146.937 |
1.10 |
141.6 |
4.80 |
306.823 |
1.40 |
86.1 |
4.12 |
325.08 |
1.10 |
127.3 |
4.33 |
245.216 |
0.90 |
106.3 |
4.10 |
196.662 |
1.80 |
80.4 |
7.10 |
253.778 |
3.20 |
77.2 |
5.63 |
195.594 |
1.70 |
126.6 |
4.31 |
243.75 |
2.10 |
93.3 |
5.33 |
234.323 |
1.40 |
59.8 |
5.33 |
151.306 |
2.20 |
96.1 |
5.66 |
251.27 |
1.20 |
110.9 |
5.36 |
274.873 |
1.70 |
118.9 |
5.00 |
280.189 |
0.90 |
80.5 |
5.64 |
211.163 |
3.20 |
100 |
5.37 |
256.202 |
1.40 |
83.2 |
3.80 |
148.986 |
1.10 |
80.6 |
5.30 |
202.274 |
1.10 |
82.6 |
3.79 |
150.481 |
1.30 |
97.2 |
5.35 |
244.131 |
1.10 |
100.1 |
5.22 |
243.838 |
1.60 |
86 |
6.51 |
262.418 |
2.00 |
101.2 |
6.86 |
332.057 |
2.30 |
Take dynamic stiffness K as ordinate, with the shake self-vibration angular frequency of body of ginseng
-2Set up plane coordinate system, to (K, the ω that sets up
-2, V
0, sample set η) carries out the precision correction according to the algorithm shown in the accompanying drawing 1, and utilizes the revised sample set of precision to set up the triangulation network, as shown in Figure 2.
Utilize the triangulation network to carry out contour tracing, generate vibratory volume isogram and water percentage isogram, such as accompanying drawing 3.After isogram carried out gridding, obtain the digital template shown in the accompanying drawing 4.
After other some positions of test site utilize the additional mass method collection in worksite, also analyze and obtain mass of vibration
With dynamic stiffness K
0=100 (MN/m).Utilize formula
Calculate the self-vibration angular frequency
0 -2=5.1 (kg.m/MN).Search a p (K by digital template
0, ω
0 -2* 40) triangle at place, according in the accompanying drawing 5 with triangular apex V
E0Cross some p and limit V
E1V
E2Intersect at a p
t, utilize formula (1) (2) interpolation calculation to go out a p
tVibratory volume V
PtWith water percentage η
t, recycling formula (3) (4) interpolation calculation goes out the vibratory volume of a p
Water percentage η
p=1.62%.Calculate wet density
Dry density
The expansion of embodiment two digital templates
Utilize additional mass method to record the mass of vibration M of ground rockfill
0=247kg and dynamic stiffness of foundation K=69 (MN/m).The method of testing pits in additional mass method institute measuring point position, the method for taking to dig pit, take a sample, weigh, measure volume obtains the wet density ρ of ground rockfill point position
Wet=1.98 * 10
3Kg.m
-3, dry density ρ
Do=1.94 * 10
3Kg.m
-3And water percentage η=1.8%.Calculate the volume V that joins the body that shakes
0=M
0/ ρ
Wet=247/1.98=124.75m
3Utilize formula to calculate the self-vibration angular frequency
-2=M
0/ K sets up data sample K, ω
-2, V
0, η, and the new data sample is added in the sample set of digital template.Revise the sample set precision, rebulid the triangulation network, carry out contour tracing and gridding, form new digital template.
In the rock master windrow rockfill Density Detection of cattle pen river diversion works as example, table 3 has represented that the digital template that utilizes glutinous bundle to cross the Dam rockfill calculates the Density Detection result of cattle pen river surface sheetpile masonry dam, showed simultaneously the compare error of surveying method with cattle pen Jiang Keng, because the rockfill of cattle pen river surface sheetpile masonry dam numeral template is not set up, what adopt is the quantity numeral template that glutinous bundle crosses the Dam rockfill, so bigger error.By progressively the corresponding data sample of setting up additional mass method and hole survey method on cattle pen river surface sheetpile masonry dam, live the digital template of just setting up cattle pen river surface sheetpile masonry dam after first experiment contrast finishes.Simultaneously, along with the 2nd, 3 ... contrast, progressively expand and revise digital template, so that digital template precision progressively improves, the final the result of calculation of additional mass method during to the 5th, 6 and hole survey method error calculated very little (such as table 4), just there is larger difference in indivedual some position.
The 1st, 2 additional mass method test parameters of table 3 cattle pen river engineering master windrow and hole survey method comparing result
The 5th, 6 additional mass method test parameters of table 4 cattle pen river engineering master windrow and hole survey method comparing result
The portability of embodiment three method achievements
Use respectively the identical test parameter dynamic stiffness of foundation K in glutinous Zha Du power station, mass of vibration M
0In the template 1 (such as Fig. 6) of different characteristic dam material, template 2 (such as Fig. 7), calculate, wherein sample (K, M
0) being in the template computer capacity, the data list that calculates is as follows:
Table 5 is set up the template that expands
K | M
0 | Template one bulk density | Template two bulk densities |
104 | 580 | 2.15 | 2.21 |
140 | 573 | 2.24 | 2.27 |
127 | 588 | 2.19 | 2.20 |
90 | 509 | 2.19 | 2.22 |
105 | 580 | 2.14 | 2.22 |
113 | 563 | 2.07 | 2.11 |
110 | 594 | 2.13 | 2.22 |
87 | 462 | 2.15 | 2.21 |
111 | 611 | 2.14 | 2.22 |
123 | 581 | 2.14 | 2.16 |
94 | 470 | 2.18 | 2.12 |
112 | 511 | 2.21 | 2.23 |
105 | 477 | 2.21 | 2.15 |
145 | 681 | 2.17 | 2.23 |
105 | 402 | 2.16 | 2.20 |
116 | 648 | 2.17 | 2.23 |
86 | 378 | 2.04 | 2.21 |
77 | 448 | 2.23 | 2.23 |
86 | 557 | 2.18 | 2.16 |
132 | 581 | 2.19 | 2.23 |
90 | 440 | 2.15 | 2.18 |
106 | 557 | 2.14 | 2.14 |
95 | 514 | 2.18 | 2.20 |
102 | 424 | 2.17 | 2.20 |
111 | 610 | 2.16 | 2.22 |
86 | 387 | 2.10 | 2.20 |
101 | 437 | 2.15 | 2.15 |
96 | 418 | 2.13 | 2.13 |
109 | 533 | 2.12 | 2.05 |
99 | 534 | 2.11 | 2.16 |
81 | 493 | 2.12 | 2.20 |
85 | 420 | 2.10 | 2.20 |
80 | 473 | 2.13 | 2.21 |
78 | 450 | 2.13 | 2.18 |
112 | 693 | 2.15 | 2.23 |
Obtain following table 6 behind Data induction and the rejecting abnormalities point.
The tables of data of table 6 arrangement
Template one bulk density |
Template two bulk densities |
2.07 |
2.12 |
2.11 |
2.16 |
2.12 |
2.17 |
2.14 |
2.17 |
2.15 |
2.19 |
2.16 |
2.21 |
2.17 |
2.21 |
2.19 |
2.23 |
2.21 |
2.24 |
2.24 |
2.27 |
As shown in Figure 8, utilize Least Square Regression Analysis to draw, equation of linear regression is: y=0.8725x+0.3158.Be that the density value that template 1 calculates is x, the as a result y that utilizes regression equation calculation to go out is similar to the density value that template 2 calculates, thereby has realized the relation between the template.