Intelligent side slope monitoring method and system based on digital twinning and storage medium

Abstract

The invention discloses a side slope intelligent monitoring method and system based on digital twinning and a storage medium, and belongs to the technical field of digital twinning, and the method comprises the steps: obtaining point cloud data of a side slope surface, and creating a digital twinning model; dividing the digital twin model to generate a plurality of first grids, and calculating displacement change values at the first grids within the target duration; defining the first grid of which the displacement change value is greater than a first threshold value as a first area, and dividing the target duration into a plurality of first durations with unequal time; setting a plurality of rainfall modes, screening out the optimal mode of each first duration, calculating the corresponding representative rainfall capacity, constructing a prediction model, and predicting and generating the representative rainfall capacity in a future time period according to the optimal mode and the representative rainfall capacity of the historical rainfall data in the first duration; based on the representative rainfall and the pore water pressure of the soil, the safety coefficient of the side slope in the future time period of the first area is predicted and generated, and the accuracy of side slope stability monitoring is improved.

Description

Technical Field

[0001] The present invention belongs to the field of digital twin technology, and specifically relates to a slope intelligent monitoring method, system and storage medium based on digital twin. Background Art

[0002] In slope engineering, intelligent monitoring using digital twin technology is an emerging approach that combines advanced sensor technology, data acquisition, data processing, and numerical simulation technology to improve the accuracy and real-time performance of slope stability assessment.

[0003] For example, Chinese patent application publication number CN117763969A proposes a digital twin-based slope monitoring method, apparatus, device, and storage medium. The method includes: obtaining monitoring data sets from multiple types of sensors installed on the slope to be monitored; performing sensitivity analysis on the monitoring data in the monitoring data sets, and selecting sensitive data sets that meet preset sensitivities; inputting the sensitive data sets into a pre-built digital twin model to invert and obtain an inversion data set for the slope to be monitored; fusing the inversion data sets to obtain a comprehensive monitoring information curve for the slope to be monitored; and determining the state of the slope to be monitored based on the inclination of the comprehensive monitoring information curve. Another example is Chinese patent application publication number CN117828732A, which discloses a digital twin-based slope stability determination method, system, medium, and terminal. This method relates to the fields of digital twin technology and ionic rare earth minerals. Its primary purpose is to prevent changes in the stress conditions of mine slopes caused by the injection of liquid leachate into the ore body, which in turn affects the stability of the mine slope. The method includes: receiving a determination instruction for the target mine slope stability parameters, and determining that the instruction carries a target time; determining the groundwater level data of the target mine at the target time from the water potential field time distribution data of the target mine; performing a three-dimensional slope stability calculation based on the groundwater level data based on the strength reduction method, and obtaining the mine slope stability prediction result of the target mine at the target time.

[0004] However, the above-mentioned existing technologies for intelligent slope monitoring may not be able to fully integrate the impact of multiple factors on slope stability, such as displacement changes, weather conditions and soil properties, which limits the monitoring system's ability to comprehensively evaluate slope stability. Summary of the Invention

[0005] To solve the above problems, the present invention provides a slope intelligent monitoring method, system and storage medium based on digital twins to solve the problems in the prior art.

[0006] In order to achieve the above-mentioned purpose of the invention, the present invention proposes a slope intelligent monitoring method based on digital twin, comprising: S1: Acquire first point cloud data of the slope surface at a preset time point in the past, and create a digital twin model based on the first point cloud data; S2: Acquire second point cloud data of the slope surface at the current time point, convert the first point cloud data, the second point cloud data, and the digital twin model into the same coordinate system, mesh the digital twin model based on a preset size to generate a plurality of first grids, define the time period between the preset time point and the current time point as a target time period, calculate the displacement change value at each first grid within the target time period, and display it in the digital twin model; S3: defining the first grid where the displacement change value is greater than a first threshold as a first area, obtaining the actual rainfall in the first area within the target duration, establishing a first coordinate system with time as the horizontal axis and the actual rainfall as the vertical axis, plotting a first variation curve of the actual rainfall over time in the first coordinate system, obtaining extreme points of the first variation curve, defining the time point corresponding to each extreme point as a first time point, and dividing the target duration into multiple first durations of varying lengths based on the first time points; S4: setting multiple rainfall patterns, each rainfall pattern corresponding to multiple historical rainfall data, selecting a rainfall pattern from the rainfall patterns as an optimal pattern based on the actual rainfall within the first time period, obtaining an optimal time point within each first time period based on the optimal pattern and the historical rainfall data, using the actual rainfall corresponding to the optimal time point as a representative rainfall within the first time period, building a prediction model based on the historical rainfall data, predicting and generating the optimal pattern and the representative rainfall in a future time period based on the optimal pattern and the representative rainfall in each first time period, and defining the representative rainfall in the future time period as predicted rainfall; S5: Obtain historical rainfall within the target time period, obtain a first pore water pressure distribution of the soil in the first area based on the historical rainfall, predict and generate a first safety factor of the first area in the future time period based on the predicted rainfall and the first pore water pressure distribution, take the average of all the first safety factors of the slope as the final safety factor of the slope, and display it in the digital twin model.

[0007] Furthermore, calculating the displacement change value at each of the first grids includes the following steps: defining a vertex of the first mesh in the digital twin model as a first node, and defining the three-dimensional position coordinates after conversion of the first point cloud data and the second point cloud data as first position coordinates and second position coordinates, respectively; Obtain a first number of the first point cloud data and a second number of the second point cloud data within the first grid. If both the first number and the second number are less than or equal to 1, define a vector formed by taking the first node as the starting point and the first position coordinate of the first point cloud data as the end point as a first vector, and define a vector formed by taking the first node as the starting point and the second position coordinate of the second point cloud data as the end point as a second vector. If the first number is greater than 1, then calculating the barycentric coordinates of all the first point cloud data in the first grid, and using the first node as the starting point and the barycentric coordinates of the first point cloud data as the end point to form a vector as the first vector; if the second number is greater than 1, then calculating the barycentric coordinates of all the second point cloud data in the first grid, and using the first node as the starting point and the barycentric coordinates of the second point cloud data as the end point to form a vector as the second vector; Get the unit vector of the first grid in the normal direction , calculate the first component of the first vector and the second vector in the normal direction based on the first formula and the second formula respectively and the second component , the first formula is: , the second formula is: ,in, and The first vector and the second vector are respectively defined, and the difference between the first component and the second component is defined as the displacement change value of the first grid.

[0008] Furthermore, selecting the optimal pattern from the rainfall patterns includes the following steps: Calculate the average value of the actual rainfall of the first change curve at multiple first time points. If the first duration is greater than the second threshold, the average value is greater than the third threshold, and the first time point is a maximum point, match the optimal mode to the first mode. If the first duration is greater than the second threshold, the average value is less than the third threshold, and the first time point is a minimum point, match the optimal mode to the second mode. If the first duration is less than the second threshold, the average value is greater than the third threshold, and the first time point is a maximum point, match the optimal mode to the third mode. If the first duration is less than the second threshold, the average value is less than the third threshold, and the first time point is a minimum point, match the optimal mode to the fourth mode.

[0009] Furthermore, obtaining the optimal time point within each of the first durations includes the following steps: Obtain the average rainfall of multiple pieces of historical rainfall data corresponding to the optimal mode at the first time point, calculate the difference between the actual rainfall at each of the first time points and the corresponding average rainfall, and define the time point when the difference is less than or equal to a fourth threshold as the optimal time point.

[0010] Furthermore, building a prediction model includes the following steps: The time period in which each piece of the historical rainfall data is located is defined as a historical time length, and the historical time length is divided into multiple first time lengths. The corresponding optimal mode is matched for each first time length, and the representative rainfall under the optimal mode is calculated. A division point is set within the historical time length, and the optimal mode and the representative rainfall of each first time length in the time period before the division point are used as input features. After the division point, one of the first time lengths is selected, and the corresponding optimal mode and the representative rainfall are used as output features. This step is repeated to process each piece of the historical rainfall data, and a neural network model is constructed. The neural network model is trained based on the input features and the corresponding output features to obtain the prediction model.

[0011] Furthermore, generating a first safety factor of the first area at the future time period includes the following steps: If the predicted rainfall in the future time period is greater than the fifth threshold, the first pore water pressure distribution is used as the initial value and the predicted rainfall is used as the boundary condition. The second pore water pressure distribution of the soil in the future time period is predicted based on the one-dimensional permeability model. The second pore water pressure of the soil in the first area is obtained based on the second pore water pressure distribution. The effective stress of the soil in the first area is calculated based on the second pore water pressure. The effective stress is the difference between the total stress of the soil in the first area and the second pore water pressure. The shear strength of the soil under the second pore water pressure is obtained, and the ratio of the shear strength to the effective stress is used as the first safety factor of the first area in the future time period.

[0012] The present invention also provides a digital twin-based intelligent slope monitoring system, which is used to implement the above-mentioned method. The system mainly includes: A creation module is used to obtain first point cloud data of the slope surface at a preset time point in the past, and create a digital twin model based on the first point cloud data; a calculation module, configured to obtain second point cloud data of the slope surface at a current time point, convert the first point cloud data, the second point cloud data, and the digital twin model into the same coordinate system, mesh the digital twin model based on a preset size to generate a plurality of first grids, define the duration between the preset time point and the current time point as a target duration, calculate a displacement change value at each of the first grids within the target duration, and display the calculated displacement change value in the digital twin model; a partitioning module, configured to define the first grid having a displacement change value greater than a first threshold as a first region, obtain actual rainfall in the first region within the target duration, establish a first coordinate system with time as the horizontal axis and the actual rainfall as the vertical axis, plot a first variation curve of the actual rainfall over time in the first coordinate system, obtain extreme value points of the first variation curve, define a time point corresponding to each extreme value point as a first time point, and divide the target duration into a plurality of first durations of varying lengths based on the first time points; a prediction module for setting a plurality of rainfall patterns, each rainfall pattern corresponding to a plurality of historical rainfall data, selecting a rainfall pattern from the rainfall patterns as an optimal pattern based on the actual rainfall within the first time period, obtaining an optimal time point within each first time period based on the optimal pattern and the historical rainfall data, using the actual rainfall corresponding to the optimal time point as a representative rainfall within the first time period, constructing a prediction model based on the historical rainfall data, predicting and generating the optimal pattern and the representative rainfall within a future time period based on the optimal pattern and the representative rainfall within each first time period, and defining the representative rainfall within the future time period as the predicted rainfall; A monitoring module is used to obtain historical rainfall within the target time period, obtain a first pore water pressure distribution of the soil in the first area based on the historical rainfall, and predict and generate a first safety factor for the first area in the future time period based on the predicted rainfall and the first pore water pressure distribution. The average of all the first safety factors of the slope is used as the final safety factor of the slope and is displayed in the digital twin model.

[0013] The present invention also provides a computer storage medium storing program instructions, wherein when the program instructions are executed, the device where the computer storage medium is located is controlled to execute the above-mentioned method.

[0014] Compared with the prior art, the beneficial effects of the present invention are at least as follows: By acquiring the first point cloud data at a preset time point in the past, the present invention can create an accurate digital twin model of the slope. This model can reflect the initial state of the slope in detail, providing an accurate basis for subsequent analysis and prediction; the digital twin model is gridded to generate multiple first grids, and by calculating the displacement change of the slope in each first grid from the preset time point to the current time point, the displacement change of the slope in each local area is analyzed in more detail; the first grid with a displacement change value greater than a preset threshold is defined as the first area. The present invention focuses on the first area with a large landslide risk or strong instability, and analyzes and predicts the rainfall in the first area.

[0015] The present invention also divides the target time into multiple first time periods of varying lengths, matches a corresponding rainfall pattern for each first time period, and each rainfall pattern includes multiple historical rainfall data. A prediction model is constructed based on the rainfall data, and the prediction model is a neural network model. The historical time period in which each historical rainfall data is located is also divided into multiple first time periods, and the optimal pattern is assigned to each first time period, and the representative rainfall under the optimal pattern is calculated. The neural network model learns the rainfall change trend under each rainfall pattern corresponding to each first time period in the historical time period, and predicts and generates the optimal rainfall pattern and predicted rainfall in the future time period, thereby improving the accuracy of predicting future rainfall trends; the first pore water pressure distribution of the slope soil in the first area of ​​the future time period is determined based on the predicted rainfall in the future time period and the historical rainfall in the target time period, and the first safety factor of the slope is obtained based on the first pore water pressure distribution. By predicting the slope stability under future rainfall conditions, preventive measures can be taken in advance to reduce disasters such as landslides and mudslides. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] Figure 1 This is a flowchart of the steps of the digital twin-based intelligent slope monitoring method of the present invention; Figure 2 is the first grid diagram of the present invention; Figure 3 is a first change curve diagram of the present invention; Figure 4 This is the system structure diagram of the digital twin-based intelligent slope monitoring system of the present invention. DETAILED DESCRIPTION

[0017] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0018] It is understood that the terms "first," "second," etc., used herein may be used to describe various elements, but unless otherwise specified, these elements are not limited by these terms. These terms are only used to distinguish a first element from another element. For example, a first xx script may be referred to as a second xx script, and similarly, a second xx script may be referred to as a first xx script without departing from the scope of this application.

[0019] like Figure 1 As shown in FIG, a digital twin-based intelligent slope monitoring method includes: S1: Acquire first point cloud data of the slope surface at a preset time point in the past, and create a digital twin model based on the first point cloud data.

[0020] Specifically, in order to clearly monitor the displacement changes and stability of the slope, it is first necessary to construct a digital twin model of the slope based on the three-dimensional point cloud data of the slope to obtain the entire regional range of the entire slope, and then analyze the displacement changes of the slope; in this embodiment, the slope is not limited to natural slopes, clay slopes, coastal or riverside slopes, etc.; a three-dimensional laser scanner or a scanning radar of a drone is used to obtain point cloud data of the slope surface at a preset time point in the past (0:00 am), and the point cloud data is preprocessed, including noise removal, filtering and data alignment, etc., and the processed point cloud data is converted into a three-dimensional model using three-dimensional modeling software such as AutoCAD, and the key geological features of the slope, such as structural surfaces, cracks, slopes, etc., are extracted from the three-dimensional model. Finally, these multi-source data are integrated to construct a digital twin model.

[0021] S2: Obtain the second point cloud data of the slope surface at the current time point, convert the first point cloud data, the second point cloud data and the digital twin model into the same coordinate system, mesh the digital twin model based on the preset size to generate multiple first grids, define the time length between the preset time point and the current time point as the target time length, calculate the displacement change value at each first grid within the target time length, and display it in the digital twin model.

[0022] Specifically, coordinate transformation operations such as translation, rotation, and scaling are performed on the first point cloud data and the second point cloud data on the slope surface, and the first point cloud data and the second point cloud data are aligned with the digital twin model to ensure that they are in the same coordinate system. In order to make the calculation of the slope displacement change more accurate, the model is gridded according to a predetermined size, such as 2m*2m. The slope can be divided into smaller parts through grid division. Each grid can represent a small block of the slope, and its displacement can be calculated independently. The displacement change value of each first grid within the target time length is calculated and finally displayed in the digital twin model. The specific calculation method will be introduced later.

[0023] S3: Define the first grid whose displacement change value is greater than the first threshold as the first area, obtain the actual rainfall in the first area within the target time length, establish a first coordinate system with time as the horizontal axis and the actual rainfall as the vertical axis, plot the first change curve of the actual rainfall over time in the first coordinate system, obtain the extreme value points of the first change curve, define the time point corresponding to each extreme value point as the first time point, and divide the target time length into multiple first time lengths of unequal time based on the first time point.

[0024] Specifically, the first area is a grid area where the slope displacement change value is large, the instability is strong, and landslides are prone to occur. The first threshold is set based on the actual situation. The important factors causing landslides in the first area include the rainfall in the area. Heavy rainfall or continuous rainfall will increase the pore water pressure of the soil and reduce the shear strength of the soil, thereby increasing the risk of landslides. This embodiment mainly analyzes and predicts the stability of the slope based on precipitation. The purpose of setting the first change curve is to intuitively show the changes in the actual rainfall in the first area over time during the duration of rainfall. Figure 3 , obtain extreme points in the first change curve, such as s1, s2, s3, s4, s5, etc., and divide 1 day into multiple first time periods of varying lengths, including 0-4h, 4-8h, 12-20h, 20-22h, and 22-24h.

[0025] S4: Set multiple rainfall patterns, each rainfall pattern corresponds to multiple historical rainfall data, select a rainfall pattern from the rainfall patterns as the optimal pattern based on the actual rainfall in the first time period, obtain the optimal time point in each first time period based on the optimal pattern and the historical rainfall data, and use the actual rainfall corresponding to the optimal time point as the representative rainfall in the first time period. Build a prediction model based on the historical rainfall data, and predict and generate the optimal pattern and representative rainfall in the future time period based on the optimal pattern and representative rainfall in each first time period, and define the representative rainfall in the future time period as the predicted rainfall.

[0026] Specifically, according to the rainfall change trend of the slope in the historical time period, a variety of rainfall patterns are obtained, including but not limited to long-term heavy rainfall pattern, long-term weak rainfall pattern, short-term heavy rainfall pattern and short-term weak rainfall pattern. The rainfall pattern that best matches the rainfall pattern in each first time period is screened out from the multiple rainfall patterns as the optimal pattern. In this way, the change of the rainfall pattern in the target time period can be obtained, and a prediction model is constructed. The optimal pattern of each first time period in the target time period and the average rainfall of the optimal pattern are input into the prediction model to obtain the predicted rainfall in the future time period. This step splits the target time period into multiple first time periods of different time periods, so as to better analyze the rainfall in each stage within the target time period in a more detailed manner. The rainfall change trend of each unequal time period is different, and the corresponding rainfall pattern is matched for each unequal first time period based on the change trend.

[0027] By matching the corresponding rainfall patterns to multiple first time periods within the target time period, the corresponding representative rainfall is calculated. Based on these data, a prediction model is used to predict the rainfall pattern and representative rainfall in future time periods, so that countermeasures can be taken in advance.

[0028] Step S5: Obtain historical rainfall within the target time period, obtain the first pore water pressure distribution of the soil in the first area based on the historical rainfall, and predict and generate the first safety factor of the first area in the future time period based on the predicted rainfall and the first pore water pressure distribution. The average of all the first safety factors of the slope is used as the final safety factor of the slope and displayed in the digital twin model.

[0029] Specifically, it is defined that if there is no rainfall for a long time before this rainfall, the soil is in a dry state, otherwise it is in a wet state, and the first absorption capacity of the soil in the first area is obtained. The first absorption capacity is the maximum value of water absorbed by the soil in a dry state. The standard for defining the dry state is based on the one-dimensional infiltration model combined with the first absorption capacity to obtain the initial pore water pressure distribution of the first area in a dry state, and the initial pore water pressure distribution is defined as the first pore water pressure distribution. If the historical rainfall is greater than or equal to a preset threshold, the first absorption capacity of the soil in the first area in the historical time period is obtained, and the initial pore water pressure distribution of the first area is obtained based on the first absorption capacity. With the initial pore water pressure distribution as the initial value and the historical rainfall as the boundary condition, the first pore water pressure distribution of the first area for the target time length is calculated based on the one-dimensional infiltration model.

[0030] Based on the predicted rainfall and the first pore water pressure distribution, a one-dimensional permeability model is used to predict and generate the first safety factor of the first area in the future time period. Each first area has a first safety factor, and the average value of all first safety factors is used as the final safety factor of the slope. If the final safety factor is less than the preset threshold, it is determined that the slope has a landslide hazard. All first areas are located, and the slope of each area is repaired in advance according to the first safety factor of each area.

[0031] By acquiring the first point cloud data at a preset time point in the past, the present invention can create an accurate digital twin model of the slope. This model can reflect the initial state of the slope in detail, providing an accurate basis for subsequent analysis and prediction; the digital twin model is gridded to generate multiple first grids, and by calculating the displacement change of the slope in each first grid from the preset time point to the current time point, the displacement change of the slope in each local area is analyzed in more detail; the first grid with a displacement change value greater than a preset threshold is defined as the first area. The present invention focuses on the first area with a large landslide risk or strong instability, and analyzes and predicts the rainfall in the first area.

[0032] The present invention also divides the target time into multiple first time periods of varying lengths, matches a corresponding rainfall pattern for each first time period, and each rainfall pattern includes multiple historical rainfall data. A prediction model is constructed based on the rainfall data, and the prediction model is a neural network model. The historical time period in which each historical rainfall data is located is also divided into multiple first time periods, and the optimal pattern is assigned to each first time period, and the representative rainfall under the optimal pattern is calculated. The neural network model learns the rainfall change trend under each rainfall pattern corresponding to each first time period in the historical time period, and predicts and generates the optimal rainfall pattern and predicted rainfall in the future time period, thereby improving the accuracy of predicting future rainfall trends; the first pore water pressure distribution of the slope soil in the first area of ​​the future time period is determined based on the predicted rainfall in the future time period and the historical rainfall in the target time period, and the first safety factor of the slope is obtained based on the first pore water pressure distribution. By predicting the slope stability under future rainfall conditions, preventive measures can be taken in advance to reduce disasters such as landslides and mudslides.

[0033] As a preferred technical solution of the present invention, calculating the displacement change value at each first grid includes the following steps: A vertex of the first mesh in the digital twin model is defined as a first node, and three-dimensional position coordinates after conversion of the first point cloud data and the second point cloud data are defined as first position coordinates and second position coordinates, respectively; Obtain a first number of first point cloud data and a second number of second point cloud data within the first grid. If both the first number and the second number are less than or equal to 1, define a vector formed by taking the first node as the starting point and the first position coordinate of the first point cloud data as the end point as the first vector, and define a vector formed by taking the first node as the starting point and the second position coordinate of the second point cloud data as the end point as the second vector. If the first number is greater than 1, the barycentric coordinates of all first point cloud data in the first grid are calculated, and a vector formed by taking the first node as the starting point and the barycentric coordinates of the first point cloud data as the end point is used as the first vector. If the second number is greater than 1, the barycentric coordinates of all second point cloud data in the first grid are calculated, and a vector formed by taking the first node as the starting point and the barycentric coordinates of the second point cloud data as the end point is used as the second vector. Get the unit vector of the first mesh in the normal direction , calculate the first component of the first vector and the second vector in the normal direction based on the first formula and the second formula respectively and the second component , the first formula is: , the second formula is: ,in, and The first and second vectors are defined as the displacement change value of the first grid.

[0034] Specifically, if Figure 2 As shown in the figure, it is the first grid diagram. The K point in the first grid is defined as the first node. The first grid has b1, b2, c1, c2 and c3, where b1 and b2 are the first point cloud data, and c1, c2 and c3 are the second point cloud data. Since the number of the first point cloud data and the second point cloud data is greater than 1, the centroid coordinates of b1 and b2, point B, and the centroid coordinates of c1, c2 and c3 are calculated and obtained. The K point and the B point are connected as the first vector. , the coordinates of point K and point C are used as the second vector , calculate the first component of the first vector in the normal direction , and the second component of the second vector in the normal direction , the difference between the first and second components As the displacement change value of the first grid, the components of the first vector and the second vector in the normal direction are calculated to obtain the displacement of the first grid in the vertical direction of the slope at different time points. By focusing only on the normal component, the error caused by terrain undulation, measurement error or data noise can be reduced, thereby improving the accuracy of displacement measurement.

[0035] Screening the optimal pattern from rainfall patterns includes the following steps: Calculate the average value of the actual rainfall of the first change curve at multiple first time points. If the first duration is greater than the second threshold, the average value is greater than the third threshold, and the first time point is a maximum point, the optimal mode is matched to the first mode. If the first duration is greater than the second threshold, the average value is less than the third threshold, and the first time point is a minimum point, the optimal mode is matched to the second mode. If the first duration is less than the second threshold, the average value is greater than the third threshold, and the first time point is a maximum point, the optimal mode is matched to the third mode. If the first duration is less than the second threshold, the average value is less than the third threshold, and the first time point is a minimum point, the optimal mode is matched to the fourth mode.

[0036] Specifically, if Figure 3 As shown, the actual rainfall at the first time points 0:30, 1:00, 1:30, 2:00, 2:30, 3:00, 3:30, and 4:00 is obtained at intervals of 30 minutes, which are 22, 24, 26, 28, 30, 34, 36, and 38 respectively. The average actual rainfall of the first duration 0-4h is 29.75, because the first duration of 4h is greater than the second threshold, assuming it is 2 hours, and the average value of 29.75 is greater than the third threshold, assuming it is 25, and the first time point is at a maximum at 4h. value point, it is determined that the optimal prediction mode match of the first time length 0-s1 is the first mode, that is, the long-term heavy rainfall mode. According to this step, the rainfall mode in each first time length is obtained. The rainfall mode of the first time length 4-8h is the second mode, that is, the long-term weak rainfall mode. The rainfall mode of the first time length 8-12h is the first mode, that is, the long-term heavy rainfall mode. The rainfall mode of the first time length 12-20h is the second mode, that is, the long-term weak rainfall mode. The rainfall mode of the first time length 20-22h is the third mode, that is, the short-term heavy rainfall mode.

[0037] Obtaining the optimal time point within each first duration includes the following steps: The average rainfall at the first time point of multiple historical rainfall data corresponding to the optimal mode is obtained, the difference between the actual rainfall at each first time point and the corresponding average rainfall is calculated, and the time point when the difference is less than or equal to the fourth threshold is defined as the optimal time point.

[0038] Specifically, the optimal mode corresponding to the first duration of 0-4h is obtained: in the long-term heavy rainfall mode, the first time points in this rainfall mode include 0:30, 1:00, 1:30, 2:00, 2:30, 3:00, 3:30, and 4:00, and the average rainfall at these first time points is obtained, which is 25, 26, 27, 28, 29, 30, 31, and 33 respectively. The error between the average rainfall and the actual rainfall at each first time point is calculated, and the difference at the first time point 2:00 is 0 by calculation, which is less than or equal to the fourth threshold value 0. The moment of 2:00 is defined as the optimal time point. By confirming the optimal time point, a rainfall with the highest accuracy and the smallest error can be obtained as a representative rainfall. Taking the representative rainfall as a benchmark, the prediction accuracy of future rainfall events can be improved.

[0039] Building a predictive model involves the following steps: The time period of each historical rainfall data is defined as the historical duration, and the historical duration is divided into multiple first durations. The corresponding optimal mode is matched for each first duration, and the representative rainfall under the optimal mode is calculated. A division point is set within the historical duration, and the optimal mode and representative rainfall of each first duration in the time period before the division point are used as input features. After the division point, a first duration is selected, and the corresponding optimal mode and representative rainfall are used as output features. This step is repeated to process each historical rainfall data, and a neural network model is constructed. The neural network model is trained based on the input features and the corresponding output features to obtain a prediction model.

[0040] Specifically, based on the method of dividing the first time period within the target time period, and matching the optimal pattern for each first time period and calculating the representative rainfall, each historical time period is divided into multiple first time periods (the historical time period also includes the target time period, which can be rainfall data for multiple days). A point is selected in the time series of the historical time period, and the historical rainfall data is divided into two parts: before the division point and after the division point. This division point is used to distinguish between model training data and model verification data. The optimal patterns and representative rainfall of all first time periods before the division point are used as input features (i.e., feature set), and the rainfall pattern and representative rainfall of the first "first time period" after the division point are used as output features (i.e., target value). The purpose of this is to train the model so that it can predict the optimal pattern and representative rainfall in the future time period based on the historical rainfall data. A neural network model is constructed using the above input features and output features. Through multiple iterations, the model parameters are adjusted. After training, the neural network model can capture the time series characteristics and patterns of the rainfall data, thereby predicting the optimal pattern and representative rainfall in the future time period. Based on the representative precipitation in the future time period, the pore water pressure distribution of the soil in the future time period is obtained in advance, and the stability of the slope is further analyzed.

[0041] Generating a first safety factor for a first region at a future time period includes the following steps: If the predicted rainfall in the future time period is greater than the fifth threshold, the first pore water pressure distribution is used as the initial value and the predicted rainfall is used as the boundary condition. The second pore water pressure distribution of the soil in the future time period is predicted based on the one-dimensional permeability model. The second pore water pressure of the soil in the first area is obtained based on the second pore water pressure distribution. The effective stress of the soil in the first area is calculated based on the second pore water pressure. The effective stress is the difference between the total stress of the soil in the first area and the second pore water pressure. The shear strength of the soil under the second pore water pressure is obtained, and the ratio of the shear strength to the effective stress is used as the first safety factor of the first area in the future time period.

[0042] Specifically, the predicted rainfall for a future time period is obtained and compared with a preset fifth threshold. If the predicted rainfall exceeds this threshold, it indicates that a large amount of water may infiltrate the soil, affecting the stability of the slope. Using a one-dimensional infiltration model, the current pore water pressure distribution is used as the initial value and the predicted rainfall as the boundary condition to predict the soil pore water pressure distribution for the future time period. Based on the predicted pore water pressure distribution, the effective stress of the soil in the first region is calculated. Effective stress refers to the difference between the total soil stress and the pore water pressure and is an important parameter for assessing soil stability. Given the second pore water pressure, the shear strength of the soil is obtained. Shear strength is the ability of soil to resist shear deformation and is crucial for assessing slope stability. The ratio of the soil shear strength to the effective stress is used as the first safety factor for the first region at the predicted future time period. The safety factor is a direct indicator of slope stability; a higher factor indicates a more stable slope. By predicting the safety factor for future time periods, potential slope instability risks can be identified in advance, providing a time window for taking preventive measures.

[0043] like Figure 4 As shown, the present invention also provides a digital twin-based intelligent slope monitoring system, which is used to implement the above method. The system mainly includes: A creation module is used to obtain first point cloud data of the slope surface at a preset time point in the past and create a digital twin model based on the first point cloud data.

[0044] The calculation module is used to obtain the second point cloud data of the slope surface at the current time point, convert the first point cloud data, the second point cloud data and the digital twin model into the same coordinate system, mesh the digital twin model based on the preset size to generate multiple first grids, define the time between the preset time point and the current time point as the target time length, calculate the displacement change value at each first grid within the target time length, and display it in the digital twin model.

[0045] A division module is used to define the first grid whose displacement change value is greater than the first threshold as a first area, obtain the actual rainfall in the first area within the target time length, establish a first coordinate system with time as the horizontal axis and the actual rainfall as the vertical axis, plot a first change curve of the actual rainfall over time in the first coordinate system, obtain the extreme value points of the first change curve, define the time point corresponding to each extreme value point as a first time point, and divide the target time length into multiple first time lengths of unequal time based on the first time point.

[0046] The prediction module is used to set multiple rainfall patterns, each rainfall pattern corresponds to multiple historical rainfall data, and a rainfall pattern is screened out from the rainfall patterns as the optimal pattern based on the actual rainfall in the first time period. The optimal time point in each first time period is obtained based on the optimal pattern and the historical rainfall data. The actual rainfall corresponding to the optimal time point is used as the representative rainfall in the first time period. A prediction model is constructed based on the historical rainfall data. According to the optimal pattern and representative rainfall in each first time period, the optimal pattern and representative rainfall in the future time period are predicted and generated, and the representative rainfall in the future time period is defined as the predicted rainfall.

[0047] The monitoring module is used to obtain historical rainfall within the target time period, obtain the first pore water pressure distribution of the soil in the first area based on the historical rainfall, and predict and generate the first safety factor of the first area in the future time period based on the predicted rainfall and the first pore water pressure distribution. The average of all the first safety factors of the slope is used as the final safety factor of the slope and displayed in the digital twin model.

[0048] The present invention also provides a computer storage medium, which stores program instructions. When the program instructions are executed, the device where the computer storage medium is located is controlled to execute the above method.

[0049] It should be understood that, although the various steps in the flow chart of each embodiment of the present invention are shown in sequence according to the indication of the arrows, these steps are not necessarily performed in sequence according to the order indicated by the arrows. Unless otherwise specified herein, the execution of these steps is not strictly limited in order, and these steps can be performed in other orders. Moreover, at least a portion of the steps in each embodiment may include a plurality of sub-steps or a plurality of stages, and these sub-steps or stages are not necessarily performed at the same time, but can be performed at different times, and the execution order of these sub-steps or stages is not necessarily performed in sequence, but can be performed in turn or alternately with at least a portion of other steps or sub-steps or stages of other steps.

[0050] Those skilled in the art will appreciate that all or part of the processes in the above-described method embodiments can be implemented by instructing the relevant hardware through a computer program. The above-described program can be stored in a non-volatile computer-readable storage medium. When executed, the program can include the processes of the above-described method embodiments. Any reference to memory, storage, database, or other media used in the various embodiments provided herein may include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), RAMbus direct RAM (RDRAM), direct RAMbus dynamic RAM (DRDRAM), and RAMbus dynamic RAM (RDRAM).

[0051] The technical features of the above-mentioned embodiments can be combined arbitrarily. In order to make the description concise, not all possible combinations of the technical features in the above-mentioned embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0052] The above embodiments merely illustrate several implementations of the present invention, and while their descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that a person skilled in the art would be able to make numerous variations and improvements without departing from the spirit of the present invention, all of which fall within the scope of protection of the present invention. Therefore, the scope of protection of the present invention shall be determined by the appended claims.

[0053] The above are only preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A slope intelligent monitoring method based on digital twins, characterized in that: The method comprises the following steps: S1: Acquire first point cloud data of the slope surface at a preset time point in the past, and create a digital twin model based on the first point cloud data; S2: Acquire the second point cloud data of the slope surface at the current time point, convert the first point cloud data, the second point cloud data and the digital twin model into the same coordinate system, mesh the digital twin model based on a preset size to generate a plurality of first grids, define the duration from the preset time point to the current time point as a target duration, calculate the displacement change value at each of the first grids within the target duration, and display it in the digital twin model; S3: defining the first grid whose displacement change value is greater than the first threshold as a first area, obtaining the actual rainfall in the first area within the target duration, establishing a first coordinate system with time as the horizontal axis and the actual rainfall as the vertical axis, plotting a first variation curve of the actual rainfall over time in the first coordinate system, obtaining the extreme value points of the first variation curve, defining the time point corresponding to each extreme value point as a first time point, and dividing the target duration into a plurality of first durations of unequal time based on the first time point; S4: setting a plurality of rainfall modes, each of which corresponds to a plurality of historical rainfall data, selecting a rainfall mode from the rainfall modes as an optimal mode based on the actual rainfall within the first time period, obtaining an optimal time point within each of the first time period based on the optimal mode and the historical rainfall data, using the actual rainfall corresponding to the optimal time point as the representative rainfall within the first time period, building a prediction model based on the historical rainfall data, predicting and generating the optimal mode and the representative rainfall within a future time period based on the optimal mode and the representative rainfall within each of the first time periods, and defining the representative rainfall within the future time period as predicted rainfall; S5: Obtain historical rainfall within the target time period, obtain a first pore water pressure distribution of the soil in the first area based on the historical rainfall, predict and generate a first safety factor of the first area in the future time period based on the predicted rainfall and the first pore water pressure distribution, take the average of all the first safety factors of the slope as the final safety factor of the slope, and display it in the digital twin model.

2. The method according to claim 1, characterized in that Calculating the displacement change value at each of the first grids comprises the following steps: Define a vertex of the first mesh in the digital twin model as a first node, and define the three-dimensional position coordinates after the first point cloud data and the second point cloud data are converted as first position coordinates and second position coordinates, respectively; Obtain a first number of the first point cloud data and a second number of the second point cloud data in the first grid. If both the first number and the second number are less than or equal to 1, a vector formed by taking the first node as a starting point and the first position coordinate of the first point cloud data as an end point is defined as a first vector, and a vector formed by taking the first node as a starting point and the second position coordinate of the second point cloud data as an end point is defined as a second vector; If the first number is greater than 1, the barycentric coordinates of all the first point cloud data in the first grid are calculated, and a vector formed by taking the first node as a starting point and the barycentric coordinates of the first point cloud data as an end point is used as the first vector; if the second number is greater than 1, the barycentric coordinates of all the second point cloud data in the first grid are calculated, and a vector formed by taking the first node as a starting point and the barycentric coordinates of the second point cloud data as an end point is used as the second vector; Get the unit vector of the first grid in the normal direction , respectively calculate the first component of the first vector and the second vector in the normal direction based on the first formula and the second formula and the second component , the first formula is: , the second formula is: ,in, and are respectively the first vector and the second vector, and the difference between the first component and the second component is defined as the displacement change value of the first grid.

3. The method according to claim 1, characterized in that Screening out the optimal mode from the rainfall modes comprises the following steps: Calculate the average value of the actual rainfall of the first change curve at multiple first time points. If the first duration is greater than the second threshold, the average value is greater than the third threshold, and the first time point is a maximum point, match the optimal mode to the first mode. If the first duration is greater than the second threshold, the average value is less than the third threshold, and the first time point is a minimum point, match the optimal mode to the second mode. If the first duration is less than the second threshold, the average value is greater than the third threshold, and the first time point is a maximum point, match the optimal mode to the third mode. If the first duration is less than the second threshold, the average value is less than the third threshold, and the first time point is a minimum point, match the optimal mode to the fourth mode.

4. The method according to claim 3, characterized in that Obtaining the optimal time point within each of the first durations comprises the following steps: Obtain the average rainfall of multiple pieces of historical rainfall data corresponding to the optimal mode at the first time point, calculate the difference between the actual rainfall at each of the first time points and the corresponding average rainfall, and define the time point when the difference is less than or equal to a fourth threshold as the optimal time point.

5. The method according to claim 3, characterized in that: Building a predictive model involves the following steps: The time period in which each piece of the historical rainfall data is located is defined as a historical time length, the historical time length is divided into multiple first time lengths, the corresponding optimal mode is matched for each first time length, and the representative rainfall under the optimal mode is calculated, a division point is set in the historical time length, the optimal mode and the representative rainfall of each first time length in the time period before the division point are used as input features, one of the first time lengths is selected after the division point, and the corresponding optimal mode and the representative rainfall are used as output features, this step is repeated to process each piece of the historical rainfall data, a neural network model is constructed, and the neural network model is trained based on the input features and the corresponding output features to obtain the prediction model.

6. The method according to claim 1, characterized in that Generating a first safety factor for the first area in the future time period comprises the following steps: If the predicted rainfall in the future time period is greater than the fifth threshold value, the first pore water pressure distribution is used as the initial value and the predicted rainfall is used as the boundary condition, and the second pore water pressure distribution of the soil in the future time period is predicted based on the one-dimensional permeability model, and the second pore water pressure of the soil in the first area is obtained based on the second pore water pressure distribution, and the effective stress of the soil in the first area is calculated based on the second pore water pressure, and the effective stress is the difference between the total stress of the soil in the first area and the second pore water pressure, and the shear strength of the soil under the second pore water pressure is obtained, and the ratio of the shear strength to the effective stress is used as the first safety factor of the first area in the future time period.

7. A digital twin-based intelligent slope monitoring system, used to implement the method according to any one of claims 1 to 6, characterized in that: The system includes the following modules: A creation module, used to obtain first point cloud data of the slope surface at a preset time point in the past, and create a digital twin model based on the first point cloud data; a calculation module, used to obtain second point cloud data of the slope surface at the current time point, convert the first point cloud data, the second point cloud data and the digital twin model into the same coordinate system, mesh the digital twin model based on a preset size to generate a plurality of first grids, define the duration from the preset time point to the current time point as a target duration, calculate the displacement change value at each of the first grids within the target duration, and display it in the digital twin model; A division module, used for defining the first grid whose displacement change value is greater than a first threshold as a first area, obtaining the actual rainfall in the first area within the target duration, establishing a first coordinate system with time as the horizontal axis and the actual rainfall as the vertical axis, plotting a first variation curve of the actual rainfall over time in the first coordinate system, obtaining extreme value points of the first variation curve, defining a time point corresponding to each extreme value point as a first time point, and dividing the target duration into a plurality of first durations of unequal time based on the first time point; A prediction module is used to set a plurality of rainfall modes, each of which corresponds to a plurality of historical rainfall data, select a rainfall mode from the rainfall modes as an optimal mode based on the actual rainfall in the first time period, obtain an optimal time point in each of the first time period based on the optimal mode and the historical rainfall data, use the actual rainfall corresponding to the optimal time point as a representative rainfall in the first time period, build a prediction model based on the historical rainfall data, predict and generate the optimal mode and the representative rainfall in a future time period according to the optimal mode and the representative rainfall in each of the first time periods, and define the representative rainfall in the future time period as predicted rainfall; A monitoring module is used to obtain historical rainfall within the target time period, obtain a first pore water pressure distribution of the soil in the first area based on the historical rainfall, predict and generate a first safety factor of the first area in the future time period based on the predicted rainfall and the first pore water pressure distribution, take the average of all the first safety factors of the slope as the final safety factor of the slope, and display it in the digital twin model.

8. A computer storage medium, characterized in that: The computer storage medium stores program instructions, wherein when the program instructions are executed, the device where the computer storage medium is located is controlled to execute the method according to any one of claims 1 to 6.

Images