CN115906256A - Reservoir landslide surge numerical simulation method and system - Google Patents
Reservoir landslide surge numerical simulation method and system Download PDFInfo
- Publication number
- CN115906256A CN115906256A CN202211547962.XA CN202211547962A CN115906256A CN 115906256 A CN115906256 A CN 115906256A CN 202211547962 A CN202211547962 A CN 202211547962A CN 115906256 A CN115906256 A CN 115906256A
- Authority
- CN
- China
- Prior art keywords
- landslide
- reservoir
- simulation
- fluid
- slope
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 83
- 238000004088 simulation Methods 0.000 title claims abstract description 72
- 239000002245 particle Substances 0.000 claims abstract description 96
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 54
- 239000012530 fluid Substances 0.000 claims abstract description 49
- 230000008569 process Effects 0.000 claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 33
- 239000011435 rock Substances 0.000 claims abstract description 25
- 239000002689 soil Substances 0.000 claims abstract description 21
- 238000010168 coupling process Methods 0.000 claims abstract description 20
- 230000008878 coupling Effects 0.000 claims abstract description 10
- 238000005859 coupling reaction Methods 0.000 claims abstract description 10
- 230000000704 physical effect Effects 0.000 claims abstract description 7
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims description 9
- 230000003993 interaction Effects 0.000 claims description 8
- 230000000007 visual effect Effects 0.000 claims description 6
- 238000005096 rolling process Methods 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000001808 coupling effect Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 230000007547 defect Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000008187 granular material Substances 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/23—Dune restoration or creation; Cliff stabilisation
Landscapes
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
The invention discloses a reservoir landslide surge numerical simulation method and a reservoir landslide surge numerical simulation system, which relate to the technical field of slope simulation analysis, and have the technical scheme that: acquiring topographic data before and after occurrence of reservoir landslide and surge disasters, and constructing a geometric model; dividing a reservoir flow field area by utilizing OpenFOAM software, and configuring fluid physical mechanical parameters; establishing a bank slope model on the geometric model by using LIGGGHTS software, and configuring physical and mechanical parameters and particle grading of particles according to physical properties of the bank slope rock-soil mass; performing analog simulation on the bank slope model by a CFDEM coupling method to obtain a simulation calculation result; and analyzing the simulation calculation result to obtain the slope movement process, the surge height and/or the influence range. According to the invention, through the coupling of the CFD method and the DEM method, the coupling effect among landslide bodies, reservoir water and air is fully considered, and two reservoir landslide and surge modes of granular break and dangerous rock collapse can be simulated.
Description
Technical Field
The invention relates to the technical field of slope simulation analysis, in particular to a reservoir landslide surge numerical simulation method and system.
Background
In recent years, the development of the water conservancy and hydropower industry in China is rapid, more and more hydraulic engineering is built, and the hidden danger of some geological disasters is gradually highlighted. Due to the influence of factors such as strong rainfall in reservoir areas, reservoir operation and the like, severe slope stability and landslide surge problems exist in a plurality of reservoir areas. When the reservoir is landslide and swells, huge swell is possibly lifted, and the threat to the life and property safety of people, the safety of ships, the operation of hydropower stations and the like is formed.
The numerical simulation method is an important means for developing research on the landslide and swell of the related reservoir at home and abroad, and can quickly, accurately and comprehensively analyze the dynamic process of the landslide and swell of the reservoir. However, the process of surging caused by unstable water entering of the side slope is complex, along with the fluid-solid coupling effect among rock-soil bodies, water bodies and air, and the space scale is large, so that the solution is difficult to be carried out by using a single medium model; in addition, rock-soil bodies and fluids have different mechanical properties, and different numerical simulation methods are required; in addition, different instability mechanisms exist on the bank side slopes under different geological environment conditions, such as granular landslide and dangerous rock body collapse, the accumulation forms and the surge propagation characteristics of the bank side slopes have certain differences, the traditional numerical simulation method is difficult to simulate large-scale slope deformation and simulate two reservoir landslide and surge modes of the granular landslide and the dangerous rock body collapse, and the application condition is limited to a certain extent.
Therefore, how to research and design a reservoir landslide surge numerical simulation method and system capable of overcoming the defects is a problem which is urgently needed to be solved at present.
Disclosure of Invention
In order to solve the defects in the prior art, the invention aims to provide a reservoir landslide surge numerical simulation method and a reservoir landslide surge numerical simulation system, which can simulate two reservoir landslide surge modes of granular break and dangerous rock collapse by coupling a CFD method and a DEM method and fully considering the coupling effect among landslides, reservoir water and air, obtain the calculation results of the motion process, the surge height, the influence range and the like after the reservoir landslide surge, and provide technical support for reservoir bank slope risk evaluation, disaster prevention and reduction.
The technical purpose of the invention is realized by the following technical scheme:
in a first aspect, a reservoir landslide swell numerical simulation method is provided, which comprises the following steps:
acquiring topographic data before and after occurrence of reservoir landslide and surge disasters, and constructing a geometric model;
dividing a reservoir flow field area in a geometric model by utilizing OpenFOAM software, and configuring fluid physical mechanical parameters;
establishing a bank slope model on the geometric model by using LIGGGHTS software, and configuring physical and mechanical parameters and particle grading of particles according to physical properties of the bank slope rock-soil mass;
performing analog simulation on the bank slope model by a CFDEM coupling method to obtain a simulation calculation result;
and analyzing the derived simulation calculation result by using visual post-processing software to obtain the slope movement process, the surge height and/or the influence range.
Further, the topographic data acquisition process before and after the occurrence of the reservoir landslide and surge disaster specifically comprises the following steps:
collecting terrain elevation data by adopting an unmanned aerial vehicle;
converting the terrain elevation data into terrain raster data;
and subtracting the terrain raster data before and after the disaster occurs to obtain the slope instability range.
Further, the construction process of the geometric model specifically includes:
according to terrain raster data before and after a disaster occurs and a slope instability range, a pre-processing script pre.sh is operated at a terminal, and a three-dimensional stl geometric model of an integral area and a instability area after the landslide and the swell of a reservoir is generated through a format conversion program;
and for generalized regular terrain, directly generating stl geometric models by adopting three-dimensional drawing software, and then importing and calculating.
Further, the dividing process of the reservoir flow field area specifically comprises the following steps:
determining the water level elevation, wherein the area below the water level elevation in the geometric model is a water body, and the area above the water level elevation is air;
the motion states of water and air in the landslide and swell process of the reservoir are described by solving a Navier-Stokes equation;
solving an interface between two-phase fluids by adopting a VOF model;
and, accounting for the effect of porosity in the fluid control equation to account for the volume occupied by the landslide particles.
Further, the expression of the fluid control equation is specifically as follows:
wherein alpha is p Represents porosity; alpha is alpha i Representing the volume fraction of air or water; u. of f Representing the fluid velocity; u. of c Represents the compression speed; rho f Represents the fluid density; r pf Representing the calculation of momentum exchange items based on a drag force model in each grid unit; τ represents the stress tensor of the fluid; f σ Represents the surface tension; when water and air are simultaneously present in the grid cells, the average characteristics of the two-phase fluid are used to characterize the cells.
Further, the fluid physical mechanical parameters comprise air density, air kinematic viscosity, water body density and water body kinematic viscosity.
Further, the configuration process of the physical and mechanical parameters and the particle grading of the particles specifically comprises the following steps:
obtaining the density, young modulus, poisson ratio, friction coefficient, coefficient of restitution and rolling friction coefficient of the rock-soil mass of the bank side slope through a geotechnical test or trial calculation mode, and analyzing to obtain the particle grading of the rock-soil mass of the bank side slope;
according to the particle grading of the rock-soil mass of the bank side slope, setting particles with different particle sizes according to a preset similarity ratio for simulation;
and calling different operation commands to respectively simulate granular landslide and dangerous rock collapse.
Further, the simulation of the granular body landslide and the dangerous rock collapse is specifically as follows:
for the granular landslide, setting loose particles with different particle sizes for simulation;
for dangerous rock body collapse, a plurality of small particles are subjected to rigid bonding to form a large-size block body without overlapping quantity for simulation.
Further, the process of performing simulation on the bank slope model by the CFDEM coupling method specifically includes:
the DEM solver calculates the forces acting on the particles and solves the particle momentum equation to update the position and velocity of each particle;
after DEM circulation is completed, calculating the porosity and fluid-particle interaction force in the CFD grid unit, and transmitting the interaction force to the CFD circulation;
solving a fluid momentum equation and a continuity control equation to obtain a fluid velocity field and a fluid pressure field;
after the flow field is converged, calculating the fluid acting force acting on the particles and returning the fluid acting force to the DEM solver;
and circulating the CFD-DEM calculation until the calculation end time is reached.
In a second aspect, a reservoir landslide surge numerical simulation system is provided, comprising:
the geometric modeling module is used for acquiring topographic data before and after a reservoir landslide surge disaster occurs and constructing a geometric model;
the region division module is used for dividing a reservoir flow field region in the geometric model by utilizing OpenFOAM software and configuring fluid physical mechanical parameters;
the side slope modeling module is used for establishing a bank side slope model on the geometric model by utilizing LIGGGHTS software and configuring particle physical mechanical parameters and particle gradation according to the physical properties of the bank side slope rock and soil mass;
the coupling analysis module is used for carrying out analog simulation on the bank slope model by a CFDEM coupling method to obtain a simulation calculation result;
and the simulation processing module is used for analyzing the derived simulation calculation result by utilizing the visual post-processing software to obtain the slope movement process, the surge height and/or the influence range.
Compared with the prior art, the invention has the following beneficial effects:
1. according to the reservoir landslide surge numerical simulation method provided by the invention, the coupling effect among landslides, reservoir water and air is fully considered through the coupling of the CFD method and the DEM method, two reservoir landslide surge modes of granular break and dangerous rock collapse can be simulated, the calculation results of the motion process, the surge height, the influence range and the like after the reservoir landslides and surges are obtained, and the technical support is provided for risk evaluation of reservoir bank slopes and disaster prevention and reduction.
2. The method can overcome the defects of the traditional simulation method in the research of the fluid-solid coupling phenomenon, and reduce the real evolution process of the reservoir landslide under the influence of the environment; the control equation of the water body and the air is a Navier-Stokes equation, contains momentum exchange terms acting with the particles, and simultaneously considers the influence of porosity in the fluid control equation to explain the volume occupied by the landslide particles; the control equation of the landslide particles is Newton's second law, and the resultant force of water and air acting on the particles is added in the control equation.
3. The invention can respectively simulate two reservoir landslide and surge modes of granular material landslide and dangerous rock body collapse, and solves the defects that other numerical simulation methods can only simulate one instability mode and have a small application range in practical situations. Different reservoir bank side slope models are established for the two reservoir landslide surge modes, different calculation modes are set, the process of propagation of reservoir landslide surge disasters can be reflected more truly, and the application range is expanded.
Drawings
The accompanying drawings, which are included to provide a further understanding of the embodiments of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. In the drawings:
FIG. 1 is a flow chart in an embodiment of the invention;
FIG. 2 is a graph showing the results of numerical simulation of granular landslide in an embodiment of the present invention;
FIG. 3 is a diagram showing the result of numerical simulation of collapse of dangerous rock mass in the embodiment of the invention;
fig. 4 is a block diagram of a system in an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to examples and accompanying drawings, and the exemplary embodiments and descriptions thereof are only used for explaining the present invention and are not meant to limit the present invention.
Example 1: a numerical simulation method for landslide and swell of a reservoir is shown in figure 1 and is specifically realized by the following steps.
Step S1: and acquiring topographic data before and after occurrence of reservoir landslide and surge disasters, and constructing a geometric model.
The topographic data acquisition process before and after the occurrence of the reservoir landslide surge disaster specifically comprises the following steps: collecting terrain elevation data by adopting an unmanned aerial vehicle; converting the terrain elevation data into terrain raster data; and subtracting the terrain raster data before and after the disaster occurs to obtain the slope instability range.
The construction process of the geometric model specifically comprises the following steps: according to terrain raster data before and after a disaster occurs and a slope instability range, a preprocessing script pre.sh is operated at a terminal, and a three-dimensional stl geometric model of an integral area and a instability area after the landslide and the surge of the reservoir is generated through a format conversion program; and for generalized regular terrain, directly generating stl geometric models by adopting three-dimensional drawing software, and then importing and calculating.
Step S2: reservoir flow field areas are divided in the geometric model by utilizing OpenFOAM software, and fluid physical mechanical parameters are configured.
The dividing process of the reservoir flow field area specifically comprises the following steps: determining the water level elevation, wherein the area below the water level elevation in the geometric model is a water body, and the area above the water level elevation is air; the motion states of water and air in the process of reservoir landslide and surge are described by solving a Navier-Stokes equation; solving an interface between two-phase fluids by adopting a VOF model; and, accounting for the effect of porosity in the fluid control equation to account for the volume occupied by the landslide particles.
The expression of the fluid control equation is specifically:
wherein alpha is p Represents porosity;α i represents the volume fraction of air or water; u. u f Representing the fluid velocity; u. of c Represents the compression speed; ρ is a unit of a gradient f Represents the fluid density; r pf Representing the calculation of momentum exchange items based on a drag force model in each grid unit; τ represents the stress tensor of the fluid; f σ Represents the surface tension; when water and air are simultaneously present in the grid cells, the average characteristics of the two-phase fluid are used to characterize the cells.
In addition, the fluid physical mechanical parameters include, but are not limited to, air density, air kinematic viscosity, water density, and water kinematic viscosity.
And step S3: and establishing a bank side slope model on the geometric model by utilizing LIGGGHTS software, and configuring particle physical mechanical parameters and particle grading according to the physical properties of the bank side slope rock-soil body.
The configuration process of the physical and mechanical parameters and the grain composition of the grains is specifically as follows: obtaining the density, young modulus, poisson ratio, friction coefficient, coefficient of restitution and rolling friction coefficient of the rock-soil mass of the bank side slope through a geotechnical test or trial calculation mode, and analyzing to obtain the particle grading of the rock-soil mass of the bank side slope; setting particles with different particle sizes according to a preset similarity ratio for simulation according to the particle composition of the rock-soil mass of the bank slope; and calling different operation commands to respectively simulate granular landslide and dangerous rock body collapse.
Specifically, loose particles with different particle sizes are set for simulating the landslide of the loose particles; for dangerous rock body collapse, a plurality of small particles are subjected to rigid bonding to form a large-size block body without overlapping quantity for simulation.
The commands for setting different bank slope models are as follows:
1. a three-dimensional geometric model of the destabilized region is imported.
2. The particle properties are defined.
3. Setting the grain composition. When simulating the sliding slope of the granular body, n is the number of the type of the granule, c i The number of particle types is a percentage of the total number of particles. When dangerous rock mass collapse is simulated, n =1, c =1 without considering the grain composition.
4. Particle packing is provided.
5. A calculation mode is set.
For bulk landslide, considering the landslide source as an independent bulk particle, NVE integration was performed to update the position, velocity and angular velocity of all particles in the simulation.
For dangerous rock body collapse, all object source particles are combined into an independent rigid body, the mass center and the moment of inertia of the formed rigid body can be calculated after the size and the position information of all the particles are counted, the force and the moment applied to the rigid body are the sum of the force and the moment on the particles formed by the rigid body, and the rigid body is used as a single entity to move and rotate.
In the process of reservoir landslide and surge, the motion state of a landslide body is described by solving a Newton's second law, the motion state comprises a translation part and a rotation part, besides acting force among particles, resultant force of water and air acting on the particles is added into a control equation, and the control equation is as follows:
wherein u is i Represents the translational velocity of particle i; f ij Representing the interaction force of particle i and particle j; f i f Representing the resultant force of the water and air acting on the particles i. I is i Is the moment of inertia; theta i Represents the angular velocity of rotation of the particle i; m ij Represents the moment of particle j acting on particle i;
representing the moment of the water and air acting on the particle i.
And step S4: and performing analog simulation on the bank slope model by a CFDEM coupling method to obtain a simulation calculation result.
The process of performing analog simulation on the bank slope model by the CFDEM coupling method specifically comprises the following steps: (1) The DEM solver calculates the forces acting on the particles and solves the particle momentum equation to update the position and velocity of each particle; (2) After DEM circulation is completed, calculating the porosity and fluid-particle interaction force in the CFD grid unit, and transmitting the interaction force to the CFD circulation; (3) Solving a fluid momentum equation and a continuity control equation to obtain a fluid velocity field and a fluid pressure field; (4) After the flow field is converged, calculating the fluid acting force acting on the particles and returning the fluid acting force to the DEM solver; (5) And circulating the CFD-DEM calculation until the calculation end time is reached.
Step S5: and analyzing the derived simulation calculation result by using visual post-processing software to obtain the slope movement process, the surge height and/or the influence range.
And specifically, after the terminal runs, processing a script post.sh, and generating a particle information file and a water surface information file at each moment by using a written format conversion program.
1. The reservoir landslide and swell mode is taken as the granular landslide for example.
Setting the water level elevation to 1840m, wherein the area below the water level elevation in the calculation model is a water body, and the area above the water level elevation is air. The method specifically comprises the following steps of: air density 1Kg/m 3 Air kinematic viscosity 1.48X 10 -5 m 2 S, water density 1000Kg/m 3 The kinematic viscosity of the water body is 1.0 multiplied by 10 -6 m 2 /s。
For the granular landslide, loose particles with different particle sizes are set for simulation. The values of the physical and mechanical parameters of the particles are as follows: density 1730kg/m of rock-soil mass of bank side slope 3 Young's modulus 2.5X 10 7 Pa, poisson's ratio of 0.35, coefficient of friction of 0.5, coefficient of restitution of 0.8, and coefficient of rolling friction of 0.01. According to the particle grading of rock-soil mass of bank side slope, setting particles with different particle sizes (1-4 m) according to a certain similarity ratio for simulation.
Fig. 2 shows the simulation results of the dynamic evolution process of granular landslide motion and surge propagation, which show the landslide and the surge morphology when t =0s, t =5s, t =10s, and t =15s, respectively. When t =5s or so, under the impact action of landslide particles, surge begins to be gradually formed; t =10s, most landslide particles slide down, the surge height is further increased, and the particles are diffused outwards in an arc shape along the landslide occurrence point; t =15s, the landslide particles all enter the reservoir, the water body near the bank sinks along the movement direction of the landslide particles, and the swell begins to climb on the opposite bank. As can be seen from the figure, the method provided by the invention can reliably simulate the process of generating surge due to instability of the landslide.
2. The reservoir landslide and surge mode is taken as an example for dangerous rock collapse.
Setting the water level elevation to be 0.2m, wherein the area below the water level elevation in the calculation model is a water body, and the area above the water level elevation is air. The method specifically comprises the following steps of: air density 1Kg/m 3 Air kinematic viscosity 1.48X 10 -5 m 2 S, water density 1000Kg/m 3 The kinematic viscosity of the water body is 1.0 multiplied by 10 -6 m 2 /s。
For dangerous rock body collapse, a plurality of small particles are subjected to rigid bonding to form a large-size block body without overlapping quantity for simulation. The values of the physical and mechanical parameters of the particles are as follows: bank side slope rock-soil mass density 1730kg/m 3 Young's modulus of 2.5X 10 7 Pa, poisson's ratio of 0.35, coefficient of friction of 0.5, coefficient of restitution of 0.8, and coefficient of rolling friction of 0.01.
The simulation result of the collapse movement of the dangerous rock mass and the dynamic evolution process of the surge propagation is shown in fig. 3, which shows the shapes of the bank side slope and the surge when the bank side slope rock slides into the water in a whole block form, and t =0s, t =1s, t =1.3s and t =1.6s in the figure. By adopting the method provided by the invention, the process of surging caused by collapse of dangerous rock bodies can be simulated more reliably.
Example 2: a reservoir landslide swell numerical simulation system for implementing the reservoir landslide swell numerical simulation method described in embodiment 1, as shown in fig. 4, includes a geometric modeling module, a region division module, a slope modeling module, a coupling analysis module, and a simulation processing module.
The geometric modeling module is used for acquiring topographic data before and after a reservoir landslide and surge disaster occurs and constructing a geometric model; the region division module is used for dividing a reservoir flow field region in the geometric model by utilizing OpenFOAM software and configuring fluid physical mechanical parameters; the side slope modeling module is used for establishing a bank side slope model on the geometric model by utilizing LIGGGHTS software and configuring particle physical mechanical parameters and particle gradation according to the physical properties of the bank side slope rock and soil mass; the coupling analysis module is used for performing analog simulation on the bank slope model by a CFDEM coupling method to obtain a simulation calculation result; and the simulation processing module is used for analyzing the derived simulation calculation result by using the visual post-processing software to obtain the slope movement process, the surge height and/or the influence range.
The working principle is as follows: when discrete particles are processed, a discrete unit method (DEM) has good flexibility and can well reflect the processes of unstable sliding and fragmentation of a side slope, so that the DEM is widely applied to numerical simulation of a geological disaster dynamic process. The Computational Fluid Dynamics (CFD) method has its unique advantages in terms of computational accuracy of the fluid evolution process. According to the method, through the coupling solution between the two methods, the defects of the traditional simulation method in the research of the fluid-solid coupling phenomenon can be overcome, the interaction among a landslide body, reservoir water and air in the landslide motion process is considered, and the real evolution process of the reservoir landslide under the influence of the environment is restored. By establishing different reservoir bank side slope models, numerical calculation of two instability modes of granular break and dangerous rock collapse can be realized, and certain reliability is achieved when real landslide surge disasters with large simulation time and space scales are achieved.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and so forth) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above embodiments, objects, technical solutions and advantages of the present invention are further described in detail, it should be understood that the above embodiments are only examples of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A reservoir landslide surge numerical simulation method is characterized by comprising the following steps:
acquiring topographic data before and after occurrence of reservoir landslide and surge disasters, and constructing a geometric model;
dividing a reservoir flow field area in a geometric model by utilizing OpenFOAM software, and configuring fluid physical mechanical parameters;
establishing a bank side slope model on the geometric model by using LIGGGHTS software, and configuring particle physical mechanical parameters and particle grading according to the physical properties of the bank side slope rock-soil body;
performing analog simulation on the bank slope model by a CFDEM coupling method to obtain a simulation calculation result;
and analyzing the derived simulation calculation result by using visual post-processing software to obtain the slope movement process, the surge height and/or the influence range.
2. The method according to claim 1, wherein the topographic data acquisition process before and after occurrence of the reservoir landslide and surge disaster is as follows:
collecting terrain elevation data by adopting an unmanned aerial vehicle;
converting the terrain elevation data into terrain raster data;
and subtracting the terrain raster data before and after the disaster occurs to obtain the slope instability range.
3. The method according to claim 2, wherein the geometric model is constructed by the following steps:
according to terrain raster data before and after a disaster occurs and a slope instability range, a pre-processing script pre.sh is operated at a terminal, and a three-dimensional stl geometric model of an integral area and a instability area after the landslide and the swell of a reservoir is generated through a format conversion program;
and for generalized regular terrain, directly generating stl geometric models by adopting three-dimensional drawing software, and then importing and calculating.
4. The method for simulating the landslide and swell value of the reservoir according to claim 1, wherein the dividing process of the reservoir flow field area is as follows:
determining the water level elevation, wherein the area below the water level elevation in the geometric model is a water body, and the area above the water level elevation is air;
the motion states of water and air in the landslide and swell process of the reservoir are described by solving a Navier-Stokes equation;
solving an interface between two-phase fluids by adopting a VOF model;
and, accounting for the effect of porosity in the fluid control equation to account for the volume occupied by the landslide particles.
5. The method according to claim 4, wherein the fluid control equation has an expression as follows:
wherein alpha is p Represents porosity; alpha is alpha i Represents the volume fraction of air or water; u. u f Representing the fluid velocity; u. of c Represents the compression speed; rho f Represents the fluid density; r is pf Representing the calculation of momentum exchange items based on a drag force model in each grid unit; τ represents the stress tensor of the fluid; f σ Represents a surface tension; when water and air exist in the grid cells at the same time, the average characteristics of the two-phase fluid are used for characterizing the cells.
6. The method of claim 1, wherein the hydrophysical parameters include air density, air kinematic viscosity, water density, and water kinematic viscosity.
7. The method for simulating the numerical value of the swell of the landslide of the reservoir as claimed in claim 1, wherein the configuration process of the physical and mechanical parameters and the particle composition of the particles is specifically as follows:
obtaining the density, young modulus, poisson ratio, friction coefficient, coefficient of restitution and rolling friction coefficient of the rock-soil mass of the bank slope through a geotechnical test or trial calculation mode, and analyzing to obtain the grain composition of the rock-soil mass of the bank slope;
setting particles with different particle sizes according to a preset similarity ratio for simulation according to the particle composition of the rock-soil mass of the bank slope;
and calling different operation commands to respectively simulate granular landslide and dangerous rock collapse.
8. The method for simulating the numerical value of the swell of the landslide of the reservoir according to claim 7, wherein the simulation of the landslide of the dispersed particles and the collapse of the dangerous rock mass is specifically as follows:
for the granular body landslide, loose particles with different particle sizes are set for simulation;
for dangerous rock body collapse, a plurality of small particles are subjected to rigid bonding to form a large-size block body without overlapping quantity for simulation.
9. The method for simulating the numerical value of the swell of the landslide of the reservoir according to claim 1, wherein the simulation process of the bank slope model by the CFDEM coupling method specifically comprises the following steps:
the DEM solver calculates each force acting on the particles and solves a particle momentum equation to update the position and the speed of each particle;
after DEM circulation is completed, calculating the porosity and fluid-particle interaction force in the CFD grid unit, and transmitting the interaction force to the CFD circulation;
solving a fluid momentum equation and a continuity control equation to obtain a fluid velocity field and a fluid pressure field;
after the flow field is converged, calculating the fluid acting force acting on the particles and returning the fluid acting force to the DEM solver;
and circulating the CFD-DEM calculation until the calculation end time is reached.
10. The utility model provides a reservoir landslide numerical simulation system that surges, characterized by includes:
the geometric modeling module is used for acquiring topographic data before and after a reservoir landslide and surge disaster occurs and constructing a geometric model;
the region division module is used for dividing a reservoir flow field region in the geometric model by utilizing OpenFOAM software and configuring fluid physical mechanical parameters;
the side slope modeling module is used for establishing a bank side slope model on the geometric model by utilizing LIGGGHTS software and configuring particle physical mechanical parameters and particle grading according to the physical properties of the bank side slope rock-soil body;
the coupling analysis module is used for performing analog simulation on the bank slope model by a CFDEM coupling method to obtain a simulation calculation result;
and the simulation processing module is used for analyzing the derived simulation calculation result by using the visual post-processing software to obtain the slope movement process, the surge height and/or the influence range.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211547962.XA CN115906256A (en) | 2022-12-05 | 2022-12-05 | Reservoir landslide surge numerical simulation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211547962.XA CN115906256A (en) | 2022-12-05 | 2022-12-05 | Reservoir landslide surge numerical simulation method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115906256A true CN115906256A (en) | 2023-04-04 |
Family
ID=86496771
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211547962.XA Pending CN115906256A (en) | 2022-12-05 | 2022-12-05 | Reservoir landslide surge numerical simulation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115906256A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117150962A (en) * | 2023-08-30 | 2023-12-01 | 中国地质大学（武汉） | Method and device for calculating surging probability of random sampling and storage device |
| CN117195766A (en) * | 2023-08-29 | 2023-12-08 | 中国地质大学（武汉） | Landslide-surge climbing disaster override probability evaluation method, landslide-surge climbing disaster override probability evaluation equipment and storage equipment |
| CN117272703A (en) * | 2023-11-23 | 2023-12-22 | 中国电建集团西北勘测设计研究院有限公司 | Dynamic analysis method and system for landslide surge in reservoir area |
| CN117422006A (en) * | 2023-06-26 | 2024-01-19 | 中国矿业大学(北京) | Road collapse disaster deduction simulation method based on CFD-DEM and FEM combined algorithm |
| CN118036420A (en) * | 2024-02-26 | 2024-05-14 | 武汉大学 | Landslide surge impact floating vessel disaster chain estimation method and system |
| CN118036420B (en) * | 2024-02-26 | 2024-07-16 | 武汉大学 | Landslide surge impact floating vessel disaster chain estimation method and system |
-
2022
- 2022-12-05 CN CN202211547962.XA patent/CN115906256A/en active Pending
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117422006A (en) * | 2023-06-26 | 2024-01-19 | 中国矿业大学(北京) | Road collapse disaster deduction simulation method based on CFD-DEM and FEM combined algorithm |
| CN117195766A (en) * | 2023-08-29 | 2023-12-08 | 中国地质大学（武汉） | Landslide-surge climbing disaster override probability evaluation method, landslide-surge climbing disaster override probability evaluation equipment and storage equipment |
| CN117195766B (en) * | 2023-08-29 | 2024-05-10 | 中国地质大学（武汉） | Landslide-surge climbing disaster override probability evaluation method, landslide-surge climbing disaster override probability evaluation equipment and storage equipment |
| CN117150962A (en) * | 2023-08-30 | 2023-12-01 | 中国地质大学（武汉） | Method and device for calculating surging probability of random sampling and storage device |
| CN117150962B (en) * | 2023-08-30 | 2024-04-02 | 中国地质大学（武汉） | Method and device for calculating surging probability of random sampling and storage device |
| CN117272703A (en) * | 2023-11-23 | 2023-12-22 | 中国电建集团西北勘测设计研究院有限公司 | Dynamic analysis method and system for landslide surge in reservoir area |
| CN117272703B (en) * | 2023-11-23 | 2024-02-20 | 中国电建集团西北勘测设计研究院有限公司 | Dynamic analysis method and system for landslide surge in reservoir area |
| CN118036420A (en) * | 2024-02-26 | 2024-05-14 | 武汉大学 | Landslide surge impact floating vessel disaster chain estimation method and system |
| CN118036420B (en) * | 2024-02-26 | 2024-07-16 | 武汉大学 | Landslide surge impact floating vessel disaster chain estimation method and system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115906256A (en) | Reservoir landslide surge numerical simulation method and system | |
| Lin et al. | Evolution of the large landslide induced by Typhoon Morakot: A case study in the Butangbunasi River, southern Taiwan using the discrete element method | |
| Pitman et al. | Computing granular avalanches and landslides | |
| Cuomo et al. | Modelling the spatio-temporal evolution of a rainfall-induced retrogressive landslide in an unsaturated slope | |
| CN102708227B (en) | SPH (smoothed particle hydrodynamics) algorithm-based simulation method and simulation system of process of breaking dam by flood | |
| CN109033537B (en) | Calculation method and system for numerical simulation in rock-fill concrete pouring process | |
| Liu et al. | Numerical simulation of the 1995 rainfall-induced Fei Tsui Road landslide in Hong Kong: new insights from hydro-mechanically coupled material point method | |
| CN112417715B (en) | Rock mass fracture simulation method and system under true triaxial servo loading state | |
| CN110750864B (en) | Dam break flood demonstration early warning method and platform based on BIM technology | |
| CN114021487B (en) | Early warning method, device and equipment for landslide collapse and readable storage medium | |
| CN108984829B (en) | Calculation method and system for stacking process of rock-fill concrete rock-fill body | |
| AU2017227323B2 (en) | Particle simulation device, particle simulation method, and particle simulation program | |
| Al Ibrahim et al. | Particula: A simulator tool for computational rock physics of granular media | |
| Hu et al. | 3D modelling of soil-rock mixtures considering the morphology and fracture characteristics of breakable blocks | |
| Wang et al. | Numerical simulation of the fluid–solid coupling mechanism of internal erosion in granular soil | |
| Zhao et al. | Numerical investigation of the failure mechanisms of soil–rock mixture slopes by material point method | |
| CN111931369A (en) | Rainfall type landslide stability analysis and movement distance measuring and calculating method, equipment and medium | |
| An et al. | Landslide dynamic process and parameter sensitivity analysis by discrete element method: the case of Turnoff Creek rock avalanche | |
| CN109522589B (en) | Non-analytic method for simulating two-phase flow of pipeline particles and electronic equipment | |
| CN113255197B (en) | Physical mud flow erosion and deposition distance simulation method based on discrete element method | |
| CN114707376A (en) | Method for simulating occurrence distribution of seabed polymetallic tuberculosis | |
| CN112241603B (en) | Numerical simulation method for high-order landslide impact scraping and underlayer converging process | |
| CN110378023A (en) | Roofing rainfall runoff simulation method, apparatus, computer equipment and storage medium | |
| Su et al. | A new GPU-accelerated coupled discrete element and depth-averaged model for simulation of flow-like landslides | |
| Das et al. | Modeling granular particle shape using discrete element method |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Xiao Huabo Inventor after: Shi Weiming Inventor after: Wang Zehao Inventor after: Wang Dongpo Inventor after: Shi Dingguo Inventor after: Liu Shiyong Inventor before: Xiao Huabo Inventor before: Shi Weiming Inventor before: OuYang Chaojun Inventor before: Wang Zehao Inventor before: Wang Dongpo Inventor before: Shi Dingguo Inventor before: Liu Shiyong |