Slope disaster early warning method and system

Abstract

The invention belongs to the technical field of geological monitoring, and particularly discloses a slope disaster early warning method and system, and the method comprises the steps: carrying out the multi-dimensional data analysis through obtaining the position detection data of each monitoring point in a target slope region; determining a relative displacement vector, a dynamic strain tensor, a displacement rate cooperation parameter and a spatial correlation parameter between each pair of monitoring points; and then a comprehensive early warning index of the target slope area at the current detection time point is calculated by using the relative displacement vector, the dynamic strain tensor, the displacement rate cooperation parameter and the spatial correlation parameter between each pair of monitoring points, and disaster early warning of the target slope area is carried out according to the comprehensive early warning index. According to the invention, the deformation degree of the slope can be accurately grasped, the dynamic monitoring and evaluation of the slope condition are realized, the slope disaster monitoring accuracy and the early warning response efficiency are improved, and the slope disaster monitoring and early warning capability is enhanced.

Description

Technical Field

[0001] The invention belongs to the technical field of geological monitoring, and in particular relates to a slope disaster early warning method and system. Background Art

[0002] Slope disaster is a common geological disaster, which is the result of the combined effects of multiple factors such as geography, geology, climate, and human activities. Generally speaking, the formation mechanism of slope disasters is mainly shear failure, rupture failure, and liquefaction failure. According to the form and development process of landslides, slope disasters can be divided into landslides, collapses, mudslides, etc. Slope disaster monitoring and early warning is the process of monitoring and evaluating the deformation degree and stability of the slope, and timely warning the occurrence of slope disasters, so as to take response and control measures.

[0003] Existing slope disaster monitoring and early warning methods can be divided into static monitoring methods and dynamic monitoring methods. Static monitoring methods mainly include ground potential difference method, ground resistivity method, seismic exploration method and gravity measurement method, while dynamic monitoring methods mainly use tools such as inclinometers, strain gauges, level gauges and pressure gauges for dynamic measurement and analysis. However, these existing slope disaster monitoring and early warning methods generally have the problems of poor measurement timeliness, insufficient dynamic response capability, single analysis dimension and insufficient monitoring accuracy. Summary of the invention

[0004] The purpose of the present invention is to provide a slope disaster early warning method and system to solve the above problems existing in the prior art.

[0005] In order to achieve the above object, the present invention adopts the following technical solutions: In a first aspect, a slope disaster early warning method is provided, which is applied to a cloud platform and includes: Obtain the position detection data of each monitoring point in the target slope area at the current detection time point; Determine the three-dimensional coordinates of each monitoring point at the current detection time point according to the position detection data of each monitoring point at the current detection time point; Determine the relative displacement vector and dynamic strain tensor between any two monitoring points at the current detection time point based on the three-dimensional coordinates of each monitoring point at the current detection time point; Retrieving the three-dimensional coordinates of each monitoring point at the previous detection time point, and determining the single-point displacement rate of each monitoring point within the time window from the previous detection time point to the current detection time point based on the three-dimensional coordinates of each monitoring point at the previous detection time point and the three-dimensional coordinates of each monitoring point at the current detection time point; According to the relative displacement vector between a pair of monitoring points at the current detection time point and the single-point displacement rate of the two monitoring points in the time window, the displacement rate coordination parameter between the two monitoring points at the current detection time point is calculated; Determine the angle parameter of the line connecting the two monitoring points according to the relative displacement vector between a pair of monitoring points at the current detection time point, and calculate the spatial correlation parameter between the two monitoring points at the current detection time point according to the angle parameter of the line connecting the two monitoring points; According to the relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points at the current detection time point, the comprehensive early warning index of the target slope area at the current detection time point is calculated; Disaster warning for the target slope area is carried out based on the comprehensive warning index of the target slope area at the current detection time point.

[0006] In a possible design, determining the three-dimensional coordinates of each monitoring point at the current detection time point according to the position detection data of each monitoring point at the current detection time point includes: The position detection data of each monitoring point at the current detection time point is projected into the set three-dimensional coordinate system to obtain the three-dimensional coordinates of each monitoring point at the current detection time point in the set three-dimensional coordinate system.

[0007] In a possible design, the determining of the relative displacement vector and the dynamic strain tensor between any two monitoring points at the current detection time point based on the three-dimensional coordinates of each monitoring point at the current detection time point includes: Substitute the three-dimensional coordinates of any two monitoring points at the current detection time point into the preset relative displacement vector formula for calculation to obtain the relative displacement vector between any two monitoring points at the current detection time point. The relative displacement vector formula is:

[0008] Where t represents the current detection time point, i and j are the numbers of any two monitoring points, and i≠j, ΔD ij (t) represents the relative displacement vector between monitoring points i and j at the current detection time point t, x i (t)、y i (t) and z i (t) is the three-dimensional coordinate of monitoring point i at the current detection time point, x j (t)、y j (t) and z j (t) is the three-dimensional coordinate of monitoring point j at the current detection time point; Substitute the relative displacement vector between any two monitoring points at the current detection time point into the preset dynamic strain tensor formula for calculation to obtain the dynamic strain tensor between any two monitoring points at the current detection time point. The dynamic strain tensor formula is:

[0009] Among them, ε ij(t) represents the dynamic strain tensor between monitoring points i and j at the current detection time point t, L ij is the initial distance between the given monitoring points i and j, and ‖‖ represents the vector modulo operation.

[0010] In a possible design, the method of determining the single-point displacement rate of each monitoring point within a time window from the previous detection time point to the current detection time point based on the three-dimensional coordinates of each monitoring point at the previous detection time point and the three-dimensional coordinates of each monitoring point at the current detection time point includes: Substitute the three-dimensional coordinates of a monitoring point at the previous detection time point and the three-dimensional coordinates of the current detection time point into the single-point displacement rate formula for calculation to obtain the single-point displacement rate of the monitoring point in the time window from the previous detection time point to the current detection time point. The single-point displacement rate formula is:

[0011] Among them, Δt represents the time interval from the previous detection time point to the current detection time point, t-Δt represents the previous detection time point, and v i (Δt) represents the single-point displacement rate of monitoring point i within the time window from the previous detection time point to the current detection time point.

[0012] In a possible design, the calculation of the displacement rate coordination parameter between the two monitoring points corresponding to the current detection time point according to the relative displacement vector between a pair of monitoring points at the current detection time point and the single-point displacement rate of the two monitoring points in the time window includes: Substitute the relative displacement vector between two monitoring points at the current detection time point and the single-point displacement rate of the two monitoring points in the time window into the preset displacement rate coordination parameter formula for calculation to obtain the displacement rate coordination parameter between the two monitoring points at the current detection time point. The displacement rate coordination parameter formula is:

[0013] Among them, δ ij (t) represents the displacement rate coordination parameter between monitoring point i and monitoring point j at the current detection time point, v i (Δt) represents the displacement rate of a single point at monitoring point i in the time window from the previous detection time point to the current detection time point, v j (Δt) represents the single-point displacement rate of monitoring point j within the time window from the previous detection time point to the current detection time point.

[0014] In a possible design, the method of determining the angle parameter of the line connecting the two monitoring points according to the relative displacement vector between a pair of monitoring points at the current detection time point, and calculating the spatial correlation parameter between the two monitoring points at the current detection time point according to the angle parameter of the line connecting the two monitoring points, includes: The angle between the relative displacement vector of a pair of monitoring points and the set main structural surface is used as the angle parameter of the line connecting the two monitoring points; Substitute the angle parameter of the line connecting the two corresponding monitoring points into the preset spatial correlation parameter formula for calculation to obtain the spatial correlation parameter between the two corresponding monitoring points. The spatial correlation parameter formula is:

[0015] Among them, i and j are the numbers of the corresponding two monitoring points, ω ij is the spatial correlation parameter between monitoring point i and monitoring point j, L ij is the initial distance between the given monitoring point i and the monitoring point j, θ ij is the angle parameter of the line connecting monitoring point i and monitoring point j, θ0 is the set reference angle constant, L0 is the set characteristic attenuation distance, and exp represents the natural exponential function.

[0016] In a possible design, the calculation of the comprehensive early warning index of the target slope area at the current detection time point according to the relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points at the current detection time point includes: The relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points at the current detection time point are substituted into the preset comprehensive warning index formula for calculation to obtain the comprehensive warning index of the target slope area at the current detection time point. The comprehensive warning index formula is:

[0017] Among them, S(t) is the comprehensive warning index of the target slope area at the current detection time point, i and j are the numbers of the corresponding two monitoring points, represents the combination of monitoring point i and monitoring point j, N is the number of all possible combinations of monitoring points, ω ij is the spatial correlation parameter between monitoring point i and monitoring point j, ΔD ij (t) is the relative displacement vector between monitoring points i and j at the current detection time point t, ‖‖ represents the vector modulus operation, ε ij (t) is the dynamic strain tensor between monitoring points i and j at the current detection time point t, δ ij(t) is the displacement rate coordination parameter between monitoring point i and monitoring point j at the current detection time point t, α is the set first adjustment coefficient, β is the set second adjustment coefficient, and γ is the set third adjustment coefficient.

[0018] In a possible design, the disaster warning of the target slope area based on the comprehensive warning index of the target slope area at the current detection time point includes: When it is determined that the comprehensive warning index of the target slope area at the current detection time point exceeds the set first threshold, the yellow disaster warning information of the target slope area is output; When it is determined that the comprehensive warning index of the target slope area at the current detection time point exceeds a set second threshold, red disaster warning information of the target slope area is output, and the second threshold is greater than the first threshold.

[0019] In a second aspect, a slope disaster warning system is provided, comprising a data acquisition unit, a coordinate mapping unit, a first determination unit, a second determination unit, a third determination unit, a fourth determination unit, a comprehensive calculation unit and a disaster warning unit, wherein: A data acquisition unit, used to acquire position detection data of each monitoring point in the target slope area at a current detection time point; A coordinate mapping unit, used to determine the three-dimensional coordinates of each monitoring point at the current detection time point according to the position detection data of each monitoring point at the current detection time point; A first determination unit, configured to determine a relative displacement vector and a dynamic strain tensor between any two monitoring points at a current detection time point based on the three-dimensional coordinates of each monitoring point at a current detection time point; The second determination unit is used to retrieve the three-dimensional coordinates of each monitoring point at the previous detection time point, and determine the single-point displacement rate of each monitoring point in the time window from the previous detection time point to the current detection time point based on the three-dimensional coordinates of each monitoring point at the previous detection time point and the three-dimensional coordinates of each monitoring point at the current detection time point; The third determination unit is used to calculate the displacement rate coordination parameter between the two monitoring points corresponding to the current detection time point according to the relative displacement vector between a pair of monitoring points at the current detection time point and the single-point displacement rate of the two monitoring points within the time window; The fourth determination unit is used to determine the angle parameter of the line connecting the two monitoring points according to the relative displacement vector between a pair of monitoring points at the current detection time point, and calculate the spatial correlation parameter between the two monitoring points at the current detection time point according to the angle parameter of the line connecting the two monitoring points; A comprehensive calculation unit is used to calculate the comprehensive early warning index of the target slope area at the current detection time point according to the relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points at the current detection time point; The disaster warning unit is used to provide disaster warning for the target slope area based on the comprehensive warning index of the target slope area at the current detection time point.

[0020] In a third aspect, a slope disaster early warning system is provided, comprising: A memory for storing instructions; A processor is used to read the instructions stored in the memory and execute any one of the methods described in the first aspect according to the instructions.

[0021] In a fourth aspect, a computer-readable storage medium is provided, wherein instructions are stored on the computer-readable storage medium, and when the instructions are executed on a computer, the computer executes any one of the methods described in the first aspect. In addition, a computer program product is provided, and when the computer program product is executed on a computer, the computer executes any one of the methods described in the first aspect.

[0022] Beneficial effects: The present invention obtains the position detection data of each monitoring point in the target slope area for multi-dimensional data analysis, determines the relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points, and then uses the relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points to calculate the comprehensive early warning index of the target slope area at the current detection time point, and then performs disaster early warning of the target slope area according to the comprehensive early warning index to achieve efficient and intelligent slope disaster early warning. The present invention can accurately grasp the deformation degree of the slope by performing multi-dimensional reference point data analysis on the slope deformation process, realize dynamic monitoring and evaluation of the slope condition, improve the accuracy of slope disaster monitoring and early warning response efficiency, and enhance the slope disaster monitoring and early warning capabilities. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0024] Figure 1 Schematic diagram of the steps of the method in Example 1 of the present invention; Figure 2 This is a schematic diagram of the structure of the system in Example 2 of the present invention; Figure 3 This is a schematic diagram of the hardware system in Example 3 of the present invention. DETAILED DESCRIPTION

[0025] It should be noted that the description of these embodiments is used to help understand the present invention, but does not constitute a limitation of the present invention. The specific structures and functional details disclosed herein are only used to describe the exemplary embodiments of the present invention. However, the present invention can be embodied in many alternative forms, and it should not be understood that the present invention is limited to the embodiments set forth herein.

[0026] It should be understood that unless otherwise clearly specified and limited, the corresponding terms should be understood in a broad sense. For example, "connection" can be a fixed connection, a detachable connection, or an integral connection; it can be an electrical connection, a direct connection, an indirect connection through an intermediate medium, or the internal connection of two elements. For ordinary technicians in this field, the specific meanings of the above terms in the embodiments can be understood according to specific circumstances.

[0027] In the following description, certain details are provided to facilitate a complete understanding of the example embodiments. However, it will be appreciated by those of ordinary skill in the art that the example embodiments may be implemented without these certain details. For example, devices may be shown in block diagrams to avoid obscuring the examples with unnecessary details. In other embodiments, well-known processes, structures, and techniques may not be shown in unnecessary detail to avoid obscuring the embodiments.

[0028] Embodiment 1: This embodiment provides a slope disaster early warning method, which can be applied to corresponding early warning monitoring servers, such as Figure 1 As shown, the method comprises the following steps: S1. Obtain the position detection data of each monitoring point in the target slope area at the current detection time point.

[0029] In specific implementation, a number of monitoring points (such as reference anchor points) can be preset in the target slope area so that the monitoring points are spread throughout the target slope area. Then, the position detection data of each monitoring point in the target slope area can be collected at the current detection time point. For example, the position detection data of the corresponding monitoring point can be collected through the detection instruments pre-buried at each monitoring point, or the position detection data of each monitoring point can be collected by using unmanned aerial vehicles carrying optical detection instruments, radar detection instruments or infrared detection instruments for high-altitude detection.

[0030] S2. Determine the three-dimensional coordinates of each monitoring point at the current detection time point based on the position detection data of each monitoring point at the current detection time point.

[0031] In specific implementation, after obtaining the position detection data of each monitoring point in the target slope area at the current detection time point, the server can project the position detection data of each monitoring point at the current detection time point into a fixed three-dimensional coordinate system (the three-dimensional coordinate system uses the initial position of a set reference point as the coordinate origin) to obtain the three-dimensional coordinates of each monitoring point at the current detection time point in the set three-dimensional coordinate system.

[0032] S3. Determine the relative displacement vector and dynamic strain tensor between any two monitoring points at the current detection time point based on the three-dimensional coordinates of each monitoring point at the current detection time point.

[0033] In specific implementation, the server substitutes the three-dimensional coordinates of any two monitoring points at the current detection time point into the preset relative displacement vector formula for calculation to obtain the relative displacement vector between any two monitoring points at the current detection time point. The relative displacement vector formula is:

[0034] Where t represents the current detection time point, i and j are the numbers of any two monitoring points, and i≠j, ΔD ij (t) represents the relative displacement vector between monitoring points i and j at the current detection time point t, x i (t)、y i (t) and z i (t) is the three-dimensional coordinate of monitoring point i at the current detection time point, x j (t)、y j (t) and z j (t) is the three-dimensional coordinate of monitoring point j at the current detection time point; Substitute the relative displacement vector between any two monitoring points at the current detection time point into the preset dynamic strain tensor formula for calculation to obtain the dynamic strain tensor between any two monitoring points at the current detection time point. The dynamic strain tensor formula is:

[0035] Among them, ε ij (t) represents the dynamic strain tensor between monitoring points i and j at the current detection time point t, L ij is the initial distance between the given monitoring points i and j, and ‖‖ represents the vector modulo operation.

[0036] S4. Retrieve the three-dimensional coordinates of each monitoring point at the previous detection time point, and determine the single-point displacement rate of each monitoring point within the time window from the previous detection time point to the current detection time point based on the three-dimensional coordinates of each monitoring point at the previous detection time point and the three-dimensional coordinates of each monitoring point at the current detection time point.

[0037] In specific implementation, the server can retrieve the pre-stored three-dimensional coordinates of each monitoring point at the previous detection time point. Then, the three-dimensional coordinates of a monitoring point at the previous detection time point and the three-dimensional coordinates of the current detection time point are substituted into the single-point displacement rate formula for calculation to obtain the single-point displacement rate of the monitoring point in the time window from the previous detection time point to the current detection time point. The single-point displacement rate formula is:

[0038] Among them, Δt represents the time interval from the previous detection time point to the current detection time point, t-Δt represents the previous detection time point, and v i (Δt) represents the single-point displacement rate of monitoring point i within the time window from the previous detection time point to the current detection time point.

[0039] S5. Calculate the displacement rate coordination parameter between the two monitoring points corresponding to the current detection time point according to the relative displacement vector between a pair of monitoring points at the current detection time point and the single-point displacement rate of the two monitoring points within the time window.

[0040] In specific implementation, the server substitutes the relative displacement vector between two monitoring points at the current detection time point and the single-point displacement rate of the two monitoring points in the time window into the preset displacement rate coordination parameter formula for calculation, and obtains the displacement rate coordination parameter between the two monitoring points at the current detection time point. The displacement rate coordination parameter formula is:

[0041] Among them, δ ij (t) represents the displacement rate coordination parameter between monitoring point i and monitoring point j at the current detection time point, v i (Δt) represents the displacement rate of a single point at monitoring point i in the time window from the previous detection time point to the current detection time point, v j (Δt) represents the single-point displacement rate of monitoring point j within the time window from the previous detection time point to the current detection time point.

[0042] S6. Determine the angle parameter of the line connecting the two monitoring points according to the relative displacement vector between a pair of monitoring points at the current detection time point, and calculate the spatial correlation parameter between the two monitoring points at the current detection time point according to the angle parameter of the line connecting the two monitoring points.

[0043] In specific implementation, the server uses the angle between the relative displacement vector between a pair of monitoring points and the set main structural surface as the angle parameter of the line connecting the two monitoring points; Substitute the angle parameter of the line connecting the two corresponding monitoring points into the preset spatial correlation parameter formula for calculation to obtain the spatial correlation parameter between the two corresponding monitoring points. The spatial correlation parameter formula is:

[0044] Among them, i and j are the numbers of the corresponding two monitoring points, ω ij is the spatial correlation parameter between monitoring point i and monitoring point j, L ij is the initial distance between the given monitoring point i and the monitoring point j, θ ij is the angle parameter of the line connecting monitoring point i and monitoring point j, θ0 is the set reference angle constant, L0 is the set characteristic attenuation distance, and exp represents the natural exponential function.

[0045] S7. Calculate the comprehensive early warning index of the target slope area at the current detection time point based on the relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points at the current detection time point.

[0046] In specific implementation, the server substitutes the relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points at the current detection time point into the preset comprehensive warning index formula for calculation, and obtains the comprehensive warning index of the target slope area at the current detection time point. The comprehensive warning index formula is:

[0047] Among them, S(t) is the comprehensive warning index of the target slope area at the current detection time point, i and j are the numbers of the corresponding two monitoring points, represents the combination of monitoring point i and monitoring point j, N is the number of all possible combinations of monitoring points, ω ij is the spatial correlation parameter between monitoring point i and monitoring point j, ΔD ij (t) is the relative displacement vector between monitoring points i and j at the current detection time point t, ‖‖ represents the vector modulus operation, ε ij (t) is the dynamic strain tensor between monitoring points i and j at the current detection time point t, δ ij (t) is the displacement rate coordination parameter between monitoring point i and monitoring point j at the current detection time point t, α is the set first adjustment coefficient, β is the set second adjustment coefficient, and γ is the set third adjustment coefficient.

[0048] S8. Conduct disaster warning for the target slope area based on the comprehensive warning index of the target slope area at the current detection time point.

[0049] In specific implementation, when the server determines that the comprehensive warning index of the target slope area at the current detection time point exceeds the set first threshold, the server outputs yellow disaster warning information of the target slope area; when the server determines that the comprehensive warning index of the target slope area at the current detection time point exceeds the set second threshold, the server outputs red disaster warning information of the target slope area, and the second threshold is greater than the first threshold; so as to achieve different levels of warnings according to actual conditions.

[0050] This method can accurately grasp the deformation degree of the slope by performing multi-dimensional reference point data analysis on the slope deformation process, realize dynamic monitoring and evaluation of the slope condition, improve the accuracy of slope disaster monitoring and early warning response efficiency, and enhance the slope disaster monitoring and early warning capabilities.

[0051] Embodiment 2: This embodiment provides a slope disaster early warning system. Figure 2 As shown, it includes a data acquisition unit, a coordinate mapping unit, a first determination unit, a second determination unit, a third determination unit, a fourth determination unit, a comprehensive calculation unit and a disaster warning unit, wherein: A data acquisition unit, used to acquire position detection data of each monitoring point in the target slope area at a current detection time point; A coordinate mapping unit, used to determine the three-dimensional coordinates of each monitoring point at the current detection time point according to the position detection data of each monitoring point at the current detection time point; A first determination unit, configured to determine a relative displacement vector and a dynamic strain tensor between any two monitoring points at a current detection time point based on the three-dimensional coordinates of each monitoring point at a current detection time point; The second determination unit is used to retrieve the three-dimensional coordinates of each monitoring point at the previous detection time point, and determine the single-point displacement rate of each monitoring point in the time window from the previous detection time point to the current detection time point based on the three-dimensional coordinates of each monitoring point at the previous detection time point and the three-dimensional coordinates of each monitoring point at the current detection time point; The third determination unit is used to calculate the displacement rate coordination parameter between the two monitoring points corresponding to the current detection time point according to the relative displacement vector between a pair of monitoring points at the current detection time point and the single-point displacement rate of the two monitoring points within the time window; The fourth determination unit is used to determine the angle parameter of the line connecting the two monitoring points according to the relative displacement vector between a pair of monitoring points at the current detection time point, and calculate the spatial correlation parameter between the two monitoring points at the current detection time point according to the angle parameter of the line connecting the two monitoring points; A comprehensive calculation unit is used to calculate the comprehensive early warning index of the target slope area at the current detection time point according to the relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points at the current detection time point; The disaster warning unit is used to provide disaster warning for the target slope area based on the comprehensive warning index of the target slope area at the current detection time point.

[0052] Embodiment 3: This embodiment provides a slope disaster early warning system. Figure 3 As shown, at the hardware level, it includes: Data interface, used to establish data connection between the processor and the external data acquisition terminal; A memory for storing instructions; The processor is used to read the instructions stored in the memory and execute the slope disaster warning method in Example 1 according to the instructions.

[0053] Optionally, the system further includes an internal bus, through which the processor, the memory and the data interface can be interconnected, and the internal bus can be an ISA (Industry Standard Architecture) bus, a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc.

[0054] The memory may include, but is not limited to, random access memory (RAM), read only memory (ROM), flash memory, first input first output (FIFO) and / or first in last out (FILO), etc. The processor may be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it may also be a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field programmable gate array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components.

[0055] Embodiment 4: This embodiment provides a computer-readable storage medium, on which instructions are stored, and when the instructions are executed on a computer, the computer executes the slope disaster early warning method in Embodiment 1. The computer-readable storage medium refers to a carrier for storing data, which may include but is not limited to a floppy disk, an optical disk, a hard disk, a flash memory, a USB flash drive, and / or a memory stick, etc. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices.

[0056] This embodiment also provides a computer program product, which, when executed on a computer, executes the slope disaster early warning method in Embodiment 1. The computer may be a general-purpose computer, a special-purpose computer, a computer network or other programmable devices.

[0057] Finally, it should be noted that the above description is only a preferred embodiment of the present invention and is not intended to limit the protection scope of the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A slope disaster early warning method, characterized in that: include: Obtain the position detection data of each monitoring point in the target slope area at the current detection time point; Determine the three-dimensional coordinates of each monitoring point at the current detection time point according to the position detection data of each monitoring point at the current detection time point; Determine the relative displacement vector and dynamic strain tensor between any two monitoring points at the current detection time point based on the three-dimensional coordinates of each monitoring point at the current detection time point; Retrieving the three-dimensional coordinates of each monitoring point at the previous detection time point, and determining the single-point displacement rate of each monitoring point within the time window from the previous detection time point to the current detection time point based on the three-dimensional coordinates of each monitoring point at the previous detection time point and the three-dimensional coordinates of each monitoring point at the current detection time point; According to the relative displacement vector between a pair of monitoring points at the current detection time point and the single-point displacement rate of the two monitoring points in the time window, the displacement rate coordination parameter between the two monitoring points at the current detection time point is calculated; Determine the angle parameter of the line connecting the two monitoring points according to the relative displacement vector between a pair of monitoring points at the current detection time point, and calculate the spatial correlation parameter between the two monitoring points at the current detection time point according to the angle parameter of the line connecting the two monitoring points; According to the relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points at the current detection time point, the comprehensive early warning index of the target slope area at the current detection time point is calculated; Disaster warning for the target slope area is carried out based on the comprehensive warning index of the target slope area at the current detection time point.

2. A slope disaster early warning method according to claim 1, characterized in that: Determining the three-dimensional coordinates of each monitoring point at the current detection time point according to the position detection data of each monitoring point at the current detection time point includes: The position detection data of each monitoring point at the current detection time point is projected into the set three-dimensional coordinate system to obtain the three-dimensional coordinates of each monitoring point at the current detection time point in the set three-dimensional coordinate system.

3. A slope disaster early warning method according to claim 1, characterized in that: The method of determining the relative displacement vector and the dynamic strain tensor between any two monitoring points at the current detection time point based on the three-dimensional coordinates of each monitoring point at the current detection time point includes: Substitute the three-dimensional coordinates of any two monitoring points at the current detection time point into the preset relative displacement vector formula for calculation to obtain the relative displacement vector between any two monitoring points at the current detection time point. The relative displacement vector formula is: Where t represents the current detection time point, i and j are the numbers of any two monitoring points, and i≠j, ΔD ij (t) represents the relative displacement vector between monitoring points i and j at the current detection time point t, x i (t)、y i (t) and z i (t) is the three-dimensional coordinate of monitoring point i at the current detection time point, x j (t)、y j (t) and z j (t) is the three-dimensional coordinate of monitoring point j at the current detection time point; Substitute the relative displacement vector between any two monitoring points at the current detection time point into the preset dynamic strain tensor formula for calculation to obtain the dynamic strain tensor between any two monitoring points at the current detection time point. The dynamic strain tensor formula is: Among them, ε ij (t) represents the dynamic strain tensor between monitoring points i and j at the current detection time point t, L ij is the initial distance between the given monitoring points i and j, and ‖‖ represents the vector modulo operation.

4. A slope disaster early warning method according to claim 3, characterized in that: The method of determining the single-point displacement rate of each monitoring point within a time window from the previous detection time point to the current detection time point based on the three-dimensional coordinates of each monitoring point at the previous detection time point and the three-dimensional coordinates of each monitoring point at the current detection time point includes: Substitute the three-dimensional coordinates of a monitoring point at the previous detection time point and the three-dimensional coordinates of the current detection time point into the single-point displacement rate formula for calculation to obtain the single-point displacement rate of the monitoring point in the time window from the previous detection time point to the current detection time point. The single-point displacement rate formula is: Among them, Δt represents the time interval from the previous detection time point to the current detection time point, t-Δt represents the previous detection time point, and v i (Δt) represents the single-point displacement rate of monitoring point i within the time window from the previous detection time point to the current detection time point.

5. A slope disaster early warning method according to claim 4, characterized in that: The method of calculating the displacement rate coordination parameter between the two monitoring points corresponding to the current detection time point according to the relative displacement vector between a pair of monitoring points at the current detection time point and the single-point displacement rate of the two monitoring points in the time window includes: Substitute the relative displacement vector between two monitoring points at the current detection time point and the single-point displacement rate of the two monitoring points in the time window into the preset displacement rate coordination parameter formula for calculation to obtain the displacement rate coordination parameter between the two monitoring points at the current detection time point. The displacement rate coordination parameter formula is: Among them, δ ij (t) represents the displacement rate coordination parameter between monitoring point i and monitoring point j at the current detection time point, v i (Δt) represents the displacement rate of a single point at monitoring point i in the time window from the previous detection time point to the current detection time point, v j (Δt) represents the single-point displacement rate of monitoring point j within the time window from the previous detection time point to the current detection time point.

6. A slope disaster early warning method according to claim 1, characterized in that: The method of determining the angle parameter of the line connecting the two monitoring points according to the relative displacement vector between a pair of monitoring points at the current detection time point, and calculating the spatial correlation parameter between the two monitoring points at the current detection time point according to the angle parameter of the line connecting the two monitoring points, includes: The angle between the relative displacement vector of a pair of monitoring points and the set main structural surface is used as the angle parameter of the line connecting the two monitoring points; Substitute the angle parameter of the line connecting the two corresponding monitoring points into the preset spatial correlation parameter formula for calculation to obtain the spatial correlation parameter between the two corresponding monitoring points. The spatial correlation parameter formula is: Among them, i and j are the numbers of the corresponding two monitoring points, ω ij is the spatial correlation parameter between monitoring point i and monitoring point j, L ij is the initial distance between the given monitoring point i and the monitoring point j, θ ij is the angle parameter of the line connecting monitoring point i and monitoring point j, θ0 is the set reference angle constant, L0 is the set characteristic attenuation distance, and exp represents the natural exponential function.

7. A slope disaster early warning method according to claim 1, characterized in that: The method of calculating the comprehensive early warning index of the target slope area at the current detection time point according to the relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points at the current detection time point includes: The relative displacement vector, dynamic strain tensor, displacement rate coordination parameter and spatial correlation parameter between each pair of monitoring points at the current detection time point are substituted into the preset comprehensive warning index formula for calculation to obtain the comprehensive warning index of the target slope area at the current detection time point. The comprehensive warning index formula is: Among them, S(t) is the comprehensive warning index of the target slope area at the current detection time point, i and j are the numbers of the corresponding two monitoring points, represents the combination of monitoring point i and monitoring point j, N is the number of all possible combinations of monitoring points, ω ij is the spatial correlation parameter between monitoring point i and monitoring point j, ΔD ij (t) is the relative displacement vector between monitoring points i and j at the current detection time point t, ‖‖ represents the vector modulus operation, ε ij (t) is the dynamic strain tensor between monitoring points i and j at the current detection time point t, δ ij (t) is the displacement rate coordination parameter between monitoring point i and monitoring point j at the current detection time point t, α is the set first adjustment coefficient, β is the set second adjustment coefficient, and γ is the set third adjustment coefficient.

8. A slope disaster early warning method according to claim 1, characterized in that: The disaster warning of the target slope area based on the comprehensive warning index of the target slope area at the current detection time point includes: When it is determined that the comprehensive warning index of the target slope area at the current detection time point exceeds the set first threshold, the yellow disaster warning information of the target slope area is output; When it is determined that the comprehensive warning index of the target slope area at the current detection time point exceeds a set second threshold, red disaster warning information of the target slope area is output, and the second threshold is greater than the first threshold.

9. A slope disaster early warning system, characterized in that: include: A memory for storing instructions; A processor is used to read the instructions stored in the memory and execute the slope disaster warning method described in any one of claims 1-8 according to the instructions.

Images