CN114662950A - Slope unit thinning method suitable for regional earthquake landslide risk assessment - Google Patents
Slope unit thinning method suitable for regional earthquake landslide risk assessment Download PDFInfo
- Publication number
- CN114662950A CN114662950A CN202210328414.1A CN202210328414A CN114662950A CN 114662950 A CN114662950 A CN 114662950A CN 202210328414 A CN202210328414 A CN 202210328414A CN 114662950 A CN114662950 A CN 114662950A
- Authority
- CN
- China
- Prior art keywords
- slope
- data
- positive
- negative
- unit
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000012502 risk assessment Methods 0.000 title claims abstract description 22
- 239000011435 rock Substances 0.000 claims abstract description 53
- 238000004458 analytical method Methods 0.000 claims abstract description 49
- 230000001133 acceleration Effects 0.000 claims abstract description 38
- 238000011160 research Methods 0.000 claims abstract description 7
- 230000011218 segmentation Effects 0.000 claims description 18
- 238000007670 refining Methods 0.000 claims description 16
- 230000005484 gravity Effects 0.000 claims description 8
- 230000003068 static effect Effects 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000004927 fusion Effects 0.000 claims description 3
- 239000003673 groundwater Substances 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 238000013507 mapping Methods 0.000 abstract description 6
- 230000007246 mechanism Effects 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 238000000605 extraction Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 230000001488 breeding effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007477 logistic regression Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Economics (AREA)
- Operations Research (AREA)
- Strategic Management (AREA)
- Pure & Applied Mathematics (AREA)
- Geometry (AREA)
- Mathematical Physics (AREA)
- Entrepreneurship & Innovation (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Development Economics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Educational Administration (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Game Theory and Decision Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marketing (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
The invention discloses a slope unit thinning method suitable for regional earthquake landslide risk assessment, which comprises the following steps: s1, collecting and sorting topographic data and geological data of the research area, and acquiring key rock mass mechanical parameters based on the topographic data and the geological data; s2, extracting a slope unit by adopting a traditional hydrological analysis method based on the topographic data; s3, calculating the slope yield acceleration based on the grid unit according to the earthquake slope stability analysis principle; s4, the slope units extracted by the traditional hydrological analysis method are used as limiting layers, and the traditional slope units are refined based on the slope yield acceleration difference obtained in S3, so that the refined slope units are obtained. The slope unit refined by the method disclosed by the invention not only retains the boundary information of the traditional slope unit, but also integrates a related physical mechanism formed by the earthquake landslide, can be effectively used for a mapping unit for regional earthquake landslide risk assessment, and is beneficial to improving the reliability of an assessment result.
Description
Technical Field
The invention relates to the technical field of regional earthquake geological disaster risk assessment, in particular to a slope unit thinning method suitable for regional earthquake landslide risk assessment.
Background
Landslide is a sliding geological phenomenon which occurs along a weak sliding surface of a slope rock-soil body, and a large amount of casualties and property loss are often caused. The earthquake landslide is an earthquake secondary disaster phenomenon directly triggered by an earthquake, has the characteristics of strong burst property, wide distribution range, large quantity, long duration and the like, and is a main geological disaster type in the southwest mountain areas of China. In order to prevent and alleviate the influence of earthquake landslide disasters, carry out regional earthquake landslide risk assessment, accurately divide the spatial distribution of potential earthquake landslide hidden dangers, and take prevention and control measures in time, the method is the most effective way at present.
The selection of the drawing unit is an important precondition for carrying out regional earthquake landslide risk assessment. Grid elements and ramp elements are the most commonly used mapping elements at present. The grid unit is a regular mapping unit, is simple to operate, and is widely used as a mapping unit for earthquake landslide risk assessment. However, for regional earthquake landslide risk assessment, the data processing amount is large, the operation time is long, the requirements on software and hardware are high, and the rapid emergency assessment of a large region is difficult to realize. In addition, the regular grid cells correspond to independent landform cells, which do not reflect the geometric features of the regional ground. In contrast, the slope unit is a set of grid units, which has abundant geometric shape information and can more effectively reflect the surface shape characteristics of the area. Meanwhile, landslide is used as a geological disaster developing on the slope unit, earthquake landslide risk assessment is carried out based on the slope unit, the assessment result has indication significance, and the landslide risk assessment method is more suitable for being used as an important reference for subsequent disaster reduction and prevention work.
At present, various methods exist for extracting the slope unit, including a forward and reverse DEM-based hydrological analysis extraction method, a forward and reverse terrain curvature hydrological analysis-based extraction method, and an r.slope units method proposed by alvioi. The extraction methods can quickly and effectively extract the slope units. However, it is worth pointing out that these extraction methods are based on the hydrological angle establishment, and the extracted slope units, although they can effectively correspond to the actual side slope, tend to have a larger area compared to the landslide. In addition, the traditional extraction method mainly uses the terrain as a dividing basis of the slope unit, and often neglects the difference of other geological environment factors, so that the geological environment attributes corresponding to the same slope unit do not have homogeneity. The spatial statistical analysis according to a large number of historical earthquake landslides shows that the breeding development of the landslide is closely related to regional geological environment attributes, such as engineering rock mass, terrain gradient and the like. Therefore, for landslide, the slope unit should be subdivided by integrating the disaster recovery mechanism of landslide. Obviously, the existing ramp unit extraction method does not adopt a corresponding refinement strategy for a specific application scenario. In summary, aiming at the defects of the existing slope unit extraction method, the invention provides a slope unit refining method suitable for regional earthquake landslide risk assessment based on earthquake slope stability analysis by comprehensively considering regional geological environment information, and the reliability of earthquake landslide risk assessment can be effectively improved.
Disclosure of Invention
The invention aims to provide a slope unit refining method suitable for regional earthquake landslide risk assessment, and aims to solve the problems in the background art.
In order to achieve the purpose, the invention provides a slope unit refining method suitable for regional earthquake landslide risk assessment, which comprises the following steps of:
s1, collecting and sorting topographic data and geological data of the research area, and acquiring key rock mass mechanical parameters based on the topographic data and the geological data; the topographic data mainly comprises a Digital Elevation Model (DEM) and topographic gradient information calculated based on the DEM, and the geological data mainly comprises engineering geological maps with various scales; the key rock mass mechanics parameters mainly comprise effective cohesion, effective internal friction angle and severe parameters corresponding to different rock mass types in geological data.
S2, extracting slope units by adopting a traditional hydrological analysis method based on the topographic data acquired in the step S1;
s3, calculating slope yield acceleration based on the grid unit according to the terrain data and the key rock mass mechanical parameters acquired in the step S1 and the earthquake slope stability analysis principle;
s4, taking the slope unit extracted in the step S2 as a limiting layer, and refining the traditional slope unit based on the difference of the yield acceleration of the slope acquired in the step S3, so as to acquire the refined slope unit.
Further, the specific steps of obtaining the key rock mass mechanical parameters in step S1 are as follows:
s1.1, extracting terrain gradient information based on the terrain data;
s1.2, counting the mean value mu, the standard deviation sigma and the maximum value M of the gradient in different rock mass ranges based on rock mass information in the geological data;
s1.3, respectively calculating critical slope values a corresponding to different rock masses according to the mean value mu, the standard deviation sigma and the maximum value M of the slopes in different rock mass ranges counted in the step S1.2, wherein when a is equal to mu +3 sigma and is less than or equal to M, a is equal to mu +3 sigma; when a ═ μ +3 σ > M, a ═ M;
s1.4, endowing all rock masses with uniform internal friction angles and weights of the rock masses, and respectively calculating the minimum cohesive force of different rock masses which are damaged under the condition of the critical slope values according to the Moore-coulomb damage criterion and the critical slope values of different rock masses obtained in the step S1.3;
s1.5, determining the internal friction angle of the rock mass, the weight of the rock mass and the minimum cohesion calculated based on the critical gradient value a as key rock mass mechanics parameters.
Further, the topographic data utilized in step S2 is mainly DEM; in step S2, a slope unit is extracted by using a conventional hydrological analysis method, the main purpose is to divide an initial slope unit, a hydrological analysis extraction method based on a positive-negative DEM is used, and the concrete steps of extracting the slope unit are as follows:
s2.1, positive and negative DEM data acquisition and terrain depression filling: the positive DEM data represents original DEM data, the reverse DEM data is obtained by differentiating an elevation maximum value in the DEM data and the original DEM data, and terrain filling is carried out on the DEM data by using a filling tool in the GIS hydrological analysis module, so that the influence of terrain cavities is eliminated;
s2.2, acquiring positive and negative flow data: based on the positive and negative DEM data after the depression is filled, respectively extracting the positive and negative flow direction data by adopting a flow direction analysis tool in a GIS hydrological analysis module;
s2.3, acquiring positive and negative accumulated flow data: acquiring positive and negative accumulated flow data by adopting an accumulated flow statistical tool in a GIS hydrological analysis module based on the positive and negative flow data;
s2.4, network division of positive and negative river flows: determining a proper flow threshold value according to the positive and negative DEM accumulated flow data and the water system distribution on the remote sensing image, carrying out binarization on the positive and negative accumulated flow data by adopting grid calculation in a GIS hydrological analysis module, and acquiring positive and negative river vector data by using a grid line-turning tool;
s2.5, extracting the positive and negative catchment watersheds: according to positive and negative river vector data, based on a basin analysis tool in the GIS hydrological analysis module, extracting positive and negative catchment basins: wherein, the boundary of the positive catchment basin corresponds to a ridge line, and the boundary of the reverse catchment basin corresponds to a valley line;
and S2.6, carrying out vector fusion on the positive and negative catchment watersheds to obtain an initial slope unit consisting of a valley and a ridge line.
Further, the concrete steps of calculating the yield acceleration of the side slope in the step S3 are as follows;
s3.1, calculating slope static safety factor F based on slope stability analysis methods：
In formula 1), c' and
effective cohesion and effective internal friction angle of the rock mass, gamma and gamma, respectivelywRespectively representing the rock mass gravity and the groundwater gravity, and m represents the landslide saturation of the slope, namely the landslideThe ratio of the middle saturated part to the thickness of the sliding body, t is the potential thickness of the sliding body of the side slope, and theta is the inclination angle of the sliding surface, namely the slope of the regional terrain;
s3.2, according to the static safety factor FsThe method adopts a quasi-static analysis method to analyze the stability of the earthquake slope and calculate the yield acceleration a of the slopec：
ac＝(Fs-1)*gsinθ2)；
In the formula 2), g is the gravity acceleration, and θ is the inclination angle of the sliding surface, namely the slope of the regional terrain.
Further, the step S4 of refining the slope unit specifically includes:
s4.1, taking the slope unit extracted in the step S2 as a limiting layer, and performing initial segmentation on the slope yield acceleration obtained in the step S3;
and S4.2, merging the grid units in the segmentation object based on the yield acceleration information of the grid units in the segmentation object according to the obtained initial segmentation result.
Further, the step S4.2 of merging the grid cells inside the segmentation object specifically includes:
s4.2.1, taking the grid cell corresponding to the boundary of the segmentation object as an initial cell, and calculating the yield acceleration difference value of the initial cell and the neighborhood cell;
s4.2.2, comparing the calculated yield acceleration difference with a set threshold, and if the yield acceleration difference is smaller than the set threshold, combining the yield acceleration difference into a unit; otherwise, no merging is performed.
Compared with the prior art, the invention has the following beneficial effects:
the invention relates to a slope unit refining method suitable for regional earthquake landslide risk assessment, which comprises the steps of firstly utilizing regional topographic data and extracting slope units by adopting a traditional hydrological analysis method; then, integrating regional topographic data and geological data, and calculating yield acceleration according to the earthquake slope stability analysis principle; and finally, taking the slope units extracted by the traditional hydrological analysis method as limiting layers, and combining the limiting layers based on the slope yield acceleration differences of the grid units in the slope units to obtain the refined slope units. The slope unit refined by the slope unit refining method not only retains the boundary information of the traditional slope unit, but also considers the relevant physical mechanism of landslide formation, so that the correlation between the refined slope unit and the landslide is tighter, the method can be effectively used for a mapping unit for regional earthquake landslide risk assessment, and is beneficial to improving the reliability of the assessment result.
In addition to the objects, features and advantages described above, other objects, features and advantages of the present invention are also provided. The present invention will be described in further detail below with reference to the drawings.
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 specification, illustrate embodiments of the invention and together with the description serve to explain the embodiments of the invention without limiting the embodiments of the invention. In the drawings:
FIG. 1 is a schematic flow chart of a slope unit refining method suitable for regional earthquake landslide risk assessment according to the invention;
FIG. 2(a) is a slope unit extracted by the forward and reverse DEM-based hydrological analysis method in the embodiment of the present invention;
FIG. 2(b) is an enlarged schematic view of A in FIG. 2 (a);
FIG. 2(c) is an enlarged schematic view of B in FIG. 2 (a);
FIG. 3(a) is the result of the initial segmentation of the yield acceleration based on the initial ramp unit in an embodiment of the present invention;
FIG. 3(b) is an enlarged schematic view of the structure at C in FIG. 3 (a);
FIG. 3(c) is an enlarged schematic view of D in FIG. 3 (a);
FIG. 4(a) is a schematic diagram comparing a refined ramp unit with a conventional ramp unit based on the method of the present invention;
FIG. 4(b) is a schematic diagram of a partial enlarged structure based on conventional ramp unit extraction;
FIG. 4(c) is a schematic diagram of a partial enlarged structure (labeled in FIG. 4 (a)) of the slope unit extraction after the thinning of the present invention in FIG. 4 (a);
FIG. 5(a) is a regional seismic landslide hazard map for a study area based on refined ramp units in an embodiment of the invention;
FIG. 5(b) is a schematic diagram of a partial enlargement of the regional landslide hazard mapping of FIG. 5(a) based on the refined ramp units of the present invention (marked in FIG. 5 (a));
fig. 5(c) is a schematic view of a partially enlarged structure (at the same position as the partially enlarged structure in fig. 5 (b)) based on a regional landslide risk map of a conventional slope unit;
fig. 5(d) is a schematic view of a partially enlarged structure (at the same position as the partially enlarged structure in fig. 5 (b)) of the grid cell-based regional landslide risk map.
Detailed Description
Embodiments of the invention will be described in detail below with reference to the drawings, but the invention can be implemented in many different ways, which are defined and covered by the claims.
The study area of one embodiment of the present invention was located in southeast of the Atha Qiang autonomous State of Kagawa, Sichuan province. Fig. 1 is a schematic flow diagram of a slope unit refining method suitable for regional earthquake landslide risk assessment according to an embodiment of the present invention, where the slope unit refining method specifically includes the following steps:
and S1, collecting and sorting the topographic data and the geological data of the research area, and acquiring key rock mass mechanical parameters based on the topographic data and the geological data. Specifically, the key rock mass mechanical parameters are obtained by the following steps:
and S1.1, extracting terrain gradient information based on the terrain data.
S1.2, counting the mean value mu, the standard deviation sigma and the maximum value M of the gradient in different rock mass ranges based on rock mass information in geological data.
S1.3, respectively calculating critical slope values a corresponding to different rock masses according to the mean value mu, the standard deviation sigma and the maximum value M of the slopes in different rock mass ranges obtained in the step S1.2, wherein when a is equal to mu +3 sigma and is less than or equal to M, a is equal to mu +3 sigma; when a ═ μ +3 σ > M, a ═ M.
And S1.4, endowing all rock masses with uniform internal friction angles and weights of the rock masses, and respectively calculating the minimum cohesive force of different rock masses which are damaged under the condition of the critical slope values according to the molar-coulomb damage criterion and the critical slope values of different rock masses obtained in the step S1.3. In one embodiment of the invention, the internal friction angles corresponding to all rock masses are set to be 30 degrees, and the gravity is 22kN/m3。
S1.5, determining the internal friction angle of the rock mass, the weight of the rock mass and the minimum cohesion calculated based on the critical gradient value a as key rock mass mechanical parameters.
In step S1, the terrain data is digital elevation models with different resolutions or digital terrain maps, and the geological data corresponds to engineering geological maps with different scales. In the embodiment of the invention, the topographic data adopts SRTM DEM with 90m resolution, and the geological data adopts the following formula 1: A50W engineering geological map, a research area is one of Wenchuan earthquake severe disaster areas.
S2, extracting the slope unit by using the conventional hydrological analysis method based on the topographic data acquired in step S1, i.e. extracting the conventional slope unit. In one embodiment of the invention, the slope unit is extracted by adopting a forward and reverse DEM-based hydrological analysis method according to DEM data of a research area. The method comprises the following specific steps:
s2.1, positive and negative DEM data acquisition and terrain depression filling: the positive DEM data represent original DEM data, the reverse DEM data are obtained by differentiating elevation maximum values in the DEM data and the original DEM data, and terrain filling is carried out on the DEM data by using a filling tool in the GIS hydrological analysis module, so that the influence of terrain cavities is eliminated.
S2.2, acquiring positive and negative flow data: and based on the positive and negative DEM data after the hollow filling, respectively extracting the positive and negative flow data by adopting a flow direction analysis tool in the GIS hydrological analysis module.
S2.3, acquiring positive and negative accumulated flow data: and acquiring the positive and negative accumulated flow data by adopting an accumulated flow statistical tool in the GIS hydrological analysis module based on the positive and negative flow data.
S2.4, network division of positive and negative river flows: and determining a proper flow threshold value according to the positive and negative DEM accumulated flow data and the water system distribution on the remote sensing image, carrying out binarization on the positive and negative accumulated flow data by adopting grid calculation in a GIS hydrological analysis module, and acquiring positive and negative river vector data by using a grid line-turning tool. In one embodiment of the invention, the flow threshold is set at 2500.
S2.5, extracting the positive and negative catchment watersheds: according to positive and negative river vector data, based on a basin analysis tool in the GIS hydrological analysis module, extracting positive and negative catchment basins: wherein, the positive catchment basin boundary is corresponding to the crest line, and the negative catchment basin boundary is corresponding to the valley line.
And S2.6, carrying out vector fusion on the positive and negative catchment watersheds to obtain an initial slope unit consisting of a valley and a ridge line.
The slope units extracted by the hydrological analysis method based on the positive DEM and the negative DEM in step S2 are shown in fig. 2(a), 2(b), and 2 (c).
S3, calculating the slope yield acceleration based on the grid unit according to the terrain data and the key rock mass mechanical parameters acquired in the step S1 and the earthquake slope stability analysis principle; the method comprises the following specific steps:
s3.1, calculating slope static safety factor F based on slope stability analysis methods：
In formula 1), c' and
effective cohesion and effective internal friction angle of the rock mass, gamma and gamma, respectivelywRespectively representing the rock mass weight and the groundwater weight, m representing the landslide saturation of the side slope, namely the ratio of the saturated part in the landslide to the thickness of the landslide, t representing the potential thickness of the side slope, and theta representing the inclination angle of the sliding surface, namely the slope of the regional terrain. The slip saturation m is generally a global uniform empirical value, and in the embodiment of the invention, the slip saturation m is set to be 0, namely the influence of underground water is not considered; the thickness t of the slider is set to 2.5 m.
S3.2 according to the static safety systemNumber FsAnalyzing the stability of the earthquake slope by adopting a quasi-static analysis method, and calculating the yield acceleration a of the slopec：
ac＝(Fs-1)*gsinθ2)；
In the formula 2), g is the gravity acceleration, and θ is the inclination angle of the sliding surface, namely the slope of the regional terrain.
And S4, taking the traditional slope unit extracted in the step S2 as a limiting image layer, and refining the traditional slope unit based on the difference of the yield acceleration of the slope acquired in the step S3, so as to acquire the refined slope unit.
The basic idea of the slope unit refining method provided by the embodiment of the invention is derived from the application scene, namely the related slope unit is mainly applied to landslide risk evaluation, especially regional earthquake landslide risk evaluation. Therefore, physical mechanism analysis related to earthquake landslide is integrated in the process of thinning the slope units, and the closer connection between the thinned slope units and the earthquake landslide is ensured. The method comprises the following two steps:
and S4.1, taking the slope unit extracted in the step S2 as a limiting layer, and initially dividing the slope yield acceleration obtained in the step S3.
And S4.2, merging the grid cells in the segmentation object based on the yield acceleration information of the grid cells in the segmentation object according to the obtained initial segmentation result. Further, the specific steps of merging the grid cells inside the segmentation object are as follows:
s4.2.1, taking the grid cell corresponding to the boundary of the segmentation object as an initial cell, calculating the yield acceleration difference value of the initial cell and the cell in the neighborhood of the initial cell.
S4.2.2, comparing the calculated yield acceleration difference with a set threshold, and combining the calculated yield acceleration difference with a unit if the yield acceleration difference is smaller than the set threshold; otherwise, no combination is performed. The smaller the threshold setting, the finer the ramp cell refinement result. In one embodiment of the present invention, it sets the threshold value to 0.2.
The result of the initial division of the yield acceleration based on the initial slope unit using step S4 is shown in fig. 3(a), 3(b), and 3 (c).
Referring to fig. 4(a), 4(b) and 4(c), it is apparent that the slope unit after being refined is subdivided on the basis of the slope unit extracted by the conventional method, and the corresponding slope unit is more refined. Meanwhile, it is worth to be noted that the slope unit extracted by the conventional method is difficult to effectively distinguish flat ground from a slope, and a flat ground area is usually directly merged with an adjacent slope unit, so that a terrain abrupt change phenomenon exists in the conventional slope unit. The slope unit after thinning can effectively distinguish flat ground and side slope, and the corresponding topographic features are more homogeneous.
In addition, in order to verify the reliability of the refined slope unit, the refined slope unit is applied to risk assessment of landslide of the Wenchuan earthquake, and the assessment results of the slope unit and the grid unit extracted based on the traditional method are compared. Fig. 5(a), 5(b), 5(c) and 5(d) are graphs of the risk of earthquake landslide evaluated by applying a logistic regression method based on different graph units in the research area. As can be seen from the figure, for the conventional slope unit, although the conventional slope unit can better reflect the actual slope form, the same slope only has one earthquake landslide risk, and the evaluation result is easy to overestimate or underestimate the actual earthquake landslide risk; for the grid unit, the evaluation result is more discrete, and the phenomenon of overestimating or underestimating the earthquake landslide risk of the area also exists; for the slope unit after the thinning, the earthquake landslide risk evaluated by the slope unit is more focused, and the earthquake landslide risk partition result is more consistent with the actual earthquake landslide distribution.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (6)
1. A slope unit thinning method suitable for regional earthquake landslide risk assessment is characterized by comprising the following steps:
s1, collecting and sorting topographic data and geological data of the research area, and acquiring key rock mass mechanical parameters based on the topographic data and the geological data;
s2, extracting slope units by adopting a hydrological analysis method based on the topographic data acquired in the step S1;
s3, calculating slope yield acceleration based on grid units according to the terrain data and the key rock mass mechanical parameters obtained in the step S1 and the earthquake slope stability analysis principle;
s4, taking the slope unit extracted in the step S2 as a limiting layer, and thinning the slope unit based on the difference of the yield acceleration of the slope acquired in the step S3, so as to acquire a thinned slope unit.
2. The slope unit refining method according to claim 1, wherein the concrete steps of obtaining key rock mechanical parameters in the step S1 are as follows:
s1.1, extracting terrain gradient information based on the terrain data;
s1.2, counting the mean value mu, the standard deviation sigma and the maximum value M of the slopes in different rock mass ranges based on rock mass information in the geological data;
s1.3, respectively calculating critical slope values a corresponding to different rock masses according to the mean value mu, the standard deviation sigma and the maximum value M of the slopes in different rock mass ranges counted in the step S1.2, wherein when a is equal to mu +3 sigma and is less than or equal to M, a is equal to mu +3 sigma; when a ═ μ +3 σ > M, a ═ M;
s1.4, endowing all rock masses with uniform internal friction angles and weights of the rock masses, and respectively calculating the minimum cohesive force of different rock masses which are damaged under the condition of the critical slope values according to the Moore-coulomb damage criterion and the critical slope values of different rock masses obtained in the step S1.3;
s1.5, determining the internal friction angle of the rock mass, the weight of the rock mass and the minimum cohesion calculated based on the critical gradient value a as key rock mass mechanics parameters.
3. The slope unit refining method according to claim 1, wherein the concrete steps of extracting the slope unit by using the hydrologic analysis method in the step S2 are as follows:
s2.1, positive and negative DEM data acquisition and terrain depression filling: the positive DEM data represents original DEM data, the reverse DEM data is obtained by differentiating elevation maximum values in the DEM data and the original DEM data, and terrain filling is carried out on the DEM data by using a filling tool in the GIS hydrological analysis module, so that the influence of terrain cavities is eliminated;
s2.2, acquiring positive and negative flow data: based on the positive and negative DEM data after the hollow filling, respectively extracting the positive and negative flow data by adopting a flow direction analysis tool in a GIS hydrological analysis module;
s2.3, acquiring positive and negative accumulated flow data: acquiring positive and negative accumulated flow data by adopting an accumulated flow statistical tool in a GIS hydrological analysis module based on the positive and negative flow data;
s2.4, network division of positive and negative river flows: determining a proper flow threshold value according to the positive and negative DEM accumulated flow data and the water system distribution on the remote sensing image, carrying out binaryzation on the positive and negative accumulated flow data by adopting the grid calculation in the GIS hydrological analysis module, and obtaining the vector data of the positive and negative rivers by using a grid line-turning tool;
s2.5, extracting the positive and negative catchment watersheds: according to positive and negative river vector data, based on a basin analysis tool in the GIS hydrological analysis module, extracting positive and negative catchment basins: wherein, the boundary of the positive catchment basin is corresponding to a ridge line, and the boundary of the reverse catchment basin is corresponding to a valley line;
and S2.6, carrying out vector fusion on the positive and negative catchment watersheds to obtain an initial slope unit consisting of a valley and a ridge line.
4. The slope unit refinement method according to claim 1, wherein the concrete steps of calculating the yield acceleration of the slope in the step S3 are as follows:
s3.1, calculating slope static safety factor F based on slope stability analysis methods：
In formula 1), c' and
effective cohesion and effective internal friction angle of the rock mass, gamma and gamma, respectivelywRespectively representing the rock mass gravity and the groundwater gravity, m represents the landslide saturation of the side slope, namely the ratio of the saturated part in the landslide to the thickness of the landslide, t is the potential thickness of the side slope, and theta is the dip angle of the sliding surface, namely the regional terrain gradient;
s3.2, according to the static safety factor FsThe method adopts a quasi-static analysis method to analyze the stability of the earthquake slope and calculate the yield acceleration a of the slopec：
ac＝(Fs-1)*gsinθ 2)；
In the formula 2), g is the gravity acceleration, and θ is the inclination angle of the sliding surface, namely the slope of the regional terrain.
5. The ramp unit refinement method according to claim 1, wherein the concrete steps of refining the ramp unit in step S4 are as follows:
s4.1, taking the slope unit extracted in the step S2 as a limiting layer, and performing initial segmentation on the slope yield acceleration obtained in the step S3;
and S4.2, merging the grid units in the segmentation object based on the yield acceleration information of the grid units in the segmentation object according to the obtained initial segmentation result.
6. The ramp element refinement method according to claim 5, characterized in that the concrete steps of merging the grid elements inside the segmentation object in the step S4.2 are as follows:
s4.2.1, taking the grid cell corresponding to the boundary of the segmentation object as an initial cell, and calculating the yield acceleration difference value of the initial cell and the neighborhood cell;
s4.2.2, comparing the calculated yield acceleration difference with a set threshold, and if the yield acceleration difference is smaller than the set threshold, combining the yield acceleration difference into a unit; otherwise, no merging is performed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210328414.1A CN114662950B (en) | 2022-03-31 | 2022-03-31 | Slope unit refining method suitable for regional earthquake landslide risk assessment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210328414.1A CN114662950B (en) | 2022-03-31 | 2022-03-31 | Slope unit refining method suitable for regional earthquake landslide risk assessment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114662950A true CN114662950A (en) | 2022-06-24 |
| CN114662950B CN114662950B (en) | 2024-06-14 |
Family
ID=82033869
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210328414.1A Active CN114662950B (en) | 2022-03-31 | 2022-03-31 | Slope unit refining method suitable for regional earthquake landslide risk assessment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114662950B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118397378A (en) * | 2024-06-27 | 2024-07-26 | 湖南科技大学 | Method, system, equipment and storage medium for classifying and extracting slope units |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103455711A (en) * | 2013-08-15 | 2013-12-18 | 广州地理研究所 | Small watershed region-oriented landslide hazard risk division method based on mechanism analysis |
| CN104240034A (en) * | 2014-09-23 | 2014-12-24 | 北京师范大学 | Method for quickly evaluating economic losses generated by landslides caused by earthquakes |
| CN107463991A (en) * | 2017-06-28 | 2017-12-12 | 西南石油大学 | A kind of Regional Landslide method for evaluating hazard based on slopes unit and machine learning |
| CN110390169A (en) * | 2019-07-25 | 2019-10-29 | 西南交通大学 | A kind of Seismic Landslide Hazard quantitative evaluation method based on mechanical model |
| KR102091758B1 (en) * | 2019-04-08 | 2020-03-20 | 주식회사 아이자랩 | Method for predicting collapse of slope land by using acceleration and acoustic emission wave |
| WO2021008282A1 (en) * | 2019-07-12 | 2021-01-21 | 清华大学 | Seismic landslide quick report analysis method and apparatus based on actually-measured seismic motion |
| CN112686522A (en) * | 2020-12-25 | 2021-04-20 | 中南大学 | Newmark correction model earthquake landslide risk assessment method |
| CN113268852A (en) * | 2021-04-14 | 2021-08-17 | 西南交通大学 | Monte Carlo simulation-based earthquake landslide probability risk analysis method |
| CN113850822A (en) * | 2021-09-18 | 2021-12-28 | 四川大学 | Automatic slope unit dividing method based on confluence segmentation |
-
2022
- 2022-03-31 CN CN202210328414.1A patent/CN114662950B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103455711A (en) * | 2013-08-15 | 2013-12-18 | 广州地理研究所 | Small watershed region-oriented landslide hazard risk division method based on mechanism analysis |
| CN104240034A (en) * | 2014-09-23 | 2014-12-24 | 北京师范大学 | Method for quickly evaluating economic losses generated by landslides caused by earthquakes |
| CN107463991A (en) * | 2017-06-28 | 2017-12-12 | 西南石油大学 | A kind of Regional Landslide method for evaluating hazard based on slopes unit and machine learning |
| KR102091758B1 (en) * | 2019-04-08 | 2020-03-20 | 주식회사 아이자랩 | Method for predicting collapse of slope land by using acceleration and acoustic emission wave |
| WO2021008282A1 (en) * | 2019-07-12 | 2021-01-21 | 清华大学 | Seismic landslide quick report analysis method and apparatus based on actually-measured seismic motion |
| CN110390169A (en) * | 2019-07-25 | 2019-10-29 | 西南交通大学 | A kind of Seismic Landslide Hazard quantitative evaluation method based on mechanical model |
| CN112686522A (en) * | 2020-12-25 | 2021-04-20 | 中南大学 | Newmark correction model earthquake landslide risk assessment method |
| CN113268852A (en) * | 2021-04-14 | 2021-08-17 | 西南交通大学 | Monte Carlo simulation-based earthquake landslide probability risk analysis method |
| CN113850822A (en) * | 2021-09-18 | 2021-12-28 | 四川大学 | Automatic slope unit dividing method based on confluence segmentation |
Non-Patent Citations (4)
| Title |
|---|
| 王佳佳;殷坤龙;杜娟;王轶力;: "基于GIS考虑准动态湿度指数的滑坡危险性预测水文力学耦合模型研究", 岩石力学与工程学报, no. 2, 15 July 2013 (2013-07-15) * |
| 王凯等: ""斜坡单元提取方法研究进展和展望"", 《长江科学学院院报》, vol. 37, no. 6, 30 June 2020 (2020-06-30), pages 86 - 90 * |
| 王涛;吴树仁;石菊松;辛鹏;: "基于简化Newmark位移模型的区域地震滑坡危险性快速评估――以汶川M_S8.0级地震为例", 工程地质学报, no. 01, 15 February 2013 (2013-02-15) * |
| 赵洪涛: "国道212线甘肃境内段沿线滑坡泥石流活动强度评价与区划研究", 兰州大学学报(自然科学版), no. 05, 28 October 2004 (2004-10-28) * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118397378A (en) * | 2024-06-27 | 2024-07-26 | 湖南科技大学 | Method, system, equipment and storage medium for classifying and extracting slope units |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114662950B (en) | 2024-06-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fan et al. | Earthquake‐induced chains of geologic hazards: Patterns, mechanisms, and impacts | |
| Cheng et al. | UAV photogrammetry-based remote sensing and preliminary assessment of the behavior of a landslide in Guizhou, China | |
| 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 | |
| Park et al. | Assessment of rainfall-induced shallow landslide susceptibility using a GIS-based probabilistic approach | |
| Malamud et al. | Landslide inventories and their statistical properties | |
| Lee et al. | Assessment of shallow landslide susceptibility using the transient infiltration flow model and GIS-based probabilistic approach | |
| Dong et al. | The role of the sediment budget in understanding debris flow susceptibility | |
| Robinson et al. | Coseismic landsliding estimates for an Alpine Fault earthquake and the consequences for erosion of the Southern Alps, New Zealand | |
| Li et al. | Geological survey and unstable rock block movement monitoring of a post-earthquake high rock slope using terrestrial laser scanning | |
| Weng et al. | Evaluating failure mechanisms of dip slope using a multiscale investigation and discrete element modelling | |
| CN111337997B (en) | Method for quickly identifying potential collapse disaster body | |
| Khal et al. | Evaluation of open Digital Elevation Models: estimation of topographic indices relevant to erosion risk in the Wadi M’Goun watershed, Morocco | |
| Mora et al. | Small landslide susceptibility and hazard assessment based on airborne lidar data | |
| Maharjan et al. | Rapid terrain assessment for earthquake-triggered landslide susceptibility with high-resolution DEM and critical acceleration | |
| CN114662950A (en) | Slope unit thinning method suitable for regional earthquake landslide risk assessment | |
| Jiang et al. | Debris flow assessment in the Gaizi-Bulunkou section of Karakoram highway | |
| Li et al. | Detection and analysis of potential landslides based on SBAS-InSAR technology in alpine canyon region | |
| Baltakova et al. | Analysis of debris flows by application of GIS and remote sensing: case study of western foothills of Pirin Mountains (Bulgaria) | |
| Baek et al. | A study on the use of planarity for quick identification of potential landslide hazard | |
| Pei et al. | Debris flow susceptibility mapping in Colorado Front Range, USA: A comparison of physics-based and data-driven approaches | |
| Amanzio et al. | Integration of terrestrial laser scanning and GIS analysis for multi-temporal landslide monitoring: A case study of the Mont de La Saxe (Aosta Valley, NW Italy) | |
| Imamura et al. | A practical method of nation wide Vs30 mapping with 250m grids for developing countries by automated topographical classification and published data in Japan | |
| CN118467892B (en) | Groove type avalanche debris flow movement distance and dangerous range prediction method | |
| Yang et al. | Characteristics and particle flow simulation of failure mechanism of Ludaping landslide in Renhuai red bed area, China | |
| OVERBECK et al. | Enhancing evaluation of post-storm morphologic response using aerial orthoimagery from Hurricane Sandy |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |