Barrier dam post-construction reconstruction method based on physical mechanism and survey monitoring data

Abstract

The invention relates to a barrier dam post-construction reconstruction method and device based on a physical mechanism and survey monitoring data and computer equipment. The method comprises the steps that a barrier dam model of the barrier dam is established, and a material parameter set, a drilling data set and a real dam body deformation data set of the barrier dam are obtained; based on the material parameter set, the drilling data set and a finite element algorithm, determining a stratum boundary change range in the barrier dam model, and determining an initial stratum change curve in the stratum boundary change range; according to the initial stratum change curve and a finite element algorithm, dam body deformation prediction is carried out, a predicted dam body deformation data set is obtained, and the loss between the predicted dam body deformation data set and the real dam body deformation data set is calculated until the loss meets a preset loss condition; and determining the initial stratum change curve corresponding to the dam body deformation prediction data set in the current round as a stratum change curve. By adopting the method, the accuracy of the barrier dam post-construction reconstruction method based on a physical mechanism and survey monitoring data can be improved.

Description

Technical Field

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

[0002] A barrier lake is a natural dam formed by geological disasters blocking a river channel, and its stability is directly related to the safety of the downstream area. To ensure the safety of the downstream area, it is necessary to construct the stratum structure of the barrier lake and predict whether the barrier lake is in a safe state based on the stratum structure of the barrier lake.

[0003] In traditional technologies, limited borehole data of each borehole in the barrier lake are obtained. For two adjacent boreholes, the stratum structure where the boreholes are located is determined based on the borehole data of the two adjacent boreholes, and the geological statistics interpolation method is used to predict the stratum change curve in the area between the two adjacent boreholes. Then, the stratum structure of the barrier lake is constructed based on each stratum change curve and the stratum structure of each borehole.

[0004] However, in traditional technologies, due to the complex geological structure of the barrier lake, it is difficult to accurately predict the stratum structure in the area between two adjacent boreholes only relying on the geological statistics interpolation method. Therefore, the accuracy of the current barrier lake construction method is relatively low. Summary of the Invention

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

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

[0007] Establish a barrier lake model of the barrier lake, and obtain a material parameter set, a borehole data set, and a real dam body deformation data set of the barrier lake;

[0008] Based on the material parameter set, the borehole data set, and the finite element algorithm, determine the stratum boundary change range in the barrier lake model, and determine an initial stratum change curve in the stratum boundary change range;

[0009] According to the initial stratum change curve and the finite element algorithm, perform dam body deformation prediction to obtain a predicted dam body deformation data set, and calculate 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;

[0010] Determine the initial stratum change curve corresponding to the predicted dam body deformation data set in the current round as the stratum change curve.

[0011] In one embodiment, determining the range of formation boundary changes in the barrier lake dam model based on the material parameter set, the borehole data set, and the finite element algorithm includes:

[0012] Determining an influence area in the barrier lake 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 barrier lake dam.

[0013] Determining the range of formation boundary changes in the influence area according to the borehole data set, the material parameter set, a 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 barrier lake dam.

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

[0015] Dividing the barrier lake dam model into an initial influence area and an initial non-influence area;

[0016] Determining 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;

[0017] Adjusting the initial influence area and the initial non-influence area, and determining 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;

[0018] Determining the influence area in the barrier lake dam model based on the first stress, the second stress, and a preset first difference threshold.

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

[0020] Performing a difference operation 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 influence area as the influence area;

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

[0023] Based on the adjusted initial affected area and the adjusted initial non - affected area, perform the steps of adjusting the initial affected area and the initial non - affected area until the first difference exceeds the first difference threshold, and determine the adjusted initial affected area as the affected area.

[0024] In one embodiment, according to the borehole dataset, the material parameter set, a preset second difference threshold, and the finite element algorithm, determining the range of formation boundary changes in the affected area includes:

[0025] Determine the initial range of formation boundary changes in the affected area, and determine formation demarcation points in the initial range of formation boundary changes according to the borehole dataset;

[0026] Based on the formation demarcation points and the initial range of formation boundary changes, generate an initial formation demarcation curve;

[0027] According to the initial formation demarcation curve, the affected area, the material parameter set, a preset second difference threshold, and the finite element algorithm, determine the range of formation boundary changes in the initial range of formation boundary changes.

[0028] In one embodiment, the determining the range of formation boundary changes in the initial range of formation boundary changes according to the initial formation demarcation curve, the affected area, the material parameter set, a preset second difference threshold, and the finite element algorithm includes:

[0029] According to the initial formation demarcation curve, the affected area, the material parameter set, and the finite element algorithm, determine the third stress of the anti - seepage structure in the barrier dam;

[0030] Adjust the initial formation demarcation curve in the direction of the formation demarcation points, and according to the adjusted initial formation demarcation curve, the affected area, the material parameter set, and the finite element algorithm, determine the fourth stress of the anti - seepage structure;

[0031] Based on the third stress, the fourth stress, the adjusted initial formation demarcation curve, and a preset second difference threshold, determine the range of formation boundary changes in the initial range of formation boundary changes.

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

[0033] Perform a difference operation on the third stress and the fourth stress to obtain a second difference, and determine whether the second difference exceeds a preset second difference threshold;

[0034] If the second difference exceeds the second difference threshold, determine a formation boundary change range within the initial formation boundary change range based on the adjusted initial formation demarcation curve.

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

[0036] Execute the step of adjusting the initial formation demarcation curve towards the formation demarcation point according to the adjusted initial formation demarcation curve until the second difference exceeds the second difference threshold, and determine a formation boundary change range within the initial formation boundary change range based on the adjusted initial formation demarcation curve.

[0037] In one embodiment, the initial formation change curve includes various control parameters. According to the initial formation change curve and the finite element algorithm, predict the dam deformation to obtain a 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 a preset loss condition, including:

[0038] 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 prediction according to the optimized formation change curve to obtain a predicted dam deformation data set.

[0039] 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.

[0040] If the loss does not exceed the loss threshold, determine that the loss meets the preset loss condition.

[0041] 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.

[0042] In one embodiment, the step of optimizing 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 includes:

[0043] Determine the perturbation values of each of the control parameters, and perturb each of the control parameters according to the perturbation values of each of the control parameters to generate a perturbed initial formation change curve.

[0044] Based on the finite element algorithm and the perturbed initial formation change curve, dam deformation prediction is carried out to obtain a perturbed dam deformation data set;

[0045] Calculate the perturbation error between the perturbed dam deformation data set and the true dam deformation data set;

[0046] Calculate the perturbation gradient according to the perturbation error, and optimize each control parameter according to the perturbation gradient and the optimization algorithm to generate an optimized initial formation change curve.

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

[0048] An acquisition module for establishing a barrier dam model of the barrier dam and acquiring the material parameter set, borehole data set and true dam deformation data set of the barrier dam;

[0049] A first determination module for determining the formation boundary change range in the barrier dam model based on the material parameter set, the borehole data set and the finite element algorithm, and determining an initial formation change curve within the formation boundary change range;

[0050] A prediction module for performing dam deformation prediction according to the initial formation change curve and the finite element algorithm to obtain a predicted dam deformation data set, and calculating the loss between the predicted dam deformation data set and the true dam deformation data set until the loss meets a preset loss condition;

[0051] A second determination module for determining the initial formation change curve corresponding to the predicted dam deformation data set in the current round as the formation change curve.

[0052] In a third aspect, the present application provides a computer device, including a memory and a processor. A computer program is stored in the memory, and when the processor executes the computer program, the steps in the above method embodiments are implemented.

[0053] In a fourth aspect, the present application provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0054] In a fifth aspect, the present application provides a computer program product, including a computer program. When the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0055] The above post-construction reconstruction method, device, computer equipment and computer-readable storage medium for barrier lakes based on physical mechanisms and survey and monitoring data establish a barrier lake model, and obtain the material parameter set, borehole data set and real dam body deformation data set of the barrier lake; based on the material parameter set, the borehole data set and the finite element algorithm, determine the range of formation boundary changes in the barrier lake model, and determine the initial formation change curve within the range of formation boundary changes; according to the initial formation change curve and the finite element algorithm, perform dam body deformation prediction to obtain a predicted dam body deformation data set, and calculate 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; determine the initial formation change curve corresponding to the predicted dam body deformation data set in the current round as the formation change curve. By using this method, the range of formation boundary changes of the barrier lake is initially determined through the material parameter set, borehole data set and finite element algorithm. Then, with the help of the dam body deformation data set, the initial formation change curve is continuously iteratively optimized, making the initial formation change curve continuously approach the true value until the loss between the predicted dam body deformation data set and the real dam body deformation data set meets the preset loss condition, obtaining an accurate formation change curve, thereby improving the accuracy of the barrier lake formation 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 will briefly introduce the drawings required for use in the description of the embodiments of the present application or related technologies. Obviously, the drawings in the following description are only some embodiments of the present application. For those of ordinary skill in the art, without creative efforts, other related drawings can also be obtained based on these drawings.

[0057] Figure 1 It is a schematic flowchart of a post-construction reconstruction method for a barrier lake based on physical mechanisms and survey and monitoring data in an embodiment;

[0058] Figure 2 It is a schematic flowchart of determining the range of formation boundary changes in an embodiment;

[0059] Figure 3 It is a schematic flowchart of determining the influence area in an embodiment;

[0060] Figure 4 It is a schematic flowchart of determining the influence area according to the first difference in an embodiment;

[0061] Figure 5 It is a schematic diagram of the zoning of the influence level of the formation boundary in an embodiment;

[0062] Figure 6Flow chart for determining influence zones in a barrier dam model in an exemplary embodiment;

[0063] Figure 7 Schematic flow chart for determining the range of formation boundary changes in an influence area in an embodiment;

[0064] Figure 8 Schematic flow chart for determining the range of formation boundary changes in an initial range of formation boundary changes in an embodiment;

[0065] Figure 9 Schematic flow chart for determining the range of formation boundary changes according to a second difference in an embodiment;

[0066] Figure 10 Schematic flow chart for determining the range of formation boundary changes in an exemplary embodiment;

[0067] Figure 11 Schematic flow chart for performing dam body prediction in an embodiment;

[0068] Figure 12 Schematic flow chart for optimizing an initial formation change curve in an embodiment;

[0069] Figure 13 Schematic flow chart for determining a formation change curve in an exemplary embodiment;

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

[0071] Figure 15 Structure block diagram of a post - construction reconstruction device for a barrier dam based on physical mechanisms and survey and monitoring data in an embodiment;

[0072] Figure 16 Internal structure diagram of a computer device in an embodiment. Detailed implementation manners

[0073] In order to make the objectives, technical solutions and advantages of the present application clearer, the present application will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0074] A barrier dam is a natural dam formed by geological disasters blocking a river channel, and its stability is directly related to the safety of the downstream area. To ensure the safety of the downstream area, it is necessary to construct the formation structure of the barrier dam and predict whether the barrier dam is in a safe state based on the formation structure of the barrier dam.

[0075] In the traditional technology, limited drilling data of each borehole in the barrier lake dam is obtained. For two adjacent boreholes, the formation structure where the boreholes are located is determined according to the drilling data of the two adjacent boreholes, and the geostatistical interpolation method is used to predict the formation change curve in the area between the two adjacent boreholes. Then, the formation structure of the barrier lake dam is constructed according to each formation change curve and the formation structure of each borehole.

[0076] However, in the traditional technology, due to the complex geological structure of the barrier lake dam, it is difficult to accurately predict the formation structure in the area between two adjacent boreholes only relying on the geostatistical interpolation method. Therefore, the accuracy of the current method for constructing the barrier lake dam is relatively low. Moreover, the drilling data is spatially sparse, and the formation structure determined based on the drilling data and the geostatistical interpolation method is inaccurate. In addition, the number of monitoring points in the real dam body deformation data set is small and most of them are located on the surface of the dam body, and due to the limitation of the sensor accuracy, there are certain errors in the real dam body deformation data set. It is easy to have the problem of multiple solutions when relying on the real dam body deformation data set to invert the formation structure, that is, there are multiple formation curves that can fit the monitoring data well, and the response of the anti-seepage structure is not considered when relying on the real dam body deformation data set to invert the formation structure.

[0077] Therefore, the present application provides a method for post-construction reconstruction of a barrier lake dam based on physical mechanisms and survey and monitoring data. The change range of the formation boundary of the barrier lake dam is initially determined through a material parameter set, a drilling data set, and a finite element algorithm. Then, with the help of the dam body deformation data set, the initial formation change curve is continuously iteratively optimized, so that the initial formation change curve continuously approaches the true value until the loss between the predicted dam body deformation data set and the real dam body deformation data set meets the preset loss condition, and an accurate formation change curve is obtained, thereby improving the accuracy of the method for post-construction reconstruction of the barrier lake dam based on physical mechanisms and survey and monitoring data. Moreover, based on the stress change of the anti-seepage structure in the barrier lake dam, the material data set, the drilling data set, and the finite element data set, the present application determines the change range of the formation of the barrier lake dam, solves the problem of multiple solutions that is prone to occur when relying on the real dam body deformation data set to invert the formation structure in the prior art. Through the formation change curve in the change range of the formation, not only can the real dam body deformation data set obtained by real survey be fitted, but also the formation change curve is sensitive to the response of the anti-seepage structure, further improving the accuracy of the formation change curve.

[0078] In one embodiment, as Figure 1 shown, a method for post-construction reconstruction of a barrier lake dam based on physical mechanisms and survey and monitoring data is provided. In the embodiments of the present application, taking the application of this method to a computer device as an example for illustration, the embodiments of the present application do not limit the execution device of the method for post-construction reconstruction of the barrier lake dam based on physical mechanisms and survey and monitoring data, and include the following steps 102 to step 108:

[0079] Step 102: Establish a barrier dam model of the barrier lake dam, and obtain the material parameter set, borehole data set, and real dam body deformation data set of the barrier lake dam.

[0080] In implementation, the computer device obtains the attribute data set of the barrier lake dam and establishes the barrier dam model of the barrier lake dam according to the attribute data set of the barrier lake dam. At the same time, the computer device obtains the borehole data set, material parameter set, and real dam body deformation data set of the barrier lake dam.

[0081] Specifically, the computer device obtains the attribute data set of the barrier lake dam, and this attribute data set includes the terrain data, location data, and size data of the barrier lake dam. The terrain data characterizes the detailed undulations of the terrain of the barrier lake dam. The location data includes, but is not limited to, the coordinate data and contour line data of the barrier lake dam. In addition, since an anti-seepage structure is also provided in the barrier lake dam, the computer device obtains the attribute data set of the anti-seepage structure. This attribute data set includes the location data and size data of the anti-seepage structure. The computer device establishes the barrier dam model of the barrier lake dam through the finite element algorithm, the attribute data set of the anti-seepage structure, and the attribute data set of the barrier lake dam. At the same time, the computer device obtains the borehole data set of the barrier lake dam input by the user. Among them, this borehole data set is obtained by the user's on-site inspection and drilling the barrier lake dam at a preset distance to collect the borehole data of each borehole. The borehole data set reflects the stratum structure of the barrier lake dam at the location of the borehole, that is, the type of geological stratum and the stratification position of the geological stratum of the barrier lake dam at the location of the borehole. Sensors monitor the displacement of the dam body of the barrier lake dam to obtain the real dam body deformation data set. The computer device obtains the first material parameter set and the second material parameter set input by the user. Among them, the material parameters in the first material parameter set are the material parameters of the affected area of the barrier lake dam, while the material parameters in the second material parameter set are the material parameters of the non-affected area of the barrier lake dam.

[0082] Optionally, the barrier dam model can be, but is not limited to, a three-dimensional barrier dam model or a two-dimensional barrier dam model. If the user has a high requirement for the accuracy of constructing the stratum structure of the barrier lake dam, a three-dimensional barrier dam model can be established. If the user has a high requirement for the efficiency of constructing the stratum structure of the barrier lake dam, a two-dimensional barrier dam model can be established, that is, the dimension of the established barrier dam model is determined according to the construction requirements of the stratum of the barrier lake dam. The embodiments of the present application do not limit the dimension of the barrier dam model. The finite element algorithm is implemented through finite element software, and the finite element software can be, but is not limited to, abaqus software (a set of powerful finite element software for engineering simulation). The embodiments of the present application do not limit the finite element software.

[0083] Step 104: Based on the material parameter set, borehole data set, and finite element algorithm, determine the change range of the stratum boundary in the barrier dam model, and determine the initial stratum change curve within the change range of the stratum boundary.

[0084] Among them, the range of formation boundary change is the envelope range of the formation boundary and also the interval where the formation boundary is located.

[0085] In implementation, the computer device determines the influence area in the barrier lake dam model according to the material parameter set, the finite element algorithm, and a preset first difference threshold, and determines the range of formation boundary change in the influence area according to the finite element algorithm, the borehole dataset, and a preset second difference threshold. Then, the computer device generates an initial formation change curve in the range of formation boundary change by a curve generation algorithm (a Python algorithm).

[0086] Specifically, in order to reinforce the barrier lake dam, an anti-seepage structure is set in the barrier lake dam. The computer device determines the dam body area in the barrier lake dam model that has a greater influence on the stress response of the anti-seepage structure according to the material parameter set, the finite element algorithm, and the first difference threshold, and determines this dam body area as the influence area. Then, the computer device determines the initial range of formation boundary change in the influence area, and narrows down the initial range of formation boundary change according to the borehole dataset, the finite element algorithm, and a preset second difference threshold to obtain the range of formation boundary change. The computer device determines each control parameter according to the range of formation boundary change and the simulation requirements, and generates an initial formation boundary change curve according to each control parameter and the curve generation algorithm.

[0087] In an exemplary embodiment, the computer device divides the barrier lake dam model into an initial influence area and an initial non-influence area, and determines 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. Then, the computer device adjusts the initial influence area and the initial non-influence area to obtain the adjusted initial influence area and the adjusted initial non-influence area. The computer device determines the second stress of the anti-seepage structure according to the adjusted initial influence area, the adjusted initial non-influence area, the material parameter set, and the finite element algorithm. The computer device calculates the first difference between the first stress and the second stress until the first difference meets the preset first difference threshold, and the computer device determines the adjusted initial influence area of the current round as the influence area. The computer device determines the initial range of formation boundary change in the influence area, and generates an initial formation boundary curve in the initial range of formation boundary change according to the borehole dataset. The computer device narrows down the initial range of formation boundary change according to the initial formation boundary curve, the influence area, the material parameter set, the second difference threshold, and the finite element algorithm to obtain the range of formation boundary change. The computer device determines each control parameter according to the range of formation boundary change and the simulation requirements, and generates an initial formation boundary change curve according to each control parameter.

[0088] Optionally, the number of control parameters is determined according to the accuracy requirement of the stratum structure of the barrier dam. The more the number of control parameters, the more accurate the simulation of the stratum boundary change curve, and the higher the accuracy of the stratum structure of the barrier dam; the fewer the number of control parameters, the rougher the simulation of the stratum boundary change curve, and the lower the accuracy of the stratum structure of the barrier dam.

[0089] Step 106: According to the initial stratum change curve and the finite element algorithm, perform dam body deformation prediction to obtain a predicted dam body deformation data set, and calculate the loss between the predicted dam body deformation data set and the real dam body deformation data set until the loss meets the preset loss condition.

[0090] Among them, the initial stratum change curve contains each control parameter. The predicted dam body deformation data sets that do not meet the loss condition are all predicted from the initial stratum change curve.

[0091] In implementation, the loss condition is preset in the computer device. The computer device optimizes the initial stratum change curve according to each control parameter, and performs dam body deformation prediction according to the optimized initial stratum change curve and the finite element algorithm to obtain a predicted dam body deformation data set. The computer device calculates the loss between the predicted dam body deformation data set and the real dam body 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 elaborated in detail in the following embodiments, and will not be repeated in the embodiments of the present application. If the loss does not meet the loss condition, the computer device continues to optimize the initial stratum change curve until the loss meets the loss condition.

[0092] Specifically, the loss condition preset 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 stratum change curve according to each control parameter, the finite element algorithm and the optimization algorithm to obtain an optimized initial stratum change curve. Then, the computer device performs dam body deformation prediction according to the optimized initial stratum change curve and the finite element algorithm to obtain a predicted dam body deformation data set, and calculates the error between the predicted dam body deformation data set and the real dam body deformation data set. This error is also the loss between the predicted dam body deformation data set and the real dam body deformation data set. Then, the computer device determines whether the error exceeds the error threshold. The computer device determines that the error does not meet the loss condition. In the case where the error does not meet the loss condition, the computer device executes the step of optimizing each control parameter until the loss meets the loss condition, and the computer device obtains an optimized initial stratum change curve.

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

[0094] In implementation, the computer device determines the optimized initial formation change curve corresponding to the predicted dam deformation data set of the current round as the formation change curve. Then, the computer device constructs the formation structure of the barrier lake dam based on the borehole data set, the formation change curve, and the barrier lake dam model.

[0095] Specifically, the computer device determines the optimized initial formation change curve corresponding to the predicted dam deformation data set of the current round as the formation change curve. Then, the computer device determines two geological layers corresponding to the formation change curve according to the borehole data set, and constructs the formation structure of the barrier lake dam based on the two geological layers, the formation change curve, and the barrier lake dam model.

[0096] In the above method for post-construction reconstruction of a barrier lake dam based on physical mechanisms and survey and monitoring data, the change range of the formation boundary of the barrier lake dam is initially determined through the material parameter set, the borehole data set, and the finite element algorithm. Then, the initial formation change curve is continuously iteratively optimized with the help of the dam deformation data set, making the initial formation change curve continuously approach the true value until the loss between the predicted dam deformation data set and the true dam deformation data set meets the preset loss condition, obtaining an accurate formation change curve, thereby improving the accuracy of the formation reconstruction method for the barrier lake dam.

[0097] In an exemplary embodiment, as Figure 2 shown, the specific processing procedure for determining the change range of the formation boundary in the barrier lake 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. Among them:

[0098] Step 202, based on the material parameter set, the finite element algorithm, and a preset first difference threshold, determines the affected area in the barrier lake dam model.

[0099] Among them, the first difference threshold is the stress change threshold of the anti-seepage structure in the barrier lake 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 implementation, the affected area is the area in the barrier dam body that significantly affects the stress response of the anti-seepage structure, and the non-affected area is the area in the barrier dam body that has an insignificant effect on the stress response of the anti-seepage structure. Since both the affected area and the non-affected area are related to the stress of the anti-seepage structure, therefore, this application uses the change in the stress of the anti-seepage structure (the first difference threshold) to verify whether the affected area and the non-affected area are correct. Specifically, the computer device divides the barrier dam model to obtain the initial affected area and the initial non-affected area. Then, the computer device determines the first stress of the anti-seepage structure in the barrier dam body according to the material parameter set, the finite element algorithm, the initial affected area, and the initial non-affected area. The computer device reduces the initial affected area to obtain the adjusted initial affected area, and adjusts the initial non-affected area according to the adjusted initial affected area to obtain the adjusted initial non-affected area. The computer device determines the second stress of the anti-seepage structure in the barrier dam body according to the material parameter set, the finite element algorithm, the adjusted initial affected area, and the adjusted initial non-affected area. Then, the computer device determines the 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 affected area of the current round as the affected area.

[0101] Step 204, determine the change range of the formation boundary in the affected area according to the borehole dataset, the material parameter set, the preset second difference threshold, and the finite element algorithm.

[0102] Among them, the second difference threshold is the stress change threshold of the anti-seepage structure in the barrier dam.

[0103] In implementation, the computer device determines the initial change range of the formation boundary in the affected area, and generates an initial formation boundary curve according to the initial change range of the formation boundary and the borehole dataset. Then, the computer device reduces the initial change range of the formation boundary according to the initial formation boundary curve, the affected area, the material parameter set, the second difference threshold, and the finite element algorithm to obtain the change range of the formation boundary.

[0104] In this embodiment, by using the material dataset, the borehole dataset, the finite element algorithm, and the stress change threshold of the anti-seepage structure, the change range of the formation boundary is determined in the barrier dam model, which provides a constraint condition for the subsequent optimization of the initial formation change curve, reduces the computational complexity during optimization, and improves the accuracy of the barrier dam formation reconstruction method. Moreover, by determining the change range of the formation of the barrier dam, the problem of multiple solutions easily occurring in the inversion of the formation structure relying on the real dam body deformation dataset in the prior art is solved. Through the formation change curve in the change range of the formation, not only can the real dam body deformation dataset obtained by real survey be fitted, but also the formation change curve is sensitive to the response of the anti-seepage structure, further improving the accuracy of the formation change curve.

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

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

[0107] In implementation, the computer device divides the barrier dam model into an initial affected area and an initial non-affected area according to the division instruction.

[0108] Specifically, the user operates the computer device based on their own engineering experience to send a division instruction to the computer device. The division instruction includes each control point and the position information of each control point. The computer device receives the division instruction and divides the entire barrier dam model into an initial affected area and an initial non-affected area according to the control points and the position information of each control point in the division instruction.

[0109] Step 304: Determine 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 through finite element software. The material parameter set includes a first material parameter set of the affected area and a second material parameter set of the non-affected area. A stress algorithm is set in the finite element software. The computer device divides the initial affected area into each initial affected area grid according to the first division size corresponding to the affected area. The computer device divides the initial non-affected area into each initial non-affected area grid according to the second division size corresponding to the non-affected area. The computer device performs data processing on each initial affected area grid, each initial non-affected 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 affected area grid is smaller than the initial non-affected area grid. That is, the computer device divides the initial affected area with a higher density and divides the initial non-affected area with a lower density.

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

[0112] In implementation, a stress algorithm is set in the finite element software. The computer device reduces the initial influence area by adjusting the positions of each control point to obtain the adjusted initial influence area. Then, the computer device determines the initial non-influence area as the area outside the adjusted initial influence area in the barrier lake dam model. The computer device divides the adjusted initial influence area into each adjusted initial influence area grid according to the first division size corresponding to the influence area. The computer device divides the adjusted initial non-influence area into each adjusted initial non-influence area grid according to the second division size corresponding to the non-influence area. The computer device performs data processing on each adjusted initial influence area grid, each adjusted initial non-influence 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, generally three control points are set. Each time the initial influence area is adjusted, at least the position of one control point is adjusted, which is determined according to the adjustment requirements. The embodiments of the present application do not limit the number of control points and the number of adjusted control points.

[0114] Step 308: Based on the first stress, the second stress, and a preset first difference threshold, determine the influence area in the barrier lake dam model.

[0115] In implementation, the computer device determines the first difference between the first stress and the second stress and determines whether the first difference threshold meets the preset first difference threshold. If the first difference meets the first difference threshold, the computer device determines the adjusted initial influence area as the influence area in the barrier lake dam model. If the first difference does not meet the first difference threshold, the computer device updates the first stress according to the second stress and executes the above step 306 until the first difference meets the first difference threshold, and the computer device determines the influence area in the barrier lake dam model.

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

[0117] Step 402: Perform a difference process on the second stress and the first stress to obtain the first difference and determine whether the first difference exceeds the preset first difference threshold.

[0118] In implementation, a first difference threshold is preset in the computer device. The computer device performs a difference process on the second stress and the first stress to obtain a first difference. This first difference characterizes the change in the stress distribution of the anti-seepage structure. The computer device determines whether the first difference exceeds the first difference threshold, that is, the computer device 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 according to the attribute information of the anti-seepage structure and the barrier lake dam, and the embodiments of the present application do not limit the first difference threshold.

[0120] Step 404, if the first difference exceeds the first difference threshold, determine the adjusted initial influence area as the influence area.

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

[0122] In an exemplary embodiment, Figure 5 is a schematic diagram of the influence level zoning of the formation boundary in an embodiment. As Figure 5 shown, the left figure shows the barrier lake dam model of the barrier lake dam, and the anti-seepage structure in the barrier lake dam is the anti-seepage wall. The right figure shows the barrier lake dam model in which the influence area and the non-influence area have been 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 the 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 influence area and the adjusted initial non-influence area, execute the steps of adjusting the initial influence area and the initial non-influence area until the first difference exceeds the first difference threshold, and determine the adjusted initial influence area as the influence area.

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

[0127] In one exemplary embodiment, to improve the optimization efficiency and provide a basis for subsequent exploration and monitoring layout, the present application introduces a method for zoning the formation influence level based on the response of the anti - seepage structure. Through this method and two material parameters with large differences in properties, the dam body area with a greater influence on the stress response of the anti - seepage structure is identified as the affected area. Specifically, Figure 6 is a flowchart for determining the influence zoning in the barrier lake dam model in one exemplary embodiment. As Figure 6 shown, the process of determining the influence zoning in the barrier lake dam model includes:

[0128] Step 601, generate the boundary line of the affected area, and divide the initial affected area and the initial non - affected area in the barrier lake dam model according to engineering experience and the boundary line of the affected area; the initial affected area corresponds to the first set of material parameters; the initial non - affected area corresponds to the second set of material parameters;

[0129] Step 602, determine the first stress of the anti - seepage structure through the first set of material parameters, the second set of material parameters, the initial affected area, the initial non - affected area, and the finite element algorithm; the stress of the anti - seepage structure characterizes the stress distribution of the anti - seepage structure;

[0130] Step 603, reduce the initial affected area by adjusting the boundary of the affected area to obtain the adjusted initial affected area, and determine the adjusted initial non - affected area according to the adjusted initial affected area;

[0131] Step 604, determine the second stress of the anti - seepage structure through the first set of material parameters, the second set of material parameters, the adjusted initial affected area, the adjusted initial non - affected area, and the finite element algorithm;

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

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

[0134] Step 607, determine the second stress as the first stress, and execute the above - mentioned step 603;

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

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

[0137] Step 702: Determine the change range of the initial formation boundary in the affected area, and determine the formation demarcation points in the change range of the initial formation boundary according to the borehole dataset.

[0138] Among them, the borehole dataset includes the positions of each borehole and the geological stratification positions corresponding to the positions of each borehole.

[0139] In implementation, the computer device determines the change range of the initial formation boundary in the affected area. The computer device determines the formation demarcation points in the change range of the initial formation boundary according to the geological stratification positions.

[0140] Specifically, a mapping relationship between the geodetic coordinates and the model coordinates of the barrier dam model is set in the computer device. The user operates the computer device according to engineering experience and sends a range determination instruction to the computer device. The range determination instruction includes each range coordinate. The computer device receives the range determination instruction and preliminarily determines the change range of the initial formation boundary in the affected area according to each range coordinate in the range determination instruction. The change range of the initial formation boundary is the range or interval where the formation may change preliminarily determined. The computer device determines the formation demarcation points in the change range of the initial formation boundary according to the geological stratification positions corresponding to the positions of each borehole and the mapping relationship.

[0141] Step 704: Generate an initial formation demarcation curve based on the formation demarcation points and the change range of the initial formation boundary.

[0142] Among them, the formation demarcation points include the first formation demarcation point and the second formation demarcation point.

[0143] In implementation, the computer device obtains each formation demarcation curve control parameter. Then, the computer device takes the first formation demarcation point as the starting point and the second formation demarcation point as the ending point, and generates an initial formation demarcation curve in the change range of the initial formation boundary through the finite element algorithm and each formation demarcation curve control parameter. The initial formation demarcation curve is a B - spline curve, and the initial formation demarcation curve is used to simulate the formation boundary of the barrier dam.

[0144] In an exemplary embodiment, since the initial formation boundary change range has an upper limit of change in the y - direction and a lower limit of change in the y - direction, in order to construct the formation boundary change range through the initial formation demarcation curve, two initial formation demarcation curves need to be generated. Therefore, the computer device acquires the formation demarcation curve control parameters. Among them, the formation demarcation curve control parameters include the first demarcation control curve parameters and the second demarcation control curve parameters. The computer device takes the first formation demarcation point as the starting point and the second formation demarcation point as the ending point, and generates the first initial formation demarcation curve in the initial formation boundary change range through the finite - element algorithm and each first formation demarcation curve control parameter. The first initial formation demarcation curve coincides with the upper limit of change in the y - direction of the initial formation boundary change range. The computer device takes the first formation demarcation point as the starting point and the second formation demarcation point as the ending point, and generates the second initial formation demarcation curve in the initial formation boundary change range through the finite - element algorithm and each second formation demarcation curve control parameter. The second initial formation demarcation curve coincides with the lower limit of change in the y - direction of the initial formation boundary change range.

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

[0146] Step 706, determine the formation boundary change range in the initial formation boundary change range according to the initial formation demarcation curve, the influence area, the material parameter set, the preset second difference threshold, and the finite - element algorithm.

[0147] In implementation, the computer device determines the third stress of the anti - seepage structure according to the initial formation demarcation curve, the influence area, the material parameter set, and the finite - element algorithm. Then, the computer device adjusts the initial formation demarcation curve in the direction of the formation demarcation point to obtain the adjusted initial formation demarcation curve. The computer device determines the fourth stress of the anti - seepage structure according to the adjusted initial formation demarcation curve, the influence area, the material parameter set, and the finite - element algorithm. Then, the computer device performs a difference process on the third stress and the fourth stress to obtain a second difference, and determines whether the second difference exceeds the second difference threshold. If the second difference exceeds the second difference threshold, the computer device determines the formation boundary change range in the initial formation boundary change range according to the adjusted initial formation demarcation curve.

[0148] In an exemplary embodiment, since the barrier dam model includes an affected area and a non-affected area, the computer device needs to determine the formation boundary of the barrier dam in the affected area and the non-affected area respectively. The computer device determines the formation boundary in the non-affected area according to the borehole dataset and the linear interpolation algorithm. Then, the computer device determines the initial formation boundary change range in the affected area, and determines the formation boundary change range in the initial formation boundary change range according to the change of the anti-seepage structure stress. Specifically, the computer device determines the third stress of the anti-seepage structure according to the initial formation boundary curve, the affected area, the material parameter set and the finite element algorithm. Then, the computer device adjusts the initial formation boundary curve in the direction of the formation boundary point to obtain the adjusted initial formation boundary curve. The computer device determines the fourth stress of the anti-seepage structure according to the adjusted initial formation boundary curve, the affected area, the material parameter set and the finite element algorithm. Then, the computer device performs a difference process on the third stress and the fourth stress 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 formation boundary change range in the initial formation boundary change range according to the adjusted initial formation boundary curve.

[0149] Optionally, the computer device can process each initial formation boundary curve one by one to obtain the formation boundary curve. It can also process each initial formation boundary curve simultaneously to obtain the formation boundary curve.

[0150] In this embodiment, by using the material dataset, borehole data and finite element algorithm, the formation boundary change range is determined in the affected area, which provides constraint conditions for subsequent optimization of the initial formation change curve, reduces the computational complexity during optimization, reduces the execution time of the post-construction reconstruction method for barrier dams based on physical mechanisms and survey and monitoring data, and improves the efficiency of the post-construction reconstruction method for barrier dams based on physical mechanisms and survey and monitoring data.

[0151] In an exemplary embodiment, as Figure 8 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 barrier dam according to the initial formation boundary curve, the affected area, the material parameter set and the finite element algorithm.

[0153] Among them, the material parameter set includes, but is not limited to, a first material parameter set and a second material parameter set.

[0154] In implementation, the computer device determines the area other than the affected area as the non - affected area in the barrier lake dam model according to the affected area. Then, the computer device imports the material parameter set, the affected area, the non - affected area, and the initial stratigraphic boundary curve into the finite - element software corresponding to the finite - element algorithm. Through the finite - element software, data processing is performed on the material parameter set, the affected area, the non - affected area, and the initial stratigraphic boundary curve, and the stress of the anti - seepage structure is solved to obtain the third stress of the anti - seepage structure. The third stress characterizes the stress distribution of the anti - seepage structure in the barrier lake dam when the geological stratification of the barrier lake dam is the initial stratigraphic boundary curve.

[0155] In an exemplary embodiment, the initial stratigraphic boundary curve includes a first initial stratigraphic boundary curve and a second initial stratigraphic boundary curve. The computer device determines the area other than the affected area as the non - affected area in the barrier lake dam model according to the affected area. Then, the computer device imports the material parameter set, the affected area, the non - affected area, and the first initial stratigraphic boundary curve into the finite - element software corresponding to the finite - element algorithm. Through the finite - element software, data processing is performed on the material parameter set, the affected area, the non - affected area, and the first initial stratigraphic boundary curve, and the stress of the anti - seepage structure is solved to obtain the third stress corresponding to the first initial stratigraphic boundary curve. The computer device imports the material parameter set, the affected area, the non - affected area, and the second initial stratigraphic boundary curve into the finite - element software corresponding to the finite - element algorithm. Through the finite - element software, data processing is performed on the material parameter set, the affected area, the non - affected area, and the second initial stratigraphic boundary curve, and the stress of the anti - seepage structure is solved to obtain the third stress corresponding to the second initial stratigraphic boundary curve.

[0156] Step 804: Adjust the initial stratigraphic boundary curve in the direction of the stratigraphic boundary point, and determine the fourth stress of the anti - seepage structure according to the adjusted initial stratigraphic boundary curve, the affected area, the material parameter set, and the finite - element algorithm.

[0157] In implementation, the computer device adjusts the initial stratigraphic boundary curve in the direction of the stratigraphic boundary point, so that the adjusted initial stratigraphic boundary curve is closer to the stratigraphic boundary point as a whole. The computer device imports the material parameter set, the affected area, the non - affected area, and the adjusted initial stratigraphic boundary curve into the finite - element software corresponding to the finite - element algorithm. Through the finite - element software, data processing is performed on the material parameter set, the affected area, the non - affected area, and the adjusted initial stratigraphic boundary curve, and the stress of the anti - seepage structure is solved to obtain the fourth stress of the anti - seepage structure. The fourth stress characterizes the stress distribution of the anti - seepage structure in the barrier lake dam when the geological stratification of the barrier lake dam is the adjusted initial stratigraphic boundary curve.

[0158] In an exemplary embodiment, the initial formation boundary curve includes a first initial formation boundary curve and a second initial formation boundary curve. The first initial formation boundary curve coincides with the upper limit of the change in the y-direction of the initial formation boundary change range. The second initial formation boundary curve coincides with the lower limit of the change in the y-direction of the initial formation boundary change range. The computer device adjusts the first initial formation boundary curve in the direction of the formation boundary point and adjusts the second initial formation boundary curve in the direction of the formation boundary point. At this time, the interval formed by the adjusted first initial formation boundary curve and the adjusted second initial formation boundary curve is smaller than the interval formed by the first initial formation boundary curve and the second initial formation boundary curve. Therefore, by continuously adjusting the initial formation boundary curve in the direction of the formation boundary point, the initial formation boundary change range can be continuously reduced. 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 formation boundary curve into the finite element software corresponding to the finite element algorithm, and the finite element software processes the data of the first material parameter set, the second material parameter set, the affected area, the non-affected area, and the adjusted first initial formation boundary curve to solve the stress of the anti-seepage structure and obtain the fourth stress corresponding to the adjusted first initial formation boundary curve. The computer device imports the material parameter set, the affected area, the non-affected area, and the adjusted second initial formation boundary curve into the finite element software corresponding to the finite element algorithm, and the finite element software processes the data of the material parameter set, the affected area, the non-affected area, and the adjusted second initial formation boundary curve to solve the stress of the anti-seepage structure and obtain the fourth stress of the adjusted second initial formation boundary curve.

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

[0160] In implementation, the computer device performs a difference process on the third stress and the fourth stress to obtain a second difference. The computer device determines whether the second difference exceeds the preset second difference threshold. If the second difference exceeds the second difference threshold, the computer device determines the formation boundary change range within the initial formation boundary change 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 according to the fourth stress and continues to execute the above step 804 until the second difference exceeds the second difference threshold, and then the computer device determines the formation boundary change range within the initial formation boundary change range according to the adjusted initial formation boundary curve.

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

[0162] Step 902: Perform a difference operation on the third stress and the fourth stress to obtain a second difference, and determine whether the second difference exceeds a preset second difference threshold.

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

[0164] Optionally, the second difference threshold is determined according to the attribute information of the barrier dam, and the 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 change range within the initial formation boundary change range based on the adjusted initial formation boundary curve.

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

[0167] Specifically, the initial formation boundary curve includes the first initial formation boundary curve and the second initial formation boundary curve. If the second difference exceeds the second difference threshold, the computer device determines the adjusted first initial formation boundary curve as the first formation boundary curve and determines the adjusted second initial formation boundary curve as the second formation boundary curve. Since the starting point and the ending point of the first formation boundary curve and the second formation boundary curve are both formation boundary points, and the first formation boundary curve and the second formation boundary curve are both within the initial formation boundary change range, therefore, the first formation boundary curve and the second formation boundary curve can enclose an area within the initial formation boundary change range. The computer device determines the area enclosed by the first formation boundary curve and the second formation boundary curve as the formation boundary change 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] In 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, and obtains 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, and obtains the updated third stress corresponding to the second initial formation boundary curve.

[0170] Step 908: According to the adjusted initial formation boundary curve, perform the step of adjusting the initial formation boundary curve towards the formation boundary point until the second difference exceeds the second difference threshold, and based on the adjusted initial formation boundary curve, determine the formation boundary change range within the initial formation boundary change range.

[0171] In implementation, the computer device continues to execute the above-mentioned step 804 according to the adjusted initial formation boundary curve until the second difference exceeds the second difference threshold. The computer device determines the adjusted initial formation boundary curve as the formation boundary curve. Then, the computer device determines the area enclosed by the formation boundary curve as the formation boundary change range within the initial formation boundary change range. Among them, the specific processing process of step 804 has been elaborated in detail in the above-mentioned embodiments, and will not be repeated in the embodiments of the present application.

[0172] In an exemplary embodiment, Figure 10 is a schematic flowchart for determining the formation boundary change range in an exemplary embodiment. As Figure 10As shown, the computer device preliminarily determines the possible change range of the formation boundary within the influence area range to obtain the initial formation boundary change range. Then, the computer device determines formation demarcation points within the initial formation boundary change range through the borehole dataset. The computer device uses the formation demarcation points as the two ends of the initial formation demarcation curve and simulates the formation boundary within the initial formation boundary change range using a B-spline curve to obtain the initial formation demarcation curve. Then, the computer device determines the third stress of the anti-seepage structure in the barrier dam according to the initial formation demarcation curve, the influence area, the material parameter set, and the finite element algorithm. The computer device adjusts the initial formation demarcation curve in the y direction towards the formation demarcation point, making the adjusted initial formation demarcation curve continuously approach the formation demarcation point, and determines the fourth stress of the anti-seepage structure according to the adjusted initial formation demarcation curve, the influence area, the material parameter set, and the finite element algorithm. The computer device determines the second difference between the third stress and the fourth stress. This second difference characterizes the influence of the adjusted initial formation demarcation curve on the stress response of the anti-seepage structure. The computer device determines whether the second difference exceeds a preset second difference threshold to determine whether the adjusted initial formation demarcation curve has a significant impact on the stress of the anti-seepage structure. If the second difference exceeds the second difference threshold, the computer device determines the adjusted initial formation demarcation curve as the formation demarcation curve, which is the formation boundary curve that first has a significant impact on the stress response of the anti-seepage structure. Then, the computer device determines the interval enclosed by the formation demarcation curve within the initial formation boundary change range as the formation boundary change range (formation boundary envelope range). 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 the computer device adjusting the initial formation demarcation curve in the y direction towards the formation demarcation point until the second difference exceeds the second difference threshold, and the computer device determines the interval enclosed by the formation demarcation curve within the initial formation boundary change range as the formation boundary change 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, optimize the initial formation change curve according to various 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 the predicted dam body deformation dataset.

[0175] Among them, the real dam body deformation dataset includes the real displacement data of the barrier dam during the historical time period.

[0176] In implementation, the computer device obtains each initial control parameter and the disturbance value corresponding to each control parameter, and perturbs each control parameter according to each disturbance value to obtain the perturbed initial formation change curve. Then, the computer device optimizes the control parameters based on the perturbed initial formation change curve and the true predicted dam deformation data set to obtain each optimized control parameter, and generates an optimized initial formation change curve according to each optimized control parameter. The computer device processes the optimized initial formation change curve, the material parameter set, the borehole data set, the affected area and the non-affected area according to the finite element software corresponding to the finite element algorithm, so as to predict the displacement of the barrier lake dam and obtain the predicted dam deformation data.

[0177] Step 1104, calculate the loss between the predicted dam deformation data set and the true dam deformation data set, and determine whether the loss exceeds a preset loss threshold.

[0178] Among them, the loss condition is that the loss between the predicted dam deformation data set and the true dam deformation data set does not exceed the loss threshold.

[0179] In implementation, a loss threshold is preset in the computer device. The computer device performs data operations on the predicted dam deformation data set and the true dam deformation data set according to the preset loss algorithm to obtain the loss between the predicted dam deformation data set and the true dam deformation data set. Then, the computer device determines whether the current loss exceeds the loss threshold. If the loss exceeds the loss threshold, the computer device executes the following step 1108. If the loss does not exceed the loss threshold, the computer device executes the following step 1106.

[0180] Optionally, the loss threshold is determined according to the simulation accuracy requirement of the formation change curve. If the simulation accuracy requirement of the formation change curve is high, the loss threshold is set small; if the simulation accuracy requirement of the formation change curve is low, the loss threshold is set large. The embodiments of the present application do not limit the loss threshold.

[0181] Step 1106, if the loss does not exceed the loss threshold, determine 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 each control parameter, 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.

[0184] In implementation, if the loss exceeds the loss threshold, the computer device determines that the loss does not meet the preset loss condition. When the loss does not meet the loss condition, the computer device continues to execute step 1102 above until the loss does not exceed the loss threshold, that is, the loss meets the preset loss condition. Among them, the specific processing process of step 1102 has been elaborated in detail in the above embodiments, and will not be repeated in the embodiments of the present application.

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

[0186] Step 1202: Determine the perturbation values of the various control parameters, and perturb the various control parameters according to the perturbation values of the various control parameters to generate a perturbed initial formation change curve.

[0187] In implementation, the computer device determines the perturbation values of the various control parameters. For example, the user inputs the perturbation values of the various control parameters to the computer device by operating the computer device. Among them, the perturbation values include positive perturbation values and negative perturbation values. The computer device updates the control parameters according to the perturbation value of each control parameter to obtain the perturbed control parameters. The perturbed control parameters include the various control parameters after positive perturbation and the various control parameters after negative perturbation. Then, the computer device uses a B-spline curve to simulate the formation boundary of the barrier dam, and generates a perturbed initial formation change curve representing the formation boundary based on the perturbed control parameters. The perturbed initial formation change curve includes a positively perturbed initial formation change curve and a negatively perturbed initial formation change curve.

[0188] Optionally, the perturbation values of each control parameter can be different or the same, which is determined according to the adjustment requirements of the initial formation change curve, and the perturbation values of the various control parameters of the present application are not limited.

[0189] Step 1204: Based on the finite element algorithm and the perturbed initial formation change curve, perform dam body deformation prediction to obtain a perturbed dam body deformation data set.

[0190] In implementation, the computer device imports the initial formation change curve after positive perturbation, the material parameter set, the borehole data set, the affected area and the non - affected area into the finite - element software corresponding to the finite - element algorithm. Through the finite - element software, data processing is performed on the initial formation change curve after positive perturbation, the material parameter set, the borehole data set, the affected area and the non - affected area to predict the displacement of the barrier dam within the historical time period, and a positive - perturbation dam - body deformation data set is obtained. The computer device imports the initial formation change curve after negative perturbation, the material parameter set, the borehole data set, the affected area and the non - affected area into the finite - element software corresponding to the finite - element algorithm. Through the finite - element software, data processing is performed on the initial formation change curve after negative perturbation, the material parameter set, the borehole data set, the affected area and the non - affected area to predict the displacement of the barrier dam within the historical time period, and a negative - perturbation dam - body deformation data set is obtained.

[0191] Step 1206: Calculate the perturbation error between the perturbed dam - body deformation data set and the real dam - body deformation data set.

[0192] In implementation, the computer device performs data processing on the perturbed dam - body deformation data set and the real dam - body deformation data set according to a preset loss algorithm to obtain the perturbation error between the perturbed dam - body deformation data set and the real dam - body deformation data set.

[0193] Specifically, the perturbed dam - body deformation data set includes the positive - perturbation dam - body deformation data set and the negative - perturbation dam - body deformation data set. The computer device performs data processing on the positive - perturbation dam - body deformation data set and the real dam - body deformation data set according to a preset loss algorithm to obtain the positive - perturbation error between the positive - perturbation dam - body deformation data set and the real dam - body deformation data set. The computer device performs data processing on the negative - perturbation dam - body deformation data set and the real dam - body deformation data set according to a preset loss algorithm to obtain the negative - perturbation error between the negative - perturbation dam - body deformation data set and the real dam - body deformation data set.

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

[0195] In implementation, the optimization algorithm is the gradient - descent method. The computer device calculates the perturbation gradient of each control parameter according to the central - difference method and the perturbation error, and updates the control parameters according to the perturbation gradient and the gradient - descent method to obtain updated control parameters. The computer device determines whether the updated control parameters are within the formation - boundary change range to obtain a judgment result. Then, the computer device generates an optimized formation change curve based on each judgment result, the perturbation values of each control parameter, and each updated control parameter.

[0196] Specifically, the computer device calculates the approximate value of the partial derivative of each control parameter according to the central difference method, and calculates the perturbation gradient of the control parameter according to the approximate value of the partial derivative. Then, the computer device updates the control parameter according to the gradient descent method and the perturbation gradient to obtain the updated control parameter. The computer device determines whether the updated control parameter is within the formation boundary change range to obtain a judgment result. If the judgment result indicates that the updated control parameter is within the formation boundary change 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 change range, the computer device adjusts the perturbation value corresponding to the control parameter, and re-updates the control parameter according to the adjusted perturbation value so that the updated control parameter is within the formation boundary change range. Then, the computer device generates an optimized formation change curve within the formation boundary change range according to each updated control parameter.

[0197] In one exemplary embodiment, Figure 13 is a schematic flowchart for determining a formation change curve in an exemplary embodiment. As Figure 13 shown, the specific processing process for 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 parameter is represented by ;

[0199] Step 1302, obtain the perturbation value of each control parameter;

[0200] Step 1303, perturb the control parameter according to the perturbation value of the control parameter to obtain the perturbed control parameter; where the perturbed control parameter is ; is the th control parameter, is the perturbation value;

[0201] Step 1304, generate an optimized perturbed formation change curve according to each perturbed control parameter;

[0202] Step 1305, based on the finite element algorithm and the perturbed initial formation change curve, perform dam deformation prediction to obtain a perturbed dam deformation data set;

[0203] Step 1306, calculate the perturbation error and between the perturbed dam deformation data set and the real dam deformation data set;

[0204] Step 1307, calculate the perturbation gradient of each control parameter according to the perturbation error;

[0205] Step 1308: Update the control parameters according to the perturbation gradient to obtain the updated control parameters; the updated control parameters ; where is the th updated control parameter, is the th control parameter, is the perturbation gradient, is the learning rate;

[0206] Step 1309: Determine whether the updated control parameters are within the formation boundary change range; if the updated control parameters are within the formation boundary change range, then execute Step 1310, if the updated control parameters are not within the formation boundary change range, then execute Step 1302;

[0207] Step 1310: Generate an optimized formation change curve based on each updated control parameter;

[0208] Step 1311: Based on the finite element algorithm and the optimized initial formation change curve, perform dam deformation prediction to obtain a predicted dam deformation data set;

[0209] Step 1312: Calculate the loss between the predicted dam deformation data set and the real dam deformation data set;

[0210] Step 1313: Determine whether the loss exceeds a preset loss threshold. If the loss exceeds the loss threshold, then execute Step 1302. If the error does not exceed the preset error threshold, then 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 is a flowchart of a post-construction reconstruction method for a barrier dam based on physical mechanisms and survey and monitoring data in an exemplary embodiment, as Figure 14 shown. The specific processing process of the post-construction reconstruction method for a barrier dam based on physical mechanisms and survey and monitoring data includes;

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

[0214] Step 1402: Determine the formation boundary change range in the influence area based on the borehole information in the borehole data set;

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

[0216] Among them, the specific processing process of step 1401 includes establishing a barrier dam model and performing a sensitivity analysis on the response of the anti-seepage structure to determine the affected area and non-affected area in the barrier dam model. The specific processing process of step 1402 includes determining the initial formation boundary change range (formation boundary variable range) in the affected area according to the borehole dataset, and continuously changing the boundary of the formation boundary variable range through the anti-seepage structure response analysis to obtain the formation boundary change range (formation boundary envelope range). The specific processing process of step 1403 includes continuously changing the formation boundary within the formation boundary envelope range based on the initial formation change curve and performing finite element calculation of the barrier dam, so that the predicted dam body deformation dataset fits the real dam body deformation dataset (monitoring data) to obtain the formation change curve (optimal formation boundary).

[0217] It should be understood that although the steps in the flowcharts involved in the above-described embodiments are sequentially shown according to the arrows, these steps are not necessarily executed in the order indicated by the arrows.

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

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

[0220] The first determination module 1502 is used to determine the formation boundary change range in the barrier dam model based on the material parameter set, borehole dataset, and finite element algorithm, and determine the initial formation change curve in the formation boundary change range.

[0221] The prediction module 1503 is used to perform dam body deformation prediction according to the initial formation change curve and the finite element algorithm to obtain the predicted dam body deformation dataset, and calculate the loss between the predicted dam body deformation dataset and the real dam body deformation dataset until the loss meets the preset loss condition.

[0222] The second determination module 1504 is used to determine the initial formation change curve corresponding to the predicted dam body deformation dataset in the current round as the formation change curve.

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

[0224] The first determination sub-module is configured to determine the influence area in the barrier lake dam model based on a set of material parameters, a 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 barrier lake dam.

[0225] The second determination sub-module is configured to determine the change range of the formation boundary in the influence area according to a borehole data set, a set of material parameters, a preset second difference threshold, and a finite element algorithm; the second difference threshold is the stress change threshold of the anti-seepage structure in the barrier lake dam.

[0226] In an exemplary embodiment, the barrier lake dam includes an anti-seepage structure, and the first determination sub-module includes:

[0227] The first division sub-module is configured to divide the barrier lake dam model into an initial influence area and an initial non-influence area.

[0228] The third determination sub-module is configured to determine the first stress of the anti-seepage structure according to a set of material parameters, the initial influence area, the initial non-influence area, and a finite element algorithm.

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

[0230] The fourth determination sub-module is configured to determine the influence area in the barrier lake dam model based on the first stress, the second stress, and a preset first difference threshold.

[0231] In an exemplary embodiment, the fourth determination sub-module is specifically configured to: the first processing sub-module is configured to perform a difference processing on the second stress and the first stress to obtain a first difference, and determine whether the first difference exceeds a preset first difference threshold. The fifth determination sub-module is configured to, if the first difference exceeds the first difference threshold, determine the adjusted initial influence area as the influence area; the first update sub-module is configured to, if the first difference does not exceed the first difference threshold, update the first stress according to the second stress; the first execution sub-module is configured to perform the steps of adjusting the initial influence area and the initial non-influence area according to the adjusted initial influence area and the adjusted initial non-influence area until the first difference exceeds the first difference threshold, and determine the adjusted initial influence area as the influence area.

[0232] In an exemplary embodiment, the second determination sub-module includes:

[0233] The sixth determination sub-module is configured to determine the initial change range of the formation boundary in the influence area, and determine the formation demarcation point in the initial change range of the formation boundary according to the borehole data set.

[0234] The first generation sub-module is used to generate an initial formation demarcation curve based on the formation demarcation point and the initial formation boundary change range.

[0235] The seventh determination sub-module is used to determine the formation boundary change range within the initial formation boundary change range according to the initial formation demarcation 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 determination sub-module includes:

[0237] The eighth determination sub-module is used to determine the third stress of the anti-seepage structure in the barrier dam according to the initial formation demarcation curve, the influence area, the material parameter set, and the finite element algorithm.

[0238] The second adjustment sub-module is used to adjust the initial formation demarcation curve towards the formation demarcation point, and determine the fourth stress of the anti-seepage structure according to the adjusted initial formation demarcation curve, the influence area, the material parameter set, and the finite element algorithm.

[0239] The ninth determination sub-module is used to determine the formation boundary change range within the initial formation boundary change range based on the third stress, the fourth stress, the adjusted initial formation demarcation curve, and the preset second difference threshold.

[0240] In an exemplary embodiment, the ninth determination sub-module specifically includes: The second processing sub-module is used to perform a difference processing on the third stress and the fourth stress to obtain a second difference, and determine whether the second difference exceeds the preset second difference threshold; the tenth determination sub-module is used to, if the second difference exceeds the second difference threshold, determine the formation boundary change range within the initial formation boundary change range based on the adjusted initial formation demarcation curve; the second update sub-module is used to, if the second difference does not exceed the second difference threshold, update the third stress based on the fourth stress; the second execution sub-module is used to execute the step of adjusting the initial formation demarcation curve towards the formation demarcation point according to the adjusted initial formation demarcation curve until the second difference exceeds the second difference threshold, and determine the formation boundary change range within the initial formation boundary change range based on the adjusted initial formation demarcation 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 sub-module is used to optimize the initial formation change curve according to the various control parameters, the finite element algorithm, and the optimization algorithm to obtain an 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.

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

[0244] The eleventh determination sub-module is used to determine that the loss meets the preset loss condition if the loss does not exceed the loss threshold.

[0245] The twelfth determination sub-module is used to, if the loss exceeds the loss threshold, perform the step of optimizing the initial formation change curve according to each control parameter, 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.

[0246] In an exemplary embodiment, the first optimization sub-module includes a second optimization sub-module and a first prediction sub-module. Among them, the first optimization sub-module is specifically used for the perturbation sub-module to determine the perturbation values of each control parameter, and perturb each control parameter according to the perturbation values of each control parameter to generate a perturbed initial formation change curve; the second prediction sub-module is used to perform dam deformation prediction based on the finite element algorithm and the perturbed initial formation change curve to obtain a perturbed dam deformation data set; the second calculation sub-module is used to calculate the perturbation error between the perturbed dam deformation data set and the real dam deformation data set; the second generation sub-module is used to calculate the perturbation gradient according to the perturbation error, and optimize each control parameter according to the perturbation gradient and the optimization algorithm to generate an optimized initial formation change curve.

[0247] Each module in the above-mentioned post-construction reconstruction device of the barrier dam based on physical mechanisms and survey and monitoring data can be implemented in whole or in part by software, hardware, and their combination. The above-mentioned modules can be embedded in the processor of the computer device in the form of hardware or independent of it, or stored in the memory of the computer device in the form of software, so that the processor can call and execute the operations corresponding to the above-mentioned modules.

[0248] In an exemplary embodiment, a computer device is provided. The computer device can be a terminal, and its internal structure diagram can be as Figure 16As shown in the figure. The computer device includes a processor, a memory, an input / output interface, a communication interface, a display unit, and an input device. Among them, the processor, the memory, and the input / output interface are connected through a system bus, and the communication interface, the display unit, and the input device are connected to the system bus through the input / output interface. Among them, 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 computer programs. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The input / output interface of the computer device is used to exchange information between the processor and external devices. The communication interface of the computer device is used to communicate with external terminals in a wired or wireless manner, and the wireless manner can be implemented through WIFI, a mobile cellular network, near field communication (NFC), or other technologies. When the computer program is executed by the processor, it implements a method for predicting barrier lakes.

[0249] Those skilled in the art can understand that Figure 16 the structure shown in the figure is only a block diagram of some structures 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 those shown in the figure, or combine some components, or have different component arrangements.

[0250] In one embodiment, a computer device is further provided, including a memory and a processor. A computer program is stored in the memory, and when the processor executes the computer program, the steps in the above method embodiments are implemented.

[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 the processor, the steps in the above method embodiments are implemented.

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

[0253] The technical features of the above embodiments can be combined arbitrarily. For the sake of brevity of description, all possible combinations of the technical features in the above embodiments are not described. However, as long as there is no contradiction in the combination of these technical features, it should be considered to be within the scope recorded in the present application.

[0254] The above-described embodiments merely represent several implementation manners of the present application. The description thereof is relatively specific and detailed, but it should not be construed as a limitation to the patent scope of the present application. It should be noted that for those of ordinary skill in the art, without departing from the concept of the present application, several modifications and improvements can still be made, and these all fall within the protection scope of the present application. Therefore, the protection scope of the present application shall be subject to the appended claims.

Claims

1. A post-construction reconstruction method for barrier dams based on physical mechanisms and survey and monitoring data, characterized in that The method includes: establishing a barrier dam model of the barrier dam, and obtaining a set of material parameters, a set of borehole data, and a set of real dam body deformation data of the barrier dam; determining a range of formation boundary changes in the barrier dam model based on the set of material parameters, the set of borehole data, and the finite element algorithm, and determining an initial formation change curve within the range of formation boundary changes; performing dam body deformation prediction according to the initial formation change curve and the finite element algorithm to obtain a predicted dam body deformation data set, and calculating a 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; determining the initial formation change curve corresponding to the predicted dam body deformation data set in the current round as the formation change curve.

2. The method according to claim 1, wherein The determining the range of formation boundary changes in the barrier dam model based on the set of material parameters, the set of borehole data, and the finite element algorithm includes: determining an influence area in the barrier dam model based on the set of material parameters, the finite element algorithm, and a preset first difference threshold; the first difference threshold is a stress change threshold of the anti-seepage structure in the barrier dam; determining the range of formation boundary changes in the influence area according to the set of borehole data, the set of material parameters, 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 barrier dam.

3. The method according to claim 2, wherein The barrier dam includes an anti-seepage structure, and the determining an influence area in the barrier dam model based on the set of material parameters, the finite element algorithm, and a preset first difference threshold includes: dividing the barrier dam model into an initial influence area and an initial non-influence area; determining a first stress of the anti-seepage structure according to the set of material parameters, the initial influence area, the initial non-influence area, and the finite element algorithm; adjusting the initial influence area and the initial non-influence area, and determining a second stress of the anti-seepage structure according to the set of material parameters, the adjusted initial influence area, the adjusted initial non-influence area, and the finite element algorithm; determining an influence area in the barrier dam model based on the first stress, the second stress, and a preset first difference threshold.

4. The method according to claim 3, wherein The determining an influence area in the barrier dam model based on the first stress, the second stress, and a preset first difference threshold includes: performing a difference process 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 influence area as the influence area; if the first difference does not exceed the first difference threshold, updating the first stress according to the second stress; performing the step of adjusting the initial influence area and the initial non-influence area according to the adjusted initial influence area and the adjusted initial non-influence area until the first difference exceeds the first difference threshold, and determining the adjusted initial influence area as the influence area.

5. The method according to claim 2, characterized in that Determining the range of formation boundary changes in the influence area according to the borehole dataset, the material parameter set, a preset second difference threshold, and the finite element algorithm includes: Determining the initial range of formation boundary changes in the influence area, and determining formation demarcation points in the initial range of formation boundary changes according to the borehole dataset; Generating an initial formation demarcation curve based on the formation demarcation points and the initial range of formation boundary changes; Determining the range of formation boundary changes in the initial range of formation boundary changes according to the initial formation demarcation curve, the influence area, the material parameter set, the preset second difference threshold, and the finite element algorithm.

6. The method according to claim 5, wherein The determining the range of formation boundary changes in the initial range of formation boundary changes according to the initial formation demarcation curve, the influence area, the material parameter set, the preset second difference threshold, and the finite element algorithm includes: Determining the third stress of the anti-seepage structure in the barrier dam according to the initial formation demarcation curve, the influence area, the material parameter set, and the finite element algorithm; Adjusting the initial formation demarcation curve in the direction of the formation demarcation points, and determining the fourth stress of the anti-seepage structure according to the adjusted initial formation demarcation curve, the influence area, the material parameter set, and the finite element algorithm; Determining the range of formation boundary changes in the initial range of formation boundary changes based on the third stress, the fourth stress, the adjusted initial formation demarcation curve, and the preset second difference threshold.

7. The method according to claim 6, wherein The determining the range of formation boundary changes in the initial range of formation boundary changes based on the third stress, the fourth stress, the adjusted initial formation demarcation curve, and the preset second difference threshold includes: Performing a difference process on the third stress and the fourth stress to obtain a second difference, and determining whether the second difference exceeds the preset second difference threshold; If the second difference exceeds the second difference threshold, determining the range of formation boundary changes in the initial range of formation boundary changes based on the adjusted initial formation demarcation curve; If the second difference does not exceed the second difference threshold, updating the third stress based on the fourth stress; Executing the step of adjusting the initial formation demarcation curve in the direction of the formation demarcation points according to the adjusted initial formation demarcation curve until the second difference exceeds the second difference threshold, and determining the range of formation boundary changes in the initial range of formation boundary changes based on the adjusted initial formation demarcation curve.

8. The method according to claim 1, characterized in that, The initial formation change curve includes various control parameters. According to the initial formation change curve and the finite element algorithm, predicting the dam body deformation to obtain a predicted dam body deformation dataset, and calculating the loss between the predicted dam body deformation dataset and the real dam body deformation dataset until the loss meets the preset loss condition, including: 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 based on 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 real 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.

9. The method according to claim 8, characterized in that, The step of optimizing 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 includes: Determine the perturbation values of each of the control parameters, and perturb each of the control parameters according to the perturbation values of each of the control parameters to generate a perturbed initial formation change curve; Based on the finite element algorithm and the perturbed initial formation change curve, perform dam body deformation prediction to obtain a perturbed dam body deformation data set; Calculate the perturbation error between the perturbed dam body deformation data set and the real dam body deformation data set; Calculate the perturbation gradient according to the perturbation error, and optimize each of the control parameters according to the perturbation gradient and the optimization algorithm to generate an optimized initial formation change curve.

10. A post-construction reconstruction device for barrier dams based on physical mechanisms and survey and monitoring data, characterized in that, The device includes: An acquisition module, configured to establish a barrier dam model of the barrier dam, and acquire a material parameter set, a borehole data set, and a real dam body deformation data set of the barrier dam; A first determination module, configured to determine a formation boundary change range in the barrier dam model based on the material parameter set, the borehole data set, and the finite element algorithm, and determine an initial formation change curve within the formation boundary change range; A prediction module, configured to perform dam body deformation prediction according to the initial formation change curve and the finite element algorithm to obtain a predicted dam body deformation data set, and calculate the loss between the predicted dam body deformation data set and the real dam body deformation data set until the loss meets the preset loss condition; A second determination module, configured to determine the initial formation change curve corresponding to the predicted dam body deformation data set in the current round as the formation change curve.

Images