Post-construction reconstruction method of landslide dam based on physical mechanism and survey monitoring data

Abstract

The present application relates to a method, device, and computer equipment for post-construction reconstruction of a landslide dam based on physical mechanisms and survey and monitoring data. The method comprises: establishing a landslide dam model of the landslide dam, and obtaining a material parameter set, a borehole data set, and a real dam body deformation data set of the landslide dam; determining the stratum boundary variation range in the landslide dam model based on the material parameter set, the borehole data set, and a finite element algorithm, and determining an initial stratum variation curve within the stratum boundary variation range; predicting dam body deformation based on the initial stratum variation curve and the finite element algorithm, obtaining a predicted dam body deformation data set, and calculating the loss between the predicted dam body deformation data set and the real dam body deformation data set until the loss meets a preset loss condition; and determining the initial stratum variation curve corresponding to the current round of predicted dam body deformation data set as the stratum variation curve. The present method can improve the accuracy of the post-construction reconstruction method of a landslide dam based on physical mechanisms and survey and monitoring data.

Description

Technical Field

[0001] The present application relates to the field of data processing technology, and in particular to a method, device, computer equipment and computer-readable storage medium for post-construction reconstruction of a landslide dam based on physical mechanisms and survey and monitoring data. Background Art

[0002] A landslide dam is a natural dam formed by a geological disaster blocking a river. Its stability is directly related to the safety of downstream areas. To ensure the safety of downstream areas, it is necessary to construct the stratigraphic structure of the landslide dam and predict whether the landslide dam is in a safe state based on the stratigraphic structure of the landslide dam.

[0003] Traditionally, limited borehole data is collected from each borehole in a landslide dam. For two adjacent boreholes, the stratigraphic structure of the boreholes is determined based on the data. Geostatistical interpolation methods are then used to predict the stratigraphic curves for the area between the two adjacent boreholes. The stratigraphic structure of the landslide dam is then constructed based on these stratigraphic curves and the stratigraphic structures of the boreholes.

[0004] However, due to the complex geological structure of the dam, it is difficult to accurately predict the stratigraphic structure of the area between two adjacent boreholes using only geostatistical interpolation methods. Therefore, the accuracy of current dam construction methods is relatively low. Summary of the Invention

[0005] Based on this, it is necessary to provide a method, device, computer equipment and computer-readable storage medium for post-construction reconstruction of a landslide dam based on physical mechanisms and survey and monitoring data to address the above technical problems.

[0006] In a first aspect, the present application provides a post-construction reconstruction method for a landslide dam based on physical mechanisms and survey and monitoring data, comprising:

[0007] Establishing a landslide dam model of the landslide dam, and obtaining a material parameter set, a drilling data set, and a real dam body deformation data set of the landslide dam;

[0008] Determining a stratum boundary variation range in the landslide dam model based on the material parameter set, the drilling data set, and a finite element algorithm, and determining an initial stratum variation curve within the stratum boundary variation range;

[0009] dam deformation prediction is performed based on the initial stratum change curve and the finite element algorithm to obtain a predicted dam deformation data set, and a loss between the predicted dam deformation data set and the actual dam deformation data set is calculated until the loss satisfies a preset loss condition;

[0010] The initial stratum change curve corresponding to the predicted dam deformation data set of the current round is determined as the stratum change curve.

[0011] In one embodiment, determining the stratum boundary variation range in the landslide dam model based on the material parameter set, the drilling data set, and a finite element algorithm includes:

[0012] Determining an affected area in the landslide dam model based on the material parameter set, the finite element algorithm, and a preset first difference threshold value, wherein the first difference threshold value is a stress change threshold value of the anti-seepage structure in the landslide dam;

[0013] The range of stratum boundary changes in the affected area is determined based on the drilling data set, the material parameter set, a preset second difference threshold and the finite element algorithm; the second difference threshold is a stress change threshold of the anti-seepage structure in the landslide dam.

[0014] In one embodiment, the landslide dam includes an anti-seepage structure, and determining the affected area in the landslide dam model based on the material parameter set, the finite element algorithm, and a preset first difference threshold includes:

[0015] Dividing the landslide dam model into an initial impact area and an initial non-impact area;

[0016] determining a first stress of the anti-seepage structure according to the material parameter set, the initial affected area, the initial non-affected area, and the finite element algorithm;

[0017] Adjusting the initial impact area and the initial non-impact area, and determining a second stress of the anti-seepage structure according to the material parameter set, the adjusted initial impact area, the adjusted initial non-impact area, and the finite element algorithm;

[0018] An affected area in the landslide dam model is determined based on the first stress, the second stress, and a preset first difference threshold.

[0019] In one embodiment, determining the affected area in the landslide dam model based on the first stress, the second stress, and a preset first difference threshold includes:

[0020] performing difference processing on the second stress and the first stress to obtain a first difference, and determining whether the first difference exceeds a preset first difference threshold;

[0021] If the first difference exceeds the first difference threshold, determining the adjusted initial impact area as the impact area;

[0022] If the first difference does not exceed the first difference threshold, updating the first stress according to the second stress;

[0023] According to the adjusted initial impact area and the adjusted initial non-impact area, the step of adjusting the initial impact area and the initial non-impact area is performed until the first difference exceeds the first difference threshold, and the adjusted initial impact area is determined as the impact area.

[0024] In one embodiment, determining the range of stratum boundary variation in the affected area based on the drilling data set, the material parameter set, a preset second difference threshold, and the finite element algorithm includes:

[0025] determining an initial stratigraphic boundary variation range in the affected area, and determining a stratigraphic demarcation point in the initial stratigraphic boundary variation range based on the drill hole dataset;

[0026] generating an initial stratum boundary curve based on the stratum boundary point and the initial stratum boundary variation range;

[0027] A stratum boundary variation range is determined within the initial stratum boundary variation range according to the initial stratum boundary curve, the affected area, the material parameter set, a preset second difference threshold, and the finite element algorithm.

[0028] In one embodiment, determining the formation boundary variation range within the initial formation boundary variation range based on the initial formation boundary curve, the affected area, the material parameter set, a preset second difference threshold, and the finite element algorithm includes:

[0029] determining a third stress of the anti-seepage structure in the landslide dam according to the initial stratum boundary curve, the affected area, the material parameter set, and a finite element algorithm;

[0030] Adjusting the initial stratum boundary curve toward the stratum boundary point, and determining the fourth stress of the anti-seepage structure based on the adjusted initial stratum boundary curve, the affected area, the material parameter set, and a finite element algorithm;

[0031] Based on the third stress, the fourth stress, the adjusted initial formation boundary curve and a preset second difference threshold, a formation boundary variation range is determined within the initial formation boundary variation range.

[0032] In one embodiment, determining the formation boundary variation range within the initial formation boundary variation range based on the third stress, the fourth stress, the adjusted initial formation boundary curve, and a preset second difference threshold value includes:

[0033] performing difference processing on the third stress and the fourth stress to obtain a second difference, and determining whether the second difference exceeds a preset second difference threshold;

[0034] If the second difference exceeds the second difference threshold, determining a formation boundary variation range within the initial formation boundary variation range based on the adjusted initial formation boundary curve;

[0035] If the second difference does not exceed the second difference threshold, updating the third stress based on the fourth stress;

[0036] The step of adjusting the initial stratigraphic boundary curve toward the stratigraphic boundary point is performed according to the adjusted initial stratigraphic boundary curve until the second difference exceeds the second difference threshold, and the stratigraphic boundary change range is determined within the initial stratigraphic boundary change range based on the adjusted initial stratigraphic boundary curve.

[0037] In one embodiment, the initial formation change curve includes various control parameters, and the dam deformation prediction is performed based on the initial formation change curve and the finite element algorithm to obtain a predicted dam deformation dataset, and the loss between the predicted dam deformation dataset and the actual dam deformation dataset is calculated until the loss meets a preset loss condition, including:

[0038] Optimizing the initial stratum change curve according to the control parameters, the finite element algorithm, and the optimization algorithm to obtain an optimized stratum change curve, and performing dam body prediction according to the optimized stratum change curve to obtain a predicted dam body deformation data set;

[0039] Calculating the loss between the predicted dam deformation dataset and the actual dam deformation dataset, and determining whether the loss exceeds a preset loss threshold;

[0040] If the loss does not exceed the loss threshold, determining that the loss meets the preset loss condition;

[0041] If the loss exceeds the loss threshold, the step of optimizing the initial formation change curve according to the control parameters, the finite element algorithm and the optimization algorithm is executed until the loss does not exceed the loss threshold, and it is determined that the loss meets the preset loss condition.

[0042] In one embodiment, optimizing the initial formation change curve according to the control parameters, the finite element algorithm, and the optimization algorithm to obtain the optimized formation change curve includes:

[0043] Determining a disturbance value of each of the control parameters, and perturbing each of the control parameters according to the disturbance value of each of the control parameters to generate the initial formation change curve after disturbance;

[0044] Based on the finite element algorithm and the initial stratum change curve after disturbance, dam body deformation prediction is performed to obtain a disturbed dam body deformation data set;

[0045] Calculating a disturbance error between the disturbance dam body deformation data set and the real dam body deformation data set;

[0046] A disturbance gradient is calculated according to the disturbance error, and each of the control parameters is optimized according to the disturbance gradient and an optimization algorithm to generate an optimized initial formation change curve.

[0047] In a second aspect, the present application further provides a post-construction reconstruction device for a landslide dam based on physical mechanisms and survey and monitoring data, comprising:

[0048] an acquisition module, configured to establish a landslide dam model of the landslide dam and acquire a material parameter set, a drilling data set, and a real dam body deformation data set of the landslide dam;

[0049] a first determining module, configured to determine a stratum boundary variation range in the landslide dam model based on the material parameter set, the drilling data set, and a finite element algorithm, and to determine an initial stratum variation curve within the stratum boundary variation range;

[0050] a prediction module, configured to predict dam deformation based on the initial stratum change curve and the finite element algorithm, obtain a predicted dam deformation data set, and calculate a loss between the predicted dam deformation data set and the actual dam deformation data set until the loss satisfies a preset loss condition;

[0051] The second determining module is configured to determine the initial stratum change curve corresponding to the predicted dam deformation data set in the current round as the stratum change curve.

[0052] In a third aspect, the present application provides a computer device including a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps in the above-mentioned method embodiments when executing the computer program.

[0053] In a fourth aspect, the present application provides a computer-readable storage medium having a computer program stored thereon, which implements the steps in the above-mentioned method embodiments when executed by a processor.

[0054] In a fifth aspect, the present application provides a computer program product, including a computer program, which implements the steps in the above-mentioned method embodiments when executed by a processor.

[0055] The above-mentioned post-construction reconstruction method, apparatus, computer device, and computer-readable storage medium for a landslide dam based on physical mechanisms and survey and monitoring data establishes a landslide dam model for the landslide dam and obtains a material parameter set, a borehole dataset, and an actual dam deformation dataset for the landslide dam. Based on the material parameter set, the borehole dataset, and a finite element algorithm, the stratum boundary variation range in the landslide dam model is determined, and an initial stratum variation curve is determined within the stratum boundary variation range. Dam deformation is predicted based on the initial stratum variation curve and the finite element algorithm to obtain a predicted dam deformation dataset, and the loss between the predicted dam deformation dataset and the actual dam deformation dataset is calculated until the loss meets a preset loss condition. The initial stratum variation curve corresponding to the predicted dam deformation dataset for the current round is determined as the stratum variation curve. Using this method, the stratum boundary variation range of the landslide dam is initially determined using the material parameter set, the borehole dataset, and the finite element algorithm. Then, the initial stratigraphic change curve is iteratively optimized with the help of the dam deformation dataset, so that the initial stratigraphic change curve continuously approaches the true value until the loss between the predicted dam deformation dataset and the true dam deformation dataset meets the preset loss condition. An accurate stratigraphic change curve is obtained, thereby improving the accuracy of the landslide dam stratigraphic reconstruction method. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] In order to more clearly illustrate the technical solutions in the embodiments of the present application or related technologies, the following briefly introduces the drawings required for use in the embodiments of the present application or related technical descriptions. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying any creative work.

[0057] Figure 1 1 is a flow chart of a post-construction reconstruction method for a landslide dam based on physical mechanisms and survey and monitoring data in one embodiment;

[0058] Figure 2 A schematic diagram of a process for determining a range of formation boundary variation in one embodiment;

[0059] Figure 3 A schematic diagram of a process for determining an impact area in one embodiment;

[0060] Figure 4 A schematic diagram of a process for determining an affected area according to a first difference in one embodiment;

[0061] Figure 5 A schematic diagram of stratum boundary influence level zoning in one embodiment;

[0062] Figure 6A flowchart for determining impact zones in a landslide dam model in an exemplary embodiment;

[0063] Figure 7 A schematic diagram of a process for determining a range of formation boundary changes in an impact area in one embodiment;

[0064] Figure 8 A schematic diagram of a flow chart of determining a formation boundary variation range within an initial formation boundary variation range in one embodiment;

[0065] Figure 9 Schematic diagram of a process for determining a formation boundary variation range based on a second difference in one embodiment;

[0066] Figure 10 A schematic diagram of a process for determining a range of formation boundary changes in an exemplary embodiment;

[0067] Figure 11 A schematic diagram of a process for performing dam body prediction in one embodiment;

[0068] Figure 12 A schematic diagram of a process for optimizing an initial formation change curve in one embodiment;

[0069] Figure 13 A schematic diagram of a flow chart for determining a formation change curve according to an exemplary embodiment;

[0070] Figure 14 is a flow chart of a post-construction reconstruction method of a landslide dam based on physical mechanisms and survey and monitoring data in an exemplary embodiment;

[0071] Figure 15 1. A structural block diagram of a post-construction reconstruction device for a landslide dam based on physical mechanisms and survey and monitoring data in one embodiment;

[0072] Figure 16 FIG. 1 is a diagram showing the internal structure of a computer device in one embodiment. DETAILED DESCRIPTION

[0073] In order to make the purpose, technical solutions and advantages of this application more clear, the following further describes this application in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain this application and are not intended to limit this application.

[0074] A landslide dam is a natural dam formed by a geological disaster blocking a river. Its stability is directly related to the safety of downstream areas. To ensure the safety of downstream areas, it is necessary to construct the stratigraphic structure of the landslide dam and predict whether the landslide dam is in a safe state based on the stratigraphic structure of the landslide dam.

[0075] Traditionally, limited borehole data is collected from each borehole in a landslide dam. For two adjacent boreholes, the stratigraphic structure of the boreholes is determined based on the data. Geostatistical interpolation methods are then used to predict the stratigraphic curves for the area between the two adjacent boreholes. The stratigraphic structure of the landslide dam is then constructed based on these stratigraphic curves and the stratigraphic structures of the boreholes.

[0076] However, in traditional technologies, due to the complex geological structure of the landslide dam, it is difficult to accurately predict the stratigraphic structure of the area between two adjacent boreholes by relying solely on geostatistical interpolation methods. Therefore, the accuracy of current landslide dam construction methods is low. In addition, the borehole data space is sparse, and the stratigraphic structure determined based on borehole data and geostatistical interpolation methods is inaccurate. In addition, the number of monitoring points in the real dam deformation dataset is small and most of them are located on the dam surface. Due to the limitation of sensor accuracy, the real dam deformation dataset has certain errors. Relying on the real dam deformation dataset to invert the stratigraphic structure is prone to multiple solution problems, that is, there are multiple stratigraphic curves that can fit the monitoring data well, and relying on the real dam deformation dataset to invert the stratigraphic structure does not take into account the response of the anti-seepage structure.

[0077] Therefore, the present application provides a post-construction reconstruction method for a landslide dam based on physical mechanisms and survey and monitoring data. The stratum boundary variation range of the landslide dam is initially determined using a material parameter set, a borehole dataset, and a finite element algorithm. Then, the initial stratum variation curve is iteratively optimized using the dam deformation dataset, so that the initial stratum variation curve continuously approaches the true value until the loss between the predicted dam deformation dataset and the actual dam deformation dataset meets a preset loss condition. An accurate stratum variation curve is obtained, thereby improving the accuracy of the post-construction reconstruction method for landslide dams based on physical mechanisms and survey and monitoring data. Furthermore, the present application determines the stratum variation range of the landslide dam based on the stress changes of the anti-seepage structure in the landslide dam, the material dataset, the borehole dataset, and the finite element dataset. This solves the problem of multiple solutions that occur when inverting the stratum structure based on the actual dam deformation dataset in the prior art. By using the stratum variation curve within the stratum variation range, not only can the actual dam deformation dataset obtained from the actual survey be fitted, but the stratum variation curve is also sensitive to the response of the anti-seepage structure, further improving the accuracy of the stratum variation curve.

[0078] In one embodiment, Figure 1 As shown, a method for post-construction reconstruction of a landslide dam based on physical mechanisms and survey and monitoring data is provided. This embodiment of the application takes the application of this method to a computer device as an example for explanation. This embodiment of the application does not limit the execution device of the method for post-construction reconstruction of a landslide dam based on physical mechanisms and survey and monitoring data. The method includes the following steps 102 to 108:

[0079] Step 102 : establishing a landslide dam model of the landslide dam, and obtaining a material parameter set, a drilling data set, and a real dam body deformation data set of the landslide dam.

[0080] During implementation, the computer device acquires a property dataset of the landslide dam and establishes a landslide dam model based on the property dataset of the landslide dam. Simultaneously, the computer device acquires a drilling dataset, a material parameter dataset, and a real dam body deformation dataset of the landslide dam.

[0081] Specifically, the computer device acquires a property dataset for the landslide dam, which includes topographic data, location data, and dimensional data. The topographic data characterizes the detailed relief of the landslide dam's terrain. The location data includes, but is not limited to, coordinate data and contour data for the landslide dam. Furthermore, since the landslide dam is also equipped with an anti-seepage structure, the computer device acquires a property dataset for the anti-seepage structure. This property dataset includes the location and dimensional data for the anti-seepage structure. The computer device uses a finite element algorithm, the property datasets for the anti-seepage structure, and the landslide dam's property datasets to build a landslide dam model. Simultaneously, the computer device acquires a user-entered borehole dataset for the landslide dam. This borehole dataset is obtained by the user conducting on-site field surveys and drilling holes into the landslide dam at preset distances, thereby collecting drilling data for each hole. The borehole dataset reflects the stratigraphic structure of the landslide dam at the locations of the holes, namely, the types and locations of the geological strata present in the landslide dam at the locations of the holes. Sensors monitor the displacement of the dam body, generating a true dam deformation dataset. The computer device obtains a first material parameter set and a second material parameter set input by a user, wherein the material parameters in the first material parameter set are material parameters of the landslide dam impact area, and the material parameters in the second material parameter set are material parameters of the landslide dam non-impact area.

[0082] Optionally, the landslide dam model can be, but is not limited to, a three-dimensional model of the landslide dam or a two-dimensional model of the landslide dam. If the user has high requirements for the accuracy of the stratigraphic structure of the landslide dam, a three-dimensional model of the landslide dam can be established. If the user has high requirements for the efficiency of the stratigraphic structure of the landslide dam, a two-dimensional model of the landslide dam can be established, that is, the dimensions of the established landslide dam model are determined according to the construction requirements of the landslide dam stratigraphic structure. The embodiment of the present application does not limit the dimensions of the landslide dam model. The finite element algorithm is implemented by finite element software, which can be, but is not limited to, Abaqus software (a powerful finite element software for engineering simulation). The embodiment of the present application does not limit the finite element software.

[0083] Step 104 : determining the stratum boundary variation range in the landslide dam model based on the material parameter set, the borehole data set, and the finite element algorithm, and determining an initial stratum variation curve in the stratum boundary variation range.

[0084] The range of stratum boundary variation is the envelope of the stratum boundary, which is also the interval where the stratum boundary is located.

[0085] During implementation, the computer determines the impact zone within the landslide dam model based on a material parameter set, a finite element method, and a preset first difference threshold. It also determines the stratigraphic boundary variation range within the impact zone based on the finite element method, a borehole dataset, and a preset second difference threshold. The computer's curve generation algorithm (a Python algorithm) then generates an initial stratigraphic variation curve within the stratigraphic boundary variation range.

[0086] Specifically, to reinforce the landslide dam, an anti-seepage structure is installed in the landslide dam. Based on the material parameter set, the finite element algorithm, and the first difference threshold, the computer device determines the dam body area that has a greater impact on the stress response of the anti-seepage structure in the landslide dam model, and determines this dam body area as the affected area. The computer device then determines the initial formation boundary variation range in the affected area, and based on the drilling data set, the finite element algorithm, and the preset second difference threshold, narrows the initial formation boundary variation range to obtain the formation boundary variation range. The computer device determines the various control parameters based on the formation boundary variation range and simulation requirements, and generates the initial formation boundary variation curve based on the various control parameters and the curve generation algorithm.

[0087] In one exemplary embodiment, a computer device divides a landslide dam model into an initial impact area and an initial non-impact area, and determines a first stress of the anti-seepage structure based on a material parameter set, the initial impact area, the initial non-impact area, and a finite element method. The computer device then adjusts the initial impact area and the initial non-impact area to obtain adjusted initial impact area and initial non-impact area. The computer device determines a second stress of the anti-seepage structure based on the adjusted initial impact area, the adjusted initial non-impact area, the material parameter set, and the finite element method. The computer device calculates a first difference between the first stress and the second stress until the first difference meets a preset first difference threshold. The computer device then determines the adjusted initial impact area for the current round as the impact area. The computer device determines an initial formation boundary variation range within the impact area and generates an initial formation demarcation curve within the initial formation boundary variation range based on the borehole dataset. The computer device narrows the initial formation boundary variation range based on the initial formation demarcation curve, the impact area, the material parameter set, the second difference threshold, and the finite element method to obtain a formation boundary variation range. The computer device determines various control parameters according to the formation boundary variation range and simulation requirements, and generates an initial formation boundary variation curve according to the various control parameters.

[0088] Optionally, the number of control parameters is determined based on the required accuracy of the stratum structure of the landslide dam. A greater number of control parameters results in a more accurate simulation of the stratum boundary variation curve and a higher accuracy of the stratum structure of the landslide dam. A smaller number of control parameters results in a coarser simulation of the stratum boundary variation curve and a lower accuracy of the stratum structure of the landslide dam.

[0089] Step 106 , based on the initial formation change curve and the finite element algorithm, dam deformation is predicted to obtain a predicted dam deformation data set, and the loss between the predicted dam deformation data set and the actual dam deformation data set is calculated until the loss meets the preset loss condition.

[0090] The initial stratum change curve contains various control parameters. The predicted dam deformation data sets that do not meet the loss conditions are all predicted by the initial stratum change curve.

[0091] During implementation, a loss condition is pre-set in the computer device. The computer device optimizes the initial formation change curve according to various control parameters, and predicts the dam deformation based on the optimized initial formation change curve and the finite element algorithm to obtain a predicted dam deformation data set. The computer device calculates the loss between the predicted dam deformation data set and the actual dam deformation data set, and determines whether the loss meets the preset loss condition. If the loss meets the preset loss condition, the computer device continues to execute the following step 108. The specific processing process of step 108 is described in detail in the following embodiment, and the embodiment of this application will not be repeated here. If the loss does not meet the loss condition, the computer device continues to optimize the initial formation change curve until the loss meets the loss condition.

[0092] Specifically, the loss condition pre-set in the computer device is that the error does not exceed the error threshold. The computer device determines each control parameter and optimizes the initial formation change curve based on each control parameter, the finite element algorithm, and the optimization algorithm to obtain an optimized initial formation change curve. The computer device then predicts dam deformation based on the optimized initial formation change curve and the finite element algorithm to obtain a predicted dam deformation dataset and calculates the error between the predicted dam deformation dataset and the actual dam deformation dataset. This error is also the loss between the predicted dam deformation dataset and the actual dam deformation dataset. The computer device then determines whether the error exceeds the error threshold. The computer device determines that the error does not meet the loss condition. If the error does not meet the loss condition, the computer device executes the steps of optimizing each control parameter until the loss meets the loss condition, and the computer device obtains the optimized initial formation change curve.

[0093] Step 108: Determine the initial stratum change curve corresponding to the current round of predicted dam deformation data set as the stratum change curve.

[0094] During implementation, the computer determines the optimized initial stratum change curve corresponding to the current round of predicted dam deformation dataset as the stratum change curve. The computer then constructs the stratum structure of the landslide dam based on the borehole dataset, the stratum change curve, and the landslide dam model.

[0095] Specifically, the computer determines the optimized initial stratigraphic change curve corresponding to the current round of predicted dam deformation dataset as the stratigraphic change curve. The computer then determines the two geological layers corresponding to the stratigraphic change curve based on the borehole dataset and constructs the stratigraphic structure of the landslide dam based on the two geological layers, the stratigraphic change curve, and the landslide dam model.

[0096] In this post-construction reconstruction method for landslide dams based on physical mechanisms and survey and monitoring data, the dam's stratigraphic boundary variation range is initially determined using a material parameter set, a borehole dataset, and a finite element algorithm. The initial stratigraphic variation curve is then iteratively optimized using the dam deformation dataset, continuously approaching the true value until the loss between the predicted and actual dam deformation datasets meets a preset loss condition. This results in an accurate stratigraphic variation curve, thereby improving the accuracy of the landslide dam stratigraphic reconstruction method.

[0097] In an exemplary embodiment, Figure 2 As shown, the specific process of determining the range of stratum boundary variation in the landslide dam model based on the material parameter set, the borehole data set and the finite element algorithm in step 104 includes steps 202 to 204.

[0098] Step 202 : determining an impact area in the landslide dam model based on a material parameter set, a finite element algorithm, and a preset first difference threshold.

[0099] The first difference threshold is a stress change threshold of the anti-seepage structure in the landslide dam. The material parameter set includes a first material parameter set for the affected area and a second material parameter set for the non-affected area.

[0100] In practice, the affected area is the region within the dam body that significantly affects the stress response of the anti-seepage structure, while the unaffected area is the region within the dam body that has a less significant impact on the stress response of the anti-seepage structure. Since both the affected and unaffected areas are related to the stress of the anti-seepage structure, this application uses the stress change of the anti-seepage structure (a first difference threshold) to verify the accuracy of the affected and unaffected areas. Specifically, a computer device divides the dam model to obtain an initial affected area and an initial unaffected area. The computer device then determines the first stress of the anti-seepage structure within the dam body based on the material parameter set, the finite element method, the initial affected area, and the initial unaffected area. The computer device reduces the initial affected area to obtain an adjusted initial affected area, and adjusts the initial unaffected area based on the adjusted initial affected area to obtain an adjusted initial unaffected area. The computer device then determines the second stress of the anti-seepage structure within the dam body based on the material parameter set, the finite element method, the adjusted initial affected area, and the adjusted initial unaffected area. The computer device then determines a first difference between the first stress and the second stress until the first difference meets the first difference threshold. The computer device determines the adjusted initial impact area of ​​the current round as the impact area.

[0101] Step 204 : determining the range of stratum boundary variation in the affected area based on the drilling data set, the material parameter set, the preset second difference threshold, and the finite element algorithm.

[0102] The second difference threshold is the stress change threshold of the anti-seepage structure in the landslide dam.

[0103] During implementation, the computer determines an initial stratigraphic boundary variation range within the affected area and generates an initial stratigraphic boundary curve based on the initial stratigraphic boundary variation range and the drill hole dataset. The computer then reduces the initial stratigraphic boundary variation range based on the initial stratigraphic boundary curve, the affected area, the material parameter set, the second difference threshold, and a finite element algorithm to obtain the stratigraphic boundary variation range.

[0104] In this embodiment, the stratum boundary variation range is determined within the landslide dam model by using a material dataset, a borehole dataset, a finite element algorithm, and a stress change threshold for the anti-seepage structure. This provides constraints for the subsequent optimization of the initial stratum variation curve, reduces computational complexity during optimization, and improves the accuracy of the stratum reconstruction method for the landslide dam. Furthermore, by determining the stratum variation range of the landslide dam, the problem of multiple solutions often occurring when inverting the stratum structure based on a real dam deformation dataset is resolved. The stratum variation curve within the stratum variation range not only fits the real dam deformation dataset obtained through real-world surveys, but also demonstrates sensitivity to the response of the anti-seepage structure, further improving the accuracy of the stratum variation curve.

[0105] In an exemplary embodiment, Figure 3 As shown, the dam includes an anti-seepage structure, and the specific processing process of step 202 includes steps 302 to 308. Among them:

[0106] Step 302: Divide the landslide dam model into an initial impact area and an initial non-impact area.

[0107] During implementation, the computer device divides the landslide dam model into an initial impact area and an initial non-impact area according to the division instruction.

[0108] Specifically, the user, based on their engineering experience, operates a computer device and sends a partitioning instruction to the computer device. The partitioning instruction includes information about each control point and its location. The computer device receives the partitioning instruction and, based on the information about each control point and its location, divides the entire landslide dam model into an initial impact area and an initial non-impact area.

[0109] Step 304 : determining the first stress of the anti-seepage structure according to the material parameter set, the initial affected area, the initial non-affected area, and the finite element algorithm.

[0110] In implementation, the finite element algorithm is implemented by finite element software. The material parameter set includes a first material parameter set for the impact area and a second material parameter set for the non-impact area. A stress algorithm is provided in the finite element software. The computer device divides the initial impact area into each initial impact area grid according to the first division size corresponding to the impact area. The computer device divides the initial non-impact area into each initial non-impact area grid according to the second division size corresponding to the non-impact area. The computer device processes the data of each initial impact area grid, each initial non-impact area grid, the first material parameter set and the second material parameter set through the stress algorithm to obtain the first stress of the anti-seepage structure. Since the first division size is smaller than the second division size, the initial impact area grid is smaller than the initial non-impact area grid. That is, the density of the computer device dividing the initial impact area is higher, and the density of dividing the initial non-impact area is lower.

[0111] Step 306 , adjusting the initial impact area and the initial non-impact area, and determining the second stress of the anti-seepage structure according to the material parameter set, the adjusted initial impact area, the adjusted initial non-impact area, and the finite element algorithm.

[0112] During implementation, a stress algorithm is provided in the finite element software. The computer device reduces the initial impact area by adjusting the position of each control point to obtain an adjusted initial impact area. Then, the computer device determines the area outside the adjusted initial impact area as the initial non-impact area in the landslide dam model. The computer device divides the adjusted initial impact area into each adjusted initial impact area grid according to the first division size corresponding to the impact area. The computer device divides the adjusted initial non-impact area into each adjusted initial non-impact area grid according to the second division size corresponding to the non-impact area. The computer device processes the data of each adjusted initial impact area grid, each adjusted initial non-impact area grid, the first material parameter set and the second material parameter set through the stress algorithm to obtain the second stress of the anti-seepage structure.

[0113] Optionally, the number of control points is generally set to three. Each time the initial impact area is adjusted, the position of at least one control point is adjusted according to the adjustment requirements. The embodiment of the present application does not limit the number of control points and the number of adjusted control points.

[0114] Step 308 : determining an affected area in the landslide dam model based on the first stress, the second stress, and a preset first difference threshold.

[0115] During implementation, the computer device determines a first difference between the first stress and the second stress, and determines whether the first difference threshold satisfies a preset first difference threshold. If the first difference satisfies the first difference threshold, the computer device determines the adjusted initial impact region as the impact region in the landslide dam model. If the first difference does not satisfies the first difference threshold, the computer device updates the first stress based on the second stress and executes step 306 above until the first difference satisfies the first difference threshold, at which point the computer device determines the impact region in the landslide dam model.

[0116] In an exemplary embodiment, Figure 4 As shown, the specific processing process of step 308 includes steps 402 to 408. Among them:

[0117] Step 402 : performing difference processing on the second stress and the first stress to obtain a first difference, and determining whether the first difference exceeds a preset first difference threshold.

[0118] During implementation, a first difference threshold is pre-set in the computer device. The computer device performs difference processing on the second stress and the first stress to obtain a first difference. This first difference represents the change in the stress distribution of the anti-seepage structure. The computer device determines whether the first difference exceeds the first difference threshold, which also determines whether the change in the stress distribution of the anti-seepage structure exceeds the set first difference threshold. If the first difference exceeds the first difference threshold, the computer device executes step 404 below. If the first difference does not exceed the first difference threshold, the computer device executes step 406 below.

[0119] Optionally, the first difference threshold is determined based on attribute information of the anti-seepage structure and the landslide dam. This embodiment of the application does not limit the first difference threshold.

[0120] Step 404: If the first difference exceeds the first difference threshold, the adjusted initial impact area is determined as the impact area.

[0121] In an implementation, if the first difference exceeds a first difference threshold, the computer device determines the adjusted initial impact area as the impact area, and determines the adjusted initial non-impact area as the non-impact area.

[0122] In an exemplary embodiment, Figure 5 FIG. 1 is a schematic diagram of the stratum boundary influence level zoning in one embodiment. Figure 5 As shown in the figure, the left figure shows a landslide dam model with an anti-seepage wall as the anti-seepage structure. The right figure shows a landslide dam model with the impact area and non-impact area divided.

[0123] Step 406: If the first difference does not exceed the first difference threshold, update the first stress according to the second stress.

[0124] In implementation, if the first difference does not exceed the first difference threshold, the computer device determines that the change in stress distribution of the anti-seepage structure does not meet the first difference threshold, and the computer device determines the second stress as the first stress.

[0125] Step 408 : According to the adjusted initial impact area and the adjusted initial non-impact area, the step of adjusting the initial impact area and the initial non-impact area is performed until the first difference exceeds a first difference threshold, and the adjusted initial impact area is determined as the impact area.

[0126] In implementation, the computer device performs the steps of adjusting the initial impact area and the initial non-impact area based on the adjusted initial impact area and the adjusted initial non-impact area, that is, the computer device performs the steps of adjusting the adjusted initial impact area and the adjusted initial non-impact area until the first difference exceeds the first difference threshold. The computer device determines the adjusted initial impact area corresponding to the first difference of the current round as the impact area, and determines the adjusted initial non-impact area corresponding to the first difference of the current round as the non-impact area.

[0127] In an exemplary embodiment, to improve optimization efficiency and provide a basis for subsequent survey and monitoring arrangements, this application introduces a stratum impact level zoning method based on the response of the anti-seepage structure. Through this method and two material parameters with large differences in properties, the dam area with a greater impact on the stress response of the anti-seepage structure is identified as the affected area. Specifically, Figure 6 FIG. 1 is a flow chart for determining the impact zones in a landslide dam model in an exemplary embodiment. Figure 6 As shown in Figure 2, the process of determining the impact zones in the landslide dam model includes:

[0128] Step 601: Generate an impact zone boundary line, and divide the initial impact zone and initial non-impact zone in the landslide dam model based on engineering experience and the impact zone boundary line; the initial impact zone corresponds to the first material parameter set; the initial non-impact zone corresponds to the second material parameter set;

[0129] Step 602: determining a first stress of the anti-seepage structure using the first material parameter set, the second material parameter set, the initial impact region, the initial non-impact region, and a finite element algorithm; the stress of the anti-seepage structure represents the stress distribution of the anti-seepage structure;

[0130] Step 603: The initial impact area is reduced by adjusting the impact area boundary to obtain an adjusted initial impact area, and an adjusted initial non-impact area is determined based on the adjusted initial impact area.

[0131] Step 604 , determining a second stress of the anti-seepage structure using the first material parameter set, the second material parameter set, the adjusted initial impact region, the adjusted initial non-impact region, and a finite element algorithm;

[0132] Step 605, determining a first difference between the first stress and the second stress;

[0133] Step 606, determining whether the first difference exceeds a first difference threshold; if the first difference exceeds the first difference threshold, executing step 608; if the first difference does not exceed the first difference threshold, executing step 607;

[0134] Step 607, determining the second stress as the first stress, and executing the above step 603;

[0135] Step 608: Determine the adjusted initial impact area as the impact area, and determine the adjusted initial non-impact area as the non-impact area.

[0136] In an exemplary embodiment, Figure 7 As shown, the specific processing process of step 204 includes steps 702 to 706. Among them:

[0137] Step 702: determine the initial formation boundary variation range in the affected area, and determine the formation demarcation point in the initial formation boundary variation range based on the drilling data set.

[0138] The drilling data set includes the location of each drilling hole and the geological layer position corresponding to each drilling hole location.

[0139] In implementation, the computer device determines an initial stratum boundary variation range in the affected area and determines stratum demarcation points in the initial stratum boundary variation range according to the positions of the geological layers.

[0140] Specifically, a mapping relationship between geodetic coordinates and model coordinates of the landslide dam model is set in the computer device. The user operates the computer device based on engineering experience and sends a range determination instruction to the computer device. The range determination instruction includes the coordinates of each range. The computer device receives the range determination instruction and preliminarily determines the initial stratigraphic boundary change range in the affected area based on the coordinates of each range in the range determination instruction. The initial stratigraphic boundary change range is the range or interval of the preliminarily determined possible stratigraphic changes. The computer device determines the stratigraphic boundary point within the initial stratigraphic boundary change range based on the geological stratification positions corresponding to each drilling position and the mapping relationship.

[0141] Step 704: Generate an initial stratum boundary curve based on the stratum boundary point and the initial stratum boundary variation range.

[0142] The stratigraphic boundary points include a first stratigraphic boundary point and a second stratigraphic boundary point.

[0143] During implementation, a computer device acquires control parameters for each stratum boundary curve. Then, using a finite element algorithm and the control parameters for each stratum boundary curve, the computer device generates an initial stratum boundary curve within the initial stratum boundary variation range, starting from the first stratum boundary point and ending at the second stratum boundary point. This initial stratum boundary curve is a B-spline curve and is used to simulate the stratum boundary of the landslide dam.

[0144] In an exemplary embodiment, since the initial stratigraphic boundary variation range has an upper limit and a lower limit in the y-direction, two initial stratigraphic boundary curves need to be generated to construct the stratigraphic boundary variation range using the initial stratigraphic boundary curves. Therefore, the computer device obtains stratigraphic boundary curve control parameters. The stratigraphic boundary curve control parameters include first boundary control curve parameters and second boundary control curve parameters. The computer device uses a finite element algorithm and the control parameters of each first stratigraphic boundary curve to generate a first initial stratigraphic boundary curve within the initial stratigraphic boundary variation range, starting from the first stratigraphic boundary point and ending at the second stratigraphic boundary point. The first initial stratigraphic boundary curve coincides with the upper limit of the y-direction variation of the initial stratigraphic boundary variation range. The computer device uses a finite element algorithm and the control parameters of each second stratigraphic boundary curve to generate a second initial stratigraphic boundary curve within the initial stratigraphic boundary variation range, starting from the first stratigraphic boundary point and ending at the second stratigraphic boundary point. The second initial stratigraphic boundary curve coincides with the lower limit of the y-direction variation of the initial stratigraphic boundary variation range.

[0145] Optionally, the number of stratum boundary curve control parameters is determined according to the simulation requirements. If the simulation requirements require accurate simulation of the stratum boundary, the number of stratum boundary curve control parameters is relatively large. The embodiment of the present application does not limit the number of stratum boundary curve control parameters.

[0146] Step 706 : determining a stratum boundary variation range within the initial stratum boundary variation range according to the initial stratum boundary curve, the affected area, the material parameter set, the preset second difference threshold, and the finite element algorithm.

[0147] During implementation, the computer device determines the third stress of the anti-seepage structure based on the initial stratum boundary curve, the affected area, the material parameter set, and the finite element algorithm. The computer device then adjusts the initial stratum boundary curve toward the stratum boundary point to obtain an adjusted initial stratum boundary curve. The computer device then determines the fourth stress of the anti-seepage structure based on the adjusted initial stratum boundary curve, the affected area, the material parameter set, and the finite element algorithm. The computer device then performs a difference operation on the third and fourth stresses to obtain a second difference, and determines whether the second difference exceeds a second difference threshold. If the second difference exceeds the second difference threshold, the computer device determines the stratum boundary variation range within the initial stratum boundary variation range based on the adjusted initial stratum boundary curve.

[0148] In an exemplary embodiment, since the landslide dam model includes an impact area and an unaffected area, the computer device needs to determine the stratigraphic boundaries of the landslide dam in both the impact area and the unaffected area. The computer device determines the stratigraphic boundaries in the unaffected area based on the borehole dataset and a linear interpolation algorithm. The computer device then determines the initial stratigraphic boundary variation range in the impact area and, based on the stress variation of the anti-seepage structure, determines the stratigraphic boundary variation range within the initial stratigraphic boundary variation range. Specifically, the computer device determines the third stress of the anti-seepage structure based on the initial stratigraphic boundary curve, the impact area, the material parameter set, and a finite element method. The computer device then adjusts the initial stratigraphic boundary curve toward the stratigraphic boundary point to obtain an adjusted initial stratigraphic boundary curve. The computer device then determines the fourth stress of the anti-seepage structure based on the adjusted initial stratigraphic boundary curve, the impact area, the material parameter set, and the finite element method. The computer device then performs a difference calculation on the third and fourth stresses to obtain a second difference value, and determines whether the second difference value exceeds a second difference threshold. If the second difference exceeds the second difference threshold, the computer device determines the formation boundary variation range within the initial formation boundary variation range according to the adjusted initial formation boundary curve.

[0149] Optionally, the computer device may process each initial stratigraphic boundary curve one by one to obtain the stratigraphic boundary curve, or may process each initial stratigraphic boundary curve simultaneously to obtain the stratigraphic boundary curve.

[0150] In this embodiment, the range of stratum boundary changes in the affected area is determined through material data sets, drilling data and finite element algorithms, which provides constraints for the subsequent optimization of the initial stratum change curve, reduces the computational complexity during optimization, reduces the execution time of the post-construction reconstruction method of the landslide dam based on physical mechanisms and survey and monitoring data, and improves the efficiency of the post-construction reconstruction method of the landslide dam based on physical mechanisms and survey and monitoring data.

[0151] In an exemplary embodiment, Figure 8 As shown, the specific processing process of step 706 includes steps 802 to 806. Among them:

[0152] Step 802: Determine the third stress of the anti-seepage structure in the landslide dam based on the initial stratum boundary curve, the affected area, the material parameter set, and the finite element algorithm.

[0153] The material parameter set includes but is not limited to a first material parameter set and a second material parameter set.

[0154] During implementation, the computer system identifies the areas outside the impact zone as non-impact zones within the landslide dam model, based on the impact zone. The computer then imports the material parameter set, impact zone, non-impact zone, and initial stratigraphic boundary curve into the finite element software corresponding to the finite element algorithm. The finite element software then processes the material parameter set, impact zone, non-impact zone, and initial stratigraphic boundary curve to solve for the stress associated with the anti-seepage structure, resulting in the third stress of the anti-seepage structure. This third stress represents the stress distribution within the anti-seepage structure within the landslide dam, assuming the geological stratification of the landslide dam corresponds to the initial stratigraphic boundary curve.

[0155] In an exemplary embodiment, the initial stratum boundary curve includes a first initial stratum boundary curve and a second initial stratum boundary curve. Based on the impact area, the computer device determines the area other than the impact area as a non-impact area in the landslide dam model. Then, the computer device imports the material parameter set, the impact area, the non-impact area, and the first initial stratum boundary curve into finite element software corresponding to the finite element algorithm. The finite element software performs data processing on the material parameter set, the impact area, the non-impact area, and the first initial stratum boundary curve to solve the stress related to the anti-seepage structure and obtain the third stress corresponding to the first initial stratum boundary curve. The computer device imports the material parameter set, the impact area, the non-impact area, and the second initial stratum boundary curve into finite element software corresponding to the finite element algorithm. The finite element software performs data processing on the material parameter set, the impact area, the non-impact area, and the second initial stratum boundary curve to solve the stress related to the anti-seepage structure and obtain the third stress corresponding to the second initial stratum boundary curve.

[0156] Step 804: Adjust the initial stratum boundary curve toward the stratum boundary point, and determine the fourth stress of the anti-seepage structure based on the adjusted initial stratum boundary curve, the affected area, the material parameter set, and the finite element algorithm.

[0157] During implementation, the computer adjusts the initial stratigraphic boundary curve toward the stratigraphic boundary point, bringing the adjusted initial stratigraphic boundary curve closer overall to the stratigraphic boundary point. The computer imports the material parameter set, the impact area, the non-impact area, and the adjusted initial stratigraphic boundary curve into the finite element software corresponding to the finite element algorithm. The finite element software processes the material parameter set, the impact area, the non-impact area, and the adjusted initial stratigraphic boundary curve to solve for the stress associated with the anti-seepage structure, thereby obtaining the fourth stress of the anti-seepage structure. This fourth stress represents the stress distribution within the anti-seepage structure of the landslide dam, given that the geological stratification of the landslide dam corresponds to the adjusted initial stratigraphic boundary curve.

[0158] In an exemplary embodiment, the initial stratigraphic boundary curve includes a first initial stratigraphic boundary curve and a second initial stratigraphic boundary curve. The first initial stratigraphic boundary curve coincides with the upper limit of the initial stratigraphic boundary variation range in the y direction. The second initial stratigraphic boundary curve coincides with the lower limit of the initial stratigraphic boundary variation range in the y direction. The computer device adjusts the first initial stratigraphic boundary curve toward the stratigraphic boundary point and adjusts the second initial stratigraphic boundary curve toward the stratigraphic boundary point. At this time, the interval formed by the adjusted first initial stratigraphic boundary curve and the adjusted second initial stratigraphic boundary curve is smaller than the interval formed by the first initial stratigraphic boundary curve and the second initial stratigraphic boundary curve. Therefore, by continuously adjusting the initial stratigraphic boundary curve toward the stratigraphic boundary point, the initial stratigraphic boundary variation range can be continuously narrowed. The computer device imports the first material parameter set, the second material parameter set, the affected area, the non-affected area, and the adjusted first initial stratum boundary curve into finite element software corresponding to the finite element algorithm. The finite element software processes the first material parameter set, the second material parameter set, the affected area, the non-affected area, and the adjusted first initial stratum boundary curve to solve for stress associated with the anti-seepage structure, thereby obtaining a fourth stress associated with the adjusted first initial stratum boundary curve. The computer device imports the material parameter set, the affected area, the non-affected area, and the adjusted second initial stratum boundary curve into finite element software corresponding to the finite element algorithm. The finite element software processes the material parameter set, the affected area, the non-affected area, and the adjusted second initial stratum boundary curve to solve for stress associated with the anti-seepage structure, thereby obtaining a fourth stress associated with the adjusted second initial stratum boundary curve.

[0159] Step 806 : determining a formation boundary variation range within the initial formation boundary variation range based on the third stress, the fourth stress, the adjusted initial formation boundary curve, and a preset second difference threshold.

[0160] During implementation, the computer device performs difference processing on the third stress and the fourth stress to obtain a second difference. The computer device then determines whether the second difference exceeds a preset second difference threshold. If the second difference exceeds the second difference threshold, the computer device determines the formation boundary variation range within the initial formation boundary variation range based on the adjusted initial formation boundary curve. If the second difference does not exceed the second difference threshold, the computer device updates the third stress based on the fourth stress and continues executing step 804 above until the second difference exceeds the second difference threshold, at which point the computer device determines the formation boundary variation range within the initial formation boundary variation range based on the adjusted initial formation boundary curve.

[0161] In an exemplary embodiment, Figure 9 As shown, the specific processing process of step 806 includes steps 902 to 908. Among them:

[0162] Step 902 : performing difference processing on the third stress and the fourth stress to obtain a second difference, and determining whether the second difference exceeds a preset second difference threshold.

[0163] In an embodiment, a computer device performs difference processing on the third stress corresponding to the first initial stratigraphic boundary curve and the fourth stress corresponding to the adjusted first initial stratigraphic boundary curve to obtain a second difference corresponding to the first initial stratigraphic boundary curve. The computer device also performs difference processing on the third stress corresponding to the second initial stratigraphic boundary curve and the fourth stress corresponding to the adjusted second initial stratigraphic boundary curve to obtain a second difference corresponding to the second initial stratigraphic boundary curve. The computer device determines whether the second difference corresponding to the first initial stratigraphic boundary curve exceeds a preset second difference threshold, and also determines whether the second difference corresponding to the second initial stratigraphic boundary curve exceeds the preset second difference threshold.

[0164] Optionally, the second difference threshold is determined based on attribute information of the landslide dam. This embodiment of the present application does not limit the second difference threshold.

[0165] Step 904 : If the second difference exceeds the second difference threshold, determine the formation boundary variation range within the initial formation boundary variation range based on the adjusted initial formation boundary curve.

[0166] In implementation, if the second difference exceeds the second difference threshold, the computer device determines the adjusted initial stratum boundary curve as the stratum boundary curve. Then, the computer device determines the area enclosed by the stratum boundary curve within the initial stratum boundary variation range as the stratum boundary variation range.

[0167] Specifically, the initial stratigraphic boundary curve includes a first initial stratigraphic boundary curve and a second initial stratigraphic boundary curve. If the second difference exceeds a second difference threshold, the computer device determines the adjusted first initial stratigraphic boundary curve as the first stratigraphic boundary curve, and determines the adjusted second initial stratigraphic boundary curve as the second stratigraphic boundary curve. Because the starting point and end point of the first stratigraphic boundary curve and the second stratigraphic boundary curve are both stratigraphic boundary points, and the first stratigraphic boundary curve and the second stratigraphic boundary curve are both within the initial stratigraphic boundary variation range, the first stratigraphic boundary curve and the second stratigraphic boundary curve can enclose an area within the initial stratigraphic boundary variation range. The computer device determines the area enclosed by the first stratigraphic boundary curve and the second stratigraphic boundary curve as the stratigraphic boundary variation range.

[0168] Step 906 : If the second difference does not exceed the second difference threshold, update the third stress based on the fourth stress.

[0169] During implementation, the computer device determines the fourth stress corresponding to the adjusted first initial formation boundary curve as the third stress corresponding to the first initial formation boundary curve, thereby obtaining the updated third stress corresponding to the first initial formation boundary curve. The computer device determines the fourth stress corresponding to the adjusted second initial formation boundary curve as the third stress corresponding to the second initial formation boundary curve, thereby obtaining the updated third stress corresponding to the second initial formation boundary curve.

[0170] Step 908, according to the adjusted initial stratigraphic boundary curve, execute the step of adjusting the initial stratigraphic boundary curve toward the stratigraphic boundary point until the second difference exceeds the second difference threshold, and based on the adjusted initial stratigraphic boundary curve, determine the stratigraphic boundary change range in the initial stratigraphic boundary change range.

[0171] During implementation, the computer device continues to execute step 804 based on the adjusted initial stratigraphic boundary curve until the second difference exceeds the second difference threshold, at which point the computer device determines the adjusted initial stratigraphic boundary curve as the stratigraphic boundary curve. The computer device then determines the area enclosed by the stratigraphic boundary curve within the initial stratigraphic boundary variation range as the stratigraphic boundary variation range. The specific processing of step 804 has been described in detail in the above embodiment and will not be repeated here in this embodiment.

[0172] In an exemplary embodiment, Figure 10 FIG. 1 is a flow chart of determining the range of formation boundary changes in an exemplary embodiment. Figure 10As shown, the computer device preliminarily determines the range of possible stratigraphic boundary variations within the affected area, obtaining an initial stratigraphic boundary variation range. The computer device then uses the borehole dataset to identify stratigraphic demarcation points within the initial stratigraphic boundary variation range. The computer device uses the stratigraphic demarcation points as the two ends of an initial stratigraphic boundary curve and uses a B-spline curve to simulate the stratigraphic boundary within the initial stratigraphic boundary variation range, obtaining an initial stratigraphic boundary curve. The computer device then determines the third stress of the anti-seepage structure in the landslide dam based on the initial stratigraphic boundary curve, the affected area, the material parameter set, and a finite element method. The computer device adjusts the initial stratigraphic boundary curve in the y-direction toward the stratigraphic demarcation point, such that the adjusted initial stratigraphic boundary curve continuously approaches the stratigraphic demarcation point. The computer device then determines the fourth stress of the anti-seepage structure based on the adjusted initial stratigraphic boundary curve, the affected area, the material parameter set, and the finite element method. The computer device then determines a second difference between the third and fourth stresses. This second difference represents the impact of the adjusted initial stratigraphic boundary curve on the stress response of the anti-seepage structure. The computer device determines whether the second difference exceeds a preset second difference threshold, thereby determining whether the adjusted initial formation boundary curve has a significant impact on the stress response of the anti-seepage structure. If the second difference exceeds the second difference threshold, the computer device determines the adjusted initial formation boundary curve as the formation boundary curve. This formation boundary curve is the formation boundary curve that first significantly affects the stress response of the anti-seepage structure. The computer device then determines the interval enclosed by the formation boundary curve within the initial formation boundary variation range as the formation boundary variation range (formation boundary envelope). If the second difference does not exceed the second difference threshold, the computer device determines the fourth stress as the third stress and executes the step of adjusting the initial formation boundary curve in the y-direction toward the formation boundary point until the second difference exceeds the second difference threshold. The computer device then determines the interval enclosed by the formation boundary curve within the initial formation boundary variation range as the formation boundary variation range, thereby providing a constraint condition for subsequent optimization.

[0173] In an exemplary embodiment, the initial formation change curve includes various control parameters, such as Figure 11 As shown, the specific processing process of step 106 includes steps 1102 to 1108. Among them:

[0174] Step 1102 , optimizing the initial formation change curve according to various control parameters, finite element algorithm and optimization algorithm to obtain an optimized formation change curve, and performing dam body prediction according to the optimized formation change curve to obtain a predicted dam body deformation data set.

[0175] Among them, the real dam deformation dataset contains the real displacement data of the landslide dam during the historical time period.

[0176] During implementation, the computer obtains each initial control parameter and the corresponding disturbance value, and perturbs each control parameter based on the disturbance value to obtain the initial ground change curve after the disturbance. The computer then optimizes the control parameters based on the perturbed initial ground change curve and the actual predicted dam deformation data set, obtaining optimized control parameters and generating an optimized initial ground change curve based on the optimized control parameters. Using finite element software corresponding to the finite element algorithm, the computer processes the optimized initial ground change curve, the material parameter set, the borehole data set, the affected area, and the non-affected area to predict the displacement of the landslide dam and obtain predicted dam deformation data.

[0177] Step 1104 : Calculate the loss between the predicted dam deformation dataset and the actual dam deformation dataset, and determine whether the loss exceeds a preset loss threshold.

[0178] The loss condition is that the loss between the predicted dam deformation dataset and the actual dam deformation dataset does not exceed the loss threshold.

[0179] During implementation, a loss threshold is pre-set in the computer device. Based on a pre-set loss algorithm, the computer device performs data operations on the predicted dam deformation dataset and the actual dam deformation dataset to determine the loss between the predicted and actual dam deformation datasets. The computer device then determines whether the current loss exceeds the loss threshold. If the loss exceeds the loss threshold, the computer device executes step 1108. If the loss does not exceed the loss threshold, the computer device executes step 1106.

[0180] Optionally, the loss threshold is determined based on the required accuracy of the formation change curve simulation. If the required accuracy of the formation change curve simulation is high, the loss threshold is set to be small; if the required accuracy of the formation change curve simulation is low, the loss threshold is set to be large. This embodiment of the application does not limit the loss threshold.

[0181] Step 1106: If the loss does not exceed the loss threshold, it is determined that the loss meets the preset loss condition.

[0182] In implementation, if the loss does not exceed the loss threshold, the computer device determines that the loss meets the preset loss condition and continues to execute the above step 108.

[0183] Step 1108, if the loss exceeds the loss threshold, execute the step of optimizing the initial formation change curve according to various control parameters, finite element algorithm and optimization algorithm until the loss does not exceed the loss threshold, and determine that the loss meets the preset loss condition.

[0184] In implementation, if the loss exceeds the loss threshold, the computer device determines that the loss does not meet the preset loss condition. If the loss does not meet the loss condition, the computer device continues to execute step 1102 until the loss does not exceed the loss threshold, i.e., the loss meets the preset loss condition. The specific processing process of step 1102 has been described in detail in the above embodiment and will not be repeated in this embodiment.

[0185] In an exemplary embodiment, the initial formation change curve includes various control parameters, such as Figure 12 As shown, the specific process of optimizing the initial formation change curve according to various control parameters, finite element algorithm and optimization algorithm in step 1102 to obtain the optimized formation change curve includes steps 1202 to 1208.

[0186] Step 1202: Determine the disturbance value of each control parameter, and perturb each control parameter according to the disturbance value of each control parameter to generate an initial formation change curve after disturbance.

[0187] During implementation, the computer device determines the disturbance value of each control parameter. For example, the user inputs the disturbance value of each control parameter into the computer device by operating the computer device. The disturbance value includes a positive disturbance value and a negative disturbance value. The computer device updates the control parameter according to the disturbance value of each control parameter to obtain the control parameter after disturbance. The control parameter after disturbance includes each control parameter after positive disturbance and each control parameter after negative disturbance. Then, the computer device uses a B-spline curve to simulate the stratum boundary of the landslide dam, and generates an initial stratum change curve after disturbance that represents the stratum boundary based on the control parameter after disturbance. The initial stratum change curve after disturbance includes an initial stratum change curve after positive disturbance and an initial stratum change curve after negative disturbance.

[0188] Optionally, the disturbance value of each control parameter may be different or the same, and is determined according to the adjustment requirements of the initial formation change curve. The disturbance value of each control parameter in this application is not limited.

[0189] Step 1204 : Based on the finite element algorithm and the initial stratum change curve after disturbance, the dam body deformation is predicted to obtain a disturbed dam body deformation data set.

[0190] During implementation, the computer equipment imports the initial stratum change curve, material parameter set, borehole data set, affected area, and non-affected area after the positive disturbance into the finite element software corresponding to the finite element algorithm. The finite element software processes the initial stratum change curve, material parameter set, borehole data set, affected area, and non-affected area after the positive disturbance, predicts the displacement of the landslide dam over the historical time period, and obtains the positive disturbance dam body deformation data set. The computer equipment imports the initial stratum change curve, material parameter set, borehole data set, affected area, and non-affected area after the negative disturbance into the finite element software corresponding to the finite element algorithm. The finite element software processes the initial stratum change curve, material parameter set, borehole data set, affected area, and non-affected area after the negative disturbance, predicts the displacement of the landslide dam over the historical time period, and obtains the negative disturbance dam body deformation data set.

[0191] Step 1206: Calculate the disturbance error between the disturbed dam body deformation data set and the true dam body deformation data set.

[0192] During implementation, the computer equipment processes the disturbed dam body deformation data set and the real dam body deformation data set according to a preset loss algorithm to obtain a disturbance error between the disturbed dam body deformation data set and the real dam body deformation data set.

[0193] Specifically, the disturbed dam deformation dataset includes a positive disturbance dam deformation dataset and a negative disturbance dam deformation dataset. The computer device processes the positive disturbance dam deformation dataset and the true dam deformation dataset according to a preset loss algorithm to obtain a positive disturbance error between the positive disturbance dam deformation dataset and the true dam deformation dataset. The computer device also processes the negative disturbance dam deformation dataset and the true dam deformation dataset according to a preset loss algorithm to obtain a negative disturbance error between the negative disturbance dam deformation dataset and the true dam deformation dataset.

[0194] Step 1208 : Calculate the disturbance gradient according to the disturbance error, and optimize the control parameters according to the disturbance gradient and the optimization algorithm to generate an optimized initial formation change curve.

[0195] In practice, the optimization algorithm is a gradient descent method. A computer calculates the perturbation gradient of each control parameter using the central difference method and perturbation error. The computer updates the control parameter based on the perturbation gradient and the gradient descent method, obtaining an updated control parameter. The computer then determines whether the updated control parameter is within the formation boundary variation range, obtaining a determination result. The computer then generates an optimized formation variation curve based on the determination results, the perturbation value of each control parameter, and the updated control parameter.

[0196] Specifically, the computer device calculates the approximate partial derivative of each control parameter based on the central difference method, and calculates the perturbation gradient of the control parameter based on the approximate partial derivative. Then, the computer device updates the control parameter based on the gradient descent method and the perturbation gradient to obtain an updated control parameter. The computer device determines whether the updated control parameter is within the formation boundary variation range and obtains a judgment result. If the judgment result indicates that the updated control parameter is within the formation boundary variation range, the computer device does not perform data processing on the updated control parameter. If the judgment result indicates that the updated control parameter is not within the formation boundary variation range, the computer device adjusts the perturbation value corresponding to the control parameter and re-updates the control parameter based on the adjusted perturbation value so that the updated control parameter is within the formation boundary variation range. Then, the computer device generates an optimized formation variation curve within the formation boundary variation range based on each updated control parameter.

[0197] In an exemplary embodiment, Figure 13 FIG. 1 is a flow chart of determining a formation change curve according to an exemplary embodiment. Figure 13 As shown in FIG, the specific processing process of determining the formation change curve includes:

[0198] Step 1301, obtain each control parameter, and generate an initial formation change curve according to each control parameter; the control parameters are used express;

[0199] Step 1302, obtaining a disturbance value of each control parameter;

[0200] Step 1303: Disturb the control parameter according to the disturbance value of the control parameter to obtain the disturbed control parameter; wherein the disturbed control parameter is ; For the control parameters, is the disturbance value;

[0201] Step 1304: generating an optimized disturbed formation change curve based on the control parameters after each disturbance;

[0202] Step 1305 , based on the finite element algorithm and the initial stratum change curve after disturbance, dam body deformation prediction is performed to obtain a disturbed dam body deformation data set;

[0203] Step 1306: Calculate the disturbance error between the disturbance dam deformation data set and the real dam deformation data set. and ;

[0204] Step 1307: Calculate the disturbance gradient of each control parameter based on the disturbance error. ;

[0205] Step 1308: Update the control parameters according to the disturbance gradient to obtain updated control parameters. ;in, After the update control parameters, For the control parameters, is the perturbation gradient, is the learning rate;

[0206] Step 1309: Determine whether the updated control parameter is within the formation boundary variation range; if the updated control parameter is within the formation boundary variation range, execute step 1310; if the updated control parameter is not within the formation boundary variation range, execute step 1302;

[0207] Step 1310: generating an optimized formation change curve based on the updated control parameters;

[0208] Step 1311 , based on the finite element algorithm and the optimized initial stratum change curve, dam deformation prediction is performed to obtain a predicted dam deformation data set;

[0209] Step 1312, calculating the loss between the predicted dam deformation dataset and the actual dam deformation dataset;

[0210] Step 1313, determining whether the loss exceeds a preset loss threshold. If the loss exceeds the loss threshold, execute step 1302. If the error does not exceed the preset error threshold, execute step 1314.

[0211] Step 1314: determine the optimized initial formation change curve as the formation change curve.

[0212] In an exemplary embodiment, Figure 14 FIG. 1 is a flow chart of a post-construction reconstruction method of a landslide dam based on physical mechanisms and survey monitoring data in an exemplary embodiment. Figure 14 As shown in FIG, the specific processing process of the post-construction reconstruction method of the landslide dam based on physical mechanism and survey monitoring data includes:

[0213] Step 1401: Establish a landslide dam model, and determine an impact area in the landslide dam model based on the response of the anti-seepage structure in the landslide dam;

[0214] Step 1402 , determining a stratum boundary variation range in an impact area based on the borehole information in the borehole dataset;

[0215] Step 1403: Generate a stratum change curve in the affected area based on the real dam deformation data set.

[0216] The specific processing of step 1401 includes establishing a landslide dam model and performing a sensitivity analysis on the response of the anti-seepage structure to determine the affected and non-affected areas of the landslide dam model. The specific processing of step 1402 includes determining the initial stratigraphic boundary variation range (stratigraphic boundary variable range) in the affected area based on the borehole dataset, and continuously changing the boundaries of the stratigraphic boundary variable range through anti-seepage structure response analysis to obtain the stratigraphic boundary variation range (stratigraphic boundary envelope range). The specific processing of step 1403 includes continuously changing the stratigraphic boundaries within the stratigraphic boundary envelope range based on the initial stratigraphic variation curve, and performing finite element calculations of the landslide dam to ensure that the predicted dam deformation dataset matches the actual dam deformation dataset (monitoring data) to obtain the stratigraphic variation curve (optimal stratigraphic boundary).

[0217] It should be understood that, although the steps in the flowcharts involved in the embodiments described above are displayed in sequence as indicated by the arrows, these steps are not necessarily executed in the order indicated by the arrows.

[0218] In an exemplary embodiment, Figure 15 As shown, a post-construction reconstruction device 1500 for a landslide dam based on physical mechanisms and survey and monitoring data is provided, comprising: an acquisition module 1501, a first determination module 1502, a prediction module 1503, and a second determination module 1504, wherein:

[0219] The acquisition module 1501 is used to establish a landslide dam model of the landslide dam and obtain a material parameter set, a drilling data set and a real dam body deformation data set of the landslide dam.

[0220] The first determination module 1502 is configured to determine the stratum boundary variation range in the landslide dam model based on the material parameter set, the drilling data set, and the finite element algorithm, and to determine an initial stratum variation curve within the stratum boundary variation range.

[0221] The prediction module 1503 is used to predict the dam deformation based on the initial formation change curve and the finite element algorithm, obtain the predicted dam deformation data set, and calculate the loss between the predicted dam deformation data set and the actual dam deformation data set until the loss meets the preset loss condition.

[0222] The second determining module 1504 is configured to determine the initial stratum change curve corresponding to the dam deformation prediction data set of the current round as the stratum change curve.

[0223] In an exemplary embodiment, the first determining module 1502 includes:

[0224] The first determination submodule is used to determine the affected area in the landslide dam model based on the material parameter set, the finite element algorithm and a preset first difference threshold; the first difference threshold is the stress change threshold of the anti-seepage structure in the landslide dam.

[0225] The second determination submodule is used to determine the range of changes in the stratum boundary in the affected area based on the drilling data set, the material parameter set, the preset second difference threshold and the finite element algorithm; the second difference threshold is the stress change threshold of the anti-seepage structure in the landslide dam.

[0226] In an exemplary embodiment, the barrier dam includes an anti-seepage structure, and the first determining submodule includes:

[0227] The first division submodule is used to divide the landslide dam model into an initial impact area and an initial non-impact area.

[0228] The third determination submodule is used to determine the first stress of the anti-seepage structure according to the material parameter set, the initial influence area, the initial non-influence area and the finite element algorithm.

[0229] The first adjustment submodule is used to adjust the initial influence area and the initial non-influence area, and determine the second stress of the anti-seepage structure according to the material parameter set, the adjusted initial influence area, the adjusted initial non-influence area and the finite element algorithm.

[0230] The fourth determining submodule is configured to determine an affected area in the landslide dam model based on the first stress, the second stress, and a preset first difference threshold.

[0231] In an exemplary embodiment, the fourth determination submodule is specifically configured to: a first processing submodule for performing difference processing on the second stress and the first stress to obtain a first difference, and determining whether the first difference exceeds a preset first difference threshold; a fifth determination submodule for determining the adjusted initial impact region as the impact region if the first difference exceeds the first difference threshold; a first updating submodule for updating the first stress based on the second stress if the first difference does not exceed the first difference threshold; and a first execution submodule for executing the steps of adjusting the initial impact region and the initial non-impact region based on the adjusted initial impact region and the adjusted initial non-impact region until the first difference exceeds the first difference threshold, thereby determining the adjusted initial impact region as the impact region.

[0232] In an exemplary embodiment, the second determining submodule includes:

[0233] The sixth determination submodule is configured to determine an initial stratum boundary variation range in the affected area, and determine a stratum demarcation point in the initial stratum boundary variation range according to the drilling data set.

[0234] The first generating submodule is used to generate an initial stratum boundary curve based on the stratum boundary point and the initial stratum boundary variation range.

[0235] The seventh determination submodule is used to determine the formation boundary variation range in the initial formation boundary variation range according to the initial formation boundary curve, the influence area, the material parameter set, the preset second difference threshold and the finite element algorithm.

[0236] In an exemplary embodiment, the seventh determining submodule includes:

[0237] The eighth determination submodule is used to determine the third stress of the anti-seepage structure in the landslide dam based on the initial stratum boundary curve, the affected area, the material parameter set and the finite element algorithm.

[0238] The second adjustment submodule is used to adjust the initial stratum boundary curve toward the stratum boundary point, and determine the fourth stress of the anti-seepage structure according to the adjusted initial stratum boundary curve, the affected area, the material parameter set and the finite element algorithm.

[0239] The ninth determination submodule is configured to determine the formation boundary variation range within the initial formation boundary variation range based on the third stress, the fourth stress, the adjusted initial formation boundary curve, and a preset second difference threshold.

[0240] In an exemplary embodiment, the ninth determination submodule is specifically used for: a second processing submodule, for performing difference processing on the third stress and the fourth stress to obtain a second difference, and judging whether the second difference exceeds a preset second difference threshold; a tenth determination submodule, for determining the formation boundary change range in the initial formation boundary change range based on the adjusted initial formation boundary curve if the second difference exceeds the second difference threshold; a second updating submodule, for updating the third stress based on the fourth stress if the second difference does not exceed the second difference threshold; a second execution submodule, for executing the step of adjusting the initial formation boundary curve toward the formation boundary point according to the adjusted initial formation boundary curve until the second difference exceeds the second difference threshold, and determining the formation boundary change range in the initial formation boundary change range based on the adjusted initial formation boundary curve.

[0241] In an exemplary embodiment, the initial formation change curve includes various control parameters, and the prediction module 1503 includes:

[0242] The first optimization submodule is used to optimize the initial stratum change curve according to various control parameters, finite element algorithm and optimization algorithm to obtain an optimized stratum change curve, and perform dam body prediction based on the optimized stratum change curve to obtain a predicted dam body deformation data set.

[0243] The first calculation submodule is used to calculate the loss between the predicted dam deformation data set and the actual dam deformation data set, and determine whether the loss exceeds a preset loss threshold.

[0244] The eleventh determining submodule is configured to determine whether the loss satisfies a preset loss condition if the loss does not exceed the loss threshold.

[0245] The twelfth determination submodule is used to execute the steps of optimizing the initial formation change curve according to various control parameters, finite element algorithm and optimization algorithm if the loss exceeds the loss threshold, until the loss does not exceed the loss threshold, and determine that the loss meets the preset loss condition.

[0246] In an exemplary embodiment, the first optimization submodule includes a second optimization submodule and a first prediction submodule. The first optimization submodule is specifically used in the perturbation submodule to determine the perturbation value of each control parameter and perturb each control parameter according to the perturbation value of each control parameter to generate an initial formation change curve after perturbation; the second prediction submodule is used to predict dam deformation based on the finite element algorithm and the initial formation change curve after perturbation to obtain a perturbed dam deformation data set; the second calculation submodule is used to calculate the perturbation error between the perturbed dam deformation data set and the actual dam deformation data set; and the second generation submodule is used to calculate the perturbation gradient based on the perturbation error, optimize each control parameter based on the perturbation gradient and the optimization algorithm, and generate an optimized initial formation change curve.

[0247] Each module in the aforementioned post-construction reconstruction device for landslide dams based on physical mechanisms and survey and monitoring data can be implemented in whole or in part through software, hardware, or a combination thereof. Each module can be embedded in or independent of a processor in a computer device in the form of hardware, or can be stored in a memory in the computer device in the form of software, so that the processor can call and execute the corresponding operations of each module.

[0248] In an exemplary embodiment, a computer device is provided. The computer device may be a terminal, and its internal structure diagram may be as shown in FIG. Figure 16As shown. The computer device includes a processor, memory, an input / output interface, a communication interface, a display unit, and an input device. The processor, memory, and input / output interface are connected via a system bus, and the communication interface, display unit, and input device are connected to the system bus via the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and computer program in the non-volatile storage medium. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal via wired or wireless communication, and the wireless communication can be implemented via Wi-Fi, a mobile cellular network, near field communication (NFC), or other technologies. When the computer program is executed by the processor, it implements a landslide dam prediction method.

[0249] Those skilled in the art will understand that Figure 16 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have a different component arrangement.

[0250] In one embodiment, a computer device is further provided, including a memory and a processor. The memory stores a computer program, and the processor implements the steps in the above method embodiments when executing the computer program.

[0251] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the steps in the above-mentioned method embodiments are implemented.

[0252] In one embodiment, a computer program product is provided, including a computer program, which implements the steps in the above method embodiments when executed by a processor.

[0253] The technical features of the above embodiments can be combined arbitrarily. In order to make the description concise, not all possible combinations of the technical features in the above 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 application.

[0254] The above-described embodiments merely represent several implementation methods of the present application. While the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present application. It should be noted that a person of ordinary skill in the art may make various modifications and improvements without departing from the spirit of the present application, and these modifications and improvements fall within the scope of protection of the present application. Therefore, the scope of protection of the present application shall be determined by the appended claims.

Claims

1. A post-construction reconstruction method for a landslide dam based on physical mechanisms and survey and monitoring data, characterized in that: The method comprises: Establishing a landslide dam model of the landslide dam, and obtaining a material parameter set, a drilling data set, and a real dam body deformation data set of the landslide dam; Determining a stratum boundary variation range in the landslide dam model based on the material parameter set, the drilling data set, and a finite element algorithm, and determining an initial stratum variation curve within the stratum boundary variation range; dam deformation prediction is performed based on the initial stratum change curve and the finite element algorithm to obtain a predicted dam deformation data set, and a loss between the predicted dam deformation data set and the actual dam deformation data set is calculated until the loss satisfies a preset loss condition; Determine the initial stratum change curve corresponding to the predicted dam deformation data set of the current round as the stratum change curve; The initial formation change curve includes various control parameters. The dam deformation is predicted based on the initial formation change curve and the finite element algorithm to obtain a predicted dam deformation data set, and the loss between the predicted dam deformation data set and the actual dam deformation data set is calculated until the loss meets a preset loss condition, including: Optimizing the initial stratum change curve according to the control parameters, the finite element algorithm, and the optimization algorithm to obtain an optimized stratum change curve, and performing dam body prediction according to the optimized stratum change curve to obtain a predicted dam body deformation data set; Calculating the loss between the predicted dam deformation dataset and the actual dam deformation dataset, and determining whether the loss exceeds a preset loss threshold; If the loss does not exceed the loss threshold, determining that the loss meets the preset loss condition; If the loss exceeds the loss threshold, the step of optimizing the initial formation change curve according to the control parameters, the finite element algorithm and the optimization algorithm is executed until the loss does not exceed the loss threshold, and it is determined that the loss meets the preset loss condition.

2. The method according to claim 1, characterized in that The determining of the stratum boundary variation range in the landslide dam model based on the material parameter set, the drilling data set, and the finite element algorithm includes: Determining an affected area in the landslide dam model based on the material parameter set, the finite element algorithm, and a preset first difference threshold value, wherein the first difference threshold value is a stress change threshold value of the anti-seepage structure in the landslide dam; The range of stratum boundary changes in the affected area is determined based on the drilling data set, the material parameter set, a preset second difference threshold and the finite element algorithm; the second difference threshold is a stress change threshold of the anti-seepage structure in the landslide dam.

3. The method according to claim 2, characterized in that The landslide dam includes an anti-seepage structure, and determining the affected area in the landslide dam model based on the material parameter set, the finite element algorithm, and a preset first difference threshold value includes: Dividing the landslide dam model into an initial impact area and an initial non-impact area; determining a first stress of the anti-seepage structure according to the material parameter set, the initial affected area, the initial non-affected area, and the finite element algorithm; Adjusting the initial impact area and the initial non-impact area, and determining a second stress of the anti-seepage structure according to the material parameter set, the adjusted initial impact area, the adjusted initial non-impact area, and the finite element algorithm; An affected area in the landslide dam model is determined based on the first stress, the second stress, and a preset first difference threshold.

4. The method according to claim 3, characterized in that The determining of the affected area in the landslide dam model based on the first stress, the second stress, and a preset first difference threshold includes: performing difference processing on the second stress and the first stress to obtain a first difference, and determining whether the first difference exceeds a preset first difference threshold; If the first difference exceeds the first difference threshold, determining the adjusted initial impact area as the impact area; If the first difference does not exceed the first difference threshold, updating the first stress according to the second stress; According to the adjusted initial impact area and the adjusted initial non-impact area, the step of adjusting the initial impact area and the initial non-impact area is performed until the first difference exceeds the first difference threshold, and the adjusted initial impact area is determined as the impact area.

5. The method according to claim 2, characterized in that The determining, based on the drilling data set, the material parameter set, a preset second difference threshold, and the finite element algorithm, of a stratum boundary variation range in the affected area includes: determining an initial stratigraphic boundary variation range in the affected area, and determining a stratigraphic demarcation point in the initial stratigraphic boundary variation range based on the drill hole dataset; generating an initial stratum boundary curve based on the stratum boundary point and the initial stratum boundary variation range; A stratum boundary variation range is determined within the initial stratum boundary variation range according to the initial stratum boundary curve, the affected area, the material parameter set, a preset second difference threshold, and the finite element algorithm.

6. The method according to claim 5, characterized in that The determining of the formation boundary variation range within the initial formation boundary variation range according to the initial formation boundary curve, the affected area, the material parameter set, a preset second difference threshold, and the finite element algorithm includes: determining a third stress of the anti-seepage structure in the landslide dam according to the initial stratum boundary curve, the affected area, the material parameter set, and a finite element algorithm; Adjusting the initial stratum boundary curve toward the stratum boundary point, and determining the fourth stress of the anti-seepage structure based on the adjusted initial stratum boundary curve, the affected area, the material parameter set, and a finite element algorithm; Based on the third stress, the fourth stress, the adjusted initial formation boundary curve and a preset second difference threshold, a formation boundary variation range is determined within the initial formation boundary variation range.

7. The method according to claim 6, characterized in that Determining the formation boundary variation range within the initial formation boundary variation range based on the third stress, the fourth stress, the adjusted initial formation boundary curve, and a preset second difference threshold value includes: performing difference processing on the third stress and the fourth stress to obtain a second difference, and determining whether the second difference exceeds a preset second difference threshold; If the second difference exceeds the second difference threshold, determining a formation boundary variation range within the initial formation boundary variation range based on the adjusted initial formation boundary curve; If the second difference does not exceed the second difference threshold, updating the third stress based on the fourth stress; The step of adjusting the initial stratigraphic boundary curve toward the stratigraphic boundary point is performed according to the adjusted initial stratigraphic boundary curve until the second difference exceeds the second difference threshold, and the stratigraphic boundary change range is determined within the initial stratigraphic boundary change range based on the adjusted initial stratigraphic boundary curve.

8. The method according to claim 1, characterized in that Optimizing the initial formation change curve according to the control parameters, the finite element algorithm, and the optimization algorithm to obtain the optimized formation change curve includes: Determining a disturbance value of each of the control parameters, and perturbing each of the control parameters according to the disturbance value of each of the control parameters to generate the initial formation change curve after disturbance; Based on the finite element algorithm and the initial stratum change curve after disturbance, dam body deformation prediction is performed to obtain a disturbed dam body deformation data set; Calculating a disturbance error between the disturbance dam body deformation data set and the real dam body deformation data set; A disturbance gradient is calculated according to the disturbance error, and each of the control parameters is optimized according to the disturbance gradient and an optimization algorithm to generate an optimized initial formation change curve.

9. The method according to claim 8, characterized in that The optimization algorithm is a gradient descent method, wherein the disturbance gradient is calculated according to the disturbance error, and each control parameter is optimized according to the disturbance gradient and the optimization algorithm to generate an optimized initial formation change curve, including: Calculating the disturbance gradient of each control parameter according to the central difference method and the disturbance error, and updating the control parameter according to the disturbance gradient and the gradient descent method to obtain an updated control parameter; Determining whether the updated control parameter is within the formation boundary variation range, and obtaining a determination result; Based on the judgment results, the disturbance values ​​of the control parameters and the updated control parameters, an optimized formation change curve is generated.

10. A post-construction reconstruction device for a landslide dam based on physical mechanisms and survey and monitoring data, characterized in that: The device comprises: an acquisition module, configured to establish a landslide dam model of the landslide dam and acquire a material parameter set, a drilling data set, and a real dam body deformation data set of the landslide dam; a first determining module, configured to determine a stratum boundary variation range in the landslide dam model based on the material parameter set, the drilling data set, and a finite element algorithm, and to determine an initial stratum variation curve within the stratum boundary variation range; a prediction module, configured to predict dam deformation based on the initial stratum change curve and the finite element algorithm, obtain a predicted dam deformation data set, and calculate a loss between the predicted dam deformation data set and the actual dam deformation data set until the loss satisfies a preset loss condition; A second determining module is configured to determine the initial stratum change curve corresponding to the predicted dam deformation data set of the current round as the stratum change curve; The initial formation change curve includes various control parameters. The prediction module is specifically used to optimize the initial formation change curve according to each of the control parameters, the finite element algorithm and the optimization algorithm to obtain the optimized formation change curve, and perform dam body prediction according to the optimized formation change curve to obtain a predicted dam body deformation data set; calculate the loss between the predicted dam body deformation data set and the actual dam body deformation data set, and determine whether the loss exceeds a preset loss threshold; if the loss does not exceed the loss threshold, determine that the loss meets the preset loss condition; if the loss exceeds the loss threshold, execute the step of optimizing the initial formation change curve according to each of the control parameters, the finite element algorithm and the optimization algorithm until the loss does not exceed the loss threshold, and determine that the loss meets the preset loss condition.

Images