CN113239435B - Method for determining optimal water discharge speed of reservoir - Google Patents
Method for determining optimal water discharge speed of reservoir Download PDFInfo
- Publication number
- CN113239435B CN113239435B CN202110509063.XA CN202110509063A CN113239435B CN 113239435 B CN113239435 B CN 113239435B CN 202110509063 A CN202110509063 A CN 202110509063A CN 113239435 B CN113239435 B CN 113239435B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- dam
- discharge speed
- failure probability
- speed
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000007599 discharging Methods 0.000 claims abstract description 29
- 230000008595 infiltration Effects 0.000 claims abstract description 14
- 238000001764 infiltration Methods 0.000 claims abstract description 14
- 238000013178 mathematical model Methods 0.000 claims abstract description 4
- 239000002689 soil Substances 0.000 claims description 10
- 230000001174 ascending effect Effects 0.000 claims description 4
- 238000012614 Monte-Carlo sampling Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000009736 wetting Methods 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000012502 risk assessment Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 231100000817 safety factor Toxicity 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Control Of Non-Electrical Variables (AREA)
Abstract
A method for determining the optimal water discharge speed of a reservoir relates to the technical field of safety risk assessment and comprises the following steps: step 1, establishing a reservoir dam mathematical model, and designing different water discharging speeds under the condition of the same highest water level; step 2, obtaining the positions of the saturation lines corresponding to different time steps at a certain water discharge speed; step 3, acquiring failure samples of the reservoir dam at the infiltration line positions corresponding to different time steps, and calculating corresponding failure probability; step 4, repeating the step 3 to obtain the maximum values of a plurality of failure probabilities; and 5, changing the water discharging speed, obtaining a relation curve between the water discharging speed and the maximum failure probability of the reservoir dam, and obtaining the optimal water discharging speed of the reservoir according to the expected failure probability. The invention provides a method for determining the optimal water discharge speed of a reservoir, which can evaluate the optimal water discharge speed and is beneficial to the safety and stability of a dam.
Description
Technical Field
The invention relates to the technical field of safety risk assessment, in particular to a method for determining the optimal water discharge speed of a reservoir.
Background
In the actual slope engineering, the slope instability phenomenon caused by reservoir water level reduction is often encountered, so that the influence of the dam water level reduction speed on the dam stability is especially important.
The main reason why the water level drops to cause the instability of the dam is that when the water level drops at different speeds, the infiltration line in the dam does not drop along with the water level outside the side slope at the same time, the pore water pressure of the soil body of the dam cannot be dissipated, the static water pressure provided by the water level outside the side slope of the dam is dissipated after the water level drops, the vertical effective normal stress at the bottom of the soil body is reduced, the total stress is kept unchanged, and the instability of the dam can be easily caused.
Because the research on the influence of the reservoir water level falling speed on the dam stability is not comprehensive enough, the selection of the reservoir water discharge speed still faces challenges, and a method for determining the optimal water discharge speed of the reservoir is urgently needed.
Disclosure of Invention
In order to solve the problems, the invention provides a method for determining the optimal water discharging speed of a reservoir.
In order to achieve the purpose, the technical scheme of the invention is as follows:
a method for determining the optimal water discharge speed of a reservoir comprises the following steps:
and 5, changing the water discharging speed, obtaining a relation curve between the water discharging speed and the maximum failure probability of the reservoir dam, and finally obtaining the optimal water discharging speed of the reservoir according to the expected failure probability.
Preferably, in the step 1, p reservoir discharge speed values are given and recorded as V 1 ,V 2 ,V 3 ,……Vp。
Preferably, in step 2, the water discharge speed V is set i Next, N time steps A are generated by using the Seep/w module in Geo-studio 1 ,A 2 ......A N And the position of the wetting line thereof.
Preferably, in step 3, in time step A i Under the corresponding saturation line position, M soil layer parameter samples X are generated by Matlab software based on a random method 1 ,X 2 ,……,X M Calculating X based on the Sbings method in the extreme balance method 1 ,X 2 ,……,X M Corresponding factor of safety Fs 1 ,Fs 2 ,……,Fs M (ii) a If Fs i <1, then called X i Repeating the above steps to obtain m failed samples, and calculating A i Dam failure probability P corresponding to each time step fi =m/M。
Preferably, in the step 4, the dam failure probability P corresponding to each of the N time steps of the dam is calculated by repeating the step 3 f1 ,P f2 ,P f3 ,……,P fN (ii) a Corresponding failure probability P of N time steps f1 ,P f2 ,P f3 ,……,P fN Performing ascending arrangement and taking the maximum value to obtain the water discharge speed V i Corresponding maximum probability of failure P i 。
Preferably, in the step 5, the steps 2 to 4 are repeated to obtain the maximum failure probability P corresponding to each of P water discharge speed values i Using (V) i ,P i ) And i is 1, 2 and … p, drawing a relation curve between the water discharge speed and the maximum failure probability of the reservoir dam, and inquiring the curve to determine the optimal water discharge speed of the reservoir according to the expected failure probability.
The method for determining the optimal water discharging speed of the reservoir has the beneficial effects that: according to the invention, under different water discharge speeds of the reservoir dam, the infiltration line positions of different time steps are obtained, and the relation curve between the water discharge speed of the reservoir and the maximum failure probability of the reservoir dam is obtained by combining the uncertainty of the dam parameters, so that the optimal water discharge speed of the reservoir can be effectively obtained according to the expected failure probability, and the safety and stability of the dam are facilitated.
Drawings
FIG. 1, a flow chart of the present invention;
FIG. 2 is a schematic cross-sectional view of an eastern red reservoir in Ullin city, Shanxi province;
fig. 3, a schematic diagram of the positions of the saturation lines of the dam at different time steps of 1 m/d;
fig. 4 is a schematic diagram of the positions of dam infiltration lines at different time steps of 2 m/d;
fig. 5 is a schematic diagram of the positions of the dam infiltration lines at different time steps of V-4 m/d;
FIG. 6 is a graph of maximum failure probability of a reservoir dam at different water discharge speeds;
Detailed Description
In the following, embodiments of the present invention are described in detail in a stepwise manner, which is merely a preferred embodiment of the present invention and is not intended to limit the scope of the present invention, and any modifications, equivalents, improvements, etc. made within the spirit and principle of the present invention should be included in the scope of the present invention.
In the description of the present invention, it should be noted that the terms "upper", "lower", "left", "right", "top", "bottom", "inner", "outer", and the like indicate orientations and positional relationships based on the orientations and positional relationships shown in the drawings, and are only used for describing the present invention and simplifying the description, but do not indicate or imply that the referred device or element must have a specific orientation and a specific orientation configuration and operation, and thus, the present invention is not to be construed as being limited thereto.
A method for determining the optimal water discharge speed of a reservoir comprises the following steps:
as shown in figure 1, in the step 1, p reservoir discharge speed values are given and recorded as V 1 ,V 2 ,V 3 ,……Vp;
As shown in FIG. 1, in step 2, the water discharge speed V is set i Next, N time steps A are generated by using the Seep/w module in Geo-studio 1 ,A 2 ......A N And the position of the saturation line thereof;
as shown in fig. 1, in step 3, in time step a i Under the corresponding saturation line position, M soil layer parameter samples X are generated by Matlab software based on a random method 1 ,X 2 ,……,X M Calculating X based on the Sbingo method in the limit balance method 1 ,X 2 ,……,X M Corresponding factor of safety Fs 1 ,Fs 2 ,……,Fs M (ii) a If Fs i <1, then called X i Repeating the above steps to obtain m failed samples, and calculating A i Dam failure probability P corresponding to each time step fi =m/M;
As shown in fig. 1, in the step 4, the dam failure probability P corresponding to each of the N time steps of the dam is calculated by repeating the step 3 f1 ,P f2 ,P f3 ,……,P fN (ii) a Corresponding failure probability P of N time steps f1 ,P f2 ,P f3 ,……,P fN Performing ascending arrangement and taking the maximum value to obtain the water discharge speed V i Corresponding maximum probability of failure P i ;
As shown in fig. 1, in the step 5, the steps 2 to 4 are repeated to obtain the maximum failure probability P corresponding to each of P water discharge speed values i Using (V) i ,P i ) And i is 1, 2, … p, drawing a relation curve between the water discharging speed and the maximum failure probability of the reservoir dam, and inquiring the curve to determine the optimal water discharging speed of the reservoir according to the expected failure probability.
The following example is described with reference to fig. 2.
The dam crest of the dam is 100m in length, the width is 15m, the maximum dam height is 31.5m, the stratum lithology of the dam foundation is a mud-sandstone interbed, and the basic quality level of the strongly weathered rock of the dam foundation is level v; when the dam is constructed, the residual slope soil on the slope surface of the dam abutment is removed, the dam body of the dam is silty sand, and the permeability coefficient is 2.8 multiplied by 10 -5 cm/s, internal angle of friction of shear strengthThe average value is 20 degrees, the standard deviation is 8, the average value of the cohesive force c is 9kpa, the standard deviation is 3.6, and the volume weight gamma is 21kN/m 3 。
The steps of calculating and evaluating by adopting the method provided by the invention are as follows:
according to the figure 2, a dam numerical model is established by utilizing a Seep/w module in Geo-studio software, and p is set to be 3 reservoir water discharging speeds V 1 =1m/d,V 2 =2m/d,V 3 4 m/d; generating N-3 time step samples A by using Seep/w module 1 ,A 2 ,A 3 And corresponding characteristic wetting line position, at a certain time step A i Under the condition, the cohesive force c and the internal friction angle of the silty sand are taken into consideration due to the uncertainty of the parameters of the silty sandAnd (4) regarding the soil as a lognormal random variable, and keeping the volume weight of the soil unchanged. Method based on logarithm random variables1000 soil parameter samples X are generated by utilizing Matlab software 1 ,X 2 ,……,X 1000 Calculating 1000 soil body parameter samples corresponding to M-1000 safety factors based on a Sbingo method in a limit balance method, and recording as Fs 1 ,Fs 2 ,……,Fs 1000 。
In the case of M being 1000 safety factors, if the safety factor Fs i <1, then called X i Repeating the steps to obtain m failure samples in total for the failure samples, and obtaining the failure samples according to a failure probability formula P fi Calculating the failure probability P corresponding to 3 time step samples of the dam according to M/M f1 ,P f2 ,P f3 The failure probability P corresponding to each time step f1 ,P f2 ,P f3 Performing ascending arrangement and taking the maximum value to obtain the maximum failure probability P corresponding to the water discharge speeds of different water levels fmax1 =7.7%,P fmax2 =12.5%,P fmax3 (V) 25.6%, using i ,P fmaxi ) And (3) drawing a relation curve between the water level discharge speed and the maximum failure probability of the reservoir dam according to the i-1, 2 and 3, and inquiring the curve to determine the optimal water discharge speed value of the reservoir according to the expected failure probability of the dam. For example, if the expected failure probability of a dam is 0.1, the optimal discharge speed of the reservoir is 1.5 m/d.
By comparison with the prior art, it can be found that: the prior research on the water discharging speed of the reservoir mostly depends on experience to judge, the consideration on the dam stability is not comprehensive enough, and the invention fully considers the influence of each time-step saturation line on the dam stability under different water discharging speeds, calculates the failure probability of each time-step saturation line, takes the maximum value, and obtains the optimal water discharging speed according to the expected failure probability. The failure probability of the dam is gradually increased along with the increase of the water discharging speed through the engineering example verification, finally, the curve can be inquired according to the expected failure probability to determine the optimal water discharging speed, and the effectiveness of the invention is verified through example comparison and analysis.
Claims (5)
1. A method for determining the optimal water discharge speed of a reservoir is characterized by comprising the following steps: the method comprises the following steps:
step 1, establishing a reservoir dam mathematical model, and designing different reservoir water discharging speeds under the condition of ensuring the same highest water level;
step 2, acquiring infiltration line positions corresponding to different time steps of the reservoir dam at a certain water discharge speed;
step 3, under the infiltration line positions corresponding to different time steps, considering the uncertainty of the reservoir dam material, obtaining a failure sample of the reservoir dam by utilizing a Monte Carlo sampling method and combining a Sbingo method in a limit balance method, and further calculating the failure probability corresponding to the characteristic infiltration line position of the reservoir dam;
step 4, repeating the step 3 to obtain a plurality of failure probabilities, and taking the maximum value of the failure probabilities;
step 5, changing the water discharging speed, obtaining a relation curve between the water discharging speed and the maximum failure probability of the reservoir dam, and finally obtaining the optimal water discharging speed of the reservoir according to the expected failure probability;
in the step 3, in the time step A i Under the corresponding saturation line position, M soil layer parameter samples X are generated by Matlab software based on a random method 1 ,X 2 ,……,X M Calculating X based on the Sbings method in the extreme balance method 1 ,X 2 ,……,X M Corresponding factor of safety Fs 1 ,Fs 2 ,……,Fs M (ii) a If Fs i <1, then called X i Repeating the above steps to obtain m failed samples, and calculating A i Dam failure probability P corresponding to each time step fi =m/M。
2. The method for determining the optimal discharge speed of the reservoir as claimed in claim 1, wherein: in the step 1, p reservoir water discharge speed values are given and recorded as V 1 ,V 2 ,V 3 ,……Vp。
3. A method for determining an optimal discharge speed of a reservoir as claimed in claim 2, characterized in that: in the step 2, the water discharging speed V is i Next, N time steps A are generated by using the Seep/w module in Geo-studio 1 ,A 2 ......A N And its corresponding position of the wetting line.
4. A method for determining an optimal discharge speed of a reservoir as claimed in claim 3, characterized in that: in the step 4, the step 3 is repeated to calculate dam failure probability P corresponding to N time steps of the dam respectively f1 ,P f2 ,P f3 ,……,P fN (ii) a Corresponding failure probability P of N time steps f1 ,P f2 ,P f3 ,……,P fN Performing ascending arrangement and taking the maximum value to obtain the water discharge speed V i Corresponding maximum probability of failure P i 。
5. The method for determining the optimal discharge speed of the reservoir as claimed in claim 4, wherein: in the step 5, the steps 2 to 4 are repeated to obtain the maximum failure probability P respectively corresponding to the P water discharge speed values i Using (V) i ,P i ) And i is 1, 2, … p, drawing a relation curve between the water discharging speed and the maximum failure probability of the reservoir dam, and inquiring the curve to determine the optimal water discharging speed of the reservoir according to the expected failure probability.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110509063.XA CN113239435B (en) | 2021-05-11 | 2021-05-11 | Method for determining optimal water discharge speed of reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110509063.XA CN113239435B (en) | 2021-05-11 | 2021-05-11 | Method for determining optimal water discharge speed of reservoir |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113239435A CN113239435A (en) | 2021-08-10 |
CN113239435B true CN113239435B (en) | 2022-09-20 |
Family
ID=77133448
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110509063.XA Active CN113239435B (en) | 2021-05-11 | 2021-05-11 | Method for determining optimal water discharge speed of reservoir |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113239435B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113887074B (en) * | 2021-10-25 | 2024-04-12 | 青岛理工大学 | Method for determining optimal slope angle of upstream side slope of reservoir dam |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103245391A (en) * | 2013-05-07 | 2013-08-14 | 山东大学 | Model test apparatus and test method for monitoring dam saturation line |
CN203249668U (en) * | 2013-05-07 | 2013-10-23 | 山东大学 | Dam body saturation line monitoring model testing device |
CN105279361A (en) * | 2015-04-10 | 2016-01-27 | 中国电建集团成都勘测设计研究院有限公司 | Method for detecting instability risk ratio of slope of high no-overtopped rockfill cofferdam |
CN105913184A (en) * | 2016-04-11 | 2016-08-31 | 青岛理工大学 | Tailing dam instability risk evaluation method based on real-time monitoring data |
CN106021853A (en) * | 2016-05-09 | 2016-10-12 | 青岛理工大学 | Method for developing slope stability and reliability analysis software under spatial variation characteristic |
CN107142858A (en) * | 2017-06-01 | 2017-09-08 | 中国路桥工程有限责任公司 | Single slope drainage system inside a kind of seasonal immersion roadbed |
CN111985041A (en) * | 2020-09-17 | 2020-11-24 | 青岛理工大学 | Reinforced side slope retaining wall height determination method and reinforced side slope retaining wall |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109117540B (en) * | 2018-08-02 | 2022-11-18 | 三峡大学 | Probability statistical analysis method for solving non-uniqueness of inversion of dam concrete mechanical parameters |
CN112183894B (en) * | 2020-10-29 | 2021-07-06 | 四川大学 | Gravity dam earthquake risk probability analysis method |
-
2021
- 2021-05-11 CN CN202110509063.XA patent/CN113239435B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103245391A (en) * | 2013-05-07 | 2013-08-14 | 山东大学 | Model test apparatus and test method for monitoring dam saturation line |
CN203249668U (en) * | 2013-05-07 | 2013-10-23 | 山东大学 | Dam body saturation line monitoring model testing device |
CN105279361A (en) * | 2015-04-10 | 2016-01-27 | 中国电建集团成都勘测设计研究院有限公司 | Method for detecting instability risk ratio of slope of high no-overtopped rockfill cofferdam |
CN105913184A (en) * | 2016-04-11 | 2016-08-31 | 青岛理工大学 | Tailing dam instability risk evaluation method based on real-time monitoring data |
CN106021853A (en) * | 2016-05-09 | 2016-10-12 | 青岛理工大学 | Method for developing slope stability and reliability analysis software under spatial variation characteristic |
CN107142858A (en) * | 2017-06-01 | 2017-09-08 | 中国路桥工程有限责任公司 | Single slope drainage system inside a kind of seasonal immersion roadbed |
CN111985041A (en) * | 2020-09-17 | 2020-11-24 | 青岛理工大学 | Reinforced side slope retaining wall height determination method and reinforced side slope retaining wall |
Non-Patent Citations (1)
Title |
---|
降雨作用下排土场边坡稳定性分析;陈耀华;《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》;20170215;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113239435A (en) | 2021-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113239435B (en) | Method for determining optimal water discharge speed of reservoir | |
JP6712543B2 (en) | Corrosion amount estimation device and method | |
WO2021047121A1 (en) | Method for determining particle displacement critical value after slope instability | |
Ssenyondo et al. | Effects of embedment depth on the pullout capacity of bucket foundations in sand | |
CN110377935B (en) | Method for predicting change of excess pore water pressure along depth of soil body around pile in pile sinking process | |
Zedan et al. | Performance Assessment of Shirin Earth Dam in Iraq Under Various Operational Conditions | |
CN112765803A (en) | Method and system for evaluating failure probability of anti-uplift stability of foundation pit bottom | |
CN108427654B (en) | Rapid calculation method for silted storage capacity of medium-sized or over-sized check dam | |
Intui et al. | Behavior of bearing capacity on pile foundation during fluctuating groundwater level | |
Caballero et al. | Advances in the understanding of the hydraulic behavior of wedge-shape block spillways | |
CN113312697B (en) | Method for predicting anti-sliding stability of high-pressure grouting on-shore soil blocking wall | |
CN113700057B (en) | Method for measuring foundation settlement | |
JP6432785B2 (en) | How to prevent ground liquefaction | |
CN112541257B (en) | Method for calculating service life of limestone slope based on corrosion rate | |
CN113887074A (en) | Novel method for determining optimal slope angle of upstream side slope of reservoir dam | |
CN104021280B (en) | Method for computing critical hydraulic gradient suitable for piping of tail silt | |
CN113158314A (en) | Slope stability analysis method | |
CN115186433B (en) | Analysis method for dynamic hydraulic characteristic curve of stratum caused by slurry shield | |
CN113420453A (en) | Debris flow damming dam bursting flood peak flow calculation method | |
Upomo et al. | Behaviour of levee on softsoil caused by rapid drawdown | |
Tanaka et al. | Effects of relative density of sand on seepage failure phenomena in front of sheet piles in two dimensions | |
CN110569548A (en) | Method for determining infiltration line of anisotropic rock slope | |
CN114357718A (en) | Method for determining maximum weakening coefficient of bank slope rock-soil body strength parameter | |
CN115573399A (en) | Method for predicting bearing capacity of single-disc supporting disc pile | |
CN114813500B (en) | Method for measuring loess stratum matrix suction force by using on-site immersion test |
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 |