A dam leakage electrical method and hydrological coupling modeling method, device and medium

Abstract

The application relates to the technical field of dam seepage, in particular to a dam seepage electro-hydrological coupling modeling method, equipment and medium, the method comprises the following steps: obtaining hydrological parameters of a dam research area; constructing an electro-hydrological coupling relationship to realize dynamic calculation of conductivity; mapping irregular spatial data to a calculation grid, and generating medium conductivity of the calculation grid by using interpolation calculation; performing electric field forward calculation to obtain an electric potential field; calculating a potential difference according to an electrode layout mode to obtain an electro-response; analyzing a change trend based on the electro-response, and establishing a relationship between hydrological changes and the electro-response to realize electro-hydrological coupling. Compared with a traditional method, the method can automatically update a resistivity or conductivity field from a hydrological process, reduce uncertainty caused by artificial setting of a resistivity anomaly, and improve consistency between an electro-explanation result and an actual seepage process.

Description

Technical Field

[0001] This invention relates to the field of dam seepage technology, and in particular to an electrical and hydrological coupled modeling method, equipment and medium for dam seepage. Background Technology

[0002] Existing electrical resistivity forward modeling methods rely solely on static resistivity models, which fail to reflect the impact of dynamic changes in groundwater on the electric field response. This is particularly problematic in engineering scenarios such as dam leakage monitoring, where it is difficult to predict changes in electrical resistivity response based on the evolution of hydrological conditions, thus affecting the accuracy of monitoring data interpretation and the reliability of engineering decisions.

[0003] In the field of dam leakage monitoring, resistivity method (ERT, Electrical Resistivity Tomography) is a commonly used non-destructive testing method. Its basic principle is to lay electrodes on the ground surface, inject current and measure the potential difference (ΔV) to invert the underground resistivity distribution.

[0004] Current mainstream methods include: Static resistivity forward modeling method: based on known conductivity distribution The Laplace equation is solved using the finite element or finite difference method to obtain the potential distribution U, and then the potential difference ΔV between the electrodes is calculated.

[0005] Conventional resistivity inversion method: Resistivity distribution is inverted by minimizing the difference between observed data and predicted data.

[0006] Archie empirical relation (local use): Some studies use empirical formulas to correlate water saturation with electrical conductivity, but these are usually only used to interpret the results and are not included in the forward modeling process.

[0007] The disadvantages of existing methods include: The conductivity model is artificially set (non-dynamic), thus ignoring the time evolution of the hydrological process. This leads to the independence of the electrical method and the hydrological model, which affects the accuracy of the interpretation of monitoring data and the reliability of engineering decisions.

[0008] In summary, there is an urgent need for a method, equipment, and medium for coupled electrical and hydrological modeling. Summary of the Invention

[0009] To address the aforementioned technical problems, this invention proposes a coupled electrical and hydrological modeling method, equipment, and medium for dam seepage.

[0010] In a first aspect, the present invention provides a method for coupled electrical and hydrological modeling, comprising the following steps: S1. Obtain hydrological parameters of the dam study area; S2. Establish the coupling relationship between electrical methods and hydrology to realize dynamic calculation of electrical conductivity; S3. Map irregular spatial data to a computational grid and use interpolation to calculate the dielectric conductivity of the generated computational grid. S4. Perform forward electric field calculations to obtain the electric potential field. S5. Calculate the potential difference based on the electrode arrangement to obtain the electrical response; S6. Based on the electrical response analysis of the changing trend, establish the relationship between hydrological changes and electrical response, and realize the coupling of electrical method and hydrology.

[0011] Optionally, in S1, the hydrological parameters of the dam study area include water saturation. Porosity Permeability coefficient and pressure head ,in, The horizontal coordinates of the two-dimensional calculation profile. Vertical coordinates or depth coordinates For calculating time or monitoring time, express Calculate the position within the domain at any time The state at which it is located.

[0012] Optionally, in S2, the expression for the coupling relationship between electrical methods and hydrology is as follows: ; in, Indicates the conductivity of the dielectric. Represents empirical parameters. This represents the resistivity of pore water.

[0013] Optionally, in S3, irregular spatial data is mapped to a computational grid, and interpolation is used to calculate the dielectric conductivity of the generated computational grid. The process includes: Irregular spatial data is mapped to a computational grid, and an interpolation method is used to generate... Calculate the position within the domain at any time The dielectric conductivity at a given location is used to generate the corresponding dielectric resistivity, wherein the dielectric conductivity and the dielectric resistivity are reciprocals of each other.

[0014] Optionally, in S4, the electric field forward modeling is performed using the finite element method, and the process includes: Construct the element stiffness matrix and assemble the global stiffness matrix. K Applying current source boundary conditions, we solve for the electric potential field, and the calculation expression is: ,in, U Represents the electric potential field. BThis represents the global current source vector assembled from the current source terms of the power supply electrode and the return electrode, as well as the boundary conditions.

[0015] Optionally, in S5, the electrical response... The calculation expression is as follows: ; in, , These represent the electric potentials at positions M and N, respectively, under the action of the power supply electrode pair consisting of power supply electrode A and power supply electrode B.

[0016] Optionally, in S6, the trend of change is analyzed based on electrical response, and the relationship between hydrological changes and electrical response is established. The process includes: Obtain water saturation at different times ; When the local water content or water saturation in the dam study area increases, according to the Archie relation, the resistivity of the medium in the dam study area decreases, and the corresponding conductivity of the medium increases. Dielectric conductivity Using this as input for the forward modeling of the electric field, the governing equations of the electric field are obtained by solving: ; in, For electric potential field, For current source terms, Represents a spatial differential operator; when it acts on an electric potential field hour, Represents the potential gradient; For the electrode combination ABMN, the potential values ​​at measuring electrodes M and N are extracted at each time step, and the potential difference is calculated: ; in, This represents the electrical response sequence of the electrode combination ABMN at different times. Indicates at time When current is applied to power supply electrodes A and B, the potential at electrode M is measured. Indicates at time The potential at electrode N is measured when current is applied to power supply electrodes A and B.

[0017] The relationship between hydrological changes and electrical responses is constructed based on the aforementioned electrical response sequence, as follows: Changes in water saturation lead to changes in the conductivity and resistivity of the medium, which in turn lead to changes in the electric potential field, ultimately resulting in changes in the electrical response. By comparing the changes in electrical response at different times, the development trends of high water content areas, low resistivity anomaly areas, and potential seepage channels within the dam study area are obtained.

[0018] In a second aspect, the present invention also provides a computer device, including a memory and a processor; The memory stores computer programs; When the processor executes the computer program, it implements the steps of the above-described electrical and hydrological coupled modeling method.

[0019] In a third aspect, the present invention also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-described electrical and hydrological coupled modeling method.

[0020] The application of the technical solution of the present invention has at least the following beneficial effects: This invention provides a method, equipment, and medium for coupled electrical resistivity and hydrological modeling. The method establishes a dynamic coupling model between water saturation and conductivity, incorporating the hydrological time evolution process into electrical resistivity forward modeling to construct a time-varying electrical resistivity response, thus achieving unified modeling of hydrological processes and electric field responses. The model constructed by this invention uses hydrological conditions to drive conductivity, and the constructed time-varying electrical resistivity response enhances time resolution. Furthermore, this method directly correlates changes in electrical resistivity observations with hydrological information, allowing for direct application in engineering to predict leakage trends and conduct risk assessments. Attached Figure Description

[0021] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the structures shown in these drawings without creative effort.

[0022] Figure 1 This is a flowchart illustrating the steps of a preferred embodiment of the present invention for a coupled electrical and hydrological modeling method. Figure 2 This is a schematic diagram of the mesh division of the dam seepage model in a preferred embodiment of the present invention; Figure 3 This is a graph showing the resistivity forward modeling results for the second day of the numerical simulation in a preferred embodiment of the present invention. Figure 4 This is a graph showing the saturation forward modeling results for the second day of the numerical simulation in a preferred embodiment of the present invention. Figure 5 This is a day-by-day resistivity diagram of the electric field forward modeling in a preferred embodiment of the present invention. Detailed Implementation

[0023] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of the present invention.

[0024] Example 1: like Figure 1 As shown, this embodiment provides a coupled electrical and hydrological modeling method, including the following steps (S1 to S6): S1. Obtain hydrological parameters of the dam study area.

[0025] In this embodiment, the hydrological parameters of the dam study area include water saturation. Porosity Permeability coefficient and pressure head ,in, The horizontal coordinates of the two-dimensional calculation profile. Vertical coordinates or depth coordinates For calculating time or monitoring time, express Calculate the position within the domain at any time The state at which it is located.

[0026] S2. Construct a coupling relationship between electrical methods and hydrology to achieve dynamic calculation of electrical conductivity.

[0027] In this embodiment, the expression for the coupling relationship between electrical methods and hydrology is as follows: ; in, Indicates the conductivity of the dielectric. Represents empirical parameters. This represents the resistivity of pore water.

[0028] S3. Map irregular spatial data to a computational grid and use interpolation to calculate the dielectric conductivity of the generated computational grid.

[0029] Mapping irregular spatial data to a computational grid and using interpolation to calculate the dielectric conductivity of the computational grid involves the following steps: Irregular spatial data is mapped to a computational grid, and an interpolation method is used to generate... Calculate the position within the domain at any time ( x,z The dielectric conductivity at a given point is used to generate the corresponding dielectric resistivity, wherein the dielectric conductivity and the dielectric resistivity are reciprocals of each other.

[0030] S4. Perform forward electric field calculations to obtain the electric potential field.

[0031] The electric field forward modeling calculation is performed using the finite element method, and the process includes: Construct the element stiffness matrix and assemble the global stiffness matrix. K Applying current source boundary conditions, we solve for the electric potential field, and the calculation expression is: ,in, U Represents the electric potential field. B This represents the global current source vector assembled from the current source terms of the power supply electrode and the return electrode, as well as the boundary conditions.

[0032] Optionally, in S5, the electrical response... The calculation expression is as follows: ; in, , These represent the electric potentials at positions M and N, respectively, under the action of the power supply electrode pair consisting of power supply electrode A and power supply electrode B.

[0033] S5. Calculate the potential difference based on the electrode arrangement to obtain the electrical response.

[0034] S6. Based on the electrical response analysis of the changing trend, establish the relationship between hydrological changes and electrical response, and realize the coupling of electrical method and hydrology.

[0035] In this embodiment, the trend of change is analyzed based on electrical response, and the relationship between hydrological changes and electrical response is established. The process includes: Obtain water saturation at different times ; When the local water content or water saturation in the dam study area increases, according to the Archie relation, the resistivity of the medium in the dam study area decreases, and the corresponding conductivity of the medium increases. Dielectric conductivity Using this as input for the forward modeling of the electric field, the governing equations of the electric field are obtained by solving: ; in, For electric potential field, For current source terms, Represents a spatial differential operator; when it acts on an electric potential field hour, Represents the potential gradient; For the electrode combination ABMN, the potential values ​​at measuring electrodes M and N are extracted at each time step, and the potential difference is calculated: ; in, This represents the electrical response sequence of the electrode combination ABMN at different times. Indicates at time When current is applied to power supply electrodes A and B, the potential at electrode M is measured. Indicates at time The potential at electrode N is measured when current is applied to power supply electrodes A and B.

[0036] The relationship between hydrological changes and electrical responses is constructed based on the aforementioned electrical response sequence, as follows: Changes in water saturation lead to changes in the conductivity and resistivity of the medium, which in turn lead to changes in the electric potential field, ultimately resulting in changes in the electrical response. By comparing the changes in electrical response at different times, the development trends of high water content areas, low resistivity anomaly areas, and potential seepage channels within the dam study area are obtained.

[0037] To verify the beneficial effects of the method in this embodiment, a simulation experiment is conducted below: Constructing a dam leakage model, such as Figure 2 As shown, the dam body is an isosceles trapezoid with a base width of 52 m, a top width of 24 m, and a vertical height of 12 m. The dam foundation is rectangular with a base width of 72 m and a vertical height of 10 m. Based on the actual reservoir operating environment, the model boundary conditions are set as follows: 1. Upstream boundary: Apply a first-type boundary condition (Dirichlet Boundary) to the upstream slope and upstream dam foundation surface, setting the upstream total head to H = 5.0 m to simulate the upstream reservoir water level; 2. Downstream boundary: Apply a constant head boundary at the downstream far boundary and downstream slope toe, setting the downstream total head to H = 2.0 m to simulate the downstream tailwater; 3. Escape surface boundary: Apply a mixed boundary condition (System Flux / Potential Seepage Face) in the potential seepage area on the downstream slope, allowing groundwater to escape naturally under gravity; 4. Impermeable boundary: The bottom of the model and the lateral boundaries without specified head / flow are set to zero flux boundaries by default.

[0038] Using the above-mentioned dam seepage model as the dam research area, and implementing the method of this embodiment, the change in water saturation can be converted into a change in resistivity through the Archie relation. Therefore, the formation and development of low resistivity regions in the simulation results can correspond to the changes in the water content state inside the dam body, reflecting the physical coupling relationship between the hydrological field and the electric field.

[0039] like Figure 3 and Figure 4 The simulation results shown indicate that the forward modeling data for the second day shows that changes in water saturation correspond to changes in resistivity.

[0040] like Figure 5 The resistivity map shown in the forward modeling diagram illustrates the evolution of low-resistivity anomaly regions over time. This demonstrates that the method in this embodiment does not rely solely on static resistivity modeling but rather describes the time-varying electrical response during leakage development. This allows for the determination of whether leakage anomalies are intensifying, migrating, or expanding, thus showcasing trend analysis and monitoring / early warning capabilities. Compared to traditional static electrical resistivity forward modeling that only provides a fixed resistivity field, the method in this embodiment can automatically update the resistivity or conductivity field based on hydrological processes, reducing the uncertainty caused by artificially set resistivity anomalies and improving the consistency between electrical interpretation results and actual seepage processes.

[0041] In addition, this embodiment also provides a computer device, including a memory and a processor; The memory stores computer programs; When the processor executes the computer program, it implements the steps of the above-described electrical and hydrological coupled modeling method.

[0042] It should be noted that computer program code for performing the operations of this application can be written in one or more programming languages ​​or a combination thereof. These programming languages ​​include object-oriented programming languages—such as Java, Smalltalk, and C++—as well as conventional procedural programming languages—such as the "C" language or similar programming languages. The program code can be executed entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving remote computers, the remote computer can be connected to the user's computer via any type of network—including a Local Area Network (LAN) or a Wide Area Network (WAN)—or can be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0043] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of this application. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. For example, two consecutively indicated blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in the block diagrams and / or flowcharts, and combinations of blocks in the block diagrams and / or flowcharts, can be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.

[0044] In addition, this embodiment also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the above-described electrical and hydrological coupled modeling method.

[0045] The readable storage medium provided in this application is a computer-readable storage medium, which stores computer-readable program instructions (i.e., a computer program) for executing the above-described electrical and hydrological coupled modeling method. Compared with the prior art, the beneficial effects of the computer-readable storage medium provided in this application are the same as the beneficial effects of the electrical and hydrological coupled modeling method provided in the above embodiments, and will not be repeated here.

[0046] The above description is only a preferred embodiment of the present invention and does not limit the scope of the present invention. All equivalent structural transformations made under the inventive concept of the present invention using the contents of the present invention specification and drawings, or direct / indirect applications in other related technical fields, are included within the protection scope of the present invention.

Claims

1. A coupled electrical and hydrological modeling method for dam seepage, characterized in that, Includes the following steps: S1. Obtain hydrological parameters of the dam study area; S2. Establish the coupling relationship between electrical methods and hydrology to realize dynamic calculation of electrical conductivity; S3. Map irregular spatial data to a computational grid and use interpolation to calculate the dielectric conductivity of the generated computational grid. S4. Perform forward electric field calculations to obtain the electric potential field. S5. Calculate the potential difference based on the electrode arrangement to obtain the electrical response; S6. Based on the electrical response analysis of the changing trend, establish the relationship between hydrological changes and electrical response, and realize the coupling of electrical method and hydrology.

2. The electrical resistivity and hydrological coupled modeling method according to claim 1, characterized in that, In S1, the hydrological parameters of the dam study area include water saturation. Porosity Permeability coefficient and pressure head ,in, The horizontal coordinates of the two-dimensional calculation profile. Vertical coordinates or depth coordinates For calculating time or monitoring time, express Calculate the position within the domain at any time The state at which it is located.

3. The electrical resistivity and hydrological coupled modeling method according to claim 2, characterized in that, In S2, the expression for the coupling relationship between electrical methods and hydrology is as follows: ； in, Indicates the conductivity of the dielectric. Represents empirical parameters. This represents the resistivity of pore water.

4. The electrical resistivity and hydrological coupled modeling method according to claim 3, characterized in that, In S3, irregular spatial data is mapped to a computational grid, and interpolation is used to calculate the dielectric conductivity of the generated computational grid. The process includes: Irregular spatial data is mapped to a computational grid, and an interpolation method is used to generate... Calculate the position within the domain at any time The dielectric conductivity at a given location is used to generate the corresponding dielectric resistivity, wherein the dielectric conductivity and the dielectric resistivity are reciprocals of each other.

5. The electrical resistivity and hydrological coupled modeling method according to claim 4, characterized in that, In S4, the electric field forward modeling is performed using the finite element method, and the process includes: Construct the element stiffness matrix and assemble the global stiffness matrix. K Applying current source boundary conditions, we solve for the electric potential field, and the calculation expression is: ,in, U Represents the electric potential field. B This represents the global current source vector assembled from the current source terms of the power supply electrode and the return electrode, as well as the boundary conditions.

6. The electrical resistivity and hydrological coupled modeling method according to claim 5, characterized in that, In S5, electrical response The calculation expression is as follows: ； in, , These represent the electric potentials at positions M and N, respectively, under the action of the power supply electrode pair consisting of power supply electrode A and power supply electrode B.

7. The electrical resistivity and hydrological coupled modeling method according to claim 6, characterized in that, In S6, the trend of change is analyzed based on electrical response, and the relationship between hydrological changes and electrical response is established. The process includes: Obtain water saturation at different times ; When the local water content or water saturation in the dam study area increases, according to the Archie relation, the resistivity of the medium in the dam study area decreases, and the corresponding conductivity of the medium increases. Dielectric conductivity Using this as input for the forward modeling of the electric field, the governing equations of the electric field are obtained by solving: ； in, For electric potential field, For current source terms, Represents a spatial differential operator; when it acts on an electric potential field hour, Represents the potential gradient; For the electrode combination ABMN, the potential values ​​at measuring electrodes M and N are extracted at each time step, and the potential difference is calculated: ； in, This represents the electrical response sequence of the electrode combination ABMN at different times. Indicates at time When current is applied to power supply electrodes A and B, the potential at electrode M is measured. Indicates at time The potential at electrode N is measured when current is applied to power supply electrodes A and B. The relationship between hydrological changes and electrical responses is constructed based on the aforementioned electrical response sequence, as follows: Changes in water saturation lead to changes in the conductivity and resistivity of the medium, which in turn lead to changes in the electric potential field, ultimately resulting in changes in the electrical response. By comparing the changes in electrical response at different times, the development trends of high water content areas, low resistivity anomaly areas, and potential seepage channels within the dam study area are obtained.

8. A computer device, characterized in that, Including memory and processor; The memory stores computer programs; When the processor executes the computer program, it implements the steps of the electrical and hydrological coupled modeling method as described in any one of claims 1-7.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the program is executed by the processor, it implements the electrical and hydrological coupled modeling method as described in any one of claims 1-7.

Images