CN107908849B - Landslide hazard degree quantification method - Google Patents
Landslide hazard degree quantification method Download PDFInfo
- Publication number
- CN107908849B CN107908849B CN201711091919.6A CN201711091919A CN107908849B CN 107908849 B CN107908849 B CN 107908849B CN 201711091919 A CN201711091919 A CN 201711091919A CN 107908849 B CN107908849 B CN 107908849B
- Authority
- CN
- China
- Prior art keywords
- landslide
- particles
- determining
- hazard
- rock
- 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
Images
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
The invention relates to a landslide hazard quantification method, and belongs to the technical field of data processing methods specially suitable for specific applications. The landslide hazard level under the corresponding threshold value at the corresponding moment is obtained through the following processes: (1) discretizing the landslide into a plurality of discrete particles; (2) determining mechanical parameters of landslide rock and soil mass materials; (3) determining iteration time step length, and obtaining the mass, the speed and the position of rock-soil body particles at different moments; (4) summarizing rock-soil body particles with displacement larger than a certain set threshold value, and determining a particle set with quantified landslide hazard; (5) determining the main direction of the landslide at the moment; (6) the method provided by the invention comprehensively considers a plurality of factors influencing the danger degree when the landslide occurs, provides a method for quantitatively evaluating the possible danger degree when the landslide occurs, and thus, the method can be used for pertinently treating and reinforcing.
Description
Technical Field
The invention relates to a landslide hazard quantification method, and belongs to the technical field of data processing methods specially suitable for specific applications.
Background
When landslide occurs, a large amount of rock and soil mass flow is often accompanied, the danger degree of the landslide is determined to a certain extent by the flowing speed and direction and the volume of the flowing soil mass, China is one of countries with frequent geological disasters such as landslide, in order to effectively avoid adverse effects caused by the disasters, stability analysis is carried out on the landslide in advance, the landslide with lower stability degree is treated and reinforced so as to avoid landslide disasters, the possible danger degree when the landslide occurs is also quantitatively evaluated, so that treatment and reinforcement are carried out in a targeted manner, and the landslide monitoring and reinforcing device is work with social significance and engineering value.
There are many factors that affect the degree of risk when a landslide occurs, but the speed and the amount of the landslide body are two relatively important factors. Conventional landslide stability analysis methods, such as the limit balance method and the finite element strength reduction method, cannot or cannot easily obtain the velocity of the rock-soil mass, and thus cannot quantify the degree of landslide hazard. Therefore, there is a need for a method that can quantify the risk of landslide.
Disclosure of Invention
The invention aims to overcome the defects of the conventional landslide hazard quantification method, and provides the landslide hazard quantification method.
The invention is realized by adopting the following technical scheme:
the method comprises the following steps: discretizing the landslide into a plurality of discrete particles having a discrete size Δ hx,Δhy,Δhx,ΔhyThe calculation time and the calculation precision are comprehensively considered for determining, and a trial calculation method from large to small is adopted until a certain discrete size has little influence on the result.
Step two: determining mechanical parameters of landslide rock and soil mass materials;
the soil layer rock mechanical parameters comprise mass m, position x, y and speed vx、vyThe method comprises the following steps of determining the slope of the landslide, the extension range of the slope top and the slope angle and the composition condition of a landslide soil layer by using field tests and measurement, and determining the cohesive force c, the internal friction angle phi, the elastic modulus E and the Poisson ratio u according to indoor conventional geotechnical tests.
Step three: determining iteration time step length, and determining the mass, speed and position of the particles at different moments;
the time step is determined by the following formula:
where C is the empirical coefficient and s is the speed of sound.
Step four: determining a set of particles for which the landslide hazard level is quantified;
the particle set is determined by the following method: will di>Particles of Δ d as a set of particles for landslide hazard quantification, where diAnd d is the displacement between all the particles and the initial value of the particles, and is a displacement threshold value, and the size is determined according to engineering experience.
Step five: determining the main direction of the landslide at the moment;
the main direction of the landslide is determined by the following method: selecting the uppermost and lowermost particles of the slider in the set of particles as muAnd mdWith corresponding coordinates xu,yu、xd,ydCalculating the unit vector of the main sliding direction of the landslide from the two points
Step six: and summarizing the projection of the momentum of the identified particle set in the main direction of the landslide, and summarizing to obtain the landslide risk degree under the corresponding threshold value at the corresponding moment.
Calculating the momentum of the sliding particle in the main sliding direction
A projected value R ofiSummarizing to obtain the landslide risk under the corresponding threshold value at the corresponding moment
The invention has the beneficial effects that:
the invention comprehensively considers a plurality of factors which influence the danger degree when the landslide occurs, including the quality, the speed and the square quantity of the landslide body, and provides a method for quantitatively evaluating the possible danger degree when the landslide occurs, thereby pertinently managing and reinforcing.
Drawings
FIG. 1 is a block flow diagram of the present invention.
Fig. 2 is a schematic discrete diagram of the landslide geometric model according to the present invention.
FIG. 3 is a geometric model diagram of homogeneous landslide.
Fig. 4 is a geometric model diagram of the sliding state when T is 0.6 seconds.
Fig. 5 is a geometric model diagram of the sliding state when T is 1.3 seconds.
Detailed Description
The invention will be further explained with reference to the drawings.
Fig. 1 shows a flow chart of the present invention, which is further described in detail below with reference to an example.
As shown in fig. 3, for a homogeneous landslide, the height of the slope is 10m, the width of the slope is 45 m, and the slope ratio is 1:2, the process of realizing landslide hazard quantification by the invention is as follows:
the method comprises the following steps: discretizing the landslide into a plurality of discrete particles having a discrete size Δ hx,Δhy；
Step two: determining mechanical parameters of landslide rock and soil mass materials;
as shown in FIG. 2, determining a reasonable discrete size by considering the calculation time and the calculation accuracy in combination is preferable to obtain Δ hx＝ΔhyThe landslide geometric model shown in fig. 3 is discretized into 11125 circular particles with different diameters, most of the particles have the diameter of 0.2m, a small amount of particles with the diameter of less than 0.2m exist at the slope, the mass of most of the particles is 80kg, and the mass of a small amount of the particles is less than 80 kg;
step three: determining iteration time step length, and determining the mass, speed and position of the particles at different moments;
setting delta t to be 0.0002 second to ensure iterative computation convergence, then performing iterative computation of a smooth particle fluid dynamics method, and storing the mass, speed and position information of all particles;
step four: determining a set of particles for which the landslide hazard level is quantified;
the most significant characteristic when a landslide occurs is displacement, so a displacement threshold value Δ d can be set to select a sliding particle set, the value of the threshold value is generally set according to engineering experience, and is usually 0.2m, that is, if a certain particle generates displacement of more than 0.2m, the particle is considered to have slid. For example, fig. 4 shows a selected sliding particle set in dark color when T is 0.6 seconds, the set including a total of 3089 sliding particles, and the sliding particle set includes a total of 5630 sliding particles as shown by the dark particles in fig. 5 when T is 1.3.
Step five: determining the main direction of the landslide at the moment;
when T is 0.6 seconds, the uppermost and lowermost particles of the slider are selected from the particle set and denoted by muAnd mdX of itu＝37.76,yu＝14.72；xd＝11.55,yd＝5.05；
When T is 1.3, x thereofu＝43.15,yu＝14.72，xd＝7.74,yd＝5.05，
Step six: and summarizing the projection of the momentum of the identified particle set in the main direction of the landslide, and summarizing to obtain the landslide risk degree under the corresponding threshold value at the corresponding moment.
When T is 0.6 seconds and 1.3 seconds, the momentum of the dark particles shown in fig. 4 and 5 was calculated, respectively
Then calculating the momentum in the main slip direction
A projected value R ofiAnd finally, quantitative evaluation of the degree of landslide hazard when T is 0.6 and 1.3 seconds
Wherein, r (k) and f (k) are sequence number mapping functions of dark particles in the sliding particle set in all particles respectively, and the mapping functions are automatically calculated by a program when the sliding particles are screened. For example, when T is 0.6 seconds, the following program code is used to calculate r (k) sequence number mapping function:
N＝0
Do i＝1,11125
if (displacement of i-th particle greater than Δ d) then
N＝N+1
R(N)＝i
End if
END DO
Of course, the foregoing is only a preferred embodiment of the invention and should not be taken as limiting the scope of the embodiments of the invention. The present invention is not limited to the above examples, and equivalent changes and modifications made by those skilled in the art within the spirit and scope of the present invention should be construed as being included in the scope of the present invention.
Claims (4)
1. A landslide hazard quantification method is characterized in that: the method comprises the following steps:
the method comprises the following steps: discretizing the landslide into a plurality of discrete particles having a discrete size Δ hx,Δhy；
Step two: determining mechanical parameters of landslide rock and soil mass materials; the soil layer rock mechanical parameters comprise mass m, position x, y and speed vx、vyCohesive force c, internal friction angle value phi, elastic modulus E and Poisson ratio u;
step three: determining iteration time step length, and determining the mass, speed and position of the particles at different moments; the time step is determined by the following formula:
wherein C is an empirical coefficient and s is the speed of sound;
step four: determining a set of particles for which the landslide hazard level is quantified; the particle set is determined by the following method: will di>Particles of Δ d as a set of particles for landslide hazard quantification, where diThe displacement between all the particles and the initial value is shown, delta d is a displacement threshold value, and the size is judged according to engineering experience;
step five: determining the main direction of the landslide at the moment; the main direction of landslide is determined by the following method:
selecting the uppermost and lowermost particles of the slider in the set of particles as muAnd mdWith corresponding coordinates xu,yu、xd,ydCalculating the unit vector of the main sliding direction of the landslide from the two points
Step six: and summarizing the projection of the momentum of the identified particle set in the main direction of the landslide, and summarizing to obtain the landslide risk degree under the corresponding threshold value at the corresponding moment.
2. The landslide hazard quantifying method according to claim 1, wherein: step one size Δ h of the discrete particlesx,ΔhyThe calculation time and the calculation precision are comprehensively considered for determination, and a trial calculation method from large to small is adopted until a certain discrete size has little influence on the result.
3. The landslide hazard quantifying method according to claim 1, wherein: and in the second step, the soil layer geotechnical parameters are determined by utilizing field tests and measurement to determine the slope of the landslide, the extension range of the top and the angle of the slope and the composition condition of the soil layer of the landslide and according to indoor conventional geotechnical tests.
4. The landslide hazard quantifying method according to claim 1, wherein: in the sixth step, the main direction of the landslide of the particle set is calculated
A projected value R ofiSummarizing to obtain the landslide risk under the corresponding threshold value at the corresponding moment
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711091919.6A CN107908849B (en) | 2017-11-08 | 2017-11-08 | Landslide hazard degree quantification method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711091919.6A CN107908849B (en) | 2017-11-08 | 2017-11-08 | Landslide hazard degree quantification method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107908849A CN107908849A (en) | 2018-04-13 |
| CN107908849B true CN107908849B (en) | 2020-12-11 |
Family
ID=61842710
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711091919.6A Active CN107908849B (en) | 2017-11-08 | 2017-11-08 | Landslide hazard degree quantification method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107908849B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109457739B (en) * | 2018-11-08 | 2019-06-18 | 青岛理工大学 | A kind of Slope safety level evaluation method based on the downstream structures extent of damage |
| CN111475937B (en) * | 2020-04-03 | 2020-12-11 | 中国地质科学院地质力学研究所 | Simulation method for flow-solid two-phase flow fluidization landslide |
| WO2021223198A1 (en) * | 2020-05-08 | 2021-11-11 | 青岛理工大学 | Method for quantizing landslide risk |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2489526A (en) * | 2011-04-01 | 2012-10-03 | Schlumberger Holdings | Representing and calculating with sparse matrixes in simulating incompressible fluid flows. |
| CN102819650B (en) * | 2012-08-16 | 2015-04-29 | 同济大学 | Computational simulation method of flow slide catastrophe of rock and soil material |
| CN105956317B (en) * | 2016-05-18 | 2017-03-15 | 青岛理工大学 | Risk of landslip quantization method |
-
2017
- 2017-11-08 CN CN201711091919.6A patent/CN107908849B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN107908849A (en) | 2018-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107908849B (en) | Landslide hazard degree quantification method | |
| Jin et al. | Bayesian model selection for sand with generalization ability evaluation | |
| CN109030202B (en) | Method for rapidly determining discrete element model parameters of rock brittle material | |
| Kim et al. | Miniature cone tip resistance on sand in a centrifuge | |
| Kirkpatrick et al. | Probabilistic model for the chloride-induced corrosion service life of bridge decks | |
| Chin et al. | A simplified method to estimate the soil-water characteristic curve | |
| Cornelis et al. | A simplified parametric model to describe the magnitude and geometry of soil shrinkage | |
| CN109655903B (en) | Shale layer transverse wave velocity prediction method and system | |
| CN109800459B (en) | Gravity type retaining wall design method and device | |
| Vanderborght et al. | Determining convective lognormal solute transport parameters from resident concentration data | |
| CN107013813A (en) | A kind of water supply line amount of leakage estimating system and method | |
| CN109187744A (en) | A kind of shear wave velocity evaluation method based on cone penetration test | |
| Nieto Leal et al. | Improved rotational hardening rule for cohesive soils and definition of inherent anisotropy | |
| CN111475876B (en) | Method for obtaining dynamic resilience mechanical characteristic parameters of granules | |
| CN106568647A (en) | Nerve network-based concrete strength predication method | |
| Bayat et al. | Estimation of the soil water retention curve using penetration resistance curve models | |
| Araei | Artificial neural networks for modeling drained monotonic behavior of rockfill materials | |
| CN109901238B (en) | High-stress formation resistivity correction method based on stress difference resistivity experiment | |
| CN109344534B (en) | Injection-production string critical erosion flow rate determination method and device | |
| Nesamatha et al. | Numerical modeling for prediction of compression index from soil index properties | |
| Liu et al. | Prediction of elastic compressibility of rock material with soft computing techniques | |
| Wang et al. | Experimental characterizations of an aging mechanism of sands | |
| Dastider et al. | Advancement in estimation of undrained shear strength through fall cone tests | |
| WO2021223198A1 (en) | Method for quantizing landslide risk | |
| Tran | Re-visitation of actual evaporation theories |
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 |