A method for predicting the sliding volume of a single landslide

Abstract

The present invention discloses a method for predicting the sliding volume of a single landslide, which relates to the technical field of geological disasters. The method includes: respectively inputting digital elevation models (DEMs) with arbitrary resolutions and the landslide sliding boundary range; respectively preprocessing the input data; creating the steepest slope lines within the sliding boundary range; generating normal vector planes at each grid point; calculating the intersection coordinates and elevation values of the normal vector planes and the sliding boundary; predicting the depth of each grid sliding surface; and predicting the landslide sliding volume. Applying this method to the prediction of the landslide sliding surface and sliding volume can help determine the landslide failure and influence range, identify potential landslide hazard areas, and provide scientific support for landslide disaster prevention and control, as well as the planning and management of personnel disaster prevention and avoidance strategies.

Description

Technical Field

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

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

[0003] Currently, 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 hypothesis. However, the method based on the equilibrium profile calculates the sliding surface by calculating the geometric characteristics of the landslide profile, and requires a large amount of drilling and geophysical exploration field data for auxiliary calculation, so its applicability has certain limitations; the method based on the ellipsoid hypothesis pre-assumes the shape of the sliding mass as an ellipsoid and estimates the sliding surface through geometric calculations, and its adaptability to the actual sliding terrain is not strong.

[0004] Therefore, there is an urgent need for a method for predicting the sliding volume of a single landslide 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 view of the above technical problems.

[0006] The present invention adopts the following technical solutions:

[0007] The present invention provides a method for predicting the sliding volume of a single landslide, including:

[0008] 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;

[0009] Preprocess the digital elevation model data and the sliding boundary range data of the single landslide, project the preprocessed digital elevation model data and the sliding boundary range data onto the same coordinate system, and extract the grid center points of each grid cell in the coordinate system and the elevation values of the center point digital elevation model data;

[0010] Connect the two grid center points corresponding to the highest value and the lowest value among the grid center points of all grid cells and the elevation values of the center point digital elevation model data to create the steepest slope line within the sliding boundary range; for each grid center point of each grid cell, create a normal vector plane perpendicular to the steepest slope line within the sliding boundary range; calculate the intersection coordinates of the normal vector plane and the landslide sliding boundary range and the elevation values of the digital elevation model data at the intersection points;

[0011] Calculate the distance between the intersection point of the plane of the calculated normal vector and the landslide slip boundary range and the steepest slope line within the landslide slip boundary range, and perform cubic spline interpolation on it to predict the depth of the landslide slip surface at each grid point;

[0012] Based on the predicted depth of the landslide slip surface at each grid point, determine the landslide slip volume of the individual landslide.

[0013] Preferably, extract the grid center point of each grid cell in the digital elevation model data in the coordinate system and the elevation value of the center point digital elevation model data, specifically including:

[0014] Read the.shp file through the shpread function of Matalb software and store the data content in a structure array to convert the file with the landslide slip boundary range in.shp format into structure data;

[0015] Unify the coordinates of the digital elevation model data in.tif format and the landslide slip boundary range data in.shp format, and project them under the WGS1984 UTM coordinate system;

[0016] Extract the grid center point of each grid cell and the elevation value of the center point digital elevation model data from the coordinate system.

[0017] Preferably, predict the depth of the landslide slip surface at each grid point, specifically including:

[0018] Take any two angle values within the angle range of 0° to 90° as the dip angle values of the left and right interpolation profiles;

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

[0020] Preferably, based on the predicted depth of the landslide slip surface at each grid point, determine the landslide slip volume of the individual landslide, and the formula is:

[0021] ;

[0022] In the formula, are the resolutions of the input DEM in the x, y directions respectively; i, j are the i th row and the j th grid center point respectively, is the DEM elevation value at the i th row and the j th grid center point, is the i th row and the jThe predicted value of the depth of the slip surface at the center point of each grid.

[0023] The present invention provides a device for predicting the sliding volume of a single landslide, including:

[0024] A data acquisition module, configured to obtain 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;

[0025] A data extraction module, configured to preprocess the digital elevation model data and the sliding boundary range data of the single landslide, project the preprocessed digital elevation model data and the sliding boundary range data onto the same coordinate system, and extract the grid center point of each grid cell in the digital elevation model data in the coordinate system and the elevation value of the center point digital elevation model data;

[0026] A landslide steepest slope line and normal vector plane construction module, configured to connect the two grid center points corresponding to the highest value and the lowest value among the grid center points of all grid cells and the elevation values of the center point digital elevation model data, create the steepest slope line within the sliding boundary range; pass through the grid center point of each grid cell, create a normal vector plane perpendicular to the steepest slope line within the sliding boundary range, and obtain the intersection coordinates of the normal vector plane and the landslide sliding boundary range and the elevation value of the digital elevation model data of the intersection point;

[0027] A prediction module, configured to calculate the distance between the intersection point 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; based on the predicted depth of the landslide sliding surface at each grid point, determine the sliding volume of the single landslide.

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

[0029] The present invention also provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, the above method for predicting the sliding volume of a single landslide is implemented.

[0030] At least one of the above technical solutions adopted by the present invention can achieve the following beneficial effects:

[0031] A method for predicting the sliding volume of a single landslide provided by the present invention creates the steepest slope line within the sliding boundary range by connecting the two grid center points corresponding to the highest value and the lowest value of the digital elevation model data within the landslide sliding boundary; creates a normal vector plane perpendicular to the steepest slope line within the sliding boundary range through the grid center points of each grid cell in the digital elevation model data; calculates the intersection coordinates of the normal vector plane and the landslide sliding boundary range and the elevation value of the digital elevation model data at the intersection point; calculates the distance between the intersection point of the normal vector plane and the landslide sliding boundary range and the steepest slope line within the sliding boundary range and performs cubic spline interpolation on it to predict the landslide sliding surface depth of each grid point; predicts the landslide sliding volume based on the predicted landslide sliding surface depth of each grid point.

[0032] This method only requires 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 failure and influence range of landslide sliding, identify potential landslide hazard areas, and provide scientific support for landslide disaster prevention and control, personnel disaster prevention and avoidance strategy planning and management. Brief Description of the Drawings

[0033] 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 to the present application. In the drawings:

[0034] Figure 1 It is a schematic flow chart of a method for predicting the landslide sliding surface and sliding volume of a single landslide provided by the present invention;

[0035] Figure 2 It is the digital elevation model (DEM) and sliding boundary range of the Qingliu landslide for a method for predicting the landslide sliding surface and sliding volume of a single landslide provided by the present invention;

[0036] Figure 3 It is the prediction result of the landslide sliding surface depth of the Qingliu landslide for a method for predicting the landslide sliding surface and sliding volume of a single landslide provided by the present invention;

[0037] Figure 4 It is the prediction result of the landslide sliding volume of the Qingliu landslide for a method for predicting the landslide sliding surface and sliding volume of a single landslide provided by the present invention;

[0038] Figure 5 It is a schematic diagram of a device for predicting the sliding volume of a single landslide provided by the present invention;

[0039] Figure 6 It is 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. Specific Embodiments

[0040] To make the objectives, technical solutions, and advantages of the present invention clearer, the technical solutions of this application will be clearly and completely described below in conjunction with specific embodiments of the present invention and the corresponding drawings. Obviously, the described embodiments are only a part of the embodiments of this application, rather than all of the embodiments. Based on the embodiments in the specification, all other embodiments obtained by those of ordinary skill in the art without creative efforts fall within the scope of protection of this application.

[0041] The following will describe in detail the technical solutions provided by each embodiment of this application in conjunction with the drawings.

[0042] Figure 1 It is a flow schematic diagram of a method for predicting the sliding volume of a single landslide in the present invention, specifically including the following steps:

[0043] S101: 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.

[0044] Among them, the data format of the digital elevation model data is tif; the data format of the sliding boundary range data is shp.

[0045] Both the digital elevation model data and the sliding boundary range data use the WGS1984 spatial geographic coordinate system.

[0046] S102: Preprocess the digital elevation model data and the sliding boundary range data of the single landslide, project the preprocessed digital elevation model data and the sliding boundary range data onto the same coordinate system, and extract the grid center points of each grid cell in the coordinate system and the elevation values of the center point digital elevation model data.

[0047] Read the.shp file through the shpread function of the Matalb software and store the data content in a structure array to convert the file with the.shp format of the sliding boundary range into structure data.

[0048] Unify the coordinates of the digital elevation model data in tif format and the sliding boundary range data in shp format, and project them under the WGS1984 UTM coordinate system.

[0049] Extract the grid center points of each grid cell in the coordinate system and the elevation values of the center point digital elevation model data.

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

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

[0052] Specifically, input any two angle values within the range of 0° to 90° as the dip angle values of the left and right interpolation profiles.

[0053] Generate single-segment cubic spline curves at the center of each grid point respectively, and the lowest elevation value of the digital elevation model data on each interpolation curve is the sliding surface depth of each grid center point.

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

[0055] The formula for predicting the landslide sliding volume is:

[0056] ;

[0057] In the formula, are the resolutions of the input DEM in the x, y directions respectively; i, j are the i th row and the j th grid center point respectively, is the DEM elevation value at the i th row and the j th grid center point, is the predicted value of the sliding surface depth at the i th row and the j th grid center point.

[0058] In the specific embodiment of the present invention, a certain Qingliu landslide in a province is selected to predict the sliding surface and sliding volume of a single landslide. On the basis of the Figure 1 embodiment, a method for predicting the sliding volume of a single landslide proposed by the present invention is described in detail.

[0059] See 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.

[0060] 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.

[0061] 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.

[0062] 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.

[0063] 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.

[0064] 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.

[0065] Then the sliding volume of Qingliu landslide is predicted, and the calculation formula is as follows:

[0066] ;

[0067] In the above formula, are respectively the resolutions of the input DEM in the x, y directions; i, j are respectively the i th row and the j th grid center point, is the DEM elevation value at the grid center point of the i th row and the j th grid center point, is the predicted value of the slip surface depth at the grid center point of the i th row and the j th grid center point.

[0068] In the specific embodiment of the present invention, the resolutions of the input DEM in the x, y directions 、 are 30 m and 30 m respectively. The number of rows of the generated grid center points i is 31 rows in total, and the number of columns j is 48 rows in total. The predicted slip volumes of the 1488 grid center points are accumulated. Refer to Figure 4 , which is the predicted result of the slip volume of the Qingliu landslide. The predicted total slip volume is 24.8×10 4 m 3 .

[0069] For the specific limitations on a device for predicting the slip volume of a single landslide, reference can be made to the limitations on the method for predicting the slip volume of a single landslide in the above text, which will not be elaborated here. Each module in the above device for predicting the slip volume of a single landslide can be implemented in whole or in part by software, hardware, and their combination. The above modules can be embedded in the processor of the computer device in hardware form or be independent of it, or be stored in the memory of the computer device in software form, so that the processor can call and execute the operations corresponding to the above modules.

[0070] Figure 5 FIG. is a schematic diagram of a device for predicting the slip volume of a single landslide provided by the present invention, including:

[0071] A data acquisition module 601, configured to obtain digital elevation model data of a single landslide, input it into the slip boundary range of the single landslide, and obtain slip boundary range data;

[0072] A data extraction module 602, configured to preprocess the digital elevation model data and the slip boundary range data of the single landslide, project the preprocessed digital elevation model data and the slip boundary range data onto the same coordinate system, and extract the grid center points of each grid cell in the digital elevation model data in the coordinate system and the elevation values of the center point digital elevation model data;

[0073] The landslide steepest slope line and normal vector plane construction module 603 is used to connect the grid center points of all raster cells and the two grid center points corresponding to the highest and lowest elevation values in the elevation values of the center point digital elevation model data, so as to create the steepest slope line within the sliding boundary; passing through the grid center point of each raster cell, create a normal vector plane perpendicular to the steepest slope line within the sliding boundary, and obtain the intersection coordinates of the normal vector plane and the landslide sliding boundary range and the elevation value of the digital elevation model data of the intersection points;

[0074] The prediction module 604 is used to calculate the distance between the intersection point of the normal vector plane and the landslide sliding boundary range and the steepest slope line within the sliding boundary and perform cubic spline interpolation on it, so as to predict the landslide sliding surface depth of each grid point; based on the predicted landslide sliding surface depth of each grid point, determine the landslide sliding volume of the single landslide.

[0075] For the specific limitations of a device for predicting the sliding volume of a single landslide, reference can be made to the limitations of a method for predicting the sliding volume of a single landslide in the above text, which will not be elaborated here. Each module in the above device for predicting the sliding volume of a single landslide can be implemented in whole or in part by software, hardware and their combination. The above modules can be embedded in the processor in the computer device in hardware form or independent of it, or stored in the memory in the computer device in software form, so that the processor can call and execute the operations corresponding to the above modules.

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

[0077] The present invention also provides Figure 6 the structural schematic diagram of the computer device shown, as Figure 6 shown, at the hardware level, the computer device includes a processor, an internal bus, a network interface, a memory, and a non-volatile memory. Of course, other hardware required for other services may also be included. The processor reads the corresponding computer program from the non-volatile memory into the memory and then runs it to implement the above Figure 1 provided method for predicting the sliding volume of a single landslide.

[0078] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. 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 methods. Among them, any reference to a memory, storage, database, or other medium used in the various embodiments provided in the present application can include at least one of non-volatile and volatile memories. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, or optical memory, etc. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc.

[0079] The technical features of the above embodiments can be combined arbitrarily. For the sake of concise description, 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, it should be considered as the scope recorded in the present invention.

[0080] The technical features of the above embodiments can be combined arbitrarily. For the sake of concise description, 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, it should be considered as the scope recorded in the present invention.

Claims

1. A method for predicting the sliding volume of a single landslide, characterized in that, Including: 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; Preprocess the digital elevation model data and the sliding boundary range data of the single landslide, project the preprocessed digital elevation model data and the sliding boundary range data onto the same coordinate system, and extract the grid center point of each grid cell in the coordinate system and the elevation value of the center point digital elevation model data; Connect the two grid center points corresponding to the highest value and the lowest value among the grid center points of all grid cells and the elevation values of the center point digital elevation model data to create the steepest slope line within the sliding boundary range; Pass through the grid center point of each grid cell to create a normal vector plane perpendicular to the steepest slope line within the sliding boundary range, and obtain the intersection coordinates of the normal vector plane and the landslide sliding boundary range and the elevation value of the digital elevation model data at the intersection; 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 landslide sliding surface depth of each grid point; Based on the predicted landslide sliding surface depth of each grid point, determine the landslide sliding volume of the single landslide; The prediction of the landslide sliding surface depth of each grid point specifically includes: Using any two angle values within the angle range of 0° to 90° as the inclination angle values of the left and right interpolation profiles; Generating single-segment cubic spline interpolation curves at the center of each grid point respectively, and the lowest elevation value of the digital elevation model data on each interpolation curve is the sliding surface depth of each grid center point.

2. The method for predicting the sliding volume of a single landslide according to claim 1, characterized in that, The extraction of the grid center point of each grid cell in the digital elevation model data in the coordinate system and the elevation value of the center point digital elevation model data specifically includes: Read the.shp file through the shpread function of the Matalb software and store the data content in a structure array to convert the file with the sliding boundary range in.shp format into structure data; Unify the coordinates of the digital elevation model data in.tif format and the sliding boundary range data in.shp format and project them under the WGS1984 UTM coordinate system; Extract the grid center point of each grid cell in the coordinate system and the elevation value of the center point digital elevation model data.

3. The method for predicting the sliding volume of a single landslide according to claim 1, characterized in that, The determination of the sliding volume of the single landslide based on the predicted landslide sliding surface depth of each grid point, the formula is: V = ΔX × ΔY × ∑ i,j (Z DEM (i, j) - Z H (i, j)); Where, ΔX and ΔY are the resolutions of the input DEM in the x and y directions respectively; i and j are the central points of the i-th row and the j-th grid respectively, and Z DEM (i, j) is the DEM elevation value at the central point of the i-th row and the j-th grid, and Z H (i, j) is the predicted value of the slip surface depth at the central point of the i-th row and the j-th grid.

4. A device for predicting the sliding volume of a single landslide, including: A data acquisition module, configured 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, configured to preprocess the digital elevation model data and the sliding boundary range data of the single landslide, project the preprocessed digital elevation model data and the sliding boundary range data onto the same coordinate system, and extract the grid center point of each grid cell in the digital elevation model data in the coordinate system and the elevation value of the center point digital elevation model data; The steepest slope line and normal vector plane construction module of the landslide is used to connect the grid center points of all raster pixels and the two grid center points corresponding to the highest and lowest elevation values in the elevation values of the center point digital elevation model data, and create the steepest slope line within the sliding boundary; passing through the grid center points of each raster pixel, create a normal vector plane perpendicular to the steepest slope line within the sliding boundary, and obtain the intersection coordinates of the normal vector plane and the landslide sliding boundary range and the elevation value of the digital elevation model data at the intersection points; The prediction module is used to calculate the distance between the intersection points of the normal vector plane and the landslide sliding boundary range and the steepest slope line within the sliding boundary and perform cubic spline interpolation on it to predict the landslide sliding surface depth of each grid point; based on the predicted landslide sliding surface depth of each grid point, determine the landslide sliding volume of the single landslide; the prediction of the landslide sliding surface depth of each grid point specifically includes: using any two angle values within the angle range of 0° to 90° as the inclination angle values of the left and right interpolation profiles; generating single-segment cubic spline interpolation curves at the center of each grid point respectively, and the lowest elevation value of the digital elevation model data on each interpolation curve is the sliding surface depth of each grid center point.

5. 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 described in any one of claims 1 to 3 above is implemented.

6. A computer device, characterized in that, It includes a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, the method described in any one of claims 1 to 3 above is implemented.

Images