A kind of high sand-containing flood " is taken off river bed " and is washed away method of discrimination
Technical field
The present invention relates to a kind of hydrologic forecasting method, be specifically related to a kind of high moisture flood and produce the method for discrimination that " taking off river bed " washes away, belong to the sediment dynamics technical field, the IPC international Patent classificating number is E02B1/00.
Background technology
Phenomenon that high sand-containing flood " is taken off river bed " is one view on the Yellow River, generally occurs in Little Valleys In The Middle Reaches of The Yellow Xiao Bei master stream section and Downstream of The Weihe River section.When high sand-containing flood passes through this section, early stage is through the clay piece (channel aggradation thing) that density is large, intensity is higher of certain hour alluvial, started and surface by current, several square metres even tens square metres of areas, some clods erect vertical with water (flow) direction as wall, then " thump " pour (clod that has was just flutterred in water after upright two or three minutes) in water once into, by flood, swallowed up very soon.On river surface, a large amount of clods rise one after another, downstream stream seethe and under, bloom in full river, surging, the underwater sound is deafening." taking off river bed " and washing away is a kind of phenomenon that hyper-concentration flow and the borderline phase mutual effect of special riverbed produce." take off river bed " and wash away and often cause the strong incision in riverbed, sometimes a flood peak can rush riverbed dark several meters, causes riverbed and water level to decline to a great extent, and souring often causes migrating of river course major trough simultaneously, the littoral engineering in river course is done great damage, to the Yellow River flood control, bring immense pressure.
Last century the seventies, " taking off river bed " phenomenon has just caused showing great attention to of domestic and international water conservancy working person, and this problem has been carried out research." take off river bed " at present and wash away the research method of problem and divide three classes: the first kind is to utilize the field data at prototype hydrometric station, research and analyse the husky condition of water that " taking off river bed " phenomenon occurs for Xiao Bei master stream, the Yellow River and section, the Weihe River, but because the husky condition of high silt-laden river water is larger along Cheng Bianhua, therefore with " take off river bed " field data at upstream and downstream hydrometric station, section is inaccurate as the calculated data of the husky condition of " taking off river bed " flush water, and because the condition that " taking off river bed " phenomenon forms is more special, the randomness of its generation is larger, be difficult to capture in practice stronger " taking off river bed " field data of followability, also therefore caused the statistics length limited, therefore the husky condition of water of " taking off river bed " phenomenon generation of the method proposition does not have convincingness.The Equations of The Second Kind method is to utilize the husky condition of water in conceptual approach clay piece when starting of sediment transport capacity, but sediment carrying capacity reflection is transport conditions in silt downstream transport process, dynamic mechanical relation in the time of can not reflecting that the clay piece is instantaneous and start, and the clay piece is taken off the starting of lifting with single sand grain also very large difference, therefore, can reflect " taking off river bed " problem with the concept of sediment carrying capacity, be worth discussion.The 3rd class research method is experimental study, the method utilizes flume test simulation " taking off river bed " to wash away phenomenon, and built two kinds of clay pieces and take off and lift mechanical model: a kind of is that clay piece bottom is by the current erosion, when clod occurred upwards to rotate along the continuous section with riverbed (fixed end place), clay piece front end was in cantilever position; Another kind is not consider the water eroding and flooding situation, the adhesion strength of only considering clay block boundary and riverbed builds mechanics analysis model, but these two kinds of models are all to take off at clay piece full wafer the calculating of under the situation of lifting, the husky condition of water being carried out, and all do not consider the impact of the structural strength of clay piece (rupture strength) own.
" take off river bed " at present and wash away the research of problem, the differences such as emphasis due to research means, method, object and research, therefore the criterion of " taking off river bed ", discriminant criterion, mechanism etc. are still lacked unified understanding, and because the condition that " taking off river bed " phenomenon forms is more special, the randomness of its generation is larger, be difficult to capture in practice stronger " taking off river bed " field data of followability, also caused the research of " taking off river bed " phenomenon at present to be restricted.
can channel aggradation thing block be started, can produce in other words " taking off river bed " phenomenon, depend primarily on the synthesis result of various conditions, comprising channel aggradation in early stage form and adjustment situation, deposit density and relative roughness, deposit block formational situation and lift height, deposit Block Boundary condition, the mechanical strength of alluvial piece itself, peak flood flow and sand content size, peb process duration length, the channel geometry parameter, the layer vertical pressure fluctuation meter is strong and weak, pressure fluctuation meter is at the PHASE DISTRIBUTION of deposit block upper and lower surface and Phase Stacking probability etc.
Wash away the present achievement in research of problem from " taking off river bed ", deriving " taking off river bed " while washing away critical criterion, consider the Research Literature that above-mentioned " taking off river bed " wash away influence factor and not yet find.The present invention utilizes " taking off river bed " to wash away clay piece mechanical strength experimental data, and consideration " taking off river bed " is washed away the combined influence factor, the critical discrimination formula that derive " taking off river bed " washes away.
Summary of the invention
The present invention washes away under combined influence factor condition in consideration " taking off river bed ", and structure clay piece fractures to take off and lifts mechanics analysis model, and utilizes clay piece rupture strength experimental data, calculates the critical flow condition that the clay piece fractures and causes " taking off river bed " phenomenon to occur.The technical scheme that the present invention takes is:
A kind of high sand-containing flood " is taken off river bed " and is washed away method of discrimination, comprises the following steps:
(1), " take off river bed " and wash away an easily section clay piece rupture strength test:
The test soil sample washes away for taking from " taking off river bed " the clay piece of easily sending out section, at first the test soil sample is reinvented, concrete grammar is: will wash away and easily send out the clay piece that the former horseshoe bend in section fetches and pulverize from " taking off river bed ", put into container and add the water stirring, make silt flocculation deposition, by air-dry formation clay piece sample, and utilize the different moisture content of controlling clay piece sample of air-dry time; Then,, on the same group time clay piece sample not, utilize Intensity of Breaking Tester to measure rupture strength;
Test method is as follows:
Clay piece sample is disposed across on two support cylinder of Intensity of Breaking Tester, the test specimen major axis is perpendicular to the Intensity of Breaking Tester support cylinder, with the speed of the 50N/s ± 10N/s loading cylinder by Intensity of Breaking Tester, vertical load is added on the horizontal plane of prism equably, until clay piece sample fractures; Utilize formula (1) to calculate the rupture strength M of clay piece sample,
In formula, F
fVertically be added in the load at its top planes middle part when-test specimen fractures, unit is N;
Distance between L-support cylinder, unit is mm;
The length of side of B-prism square sectional, unit is mm;
Use the average of one group three clay piece sample rupture strengths as result of the test, utilize experimental result to set up the relation equation of clay piece rupture strength and moisture content, calculate " taking off river bed " by this relation equation and wash away rupture strength M under the clay piece saturated aqueous rate of easily sending out section
s
(2), build " taking off river bed " while washing away the clay piece fracture to take off and lift stressed Mathematical Modeling:
F
DFor the positive thrust of current:
In formula, V
bFor clay piece bottom water flow flow velocity, C
DFor resistance coefficient, γ
mFor muddy water unit weight, b is clay piece width, and c is clay piece thickness.
For F on plane of rupture
DReaction force:
τ is the current drag power of clay piece upper surface:
G
sClay piece weight is: G
s=γ
sLbcg(5)
In formula, γ
sFor the unit weight of clay piece, l is the erosion jib-length, and g is acceleration of gravity, and other symbol is the same.
P
1For the suffered water pressure of clay piece: P
1=γ
mLbhg(6)
In formula, h is the clay piece surface depth of water, and other symbol is the same.
P
2Buoyancy for the clay piece: P
2=γ
mLb (h+c) g(7)
Maximum flow of water pulsation uplift force:
In formula, K is linear coefficient, and units/kg ms, span: 3~4.2, J are water-surface gradient.
(3), consider to set up the torque equilibrium equation of each power to the y axle under least favorable Load Combination condition:
(4), equation solution, draw critical flow velocity V
b:
At first, to the formula (2) in step (two) and (three)-(9) simultaneous solution, obtain formula (10):
The clay piece is adopted empirical formula by the erosion jib-length of scour within the Δ t time period:
In formula, C
1For bed scour coefficient, τ
sFor the suffered drag force of sand grain around the clay piece, τ
cCritical incipient motion shearing stress when the channel aggradation thing washes away, γ
cUnit weight for the channel aggradation thing.
In formula, ρ
mFor the density of high sand-containing flood, D is the average grain diameter of movable deposit around the clay piece, and A is the combined influence coefficient, and its value is 0.044, and coefficient ξ value is 1.35, d
mFor the average grain diameter of erosion deposit, d
iFor washed out sand grain particle diameter,
Simultaneous formula (10), (11), (13) and (14) solve, and can obtain V
b, this V
bProduce exactly the critical flow velocity that " taking off river bed " washes away;
(5), judge whether certain flood produces " taking off river bed " and wash away phenomenon.
When " taking off river bed " washes away easy upstream, section generation high sand-containing flood, according to the observation data of this upstream, section streamflow measuring station to flood flow rate of water flow and sand content, actual riverbed flow velocity, sand content when calculating this flood and arriving this section and wash away duration, then utilize the method for step (four) to calculate the critical flow velocity V that generation " is taken off river bed " and washed away
b, when this actual riverbed flow velocity greater than this critical flow velocity V
bThe time, " taking off river bed " can occur wash away phenomenon, otherwise, " taking off river bed " can not occur wash away phenomenon.
Method of discrimination of the present invention, whether the high moisture flood that can forecast with unerring accuracy can produce " taking off river bed " is washed away phenomenon, for flood control provide decision-making foundation, simultaneously, also for " take off river bed " real-time data acquisition that washes away provides may.
The specific embodiment
Committed step of the present invention is that " taking off river bed " section clay piece is carried out the rupture strength test; Consider " take off river bed " wash away the combined influence factor and build the clay piece and fracture to take off and lift model mathematics, and model is solved.Its concrete implementation step is as follows:
(1) " take off river bed " and wash away an easily section clay piece rupture strength test.
The test soil sample is that river bed " is taken off " and wash away the clay piece of easily sending out section in the Yellow River.At first wash away according to " taking off river bed " the formation mechanism of easily sending out section clay piece flocculation deposition, reinvent wash away the easy soil sample of fetching section of sending out from " taking off river bed ", concrete grammar is: will pulverize from the clay piece that former horseshoe bend is fetched, put into container and add the water stirring, make silt flocculation deposition,, as sample, and utilize the different moisture content of controlling clay piece sample of air-dry time by air-dry formation clay piece., on the same group time clay piece sample not, utilize Intensity of Breaking Tester to measure rupture strength.
Test method is as follows:
Clay piece sample is disposed across on two support cylinder of Intensity of Breaking Tester, the test specimen major axis is perpendicular to the Intensity of Breaking Tester support cylinder, with the speed of the about 50N/s ± 10N/s loading cylinder by Intensity of Breaking Tester, vertical load is added on the horizontal plane of prism equably, until clay piece sample fractures.The design formulas of rupture strength M is suc as formula (1), and unit is newton/square millimeter (MPa).
In formula, F
fVertically be added in the load at its top planes middle part when-test specimen fractures, unit is N;
Distance between L-support cylinder, unit is mm;
The length of side of B-prism square sectional, unit is mm.
Use the average of one group three clay piece sample rupture strengths as result of the test.Utilize experimental result to set up the relation equation of clay piece rupture strength and moisture content, calculate " taking off river bed " by this relation equation and wash away rupture strength M under the clay piece saturated aqueous rate of easily sending out section
s
(2) build " taking off river bed " while washing away the clay piece fracture to take off and lift stressed Mathematical Modeling, calculate clay piece load.
F
DThe positive thrust of current:
In formula, V
bFor clay piece bottom water flow flow velocity, C
DFor resistance coefficient, value is 1, γ
mFor muddy water unit weight, b is clay piece width, and c is clay piece thickness.
For F
DReaction force:
τ ' is the fracture vertical shearing of section part of clay piece, and it does not produce moment to the y axle, can ignore.
τ is the current drag power of clay piece upper surface:
τ " be the current drag power of clay piece soffit, it does not produce moment to the y axle, can ignore.
Clay piece weight G
s: G
s=γ
sLbcg(5)
In formula, γ
sFor the unit weight of clay piece, l is the scour jib-length, and g is acceleration of gravity, and other symbol is the same.
The suffered water pressure P of clay piece
1: P
1=γ
mLbhg(6)
In formula, h is the clay piece surface depth of water, and other symbol is the same.
The buoyancy P of clay piece
2: P
2=γ
mLb (h+c) g(7)
Maximum flow of water pulsation uplift force:
In formula, K is linear coefficient, and units/kg ms, span: 3~4.2, J are water-surface gradient.
(3) consider to set up the torque equilibrium equation of each power to the y axle under least favorable Load Combination condition
(4) equation solution, draw flow velocity V
b, the critical flow velocity that " taking off river bed " washes away occurs exactly:
In formula (4), (5), (6), (7), (8) substitution formula (9), arrange:
The clay piece is adopted empirical formula by the erosion jib-length of scour within the Δ t time period:
In formula, C
1For bed scour coefficient, τ
sFor the suffered drag force of sand grain around the clay piece, τ
cCritical incipient motion shearing stress when the channel aggradation thing washes away, γ
cUnit weight for the channel aggradation thing.
In formula,
For the drag force coefficient, value is 0.4, α
DFor the area coefficient of sand grain, value is
ρ
mFor the density of water, d
iFor washed out sand grain particle diameter, other symbol is the same.Replace d with the average grain diameter D of movable deposit around the clay piece
i, arrangement formula (12):
Channel aggradation thing critical incipient motion shearing stress τ
c:
In formula, through test of many times, determine, the accurate value of combined influence coefficient A is 0.044, and the span of coefficient ξ is 1.25~1.65, and its accurate value is 1.35, d
mFor the average grain diameter of erosion deposit, d
iFor washed out sand grain particle diameter.
In formula (13), (14) substitution formula (11), can obtain the design formulas of scour jib-length l:
Around consideration clay piece, but the particle diameter of erosion sand grain is 10
-5(m) magnitude left and right, the span of A is 10
-2Magnitude, and
Therefore, in formula (15)
Replace d with the average grain diameter D of movable deposit around the clay piece
i, formula (15) can be expressed as:
Will
L is as variable, simultaneous formula (10), formula (16) solve: with formula (10) substitution formula (16) and arrange:
Order
Formula (17) is converted into:
In formula,
Obtain according to the radical formula of simple cubic equation:
In l substitution formula (10), can obtain
(5) judge whether certain flood produces " taking off river bed " and wash away phenomenon.
When " taking off river bed " washes away easy upstream, section generation high sand-containing flood, according to this upstream, section streamflow measuring station to flood flow rate of water flow and sand content observation data, according to this section hydrological observation data for many years, actual riverbed flow velocity, sand content in the time of can calculating at an easy rate this flood and arrive this section and wash away duration, then utilize step (four) method to calculate the critical flow velocity V that generation " is taken off river bed " and washed away
bWhen this actual riverbed flow velocity greater than this critical flow velocity V
bThe time, " taking off river bed " can occur wash away phenomenon, otherwise, " taking off river bed " can not occur wash away phenomenon.