Foundation stability analysis early warning method and system based on underground water level change

Abstract

The invention discloses a foundation stability analysis early warning method and system based on underground water level changes, and relates to the technical field of data processing. The method comprises the following steps: acquiring underground water level data of a target area, and constructing an underground water level space-time contour surface based on a water level value, sampling time and a sampling position in the underground water level data; a plurality of monitoring points on the underground water level space-time contour surface are determined, and the underground water level gradient change rate of each monitoring point is calculated; determining a foundation stress concentration area in the underground water level space-time contour surface according to the change rate of each underground water level gradient; the foundation stress concentration area is divided into a plurality of grid areas, and the deformation rate of each grid area is received; and for each grid region, when the deformation rate is greater than a rate threshold, determining that the grid region is a foundation instability risk region, and generating early warning information of the foundation instability risk region. By implementing the technical scheme provided by the invention, the effect of improving the accuracy of the foundation stability analysis early warning based on the underground water level change is achieved.

Description

Technical Field

[0001] The present application relates to the field of data processing technology, and specifically to a foundation stability analysis and early warning method, system, electronic equipment and storage medium based on groundwater level changes. Background Art

[0002] With the acceleration of urbanization, the construction of high-rise buildings and large-scale infrastructure is increasing, and the monitoring and management of foundation stability has become particularly important. The stability of the foundation is directly related to the safety and use function of the building, and the change of the groundwater level is a key factor affecting the stability of the foundation. The increase or decrease of the groundwater level may affect the bearing capacity of the soil, thereby affecting the stability of the structure above.

[0003] At present, the existing foundation stability analysis and early warning methods based on groundwater level changes usually collect groundwater level change data at a fixed location, and use this data to analyze and warn the foundation stability. However, in practical applications, only through fixed location data analysis, there is a lack of comprehensive consideration of groundwater level changes in a large area. This localized monitoring method is difficult to capture the overall impact of groundwater level changes on the stability of the entire foundation, thereby reducing the accuracy of foundation stability analysis and early warning based on groundwater level changes. Summary of the invention

[0004] The present application provides a foundation stability analysis and early warning method and system based on groundwater level changes, which has the effect of improving the accuracy of foundation stability analysis and early warning based on groundwater level changes.

[0005] In a first aspect, the present application provides a foundation stability analysis and early warning method based on groundwater level changes, comprising: Acquire groundwater level data of a target area, and construct a spatiotemporal isosurface of groundwater level in the target area based on the water level value, sampling time and sampling location in the groundwater level data; Determine a plurality of monitoring points on the groundwater level time-space isosurface, and calculate the groundwater level gradient change rate of each of the monitoring points in the groundwater level time-space isosurface; Determining the foundation stress concentration area in the groundwater level time-space isosurface according to each of the groundwater level gradient change rates; Dividing the foundation stress concentration area into a plurality of grid areas, and receiving the deformation rate of each of the grid areas; For each of the grid areas, when the deformation rate is greater than a rate threshold, the grid area is determined to be a foundation instability risk area, and early warning information of the foundation instability risk area is generated.

[0006] In a second aspect of the present application, a foundation stability analysis and early warning system based on groundwater level changes is provided, the system comprising: A data acquisition module, used to acquire groundwater level data of a target area, and construct a spatiotemporal isosurface of groundwater level in the target area based on the water level value, sampling time and sampling location in the groundwater level data; A gradient change rate determination module is used to determine a plurality of monitoring points on the groundwater level time-space isosurface, and calculate the groundwater level gradient change rate of each of the monitoring points in the groundwater level time-space isosurface; An area determination module, used to determine the foundation stress concentration area in the groundwater level time-space isosurface according to each groundwater level gradient change rate; The foundation analysis and early warning module is used to divide the foundation stress concentration area into multiple grid areas and receive the deformation rate of each grid area; for each grid area, when the deformation rate is greater than a rate threshold, the grid area is determined to be a foundation instability risk area, and early warning information of the foundation instability risk area is generated.

[0007] In a third aspect of the present application, an electronic device is provided, comprising a memory, a processor, and a program stored in the memory and executable on the processor, wherein the program can implement a foundation stability analysis and early warning method based on groundwater level changes when loaded and executed by the processor.

[0008] In a fourth aspect of the present application, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the processor implements a foundation stability analysis and early warning method based on groundwater level changes.

[0009] In summary, one or more technical solutions provided in the embodiments of the present application have at least the following technical effects or advantages: By adopting the above technical scheme, the spatiotemporal contour surface of the groundwater level in the target area is constructed in combination with the acquired groundwater level data in the target area, which realizes the overall expression of the groundwater level change characteristics in a large range, breaking through the spatial limitations of the traditional fixed point monitoring method. On the basis of the spatiotemporal contour surface, by calculating the groundwater level gradient change rate of multiple monitoring points, the dynamic change characteristics of groundwater flow can be fully captured, providing more complete data support for foundation stability analysis. The identified foundation stress concentration area is divided into grid areas and the deformation rate is monitored. The correlation analysis mechanism between the overall change of groundwater level and the local foundation deformation is established, and the system monitoring from macro water level change to micro deformation response is realized. When the deformation rate of the grid area exceeds the threshold, the system automatically marks the foundation instability risk area and generates warning information, forming a precise warning system based on comprehensive analysis. By conducting multi-point analysis and warning of large-scale groundwater level changes, the impact of groundwater level changes on foundation stability can be more accurately evaluated, and the accuracy of foundation stability analysis and warning based on groundwater level changes can be improved. BRIEF DESCRIPTION OF THE DRAWINGS

[0010] Figure 1 It is a flow chart of a foundation stability analysis and early warning method based on groundwater level changes provided in an embodiment of the present application; Figure 2 It is a structural schematic diagram of a foundation stability analysis and early warning system based on groundwater level changes provided in an embodiment of the present application; Figure 3 It is a structural schematic diagram of an electronic device provided in an embodiment of the present application.

[0011] Description of reference numerals: 300, electronic device; 301, processor; 302, communication bus; 303, user interface; 304, network interface; 305, memory. DETAILED DESCRIPTION

[0012] In order to enable technicians in this field to better understand the technical solutions in this specification, the technical solutions in the embodiments of this specification will be clearly and completely described below in conjunction with the drawings in the embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this application, not all of the embodiments.

[0013] In the description of the embodiments of the present application, words such as "for example" or "for example" are used to indicate examples, illustrations or explanations. Any embodiment or design described as "for example" or "for example" in the embodiments of the present application should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of words such as "for example" or "for example" is intended to present related concepts in a specific way.

[0014] In the description of the embodiments of the present application, the meaning of the term "multiple" refers to two or more. For example, multiple systems refer to two or more systems, and multiple screen terminals refer to two or more screen terminals. In addition, the terms "first" and "second" are used for descriptive purposes only and cannot be understood as indicating or implying relative importance or implicitly indicating the indicated technical features. Thus, the features defined as "first" and "second" may explicitly or implicitly include one or more of the features. The terms "include", "comprise", "have" and their variations all mean "including but not limited to", unless otherwise specifically emphasized.

[0015] The present application embodiment provides a foundation stability analysis and early warning method based on groundwater level changes. In one embodiment, please refer to Figure 1 , Figure 1 This is a flow chart of a foundation stability analysis and early warning method based on groundwater level changes provided by an embodiment of the present application. The method can be implemented by a computer program, which can be integrated into an application or run as an independent tool application. The method can also be implemented by a single-chip microcomputer or run on a foundation stability analysis and early warning system based on groundwater level changes based on a von Neumann system. Specifically, the method can include the following steps: Step 101: Obtain groundwater level data of a target area, and construct a spatiotemporal isosurface of groundwater level of the target area based on the water level value, sampling time and sampling location in the groundwater level data.

[0016] Among them, groundwater level data refers to a comprehensive data set collected by water level monitoring devices in the target area, reflecting the spatial distribution and temporal changes of groundwater levels. The data includes three basic elements: water level value, sampling time, and sampling location. Among them, the water level value represents the elevation data of the groundwater level relative to a certain reference surface; the sampling time records the specific moment of data collection; and the sampling location determines the specific spatial location of the monitoring point through the longitude and latitude coordinates.

[0017] The spatiotemporal isosurface of groundwater level refers to a set of surfaces formed by continuously processing the groundwater level data at different times and locations in the target area through mathematical interpolation methods in a three-dimensional spatial coordinate system. This surface reflects the continuous change characteristics of the groundwater level in both space and time dimensions.

[0018] Specifically, firstly, the groundwater level data of the target area is obtained through multiple water level monitoring devices. Each water level monitoring device is equipped with a water level sensor, a time recording unit and a position positioning module to collect water level value, sampling time and sampling location information in real time. In order to construct a complete spatiotemporal isosurface of groundwater level, the target area needs to be divided into multiple spatial units with preset spatial specifications. The size of the spatial unit can be set according to the actual monitoring needs. For example, the target area can be divided into a square grid of 10 meters × 10 meters. For each spatial unit, based on the acquired groundwater level data, the target water level value, target sampling time and target sampling location corresponding to the spatial unit are determined by the spatiotemporal interpolation algorithm. Specifically, the sampling time, target water level value and target sampling location of each spatial unit can be interpolated by the Kriging interpolation method, so as to obtain the groundwater level distribution characteristics at multiple time points. Subsequently, the groundwater level distribution data at each time point are connected in chronological order to generate the spatiotemporal isosurface of the groundwater level in the target area. The spatiotemporal contour surface of groundwater level not only contains the spatial distribution information of groundwater level, but also reflects the change characteristics of groundwater level over time, providing comprehensive data support for subsequent foundation stability analysis. By constructing the spatiotemporal contour surface of groundwater level, the limitations of traditional fixed monitoring points are broken, and the overall grasp of groundwater level changes in the target area is achieved. At the same time, the influencing factors of time and space are considered, laying the foundation for accurate assessment of foundation stability. Practice has shown that analysis based on the spatiotemporal contour surface of groundwater level can effectively identify areas with abnormal groundwater level changes and improve the accuracy of foundation stability warning.

[0019] Based on the above embodiment, as an optional embodiment, in step 101: constructing a spatiotemporal isosurface of groundwater level in the target area based on the water level value, sampling time and sampling location in the groundwater level data, this step may also include the following steps: Step 201: Divide the target area into a plurality of spatial units of preset spatial specifications; determine the target water level value, target sampling time and target sampling position corresponding to each spatial unit based on the water level value, sampling time and sampling position in the groundwater level data.

[0020] Among them, the spatial unit of the preset spatial specification refers to a regular grid unit formed by dividing the target area according to a preset uniform size standard, and each spatial unit has the same geometric shape and size. In the embodiment of the present application, the spatial unit of the preset spatial specification can be understood as a standardized square area, and its side length can be set according to actual monitoring needs, such as 10 meters × 10 meters or 20 meters × 20 meters, forming a grid monitoring system covering the entire target area.

[0021] Specifically, in order to achieve comprehensive monitoring of groundwater level changes in the entire region, the target area needs to be divided into regular grid-shaped spatial units. According to the monitoring accuracy requirements and actual engineering needs, the target area is divided into several square spatial units of equal size, and the specifications of each spatial unit can be set to 10 meters × 10 meters. The spatial units are divided in a grid-like manner to ensure that there is no overlap or omission between the spatial units, thereby forming a complete monitoring network. For each spatial unit, it is necessary to determine its corresponding target water level value, target sampling time and target sampling location. First, the existing groundwater level monitoring data in the spatial unit is obtained. If there is an actual monitoring point in the spatial unit, the water level value, sampling time and sampling location of the monitoring point are directly used as the target value of the spatial unit; if there is no actual monitoring point in the spatial unit, the target water level value of the spatial unit is calculated based on the monitoring data of the adjacent spatial units, using a spatial interpolation method, such as the Kriging interpolation method, and the interpolation calculation time is used as the target sampling time, and the center point coordinates of the spatial unit are used as the target sampling location. In this way, discrete monitoring data can be transformed into a continuous data field covering the entire target area, providing a data basis for constructing spatiotemporal isosurfaces of groundwater levels.

[0022] Step 202: interpolate the sampling time, target water level value and target sampling position of each spatial unit to obtain groundwater level distribution characteristics at multiple time points.

[0023] Specifically, since the target water level value acquisition time of each spatial unit may be different, and the data of some spatial units may be obtained through spatial interpolation, in order to accurately describe the spatiotemporal variation characteristics of the groundwater level, it is necessary to perform spatiotemporal interpolation processing on the data of each spatial unit. In this embodiment, the interpolation time series is first determined, and multiple time points can be selected according to fixed time intervals, such as every hour or every day. For each selected time point, the target water level value of each spatial unit is interpolated and calculated using the Kriging interpolation method. Specifically, a variogram model is established based on the spatial coordinates of the target sampling position and the corresponding target water level value, and the water level estimation value of each spatial unit at this time point is calculated considering the correlation of the spatial position. At the same time, for the case where the target sampling moment is inconsistent with the selected time point, the water level value is corrected by the interpolation method of the time dimension to ensure that the interpolation result can accurately reflect the change trend of the groundwater level over time. Through this spatiotemporal interpolation processing, the groundwater level distribution characteristics of the target area at each selected time point can be obtained, and these distribution characteristics not only contain the water level change law of the spatial dimension, but also reflect the dynamic characteristics of the time dimension. This interpolation processing method effectively solves the problem of spatiotemporal discontinuity of monitoring data, provides a reliable data basis for constructing a complete spatiotemporal isosurface of groundwater level, and thus improves the accuracy of foundation stability analysis. Practice shows that the water level distribution characteristics obtained based on spatiotemporal interpolation processing can better reflect the changing trend of groundwater level and help to timely discover abnormal change areas.

[0024] Step 203: Connect the groundwater level distribution data at each time point in chronological order to generate a spatiotemporal isosurface of the groundwater level in the target area.

[0025] Specifically, the groundwater level distribution characteristics of multiple time points obtained by spatiotemporal interpolation are first arranged in chronological order to form a time series data set. For the water level changes between adjacent time points, the time interpolation algorithm is used for transition processing to ensure that the changes in the groundwater level present continuous and smooth characteristics. The cubic spline interpolation method is used to interpolate the time dimension, so that the changes in the groundwater level on the time axis are more in line with the actual situation. At the same time, considering the continuity requirements of the spatial dimension, the water level differences between adjacent grid units in space are also smoothed to avoid unreasonable mutations. Through this spatiotemporal connection processing, a four-dimensional spatiotemporal isosurface of the groundwater level is finally generated. The isosurface not only reflects the spatial distribution state of the groundwater level at any time, but also reflects the trend of the water level at any position over time. This continuous spatiotemporal isosurface expression provides a complete data model for the subsequent analysis of the impact of groundwater level changes on foundation stability, and can more accurately identify and predict water level changes that may lead to foundation instability. Practice has proved that the analysis method based on spatiotemporal isosurfaces has significantly improved the accuracy of foundation stability warning.

[0026] Step 102: Determine a plurality of monitoring points on the groundwater level spatiotemporal contour surface, and calculate the groundwater level gradient change rate of each monitoring point in the groundwater level spatiotemporal contour surface.

[0027] Among them, the groundwater level gradient change rate refers to the ratio of the change in water level gradient per unit time at the monitoring point location on the groundwater level time-space isosurface to the initial water level gradient. This indicator reflects the comprehensive change rate of groundwater level in time and space dimensions.

[0028] Specifically, firstly, on the spatiotemporal isosurface of groundwater level, areas with large water level change rates, areas with large terrain slopes, and areas with concentrated building loads are identified. These areas are usually locations where foundation stability is more sensitive. Then, dense monitoring points are arranged in these key areas, and a sparser distribution of monitoring points can be used in other areas. For each determined monitoring point, its groundwater level gradient change rate needs to be calculated. When calculating, first determine the water level value of the monitoring point at different times, and calculate the water level change between adjacent time points; secondly, determine the spatial water level distribution around the monitoring point, and calculate the water level gradient in the spatial direction; finally, comprehensively consider the change characteristics of the two dimensions of time and space, and obtain the groundwater level gradient change rate of the monitoring point. The groundwater level gradient change rate can fully reflect the dynamic characteristics of groundwater level changes, especially the rapid change process that may lead to foundation instability.

[0029] Based on the above embodiment, as an optional embodiment, in step 102: determining multiple monitoring points on the spatiotemporal isosurface of groundwater level, this step may also include the following steps: Step 301: Select multiple time sections on the groundwater level time-space isosurface; use the terrain slope mutation area, building distribution area and groundwater level mutation area in each time section as the point area on the groundwater level time-space isosurface.

[0030] Among them, the time section refers to the three-dimensional spatial water level distribution data corresponding to a specific moment selected on the groundwater level space-time isosurface, which reflects the spatial distribution status of the groundwater level in the target area at that moment.

[0031] The terrain slope mutation area refers to the area within the target area where the surface elevation changes significantly in spatial distribution. Specifically, it is manifested as an area where the slope difference between adjacent spatial locations exceeds a preset threshold, such as 15 degrees, or the slope change rate is greater than a preset rate.

[0032] The building distribution area refers to the spatial range within the target area where building loads exist, including the floor area of ​​the building itself and its impact range.

[0033] The groundwater level mutation area refers to the area where the water level shows significant drop or abnormal change gradient in spatial distribution on the groundwater level spatiotemporal isosurface, specifically, the water level difference between adjacent spatial units exceeds the preset threshold, such as 1 meter, or the water level gradient change rate exceeds the preset standard.

[0034] The point area refers to the spatial range that requires key monitoring, which is determined on the groundwater level spatiotemporal isosurface through comprehensive analysis of terrain slope mutation areas, building distribution areas, and groundwater level mutation areas.

[0035] Specifically, since the stability of the foundation is affected by many factors, it is necessary to scientifically select monitoring points on the groundwater level spatiotemporal contour surface. In this embodiment, first, multiple time sections are selected on the groundwater level spatiotemporal contour surface at fixed time intervals, such as every day or every week. These time sections can reflect the groundwater level distribution characteristics in different periods. For each time section, three types of characteristic areas are focused on: terrain slope mutation area, building distribution area and groundwater level mutation area. The terrain slope mutation area refers to the area where the terrain slope changes more than the preset threshold, such as the area where the slope changes more than 15 degrees. These areas are more prone to foundation instability due to the drastic terrain changes; the building distribution area refers to the area where there is a building load, and the gravity load of the building will have a significant impact on the stability of the foundation; the groundwater level mutation area refers to the area where the groundwater level shows a significant mutation in spatial distribution, such as the water level drop exceeds 1 meter. This mutation may cause groundwater flow to intensify and affect the stability of the foundation. These three types of characteristic areas are identified by numerical analysis methods and used as point areas to ensure that the layout of monitoring points can cover these potential high-risk areas. This method comprehensively considers multiple influencing factors such as terrain, building load and groundwater level changes, making the layout of monitoring points more targeted.

[0036] Step 302: Determine a plurality of monitoring points on the spatiotemporal isosurface of the groundwater level according to the locations of the point areas.

[0037] Specifically, the layout of monitoring points needs to comprehensively consider the spatial distribution characteristics and risk level of the point area. First, the risk level of each point area is evaluated. The evaluation indicators include the degree of change of terrain slope, the size of building load, the amplitude of groundwater level mutation, etc. These indicators are quantified and scored and weighted to obtain the comprehensive risk level of each point area. Secondly, the density of monitoring points is determined according to the risk level of the point area. The higher the risk level, the denser the monitoring points are. Specifically, for high-risk point areas with multiple risk factors (such as both terrain slope mutation and groundwater level mutation), the distance between monitoring points can be set to 20-30 meters; for medium-risk point areas with only a single risk factor, the distance between monitoring points can be set to 40-50 meters; for low-risk point areas, the distance between monitoring points can be appropriately increased. When determining the location of specific monitoring points, priority is given to locations where risk factors overlap, such as the junction of the building load area and the groundwater level mutation area. At the same time, considering the operability of monitoring, the layout of monitoring points also needs to avoid underground pipelines, green belts and other areas that are not conducive to the installation of monitoring equipment. Through this risk-level-based monitoring point layout method, a monitoring network with prominent focus and reasonable coverage has been formed. This monitoring point layout scheme can effectively capture early signs of foundation deformation and improve monitoring efficiency and early warning accuracy.

[0038] Based on the above embodiment, as an optional embodiment, in step 102: calculating the groundwater level gradient change rate of each monitoring point in the groundwater level spatiotemporal isosurface, this step may also include the following steps: Step 303: Obtain the water level value of each monitoring point in the spatiotemporal contour surface of the groundwater level at multiple sampling times; for each monitoring point, calculate the relative water level difference between the monitoring point and each monitoring point within a preset radius based on the water level value of the monitoring point at multiple sampling times.

[0039] Specifically, in order to accurately evaluate the impact of groundwater level changes on foundation stability, it is necessary to obtain and analyze the water level data of each monitoring point and its spatial correlation characteristics. First, the water level values ​​of each monitoring point at different sampling times are extracted from the spatiotemporal contour surface of the groundwater level. The sampling time interval can be set according to the actual monitoring needs, such as sampling once a day or once a week. After the data is acquired, for each monitoring point, it is necessary to analyze the water level relationship between it and the surrounding monitoring points. A preset radius (such as 50 meters) is set as the influence range, and the set of adjacent monitoring points is determined within this range. Then, for each sampling moment, the relative water level difference between the target monitoring point and each monitoring point within its influence range is calculated. The relative water level difference is calculated by subtracting the water level values ​​of other monitoring points within the influence range from the water level value of the target monitoring point to obtain a set of data reflecting the spatial variation characteristics of the local water level. This analysis method based on relative water level difference can effectively reflect the spatial gradient characteristics of the groundwater level in the local area. When the relative water level difference is large, it indicates that there is a significant hydraulic gradient in the area, which may lead to strong groundwater flow; when the relative water level difference changes significantly over time, it indicates that the groundwater flow field in the area is unstable. By analyzing the relative water level difference between monitoring points and its time evolution characteristics, abnormal changes in the local groundwater flow field can be discovered in a timely manner, providing a basis for assessing foundation stability risks.

[0040] Step 304: Determine the water level gradient of the monitoring point according to each relative water level difference and the relative distance between the monitoring point and each monitoring point within a preset radius.

[0041] The water level gradient refers to the rate of change of the groundwater level in space, indicating the amount of change in the water level per unit distance, and its value is equal to the ratio of the relative water level difference between adjacent monitoring points to the relative distance. In the embodiment of the present application, the water level gradient can be understood as a vector indicator that comprehensively reflects the characteristics of the local groundwater flow field.

[0042] Specifically, for each monitoring point, it is first necessary to determine the relative distance between it and each monitoring point within the preset radius. The relative distance can be calculated by the spatial coordinates of the monitoring point. Then, each relative water level difference is divided by the corresponding relative distance to obtain the water level gradient value in different directions. Considering the complexity of groundwater flow, the water level gradient in a single direction may not fully reflect the local water flow characteristics, so it is necessary to comprehensively consider the information of all monitoring points within the preset radius. The specific calculation method adopts the weighted average method: for each monitoring point within the preset radius, a weight coefficient is set according to its relative distance from the target monitoring point. The closer the distance, the greater the weight. The water level gradient in each direction is weighted averaged to obtain the comprehensive water level gradient of the target monitoring point. This water level gradient calculation method based on spatial weighting takes into account both the dominant influence of close-range monitoring points and the secondary influence of distant monitoring points, and can more accurately reflect the spatial variation characteristics of local groundwater levels. Practice has shown that the water level gradient calculated by this method is well representative, can reflect the intensity and direction characteristics of the groundwater flow field at the monitoring point location, and provides reliable basic data for subsequent evaluation of foundation stability.

[0043] Step 305: The changing rate of the water level gradient at each sampling time is taken as the groundwater level gradient changing rate of the monitoring point.

[0044] Specifically, in order to fully grasp the dynamic characteristics of the groundwater flow field, it is necessary to further analyze the time evolution law of the water level gradient. By calculating the rate of change of the water level gradient between adjacent sampling times, the groundwater level gradient change rate reflecting the dynamic characteristics of the groundwater flow field is obtained. The specific calculation method is: for each monitoring point, first obtain its water level gradient value at each sampling time, then calculate the difference between the water level gradients of two adjacent sampling times, and divide the difference by the corresponding time interval to obtain the groundwater level gradient change rate within the time period. Considering the continuity of the groundwater flow field change, the moving average method can be used to smooth the calculation results to eliminate the influence of accidental factors. The calculation results of the groundwater level gradient change rate can reflect the changing trend of the local groundwater flow field intensity. A large change rate indicates that the groundwater flow field in the area is in a state of rapid change, which may have an adverse effect on the stability of the foundation.

[0045] Step 103: Determine the foundation stress concentration area in the groundwater level time-space isosurface according to the groundwater level gradient change rate.

[0046] The foundation stress concentration area refers to the area in the spatiotemporal isosurface of the groundwater level where the foundation stress is significantly redistributed due to the drastic changes in the groundwater flow field. It is specifically manifested as a range where the groundwater level gradient change rate exceeds the preset threshold and has spatial continuity.

[0047] Specifically, in order to accurately identify the foundation stability risk area, it is necessary to determine the foundation stress concentration area based on the groundwater level gradient change rate. In this embodiment, a threshold value of the groundwater level gradient change rate (such as 0.01 / day) is first set. When the groundwater level gradient change rate of the monitoring point exceeds the threshold, it indicates that the groundwater flow field in the area where the point is located is undergoing significant changes. Then, the groundwater level gradient change rate of all monitoring points on the groundwater level spatiotemporal isosurface is spatially analyzed, and the areas between the monitoring points are interpolated using the Kriging interpolation method to obtain a continuous groundwater level gradient change rate distribution field. On this basis, the cluster analysis method is used to identify areas with large groundwater level gradient change rates and spatial continuity, and these areas are identified as potential foundation stress concentration areas. At the same time, considering the influence of geological conditions and building loads, it is necessary to superimpose the identified areas with the terrain slope mutation area and the building distribution area for analysis, focusing on the areas where multiple unfavorable factors are superimposed. This method for identifying foundation stress concentration areas based on the groundwater level gradient change rate can effectively capture the foundation stress redistribution phenomenon caused by changes in the groundwater flow field. Practice has shown that in these areas of foundation stress concentration, there are often changes in soil structure and a reduction in bearing capacity, which are high-risk areas for foundation instability. By timely identifying these areas, spatial references can be provided for foundation stability monitoring and early warning.

[0048] Based on the above embodiment, as an optional embodiment, in step 103: determining the foundation stress concentration area in the spatiotemporal isosurface of the groundwater level according to the groundwater level gradient change rate, this step may also include the following steps: Step 401: Calculate the time cumulative value of the gradient change rate of each monitoring point, and select the monitoring points whose time cumulative value is greater than the cumulative threshold as target monitoring points.

[0049] Specifically, in order to accurately evaluate the long-term impact of groundwater flow field changes on foundation stability, it is necessary to conduct a time accumulation analysis on the groundwater level gradient change rate at the monitoring point. First, it is necessary to determine the calculation period of time accumulation (such as 30 days). During this period, the groundwater level gradient change rate of each monitoring point is time-integrated to obtain the time accumulation value reflecting the intensity of groundwater flow field changes. The time accumulation value is calculated using the trapezoidal integration method. The groundwater level gradient change rate at each sampling time in the monitoring period is connected in chronological order, and the area between the formed curve and the time axis is calculated. A larger time accumulation value indicates that the monitoring point has experienced continuous or repeated groundwater flow field changes during the calculation period. This long-term accumulation effect is more likely to cause foundation stability problems. Therefore, it is necessary to set a reasonable accumulation threshold (such as 0.3). When the time accumulation value of a monitoring point exceeds this threshold, it is marked as a target monitoring point. This monitoring point screening method based on the time accumulation effect not only considers the intensity of groundwater flow field changes, but also the duration of changes, and can more comprehensively reflect the foundation stability risk.

[0050] Based on the above embodiment, as an optional embodiment, in step 401: calculating the time cumulative value of the gradient change rate of each monitoring point, this step may also include the following steps: Step 411: For each monitoring point, divide the gradient change rate into sub-gradient change rates corresponding to a plurality of time windows.

[0051] Specifically, first determine the size of the time window (such as 7 days) and the sliding step size (such as 1 day) to achieve a refined analysis of the gradient change rate. For each monitoring point, its complete gradient change rate time series is divided according to the set time window to obtain a series of continuous sub-gradient change rate sequences. The sliding window technology is used in the division process to ensure a certain overlap between adjacent time windows, which can ensure the continuity of the data and avoid the omission of important information. Each sub-gradient change rate sequence contains the groundwater flow field change characteristics of the monitoring point within the time window, and can reflect the change law on the local time scale. This time window-based division method enables us to observe the time evolution process of groundwater flow field changes more carefully and capture possible periodic changes or sudden anomalies.

[0052] Step 421: Perform differential operation on the sub-gradient change rates of adjacent time windows to obtain the change amount between each time window.

[0053] Specifically, in order to accurately quantify the dynamic change characteristics of the groundwater level gradient change rate, it is necessary to perform differential analysis on the sub-gradient change rates of adjacent time windows. The first-order difference method is used to calculate the difference in sub-gradient change rates between adjacent time windows, that is, the sub-gradient change rate of the next time window is subtracted from the sub-gradient change rate of the previous time window to obtain the change amount reflecting the change trend of the groundwater flow field. Considering that the actual monitoring data may have noise and fluctuations, before performing the differential operation, the sub-gradient change rate data in each time window is first statistically averaged to obtain a more representative window characteristic value. The change amount obtained by the differential operation can intuitively reflect the acceleration or deceleration trend of the groundwater flow field change. A positive value indicates that the change rate is increasing, a negative value indicates that the change rate is decreasing, and the absolute value of the change amount reflects the severity of the change. This analysis method based on time window difference can not only identify the trend of groundwater flow field change, but also quantitatively describe the rate of change, providing a basis for assessing foundation stability risks.

[0054] Step 431: sum the products of each change amount and the duration of the time window to obtain the time cumulative value of the gradient change rate of the monitoring point.

[0055] Specifically, in order to comprehensively evaluate the cumulative effect of the groundwater level gradient change rate, it is necessary to consider the combined influence of the change amount and time factors. First, the change amount corresponding to each time window is multiplied by the duration of the window to obtain the product value reflecting the cumulative effect within a single time window. The duration of the time window (such as 7 days) reflects the time scale of the change, while the change amount reflects the intensity of the change. The product of the two can comprehensively reflect the cumulative impact of the change in the groundwater flow field. The product values ​​of all time windows are summed to obtain the time cumulative value of the gradient change rate of the monitoring point, which comprehensively considers the intensity and duration of the groundwater flow field change during the entire monitoring period. In the summation process, a time weight factor can be introduced to make the recent changes contribute more to the cumulative value, thereby better reflecting the current risk status. This calculation method based on weighted accumulation can effectively identify the high-risk state caused by the superposition of persistent anomalies and periodic changes.

[0056] Step 402: Taking each target monitoring point as a starting point, expanding outward by a preset distance to form an expansion area of ​​each starting point.

[0057] Specifically, in order to comprehensively evaluate the foundation stability risk in the surrounding areas of the target monitoring points, it is necessary to conduct a spatial expansion analysis for each target monitoring point. Considering the spatial continuity and impact range of the groundwater flow field changes, each target monitoring point is taken as the starting point, and a preset distance (such as 50 meters) is evenly extended to form a circular expansion area. The determination of the expansion distance needs to consider geological conditions, soil properties and engineering experience to ensure that the main impact range of the groundwater flow field changes can be covered. When determining the expansion area, the buffer zone analysis method is used to draw a circular area with the target monitoring point as the center and the preset distance as the radius. When the expansion areas of multiple target monitoring points overlap, the overlapping parts need to be focused on as areas of concern, because these areas often reflect groundwater flow field anomalies in a larger range. This analysis method based on spatial expansion can effectively identify the impact range of groundwater flow field changes and provide a spatial basis for the subsequent division of foundation stress concentration areas.

[0058] Step 403: In each expansion area, multiple overlapping expansion areas are merged into one candidate stress area to obtain multiple candidate stress areas; and the candidate stress areas whose area is greater than the area threshold are used as foundation stress concentration areas in the spatiotemporal contour surface of the groundwater level.

[0059] Specifically, all the extended areas are first spatially analyzed to identify the extended areas with overlap. When two or more extended areas overlap in space, they are merged into a complete candidate stress area. The merging process uses a spatial connection algorithm to ensure that the boundaries of the merged area are smooth and continuous. The principle of this merging method is based on the spatial continuity characteristics of groundwater flow field changes. The overlap of adjacent areas often means a larger range of foundation stability risks. After the merger is completed, the area of ​​each candidate stress area needs to be calculated and compared with a pre-set area threshold (such as 5000 square meters). The purpose of setting the area threshold is to screen out foundation stress concentration areas with practical engineering significance, because areas that are too small may be caused by local anomalies or monitoring errors and do not have typical foundation stability risk characteristics. Candidate stress areas with an area greater than the threshold are determined as the final foundation stress concentration areas. These areas often represent a larger range of groundwater flow field anomalies and need to be monitored and protected.

[0060] Step 104: Divide the foundation stress concentration area into a plurality of grid areas, and receive the deformation rate of each grid area.

[0061] The grid area refers to the square area units of equal size formed by dividing the foundation stress concentration area at regular intervals, and each grid area has independent spatial coordinates and deformation monitoring data.

[0062] The deformation rate refers to the speed at which the surface displacement in the grid area changes over time, expressed as the amount of displacement per unit time. In the embodiments of the present application, the deformation rate can be understood as a quantitative indicator reflecting the deformation trend of the foundation, which is obtained by obtaining the displacement data of a continuous time series through ground settlement monitoring, and is calculated after data preprocessing and time series analysis. The deformation rate is mainly used to quantitatively evaluate the deformation state of the foundation.

[0063] Specifically, a regular grid division method is used to divide each foundation stress concentration area into several square grid areas of equal size. The selection of the grid size (such as 5 meters × 5 meters) needs to comprehensively consider the monitoring accuracy requirements and the actual engineering characteristics. The grid division uses GIS spatial analysis technology to ensure that the grid boundary is accurately connected with the boundary of the foundation stress concentration area. For irregular areas at the boundary, the minimum area principle is used for merging. After the division is completed, the deformation rate data of each grid area is obtained through ground settlement monitoring. The deformation rate reflects the change characteristics of the surface displacement in the grid area over time. In the data acquisition process, the time series analysis method is used to process the original monitoring data, eliminate seasonal effects and outliers, and obtain the rate value reflecting the actual deformation trend. This grid-based refined monitoring method can accurately describe the deformation distribution characteristics inside the foundation stress concentration area and help identify local abnormal areas.

[0064] Step 105: For each grid area, when the deformation rate is greater than the rate threshold, the grid area is determined to be a foundation instability risk area, and early warning information of the foundation instability risk area is generated.

[0065] Among them, the foundation instability risk area refers to the regional unit in the grid area where the deformation rate exceeds the preset rate threshold, indicating that the foundation deformation in this area has reached a level that may affect the safety of the project.

[0066] Specifically, in order to timely discover and warn of the risk of foundation instability, it is necessary to perform threshold judgment on the deformation rate of the grid area and generate warning information. First, a threshold of deformation rate is set (such as -5 mm / day, negative values ​​indicate settlement), which reflects the safety limit of foundation deformation. For each grid area, the measured deformation rate is compared with the threshold. When the absolute value of the deformation rate exceeds the threshold, it indicates that the foundation deformation of the grid area has reached the warning level and needs to be marked as a foundation instability risk area. In order to achieve effective risk warning, the system automatically generates warning information containing information such as the spatial location of the risk area, the deformation rate value, and the degree of over-limit, and sends it to the relevant responsible persons through preset communication channels (such as SMS, email, etc.). The generation of warning information adopts a hierarchical warning mechanism, and different warning levels (such as yellow warning, orange warning, and red warning) are set according to the degree to which the deformation rate exceeds the threshold. Each warning level corresponds to different disposal requirements. This warning method based on threshold judgment can achieve rapid identification and timely response to the risk of foundation instability.

[0067] Based on the above embodiment, as an optional embodiment, in step 101: obtaining groundwater level data of the target area, after this step, the following steps may also be included: Step 501: Acquire meteorological data of the target area, the meteorological data including rainfall amount and rainfall duration; when the rainfall amount is greater than a rainfall amount threshold or the rainfall duration is greater than a duration threshold, issue an early warning to the target area.

[0068] Specifically, in order to prevent the adverse effects of rainfall on foundation stability, it is necessary to dynamically monitor and warn the rainfall characteristics of the target area. The rainfall and rainfall duration data of the target area are obtained in real time through automatic weather stations or meteorological departments. According to engineering experience and local hydrogeological conditions, the rainfall threshold (such as 50 mm cumulative rainfall in 24 hours) and the rainfall duration threshold (such as continuous rainfall for 48 hours) are pre-set. These thresholds reflect the critical rainfall conditions that may cause significant changes in groundwater levels. The system continuously monitors rainfall data. When the 24-hour cumulative rainfall exceeds the rainfall threshold, or the rainfall duration exceeds the duration threshold, it indicates that the rainfall intensity or persistence has reached a level that may affect the stability of the foundation. At this time, the early warning mechanism is automatically triggered. The early warning information contains key information such as rainfall, rainfall duration, and the degree of exceeding the limit, and different early warning levels are set according to the degree of exceeding the limit. This early warning method based on meteorological thresholds can issue early warnings in time before the foundation stability is affected by rainfall. By monitoring and warning rainfall characteristics, the risk of foundation instability induced by rainfall can be effectively prevented, and sufficient response time can be reserved for taking protective measures.

[0069] Reference Figure 2, is a foundation stability analysis and early warning system based on groundwater level changes provided by an embodiment of the present application, the system comprises: a data acquisition module, a gradient change rate determination module, an area determination module, and a foundation analysis and early warning module, wherein: A data acquisition module is used to acquire groundwater level data of a target area, and to construct a spatiotemporal isosurface of groundwater level in the target area based on the water level value, sampling time and sampling location in the groundwater level data; A gradient change rate determination module is used to determine multiple monitoring points on the groundwater level time-space isosurface and calculate the groundwater level gradient change rate of each monitoring point in the groundwater level time-space isosurface; The area determination module is used to determine the foundation stress concentration area in the spatiotemporal isosurface of the groundwater level according to the gradient change rate of each groundwater level; The foundation analysis and early warning module is used to divide the foundation stress concentration area into multiple grid areas and receive the deformation rate of each grid area; for each grid area, when the deformation rate is greater than the rate threshold, the grid area is determined to be a foundation instability risk area, and early warning information of the foundation instability risk area is generated.

[0070] It should be noted that: when the device provided in the above embodiment realizes its function, only the division of the above functional modules is used as an example. In actual application, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the device is divided into different functional modules to complete all or part of the functions described above. In addition, the device and method embodiments provided in the above embodiment belong to the same concept, and the specific implementation process is detailed in the method embodiment, which will not be repeated here.

[0071] The present application also discloses an electronic device. Figure 3 , Figure 3 The electronic device 300 may include: at least one processor 301 , at least one network interface 304 , a user interface 303 , a memory 305 , and at least one communication bus 302 .

[0072] The communication bus 302 is used to realize the connection and communication between these components.

[0073] The user interface 303 may include a display interface and a camera interface. Optionally, the user interface 303 may also include a standard wired interface and a wireless interface.

[0074] The network interface 304 may optionally include a standard wired interface or a wireless interface (such as a WI-FI interface).

[0075] Among them, the processor 301 may include one or more processing cores. The processor 301 uses various interfaces and lines to connect various parts in the entire server, and executes various functions of the server and processes data by running or executing instructions, programs, code sets or instruction sets stored in the memory 305, and calling data stored in the memory 305. Optionally, the processor 301 can be implemented in at least one hardware form of digital signal processing (Digital Signal Processing, DSP), field programmable gate array (Field-Programmable Gate Array, FPGA), and programmable logic array (Programmable Logic Array, PLA). The processor 301 can integrate one or a combination of a central processing unit (Central Processing Unit, CPU), a graphics processing unit (Graphics Processing Unit, GPU) and a modem. Among them, the CPU mainly processes the operating system, user interface diagrams and applications, etc.; the GPU is responsible for rendering and drawing the content to be displayed on the display screen; the modem is used to process wireless communications. It can be understood that the above-mentioned modem may not be integrated into the processor 301, and it can be implemented separately through a chip.

[0076] Among them, the memory 305 may include a random access memory (Random Access Memory, RAM) and may also include a read-only memory (Read-Only Memory). Optionally, the memory 305 includes a non-transitory computer-readable storage medium. The memory 305 can be used to store instructions, programs, codes, code sets or instruction sets. The memory 305 may include a program storage area and a data storage area, wherein the program storage area may store instructions for implementing an operating system, instructions for at least one function (such as a touch function, a sound playback function, an image playback function, etc.), instructions for implementing the above-mentioned various method embodiments, etc.; the data storage area may store data involved in the above-mentioned various method embodiments, etc. The memory 305 may optionally also be at least one storage device located away from the aforementioned processor 301. Refer to Figure 3 The memory 305 as a computer storage medium may include an operating system, a network communication module, a user interface module, and an application program of a foundation stability analysis and early warning method based on groundwater level changes.

[0077] exist Figure 3In the electronic device 300 shown, the user interface 303 is mainly used to provide an input interface for the user and obtain the data input by the user; and the processor 301 can be used to call the memory 305 to store an application program for a foundation stability analysis and early warning method based on groundwater level changes. When executed by one or more processors 301, the electronic device 300 executes one or more methods in the above-mentioned embodiments. It should be noted that for the aforementioned method embodiments, for the sake of simple description, they are all expressed as a series of action combinations, but those skilled in the art should be aware that the present application is not limited by the described order of actions, because according to the present application, certain steps can be performed in other orders or simultaneously. Secondly, those skilled in the art should also be aware that the embodiments described in the specification are all preferred embodiments, and the actions and modules involved are not necessarily required for the present application.

[0078] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0079] In the several embodiments provided in this application, it should be understood that the disclosed devices can be implemented in other ways. For example, the device embodiments described above are only schematic, such as the division of units, which is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some service interfaces, and the indirect coupling or communication connection of devices or units can be electrical or other forms.

[0080] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0081] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0082] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable memory. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, which is stored in a memory and includes several instructions for a computer device (which can be a personal computer, server or network device, etc.) to execute all or part of the steps of the various embodiments of the present application. The aforementioned memory includes: various media that can store program codes, such as USB flash drives, mobile hard drives, magnetic disks or optical disks.

[0083] The above are only exemplary embodiments of the present disclosure and cannot be used to limit the scope of the present disclosure. That is, any equivalent changes and modifications made according to the teachings of the present disclosure are still within the scope of the present disclosure. After considering the disclosure of the specification and practice, those skilled in the art will easily think of other embodiments of the present disclosure.

[0084] This application is intended to cover any variation, use or adaptation of the present disclosure, which follows the general principles of the present disclosure and includes common knowledge or customary technical means in the art not recorded in the present disclosure. The description and examples are to be regarded as exemplary only.

Claims

1. A foundation stability analysis and early warning method based on groundwater level changes, characterized in that: include: Acquire groundwater level data of a target area, and construct a spatiotemporal isosurface of groundwater level in the target area based on the water level value, sampling time and sampling location in the groundwater level data; Determine a plurality of monitoring points on the groundwater level time-space isosurface, and calculate the groundwater level gradient change rate of each of the monitoring points in the groundwater level time-space isosurface; Determining the foundation stress concentration area in the groundwater level time-space isosurface according to each of the groundwater level gradient change rates; Dividing the foundation stress concentration area into a plurality of grid areas, and receiving the deformation rate of each of the grid areas; For each of the grid areas, when the deformation rate is greater than a rate threshold, the grid area is determined to be a foundation instability risk area, and early warning information of the foundation instability risk area is generated.

2. The foundation stability analysis and early warning method based on groundwater level changes according to claim 1 is characterized in that: The step of constructing a spatiotemporal isosurface of groundwater level in the target area based on the water level value, sampling time and sampling location in the groundwater level data comprises: Dividing the target area into a plurality of spatial units with preset spatial specifications; Based on the water level value, sampling time and sampling location in the groundwater level data, determining the target water level value, target sampling time and target sampling location corresponding to each of the spatial units; Performing interpolation processing on the sampling time, target water level value and target sampling position of each of the spatial units to obtain groundwater level distribution characteristics at multiple time points; The groundwater level distribution data at each of the time points are sequentially connected in chronological order to generate a spatiotemporal isosurface of the groundwater level in the target area.

3. The foundation stability analysis and early warning method based on groundwater level changes according to claim 1 is characterized in that: The step of determining a plurality of monitoring points on the spatiotemporal contour surface of the groundwater level comprises: Selecting multiple time sections on the groundwater level time-space isosurface; The terrain slope mutation area, the building distribution area and the groundwater level mutation area in each of the time sections are respectively used as point areas on the groundwater level spatiotemporal isosurface; According to the location of each of the point areas, a plurality of monitoring points on the spatiotemporal isosurface of the groundwater level are determined.

4. The foundation stability analysis and early warning method based on groundwater level changes according to claim 1 is characterized in that: The step of calculating the groundwater level gradient change rate at each monitoring point in the groundwater level spatiotemporal isosurface comprises: Obtaining the water level values ​​of each monitoring point in the groundwater level spatiotemporal isosurface at multiple sampling times; For each of the monitoring points, based on the water level values ​​of the monitoring point at multiple sampling moments, the relative water level difference between the monitoring point and each monitoring point within a preset radius is calculated; Determine the water level gradient of the monitoring point according to each of the relative water level differences and the relative distances between the monitoring point and each of the monitoring points within a preset radius; The changing rate of the water level gradient at each sampling time is taken as the groundwater level gradient changing rate of the monitoring point.

5. The foundation stability analysis and early warning method based on groundwater level changes according to claim 1 is characterized in that: Determining the foundation stress concentration area in the groundwater level time-space isosurface according to each groundwater level gradient change rate includes: Calculating the time cumulative value of the gradient change rate of each monitoring point, and selecting the monitoring points whose time cumulative value is greater than the cumulative threshold as target monitoring points; Taking each target monitoring point as a starting point, expanding outward by a preset distance to form an expansion area of ​​each starting point; In each of the expansion regions, a plurality of overlapping expansion regions are merged into one candidate stress region to obtain a plurality of candidate stress regions; The candidate stress regions whose area is greater than the area threshold are taken as the foundation stress concentration regions in the spatiotemporal contour surface of the groundwater level.

6. The foundation stability analysis and early warning method based on groundwater level changes according to claim 5 is characterized in that: The calculating the time accumulation value of the gradient change rate of each monitoring point comprises: For each of the monitoring points, dividing the gradient change rate into sub-gradient change rates corresponding to a plurality of time windows; Performing a differential operation on the sub-gradient change rates of adjacent time windows to obtain the change amount between the time windows; The products of each of the change amounts and the duration of the time window are summed to obtain a time-accumulated value of the gradient change rate of the monitoring point.

7. The foundation stability analysis and early warning method based on groundwater level changes according to claim 1 is characterized in that: After obtaining the groundwater level data of the target area, the method further includes: Acquiring meteorological data of the target area, wherein the meteorological data includes rainfall amount and rainfall duration; When the rainfall amount is greater than a rainfall amount threshold or the rainfall duration is greater than a duration threshold, an early warning is issued to the target area.

8. A foundation stability analysis and early warning system based on groundwater level changes, characterized in that: The system comprises: A data acquisition module, used to acquire groundwater level data of a target area, and construct a spatiotemporal isosurface of groundwater level in the target area based on the water level value, sampling time and sampling location in the groundwater level data; A gradient change rate determination module is used to determine a plurality of monitoring points on the groundwater level time-space isosurface, and calculate the groundwater level gradient change rate of each of the monitoring points in the groundwater level time-space isosurface; An area determination module, used to determine the foundation stress concentration area in the groundwater level time-space isosurface according to each groundwater level gradient change rate; The foundation analysis and early warning module is used to divide the foundation stress concentration area into multiple grid areas and receive the deformation rate of each grid area; for each grid area, when the deformation rate is greater than a rate threshold, the grid area is determined to be a foundation instability risk area, and early warning information of the foundation instability risk area is generated.

9. An electronic device, characterized in that: It includes a processor, a memory, a user interface and a network interface, the memory is used to store instructions, the user interface and the network interface are used to communicate with other devices, and the processor is used to execute the instructions stored in the memory so that the electronic device executes the foundation stability analysis and early warning method based on groundwater level changes as described in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores instructions, and when the instructions are executed, the foundation stability analysis and early warning method based on groundwater level changes as described in any one of claims 1 to 7 is executed.

Images