A fully coupled simulation method and system for surface water and groundwater

Abstract

This invention discloses a fully coupled simulation method and system for surface water and groundwater. The method includes: acquiring the spatial distribution characteristics of surface water and groundwater, dividing the surface water and groundwater into separate grids, and defining the interface between them; constructing numerical models of the surface water and groundwater based on the grid division results and the interface; subsequently, using the constructed numerical models to numerically solve the hydrodynamic control equations and solute transport equations for their respective regions, calculating the water flux and solute flux at the interface, and synchronously transmitting the flux information to the control equations of the other region in real time, thus completing the fully coupled simulation calculation. This invention achieves real-time interaction and synchronous feedback between surface water and groundwater by setting bidirectional conserved flux boundary conditions at the interface, overcoming the coupling bias caused by inaccurate interface flux calculation in existing technologies, and improving the overall accuracy of the simulation calculation.

Description

Technical Field

[0001] This invention relates to the field of water conservancy engineering technology, and in particular to a fully coupled simulation method and system for surface water and groundwater. Background Technology

[0002] In ocean, river, and lake systems, surface water and groundwater are typically connected by irregular bed morphologies, and their coupling process constitutes a dynamically interconnected system. Along with water exchange, most solutes can migrate and transform between surface and groundwater systems through multiple pathways, increasing the complexity of surface and groundwater interactions and thus affecting water quality and aquatic ecosystem functions. Therefore, to deeply understand the mechanisms of water transport and solute transformation during surface and groundwater interactions, it is necessary to study them as a unified, coupled system. However, traditional numerical simulation methods for surface and groundwater coupling mostly employ unidirectional sequential coupling: that is, first calculating the flow field information of the surface water, using pressure as a boundary condition to drive the groundwater model, and then separately calculating the water flow and solute transport in the groundwater. Such methods have significant limitations:

[0003] (1) The surface water model and the underground water model are independent of each other, and no real-time synchronization and mutual feedback mechanism for key variables has been established at the interface.

[0004] (2) Such methods only consider the one-way effect of surface water on groundwater, which leads to the overall response process being underestimated or misjudged, resulting in inaccurate simulation results.

[0005] Therefore, how to effectively consider the feedback effect of groundwater on surface water and accurately simulate the interaction between surface water and groundwater is a technical problem that urgently needs to be solved. Summary of the Invention

[0006] This invention provides a fully coupled simulation method and system for surface water and groundwater, aiming to overcome the shortcomings of traditional coupling methods, ensure the conservation of flux at the interface, and realize the overall study of the interaction between surface water and groundwater.

[0007] Firstly, this invention provides a fully coupled simulation method for surface water and groundwater. This method is executed by a computing device, which can be understood as a computer or similar device; however, this invention does not limit the scope of the computing device. The method includes:

[0008] The spatial distribution characteristics of surface water and underground water bodies are obtained, and the surface water and underground water bodies are divided into separate grids, with interfaces established between them. The grid division and interface settings determine the location of the interfaces and the spatial correspondence for bidirectional data exchange. Based on the grid division results and interfaces, numerical models of the surface water and underground water bodies are constructed. The hydrodynamic control equations and solute transport equations for the surface water and underground water bodies are solved using these models. During the calculation of these equations, the water flux and solute flux at the interfaces are also calculated. These equations describe the flow process of the water body and the solute transport process, respectively. The water flux and solute flux are synchronized in real time to the control equations of the other region to complete the simulation calculation. The control equations of the other region are either the hydrodynamic control equations or the solute transport equations for the surface water or underground water bodies.

[0009] Through the above-described method, this invention divides the surface water and groundwater into separate grids and sets the interface between them, thus dividing the overall computational domain into two subdomains: surface water and groundwater. Numerical models for both are then constructed. Subsequently, based on the hydrodynamic control equations and solute transport equations for both surface and groundwater, the water flux and solute flux at the interface are calculated. These fluxes are dynamically updated over time, and the hydrodynamic control equations and solute transport equations are adjusted accordingly. This approach accurately reflects the bidirectional recharge relationship between surface water and groundwater, significantly improving the stability and accuracy of the simulation results. It not only supports the simulation of the migration and transformation of pollutants during the exchange between surface and groundwater but also provides more reliable technical support for numerical research on complex water environment problems.

[0010] In the aforementioned fully coupled simulation method for surface water and groundwater, based on the partitioned grid division results and interfaces, numerical models of surface water and groundwater are constructed, including:

[0011] A two-way conserved flux boundary condition is set at the interface. The two-way conserved flux boundary condition is used to realize synchronous feedback coupling calculation between surface water and underground water.

[0012] In this way, the present invention sets bidirectional conservation flux boundary conditions. Through the conservation pairing of "one outflow is the other inflow", the global conservation of water and solute is strictly guaranteed numerically, avoiding mass loss or cumulative error caused by improper boundary treatment in traditional coupling methods.

[0013] In the aforementioned fully coupled simulation method for surface water and groundwater, the hydrodynamic governing equations for surface water include the continuity equation and the incompressible Navier-Stokes equations:

[0014] ;

[0015] ;

[0016] Where u represents the fluid velocity, S SW This represents the surface water volume source term caused by the exchange between surface water and groundwater interfaces, where t represents time and p represents fluid pressure. μ represents fluid density, and μ represents fluid dynamic viscosity.

[0017] The hydrodynamic governing equations for underground water bodies include Darcy's equation and the head equation:

[0018] ;

[0019] ;

[0020] Where q represents the Darcy velocity, K represents the permeability coefficient, h represents the total head, and p represents the fluid pressure. Let z represent the fluid density and z represent the elevation head. The storage rate is represented by t, time is represented by ∇·q, the net outflow of water through the boundary of the unit cell is represented by Q, and the amount of water added or removed from the outside is represented by S. GW This indicates the groundwater volume source term caused by the exchange between surface water and groundwater interfaces;

[0021] In the aforementioned fully coupled simulation method for surface water and groundwater, the solute transport equations include the convection-diffusion-reaction equations:

[0022] ;

[0023] Where θ represents porosity, v represents the velocity vector, and the surface water is taken as... underground water C represents concentration, D represents hydrodynamic dispersion coefficient, R(C) represents solute reaction term, and S c This indicates the solute source term caused by the exchange between surface water and groundwater interfaces;

[0024] In the aforementioned fully coupled simulation method for surface water and groundwater, the calculation of water flux and solute flux includes:

[0025] Calculating water flux includes:

[0026] ;

[0027] in, For riverbed sediment permeability, For equivalent thickness, For surface water head, Groundwater head; when Time indicates surface water infiltration, when This indicates groundwater recharge;

[0028] The conservation relationship of water flux is:

[0029] ;

[0030] in, The interface area is used to ensure that the amount of surface water loss equals the amount of groundwater gain.

[0031] Calculate the solute flux at the interface, including:

[0032] ;

[0033] Where C represents the solute concentration, D b The proportionality coefficient representing the diffusion process, This represents the concentration gradient along the interface normal.

[0034] The conservation condition for interfacial solute flux is:

[0035] ;

[0036] Among them, J SW and J GW These represent the changes in solutes in surface water and groundwater, respectively.

[0037] Through the above methods, this invention employs a source term transfer and synchronous feedback mechanism to dynamically update water flux and solute flux during the time progression and transmit them to the other equation in real time. This can truly reflect the bidirectional recharge relationship between surface water and groundwater, significantly improving the stability and accuracy of the simulation results.

[0038] In the aforementioned fully coupled simulation method for surface water and groundwater, the following conservation relationship is satisfied:

[0039] ;

[0040] .

[0041] Through the above method, the present invention can strictly guarantee the global conservation of water volume and solute by the conservation pairing of "one outflow is the other inflow", thus avoiding the mass loss or cumulative error caused by improper boundary treatment in traditional coupling methods.

[0042] In the aforementioned fully coupled simulation method for surface water and groundwater, the source terms include interfacial water quantity source terms and interfacial solute source terms; the transmission process of the interfacial water quantity source terms is as follows:

[0043] For surface water bodies:

[0044] ;

[0045] in, This corresponds to the volume of a surface water unit. When This indicates surface water loss, with the source term being negative; if This indicates that surface water is replenished, and the source term is positive.

[0046] For underground water bodies:

[0047] ;

[0048] in, To correspond to the groundwater unit volume, the amount of groundwater gained is strictly equal to the amount of surface water lost;

[0049] The process of solute source term transfer at the interface is as follows:

[0050] For surface water bodies:

[0051] ;

[0052] For underground water bodies:

[0053] .

[0054] By adding a term to the surface water control equation to represent the flow loss caused by the interface, and adding a term to the groundwater control equation to represent the flow gain caused by the interface, the amount of groundwater gained is strictly equal to the amount of surface water lost. Simultaneously, for both surface and groundwater water systems, the increase or decrease in solute caused by interface exchange is calculated to ensure the conservation of solute exchange on both sides of the interface.

[0055] Secondly, the present invention provides a fully coupled simulation system for surface water and groundwater, which is used to implement the relevant content of the fully coupled simulation method for surface water and groundwater provided by the present invention, including a feature acquisition unit, a model construction unit, a flux calculation unit and a synchronous calculation unit;

[0056] The system comprises the following components: a feature acquisition unit, used to acquire the spatial distribution characteristics of surface and underground water bodies, dividing them into grids and setting the interfaces; the grid division and interface setting are used to determine the location of the interfaces and the spatial correspondence for bidirectional data exchange between surface and underground water bodies; a model building unit, used to construct numerical models of surface and underground water bodies based on the grid division results and interfaces; a flux calculation unit, used to solve the hydrodynamic control equations and solute transport equations of surface and underground water bodies based on the numerical models, obtaining water flux and solute flux; the hydrodynamic control equations and solute transport equations are used to describe the flow process of water bodies and the solute transport process, respectively; and a synchronization calculation unit, used to synchronize the water flux and solute flux to the control equations of the other region, completing the simulation calculation; the control equations of the other region are the hydrodynamic control equations or solute transport equations of surface or underground water bodies.

[0057] Thirdly, the present invention also provides a computing device, comprising: a memory for storing program instructions; and a processor for calling the program instructions stored in the memory and executing the method described in the first aspect according to the obtained program instructions.

[0058] Fourthly, the present invention also provides a computer-readable storage medium storing computer-readable instructions, which, when read and executed by a computer, implement the method of the first aspect described above.

[0059] Fifthly, the present invention provides a computer program product comprising a computer program executable by a computer device, wherein when the program is run on the computer device, the computer device performs the method described in the first aspect.

[0060] Beneficial Effects: Through the above-described method, this invention divides the surface water and groundwater into separate grids and sets the interface between them, thus dividing the overall computational domain into two subdomains: surface water and groundwater. Numerical models for surface water and groundwater are then constructed separately. Subsequently, based on the hydrodynamic control equations and solute transport equations for surface and groundwater, the water flux and solute flux for both are calculated. These fluxes are dynamically updated over time, and the hydrodynamic control equations are adjusted accordingly. This approach accurately reflects the bidirectional recharge relationship between surface water and groundwater, significantly improving the stability and accuracy of the simulation results. It not only supports the simulation of pollutant migration and transformation during the exchange between surface and groundwater but also provides more reliable technical support for numerical research on complex water environment problems. Attached Figure Description

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

[0062] Figure 1 This is a schematic diagram of the surface water area and the groundwater area in a fully coupled simulation method for surface water and groundwater provided in Embodiment 1 of the present invention;

[0063] Figure 2 This is a flowchart illustrating a fully coupled simulation method for surface water and groundwater provided in Embodiment 1 of the present invention.

[0064] Figure 3 This is a flowchart illustrating a fully coupled simulation method for surface water and groundwater provided in Embodiment 2 of the present invention.

[0065] Figure 4 This is a schematic diagram of a fully coupled simulation system for surface water and groundwater provided in Embodiment 3 of the present invention;

[0066] Figure 5 This is a schematic diagram of the structure of a computing device provided in Embodiment 4 of the present invention. Detailed Implementation

[0067] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings.

[0068] In the following embodiments of the present invention, "and / or" describes the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist simultaneously, and B exists alone, where A and B can be singular or plural. The character " / " generally indicates that the associated objects before and after are in an "or" relationship. "At least one (item) below" or similar expressions refer to any combination of these items, including any combination of single item (item) or plural items (items). For example, at least one (item) of a, b, or c can represent: a, b, c, a - b, a - c, b - c, or a - b - c, where a, b, and c can be single or multiple. The singular forms of expression "one", "a kind of", "", "the above", "the", and "this" are also intended to include expressions such as "one or more", unless there is a clear indication to the contrary in the context. Also, unless otherwise stated, the ordinal numbers such as "first", "second", etc. mentioned in the embodiments of the present invention are used to distinguish multiple objects and are not used to limit the order, time sequence, priority, or importance of multiple objects.

[0069] Describing reference to "one embodiment" or "some embodiments" etc. in the specification of the present invention means that in one or more embodiments of the present invention, specific features, structures, or characteristics described in combination with this embodiment are included. Thus, statements such as "in one embodiment", "in some embodiments", "in other some embodiments", "in still other embodiments", etc. that appear in different places in this specification do not necessarily refer to the same embodiment, but mean "one or more but not all embodiments", unless otherwise specifically emphasized in other ways. The terms "comprise", "include", "have" and their variants all mean "include but not limited to", unless otherwise specifically emphasized in other ways.

[0070] Embodiment 1

[0071] Embodiment 1 of the present invention provides a full - coupling simulation method for surface water and groundwater, as follows Figure 1 As shown in the schematic diagram of the two - way exchange between the surface water area and the groundwater area, after determining the interface between the surface water area and the groundwater area, the two - way exchange between the surface water area and the groundwater area can be judged. Embodiment 1 of the present invention provides a full - coupling simulation method for surface water and groundwater, and the specific steps are as follows Figure 2 As shown. This method is executed by a computing device, and the computing device can be understood as a device such as a computer, which is not limited in the present invention. The steps of this method are as follows:

[0072] Step 101, obtain the spatial distribution characteristics of the surface water area and the groundwater area, respectively perform zoning grid division on the surface water area and the groundwater area, and set the interface between the surface water area and the groundwater area.

[0073] Among them, the partitioning of grids and the setting of interfaces are used to determine the location of the interface between surface water and underground water and the spatial correspondence of two-way data exchange.

[0074] Specifically, after obtaining the spatial distribution characteristics of surface water and underground water, the entire computational domain is divided into two subdomains: surface water and underground water, with independent grids established for each. At the interface Γ between the two domains, a shared interface cell is defined as an information exchange area to enable data transfer between the two domains.

[0075] If the grids of surface water and underground water are inconsistent, information can be conveyed through interface mapping or interpolation.

[0076] Step 102: Based on the partitioned grid division results and the interface, construct the surface water numerical model and the groundwater numerical model.

[0077] Specifically, after dividing the overall computational domain into two subdomains—surface water and underground water—independent governing equations are established for each. For surface water, the continuity equation and the incompressible Navier-Stokes equations are used as the hydrodynamic governing equations; for underground water, the Darcy equation and the head equation are used.

[0078] In the aforementioned fully coupled simulation method for surface water and groundwater, the hydrodynamic governing equations for surface water are constructed using the continuity equation and the incompressible Navier-Stokes equations. These equations are as follows:

[0079] ;

[0080] ;

[0081] Where u represents the fluid velocity, S SW This represents the surface water volume source term caused by the exchange between surface water and groundwater interfaces, where t represents time and p represents fluid pressure. μ represents fluid density, and μ represents fluid dynamic viscosity.

[0082] The Darcy equation and the head equation are used as the hydrodynamic governing equations for underground water bodies. The specific equations are as follows:

[0083] ;

[0084] ;

[0085] Where q represents the Darcy velocity, K represents the permeability coefficient, h represents the total head, and p represents the fluid pressure. Let z represent the fluid density and z represent the elevation head. The storage rate is represented by t, time is represented by ∇·q, the net outflow of water through the boundary of the unit cell is represented by Q, and the amount of water added or removed from the outside is represented by S. GW This indicates the groundwater volume source term caused by the exchange between surface water and groundwater interfaces.

[0086] Step 103: Based on the numerical models of surface water bodies and underground water bodies, numerically solve the hydrodynamic control equations and solute transport equations for surface water bodies and underground water bodies.

[0087] Specifically, based on the water exchange between the two regions, the water flux is obtained according to the hydrodynamic control equations of surface water and groundwater, and the solute flux is obtained according to the solute transport equations.

[0088] The solute transport equation specifically includes the solute convection-diffusion-reaction equation, which is a universal equation for both surface and underground water bodies. This universality is achieved by changing the porosity of the solution in both surface and underground water environments. The equation is as follows:

[0089] ;

[0090] Where θ represents porosity, v represents the velocity vector, C represents concentration, and D represents the hydrodynamic dispersion coefficient. (The last part, "Surface water," appears to be a typo and can be omitted.) underground water R(C) represents the solute reaction term, S c This represents the solute source term caused by the exchange at the surface water-groundwater interface. Discrete solutions are applied to the above equations to simulate the migration and transformation of pollutants during the surface water-groundwater coupling process.

[0091] Specifically, the convection-diffusion-reaction equation is the equation for the transport process of solutes in surface water and groundwater. It carries the function of predicting the transport of solutes in groundwater and surface water. In the simulation of solute transport, the solute concentration in surface water is usually regarded as a constant, which cannot simulate the actual convection-diffusion-reaction process, resulting in the incorrect estimation of solute concentration and thus causing simulation errors. This invention calculates the exchange relationship between the two domains and obtains a universal transport process equation for both domains based on the convection-diffusion-reaction equation, thereby realizing the exchange of solutes between surface water and groundwater domains.

[0092] Step 104: In the process of numerically solving the hydrodynamic control equations and solute transport equations for surface water and underground water, calculate the water flux and solute flux at the interface, and transmit the water flux and solute flux to the control equations of the other region in real time to complete the simulation calculation.

[0093] The governing equations for the other region are the hydrodynamic governing equations or solute transport equations for the surface water or underground water.

[0094] Calculate water flux and solute flux, where water flux refers to the volume of water exchanged through the exchange surface per unit time, and solute flux represents the amount of solute exchanged through the exchange surface per unit time.

[0095] At the interface Γ, the water flux is calculated using the hydrodynamic equation, as shown below:

[0096] ;

[0097] in, For riverbed sediment permeability, For equivalent thickness, For surface water head, Groundwater head; when Time indicates surface water infiltration, when This indicates groundwater recharge;

[0098] The conservation relationship of water flux is:

[0099] ;

[0100] in, The interface area is used to ensure that the amount of surface water loss equals the amount of groundwater gain.

[0101] The solute flux at the interface is calculated using the solute transport equation, as shown below:

[0102] ;

[0103] Where C represents the solute concentration, D b The proportionality coefficient representing the diffusion process, This represents the concentration gradient along the interface normal.

[0104] The conservation condition for interfacial solute flux is:

[0105] ;

[0106] Among them, J SW and J GW These represent the changes in solutes in surface water and groundwater, respectively.

[0107] Synchronizing water flux and solute flux into the governing equations of the other region means adjusting the water flux based on the values ​​of water flux and solute flux. and solute flux The flux is added as a source term to the governing equations of the other region. Under the governing equations containing the above source term, the variable values ​​at the next time step are obtained by iterative solution. The above process is repeated, and the flux is updated in real time at each time step according to the latest head difference and solute concentration difference. The source term is transmitted in real time, realizing the coupling of bidirectional conserved flux boundary conditions and real-time synchronous feedback between surface water and groundwater.

[0108] The process of water source term transfer at the interface is as follows:

[0109] For surface water bodies, a new term is added to the governing equations to represent the flow loss caused by the interface:

[0110] ;

[0111] in, This corresponds to the volume of a surface water unit. When This indicates surface water loss, with the source term being negative; if This indicates that surface water is replenished, and the source term is positive.

[0112] For underground water bodies, a new term is added to the governing equations to represent the flow gain caused by the interface:

[0113] ;

[0114] in, To correspond to the groundwater unit volume, the amount of groundwater gained is strictly equal to the amount of surface water lost.

[0115] The process of solute source term transfer at the interface is as follows:

[0116] For surface water bodies, add a term to the solute transport equation:

[0117] ;

[0118] For underground water bodies, add a term to the solute transport equation:

[0119] ;

[0120] Water flux at each time step and solute flux The flux changes in real time with the head difference and concentration difference. Once the flux calculation is complete, it is immediately written into the corresponding source term and passed to the other equation, creating instantaneous feedback between surface water and groundwater. This synchronous feedback ensures that the numerical model not only satisfies local conservation but also dynamically reflects the recharge process of groundwater to surface water and the infiltration effect of surface water on groundwater. At any given time, the following conservation relationships are satisfied:

[0121] ;

[0122] .

[0123] Through the above methods, the present invention can achieve fully coupled simulation of surface water and underground water, realizing "outflow from one side is inflow from the other side", ensuring bidirectional conservation and synchronous feedback.

[0124] Example 2

[0125] The fully coupled simulation method for surface water and groundwater provided in Embodiment 2 of this invention is an improvement on the method provided in Embodiment 1. The specific process is as follows: Figure 3 As shown, the method is executed by a computing device, which can be understood as a computer or similar device, but is not limited to this in this invention. The steps of the method are as follows:

[0126] Step 201: Obtain the spatial distribution characteristics of the surface water model and the underground water model.

[0127] Determine the extent of surface water and underground water bodies, and obtain their spatial distribution characteristics.

[0128] Step 202: Perform grid division and interface settings.

[0129] The entire computational domain is divided into two subdomains: surface water and underground water, each with its own independent grid. At the interface Γ between the two domains, a shared interface cell is defined as the information exchange area, enabling data transfer between the two domains. The grid partitioning and interface settings determine the location of the interface between the surface water and underground water areas and the spatial correspondence for bidirectional data exchange.

[0130] Step 203: Set the boundary conditions for bidirectional conserved flux.

[0131] The hydrodynamic governing equations and solute transport equations for surface and underground water bodies are determined. At the interface Γ, the water flux and interfacial solute flux are calculated. The water flux and interfacial solute flux satisfy the conservation conditions, thus realizing the setting of bidirectional conservation flux boundary conditions.

[0132] Step 204: Perform source item transfer and synchronization feedback.

[0133] Water flux at each time step and solute flux The flux changes in real time with the head difference and concentration difference. After the flux calculation is completed, it is immediately written into the corresponding source term and passed to the other equation, so that an instantaneous feedback is formed between surface water and groundwater. This synchronous feedback ensures that the numerical model not only satisfies local conservation, but also dynamically reflects the process of groundwater recharge to surface water and the infiltration of surface water into groundwater.

[0134] Step 205 determines whether the process includes solute transport. If yes, return to step 203; otherwise, proceed to step 206.

[0135] Step 206: Perform time-progression and iterative solution.

[0136] Within each time step, the following process is executed sequentially:

[0137] i. Calculate the surface water flow field;

[0138] ii. Calculate the groundwater head field;

[0139] iii. Calculate water flux and solute flux;

[0140] iv. Use the calculated flux as a source term to propagate back to the other equation;

[0141] v. Correct the results in both domains and complete the simulation calculation.

[0142] Step 207: Determine whether the conservation error exceeds the threshold. If yes, return to step 202; otherwise, output the result.

[0143] If the conservation error exceeds the threshold, the above process is repeated within the same time step until the convergence condition is met. The convergence criterion can be expressed as:

[0144] ;

[0145] .

[0146] Based on the above approach, this invention introduces a bidirectional conservation flux boundary condition at the interface between surface water and groundwater. Through the conservation pairing of "outflow from one side is inflow from the other," it can strictly guarantee the global conservation of water quantity and solute in numerical terms, avoiding mass loss or cumulative errors caused by improper boundary treatment in traditional coupling methods. Simultaneously, by employing a source term transfer and synchronous feedback mechanism, water flux and solute flux are dynamically updated during time progression and immediately transferred to the other equation, accurately reflecting the bidirectional recharge relationship between surface water and groundwater, significantly improving the stability and accuracy of the simulation results. The provided water flux formula has a clear physical meaning and a simple calculation method, making it easy to couple with existing finite volume method or finite element method models, exhibiting good scalability and operability. Finally, this invention, combined with a solute migration-reaction transport mechanism, supports the simulation of the migration and transformation of pollutants during the exchange process between surface water and groundwater, providing more reliable technical support for numerical research on complex water environment problems.

[0147] Example 3

[0148] Embodiment 3 of the present invention provides a fully coupled simulation system for surface water and groundwater based on Embodiment 1. This system is used to implement the relevant content of the fully coupled simulation method for surface water and groundwater provided by the present invention, including a feature acquisition unit, a model construction unit, a flux calculation unit, and a synchronous calculation unit.

[0149] The system comprises the following components: a feature acquisition unit, used to acquire the spatial distribution characteristics of surface water and underground water, dividing the surface water and underground water into separate grids and setting the interfaces between them; the grid division and interface setting are used to determine the location of the interfaces between surface water and underground water and the spatial correspondence for bidirectional data exchange; a model building unit, used to construct numerical models of surface water and underground water based on the grid division results and interfaces; and a flux calculation unit, used to calculate the flux of surface water based on the numerical models of surface water and underground water. The system numerically solves the hydrodynamic control equations and solute transport equations for surface and underground water bodies. These equations describe the flow and transport processes of water bodies. A synchronous computing unit is used to calculate the water flux and solute flux at the interface during the numerical solution of the hydrodynamic control equations and solute transport equations for surface and underground water bodies. The water flux and solute flux are then synchronized to the control equations of the other region to complete the simulation calculation. The control equations of the other region are either the hydrodynamic control equations or the solute transport equations for surface or underground water bodies.

[0150] Example 4

[0151] Having introduced the fully coupled simulation system for surface water and groundwater in an exemplary embodiment of the present invention, we will now introduce a computing device in another exemplary embodiment of the present invention.

[0152] Those skilled in the art will understand that various aspects of the present invention can be implemented as systems, methods, or program products. Therefore, various aspects of the present invention can be specifically implemented in the following forms: entirely in hardware, entirely in software (including firmware, microcode, etc.), or in a combination of hardware and software, collectively referred to herein as “circuit,” “module,” or “system.”

[0153] In some possible implementations, the computing device according to the invention may include at least one processor and at least one memory. The memory stores a computer program that, when executed by the processor, causes the processor to perform the steps in the fully coupled surface water and groundwater simulation method according to various exemplary embodiments of the invention described above.

[0154] The following reference Figure 5To describe a computing device 130 according to this embodiment of the invention. Figure 5 The computing device 130 shown is merely an example and should not be construed as limiting the functionality and scope of the embodiments of the present invention. Figure 5 As shown, the computing device 130 is presented in the form of a general-purpose smart terminal (or Bluetooth headset). The components of the computing device 130 may include, but are not limited to: at least one processor 131, at least one memory 132, and a bus 133 connecting different system components (including memory 132 and processor 131).

[0155] Bus 133 represents one or more of several bus architectures, including a memory bus or memory controller, peripheral bus, processor, or local bus using any of the various bus architectures. Memory 132 may include readable media in the form of volatile memory, such as random access memory (RAM) 1321 and / or cache memory 1322, and may further include read-only memory (ROM) 1323. Memory 132 may also include a program / utility 1325 having a set (at least one) of program modules 1324, including but not limited to: an operating system, one or more application programs, other program modules, and program data, each or some combination of these examples may include an implementation of a network environment.

[0156] The computing device 130 can also communicate with one or more external devices 134 (e.g., keyboard, pointing device, etc.), and / or with any device that enables the computing device 130 to communicate with one or more other smart terminals (e.g., router, modem, etc.). This communication can be performed via input / output (I / O) interface 135. Furthermore, the computing device 130 can also communicate with one or more networks (e.g., local area network (LAN), wide area network (WAN), and / or public networks, such as the Internet) via network adapter 136. As shown, network adapter 136 communicates with other modules used in the computing device 130 via bus 133. It should be understood that, although not shown in the figures, other hardware and / or software modules can be used in conjunction with the computing device 130, including but not limited to: microcode, device drivers, redundant processors, external disk drive arrays, RAID systems, tape drives, and data backup storage systems.

[0157] In some possible implementations, various aspects of the fully coupled surface water and groundwater simulation method provided by the present invention can also be implemented in the form of a program product, which includes a computer program that, when run on a computer device, causes the computer device to perform the steps in the fully coupled surface water and groundwater simulation method according to various exemplary embodiments of the present invention as described above.

[0158] The program product may employ any combination of one or more readable media. A readable medium may be a readable signal medium or a readable storage medium. A readable storage medium may be, for example—but not limited to—an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of readable storage media (a non-exhaustive list) include: an electrical connection having one or more wires, a portable disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0159] The program product for time-domain noise processing according to embodiments of the present invention may employ a portable compact disc read-only memory (CD-ROM) and include a computer program, and may run on a smart terminal. However, the program product of the present invention is not limited thereto. In this document, the readable storage medium may be any tangible medium containing or storing a program that may be used by or in conjunction with an instruction execution system, apparatus, or device.

[0160] A readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, carrying a readable computer program. Such propagated data signals may take various forms, including—but not limited to—electromagnetic signals, optical signals, or any suitable combination thereof. A readable signal medium may also be any readable medium other than a readable storage medium, capable of sending, propagating, or transmitting a program for use by or in conjunction with an instruction execution system, apparatus, or device.

[0161] It should be noted that although several units or sub-units of the device have been mentioned in the detailed description above, this division is merely exemplary and not mandatory. In fact, according to embodiments of the present invention, the features and functions of two or more units described above can be embodied in one unit. Conversely, the features and functions of one unit described above can be further divided and embodied by multiple units.

[0162] Furthermore, although the operations of the method of the present invention are described in a specific order in the accompanying drawings, this does not require or imply that these operations must be performed in that specific order, or that all the operations shown must be performed to achieve the desired result. Additionally or alternatively, certain steps may be omitted, multiple steps may be combined into one step, and / or one step may be broken down into multiple steps.

[0163] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable access frequency prediction device to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable access frequency prediction device, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0164] These computer program instructions may also be stored in a computer-readable storage medium capable of directing a computer or other programmable access predictive device to operate in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0165] These computer program instructions can also be loaded onto a computer or other programmable access device to cause a series of operational steps to be performed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable device for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0166] Although preferred embodiments of the invention have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the invention.

[0167] Obviously, those skilled in the art can make various modifications and variations to this invention without departing from its spirit and scope. Therefore, if these modifications and variations fall within the scope of the claims of this invention and their equivalents, this invention also intends to include these modifications and variations.

Claims

1. A fully coupled simulation method for surface water and groundwater, characterized in that, include: The spatial distribution characteristics of surface water and underground water are obtained, and the surface water and underground water are divided into grids and partitions respectively, and the interface between the surface water and underground water is set. The grid division and the setting of the interface are used to determine the location of the interface between the surface water and underground water and the spatial correspondence of bidirectional data exchange. Based on the partitioned grid division results and the interface, a numerical model of surface water and a numerical model of underground water are constructed. Based on the numerical models of surface water bodies and groundwater bodies, the hydrodynamic governing equations and solute transport equations for both surface and groundwater bodies are numerically solved. The hydrodynamic governing equations describe the flow process of the water body, and the solute transport equations describe the solute transport process. The hydrodynamic governing equations for surface water bodies include the continuity equation and the incompressible Navier-Stokes equations. ； ； Wherein, u represents the fluid velocity, and S SW This represents the surface water volume source term caused by the exchange between surface water and groundwater interfaces, where t represents time, p represents fluid pressure, and so on. μ represents the fluid density, and μ represents the fluid dynamic viscosity. The hydrodynamic governing equations for the underground water area include Darcy's equation and the head equation: ； ； Wherein, q represents Darcy velocity, K represents permeability coefficient, h represents total head, and p represents fluid pressure. The fluid density is represented by z, and the elevation head is represented by z. The storage rate is represented by t, time is represented by ∇·q, the net outflow of water through the boundary of the unit cell is represented by Q, and the amount of water added or removed from the outside is represented by S. GW This indicates the groundwater volume source term caused by the exchange between surface water and groundwater interfaces; In the process of numerically solving the hydrodynamic control equations and solute transport equations of the surface water and underground water areas, the water flux and solute flux at the interface are calculated, and the water flux and solute flux are synchronously transferred to the control equations of the other region, thereby completing the fully coupled simulation calculation; wherein, the control equation of the other region is the hydrodynamic control equation or solute transport equation of the surface water or underground water area; the calculation of water flux and solute flux includes: Calculate water flux ,include: ； Among them, the The permeability of riverbed sediments, the For equivalent thickness, the For surface water head, the aforementioned Groundwater head; when Time indicates surface water infiltration, when This indicates groundwater recharge; The conservation relationship of the water flux is as follows: ； Among them, the The quantity of surface water flux, the The quantity of groundwater flux, the The interface area is used to ensure that surface water loss equals groundwater gain. Indicates water flux; Calculate the solute flux at the interface ,include: ； Wherein, C represents the solute concentration, and D... b The proportionality coefficient representing the diffusion process, the This represents the concentration gradient along the interface normal. The conservation relationship for the interfacial solute flux is as follows: ； Wherein, J SW and J GW These represent the changes in solutes in surface water and groundwater, respectively.

2. The method according to claim 1, characterized in that, The step of constructing a surface water numerical model and a groundwater numerical model based on the partitioned grid division results and the interface includes: A bidirectional conserved flux boundary condition is set at the interface. The bidirectional conserved flux boundary condition is used to realize the bidirectional transfer of water flux and solute flux between surface water and groundwater, so as to complete the synchronous feedback coupling calculation between the surface water and groundwater.

3. The method according to claim 1, characterized in that, The solute transport equations include convection-diffusion-reaction equations: ； Wherein, θ represents porosity, v represents velocity vector, and surface water is taken as... underground water source C represents the solute concentration, D represents the hydrodynamic dispersion coefficient, R(C) represents the solute reaction term, and S represents the solute reaction term. c This indicates the solute source term caused by the exchange between surface water and groundwater.

4. The method according to claim 1, characterized in that, The fully coupled simulation method for surface water and groundwater satisfies the following conservation relationship: ； 。 5. According to the method described in claim 1, water flux and solute flux are converted into source terms required by each region and simultaneously written into the control equation of the other region. Under the control equation containing the above source terms, the variable values ​​at the next time step are obtained by iterative solution. The above process is repeated, and the flux is updated in real time at each time step according to the latest head difference and solute concentration difference. The source terms are transmitted in real time, realizing the coupling of bidirectional conserved flux boundary conditions and real-time synchronous feedback between surface water and groundwater.

6. The method according to claim 5, wherein the source term includes an interfacial water source term and an interfacial solute source term; wherein, The process of transferring the interface water source term is as follows: For surface water bodies: ； Among them, the For the corresponding surface water unit volume; when the This indicates surface water loss, with a negative source term; if the stated This indicates that surface water is replenished, and the source term is positive. For underground water bodies: ； Among them, the To correspond to the groundwater unit volume, the amount of groundwater gained is strictly equal to the amount of surface water lost; The process of solute source term transfer at the interface is as follows: For the aforementioned surface water bodies: ； Regarding the aforementioned underground water area: 。 7. A fully coupled simulation system for surface water and groundwater, characterized in that, include: The feature acquisition unit is used to acquire the spatial distribution features of surface water and underground water, divide the surface water and underground water into separate grids, and set the interface between the surface water and underground water; the grid division and the setting of the interface are used to determine the location of the interface between the surface water and underground water and the spatial correspondence of bidirectional data exchange. The model building unit is used to build a numerical model of surface water and a numerical model of underground water based on the partitioned grid division results and the interface. The flux calculation unit is used to numerically solve the hydrodynamic control equations and solute transport equations for surface water and groundwater bodies based on the numerical models of the surface water and groundwater bodies. These hydrodynamic control equations and solute transport equations describe the flow process of water bodies and the solute transport process. The hydrodynamic control equations for surface water bodies include the continuity equation and the incompressible Navier-Stokes equations. ； ； Wherein, u represents the fluid velocity, and S SW This represents the surface water volume source term caused by the exchange between surface water and groundwater interfaces, where t represents time, p represents fluid pressure, and so on. μ represents the fluid density, and μ represents the fluid dynamic viscosity. The hydrodynamic governing equations for the underground water area include Darcy's equation and the head equation: ； ； Wherein, q represents Darcy velocity, K represents permeability coefficient, h represents total head, and p represents fluid pressure. The fluid density is represented by z, and the elevation head is represented by z. The storage rate is represented by t, time is represented by ∇·q, the net outflow of water through the boundary of the unit cell is represented by Q, and the amount of water added or removed from the outside is represented by S. GW This indicates the groundwater volume source term caused by the exchange between surface water and groundwater interfaces; A synchronous calculation unit is used to calculate the water flux and solute flux at the interface during the numerical solution of the hydrodynamic control equations and solute transport equations of the surface water and underground water areas, and to synchronize the water flux and solute flux to the control equations of the other region to complete the simulation calculation; wherein, the control equations of the other region are the hydrodynamic control equations or solute transport equations of the surface water or underground water areas; the calculation of water flux and solute flux includes: Calculate water flux ,include: ； Among them, the The permeability of riverbed sediments, the For equivalent thickness, the For surface water head, the aforementioned Groundwater head; when Time indicates surface water infiltration, when This indicates groundwater recharge; The conservation relationship of the water flux is as follows: ； Among them, the The quantity of surface water flux, the The quantity of groundwater flux, the The interface area is used to ensure that surface water loss equals groundwater gain. Indicates water flux; Calculate the solute flux at the interface ,include: ； Wherein, C represents the solute concentration, and D... b The proportionality coefficient representing the diffusion process, the This represents the concentration gradient along the interface normal. The conservation relationship for the interfacial solute flux is as follows: ； Wherein, J SW and J GW These represent the changes in solutes in surface water and groundwater, respectively.

8. A computing device, characterized in that, include: Memory, used to store program instructions; A processor is configured to invoke program instructions stored in the memory and execute the method as described in any one of claims 1-6 according to the obtained program instructions.

9. A computer-readable storage medium, characterized in that, Includes computer-readable instructions that, when read and executed by a computer, cause the method as described in any one of claims 1 to 6 to be implemented.

Images