Landslide local instability region boundary identification method based on surface scene deformation monitoring

Abstract

The present application relates to the technical field of geological disaster prevention, and particularly relates to a landslide local instability region boundary identification method based on surface scene deformation monitoring, S1 obtains global surface scene three-dimensional deformation history monitoring data of a landslide body of a landslide hidden danger point; S2 takes the latest surface scene deformation monitoring data as a criterion, and constructs a landslide hidden danger point local instability region target area based on a first preset method; S3 identifies a landslide hidden danger point local instability region key resolution unit in the target area region based on a second preset method; S4 clusters and groups the key resolution unit based on a third preset method, and classifies the key resolution unit into different local instability regions; S5 takes a local instability region containing a number of key resolution units less than or equal to a set number as an abnormal point and removes the abnormal point; S6 takes a peripheral connecting line containing a key resolution unit in a local instability region as a boundary of the local instability region, and obtains a boundary of one or more local instability regions on a landslide body. The present application can timely and accurately identify the boundary of a local instability region.

Description

Technical Field

[0001] This invention relates to the field of geological disaster prevention and control technology, specifically to a method for identifying the boundary of a local instability zone in a landslide based on surface scene deformation monitoring. Background Technology

[0002] Mountainous areas account for about two-thirds of my country's land area, with complex geological and geographical conditions and nearly 300,000 landslides. They are prone to landslide disasters. Effective monitoring of landslide deformation and evolution processes, and timely and accurate early warning and forecasting of local instability and overall failure of landslide bodies are important means to effectively avoid casualties and property losses.

[0003] Currently, monitoring of landslide deformation processes primarily relies on point-based monitoring equipment such as GNSS and crack gauges deployed along profile lines. This limited range of point-based equipment struggles to capture the deformation state of the entire landslide, particularly failing to meet the monitoring and early warning needs for localized landslide instability, resulting in high rates of missed and false alarms. In recent years, with the rapid development of BeiDou navigation satellite and radar technology, three-dimensional millimeter-level surface deformation technology and equipment based on BeiDou navigation satellite reflected signals have emerged and begun to be used for three-dimensional deformation monitoring of landslide surfaces. This enables the acquisition of three-dimensional deformation information across the entire surface during landslide evolution, laying the foundation for monitoring and early warning of localized landslide instability and failure events.

[0004] However, in surface-level scenarios, the local instability of a landslide is related to the surrounding environment, making it impossible to determine the boundaries of the local instability area in a timely and accurate manner. Summary of the Invention

[0005] The present invention aims to provide a method for identifying the boundary of a local instability area of ​​a landslide based on surface scene deformation monitoring, so as to solve the problem of being unable to determine the boundary of a local instability area of ​​a landslide in a timely and accurate manner.

[0006] The method for identifying the boundary of a landslide's local instability zone based on surface scene deformation monitoring in this scheme includes the following steps:

[0007] S1, acquire historical three-dimensional deformation data of the entire surface scene of the landslide body at the landslide hazard point. The historical monitoring data includes the center point coordinates of each resolution unit on the landslide body surface and deformation time series data.

[0008] S2, based on the latest surface scene deformation monitoring data, construct the target area of ​​the local instability zone of the landslide hazard point based on the first preset method;

[0009] S3, based on the second preset method, within the target area of ​​the local instability zone of the landslide hazard point constructed in step S2, identify the key distinguishing unit of the local instability zone of the landslide hazard point;

[0010] S4. Based on the third preset method, the key distinguishing units of the local instability area of ​​the landslide hazard point described in step S3 are clustered and grouped, and classified into different local instability areas.

[0011] S5, Local unstable regions containing fewer than or equal to a set number of key resolution units are removed as outliers and are not recorded as local unstable regions.

[0012] S6. For a local instability zone containing more than a set number of key resolution units, for each local instability zone, the outer line connecting the key resolution units contained in the local instability zone is taken as the boundary of the local instability zone, thereby obtaining the boundary of one or more local instability zones on the landslide body of the landslide hazard point.

[0013] The beneficial effects of this plan are:

[0014] By utilizing historical 3D deformation data of landslide hazard points, key resolving units of local instability areas at landslide hazard points are accurately identified. These units are then clustered and grouped into different local instability zones. Anomalies are removed from the categorized local instability zones, and corresponding boundaries are delineated to avoid misclassification of local instability zones due to factors such as boulder falls or small-scale human activities. Key resolving units within local instability zones are used as control points to delineate the overall boundary of the local instability zone, enabling timely and accurate identification of landslide local instability zone boundaries based on surface scene deformation monitoring.

[0015] Furthermore, in step S2, the first preset method includes the following steps:

[0016] S2.1, Based on the latest period of surface scene deformation monitoring data, obtain the cumulative deformation value of each resolution unit on the landslide surface;

[0017] S2.2 Sort all resolvable units on the surface of the landslide body in descending order of cumulative deformation value;

[0018] S2.3, select the landslide surface resolution units with larger deformation values ​​at the top of the sorting and mark them with the first identifier;

[0019] S2.4, take the set of all distinguishable units marked as the first identifier on the surface of the landslide body as the target area of ​​the local instability zone of the landslide hazard point.

[0020] The beneficial effects are as follows: By constructing target areas for local instability zones of landslide hazard points using the latest phase of surface scene deformation monitoring data, the accuracy of the target areas for local instability zones of landslide hazard points can be ensured by maximizing the use of historical 3D deformation data of the surface scene at the landslide hazard points. Furthermore, by sorting and obtaining the highest-resolution unit with a predetermined ratio as the target area, both scenarios of overall landslide body stability and localized deformation (i.e., slow overall deformation rate and fast localized deformation rate) can be considered, thus ensuring that the constructed target areas for local instability zones of landslide hazard points are more comprehensive and accurate.

[0021] Furthermore, the set ratio is set according to the calculation requirements.

[0022] Furthermore, in step S3, the second preset method includes the following steps:

[0023] S3.1, for any landslide surface resolution unit in the target area of ​​the local instability zone of the landslide hazard point described in S2.4, based on the deformation time series data obtained in S1, the deformation sequence of the resolution unit in S2.4 is set as follows: ;

[0024] S3.2, Solve for the resolvable element deformation sequence in S3.1. Cumulative deviation sequence The formula is:

[0025] ;

[0026] In the formula, The sequence length; This represents the average deformation.

[0027] S3.3, the cumulative deviation sequence is divided into sub-intervals using both sequential and reverse order methods;

[0028] S3.4, Calculate the fluctuation function value under the given subinterval length condition. :

[0029] ;

[0030] In the formula, The length of the sub-interval, The number of subintervals For the first The variance of each subinterval;

[0031] S3.5, Change the length of the subinterval Repeat steps S3.3 and S3.4 to obtain different interval lengths. The corresponding fluctuation function value Several numerical pairs were obtained. And ensure that the number of values ​​is greater than or equal to 5;

[0032] S3.6, Based on the several values ​​obtained in step S3.5, By constructing a functional relationship and fitting it using the least squares method, the classical Hurst exponent of the resolving unit is obtained. The functional relationship is expressed as:

[0033] ;

[0034] In the formula, C is a constant;

[0035] S3.7, For all surface resolution units of landslide bodies in the target area of ​​the local instability zone of the landslide hazard point described in step S2.4, repeat steps S3.1 to S3.6 to obtain the classical Hurst index of all surface resolution units of landslide bodies in the target area of ​​the local instability zone of the landslide hazard point. ;

[0036] S3.8, Classical Hurst index, a resolution unit for the surface of the landslide body in the target area of ​​the local instability zone of the landslide hazard point. All resolvable units with values ​​within a set range are labeled with a second identifier. All resolvable units labeled with the second identifier are the key resolvable units of the local instability area of ​​the landslide hazard point.

[0037] The beneficial effects are: by using operations such as cumulative deviation sequences, sub-interval division, and function construction, the accuracy of identifying key distinguishing units in the local instability area of ​​landslide hazard points can be improved.

[0038] Furthermore, in S3.3, the cumulative deviation sequence is sequentially... Perform sub-interval partitioning to obtain the number of sub-intervals. Then according to length Reverse the cumulative deviation sequence Perform sub-interval partitioning to obtain the number of sub-intervals. A total of Sub-intervals, where when If the value is not an integer, take the integer value.

[0039] The beneficial effect is that the division of sub-intervals can fully utilize the role of all statistical nodes in the resolving unit deformation sequence, thereby improving the accuracy of calculation.

[0040] Furthermore, in S3.5, the length of the sub-interval The value range is 2 to N / 2, and the length of the sub-interval changes from small to large within the value range and is evenly distributed within the value range.

[0041] The beneficial effect is that the change in the length of the sub-interval can completely capture the fluctuation function carrying the corresponding landslide information.

[0042] Furthermore, in S3.8, the set interval is... .

[0043] The beneficial effects are: by filtering the resolution units within the set interval, searching for resolution units with a continuous deformation trend in the target area of ​​the local instability zone of the landslide hazard point, the key resolution units of the local instability zone of the landslide hazard point can be accurately identified.

[0044] Furthermore, in step S4, the third preset method includes the following steps:

[0045] S4.1, the key identification units of the local instability area of ​​the landslide hazard point marked as the second identifier as described in S3.8 are numbered sequentially in ascending order;

[0046] S4.2, classify the first sequentially numbered resolving unit into the first local instability region, and record it as the resolving unit search center of the first local instability region;

[0047] S4.3, calculate the distance between the coordinates of the center point of the search center resolution cell and the coordinates of the center points of all remaining resolution cells marked with the second identifier. The calculation formula is:

[0048] ;

[0049] In the formula, To determine the coordinates of the unit center point, , To distinguish the unit number, and have ;

[0050] S4.4, classify the resolution units that meet the preset conditions into the first local instability region;

[0051] S4.5, sequentially classify the distinguishing units that meet the preset conditions described in step S4.4 into the first local instability region, and which have not previously served as the distinguishing unit search center of the first local instability region, and re-record them as the distinguishing units search center of the first local instability region. Repeat steps S4.2 to S4.4 until all distinguishing units that meet the conditions are classified into the first local instability region.

[0052] S4.6 For the remaining unclassified distinguishing units, renumber them in ascending order. Repeat steps S4.2 to S4.5 to classify the remaining unclassified distinguishing units into other local instability zones except the first local instability zone. Cluster all the key distinguishing units of the local instability zones of landslide hazard points marked with the second identifier into n local instability zones.

[0053] The beneficial effect is that, by using complete clustering and grouping, key distinguishing units belonging to different local instability zones can be segmented, enabling accurate identification of multiple local instability zones on the landslide body.

[0054] Furthermore, in S4.4, the side length of the resolving unit is denoted as... The preset conditions are or .

[0055] The beneficial effect is that by constraining from both the absolute distance and the side length of the resolution unit, the needs of radars in different resolution surface scenes can be met, thus enabling more accurate clustering. Attached Figure Description

[0056] Figure 1 This is a flowchart illustrating an embodiment of the landslide local instability zone boundary identification method based on surface scene deformation monitoring according to the present invention. Detailed Implementation

[0057] The following detailed description provides further details on specific implementation methods.

[0058] Example

[0059] A method for identifying the boundary of a local instability zone in a landslide based on surface scene deformation monitoring, such as Figure 1 As shown, it includes the following steps:

[0060] S1. Obtain historical three-dimensional deformation data of the entire surface scene of the landslide body at the landslide hazard point. The historical monitoring data includes the center point coordinates and deformation time series data of each resolution unit on the landslide body surface. This resolution unit is related to the surface scene deformation monitoring radar and is an inherent attribute of the surface scene deformation monitoring radar, which will not be elaborated here.

[0061] S2, based on the latest period's surface scene deformation monitoring data, construct the target area of ​​the local instability zone of the landslide hazard point based on the first preset method. The first preset method includes the following steps:

[0062] S2.1, Based on the latest period of surface scene deformation monitoring data, obtain the cumulative deformation value of each resolution unit on the landslide surface;

[0063] S2.2 Sort all resolvable units on the surface of the landslide body in descending order of cumulative deformation value;

[0064] S2.3, select the landslide surface resolution units with larger deformation values ​​and sort them first according to a set proportion, and mark them with the first identifier as the label. The set proportion is set according to the calculation requirements, such as 25%, 30%, 50%, etc. In this embodiment, the set proportion is set to 25%. The first identifier can be represented by letters or numbers, for example, the first identifier is set to b.

[0065] S2.4, take the set of all distinguishable units marked as the first identifier on the surface of the landslide body as the target area of ​​the local instability zone of the landslide hazard point.

[0066] S3, based on the second preset method, within the target area of ​​the local instability zone of the landslide hazard point constructed in step S2, identify the key distinguishing units of the local instability zone of the landslide hazard point. The second preset method includes the following steps:

[0067] S3.1, for any landslide surface resolution unit in the target area of ​​the local instability zone of the landslide hazard point described in S2.4, based on the deformation time series data obtained in S1, the deformation sequence of the resolution unit in S2.4 is set as follows: ;

[0068] S3.2, Solve for the resolvable element deformation sequence in S3.1. Cumulative deviation sequence The solution formula is:

[0069] ;

[0070] In the formula, The sequence length; This represents the average deformation.

[0071] S3.3, the cumulative deviation sequence is divided into sub-intervals using both sequential and reverse methods. Specifically: the cumulative deviation sequence is divided into sub-intervals sequentially. Perform sub-interval partitioning to obtain the number of sub-intervals. Then according to length Reverse the cumulative deviation sequence Perform sub-interval partitioning to obtain the number of sub-intervals. A total of Sub-intervals, where when If the value is not an integer, take the integer value.

[0072] S3.4, Calculate the fluctuation function value under the given subinterval length condition. :

[0073] ;

[0074] In the formula, The length of the sub-interval, The number of subintervals;

[0075] S3.5, Change the length of the subinterval Interval length The value range is 2 to N / 2, and the interval length is... The values ​​are changed from small to large within the range, and then evenly distributed within the range. Steps S3.3 and S3.4 are repeated to obtain different interval lengths. The corresponding fluctuation function value Several numerical pairs were obtained. And ensure that the number of values ​​is greater than or equal to 5;

[0076] S3.6, Based on the several values ​​obtained in step S3.5, By constructing a functional relationship and fitting it using the least squares method, the classical Hurst exponent of the resolving unit is obtained. The functional relationship is expressed as:

[0077] ;

[0078] In the formula, C is a constant;

[0079] S3.7, For all surface resolution units of landslide bodies in the target area of ​​the local instability zone of the landslide hazard point described in step S2.4, repeat steps S3.1 to S3.6 to obtain the classical Hurst index of all surface resolution units of landslide bodies in the target area of ​​the local instability zone of the landslide hazard point. ;

[0080] S3.8, Classical Hurst index, a resolution unit for the surface of the landslide body in the target area of ​​the local instability zone of the landslide hazard point. All resolvable units with values ​​within a set range are marked with a second identifier, and the set range is... The second identifier can be represented by uppercase or lowercase letters or numbers. For example, if the second identifier is set to B, all the distinguishing units marked with the second identifier are the key distinguishing units of the local instability area of ​​the landslide hazard point.

[0081] S4. Based on the third preset method, the key distinguishing units of the local instability area of ​​the landslide hazard point mentioned in step S3 are clustered and grouped, and classified into different local instability areas. The third preset method includes the following steps:

[0082] S4.1, The key identification units of the local instability area of ​​the landslide hazard point marked B as described in S3.8 are numbered sequentially in ascending order. The numbering can be done by letters or numbers, for example, the numbering is done in ascending order of 1, 2, 3, 4, ...

[0083] S4.2, classify the first resolving unit in the sequence number into the first local instability region, that is, classify the resolving unit numbered 1 into the first local instability region. The first local instability region is denoted as LF-1 and recorded as the search center of the resolving unit of the first local instability region LF-1.

[0084] S4.3, calculate the distance between the coordinates of the center point of the search center resolution cell and the coordinates of the center points of all remaining resolution cells labeled B. The calculation formula is:

[0085] ;

[0086] In the formula, To determine the coordinates of the unit center point, , To distinguish the unit number, and have ;

[0087] S4.4, classify the resolving elements that meet the preset conditions into the first local instability region LF-1, and denote the side length of the resolving element as... The preset conditions are or ;

[0088] S4.5, sequentially re-record the resolving units that meet the preset conditions described in step S4.4 and are classified into the first local instability zone LF-1, and which have not previously served as the search center of the resolving unit in the first local instability zone LF-1, as the search center of the resolving unit in the first local instability zone LF-1. Repeat steps S4.2 to S4.4 until all resolving units that meet the conditions are classified into the first local instability zone LF-1. That is, each resolving unit classified into the first local instability zone LF-1 in S4.4 is used as the search center of the resolving unit again for search and distance calculation to avoid missing information points.

[0089] S4.6 For the remaining unclassified distinguishable units, renumber them in ascending order, using letters or numbers, such as 1, 2, 3, 4, ... Repeat steps S4.2 to S4.5 to classify the remaining unclassified distinguishable units into other local instability zones except the first local instability zone LF-1, i.e., into LF-2, LF-3, LF-4, ..., LF-n. Cluster all the key distinguishable units of the local instability areas of landslide hazard points marked as B into n local instability zones.

[0090] S5, local unstable areas containing fewer than or equal to a set number of key resolution units are removed as outliers and are not recorded as local unstable areas. The set number is set according to actual needs. The larger the resolution units, the fewer the set number, and the smaller the resolution units, the more the set number. For example, the set number is set to 3.

[0091] S6. For a local instability zone containing more than a set number of key resolution units, for each local instability zone, the outer line connecting the key resolution units contained in the local instability zone is taken as the boundary of the local instability zone, thereby obtaining the boundary of one or more local instability zones on the landslide body of the landslide hazard point.

[0092] This embodiment acquires historical monitoring data based on a surface scene, constructs a target area for the local instability zone of a landslide hazard point using the latest surface scene monitoring data, identifies key resolution units (KJUs) within the target area, and then clusters and groups these KJUs to classify them into different local instability zones. Based on the number of KJUs in each local instability zone, the local instability zones are re-screened, and the boundaries of several local instability zones are divided according to the re-screened local instability zones. This can improve the timeliness and accuracy of identifying the boundaries of local instability zones of landslides based on surface scene deformation monitoring, while ensuring that large-scale local instability of landslides does not occur.

[0093] The above descriptions are merely embodiments of the present invention, and common knowledge regarding specific structures and characteristics is not elaborated upon here. It should be noted that those skilled in the art can make various modifications and improvements without departing from the structure of the present invention, and these should also be considered within the scope of protection of the present invention. These modifications and improvements will not affect the effectiveness of the present invention or the practicality of the patent. The scope of protection claimed in this application should be determined by the content of its claims, and the specific embodiments described in the specification can be used to interpret the content of the claims.

Claims

1. A method for identifying the boundary of a landslide's local instability zone based on surface scene deformation monitoring, characterized in that, Includes the following steps: S1, acquire historical three-dimensional deformation data of the entire surface scene of the landslide body at the landslide hazard point. The historical monitoring data includes the center point coordinates of each resolution unit on the landslide body surface and deformation time series data. S2, based on the latest surface scene deformation monitoring data, construct the target area of ​​the local instability zone of the landslide hazard point based on the first preset method; The first preset method includes the following steps: S2.1, Based on the latest period of surface scene deformation monitoring data, obtain the cumulative deformation value of each resolution unit on the landslide surface; S2.2 Sort all resolvable units on the surface of the landslide body in descending order of cumulative deformation value; S2.3, select the landslide surface resolution units with larger deformation values ​​at the top of the sorting and mark them with the first identifier; S2.4, take the set of all distinguishable units marked as the first identifier on the surface of the landslide body as the target area of ​​the local instability zone of the landslide hazard point; S3, based on the second preset method, within the target area of ​​the local instability zone of the landslide hazard point constructed in step S2, identify the key distinguishing unit of the local instability zone of the landslide hazard point; The second preset method includes the following steps: S3.1, for any landslide surface resolution unit in the target area of ​​the local instability zone of the landslide hazard point described in S2.4, based on the deformation time series data obtained in S1, the deformation sequence of the resolution unit in S2.4 is set as follows: ; S3.2, Solve for the resolvable element deformation sequence in S3.

1. Cumulative deviation sequence The formula is: ； In the formula, The sequence length; This represents the average deformation. S3.3, the cumulative deviation sequence is divided into sub-intervals using both sequential and reverse order methods; S3.4, Calculate the fluctuation function value under the given subinterval length condition. : ； In the formula, The length of the sub-interval, The number of subintervals For the first The variance of each subinterval; S3.5, Change the length of the subinterval Repeat steps S3.3 and S3.4 to obtain different interval lengths. The corresponding fluctuation function value Several numerical pairs were obtained. And ensure that the number of values ​​is greater than or equal to 5; S3.6, Based on the several values ​​obtained in step S3.5, By constructing a functional relationship and fitting it using the least squares method, the classical Hurst exponent of the resolving unit is obtained. The functional relationship is expressed as: ； In the formula, C is a constant; S3.7, For all surface resolution units of landslide bodies in the target area of ​​the local instability zone of the landslide hazard point described in step S2.4, repeat steps S3.1 to S3.6 to obtain the classical Hurst index of all surface resolution units of landslide bodies in the target area of ​​the local instability zone of the landslide hazard point. ; S3.8, Classical Hurst index, a resolution unit for the surface of the landslide body in the target area of ​​the local instability zone of the landslide hazard point. All distinguishable units within the set range are marked with a second identifier. All distinguishable units marked with the second identifier are the key distinguishable units of the local instability area of ​​the landslide hazard point. S4. Based on the third preset method, the key distinguishing units of the local instability area of ​​the landslide hazard point described in step S3 are clustered and grouped, and classified into different local instability areas. S5, Local unstable regions containing fewer than or equal to a set number of key resolution units are removed as outliers and are not recorded as local unstable regions. S6. For a local instability zone containing more than a set number of key resolution units, for each local instability zone, the outer line connecting the key resolution units contained in the local instability zone is taken as the boundary of the local instability zone, thereby obtaining the boundary of one or more local instability zones on the landslide body of the landslide hazard point.

2. The method for identifying the boundary of a landslide's local instability zone based on surface scene deformation monitoring according to claim 1, characterized in that: The set ratio is set according to the calculation requirements.

3. The method for identifying the boundary of a landslide's local instability zone based on surface scene deformation monitoring according to claim 1, characterized in that: In S3.3, the cumulative deviation sequence is sequentially paired. Perform sub-interval partitioning to obtain the number of sub-intervals. Then according to length Reverse the cumulative deviation sequence Perform sub-interval partitioning to obtain the number of sub-intervals. A total of Sub-intervals, where when If the value is not an integer, take the integer value.

4. The method for identifying the boundary of a landslide's local instability zone based on surface scene deformation monitoring according to claim 3, characterized in that: In S3.5, the length of the sub-interval The value range is 2 to N / 2, and the length of the sub-interval changes from small to large within the value range and is evenly distributed within the value range.

5. The method for identifying the boundary of a landslide's local instability zone based on surface scene deformation monitoring according to claim 3, characterized in that: In S3.8, the set interval is .

6. The method for identifying the boundary of a landslide's local instability zone based on surface scene deformation monitoring according to claim 1, characterized in that: In step S4, the third preset method includes the following steps: S4.1, the key identification units of the local instability area of ​​the landslide hazard point marked as the second identifier as described in S3.8 are numbered sequentially in ascending order; S4.2, classify the first sequentially numbered resolving unit into the first local instability region, and record it as the resolving unit search center of the first local instability region; S4.3, calculate the distance between the coordinates of the center point of the search center resolution cell and the coordinates of the center points of all remaining resolution cells marked with the second identifier. The calculation formula is: ； In the formula, To determine the coordinates of the unit center point, , To distinguish the unit number, and have ; S4.4, classify the resolution units that meet the preset conditions into the first local instability region; S4.5, sequentially classify the distinguishing units that meet the preset conditions described in step S4.4 into the first local instability region, and which have not previously served as the distinguishing unit search center of the first local instability region, and re-record them as the distinguishing units search center of the first local instability region. Repeat steps S4.2 to S4.4 until all distinguishing units that meet the conditions are classified into the first local instability region. S4.6 For the remaining unclassified distinguishing units, renumber them in ascending order. Repeat steps S4.2 to S4.5 to classify the remaining unclassified distinguishing units into other local instability zones except the first local instability zone. Cluster all the key distinguishing units of the local instability zones of landslide hazard points marked with the second identifier into n local instability zones.

7. The method for identifying the boundary of a landslide's local instability zone based on surface scene deformation monitoring according to claim 6, characterized in that: In S4.4, the side length of the resolving unit is denoted as... The preset conditions are: or .

Images