A reservoir landslide type rapid identification method and system coupled with time series InSAR technology

Abstract

The application provides a reservoir landslide type rapid identification method and system coupled with a timing InSAR technology, comprising: obtaining engineering geology condition data and a radar image set of a research area; obtaining four index categories for identifying the landslide type of the research area through the engineering geology condition data and the radar image set; obtaining four index data of different types of reservoir landslides in the research area through the engineering geology condition data, the radar image set and the four index categories; and constructing a discrimination model of different types of reservoir landslides through the four index data, and identifying the landslide type in the research area through the discrimination model. The discrimination model constructed through the index categories and the index data can correctly and rapidly identify different types of reservoir landslides, and then targeted early warning schemes are taken, so that the discrimination model has strong universality.

Description

Technical Field

[0001] This invention relates to the field of landslide disaster prediction and forecasting, and in particular to a method and system for rapid identification of reservoir landslide types coupled with time-series InSAR technology. Background Technology

[0002] my country is a major water conservancy country with abundant reservoirs. However, these reservoirs are often built in high mountains and deep valleys, with significant elevation differences and complex geological structures. This unique topographical feature provides potentially favorable conditions for the development of reservoir landslides. Surveys show that approximately 60% of landslides around reservoirs occur during periods of rapid decline in reservoir water levels, while the remaining 40% occur during periods of rising water levels. Based on the mechanism by which reservoir water levels affect landslides, they can be classified into hydrodynamic pressure landslides and buoyancy-reducing landslides.

[0003] Accurate and rapid identification of hydrodynamic pressure landslides and floating-off reduction landslides is of great significance for the rapid and accurate implementation of reasonable early warning plans, as well as for landslide disaster prevention, risk management, and the protection of people's lives and property.

[0004] However, the existing technology still has the following defects: (1) Traditional field landslide investigation methods and mass monitoring and prevention disaster prevention measures are insufficient for landslide monitoring and identification; (2) At present, there are no clear classification standards for different types of reservoir landslides, which is not conducive to the subsequent geological disaster prevention and control work. Summary of the Invention

[0005] To address the aforementioned technical problems, this invention provides a rapid identification method for reservoir landslide types coupled with time-series InSAR technology, comprising the following steps:

[0006] S1: Acquire engineering geological condition data and radar image sets for the study area;

[0007] S2: Four types of indicators are obtained through engineering geological data and radar image sets to identify landslide types in the study area, including: cumulative displacement curve, soil-rock ratio, sliding surface morphology and landslide leading edge shear exit elevation;

[0008] S3: By using engineering geological data, radar image sets, and four types of indicators, four types of indicators were obtained for different types of reservoir landslides in the study area, including: cumulative displacement curve shape, soil-rock ratio threshold, type of sliding surface morphology, and threshold of landslide leading edge shear outlet elevation.

[0009] S4: Construct identification models for different types of reservoir landslides using four indicator data, and identify landslide types in the study area using the identification models.

[0010] Preferably, the engineering geological data includes: the number and geographical location of existing landslides in the study area, the climate characteristics, geomorphological characteristics and stratigraphic lithology characteristics of the study area, and the rainfall data and reservoir water level scheduling data of the study area over the past ten years.

[0011] Preferably, the process of obtaining the cumulative displacement curve is as follows:

[0012] We acquired a set of Sentinel-1A radar images of the slopes in the study area over the past five years. We used a digital elevation model to remove the terrain phase and geocode the radar images to obtain interferograms. We used a small baseline set method to process the interferograms and calculated the surface deformation information in the study area with millimeter-level accuracy. We then used the surface deformation information to draw the cumulative displacement curve of the study area.

[0013] Preferably, the process for obtaining the soil-rock ratio is as follows:

[0014] By analyzing the lithological characteristics of the strata, the soil-to-rock ratio of the landslide deposits in the study area was calculated using the following formula:

[0015]

[0016]

[0017]

[0018] in, To determine the thickness of the soil layer sampled from the borehole, The value represents the thickness of the gravel layer sampled from the borehole, where i is the borehole number and n is the total number of boreholes. The ratio of soil to stone.

[0019] Preferably, the process of obtaining the smooth surface morphology is as follows:

[0020] Based on engineering geological data and borehole data, we analyzed and obtained information on the topographic features, lithological characteristics, rock component boundaries and geological structure of the study area, drew landslide profile maps of the study area, and obtained the sliding surface morphology based on the landslide profile maps.

[0021] Preferably, the process for obtaining the elevation of the landslide leading edge shear outlet is as follows:

[0022] Based on engineering geological data and borehole data, the topographic features, elevation data, groundwater level and geological structure information of the study area were analyzed and obtained. Landslide profile elevation maps of the study area were drawn, and the elevation of the shear outlet at the leading edge of the landslide was obtained from the landslide profile elevation maps.

[0023] Preferably, step S3 specifically includes:

[0024] S31: The reservoir landslides in the study area are divided into two types: hydrodynamic pressure type landslides and floating weight reduction type landslides. The cumulative displacement curve of the hydrodynamic pressure type landslide is a step-type cumulative displacement curve, while the cumulative displacement curve of the floating weight reduction type landslide is a straight-line cumulative displacement curve.

[0025] S32: Set the threshold for the soil-rock ratio of hydrodynamic pressure type landslides to SRR≥1, and set the threshold for the soil-rock ratio of floating weight reduction type landslides to SRR<1, where SRR is the soil-rock ratio;

[0026] S33: The sliding surface morphology of the hydrodynamic pressure type landslide is an arc-shaped sliding surface, while the sliding surface morphology of the floating weight reduction type landslide is a chair-shaped sliding surface.

[0027] S34: The threshold for the shear outlet elevation at the leading edge of the landslide for both hydrodynamic pressure type landslides and floating overload reduction type landslides is not higher than the maximum elevation of the reservoir water level drawdown zone.

[0028] Preferably, the construction process of the identification model is as follows:

[0029] The identification criteria for the identification model are set, including: a first set of conditions and a second set of conditions;

[0030] A landslide that simultaneously satisfies all conditions in the first set of conditions is a reservoir hydrodynamic pressure type landslide. The first set of conditions includes:

[0031] Step-type cumulative displacement curve; soil-to-rock ratio SRR≥1; arc-shaped sliding surface; landslide leading edge shear exit elevation not higher than the threshold;

[0032] A landslide that simultaneously satisfies all conditions in the second set of conditions is a reservoir-assisted floating reduction type landslide. The second set of conditions includes:

[0033] Linear cumulative displacement curve; soil-to-rock ratio SRR < 1; chair-shaped sliding surface; shear exit elevation at the leading edge of the landslide not higher than the threshold.

[0034] A rapid identification system for reservoir landslide types coupled with temporal InSAR technology includes:

[0035] The engineering geological condition data acquisition module is used to acquire engineering geological condition data and radar image sets for the study area;

[0036] The indicator type acquisition module is used to acquire four types of indicators for identifying landslide types in the study area through engineering geological condition data and radar image sets, including: cumulative displacement curve, soil-rock ratio, sliding surface morphology and landslide leading edge shear exit elevation.

[0037] The indicator data acquisition module is used to acquire four indicator data for different types of reservoir landslides in the study area through engineering geological condition data, radar image sets and four indicator types, including: cumulative displacement curve shape, soil-rock ratio threshold, type of sliding surface morphology and threshold of landslide leading edge shear outlet elevation.

[0038] The landslide type identification module is used to construct identification models for different types of reservoir landslides using four indicator data, and to identify the landslide types in the study area through the identification models.

[0039] The present invention has the following beneficial effects:

[0040] 1. The identification model constructed through indicator types and indicator data can accurately and quickly identify different types of reservoir landslides, and then take targeted early warning plans, which has strong versatility;

[0041] 2. The identification criteria in the judgment model are specifically set for different types of landslides, which can further improve the accuracy of landslide disaster prevention and risk management, meet the urgent needs of areas with a wide distribution of reservoir landslides for identifying different types of reservoir landslides, lay a solid foundation for regional disaster prevention and mitigation work, and are of great significance for protecting people's lives and property. Attached Figure Description

[0042] Figure 1 This is a flowchart of a method according to an embodiment of the present invention;

[0043] Figure 2 Schematic diagrams of the cumulative displacement curve shapes for step-type and linear cumulative displacement curves;

[0044] Figure 3 A simplified cross-sectional diagram of the sliding surface morphology of the arc-shaped sliding surface and the chair-shaped sliding surface;

[0045] The realization of the objective, functional features and advantages of the present invention will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0046] It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.

[0047] Reference Figure 1 This invention provides a rapid identification method for reservoir landslide types coupled with temporal InSAR technology, which can further improve the accuracy of landslide disaster prevention and risk management, meet the urgent needs of areas with a wide distribution of reservoir landslides for identifying different types of reservoir landslides, lay a solid foundation for regional disaster prevention and mitigation work, and is of great significance for protecting people's lives and property.

[0048] Includes the following steps:

[0049] S1: Acquire engineering geological condition data and radar image sets for the study area;

[0050] S2: Four types of indicators are obtained through engineering geological data and radar image sets to identify landslide types in the study area, including: cumulative displacement curve, soil-rock ratio, sliding surface morphology and landslide leading edge shear exit elevation;

[0051] S3: By using engineering geological data, radar image sets, and four types of indicators, four types of indicators were obtained for different types of reservoir landslides in the study area, including: cumulative displacement curve shape, soil-rock ratio threshold, type of sliding surface morphology, and threshold of landslide leading edge shear outlet elevation.

[0052] S4: Construct identification models for different types of reservoir landslides using four indicator data, and identify landslide types in the study area using the identification models.

[0053] Furthermore, the engineering geological data includes: the number and geographical location of existing landslides in the study area, the climate characteristics, geomorphological characteristics and stratigraphic lithology characteristics of the study area, and the rainfall data and reservoir water level scheduling data of the study area over the past ten years.

[0054] Specifically, the study area of ​​this embodiment is six typical landslide areas in the Three Gorges Reservoir area. Existing literature was extensively collected, including research papers, geological survey reports and relevant engineering construction documents that take the six landslides in the study area as research objects. Field investigations were also conducted to clarify the geographical location and scale of the six landslides in the study area, and to investigate the climate characteristics, geomorphological characteristics and stratigraphic lithology characteristics of the site, as well as the rainfall data and reservoir water level scheduling data for the past ten years.

[0055] Furthermore, time-series InSAR technology is used to acquire image data;

[0056] The process of obtaining the cumulative displacement curve is as follows:

[0057] We acquired a set of Sentinel-1A radar images of the slopes in the study area over the past five years. We used a digital elevation model to remove the terrain phase and geocode the radar images to obtain interferograms. We used a small baseline set method to process the interferograms and calculated the surface deformation information in the study area with millimeter-level accuracy. We then used the surface deformation information to draw the cumulative displacement curve of the study area.

[0058] Furthermore, based on ground surveys and data revealed by exploration boreholes, the stratigraphic and lithological characteristics of the study area were preliminarily analyzed; the thickness of the soil and gravel layers was calculated based on the borehole data; and the soil-to-rock ratio of the landslide deposits in the study area was calculated using algebraic methods.

[0059] The process for obtaining the soil-rock ratio is as follows:

[0060] By analyzing the lithological characteristics of the strata, the soil-to-rock ratio of the landslide deposits in the study area was calculated using the following formula:

[0061]

[0062]

[0063]

[0064] in, To determine the thickness of the soil layer sampled from the borehole, The value represents the thickness of the gravel layer sampled from the borehole, where i is the borehole number and n is the total number of boreholes. The ratio of soil to stone.

[0065] Furthermore, the process of obtaining the smooth surface morphology is as follows:

[0066] Based on engineering geological data and borehole data, we analyzed and obtained information on the topographic features, lithological characteristics, rock component boundaries and geological structure of the study area, drew landslide profile maps of the study area, and obtained the sliding surface morphology based on the landslide profile maps.

[0067] Furthermore, the process for obtaining the elevation of the shear outlet at the leading edge of the landslide is as follows:

[0068] Based on engineering geological data and borehole data, the topographic features, elevation data, groundwater level and geological structure information of the study area were analyzed and obtained. Landslide profile elevation maps of the study area were drawn, and the elevation of the shear outlet at the leading edge of the landslide was obtained from the landslide profile elevation maps.

[0069] Furthermore, based on reservoir water level scheduling data and InSAR image data obtained from regional surveys, the cumulative displacement curve shape of different types of landslides was determined.

[0070] Based on borehole survey data and soil and rock characteristics obtained from regional surveys, the threshold values ​​for the soil-to-rock ratio of the deposited material were determined for different types of landslides.

[0071] Based on the profile data of landslides in the study area obtained from the regional survey, the sliding surface morphology of different types of landslides was determined.

[0072] Based on the profile data of landslides in the study area obtained from the regional survey, the threshold values ​​for the leading shear exit elevation of different types of landslides were determined.

[0073] Step S3 is as follows:

[0074] S31: The reservoir landslides in the study area are divided into two types: hydrodynamic pressure type landslides and floating weight reduction type landslides. The cumulative displacement curve of the hydrodynamic pressure type landslide is a step-type cumulative displacement curve, while the cumulative displacement curve of the floating weight reduction type landslide is a straight-line cumulative displacement curve.

[0075] Specifically, the cumulative displacement curves of the step-type and linear cumulative displacement curves are as follows: Figure 2 As shown;

[0076] S32: Set the threshold for the soil-rock ratio of hydrodynamic pressure type landslides to SRR≥1, and set the threshold for the soil-rock ratio of floating weight reduction type landslides to SRR<1, where SRR is the soil-rock ratio;

[0077] S33: The sliding surface morphology of the hydrodynamic pressure type landslide is an arc-shaped sliding surface, while the sliding surface morphology of the floating weight reduction type landslide is a chair-shaped sliding surface.

[0078] Specifically, simplified cross-sectional diagrams of the sliding surface morphology of the arc-shaped sliding surface and the chair-shaped sliding surface are shown below. Figure 3 As shown;

[0079] S34: The threshold for the elevation of the shear outlet at the leading edge of the landslide is no higher than the maximum elevation of the reservoir water level drawdown zone for both dynamic water pressure type landslides and floating overload reduction type landslides. Specifically, the maximum elevation of the reservoir water level drawdown zone in the Three Gorges Reservoir area is 175m, that is, the elevation of the shear outlet at the leading edge of the landslide is no higher than 175m.

[0080] Furthermore, identification criteria for reservoir hydrodynamic pressure type landslides and floating-over-reduced type landslides are established. Specifically, based on the cumulative displacement curve shape, soil-rock ratio threshold, sliding surface morphology, and landslide leading edge shear outlet elevation threshold of reservoir landslides, identification criteria for reservoir hydrodynamic pressure type landslides and floating-over-reduced type landslides are established in combination with the geological and geomorphological characteristics and failure mechanisms of six typical landslides.

[0081] The process of constructing the identification model is as follows:

[0082] The identification criteria for the identification model are set, including: a first set of conditions and a second set of conditions;

[0083] A landslide that simultaneously satisfies all conditions in the first set of conditions is a reservoir hydrodynamic pressure type landslide. The first set of conditions includes:

[0084] Step-type cumulative displacement curve; soil-rock ratio SRR≥1; arc-shaped sliding surface; landslide leading edge shear outlet elevation not higher than the threshold (landslide leading edge shear outlet elevation cannot be higher than 175m, the maximum elevation of the Three Gorges Reservoir water level drawdown zone).

[0085] A landslide that simultaneously satisfies all conditions in the second set of conditions is a reservoir-assisted floating reduction type landslide. The second set of conditions includes:

[0086] Linear cumulative displacement curve; soil-to-rock ratio (SRR) < 1; chair-shaped sliding surface; landslide leading edge shear outlet elevation not higher than the threshold (landslide leading edge shear outlet elevation cannot be higher than 175m, the maximum elevation of the Three Gorges Reservoir water level drawdown zone).

[0087] This invention provides a rapid identification system for reservoir landslide types coupled with temporal InSAR technology, comprising:

[0088] The engineering geological condition data acquisition module is used to acquire engineering geological condition data and radar image sets for the study area;

[0089] The indicator type acquisition module is used to acquire four types of indicators for identifying landslide types in the study area through engineering geological condition data and radar image sets, including: cumulative displacement curve, soil-rock ratio, sliding surface morphology and landslide leading edge shear exit elevation.

[0090] The indicator data acquisition module is used to acquire four indicator data for different types of reservoir landslides in the study area through engineering geological condition data, radar image sets and four indicator types, including: cumulative displacement curve shape, soil-rock ratio threshold, type of sliding surface morphology and threshold of landslide leading edge shear outlet elevation.

[0091] The landslide type identification module is used to construct identification models for different types of reservoir landslides using four indicator data, and to identify the landslide types in the study area through the identification models.

[0092] It should be noted that, in this document, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or system that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or system. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or system that includes that element.

[0093] The sequence numbers of the above embodiments of the present invention are for descriptive purposes only and do not represent the superiority or inferiority of the embodiments. In the unit claims listing several devices, several of these devices may be embodied by the same hardware item. The use of the terms first, second, and third, etc., does not indicate any order and can be interpreted as identifiers.

[0094] The above are merely preferred embodiments of the present invention and do not limit the scope of the patent. Any equivalent structural or procedural transformations made based on the description and drawings of the present invention, or direct or indirect applications in other related technical fields, are similarly included within the scope of patent protection of the present invention.

Claims

1. A reservoir landslide type rapid identification method coupled with a timing sequence InSAR technology, characterized in that, The method comprises the following steps: S1: obtaining engineering geology condition data and a radar image set of a study area; S2: obtaining four index categories for identifying landslide types in the study area from the engineering geology condition data and the radar image set, including a cumulative displacement curve, a soil-rock ratio, a sliding surface morphology, and a shear outlet elevation of a landslide front edge; S3: obtaining four index data of different types of reservoir landslides in the study area from the engineering geology condition data, the radar image set, and the four index categories, including a cumulative displacement curve shape, a soil-rock ratio threshold value, a sliding surface morphology category, and a shear outlet elevation threshold value of the landslide front edge; S4: constructing a discrimination model of different types of reservoir landslides from the four index data, and identifying landslide types in the study area through the discrimination model; The construction process of the discrimination model is as follows: setting identification criteria of the discrimination model, including a first condition set and a second condition set; a reservoir dynamic water pressure type landslide is obtained when all conditions in the first condition set are met, and the first condition set includes: a step-type cumulative displacement curve, a soil-rock ratio SRR≥1, an arc-shaped sliding surface, and a shear outlet elevation of the landslide front edge not higher than a threshold value; a reservoir floatation and weight reduction type landslide is obtained when all conditions in the second condition set are met, and the second condition set includes: a straight-line type cumulative displacement curve, a soil-rock ratio SRR<1, a recliner-shaped sliding surface, and a shear outlet elevation of the landslide front edge not higher than a threshold value.

2. The reservoir landslide type rapid identification method of coupled time-series InSAR technology according to claim 1, characterized in that, The engineering geology condition data includes the number and geographical location of existing landslides in the study area, climate characteristics, landform characteristics, and stratum lithology characteristics of the study area, rainfall data in the study area in the past ten years, and reservoir water level scheduling data.

3. The reservoir landslide type rapid identification method of coupled time-series InSAR technology according to claim 1, characterized in that, The cumulative displacement curve is obtained as follows: a Sentinel-1A radar image set of slopes in the study area in the past five years is obtained, a digital elevation model is used to remove terrain phases and perform geographical coding processing on the radar image set to obtain an interferogram, a small baseline set method is used to process the interferogram, millimeter-level surface deformation information in the study area is calculated and obtained, and a cumulative displacement curve of the study area is drawn through the surface deformation information.

4. The reservoir landslide type rapid identification method of coupled time-series InSAR technology according to claim 1, characterized in that, The soil-rock ratio is obtained as follows: the soil-rock ratio of landslide accumulation layer material in the study area is calculated through analysis of stratum lithology characteristics, and the calculation formula is as follows: wherein, is the thickness of the soil layer sampled by the borehole, is the thickness of the gravel layer sampled by the borehole, i is the number of the exploration borehole, and n is the total number of the exploration boreholes, is the soil-rock ratio.

5. The reservoir landslide type rapid identification method of coupled time-series InSAR technology according to claim 1, characterized in that, The sliding surface morphology is obtained as follows: terrain and landform characteristics, lithology characteristics, rock group boundaries, and geological structure information of the study area are analyzed and obtained from the engineering geology condition data and drilling data, a landslide profile of the study area is drawn, and the sliding surface morphology is obtained from the landslide profile.

6. The reservoir landslide type rapid identification method of coupled time-series InSAR technology according to claim 1, characterized in that, The shear outlet elevation of the landslide front edge is obtained as follows: terrain and landform characteristics, elevation data, groundwater level, and geological structure information of the study area are analyzed and obtained from the engineering geology condition data and drilling data, a landslide profile elevation map of the study area is drawn, and the shear outlet elevation of the landslide front edge is obtained from the landslide profile elevation map.

7. The reservoir landslide type rapid identification method of coupled time-series InSAR technology according to claim 1, characterized in that, Step S3 is specifically as follows: S31: dividing reservoir landslides in the study area into a dynamic water pressure type landslide and a floatation and weight reduction type landslide, the cumulative displacement curve shape of the dynamic water pressure type landslide is a step-type cumulative displacement curve, and the cumulative displacement curve shape of the floatation and weight reduction type landslide is a straight-line type cumulative displacement curve; S32: setting a threshold value of the soil-rock ratio of the dynamic water pressure type landslide to SRR≥1, and setting a threshold value of the soil-rock ratio of the floatation and weight reduction type landslide to SRR<1, wherein SRR is the soil-rock ratio; S33: the sliding surface of the dynamic water pressure type landslide is an arc-shaped sliding surface, and the sliding surface of the floatation and weight reduction type landslide is a recliner-shaped sliding surface; S34: the threshold value of the shear outlet elevation of the landslide front edge of the dynamic water pressure type landslide and the floatation and weight reduction type landslide is not higher than the maximum elevation of the drawdown zone of the reservoir water level.

8. A reservoir landslide type rapid identification system coupled with a time-series InSAR technology, for implementing the method of any one of claims 1-7. Comprise: an engineering geological condition data acquisition module configured to acquire engineering geological condition data and a radar image set of a study area; an index category acquisition module configured to acquire four index categories for identifying landslide types in the study area from the engineering geological condition data and the radar image set, including: a cumulative displacement curve, a soil-rock ratio, a sliding surface morphology, and a shear outlet elevation of a landslide front edge; an index data acquisition module configured to acquire four index data of different types of reservoir landslides in the study area from the engineering geological condition data, the radar image set, and the four index categories, including: a cumulative displacement curve shape, a threshold value of a soil-rock ratio, a type of a sliding surface morphology, and a threshold value of a shear outlet elevation of a landslide front edge; a landslide type identification module configured to construct an identification model of different types of reservoir landslides by the four index data, and identify the landslide types in the study area by the identification model.

Images