Method for predicting slip volume of single landslide

Abstract

The invention discloses a slip volume prediction method for a single landslide, and relates to the technical field of geological disasters. The method comprises the following steps: respectively inputting a digital elevation model (DEM) with any resolution and a landslide slippage boundary range; respectively preprocessing the input data; establishing a steepest slope line in a slip boundary range; generating a normal vector plane at each grid point; calculating coordinates and an elevation value of an intersection point of the normal vector plane and the slippage boundary; predicting the depth of each grid slip plane; and carrying out landslide slippage volume prediction. When the method is applied to landslide slip surface and slip volume prediction, landslide damage and influence ranges can be determined, potential landslide danger areas can be discriminated, and scientific support can be provided for landslide disaster prevention and control and personnel disaster prevention and risk avoiding strategy planning and management.

Description

Technical Field

[0001] The present application relates to the technical field of geological disasters, and in particular to a method for predicting the sliding volume of a single landslide. Background Art

[0002] Single landslide disaster is one of the common geological disasters in the western mountainous areas of my country. Under the influence of various inducing factors such as rainfall, earthquakes and human engineering activities, it can become unstable and slide under certain conditions, causing related losses. Accurately predicting the depth and volume of the sliding surface of a single landslide is the key to the construction of post-disaster prevention and control projects and the calculation of construction volume.

[0003] At present, the main methods for predicting the sliding surface of a single landslide are the method based on the equilibrium profile and the method based on the ellipsoid assumption. However, the method based on the equilibrium profile obtains the sliding surface by calculating the geometric characteristics of the landslide profile, which requires a large amount of drilling and geophysical exploration field data to assist in the calculation, and its applicability is limited to a certain extent; the method based on the ellipsoid assumption pre-assumes that the shape of the sliding body is an ellipsoid, and estimates the sliding surface through geometric calculation, which is not very adaptable to the actual sliding terrain.

[0004] Therefore, a method for predicting the sliding volume of a single landslide is urgently needed to improve its adaptability to the actual sliding terrain. Summary of the invention

[0005] Based on this, it is necessary to provide a method for predicting the sliding volume of a single landslide in response to the above technical problems.

[0006] The present invention adopts the following technical solutions: The present invention provides a method for predicting the sliding volume of a single landslide, comprising: Acquire digital elevation model data of a single landslide, input it into the sliding boundary range of the single landslide, and obtain sliding boundary range data; Preprocess the digital elevation model data and sliding boundary range data of the single landslide, project the preprocessed digital elevation model data and sliding boundary range data to the same coordinate system, and extract the elevation value of the grid center point and the center point digital elevation model data of each grid pixel in the coordinate system; The grid center points of all grid pixels and the two grid center points corresponding to the highest and lowest values ​​in the elevation values ​​of the center point digital elevation model data are connected to create the steepest slope line within the sliding boundary range; a normal vector plane perpendicular to the steepest slope line within the sliding boundary range is created through the grid center point of each grid pixel; the coordinates of the intersection of the normal vector plane and the landslide sliding boundary range and the elevation value of the digital elevation model data of the intersection are calculated; The distance between the intersection of the normal vector plane and the landslide sliding boundary range and the steepest slope line within the sliding boundary range is calculated and cubic spline interpolation is performed to predict the depth of the landslide sliding surface at each grid point. Based on the predicted depth of the landslide slip surface at each grid point, the landslide slip volume of the single landslide is determined.

[0007] Preferably, extracting the grid center point of each grid pixel in the digital elevation model data in the coordinate system and the elevation value of the center point digital elevation model data specifically includes: The .shp file is read through the shpread function of Matalb software and the data content is stored in a structure array to convert the file in the .shp format with the sliding boundary range into structure data; The coordinates of the digital elevation model data in .tif format and the slip boundary range data in .shp format are unified and projected into the WGS1984 UTM coordinate system; Extract the grid center point of each raster pixel and the elevation value of the center point digital elevation model data from the coordinate system.

[0008] Preferably, the depth of the landslide slip surface at each grid point is predicted, specifically including: The angle values ​​of any two angles in the angle range of 0° to 90° are used as the inclination values ​​of the left and right interpolation sections; A single-segment cubic spline interpolation curve is generated at the center of each grid point. The lowest point of the elevation value of the digital elevation model data on each interpolation curve is the slip surface depth of each grid center point.

[0009] Preferably, the landslide volume of a single landslide is determined based on the predicted landslide slip surface depth of each grid point, using the formula: ; In the formula, The input DEM is x, y Resolution in direction; i, j Respectively i Line j The grid center points, For the i Line j The DEM elevation value at the center of the grid. For the i Line j The predicted value of the slip surface depth at the center point of each grid.

[0010] The present invention provides a sliding volume prediction device for a single landslide, comprising: The data acquisition module is used to obtain the digital elevation model data of the single landslide, input it into the sliding boundary range of the single landslide, and obtain the sliding boundary range data; A data extraction module is used to preprocess the digital elevation model data and the sliding boundary range data of the single landslide, and project the preprocessed digital elevation model data and the sliding boundary range data to the same coordinate system, and extract the grid center point of each grid pixel in the digital elevation model data in the coordinate system and the elevation value of the center point digital elevation model data; The landslide steepest slope line and normal vector plane construction module is used to connect the grid center points of all grid pixels and the two grid center points corresponding to the highest and lowest values ​​in the elevation values ​​of the center point digital elevation model data to create the steepest slope line within the sliding boundary range; through the grid center point of each grid pixel, create a normal vector plane perpendicular to the steepest slope line within the sliding boundary range, and obtain the coordinates of the intersection of the normal vector plane and the landslide sliding boundary range and the elevation value of the digital elevation model data of the intersection; The prediction module is used to calculate the distance between the intersection of the normal vector plane and the landslide slip boundary range and the steepest slope line within the slip boundary range and perform cubic spline interpolation to predict the landslide slip surface depth of each grid point; based on the predicted landslide slip surface depth of each grid point, the landslide slip volume of the single landslide is determined.

[0011] The present invention also provides a computer-readable storage medium, wherein the storage medium stores a computer program, and when the computer program is executed by a processor, the above-mentioned method for predicting the sliding volume of a single landslide is implemented.

[0012] The present invention also provides a computer device, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for predicting the sliding volume of a single landslide when executing the program.

[0013] At least one of the above technical solutions adopted by the present invention can achieve the following beneficial effects: The invention provides a method for predicting the slip volume of a single landslide. The method comprises the following steps: connecting two grid center points corresponding to the highest value and the lowest value of digital elevation model data within a landslide slip boundary range to create a steepest slope line within the slip boundary range; creating a normal vector plane perpendicular to the steepest slope line within the slip boundary range through the grid center point of each grid pixel in the digital elevation model data; calculating the coordinates of the intersection point between the normal vector plane and the landslide slip boundary range and the elevation value of the digital elevation model data of the intersection point; calculating the distance between the intersection point between the normal vector plane and the landslide slip boundary range and the steepest slope line within the slip boundary range and performing cubic spline interpolation on the distance to predict the landslide slip surface depth of each grid point; and predicting the landslide slip volume based on the predicted landslide slip surface depth of each grid point.

[0014] The method only needs the digital elevation model (DEM) data of the landslide sliding area and the vector data of the sliding range boundary to predict the landslide sliding surface and sliding volume, avoiding the need for field data and improving the terrain applicability. The present invention can help determine the damage and impact range of landslide sliding, identify potential landslide hazard areas, and provide scientific support for landslide disaster prevention and control, and personnel disaster prevention and risk avoidance strategy planning and management. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] The drawings described herein are used to provide a further understanding of the present application and constitute a part of the present application. The illustrative embodiments of the present application and their descriptions are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings: Figure 1 A schematic diagram of a flow chart of a method for predicting sliding surface and sliding volume of a single landslide provided by the present invention; Figure 2 A digital elevation model (DEM) and sliding boundary range of the Qingliu landslide for a single landslide sliding surface and sliding volume prediction method provided by the present invention; Figure 3 The prediction result of the sliding surface depth of the Qingliu landslide provided by the method for predicting the sliding surface and sliding volume of a single landslide provided by the present invention; Figure 4 The sliding volume prediction result of Qingliu landslide provided by the present invention is a single landslide sliding surface and sliding volume prediction method; Figure 5 A schematic diagram of a sliding volume prediction device for a single landslide provided by the present invention; Figure 6 A schematic diagram of a computer device for implementing a method for predicting the sliding volume of a single landslide provided by the present invention. DETAILED DESCRIPTION

[0016] In order to make the purpose, technical scheme and advantages of the present invention clearer, the technical scheme of the present application will be clearly and completely described below in combination with the specific embodiments of the present invention and the corresponding drawings. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the specification, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0017] The technical solutions provided by various embodiments of the present application are described in detail below in conjunction with the accompanying drawings.

[0018] Figure 1 The figure is a flow chart of a method for predicting the sliding volume of a single landslide in the present invention, which specifically comprises the following steps: S101: Acquire digital elevation model data of a single landslide, input the data into a sliding boundary range of the single landslide, and obtain sliding boundary range data.

[0019] The data format of the digital elevation model data is tif; the data format of the slip boundary range data is shp.

[0020] The coordinate systems used for the digital elevation model data and the slip boundary range data are both the WGS1984 spatial geographic coordinate system.

[0021] S102: Preprocessing the digital elevation model data and the sliding boundary range data of the single landslide, projecting the preprocessed digital elevation model data and the sliding boundary range data to the same coordinate system, and extracting the elevation value of the grid center point and the center point digital elevation model data of each grid pixel in the coordinate system.

[0022] The .shp file is read through the shpread function of the Matlab software and the data content is stored in a structure array to convert the file in the .shp format with the sliding boundary range into structure data.

[0023] The coordinates of the digital elevation model data in .tif format and the slip boundary range data in .shp format are unified and projected into the WGS1984 UTM coordinate system.

[0024] Extract the elevation value of the grid center point and the center point digital elevation model data of each raster pixel in the coordinate system.

[0025] S103: Connect the grid center points of all grid pixels and the two grid center points corresponding to the highest value and the lowest value in the elevation values ​​of the center point digital elevation model data to create the steepest slope line within the sliding boundary range; create a normal vector plane perpendicular to the steepest slope line within the sliding boundary range through the grid center point of each grid pixel; calculate the coordinates of the intersection of the normal vector plane and the landslide sliding boundary range and the elevation value of the digital elevation model data of the intersection.

[0026] S104: Calculate the distance between the intersection of the normal vector plane and the landslide sliding boundary range and the steepest slope line within the sliding boundary range and perform cubic spline interpolation on it to predict the depth of the landslide sliding surface at each grid point.

[0027] Specifically, any two angle values ​​in the range of 0° to 90° are input as the inclination values ​​of the left and right interpolation profiles.

[0028] A single-segment cubic spline curve is generated at the center of each grid point, and the lowest point of the elevation value of the digital elevation model data on each interpolation curve is the slip surface depth of each grid center point.

[0029] S105: Predicting the landslide sliding volume based on the predicted landslide sliding surface depth of each grid point.

[0030] The formula for predicting the landslide sliding volume is: ; In the formula, The input DEM is x, y Resolution in direction; i, j Respectively i Line j The grid center points, For the i Line j The DEM elevation value at the center of the grid. For the i Line j The predicted value of the slip surface depth at the center point of each grid.

[0031] The specific embodiment of the present invention selects Qingliu landslide in a certain province to predict the sliding surface and sliding volume of a single landslide. Figure 1 Based on the embodiments, a method for predicting the sliding volume of a single landslide proposed in the present invention is described in detail.

[0032] See also Figure 2 , for the digital elevation model (DEM) and slip boundary range of Qingliu landslide, first input the digital elevation model (DEM) data with a spatial resolution of 30m×30m (data format is .tif), the DEM data source is the Earth Observation Satellite (https: / / earth.jaxa.jp / en / ), then input the slip boundary range of Qingliu landslide (data format is .shp). The coordinate system used by the input digital elevation model (DEM) and the slip boundary range data of the landslide is the WGS1984 spatial geographic coordinate system.

[0033] The digital elevation model (DEM) data was read using the shpread function of Matalb software and the data content was stored in a structure array. The coordinates of the digital elevation model (DEM) and the landslide boundary range data were unified and both were projected to the WGS1984 UTM 48N coordinate system; then, the grid center point and center point elevation value of each grid pixel in the DEM data were extracted (( ), a total of 1489 grid center points were extracted.

[0034] Among the 1488 grid center points, the highest DEM value of the grid center points within the landslide sliding boundary is 2430m, and the lowest is 1820m. The two grid center points with the highest and lowest DEM values ​​are connected to create the steepest slope line within the Qingliu landslide sliding boundary.

[0035] Through 1488 grid center points, 1488 normal vector planes perpendicular to the steepest slope line within the sliding boundary of Qingliu landslide are generated respectively.

[0036] The coordinates of the intersections of 1488 normal vector planes and the landslide boundary range as well as the DEM elevation values ​​of the intersections were calculated.

[0037] The distances between all intersection points and the steepest slope line were calculated and cubic spline interpolation (inverse distance weighted method) was performed. Any two angle values ​​(0°-90°) were input as the inclination values ​​of the left and right interpolation profiles. 25° and 15° were input as the inclination values ​​of the left and right interpolation profiles of the Qingliu landslide slip surface, respectively. Then, 1488 single-segment cubic spline curves were interpolated at all grid center points. The lowest point of the DEM elevation value on each interpolation curve ( ) is used as the slip surface depth of each grid center point, and the combined surface of the lowest DEM value of all grid center points is the slip surface depth prediction result, see Figure 3 , which is the prediction result of the sliding surface depth of Qingliu landslide.

[0038] Then the sliding volume of Qingliu landslide is predicted, and the calculation formula is as follows: ; In the above formula, The input DEM is x, y Resolution in direction; i, j Respectively i Line j The grid center points, For the i Line j The DEM elevation value at the center of the grid. For the i Line j The predicted value of the slip surface depth at the center point of each grid.

[0039] In a specific embodiment of the present invention, the input DEM is x, y Resolution in direction 、 The number of rows of grid center points generated is 30m and 30m respectively. i Total 31 rows, columns j A total of 48 lines, the predicted Qingliu landslide sliding volume of 1488 grid center points is accumulated, see Figure 4 , is the prediction result of the sliding volume of Qingliu landslide. The predicted total sliding volume is 24.8×10 4 m 3 .

[0040] The specific definition of a slip volume prediction device for a single landslide can be found in the above definition of a slip volume prediction method for a single landslide, which will not be repeated here. Each module in the above slip volume prediction device for a single landslide can be implemented in whole or in part by software, hardware, and a combination thereof. The above modules can be embedded in or independent of a processor in a computer device in the form of hardware, or can be stored in a memory in a computer device in the form of software, so that the processor can call and execute the operations corresponding to the above modules.

[0041] Figure 5 A schematic diagram of a sliding volume prediction device for a single landslide provided by the present invention comprises: The data acquisition module 601 is used to obtain the digital elevation model data of the single landslide, input it into the sliding boundary range of the single landslide, and obtain the sliding boundary range data; The data extraction module 602 is used to pre-process the digital elevation model data and the sliding boundary range data of the single landslide, and project the pre-processed digital elevation model data and the sliding boundary range data to the same coordinate system, and extract the grid center point of each grid pixel in the digital elevation model data in the coordinate system and the elevation value of the center point digital elevation model data; The landslide steepest slope line and normal vector plane construction module 603 is used to connect the grid center points of all grid pixels and the two grid center points corresponding to the highest value and the lowest value in the elevation value of the center point digital elevation model data to create the steepest slope line within the sliding boundary range; through the grid center point of each grid pixel, create a normal vector plane perpendicular to the steepest slope line within the sliding boundary range, and obtain the coordinates of the intersection of the normal vector plane and the landslide sliding boundary range and the elevation value of the digital elevation model data of the intersection; The prediction module 604 is used to calculate the distance between the intersection of the normal vector plane and the landslide slip boundary range and the steepest slope line within the slip boundary range and perform cubic spline interpolation to predict the depth of the landslide slip surface at each grid point; based on the predicted landslide slip surface depth at each grid point, determine the landslide slip volume of the single landslide.

[0042] The specific definition of a slip volume prediction device for a single landslide can be found in the above definition of a slip volume prediction method for a single landslide, which will not be repeated here. Each module in the above slip volume prediction device for a single landslide can be implemented in whole or in part by software, hardware, and a combination thereof. The above modules can be embedded in or independent of a processor in a computer device in the form of hardware, or can be stored in a memory in a computer device in the form of software, so that the processor can call and execute the operations corresponding to the above modules.

[0043] The present invention also provides a computer-readable storage medium, which stores a computer program, which can be used to execute the above Figure 1 A method for predicting the sliding volume of a single landslide is provided.

[0044] The present invention also provides Figure 6 The structural diagram of the computer device shown in FIG. Figure 6 As shown in the figure, at the hardware level, the computer device includes a processor, an internal bus, a network interface, a memory, and a non-volatile memory, and may also include other hardware required for the business. The processor reads the corresponding computer program from the non-volatile memory into the memory and then runs it to achieve the above Figure 1 A method for predicting the sliding volume of a single landslide is provided.

[0045] Those of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), tape, floppy disk, flash memory or optical memory, etc. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0046] The technical features of the above embodiments may be arbitrarily combined. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of the present invention.

[0047] The technical features of the above embodiments may be arbitrarily combined. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of the present invention.

Claims

1. A method for predicting the sliding volume of a single landslide, characterized in that: include: Acquire digital elevation model data of a single landslide, input it into the sliding boundary range of the single landslide, and obtain sliding boundary range data; Preprocess the digital elevation model data and sliding boundary range data of the single landslide, project the preprocessed digital elevation model data and sliding boundary range data to the same coordinate system, and extract the elevation value of the grid center point and the center point digital elevation model data of each grid pixel in the coordinate system; The grid center points of all grid pixels and the two grid center points corresponding to the highest and lowest values ​​in the elevation values ​​of the center point digital elevation model data are connected to create the steepest slope line within the sliding boundary range; a normal vector plane perpendicular to the steepest slope line within the sliding boundary range is created through the grid center point of each grid pixel, and the coordinates of the intersection of the normal vector plane and the landslide sliding boundary range and the elevation value of the digital elevation model data of the intersection are obtained; The distance between the intersection of the normal vector plane and the landslide sliding boundary range and the steepest slope line within the sliding boundary range is calculated and cubic spline interpolation is performed to predict the depth of the landslide sliding surface at each grid point. Based on the predicted depth of the landslide slip surface at each grid point, the landslide slip volume of the single landslide is determined.

2. The method for predicting the sliding volume of a single landslide according to claim 1, characterized in that: The step of extracting the grid center point of each grid pixel in the digital elevation model data in the coordinate system and the elevation value of the center point digital elevation model data specifically includes: The .shp file is read through the shpread function of Matalb software and the data content is stored in a structure array to convert the file in the .shp format with the sliding boundary range into structure data; The coordinates of the digital elevation model data in .tif format and the slip boundary range data in .shp format are unified and projected into the WGS1984 UTM coordinate system; Extract the grid center point of each raster pixel and the elevation value of the center point digital elevation model data from the coordinate system.

3. The method for predicting the sliding volume of a single landslide according to claim 1, characterized in that: The prediction of the depth of the landslide sliding surface at each grid point specifically includes: The angle values ​​of any two angles in the angle range of 0° to 90° are used as the inclination values ​​of the left and right interpolation sections; A single-segment cubic spline interpolation curve is generated at the center of each grid point. The lowest point of the elevation value of the digital elevation model data on each interpolation curve is the slip surface depth of each grid center point.

4. The method for predicting the sliding volume of a single landslide according to claim 1, characterized in that: The slip volume of a single landslide is determined based on the predicted depth of the landslide slip surface at each grid point. The formula is: ； In the formula, The input DEM is x, y Resolution in direction; i, j Respectively i Line j The grid center points, For the i Line j The DEM elevation value at the center of the grid. For the i Line j The predicted value of the slip surface depth at the center point of each grid.

5. A sliding volume prediction device for a single landslide, comprising: The data acquisition module is used to obtain the digital elevation model data of the single landslide, input it into the sliding boundary range of the single landslide, and obtain the sliding boundary range data; A data extraction module is used to preprocess the digital elevation model data and the sliding boundary range data of the single landslide, and project the preprocessed digital elevation model data and the sliding boundary range data to the same coordinate system, and extract the grid center point of each grid pixel in the digital elevation model data in the coordinate system and the elevation value of the center point digital elevation model data; The landslide steepest slope line and normal vector plane construction module is used to connect the grid center points of all grid pixels and the two grid center points corresponding to the highest and lowest values ​​in the elevation values ​​of the center point digital elevation model data to create the steepest slope line within the sliding boundary range; through the grid center point of each grid pixel, create a normal vector plane perpendicular to the steepest slope line within the sliding boundary range, and obtain the coordinates of the intersection of the normal vector plane and the landslide sliding boundary range and the elevation value of the digital elevation model data of the intersection; The prediction module is used to calculate the distance between the intersection of the normal vector plane and the landslide slip boundary range and the steepest slope line within the slip boundary range and perform cubic spline interpolation to predict the depth of the landslide slip surface at each grid point; based on the predicted landslide slip surface depth at each grid point, the landslide slip volume of the single landslide is determined.

6. A computer-readable storage medium, characterized in that: The storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 4 is implemented.

7. A computer device, characterized in that: The method comprises a memory, a processor and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method described in any one of claims 1 to 4 is implemented.

Images