CN-121997430-A - Impermeable simulation method and system for dam geomembrane

CN121997430ACN 121997430 ACN121997430 ACN 121997430ACN-121997430-A

Abstract

The invention discloses an anti-seepage simulation method and system for a geomembrane of a dam body, and relates to the technical field of hydraulic engineering, wherein the method comprises the steps of integrating three-dimensional geometric data, filling layer arrangement information, a geomembrane laying form, a membrane-soil interface contact structure and water content distribution of the dam body to generate a multi-source structural feature domain; the method comprises the steps of combining a stress memory field, identifying potential seepage weak areas, dynamically dividing multi-level risk sub-areas of a dam body through intelligent disturbance signals to excite responses of seepage, stress and crack energy fields, constructing a coupling evolution field, calculating a minimum energy seepage path by combining an energy function, and constructing an anti-seepage risk field of the dam body. The method solves the technical problems that the existing dam body seepage prevention method cannot accurately identify potential seepage weak areas and evaluate seepage risks, and the seepage prevention design accuracy and reliability are insufficient, and achieves the technical effects of accurately identifying the seepage weak areas and improving the dam body seepage prevention design accuracy and reliability through comprehensive simulation and dynamic risk evaluation.

Inventors

LIU LINDONG
XIAO HUACHUN
ZHANG JINGYING

Assignees

山东天海新材料工程有限公司

Dates

Publication Date: 20260508
Application Date: 20260128

Claims (10)

1. A method for simulating seepage control of a geomembrane of a dam, the method comprising: Generating a multisource construction characteristic domain of the dam body through three-dimensional geometric data of the dam body, filling layer physical information, geomembrane laying forms, a membrane-soil interface contact structure and dam body water content distribution; after the residual stretching, the fold shape and the anchoring state of the geomembrane are called, a geomembrane stress memory field of local initial stress non-uniformity in the reaction membrane is established according to the multi-source construction characteristic domain; After extracting geometric variation gradient characteristics of the dam body, calculating and constructing a potential seepage weak area according to the geomembrane stress memory field; Injecting intelligent disturbance signals into the potential seepage weak areas to trigger initial responses of seepage fields, stress fields and crack energy fields; Dynamically dividing a multi-level risk subarea of a dam body according to the amplification tracks of the initial response in a seepage pressure field, a geomembrane stress field and an interface energy field, and forming a coupling evolution field amplified along with disturbance according to the multi-level risk subarea; And constructing a joint energy function of seepage energy, in-film strain energy and interfacial crack propagation energy based on the coupling evolution field, solving a minimum energy permeation path required by leakage damage of the geomembrane, and establishing a dam seepage prevention risk field.
2. The method for simulating seepage prevention of a geomembrane of a dam body according to claim 1, wherein after extracting gradient characteristics of geometric variation of the dam body, constructing potential seepage weak areas according to stress memory field calculation of the geomembrane comprises the following steps: After identifying the positions of the seam lines according to the multi-source structural feature domain, executing seam line segmentation processing; Taking each joint line segment as a unit, counting the rigidity and compactness difference of soil bodies at two sides segment by segment, and establishing a material non-average mark according to a counting result; Performing stress concentration association analysis on each joint line segment according to the geomembrane stress memory field, and establishing an association stress abnormal value mark; Performing joint weak area identification of the geomembrane by using the material non-average value identification and the associated stress abnormal value identification, and establishing a potential joint weak area; and constructing a potential penetration weak area according to the potential joint weak area.
3. The method of simulating seepage control of a geomembrane for a dam of claim 2, wherein constructing a potential permeable zone from the potential joint zone of weakness further comprises: Calculating curvature values point by point on the surface of the geomembrane, and constructing a neighborhood average value of the curvature; marking curvature mutation points according to the curvature value and the neighborhood average value; Tracking a continuous distribution line of residual strain along the trend of the fold by utilizing the curvature abrupt change points, and if the residual strain continuously meets a preset threshold along a linear or curve direction and crosses the curvature abrupt change points, constructing a main axis of the fold; Performing area search along the main axis of the folds, detecting local discontinuous positions of axial displacement and thickness change of the folds, and constructing a fold concentration area according to detection results; The fold concentration areas are graded according to fold superposition degree and axial strain accumulation along the major axis of the folds, and potential fold concentration weak areas are established; and constructing a potential penetration weak area according to the potential joint weak area and the potential fold concentration weak area.
4. The method for simulating seepage control of a geomembrane of a dam body according to claim 3, wherein constructing a potential seepage weakness according to the potential joint weakness and the potential fold concentration weakness comprises: Extracting film-soil interface roughness, local laminating degree and moisture distribution characteristics in a multisource structural characteristic domain; Establishing a potential interface void weak zone according to the extraction result and combining sudden drop identification and reverse distribution identification of the geomembrane stress memory field; Identifying a potential high stress concentration weak area by utilizing stress gradient peak points of a geomembrane stress memory field and combining material rigidity differences of dam body topography inflection positions and building layer abrupt change positions; And the potential joint weak area, the potential fold concentration weak area, the potential interface void weak area and the potential high stress concentration weak area construct a potential penetration weak area.
5. The method of simulating seepage control of a geomembrane of a dam of claim 1, wherein injecting intelligent perturbation signals into the potential permeable weak zone comprises: Generating a disturbance target parameter set of each weak zone based on the spatial boundary, the type label, the geomembrane stress memory field and the local material non-uniformity of the potential permeable weak zone; Performing disturbance signal type matching triggering seepage, stress and crack energy response by taking the disturbance target parameter set as a matching target, and establishing a matching result, wherein the matching result comprises local pressure gradient fluctuation, mechanical displacement, temperature disturbance or interface contact condition; And based on the initial stress peak value, the membrane-soil interface rigidity and the fold morphology of the potential osmotic weak zone, carrying out self-adaptive optimization on the amplitude, the acting direction and the time sequence on the matching result, and establishing an intelligent disturbance signal.
6. The method for simulating seepage control of a geomembrane of a dam of claim 1, wherein dynamically dividing the multi-level risk sub-zone of the dam comprises: Extracting a multidimensional response track data set in initial response, calculating seepage sensitivity, stress sensitivity and interface energy sensitivity along a dam space grid unit, and calculating local sensitivity peak value and gradient change according to a calculation result to form a local risk index; the dam body is divided into a plurality of levels of multi-level risk subareas according to the spatial distribution of the local risk indexes and the sensitivity gradient, wherein the multi-level risk subareas comprise a high risk subarea, a medium risk subarea and a low risk subarea.
7. The method for simulating seepage control of a geomembrane of a dam body according to claim 1, wherein the joint energy function is as follows: ; Wherein, the The combined energy is characterized by the fact that, Characterizing potential permeation pathways or dam unit volume areas, The osmotic pressure is characterized by the fact that, The modulus of local resistance is characterized by, The local stress of the geomembrane is represented, The elastic modulus of the geomembrane is characterized, The geomembrane unit volume is characterized by that, The critical energy release rate of the interfacial crack is characterized, The actual area of interfacial crack propagation is characterized, Is a weight coefficient.
8. The method of modeling the permeability resistance of a geomembrane of a dam of claim 7, wherein solving the minimum energy permeation path required to cause leakage damage to the geomembrane comprises: dispersing the coupling evolution place into a plurality of voxel units and interface units in a three-dimensional space in a dam body area, adaptively adjusting the unit density according to the spatial distribution of energy gradients in the coupling evolution place, and respectively calculating the energy contribution of the voxel units and the interface units by utilizing the combined energy function to form an accumulated unit energy label; Mapping the voxel center or interface unit nodes into vertexes of the graph, establishing edges of the graph by using space adjacency or unit nodes possibly communicated through cracks and joints, and configuring energy cost weights of the edges by using the accumulative unit energy labels to construct an energy weighted graph; Performing global shortest path search on a coarse grid scale based on the energy weighting graph, constructing N candidate paths, and recording total energy and path topological characteristics of each candidate path; And (3) performing grid encryption on the neighborhood of the N candidate paths, calculating joint energy for voxel units and interface units in the encrypted local area, updating corresponding side weights, re-performing path search on the encrypted grid, and constructing a minimum energy penetration path.
9. The method of modeling seepage control of a geomembrane of a dam of claim 1, wherein generating a multi-source structural feature of the dam comprises: And constructing a twin simulation model by utilizing the multi-source construction feature domain, executing multi-level risk sub-region based seepage prevention simulation based on the potential seepage weak region in the twin simulation model, and establishing a dam seepage prevention risk field.
10. An impermeable simulation system for a geomembrane of a dam, the system comprising: the multi-source construction feature domain construction module is used for generating a multi-source construction feature domain of the dam body through three-dimensional geometric data, filling layer physical information, geomembrane laying forms, a membrane-soil interface contact structure and dam body water content distribution of the dam body; The geomembrane stress memory field construction module is used for establishing a geomembrane stress memory field with local initial stress non-uniformity in the reaction membrane according to the multi-source construction characteristic field after the residual stretching, the fold morphology and the anchoring state of the geomembrane laying process are called; The potential seepage weak area identification module is used for calculating and constructing a potential seepage weak area according to the geomembrane stress memory field after extracting the geometric change gradient characteristics of the dam body; The initial response triggering module is used for injecting intelligent disturbance signals into the potential seepage weak area so as to trigger initial responses of seepage fields, stress fields and crack energy fields; The coupling evolution field generation module is used for dynamically dividing a multi-level risk subarea of the dam body according to the amplification tracks of the initial response in the seepage pressure field, the geomembrane stress field and the interface energy field, and forming a coupling evolution field amplified along with disturbance according to the multi-level risk subarea; the seepage-proofing risk field establishment module is used for establishing a joint energy function of seepage energy, in-film strain energy and interfacial crack propagation energy based on the coupling evolution field, solving a minimum energy seepage path required by seepage damage of the geomembrane, and establishing a dam seepage-proofing risk field.

Description

Impermeable simulation method and system for dam geomembrane Technical Field The invention relates to the technical field of hydraulic engineering, in particular to an anti-seepage simulation method and system for a dam geomembrane. Background The geomembrane is widely applied to dam body seepage prevention engineering and is used for preventing water seepage. However, as the service time increases, the geomembrane may be deformed or damaged due to uneven laying, wrinkles, poor joints, etc., thereby affecting the impermeable performance of the dam. The traditional dam body seepage prevention design generally depends on static analysis and simplifying assumptions, complex coupling effects of factors in the dam body are difficult to comprehensively consider, and particularly, certain limitations exist in aspects of identifying potential seepage weak areas and optimizing seepage prevention design in terms of seepage path and stress distribution, the potential seepage weak areas cannot be accurately identified and seepage risk can not be estimated, and seepage prevention design accuracy and reliability are insufficient. Disclosure of Invention The application provides an anti-seepage simulation method and system for a dam geomembrane, which are used for solving the technical problems that the existing dam anti-seepage method cannot accurately identify potential seepage weak areas and evaluate seepage risks, so that the anti-seepage design accuracy and reliability are insufficient. The application provides an impermeable simulation method of a geomembrane of a dam body, which comprises the steps of generating a multi-source structural feature field of the dam body through three-dimensional geometric data, filling layer physical information, a geomembrane laying form, a membrane-soil interface contact structure and dam body water content distribution of the dam body, after residual stretching, fold forms and anchoring states of the geomembrane laying process are called, establishing a geomembrane stress memory field with local initial stress non-uniformity in a reaction membrane according to the multi-source structural feature field, after geometrical change gradient features of the dam body are extracted, constructing a potential permeable weak region according to the geomembrane stress memory field, injecting intelligent disturbance signals into the potential permeable weak region to trigger initial responses of the permeable field, the stress field and a crack energy field, dynamically dividing a multi-level risk subarea of the dam body according to amplification tracks of the initial responses in the permeable pressure field, the geomembrane stress field and the interface energy field, constructing a combined permeable membrane evolution path based on the coupling field, solving the minimum permeable membrane expansion function, and enabling the minimum permeable membrane to be required to have an impermeable path. The application provides an impermeable simulation system of a dam geomembrane, which comprises a multisource construction feature domain construction module, an initial response triggering module, a coupling field generation module and an impermeable weak area joint analysis module, wherein the multisource construction feature domain construction module is used for generating multisource construction feature domains of the dam through three-dimensional geometric data of the dam, filling layer physical information, a geomembrane laying form, a membrane-soil interface contact structure and dam moisture content distribution, the geotechnical membrane stress memory domain construction module is used for calling residual stretching, fold forms and anchoring states of the geotechnical membrane laying process, then establishing a geomembrane stress memory field with local initial stress non-uniformity in a reaction membrane according to the multisource construction feature domains, the potential permeable weak area identification module is used for calculating and constructing a potential permeable weak area according to the geomembrane stress memory field after extracting gradient features of geometric change of the dam, the initial response triggering module is used for injecting intelligent disturbance signals into the potential permeable weak area so as to trigger the initial responses of a seepage field, a stress field and a crack energy field, the coupling field generation module is used for dynamically dividing the differential expansion tracks in the seepage field, the geotechnical membrane stress field and the interface energy field according to the initial responses, and the differential expansion risk expansion function of the differential expansion interface energy field is established according to the differential risk expansion function of the differential expansion interface stress field, and the differential risk expansion interface energy field is used for solving the minimum risk expans