Method, device, equipment and medium for simulating dynamic exchange between surface water and underground water

Abstract

The application relates to a surface water-underground water dynamic exchange simulation method, device, equipment and medium, and is applied to the hydrological technical field.The method comprises the following steps: determining a simulated surface water level according to the simulated river flow of each grid at a current time; determining simulated vertical exchange amount and simulated lateral exchange amount of each grid at a current time step according to the simulated surface water level of each grid at the current time and obtained underground water level; determining surface water-underground water exchange amount of each grid at the current time step according to the simulated vertical exchange amount and the simulated lateral exchange amount of each grid at the current time step; determining simulated river flow of each grid at a next time according to the simulated river flow of each grid at the current time, the surface water-underground water exchange amount of each grid at the current time step and surface water input amount, and taking the next time as a new current time, and repeatedly executing the above steps. The method can improve the accuracy of the exchange simulation between the surface water and the underground water.

Description

Technical Field

[0001] This application relates to the field of hydrological technology, and in particular to a method, apparatus, equipment and medium for simulating the dynamic exchange between surface water and groundwater. Background Technology

[0002] Water exchange occurs between surface water and groundwater; that is, water from surface water can flow into groundwater, and vice versa. Simulating this exchange is a prerequisite for achieving coordinated water resource management and ensuring ecological security.

[0003] Currently, surface water-groundwater coupling models are often used to simulate the exchange between surface water and groundwater.

[0004] However, current surface water-groundwater coupling models often lack detailed characterization of complex human activities and natural environments, making them difficult to apply effectively in watersheds with intense human activity and diverse natural environments, resulting in inaccurate simulations of surface water-groundwater exchange. Summary of the Invention

[0005] Based on this, this application provides a method, apparatus, equipment, and medium for simulating the dynamic exchange between surface water and groundwater, which can improve the accuracy of the simulation of the exchange between surface water and groundwater.

[0006] Firstly, this application provides a method for simulating the dynamic exchange between surface water and groundwater, the method comprising:

[0007] Based on the simulated river flow at the current moment for each grid obtained by dividing the simulated watershed into latitude and longitude grids, the simulated surface water level at the current moment for each grid is determined.

[0008] Based on the simulated surface water level and the obtained groundwater level of each grid at the current time, determine the simulated vertical exchange and simulated lateral exchange of each grid at the current time step; the current time step is the time step between the current time and the next time step;

[0009] Based on the simulated vertical and lateral exchange amounts of each grid at the current time step, determine the surface water-groundwater exchange amount of each grid at the current time step.

[0010] Based on the simulated river flow of each grid at the current time, the surface water-groundwater exchange and surface water input of each grid at the current time step, determine the simulated river flow of each grid at the next time step, and take the next time step as the new current time step, and repeat the above steps.

[0011] In some embodiments, determining the simulated surface water level of each grid at the current time based on the simulated river flow of each grid obtained by dividing the simulated watershed into latitude and longitude grids includes:

[0012] Obtain the roughness, size, and slope of the river channel covered by each grid;

[0013] Based on the simulated river flow of each grid at the current moment, as well as the roughness, size, and slope of the river covered by each grid, the simulated surface water level of each grid at the current moment is determined.

[0014] In some embodiments, the simulated vertical exchange and simulated lateral exchange of each grid at the current time step are determined based on the simulated surface water level and the acquired groundwater level of each grid at the current time, including:

[0015] Obtain the riverbed hydraulic conduction coefficient and the bank hydraulic conduction coefficient of the simulated river segment in each grid.

[0016] Based on the simulated surface water level, groundwater level and riverbed hydraulic conduction coefficient of each grid at the current time, determine the simulated vertical exchange amount of each grid at the current time step.

[0017] Based on the simulated surface water level, groundwater level, and bank hydraulic conduction coefficient of the simulated river segment in each grid at the current time step, determine the simulated lateral exchange amount of each grid at the current time step.

[0018] In some embodiments, the simulated vertical exchange amount of each grid at the current time step is determined based on the simulated surface water level, groundwater level, and riverbed hydraulic conductivity coefficient of the simulated river segment in each grid at the current time, including:

[0019] Obtain the riverbed bottom elevation for each grid cell;

[0020] Based on the simulated surface water level, groundwater level, riverbed bottom elevation, and riverbed hydraulic conduction coefficient of each grid at the current time, the simulated vertical exchange volume of each grid at the current time step is determined.

[0021] In some embodiments, the simulated vertical exchange amount of each grid at the current time step is determined based on the simulated surface water level, groundwater level, riverbed bottom elevation, and riverbed hydraulic conductivity coefficient of the simulated river segment in each grid at the current time, including:

[0022] If the groundwater level of any grid at the current time is less than or equal to the riverbed bottom elevation of any grid at the current time, the simulated vertical exchange amount of any grid at the current time is determined based on the riverbed bottom elevation of any grid at the current time, the simulated surface water level, and the riverbed hydraulic conduction coefficient of the simulated river segment in any grid.

[0023] If the groundwater level of any grid at the current time is greater than the riverbed bottom elevation of any grid at the current time, the simulated vertical exchange volume of any grid at the current time is determined based on the groundwater level of any grid at the current time, the simulated surface water level, and the hydraulic conductivity coefficient of the simulated river section in any grid.

[0024] In some embodiments, the simulated river flow of each grid at the next time step is determined based on the simulated river flow of each grid at the current time step, the surface water-groundwater exchange volume of each grid at the current time step, and the surface water input volume, including:

[0025] Obtain the surface water input for each grid at the current time step, including surface runoff input, water return flow, river surface evaporation, and river water intake.

[0026] Based on the simulated river flow of each grid at the current time step, the surface water-groundwater exchange, surface runoff input, water return flow, river surface evaporation, and river water intake of each grid at the current time step, determine the simulated river flow of each grid at the next time step.

[0027] In some embodiments, the method further includes:

[0028] Obtain the surface water input of each grid at the current time step, including surface runoff input, water return flow, river surface evaporation, and river water intake, as well as the river water storage at the current moment.

[0029] Based on the surface runoff input, water return flow, river surface evaporation, river water intake, surface water-groundwater exchange, and river water storage at the current time step for each grid, determine the river water storage at the next time step.

[0030] Secondly, this application provides a dynamic exchange simulation device for surface water and groundwater, the device comprising:

[0031] The surface water level determination module is used to determine the simulated surface water level of each grid at the current time based on the simulated river flow of each grid obtained by dividing the simulated watershed into latitude and longitude grids.

[0032] The vertical and lateral exchange quantity determination module is used to determine the simulated vertical and lateral exchange quantities of each grid at the current time step based on the simulated surface water level and the acquired groundwater level of each grid at the current time step; the current time step is the time step between the current time step and the next time step.

[0033] The surface water-groundwater exchange volume determination module is used to determine the surface water-groundwater exchange volume of each grid at the current time step based on the simulated vertical and simulated lateral exchange volumes of each grid at the current time step.

[0034] The river flow determination module is used to determine the simulated river flow of each grid at the next time step based on the simulated river flow of each grid at the current time step, the surface water-groundwater exchange volume of each grid at the current time step, and the surface water input volume, and to take the next time step as the new current time step; each module is used to cyclically execute the corresponding steps.

[0035] Thirdly, this application provides a computer device including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps of any of the methods in the first aspect.

[0036] Fourthly, this application provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the steps of the method of any one of the first aspects.

[0037] The technical solution provided in this application calculates the simulated vertical and lateral exchange volumes within the current time step grid by grid-by-grid, and determines the surface water-groundwater exchange volume accordingly. This accurately depicts the vertical and lateral runoff exchange processes of water bodies at different spatial locations, significantly improving the physical realism and spatial refinement of the watershed water cycle simulation, and enhancing the accuracy of the surface water-groundwater exchange simulation. Furthermore, based on the simulated river flow at the current time, the surface water-groundwater exchange volume at the current time step, and the surface water input, the simulated river flow at the next time step for each grid is determined, and the next time step is taken as the new current time step. This approach comprehensively considers the surface water input of each grid at each time step, enabling effective application even in watersheds with intense human activity and diverse natural environments. For example, human activities affect surface water-groundwater exchange not only through surface water extraction but also through groundwater extraction. This complex interaction and feedback between human activities and surface water-groundwater exchange ultimately affects water resource allocation. Therefore, by considering the surface water input of each grid at the current time step, this embodiment can accurately obtain the surface water-groundwater exchange volume at each moment, further improving the accuracy of the surface water-groundwater exchange simulation. Attached Figure Description

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

[0039] Figure 1 A schematic diagram of the cross-section of the river, groundwater, and river channel provided for the first embodiment;

[0040] Figure 2 A schematic diagram of the cross-section of the river, groundwater, and river channel provided for the second embodiment;

[0041] Figure 3 A flowchart illustrating a method for simulating dynamic exchange between surface water and groundwater in some embodiments;

[0042] Figure 4 A flowchart illustrating the method for determining the simulated surface water level of each grid at the current moment, provided in some embodiments;

[0043] Figure 5 A flowchart illustrating the method for determining the simulated vertical and lateral exchange amounts of each grid at the current time step, provided in some embodiments;

[0044] Figure 6 This is a flowchart illustrating the method for determining the simulated river flow of each grid at the next time step in some embodiments;

[0045] Figure 7 This is a flowchart illustrating the method for determining the river water storage at the next moment in some embodiments;

[0046] Figure 8 A flowchart illustrating the simulated exchange between surface water and groundwater for some embodiments;

[0047] Figure 9 A schematic diagram of the structure of a dynamic exchange simulation device between surface water and groundwater provided in some embodiments;

[0048] Figure 10 A schematic diagram of the structure of a computer device provided for some embodiments. Detailed Implementation

[0049] The embodiments of the technical solution of this application will now be described in detail with reference to the accompanying drawings. These embodiments are only used to more clearly illustrate the technical solution of this application and are therefore merely examples, and should not be used to limit the scope of protection of this application.

[0050] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains; the terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the application; the terms “comprising” and “having”, and any variations thereof, in the specification, claims, and foregoing description of the drawings are intended to cover non-exclusive inclusion.

[0051] In the description of the embodiments of this application, technical terms such as "first" and "second" are used only to distinguish different objects and should not be construed as indicating or implying relative importance or implicitly specifying the number, specific order, or primary and secondary relationship of the indicated technical features. In the description of the embodiments of this application, "multiple" means two or more, "multiple groups" means two or more, and "each" means each of the multiple, unless otherwise explicitly defined.

[0052] In this document, the term "embodiment" means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of this application. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.

[0053] In the description of the embodiments in this application, the term "and / or" is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, and B existing alone. Additionally, the character " / " in this document generally indicates that the preceding and following related objects have an "or" relationship.

[0054] Unless otherwise specified, the order of execution steps in the embodiments of this application is not limited. It should also be noted that any step in the embodiments of this application can be executed independently, that is, the execution of any step in the above embodiments can be performed without depending on the execution of other steps.

[0055] Figure 1 The first embodiment provides a schematic diagram of the cross-section of the river, groundwater, and river channel. Figure 1 In the embodiment shown, the river channel 11 includes low-permeability riverbed sediments 111 and the river channel body 112, the groundwater 12 is water below the ground surface, and the river 13 is surface water.

[0056] Figure 2 The second embodiment provides a schematic diagram of the river, groundwater, and cross-section of the river channel. Figure 2In the embodiment shown, the river channel 21 includes low-permeability riverbed sediments 211 and the river channel body 212, the groundwater 22 is water below the ground surface, and the river 23 is surface water.

[0057] The methods in any embodiment of this application can be applied to a computer device or a processor in a computer device. Exemplarily, the computer device may include one or a combination of at least two of the following: a server, a mobile phone, a tablet computer, a computer with transceiver capabilities, a handheld computer, a desktop computer, an industrial control computer, a virtual reality (VR) device, an augmented reality (AR) device, or a device in industrial control, etc.

[0058] Figure 3 A flowchart illustrating a method for simulating dynamic exchange between surface water and groundwater provided in some embodiments, such as... Figure 3 As shown, the method includes the following steps:

[0059] S301. Based on the simulated river flow of each grid obtained by dividing the simulated watershed into latitude and longitude grids at the current time, determine the simulated surface water level of each grid at the current time.

[0060] In some embodiments, S301 may include: simulating the watershed to obtain a simulated watershed, and dividing the simulated watershed into grids to obtain each grid; and determining the simulated surface water level of each grid at the current time based on the simulated river flow of each grid at the current time.

[0061] A watershed is the catchment area of ​​a river or lake surrounded by the surface water-groundwater watershed line, usually referring to the surface water catchment area.

[0062] For example, the watershed can be simulated in the Community Water Model (CWatM) to obtain a simulated watershed, and the simulated watershed can be divided into latitude and longitude grids in CWatM to obtain multiple grids of the simulated watershed.

[0063] In some embodiments, gridding of the simulated watershed may include dividing the simulated watershed into latitude and longitude grids. This latitude and longitude grid division involves dividing the area of ​​the simulated watershed into several grids according to preset latitude and longitude intervals (such as 0.1°×0.1° or 0.05°×0.05°). Each grid can participate in the hydrological process simulation as an independent computational unit. For example, the preset latitude and longitude intervals can be fixed values.

[0064] In other embodiments, meshing the simulated watershed may include dividing the simulated watershed into meshes according to a preset meshing method in CWatM.

[0065] In some embodiments, the simulated river flow of each grid at the current moment can be a preset value, or it can be the actual river flow that can be measured.

[0066] In some embodiments, S301 may include: determining the simulated water depth of each grid at the current time based on the simulated river flow of each grid at the current time; and determining the simulated surface water level of each grid at the current time based on the simulated water depth of each grid at the current time and the riverbed elevation corresponding to each grid.

[0067] For example, the Manning formula can be used to determine the simulated water depth of each grid at the current moment based on the simulated river flow of each grid at the current moment.

[0068] In some embodiments, the simulated surface water level of each grid can be determined in CwatM based on the simulated river flow of each grid at the current time.

[0069] For example, CwatM can transfer the simulated surface water level of each grid at the current moment to the Modular Three-Dimensional Finite-Difference Ground-Water Flow Model (MODFLOW). In MODFLOW, the simulated watershed also needs to be meshed. For example, the meshing method for the simulated watershed in MODFLOW is the same as the meshing method for the simulated watershed in CwatM.

[0070] S302. Based on the simulated surface water level and the obtained groundwater level of each grid at the current time, determine the simulated vertical exchange amount and simulated lateral exchange amount of each grid at the current time step.

[0071] The current time step is the time step between the current moment and the next moment.

[0072] Simulated vertical exchange can represent the flow exchange between surface water and the riverbed beneath it. Simulated lateral exchange can represent the flow exchange between surface water and the riverbank on the side of the surface water.

[0073] In some embodiments, the groundwater level of each grid at the current moment can change dynamically based on complex human and natural activities. For example, human extraction of groundwater can cause the groundwater level to drop. As another example, rainfall can cause the groundwater level to rise.

[0074] For example, in MODFLOW, the simulated vertical exchange and simulated lateral exchange of each grid at the current time step are determined based on the simulated surface water level and groundwater level of each grid at the current time.

[0075] S303. Based on the simulated vertical and lateral exchange amounts of each grid at the current time step, determine the surface water-groundwater exchange amount of each grid at the current time step.

[0076] Surface water-groundwater exchange represents the amount of flow exchanged between surface water and groundwater. For example, surface water-groundwater exchange represents the difference between the flow rate of surface water flowing into groundwater and the flow rate of groundwater flowing into surface water, or the difference between the flow rate of groundwater flowing into surface water and the flow rate of surface water flowing into groundwater.

[0077] In some embodiments, the surface water-groundwater exchange volume of the simulated watershed can be determined based on the surface water-groundwater exchange volume of each grid at the current time step.

[0078] For example, in MODFLOW, the surface water-groundwater exchange rate of each grid cell at the current time step is determined based on the simulated vertical and lateral exchange rates of each grid cell at the current time step. MODFLOW then transmits the surface water-groundwater exchange rates of each grid cell at the current time step to CwatM, so that CwatM can determine the simulated channel flow rate of each grid cell at the next time step.

[0079] In some embodiments, the sum of the simulated vertical exchange and simulated lateral exchange of each grid at the current time step can be determined as the surface water-groundwater exchange of each grid at the current time step.

[0080] S304. Based on the simulated river flow of each grid at the current time, the surface water-groundwater exchange volume and surface water input of each grid at the current time step, determine the simulated river flow of each grid at the next time step.

[0081] For example, the surface water input of each grid at the current time step can represent the total flow rate of surface water input into each grid within the current time step. For example, the surface water input can be caused by complex human activities and / or by natural activities.

[0082] For example, the surface water-groundwater exchange volume represents the difference between the flow rate of groundwater output to surface water and the flow rate of surface water output to groundwater. S304 may include: determining the simulated river flow rate of each grid at the next time step by summing the simulated river flow rate of each grid at the current time step, the surface water-groundwater exchange volume of each grid at the current time step, and the surface water input flow rate.

[0083] The technical solution provided in this application calculates the simulated vertical and lateral exchange volumes within the current time step grid by grid-by-grid, and determines the surface water-groundwater exchange volume accordingly. This accurately depicts the vertical and lateral runoff exchange processes of water bodies at different spatial locations, significantly improving the physical realism and spatial refinement of the watershed water cycle simulation, and enhancing the accuracy of the surface water-groundwater exchange simulation. Furthermore, based on the simulated river flow at the current time, the surface water-groundwater exchange volume at the current time step, and the surface water input, the simulated river flow at the next time step for each grid is determined, and the next time step is taken as the new current time step. This approach comprehensively considers the surface water input of each grid at each time step, enabling effective application even in watersheds with intense human activity and diverse natural environments. For example, human activities affect surface water-groundwater exchange not only through surface water extraction but also through groundwater extraction. This complex interaction and feedback between human activities and surface water-groundwater exchange ultimately affects water resource allocation. Therefore, by considering the surface water input of each grid at the current time step, this embodiment can accurately obtain the surface water-groundwater exchange volume at each moment, further improving the accuracy of the surface water-groundwater exchange simulation.

[0084] Figure 4 A flowchart illustrating the method for determining the simulated surface water level of each grid at the current moment, as provided in some embodiments, is shown below. Figure 4 As shown, this method is one interpretation of S302, and the method includes:

[0085] S401. Obtain the roughness, size, and slope of the river channel covered by each grid.

[0086] Roughness is the combined effect of riverbed and bank materials on water flow resistance. It is usually quantified using empirical or semi-empirical parameters such as the Manning coefficient. Its value depends on the composition of the river surface (such as silt, sand, bedrock, or vegetation cover). This parameter directly determines the magnitude of the frictional resistance experienced by the water flow under specific flow velocity and water depth conditions.

[0087] The dimensions of a river channel may include its length and width. For example, a grid may cover the river channel and / or its banks, and the dimensions of the area to which the river channel is covered by the grid are determined as the dimensions of the river channel covered by the grid.

[0088] Slope refers to the longitudinal slope of a river channel along the direction of water flow, that is, the rate of decline of the riverbed elevation along the direction of water flow. This parameter is used to drive the water body to move downstream under the action of gravity.

[0089] S402. Based on the simulated river flow of each grid at the current moment, and the roughness, size and slope of the river covered by each grid, determine the simulated surface water level of each grid at the current moment.

[0090] This application uses a rectangular river channel cross-section as an example to derive the process of determining the simulated surface water level of each grid at the current moment:

[0091] The expression for Manning's formula is as follows: ;in, Let be the flow velocity of each grid, and n be the Manning coefficient, used to characterize the roughness of the channel covered by each grid. Let S be the hydraulic radius, and S be the slope of the river channel covered by each grid.

[0092] Flow rate of each grid Simulated river flow at the current time for each grid The relationship between them is ;in, This represents the cross-sectional area of ​​the water passage in each grid.

[0093] Cross-sectional area of ​​water flow in each grid The width B and water depth h of the river channel covered by each grid are given. r The product of and can be used to obtain the formula. .

[0094] The hydraulic radius R is obtained by dividing the cross-sectional area A of each grid by the wetted perimeter P (the circumference of the water flow in contact with the channel wall) of each grid. For wide and shallow channels (channel width B...), water depth h r Alternatively, for each grid cell, the hydraulic radius R is approximately equal to the water depth h. r This allows us to solve for the depth h of the river channel covered by each grid. r : .

[0095] The depth h of the river channel covered by each grid r The simulated surface water level for each grid cell is obtained by adding the riverbed elevation corresponding to each grid cell. .

[0096] In the technical solution provided in this application embodiment, the simulated surface water level of each grid at the current moment is determined based on the simulated river flow of each grid at the current moment, as well as the roughness, size, and slope of the river covered by each grid. Relying on the attribute information of the river covered by each grid and the simulated river flow of each grid at the current moment, the simulated surface water level of each grid at the current moment is dynamically determined, thereby improving the accuracy of the determined simulated surface water level of each grid at the current moment.

[0097] Figure 5 A flowchart illustrating the method for determining the simulated vertical and lateral exchange quantities of each grid at the current time step, as provided in some embodiments, is shown below. Figure 5 As shown, this method is an interpretation of S302, and includes the following steps:

[0098] S501. Obtain the hydraulic conductivity coefficient of the riverbed and the hydraulic conductivity coefficient of the riverbank in each simulated river segment in each grid.

[0099] In some embodiments, the riverbed hydraulic conductivity coefficient can be obtained by collecting riverbed sediment samples and conducting permeability tests. In some embodiments, the bank hydraulic conductivity coefficient of the simulated river section can be determined by: assigning values ​​to layers or zones based on the spatial distribution characteristics of the riverbank constituent materials, combined with geological information from geological boreholes and geophysical exploration results, and determining the bank hydraulic conductivity coefficient of the simulated river section based on the layered or zoned value assignment results.

[0100] In other embodiments, the hydraulic conductivity coefficients of the riverbed and the riverbank of the simulated river section can be calculated using formulas.

[0101] For example, determining the hydraulic conductivity coefficient of a simulated river section may include: obtaining the sediment permeability coefficient, river length, river width, and riverbed thickness of the simulated river section; and determining the hydraulic conductivity coefficient of the simulated river section based on these parameters. ;in, To simulate the hydraulic conductivity coefficient of the riverbed in the river section; To simulate the sediment permeability coefficient of the riverbed in a river section; To simulate the river length of the river section; To simulate the width of the river section; To simulate the riverbed thickness of the river section.

[0102] For example, determining the bank hydraulic conductivity coefficient of a simulated river segment may include: obtaining the sediment permeability coefficient, bank length, longitudinal thickness, and transverse thickness of the simulated river segment; and determining the bank hydraulic conductivity coefficient based on these parameters. ;in, To simulate the hydraulic conduction coefficient of the riverbank in a river section; To simulate the sediment permeability coefficient of the riverbank in a river section; To simulate the length of the riverbank in a river section; To simulate the longitudinal thickness of the riverbank in a river section (the longitudinal thickness of the riverbank can be the same as the depth of the river channel); To simulate the lateral thickness of the riverbank in a river section.

[0103] In some other embodiments, the hydraulic conductivity coefficients of the riverbed and the riverbank of the simulated river section can be determined by parameter calibration. For example, the hydraulic conductivity coefficients of the riverbed and the bank are randomly initialized. Based on these coefficients and the current groundwater depth (surface elevation minus groundwater level), the groundwater depth at a subsequent moment is continuously determined. If the groundwater depth at a subsequent moment is the same as the actual groundwater depth, then the current hydraulic conductivity coefficients of the riverbed and the bank are determined as the hydraulic conductivity coefficients of the simulated river segment. If the groundwater depth at a subsequent moment is different from the actual groundwater depth, the hydraulic conductivity coefficients of the riverbed and the bank are iterated until a certain iterated hydraulic conductivity coefficient is found that makes the groundwater depth at a subsequent moment the same as the actual groundwater depth. This iterated hydraulic conductivity coefficient is then determined as the hydraulic conductivity coefficient of the simulated river segment.

[0104] S502. Based on the simulated surface water level, groundwater level and riverbed hydraulic conduction coefficient of each grid at the current time, determine the simulated vertical exchange amount of each grid at the current time step.

[0105] For example, the River Subroutine Package (RIV Package) in MODFLOW can be used to simulate the vertical exchange between rivers and groundwater.

[0106] S503. Based on the simulated surface water level, groundwater level and riverbank hydraulic conduction coefficient of each grid at the current time, determine the simulated lateral exchange amount of each grid at the current time step.

[0107] For example, the General Head Boundary (GHB) package in MODFLOW can be used to simulate the lateral exchange between the river and the riverbank.

[0108] In the technical solution provided by this application embodiment, the simulated vertical exchange volume of each grid at the current time step is determined based on the simulated surface water level, groundwater level, and riverbed hydraulic conduction coefficient of the simulated river section in each grid at the current time. Furthermore, the simulated lateral exchange volume of each grid at the current time step is determined based on the simulated surface water level, groundwater level, and riverbank hydraulic conduction coefficient of the simulated river section in each grid. This enables independent calculation of the simulated vertical and lateral exchange volumes, improving the accuracy of the determined simulated vertical and lateral exchange volumes.

[0109] The following explains how the simulated vertical swap amount for each grid cell is determined at the current time step:

[0110] In some embodiments, the simulated vertical exchange amount of each grid at the current time step is determined based on the simulated surface water level, groundwater level, and riverbed hydraulic conduction coefficient of the simulated river segment in each grid at the current time, including: obtaining the riverbed bottom elevation of each grid; and determining the simulated vertical exchange amount of each grid at the current time step based on the simulated surface water level, groundwater level, riverbed bottom elevation, and riverbed hydraulic conduction coefficient of the simulated river segment in each grid at the current time.

[0111] The groundwater level of a grid at the current moment can include the head of groundwater at the same latitude and longitude as that grid at the current moment. The riverbed elevation can be determined based on the river channel design parameters or obtained through measurement.

[0112] In the technical solution provided in this application embodiment, the simulated vertical exchange amount of each grid at the current time step is determined based on the simulated surface water level, groundwater level, riverbed bottom elevation, and riverbed hydraulic conduction coefficient of each grid at the current time, thereby improving the accuracy of the determined simulated vertical exchange amount of each grid at the current time step.

[0113] In some embodiments, the simulated vertical exchange amount of each grid at the current time step is determined based on the simulated surface water level, groundwater level, riverbed bottom elevation of each grid at the current time, and the riverbed hydraulic conduction coefficient of the simulated river segment in each grid. This includes: if the groundwater level of any grid at the current time is less than or equal to the riverbed bottom elevation of any grid at the current time, the simulated vertical exchange amount of any grid at the current time is determined based on the riverbed bottom elevation of any grid at the current time, the simulated surface water level, and the riverbed hydraulic conduction coefficient of the simulated river segment in any grid.

[0114] In some embodiments, the simulated vertical exchange amount of each grid at the current time step is determined based on the simulated surface water level, groundwater level, riverbed bottom elevation of each grid at the current time, and the riverbed hydraulic conduction coefficient of the simulated river segment in each grid. This includes: if the groundwater level of any grid at the current time is greater than the riverbed bottom elevation of any grid at the current time, the simulated vertical exchange amount of any grid at the current time is determined based on the groundwater level of any grid at the current time, the simulated surface water level, and the riverbed hydraulic conduction coefficient of the simulated river segment in any grid.

[0115] In the technical solution provided in this application embodiment, the difference in the relationship between the groundwater level of any grid at the current time and the riverbed bottom elevation of any grid at the current time is used to determine the simulated vertical exchange amount of any grid at the current time step in a differentiated manner. This makes the determined simulated vertical exchange amount of any grid at the current time step consistent with the actual situation and improves the accuracy of the determined simulated vertical exchange amount of any grid at the current time step.

[0116] For example, during the simulation, irregular river channel cross-sections in nature can be generalized as rectangles, and riverbed and bank sediments can be generalized as boundaries with specific hydraulic parameters (e.g., schematic diagrams of the generalized river, groundwater, and river channel cross-sections are shown in the figure). Figure 2 As shown in the figure, the permeability coefficient reflects the water-conducting capacity. For example, the riverbed is generalized as a horizontal boundary controlling vertical exchange, while the banks on both sides are generalized as vertical boundaries controlling horizontal exchange.

[0117] The RIV package in MODFLOW can be used to simulate the vertical exchange between river and groundwater. The amount and direction of the exchange depend on the hydraulic gradient between the river and groundwater levels and the overall hydraulic conductivity of the riverbed (i.e., the riverbed hydraulic conductivity). The riverbed hydraulic conductivity coefficient integrates the hydraulic conductivity and geometric dimensions of the riverbed and serves as the direct basis for calculating exchange fluxes in the model. To achieve dynamic interactive simulation, the river level in the model is obtained through dynamic coupling with the CWatM hydrological model. The river discharge (Q) generated by the CWatM simulation is converted into the corresponding river level in real time at each time step and passed to MODFLOW as a boundary condition.

[0118] like ,but ;like ,but Where h represents the groundwater level of each grid cell at the current moment. This represents the riverbed bottom elevation for each grid cell at the current moment. This represents the simulated vertical exchange amount for each grid cell at the current time step; a positive value indicates river seepage into groundwater. To simulate the hydraulic conductivity coefficient of the riverbed in a river section, This represents the simulated surface water level for each grid cell at the current moment.

[0119] In some embodiments, . This represents the simulated lateral exchange amount for each grid cell at the current time step; a positive value indicates river seepage into the riverbank. To simulate the hydraulic conduction coefficient of the riverbank in a river section, The head assigned to the boundary conditions (here, the simulated surface water level of each grid at the current moment) ).

[0120] For example, and The summation yields the surface water-groundwater exchange volume for each grid at the current time step. : .

[0121] Figure 6 This is a flowchart illustrating the method for determining the simulated river flow of each grid at the next time step in some embodiments, such as... Figure 6 As shown, the method includes:

[0122] S601. Obtain the surface water input of each grid at the current time step, including surface runoff input, water return flow, river surface evaporation, and river water intake.

[0123] For example, the surface runoff input of each grid at the current time step. It can be determined in the following ways: ;in The three summation symbols represent the summation of runoff for different land use / cover types (used to summarize the runoff contribution of each land cover type within the same grid), the summation of different runoff components (used to sum the runoff volumes of different runoff components under the same land type to obtain the total runoff for that land type), and the summation of runoff generated within a certain number of past time steps by time weight. Weighted summation. The total duration of the confluence process (or the number of time steps at which the peak occurs) is determined based on the surface slope and land cover. This represents the runoff volume generated by a land cover class for a runoff component at a specific point in the past. It is the time step of the simulation. The time delay is expressed in days; for example, when =1, that is ;when =2 is That is, the runoff generated yesterday; when =3 is This refers to the runoff generated the day before yesterday.

[0124] For example, surface runoff input refers to the total amount of surface water generated by events such as rainfall and snowmelt in the catchment area around the simulated river section within a time step, and which flows into the river through slope flow or near-surface runoff.

[0125] Water return flow (flow) refers to the portion of water that, within a given time step, is not completely consumed after the completion of water-consuming processes such as industrial and agricultural production and urban and rural residents' daily lives, and ultimately returns to the river channel through drainage systems or seepage.

[0126] River surface evaporation (flow) refers to the volume of water lost directly into the atmosphere from the free water surface (i.e., river surface) of a simulated river section within a time step due to physical evaporation. Its value is usually determined by the evaporation capacity of the water surface and the surface area of ​​the river.

[0127] River water intake (flow) refers to the total amount of water extracted from a river within a given time step through artificial water diversion facilities (such as pumping stations and sluice gates) or natural means to meet various water needs (such as irrigation and water supply).

[0128] S602. Based on the simulated river flow of each grid at the current time, the surface water-groundwater exchange, surface runoff input, water return flow, river surface evaporation, and river water intake of each grid at the current time step, determine the simulated river flow of each grid at the next time step.

[0129] In some embodiments, the simulated surface runoff input, surface water-groundwater exchange, surface runoff input, water return flow, river surface evaporation, and river water intake can refer to flow rate.

[0130] For example, the simulated river flow of each grid at the current time step, the simulated surface runoff input of each grid at the current time step, and the water return flow are summed, and then the surface water-groundwater exchange volume, the river surface evaporation volume, and the river water intake volume are subtracted to obtain the simulated river flow of each grid at the next time step.

[0131] Figure 7 This is a flowchart illustrating the method for determining the river water storage at the next moment in some embodiments, such as... Figure 7 As shown, the method includes:

[0132] S701. Obtain the surface water input of each grid at the current time step, including surface runoff input, water return flow, river surface evaporation, and river water intake, and obtain the river water storage at the current moment.

[0133] The current river water storage refers to the total volume of water stored in the river channel at the current specific point in time in the simulated watershed. It is a state variable that characterizes the instantaneous water storage state of the watershed and is the basis for calculating the river water balance.

[0134] S702. Based on the surface runoff input, water return flow, river surface evaporation, river water intake, surface water-groundwater exchange, and river water storage at the current time step for each grid, determine the river water storage at the next time step.

[0135] For example, river water storage refers to the total amount of water in a river that is either still or flowing at a specific moment.

[0136] In some embodiments, S702 may include the following steps: converting surface water-groundwater exchange volume into water exchange volume; converting the simulated surface runoff input, water return flow, river surface evaporation, and river water intake of each grid at the current time step into water volume; and determining the river water storage at the next time step based on the simulated river segment water volume, water exchange volume, river surface evaporation, and river water intake converted by each grid at the current time step, and the river water storage at the current time.

[0137] In obtaining Subsequently, in CWatM, surface water seepage recharge to groundwater should be considered a discharge term (negative value), while groundwater recharge to surface water should be considered a recharge term (positive value). This is exactly the opposite of MODFLOW. Therefore, [the following is a possible interpretation:] ; ;in, Indicates the duration of a time step; For the next moment, the river's water storage capacity. The current water level in the river channel; Input the simulated surface runoff for each grid at the current time step; The water return flow rate for each grid at the current time step; This refers to the amount of water exchanged. This represents the evaporation rate of the river surface at the current time step for each grid cell. This represents the water intake of each grid at the current time step.

[0138] Figure 8 A flowchart illustrating the simulated exchange process between surface water and groundwater is provided for some embodiments, such as Figure 8 As shown, precipitation / snowmelt can replenish groundwater, and the runoff from precipitation / snowmelt can also flow into surface water through surface runoff or interflow. Surface water can seep into groundwater through rivers, and groundwater can reach surface water through discharge. The sufficiency of surface water resources is determined based on water demands for domestic, industrial, and agricultural use. If sufficient, water is extracted from surface water, which consumes surface water. If surface water is insufficient, groundwater is extracted, which consumes groundwater.

[0139] Based on the same inventive concept, this application also provides a surface water-groundwater dynamic exchange simulation device for implementing the above-mentioned surface water-groundwater dynamic exchange simulation method. The solution provided by this device is similar to the solution described in the above method. Therefore, the specific limitations of one or more surface water-groundwater dynamic exchange simulation device embodiments provided below can be found in the limitations of the surface water-groundwater dynamic exchange simulation method described above, and will not be repeated here.

[0140] In one exemplary embodiment, Figure 9 A schematic diagram of the structure of a surface water-groundwater dynamic exchange simulation device provided in some embodiments, such as... Figure 9 As shown, the surface water-groundwater dynamic exchange simulation device 900 includes:

[0141] The surface water level determination module 901 is used to determine the simulated surface water level of each grid at the current time based on the simulated river flow of each grid obtained by dividing the simulated watershed into latitude and longitude grids.

[0142] The vertical and lateral exchange quantity determination module 902 is used to determine the simulated vertical exchange quantity and simulated lateral exchange quantity of each grid at the current time step based on the simulated surface water level and the obtained groundwater level of each grid at the current time; the current time step is the time step between the current time and the next time step.

[0143] The surface water-groundwater exchange volume determination module 903 is used to determine the surface water-groundwater exchange volume of each grid at the current time step based on the simulated vertical exchange volume and simulated lateral exchange volume of each grid at the current time step.

[0144] The river flow determination module 904 is used to determine the simulated river flow of each grid at the next time step based on the simulated river flow of each grid at the current time step, the surface water-groundwater exchange volume of each grid at the current time step, and the surface water input volume, and to take the next time step as the new current time step; wherein, each module is used to cyclically execute the corresponding steps.

[0145] In some embodiments, the surface water level determination module 901 includes an acquisition unit and a determination unit; the acquisition unit is used to acquire the roughness, size and slope of the river channel covered by each grid; the determination unit is used to determine the simulated surface water level of each grid at the current time based on the simulated river flow of each grid at the current time, and the roughness, size and slope of the river channel covered by each grid.

[0146] In some embodiments, the vertical and lateral exchange quantity determination module 902 includes an acquisition unit, a vertical exchange quantity determination unit, and a lateral exchange quantity determination unit; the acquisition unit is used to acquire the riverbed hydraulic conductivity coefficient and the bank hydraulic conductivity coefficient of the simulated river segment in each grid; the vertical exchange quantity determination unit is used to determine the simulated vertical exchange quantity of each grid at the current time step based on the simulated surface water level, groundwater level, and riverbed hydraulic conductivity coefficient of the simulated river segment in each grid at the current time; the lateral exchange quantity determination unit is used to determine the simulated lateral exchange quantity of each grid at the current time step based on the simulated surface water level, groundwater level, and bank hydraulic conductivity coefficient of the simulated river segment in each grid at the current time.

[0147] In some embodiments, the vertical exchange quantity determination unit is further configured to obtain the riverbed bottom elevation of each grid; and determine the simulated vertical exchange quantity of each grid at the current time step based on the simulated surface water level, groundwater level, riverbed bottom elevation of each grid at the current time and the riverbed hydraulic conduction coefficient of the simulated river segment in each grid.

[0148] In some embodiments, the vertical exchange quantity determination unit is further configured to determine the simulated vertical exchange quantity of any grid at the current time step, based on the riverbed bottom elevation of any grid at the current time, the simulated surface water level, and the riverbed hydraulic conductivity coefficient of the simulated river segment in any grid, when the groundwater level of any grid at the current time is less than or equal to the riverbed bottom elevation of any grid at the current time.

[0149] In some embodiments, the vertical exchange quantity determination unit is further configured to determine the simulated vertical exchange quantity of any grid at the current time step based on the groundwater level of any grid at the current time, the simulated surface water level, and the hydraulic conductivity coefficient of the simulated river section in any grid when the groundwater level of any grid at the current time is greater than the riverbed bottom elevation of any grid at the current time.

[0150] In some embodiments, the river flow determination module 904 includes an acquisition unit and a determination unit. The acquisition unit is used to acquire the surface water input of each grid at the current time step, including surface runoff input, water return flow, river surface evaporation, and river water intake. The determination unit is used to determine the simulated river flow of each grid at the next time step based on the simulated river flow of each grid at the current time step, the surface water-groundwater exchange of each grid at the current time step, the surface runoff input, water return flow, river surface evaporation, and river water intake.

[0151] In some embodiments, the surface water-groundwater dynamic exchange simulation device 900 further includes an acquisition module and a water storage determination module; the acquisition module is used to acquire the surface water input of each grid at the current time step, including surface runoff input, water return flow, river surface evaporation, and river water intake, and to acquire the river water storage at the current time; the water storage determination module is used to determine the river water storage at the next time step based on the surface runoff input, water return flow, river surface evaporation, river water intake, surface water-groundwater exchange volume, and the river water storage at the current time step.

[0152] The descriptions of the above device embodiments are similar to those of the above method embodiments, and have similar beneficial effects. For technical details not disclosed in the device embodiments of this application, please refer to the descriptions of the method embodiments of this application for understanding.

[0153] The modules in the aforementioned dynamic exchange simulation device between surface water and groundwater can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device, or stored in the memory of a computer device as software, so that the processor can call and execute the corresponding operations of each module.

[0154] In one exemplary embodiment, Figure 10This is a schematic diagram of the structure of a computer device provided in some embodiments. The computer device includes a processor, memory, input / output interface, communication interface, display unit, and input device. The processor, memory, and input / output interface are connected via a system bus, and the communication interface, display unit, and input device are also connected to the system bus via the input / output interface. The processor of this computer device provides computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The input / output interface of the computer device is used for exchanging information between the processor and external devices. The communication interface of the computer device is used for wired or wireless communication with external terminals. Wireless communication can be achieved through Wireless Fidelity (WIFI), mobile cellular networks, Near Field Communication (NFC), or other technologies. When the computer program is executed by the processor, it implements a dynamic exchange simulation method for surface water and groundwater. The display unit of the computer device is used to form a visually visible image and can be a display screen, a projection device, or a virtual reality imaging device. The display screen can be an LCD screen or an e-ink screen. The input device of the computer device can be a touch layer covering the display screen, or buttons, trackballs, or touchpads set on the casing of the computer device, or external keyboards, touchpads, or mice, etc.

[0155] Those skilled in the art will understand that Figure 10 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0156] For example, a computer device includes a memory and a processor, the memory storing a computer program, and the processor executing the computer program to implement the steps of the method of any of the above embodiments.

[0157] In one embodiment, a computer-readable storage medium is provided, wherein a computer program, when executed by a processor, implements the steps of the method provided in any of the above embodiments.

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

[0159] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the methods described above.

[0160] The processor, functional modules, or functional units in any embodiment of this application may include an integration of one or more of the following: a general-purpose processor, an application-specific integrated circuit (ASIC), a digital signal processor (DSP), a digital signal processing device (DSPD), a programmable logic device (PLD), a field-programmable gate array (FPGA), a central processing unit (CPU), a graphics processing unit (GPU), an embedded neural network processing unit (NPU), a controller, a microcontroller, a microprocessor, a programmable logic device, a discrete gate or transistor logic device, a discrete hardware component, a quantum computing-based data processing logic unit, an artificial intelligence (AI) processor, etc. The general-purpose processor may be a microprocessor or any conventional processor.

[0161] The memory or computer-readable storage medium in any embodiment of this application may include at least one of non-volatile memory and volatile memory. Non-volatile memory includes integration of one or more of the following: Read Only Memory (ROM), Programmable Read-Only Memory (PROM), Erasable Programmable Read-Only Memory (EPROM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Ferromagnetic Random Access Memory (FRAM), Flash Memory, Magnetic Surface Memory, Optical Disc, Compact Disc Read-Only Memory (CD-ROM), Magnetic Tape, Floppy Disk, Flash Memory, Optical Memory, High-Density Embedded Non-Volatile Memory, Resistive Random Access Memory (ReRAM), Magnetoresistive Random Access Memory (MRAM), Ferroelectric Random Access Memory (FRAM), Phase Change Memory (PCM), Graphene Memory, Volatile Memory, etc. Volatile memory includes one or more of the following: Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0162] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this application.

[0163] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A method for simulating dynamic exchange between surface water and groundwater, characterized in that, The method includes: Based on the simulated river flow of each grid obtained by dividing the simulated watershed into latitude and longitude grids at the current time, the simulated surface water level of each grid at the current time is determined. Obtain the riverbed hydraulic conduction coefficient and the bank hydraulic conduction coefficient of the simulated river segment in each grid. Based on the simulated surface water level, the obtained groundwater level, and the hydraulic conductivity coefficient of the simulated river section in each grid at the current time, the simulated vertical exchange amount of each grid at the current time step is determined. Based on the simulated surface water level, groundwater level, and bank hydraulic conduction coefficient of the simulated river segment in each grid at the current time, the simulated lateral exchange amount of each grid at the current time step is determined; the current time step is the time step between the current time and the next time step. Based on the simulated vertical and lateral exchange amounts of each grid at the current time step, determine the surface water-groundwater exchange amount of each grid at the current time step. Based on the simulated river flow of each grid at the current time, the surface water-groundwater exchange and surface water input of each grid at the current time step, determine the simulated river flow of each grid at the next time step, and use the next time step as the new current time step, repeating the above steps.

2. The method according to claim 1, characterized in that, The step of determining the simulated surface water level of each grid at the current time based on the simulated river flow rate of each grid obtained by dividing the simulated watershed into latitude and longitude grids includes: Obtain the roughness, size, and slope of the river channel covered by each of the grids; The simulated surface water level of each grid at the current moment is determined based on the simulated river flow of each grid at the current moment, as well as the roughness, size, and slope of the river covered by each grid.

3. The method according to claim 1, characterized in that, The process of obtaining the riverbed hydraulic conductivity coefficient and the bank hydraulic conductivity coefficient of the simulated river segment in each of the grids includes: Based on the sediment permeability coefficient, river length, river width, and riverbed thickness of the simulated river segment in each grid, the hydraulic conductivity coefficient of the simulated river segment in each grid is determined. Based on the sediment permeability coefficient of the simulated river section in each grid, the river length of the simulated river section, the longitudinal thickness of the simulated river section, and the transverse thickness of the simulated river section, the hydraulic conductivity coefficient of the simulated river section in each grid is determined.

4. The method according to any one of claims 1-3, characterized in that, The step of determining the simulated vertical exchange amount of each grid at the current time step based on the simulated surface water level, the obtained groundwater level, and the hydraulic conductivity coefficient of the simulated river section in each grid includes: Obtain the riverbed bottom elevation for each of the aforementioned grids; Based on the simulated surface water level, groundwater level, riverbed bottom elevation, and riverbed hydraulic conduction coefficient of the simulated river segment in each grid at the current time, the simulated vertical exchange amount of each grid at the current time step is determined.

5. The method according to claim 4, characterized in that, The step of determining the simulated vertical exchange amount of each grid at the current time step based on the simulated surface water level, groundwater level, riverbed bottom elevation, and riverbed hydraulic conduction coefficient of the simulated river segment in each grid at the current time includes: If the groundwater level of any grid at the current time is less than or equal to the riverbed bottom elevation of any grid at the current time, the simulated vertical exchange amount of any grid at the current time is determined based on the riverbed bottom elevation of any grid at the current time, the simulated surface water level, and the riverbed hydraulic conduction coefficient of the simulated river segment in any grid. If the groundwater level of any grid at the current time is greater than the riverbed bottom elevation of any grid at the current time, the simulated vertical exchange amount of any grid at the current time is determined based on the groundwater level of any grid at the current time, the simulated surface water level, and the hydraulic conductivity coefficient of the simulated river section in any grid.

6. The method according to any one of claims 1-3, characterized in that, The step of determining the simulated river flow of each grid at the next time step based on the simulated river flow of each grid at the current time step, the surface water-groundwater exchange volume of each grid at the current time step, and the surface water input volume includes: The surface water input of each grid at the current time step includes surface runoff input, water return flow, river surface evaporation, and river water intake. Based on the simulated river flow of each grid at the current time, the surface water-groundwater exchange, surface runoff input, water return flow, river surface evaporation, and river water intake of each grid at the current time step, determine the simulated river flow of each grid at the next time step.

7. The method according to any one of claims 1-3, characterized in that, The method further includes: The surface water input of each grid at the current time step includes surface runoff input, water return flow, river surface evaporation, and river water intake, as well as the river water storage at the current moment. Based on the surface runoff input, water return flow, river surface evaporation, river water intake, surface water-groundwater exchange, and river water storage at the current time step for each grid, determine the river water storage at the next time step.

8. A dynamic exchange simulation device for surface water and groundwater, characterized in that, The device includes: The surface water level determination module is used to determine the simulated surface water level of each grid at the current time based on the simulated river flow of each grid obtained by dividing the simulated watershed into latitude and longitude grids. The vertical and lateral exchange quantity determination module is used to obtain the riverbed hydraulic conductivity coefficient and the bank hydraulic conductivity coefficient of the simulated river segment in each grid; determine the simulated vertical exchange quantity of each grid at the current time step based on the simulated surface water level, the obtained groundwater level, and the riverbed hydraulic conductivity coefficient of the simulated river segment in each grid; and determine the simulated lateral exchange quantity of each grid at the current time step based on the simulated surface water level, the groundwater level, and the bank hydraulic conductivity coefficient of the simulated river segment in each grid; wherein the current time step is the time step between the current time and the next time step; The surface water-groundwater exchange volume determination module is used to determine the surface water-groundwater exchange volume of each grid at the current time step based on the simulated vertical exchange volume and simulated lateral exchange volume of each grid at the current time step. The river flow determination module is used to determine the simulated river flow of each grid at the next time step based on the simulated river flow of each grid at the current time step, the surface water-groundwater exchange volume of each grid at the current time step, and the surface water input volume, and to take the next time step as the new current time step; wherein, each module is used to cyclically execute the corresponding steps.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 7.

Images