Search

CN-121999163-A - Interactive landslide power disaster-induced rapid simulation method

CN121999163ACN 121999163 ACN121999163 ACN 121999163ACN-121999163-A

Abstract

The application discloses an interactive landslide power disaster-induced rapid simulation method, which belongs to the technical field of geological disaster modeling and comprises the following steps of obtaining remote sensing images and terrain grids, defining the range of a sliding body area, determining the maximum thickness of a sliding surface, calculating the distance proportion of each grid unit and the defined boundary of the sliding body area, constructing a thickness space distribution model of each grid unit based on an exponential function, the maximum thickness and the distance proportion, generating a three-dimensional foundation terrain by combining surface topography, establishing a landslide motion control equation by using a universal control equation, and outputting a simulation result by using the landslide motion control equation through a simulation tool, wherein the space distribution model is used for fitting a picture curve based on n times of decay functions, automatically generating the space distribution of picture depth, improving the accuracy and convenience of initial conditions of the sliding body, and effectively improving the accuracy of the simulation result of the landslide power process by introducing an improved friction resistance evolution model taking sliding into consideration of flow conversion.

Inventors

  • GUO JIAN
  • CUI YIFEI

Assignees

  • 中国科学院地理科学与资源研究所

Dates

Publication Date
20260508
Application Date
20260127

Claims (7)

  1. 1. The interactive landslide power disaster-causing rapid simulation method is characterized by comprising the following steps of: S1, acquiring a remote sensing image covering an analog region and a terrain grid formed by two-dimensional grid units, encircling a sliding body region range in the remote sensing image, determining the maximum thickness h max of a sliding surface, and projecting an encircling result into the terrain grid; S2, calculating the distance proportion d between each grid unit and the circling boundary of the sliding body area; S3, constructing a spatial distribution model of the thickness of the sliding body based on an exponential function, the maximum thickness h max and the distance proportion d, and calculating the thickness h of each grid unit: S4, after the thickness h of each grid unit is determined, generating a matched three-dimensional bed topography z by combining the surface topography S of the sliding body area; s5, based on a shallow water equation, using a general control equation to establish a landslide motion control equation: ; ; wherein u is a set velocity vector, g is a gravitational acceleration term, τ is a resistance term, Slope of the topography of the three-dimensional bed; And S6, applying a landslide motion control equation through a simulation tool, setting a simulation time t, and then running to output a simulation result.
  2. 2. The interactive landslide power disaster-induced rapid simulation method according to claim 1, wherein the terrain grid acquisition mode is that digital elevation model DEM and remote sensing image data covering the simulation area are imported, grid resampling of space projection and DEM is carried out, and a terrain grid file is generated.
  3. 3. The interactive landslide power disaster quick simulation method of claim 1, wherein the valid mark representation and the invalid mark representation are performed on all grids in the terrain grid, wherein the valid mark represents a grid unit belonging to a sliding body area, and the invalid mark represents a background area not participating in calculation.
  4. 4. The interactive landslide power disaster-induced rapid simulation method of claim 3, wherein S2 is specifically: S2.1, calculating distance values from the grid units of each effective mark of the sliding body region to the nearest circled boundary, and extracting the maximum distance from the distance values of all the grid units of the effective mark as a normalization reference; and S2.2, dividing the distance value of each effectively marked grid unit by the maximum distance value to obtain a standardized distance proportion d.
  5. 5. The interactive landslide power disaster quick simulation method of claim 4, wherein the thickness calculation formula of each grid unit is: ; where n is an attenuation index, where n=1 indicates a linear attenuation, and a larger n indicates a smoother sliding surface and a closer to a rock landslide.
  6. 6. The interactive landslide power disaster quick simulation method according to claim 1, wherein a friction resistance evolution model is introduced into the resistance term τ: ; Wherein, the For the conversion factor, c represents the cohesion, Indicating the coefficient of friction of the slide material, Representing the normal stress of the grid cells, Representing the slip mass density, v representing the current speed of movement of the grid element, Representing the turbulence coefficient; the frictional resistance evolution model shows the following characteristics in the material movement of different stages: Stage I-friction dominant stage-the mass is at a speed below the low speed threshold and the mass is structurally complete, λ≡1, where the shear stress expression is: I.e., the frictional resistance evolution model is degenerated to the Mohr-Coulomb form; The phase II-mixing behavior phase is that when the speed of a substance is increased or the structure is damaged, 0< lambda <1, the friction resistance evolution model is represented by mixing control of cohesive force, friction force and inertia resistance, and the three forces act together; phase III-inertia dominant phase when the velocity of the material is above the high velocity threshold and crushing occurs, λ≡0, the shear stress expression is: The frictional resistance evolution model gradually approaches Voellmy resistance forms.
  7. 7. The interactive landslide power disaster-causing rapid simulation method of claim 6, wherein λ is a key parameter controlling multistage evolution, influenced by speed and movement distance, and calculated by automatic iteration of the following formula: ; Wherein B r is the landslide breaking degree, the value range is 0-1, V is the maximum threshold speed, and w is the weight factor.

Description

Interactive landslide power disaster-induced rapid simulation method Technical Field The invention belongs to the technical field of geological disaster modeling, and particularly relates to an interactive landslide dynamic disaster-induced rapid simulation method. Background Quick geological disasters such as landslide and debris flow frequently occur in the global mountain region, and the method has the characteristics of strong burst, wide hazard range, complex power process and the like. In order to study the exercise mechanism and serve risk assessment and engineering prevention, various numerical simulation tools such as RAMMS, massFlow and the like have been developed at home and abroad. The software is generally based on shallow water equations or continuous medium dynamic frames, and calculates the motion evolution process of the sliding body along the surface of the terrain by inputting a Digital Elevation Model (DEM), initial sliding body thickness and friction parameters, and outputs results such as a speed field, thickness distribution, accumulation range and the like. However, the above numerical simulations have the following drawbacks: The accurate three-dimensional sliding surface morphology is an important premise for determining the landslide hazard scale and a core boundary condition for calculation, but because the landslide hazard scale is covered by the rock-soil vegetation on the ground, the accurate identification and characterization are very difficult, which becomes an important challenge for influencing the accuracy of a simulation result. In practice, the landslide is started by a solid state and gradually chipped to perform fluidization movement, the existing landslide simulation mostly adopts a single resistance model in a fixed form, and consideration of a chip fluidization and friction weakening mechanism in the landslide movement process is lacked, so that the damage range, speed and other dynamic parameters of a simulation result are inaccurate. Disclosure of Invention In order to solve the problems in the background technology, the invention provides an interactive landslide power disaster-induced rapid simulation method, which aims to solve the problems of inaccurate object source scale and inaccurate motion calculation in the existing numerical simulation process. In order to achieve the above purpose, the present invention provides the following technical solutions: The interactive landslide power disaster-causing rapid simulation method comprises the following steps: s1, acquiring a remote sensing image covering an analog region and a terrain grid formed by two-dimensional grid units, encircling a sliding body region range in the remote sensing image, determining the maximum thickness hmax of a sliding surface, and projecting an encircling result into the terrain grid; S2, calculating the distance proportion d between each grid unit and the circling boundary of the sliding body area; S3, constructing a spatial distribution model of the thickness of the sliding body based on an exponential function, the maximum thickness hmax and the distance proportion d, and calculating the thickness h of each grid unit: S4, after the thickness h of each grid unit is determined, generating a matched three-dimensional bed topography z by combining the surface topography S of the sliding body area; s5, based on a shallow water equation, using a general control equation to establish a landslide motion control equation: ; ; wherein u is a set velocity vector, g is a gravitational acceleration term, τ is a resistance term, Slope of the topography of the three-dimensional bed; And S6, applying a landslide motion control equation through a simulation tool, setting a simulation time t, and then running to output a simulation result. Preferably, the terrain grid acquisition mode is that a digital elevation model DEM and remote sensing image data covering an analog area are imported, space projection and grid resampling of the DEM are carried out, and a terrain grid file is generated. Preferably, the valid mark representation and the invalid mark representation are performed on all grids in the terrain grid, wherein the valid mark represents a grid unit belonging to the sliding body area, and the invalid mark represents a background area which does not participate in calculation. Preferably, S2 is specifically: S2.1, calculating distance values from the grid units of each effective mark of the sliding body region to the nearest circled boundary, and extracting the maximum distance from the distance values of all the grid units of the effective mark as a normalization reference; and S2.2, dividing the distance value of each effectively marked grid unit by the maximum distance value to obtain a standardized distance proportion d. Preferably, the thickness calculation formula of each grid unit is: ; Where n is an attenuation index, and when n=1, it represents linear attenuation, and when n is lar