Surface water and groundwater simulation coupling method, device, computer equipment, readable storage medium and program product

Abstract

The present application relates to a surface water and groundwater simulation coupling method, device, computer equipment, readable storage medium and program product. The method includes: in a surface water simulation model, determining a target water allocation parameter based on environmental data of a target area; in a groundwater simulation model, determining a groundwater flow function based on a finite difference simulation algorithm and target water allocation parameters, and determining a groundwater head based on the groundwater flow function; and determining a target groundwater discharge based on the groundwater head; in a surface water simulation model, determining a confluence parameter corresponding to a target area based on a preset fluid dynamics algorithm and target groundwater discharge, and the confluence parameter characterizes the degree of mutual influence between surface water and groundwater. The use of this method can improve the coupling accuracy of surface water and groundwater.

Description

Technical Field

[0001] The present application relates to the field of hydrological technology, and in particular to a surface water and groundwater simulation coupling method, device, computer equipment, readable storage medium and program product. Background Art

[0002] Under the influence of factors such as climate fluctuations, extreme weather events and excessive groundwater exploitation, the interaction and transformation relationship between surface water and groundwater has changed.

[0003] In related technologies, surface water and groundwater are studied as single resources at a smaller scale. The corresponding models used include global hydrological models, land surface models, groundwater models, etc. The above models provide feasible tools for assessing the impact of human activities on groundwater reserves. However, there is a lack of simulation of water use processes under strong human activities and the interaction process between surface water and groundwater. Therefore, the coupling accuracy of surface water and groundwater is low. Summary of the invention

[0004] Based on this, it is necessary to provide a surface water and groundwater simulation coupling method, device, computer equipment, readable storage medium and program product to address the above-mentioned technical problems, which can improve the coupling accuracy of surface water and groundwater based on more comprehensive data.

[0005] In a first aspect, the present application provides a surface water and groundwater simulation coupling method, comprising:

[0006] In the surface water simulation model, the target water allocation parameters are determined based on the environmental data of the target area;

[0007] In a groundwater simulation model, a groundwater flow function is determined based on a finite difference simulation algorithm and the target water volume allocation parameter, and a groundwater head is determined based on the groundwater flow function; and a target groundwater discharge is determined based on the groundwater head;

[0008] In the surface water simulation model, based on a preset fluid dynamics algorithm and the target groundwater discharge, a confluence parameter corresponding to the target area is determined, and the confluence parameter characterizes the degree of mutual influence between the surface water and the groundwater.

[0009] In one embodiment, the environmental data at least includes: total water demand, soil unsaturated hydraulic conductivity, soil moisture content and soil residual moisture content. In the surface water simulation model, determining the target water allocation parameters based on the environmental data includes:

[0010] In the surface water simulation model, the total water demand is divided according to the population density distribution strategy to obtain the water demand;

[0011] The product of the unsaturated hydraulic conductivity of the soil and the unit time length is determined as the water transfer amount in the soil; the difference between the soil water content and the residual water content of the soil is determined as the effective use water content of the soil;

[0012] Determine the minimum value of the water transfer amount and the effective utilization water content of the soil, and determine the product of the minimum value and the area ratio of the unsaturation of the underground aquifer as the deep leakage amount;

[0013] Determine the preferential flow based on the product value of the amount of water reaching the soil surface after being intercepted by the vegetation canopy and the relative saturation of the soil, and determine the deep infiltration amount and the preferential flow as the infiltration amount;

[0014] The water demand and the infiltration amount are determined as initial water allocation parameters, and the target water allocation parameters are determined based on the initial water allocation parameters.

[0015] In one embodiment, determining a target groundwater discharge based on the groundwater head comprises:

[0016] If the groundwater head is higher than the spatial range of the lower boundary of the soil corresponding to the target area, the difference between the groundwater head and the spatial range of the lower boundary of the soil is determined, and the product of the difference and the hydraulic conductivity of the groundwater aquifer is determined as the total discharge volume of groundwater;

[0017] The product of the total amount of discharged water and the proportion of the water area corresponding to the target area is determined as the base flow generated by the recharge from groundwater to surface water, and the difference between the total amount of discharged water and the base flow is determined as the capillary rise;

[0018] The base flow and the capillary rise are determined as an initial groundwater discharge, and a target groundwater discharge is determined based on the initial groundwater discharge.

[0019] In one embodiment, the target area includes a plurality of first grids obtained by division based on a first grid division scale, and a plurality of second grids obtained by division based on a second grid division scale; and determining the target water allocation parameter based on the initial water allocation parameter includes:

[0020] For each second grid in the groundwater simulation model, determine each first grid in the surface water simulation model corresponding to the second grid; determine a first overlapping area between the second grid and the corresponding first grids, and perform weighted summation on the product value of the first overlapping area and the initial water allocation parameter in the first grid to obtain a first weighted summation value; and determine the ratio of the first weighted summation value to the area of ​​the second grid as the target water allocation value of the second grid.

[0021] In one embodiment, determining a target groundwater discharge based on the initial groundwater discharge comprises:

[0022] For each of the first grids in the surface water simulation model, determine each of the second grids in the groundwater simulation model corresponding to the first grid; determine a second overlapping area between the first grid and each of the corresponding second grids, and perform weighted summation on the product of the second overlapping area and the initial groundwater discharge in the second grid to obtain a second weighted summation value; and determine the ratio of the second weighted summation value to the area of ​​the first network as the target groundwater discharge of the first grid.

[0023] In one embodiment, in the surface water simulation model, determining the confluence parameters based on a preset fluid dynamics algorithm and the target groundwater discharge volume includes:

[0024] In the surface water simulation model, the ratio of the flow rate of surface water to the coordinate change value of surface water is determined as the flow rate change rate of surface water, and the ratio of the cross-sectional area of ​​surface water to the unit time length is determined as the cross-sectional change rate of surface water;

[0025] A sink flow parameter is determined based on the flow rate change rate, the cross-sectional change rate, the capillary rise, and the base flow.

[0026] In one embodiment, the method further comprises:

[0027] Determining an initial surface water simulation model and an initial groundwater simulation model based on the environmental data and initial model parameters;

[0028] In the initial surface water simulation model and the initial groundwater simulation model corresponding to the first spatial resolution, calibrating the initial model parameters based on a preset calibration algorithm to obtain first model parameters;

[0029] In the initial surface water simulation model and the initial groundwater simulation model corresponding to the second spatial resolution, calibrating the first model parameters based on the preset calibration algorithm to obtain second model parameters;

[0030] The initial surface water simulation model and the initial groundwater simulation model are optimized based on the second model parameters to obtain the surface water simulation model and the groundwater simulation model.

[0031] In a second aspect, the present application also provides a surface water and groundwater simulation coupling device, comprising:

[0032] A first determination module is used to determine target water allocation parameters based on environmental data of a target area in a surface water simulation model;

[0033] A second determination module is used to determine a groundwater flow function in a groundwater simulation model based on a finite difference simulation algorithm and the target water volume allocation parameter, and to determine a groundwater head based on the groundwater flow function; and to determine a target groundwater discharge based on the groundwater head;

[0034] The third determination module is used to determine the confluence parameters corresponding to the target area in the surface water simulation model based on a preset fluid dynamics algorithm and the target groundwater discharge, wherein the confluence parameters represent the degree of mutual influence between the surface water and the groundwater.

[0035] In a third aspect, the present application further provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the following steps are implemented:

[0036] In the surface water simulation model, the target water allocation parameters are determined based on the environmental data of the target area;

[0037] In a groundwater simulation model, a groundwater flow function is determined based on a finite difference simulation algorithm and the target water volume allocation parameter, and a groundwater head is determined based on the groundwater flow function; and a target groundwater discharge is determined based on the groundwater head;

[0038] In the surface water simulation model, based on a preset fluid dynamics algorithm and the target groundwater discharge, a confluence parameter corresponding to the target area is determined, and the confluence parameter characterizes the degree of mutual influence between the surface water and the groundwater.

[0039] In a fourth aspect, the present application further provides a computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor, the following steps are implemented:

[0040] In the surface water simulation model, the target water allocation parameters are determined based on the environmental data of the target area;

[0041] In a groundwater simulation model, a groundwater flow function is determined based on a finite difference simulation algorithm and the target water volume allocation parameter, and a groundwater head is determined based on the groundwater flow function; and a target groundwater discharge is determined based on the groundwater head;

[0042] In the surface water simulation model, based on a preset fluid dynamics algorithm and the target groundwater discharge, a confluence parameter corresponding to the target area is determined, and the confluence parameter characterizes the degree of mutual influence between the surface water and the groundwater.

[0043] In a fifth aspect, the present application further provides a computer program product, including a computer program, which implements the following steps when executed by a processor:

[0044] In the surface water simulation model, the target water allocation parameters are determined based on the environmental data of the target area;

[0045] In a groundwater simulation model, a groundwater flow function is determined based on a finite difference simulation algorithm and the target water volume allocation parameter, and a groundwater head is determined based on the groundwater flow function; and a target groundwater discharge is determined based on the groundwater head;

[0046] In the surface water simulation model, based on a preset fluid dynamics algorithm and the target groundwater discharge, a confluence parameter corresponding to the target area is determined, and the confluence parameter characterizes the degree of mutual influence between the surface water and the groundwater.

[0047] The above-mentioned surface water and groundwater simulation coupling method, device, computer equipment, readable storage medium and program product can effectively simulate the interaction process between surface water and groundwater in areas with strong human activities, thereby improving the simulation accuracy of surface water and groundwater, by determining the target water allocation parameters of the target area in the surface water simulation model, calculating the groundwater head in the groundwater simulation model, determining the target groundwater discharge based on the groundwater head, and determining the evolution of surface water in the surface water simulation model based on a preset fluid dynamics algorithm and the target groundwater discharge. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] In order to more clearly illustrate the technical solutions in the embodiments of the present application or related technologies, the drawings required for use in the embodiments of the present application or related technical descriptions will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0049] FIG1 is a diagram showing an application environment of a surface water and groundwater simulation coupling method in one embodiment;

[0050] FIG2 is a schematic diagram of the coupling of a surface water simulation model and a groundwater simulation model in one embodiment;

[0051] FIG3 is a comparison diagram of the simulation and measured results of groundwater storage change anomaly in a certain area from XX to XX in one embodiment;

[0052] FIG4 is a structural block diagram of a surface water and groundwater simulation coupling device in one embodiment;

[0053] FIG. 5 is a diagram showing the internal structure of a computer device in one embodiment. DETAILED DESCRIPTION

[0054] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0055] In an exemplary embodiment, as shown in FIG1 , a surface water and groundwater simulation coupling method is provided. This embodiment uses the method applied to a terminal as an example for illustration. It is understandable that the method can also be applied to a server, and can also be applied to a system including a terminal and a server, and is implemented through the interaction between the terminal and the server. In this embodiment, the method includes the following steps:

[0056] Step 101: In a surface water simulation model, target water allocation parameters are determined based on environmental data of a target area.

[0057] The surface water simulation model can be a global surface water model (Community Water Model, CWatM), and the target area can be a specific geographical area that needs water resources management and planning, such as the Haihe River Basin or other areas of high-intensity human activities. Environmental data may include: meteorological data, surface properties, surface cover types, soil, vegetation, lakes and reservoirs, water demand of various departments and confluence related parameters, etc. The target water allocation parameter can be a water allocation parameter transmitted from the surface water simulation model to the groundwater simulation model, and the water allocation parameter can characterize the water demand of surface water and the infiltration of surface water into groundwater.

[0058] Specifically, in the surface water simulation model, the terminal can obtain the environmental data of each unit time corresponding to the target area, and calculate the environmental data to obtain the water allocation parameter, where the unit time can be one day. The terminal can perform parameter conversion on the water allocation parameter output by the surface water simulation model to obtain the target water allocation parameter, and process the target allocation parameter through the groundwater simulation model.

[0059] Step 102: In the groundwater simulation model, based on the finite difference simulation algorithm and the target water allocation parameters, determine the groundwater flow function, and determine the groundwater head based on the groundwater flow function; and determine the target groundwater discharge based on the groundwater head.

[0060] The groundwater simulation model may be a modular finite-difference groundwater flow model (MODFLOW), and the groundwater flow function is used to simulate the two-dimensional flow of groundwater. The groundwater head may be used to characterize the vertical distance of groundwater below the surface. The target groundwater discharge may be the groundwater discharge transferred from the groundwater simulation model to the surface water simulation model, and the groundwater discharge is used to characterize the recharge, flow or discharge from the groundwater system to the surface water system.

[0061] Specifically, in the groundwater simulation model, the terminal can simulate the groundwater flow through the finite difference simulation algorithm, and can construct the groundwater flow function based on the target water distribution parameters, water supply degree, and hydraulic conductivity in the horizontal and vertical directions. The groundwater flow function is used to simulate the two-dimensional flow of groundwater. The specific expression of the groundwater flow function can be:

[0062]

[0063] Among them, S y is the water supply degree, h is the groundwater head, K xx With K yy is the hydraulic conductivity in the x and y directions, W and F are the target water allocation parameters, W is the groundwater demand, and F is the infiltration rate.

[0064] In the surface water simulation model, the vertical movement of groundwater is simulated by aquifer stratification. In this embodiment, multiple layers of unconfined aquifers can be set, and the water reserves of the aquifers can be calculated according to the thickness, hydraulic conductivity, porosity and water storage rate of the aquifers. The water reserves can be determined according to the water supply degree and water storage rate. Specifically, the calculation expression can be:

[0065]

[0066]

[0067]

[0068] Where GWS is the aquifer water storage; GWS sy GWS is the water storage component of the aquifer that is proportional to the water supply degree; ss is the water storage component of the aquifer that is proportional to the water storage rate; S y Can be the water supply degree; S s can be the water storage rate, H is the thickness of the aquifer, A is the area of ​​the aquifer, h is the height of the water head relative to the bottom plate, S Fis the saturation obtained by smoothing the ratio of h to H. Smoothing is done to achieve partial derivative calculations within the model. The terminal can select whether the water storage in the aquifer requires GWS based on the user's configuration. sy If the user sets the aquifer as unconvertible, GWS sy For example, in this embodiment, two unconfined aquifers can be set, with thicknesses of 50 m and 150 m, respectively, and hydraulic conductivity of 1×10 -7 m / s and 1×10 -5 The porosity and water storage capacity of the two aquifers are 0.08 and 2.1×10 -6 m -1 In this example, both aquifers are considered GWS sy and GWS ss impact.

[0069] In the groundwater simulation model, the terminal can calculate the groundwater head by iteratively solving the groundwater flow function. The terminal can determine the specific process of groundwater action to surface water through the relationship between the groundwater head and the spatial range of the lower boundary of the underlying soil. For example, when the groundwater head is higher than the spatial range of the lower boundary of the underlying soil, natural discharge of groundwater to surface water will occur, and the groundwater discharge amount is determined based on the natural discharge of groundwater to surface water. When the terminal transfers the groundwater discharge amount from the groundwater simulation model to the surface water simulation model, it is necessary to convert the groundwater discharge amount into the target groundwater discharge amount.

[0070] Step 103: In the surface water simulation model, based on a preset fluid dynamics algorithm and a target groundwater discharge, determine the confluence parameters corresponding to the target area.

[0071] The preset fluid dynamics algorithm may be a dynamic wave method. The confluence parameter may characterize the degree of mutual influence between surface water and groundwater. The confluence time of surface water is determined according to land cover, slope, etc.

[0072] Specifically, in the surface water simulation model, the terminal can calculate the evolution process of surface water in the target area based on preset fluid dynamics algorithms, river flow, downstream water level, lateral flow velocity, pressure, target groundwater discharge and other parameters, and determine the degree of mutual influence between surface water and groundwater in the target area based on the evolution process.

[0073] The above-mentioned surface water and groundwater simulation coupling method can effectively simulate the interaction process between surface water and groundwater in areas with strong human activities, and improve the simulation accuracy of surface water and groundwater, by determining the target water allocation parameters of the target area in the surface water simulation model, calculating the groundwater head in the groundwater simulation model, determining the target groundwater discharge based on the groundwater head, and determining the evolution of surface water based on the preset fluid dynamics algorithm and the target groundwater discharge in the surface water simulation model.

[0074] In an exemplary embodiment, in a surface water simulation model, determining a target water allocation parameter based on environmental data includes:

[0075] In the surface water simulation model, the total water demand is divided according to the population density distribution strategy to obtain the water demand; the product of the soil unsaturated hydraulic conductivity and the unit time is determined as the water transfer in the soil; the difference between the soil moisture content and the soil residual moisture content is determined as the effective soil water content; the minimum value of the water transfer and the effective soil water content is determined, and the product of the minimum value and the area ratio of the unsaturation of the underground aquifer is determined as the deep infiltration; based on the product of the amount of water reaching the soil surface after being intercepted by the vegetation canopy and the relative saturation of the soil, the preferential flow is determined, and the deep infiltration and preferential flow are determined as the infiltration; the water demand and the infiltration are determined as the initial water allocation parameters, and the target water allocation parameters are determined based on the initial water allocation parameters.

[0076] Among them, environmental data include at least: total water demand, soil unsaturated hydraulic conductivity, soil moisture content and soil residual moisture content. The setting of population density distribution strategy can be based on the downscaling of population density distribution. Deep infiltration can be the process in which soil moisture exceeds the absorption capacity of plant roots and the water holding capacity of the soil during rainfall or irrigation, and then penetrates downward to the soil layer below the root activity layer, and finally replenishes groundwater. Preferential flow can be the process in which surface water bypasses the soil matrix through soil cracks, caves, etc. to quickly replenish groundwater.

[0077] Specifically, in the surface water simulation model, the terminal can downscale the total water demand in the water resources bulletin into a grid map according to the population density distribution according to the water demand of each department. Based on this, the water demand of each department can be allocated to different water sources such as lake water, reservoir water, and groundwater in the surface water simulation model.

[0078] In the surface water simulation model, surface water recharges groundwater, which can include deep seepage and preferential flow. During the simulation of the surface water simulation model, surface water infiltrates layer by layer into the soil at different depths. When the moisture content of the underlying soil is greater than its field capacity, deep seepage occurs to recharge groundwater. The deep seepage is proportional to the unsaturated hydraulic conductivity of the soil and is affected by whether the underground aquifer is saturated. The terminal can determine the product of the unsaturated hydraulic conductivity of the soil and the unit time as the water transfer in the soil; determine the difference between the soil moisture content and the residual moisture content of the soil as the effective use moisture content of the soil; determine the minimum value of the moisture transfer and the effective use moisture content of the soil, and determine the product of the minimum value and the area ratio of the unsaturation of the underground aquifer as the deep seepage.

[0079] For example, in this embodiment, three layers of soil may be provided, and the calculation formula for the deep seepage amount may be:

[0080]

[0081] Where Prec is the deep seepage, K3(θ3) is the unsaturated hydraulic conductivity of the bottom (third layer) soil, θ3 is the volumetric moisture content of the bottom soil, w3 is the moisture content of the bottom soil, and w res3 is the residual moisture content of the underlying soil, and satFrac is the area percentage of the saturated groundwater aquifer in the target area.

[0082] The expression of the saturated area fraction satFrac of the groundwater aquifer in the target area can be:

[0083]

[0084] Among them, Area_sat is the area of ​​the saturated aquifer in the target area. The area of ​​the saturated aquifer is the area where the groundwater head is higher than the lower boundary area of ​​the underlying soil. Area_tot is the total area in the target area.

[0085] The calculation formula of unsaturated hydraulic conductivity of soil can be:

[0086]

[0087] Where K(θ) is the unsaturated hydraulic conductivity of soil with volumetric moisture content θ, K s is the saturated hydraulic conductivity of soil, θ is the volumetric water content of soil, and θ s is the saturated moisture content of the soil, and β is the parameter to be determined. s and θ s is a constant. θ is the simulated soil volume moisture content at the current moment.

[0088] In the surface water simulation model, the preferential flow is determined by the amount of water reaching the soil surface after canopy interception, the relative soil saturation and the calibrated preferential flow correction coefficient. The specific formula for preferential flow can be:

[0089]

[0090] Among them, Pref is the preferential flow, availWaterInfiltration is the amount of water reaching the soil surface after being intercepted by the vegetation canopy, relSat is the relative soil saturation, and cPrefFlow is the preferential flow parameter obtained by calibration.

[0091] The terminal can calculate the sum of the precipitation in the environmental data and the water storage intercepted by the vegetation canopy, and calculate the difference between the sum and the interception capacity of the vegetation canopy in the environmental data. If the difference is greater than 0, the difference is determined as the amount of water reaching the soil surface after being intercepted by the vegetation canopy; if the difference is less than 0, the amount of water reaching the soil surface after being intercepted by the vegetation canopy is 0. The expression of the amount of water reaching the soil surface after being intercepted by the vegetation canopy can be:

[0092]

[0093] Among them, Rain(t) is the precipitation; Intcp Stor (t) is the water storage intercepted by the vegetation canopy; Intcp Cap is the interception capacity of the vegetation canopy. t is the current moment, and t-1 is the moment before the current moment.

[0094] The expression of water storage intercepted by vegetation canopy can be:

[0095]

[0096] Among them, Intcp Envap (t) is the evaporation throttled by the vegetation canopy, which can be obtained by calculating meteorological data using the Penman formula.

[0097] The expression of soil relative saturation can be:

[0098]

[0099] Among them, w1 and w2 are the water contents of the first and second layers of soil, respectively, and w s1 and w s2 are the saturated moisture contents of the first and second soil layers respectively, and the moisture contents can be environmental data.

[0100] The terminal can determine the deep seepage and preferential flow as the infiltration amount, and the water demand and infiltration amount as the initial water allocation parameters. In the process of transmitting data from the surface water simulation model to the groundwater simulation model, the initial water allocation parameters are converted to obtain the target water allocation parameters, and the target water allocation parameters are sent to the groundwater simulation model.

[0101] In an exemplary embodiment, determining a target groundwater discharge based on a groundwater head includes:

[0102] If the groundwater head is higher than the spatial range of the lower boundary of the soil corresponding to the target area, the difference between the groundwater head and the spatial range of the lower boundary of the soil is determined, and the product of the difference and the hydraulic conductivity of the groundwater aquifer is determined as the total groundwater discharge; the product of the total discharge volume and the proportion of the water area corresponding to the target area is determined as the base flow generated by the recharge from groundwater to surface water, and the difference between the total discharge volume and the base flow is determined as the capillary rise; the base flow and the capillary rise are determined as the initial groundwater discharge, and the target groundwater discharge is determined based on the initial groundwater discharge.

[0103] Specifically, the terminal can compare the relationship between the groundwater head and the spatial range of the lower boundary of the bottom soil corresponding to the target area. If the groundwater head is lower than or equal to the spatial range of the lower boundary of the bottom soil corresponding to the target area, the groundwater has no effect on the surface water, that is, the total natural discharge of groundwater is 0. If the groundwater head is higher than the spatial range of the lower boundary of the bottom corresponding to the target area, the difference between the groundwater head and the spatial range of the lower boundary of the soil is calculated, and the product of the difference and the hydraulic conductivity of the groundwater aquifer is determined as the total discharge of groundwater. The expression of the total discharge of groundwater can be:

[0104]

[0105] Where Dis is the total amount of discharged water, h is the water head of the first aquifer simulated by MODFLOW, h0 is the height range of the lower boundary of the soil, K GW1 is the hydraulic conductivity of the first aquifer.

[0106] The terminal can determine the product of the total amount of discharged water and the proportion of the water area corresponding to the target area as the base flow generated by the recharge from groundwater to surface water, and determine the difference between the total amount of discharged water and the base flow as the capillary rise. The specific expression of the base flow can be:

[0107]

[0108] Among them, Baseflow is the base flow generated by groundwater recharge to the surface, ratio SW It is the proportion of water area in the target area.

[0109] The expression for capillary rise can be:

[0110]

[0111] Among them, CapRise GW It is the capillary rise of groundwater into the underlying soil.

[0112] The terminal can determine the base flow and capillary rise as the initial groundwater discharge. When transferring the initial groundwater discharge from the groundwater simulation model to the surface water simulation model, it is necessary to convert the initial groundwater discharge into the target groundwater discharge.

[0113] In addition, when calculating the capillary rise in this embodiment, in order to ensure that the capillary rise does not exceed the water content of the lower soil, the capillary rise can be calculated by the following calculation method, and the capillary rise is corrected to ensure that the capillary rise does not exceed the difference between the water content of the lower soil and the retention water content. The formula is as follows:

[0114]

[0115] Among them, CapRise_correct is the corrected capillary rise, θ2 is the volumetric moisture content of the lower soil, and θ r2 It is the moisture content of the lower soil.

[0116] In an exemplary embodiment, determining the target water allocation parameter based on the initial water allocation parameter includes:

[0117] For each second grid in the groundwater simulation model, determine each first grid in the surface water simulation model corresponding to the second grid; determine the first overlapping area of ​​the second grid and the corresponding first grids, and perform weighted summation on the product value of the first overlapping area and the initial water allocation parameter in the first grid to obtain a first weighted summation value; and determine the ratio of the first weighted summation value to the area of ​​the second grid as the target water allocation value of the second grid.

[0118] Among them, the target area includes multiple first grids obtained by division based on the first grid division scale, and multiple second grids obtained by division based on the second grid division scale; that is, the terminal device can divide the target area by the first grid division scale corresponding to the surface water simulation model to obtain multiple first grids; the terminal can also divide the target area by the second grid division scale corresponding to the groundwater simulation model to obtain multiple second grids. The target area can be the entire divided grid, or it can be a specified area in the grid divided according to actual needs. The first grid division scale can be established based on the surface water simulation model, and the surface water simulation model can use the WGS84 geographic coordinate system (World Geodetic System 1984); the second grid division scale can be established based on the groundwater simulation model, and the groundwater simulation model can use the Universal Transverse Mercator Projection coordinate system.

[0119] Specifically, as shown in Figure 2, Figure 2 is a schematic diagram of the coupling of the surface water simulation model and the groundwater simulation model in an embodiment. The trapezoidal grid of the thick solid line in Figure 2 is the first grid, and the square grid of the light solid line is the second grid. Since the grid division scales of the surface water simulation model and the groundwater simulation model are different, the coordinate system conversion is involved when the data is transmitted between the surface water simulation model and the groundwater simulation model. When transferring data from the surface water simulation model to the groundwater simulation model, the specific initial process is: for each second grid in the groundwater simulation model, determine each first grid in the surface water simulation model corresponding to the second grid; determine the first overlapping area of ​​the second grid and the corresponding first grids, and perform weighted summation on the product value of the first overlapping area and the initial water allocation parameter in the first grid to obtain a first weighted sum value; determine the ratio of the first weighted sum value to the area of ​​the second grid as the target water allocation value of the second grid, and the specific expression formula of the target water allocation value of the second grid can be:

[0120]

[0121] Among them, MODFLOW p0,q0 represents the target water allocation value of the second grid in row p0 and column q0; CWatM i,j represents the initial water distribution value of the first grid in the i-th row and j-th column, Area i,j,p0,q0 Represents the first overlapping area of ​​the small grids after cutting that are contained in both the first grid in row i and column j and the second grid in row p0 and column q0. MODFLOW_Area p0,q0 Represents the area of ​​the second grid at row p0 and column q0.

[0122] In an exemplary embodiment, the specific implementation process of the step of “determining the target groundwater discharge volume based on the initial groundwater discharge volume” may include:

[0123] For each first grid in the surface water simulation model, determine each second grid in the groundwater simulation model corresponding to the first grid; determine the second overlapping area of ​​the first grid and the corresponding second grids, and perform weighted summation on the product value of the second overlapping area and the initial groundwater discharge in the second grid to obtain a second weighted summation value; and determine the ratio of the second weighted summation value to the area of ​​the first network as the target groundwater discharge of the first grid.

[0124] Specifically, for each first grid in the surface water simulation model, the specific expression for the target groundwater discharge volume of the first grid may be:

[0125]

[0126] Among them, CWatM i0,j0 represents the target groundwater discharge of the first grid in row i0 and column j0, CWatM_Area i0,j0 represents the area of ​​the first grid in row i0 and column j0, MODFLOW p,q represents the initial groundwater discharge of the MODFLOW grid in row p and column q, Area i0,j0p,q Represents the second overlapping area of ​​the small grids after cutting that are both included in the first grid at row i0 and column j0 and the second grid at row p and column q. To realize the above transfer process, the row and column numbers of each small grid corresponding to the first grid and the second grid and the area of ​​the small grid are made as the input of the coupling model, as shown in Figure 2.

[0127] In an exemplary embodiment, in a surface water simulation model, based on a preset fluid dynamics algorithm and a target groundwater discharge, a flow parameter is determined, including:

[0128] In the surface water simulation model, the ratio of the surface water flow rate and the coordinate change value of the surface water is determined as the flow change rate of the surface water, and the ratio of the cross-sectional area of ​​the surface water to the unit time length is determined as the cross-sectional change rate of the surface water; the confluence parameters are determined based on the flow change rate, the cross-sectional change rate, the capillary rise and the base flow.

[0129] Among them, the capillary rise will affect the moisture content of the underlying soil, and the base flow as part of the lateral flow q directly affects the confluence calculation of the surface water simulation model. The preset fluid dynamics algorithm can be implemented through the continuity equation and the momentum equation. The flow rate of surface water can be the flow rate of the river on the surface.

[0130] Specifically, in the surface water simulation model, the terminal determines the ratio of the flow rate of surface water and the coordinate change value of surface water as the flow rate change rate of surface water, and determines the ratio of the cross-sectional area of ​​surface water to the unit time as the cross-sectional change rate of surface water; the terminal can establish a continuity equation based on the flow change rate, the cross-sectional change rate and the lateral flow, and determine the cross-sectional area of ​​the surface water flow based on the momentum equation. The terminal can determine the confluence of surface water based on the flow change rate, the cross-sectional change rate, the capillary rise and the base flow, and the confluence is affected by groundwater. The specific expression of the continuity equation can be:

[0131]

[0132] Among them, Q is the river flow, A is the cross-sectional area of ​​the water flow, x is the spatial coordinate along the direction of the water flow, t is the time, and q is the lateral flow per unit length.

[0133] The specific expression of the cross-sectional area of ​​the water flow can be:

[0134]

[0135] Among them, α and β are the parameters of the momentum equation, β=0.6.

[0136] The parameter α of the momentum equation can be solved according to the Manning formula. The specific expression of the parameter α can be:

[0137]

[0138] Among them, n is the Manning coefficient, P is the wetted perimeter of the water-passing section, and S0 is the terrain gradient.

[0139] In an exemplary embodiment, the surface water and groundwater simulation coupling method further includes:

[0140] Based on environmental data and initial model parameters, an initial surface water simulation model and an initial groundwater simulation model are determined; in the initial surface water simulation model and the initial groundwater simulation model corresponding to the first spatial resolution, the initial model parameters are calibrated based on a preset calibration algorithm to obtain first model parameters; in the initial surface water simulation model and the initial groundwater simulation model corresponding to the second spatial resolution, the first model parameters are calibrated based on a preset calibration algorithm to obtain second model parameters; based on the second model parameters, the initial surface water simulation model and the initial groundwater simulation model are optimized to obtain a surface water simulation model and a groundwater simulation model.

[0141] The preset calibration algorithm may be a non-dominated sorting genetic algorithm II (NSGA-II) in a distributed evolutionary algorithm in Python (DEAP) library. The first spatial resolution is different from the second spatial resolution, and optionally, the first spatial resolution is lower than the second spatial resolution. The model parameters may affect the calculation results of the surface water simulation model and the groundwater simulation model.

[0142] Specifically, the terminal can process the environmental data to obtain multiple input files, and construct an initial surface water simulation model and an initial groundwater simulation model through the input files and output model parameters. Exemplarily, the input files can be shown in Tables 1 and 2, where Table 1 is the meteorological input files and their sources in the surface water simulation model and the groundwater simulation model, and Table 2 is the main input files in the surface water simulation model and the groundwater simulation model except for the meteorological input.

[0143] Table 1

[0144]

[0145] Table 2

[0146]

[0147] The terminal can process the meteorological data in the environmental data to obtain the specific process of the input file, including the meteorological data including precipitation, daily average temperature, daily maximum temperature, daily minimum temperature, surface air pressure, wind speed, downward shortwave radiation and downward longwave radiation. The terminal can use the Penman formula to calculate evapotranspiration, and obtain meteorological data of different spatial resolutions by inverse distance interpolation and Gaussian filtering in turn, and determine the variation parameter σ in the Gaussian filter according to the spatial resolution. Exemplarily, the meteorological data of different spatial resolutions can be 1 arc min×1 arc min and 5 arc min×5 arc min, respectively. The variation parameter σ of the Gaussian filter template is inversely proportional to the spatial resolution, that is, when 1arc min×1arc min, σ=10, and when 5arc min×5 arc min, σ=2. Among them, the inverse distance interpolation can make the range of known data points a set of known data points with the closest distance in all directions around the grid to be interpolated, and use the distance square inverse distance interpolation. Two-dimensional Gaussian filtering can smooth the meteorological input, and the Gaussian filter template obeys the two-dimensional Gaussian distribution. Due to the uncertainty of some precipitation, the ratio of the measured monthly total precipitation and the monthly total precipitation of the input data is calculated as the multiplication correction factor. Based on the multiplication factor, the daily precipitation of the input data within the month is multiplied by the correction factor to ensure the consistency of the meteorological input. Finally, the input file corresponding to the meteorological data is obtained. Other input files besides meteorological input can be produced by the terminal through publicly released data sets or water resources bulletins. The specific input files of the model are not limited here.

[0148] In order to give full play to the advantages of multi-core parallel computing of server CPU and single-core performance of desktop CPU, the terminal can preliminarily calibrate the initial model parameters in the initial surface water simulation model and the initial groundwater simulation model corresponding to the first resolution based on a preset calibration algorithm to obtain first model parameters, and then calibrate the first model parameters in the initial surface water simulation model and the initial groundwater simulation model corresponding to the second spatial resolution based on the preset calibration algorithm to obtain second model parameters. Exemplarily, the first resolution can correspond to an initial surface water simulation model with a spatial resolution of 5 arc min coupled with a spatial resolution of 5 km to obtain an initial groundwater simulation model, and the second resolution can correspond to an initial surface water simulation model with a spatial resolution of 1 arc min coupled with a spatial resolution of 1 km to obtain an initial groundwater simulation model.

[0149] In one embodiment, the specific process of calibrating the initial model parameters using a preset calibration algorithm is as follows: first, a population of initial size (model parameter set) is generated according to the set model parameters and parameter ranges. Secondly, the model is run according to each set of parameters and the objective function is calculated. Then, according to the objective function, excellent individuals in the population are selected for crossover mutation to obtain the next generation population. In the process of generating the next generation population, an elite strategy is used and crowding is considered to improve the level and diversity of the population. Finally, the process of crossover mutation and calculation of the objective function is repeated to obtain the Pareto optimal solution, which is the model parameter after calibration.

[0150] The terminal can optimize the relevant parameters in the initial surface water simulation model and the initial groundwater simulation model based on the second model parameters to obtain the surface water simulation model and the groundwater simulation model. The model parameters involved in the calibration of the coupling model are shown in Table 3.

[0151] Table 3

[0152]

[0153] In this embodiment, the measured burial depth changes of submerged wells in the target area are weighted averaged by the water supply distribution of the target area to obtain the monthly time series of submerged reserves in the target area for calibration. In order to better simulate the change trend of groundwater reserves in the target area, this embodiment preferentially selects the parameter group with the smallest difference from the measured change trend in the Pareto front, while taking into account the correlation coefficient and morphological similarity with the measured reserves.

[0154] In an exemplary embodiment, the surface water and groundwater simulation coupling method specifically includes the following steps:

[0155] First, the CWatM model and the MODFLOW model are initialized according to the input files and related model parameters, and two independent models of CWatM and MODFLOW are built to prepare for daily iterative simulation. Secondly, CWatM calculates the daily groundwater demand and infiltration, and transfers the demand and infiltration to MODFLOW. Then, MODFLOW calculates the groundwater head, capillary rise and base flow and transfers them to CWatM. Finally, CWatM calculates the confluence according to the information returned by MODFLOW to complete an iteration. In addition, the non-dominated sorting genetic algorithm II (NSGA-II) in the DEAP library is used to calibrate the model parameters according to the changes in the measured burial depth of submerged wells within the target area, and the optimized CWatM model and MODFLOW model are obtained. Based on the optimized CWatM model and MODFLOW model, a specific application scenario is set for simulation.

[0156] In this example, based on the CWatM model and the MODFLOW model, a coupling model of CWatM and MODFLOW daily exchange of vertical fluxes was built, and a "low-high spatial resolution" two-step calibration method was proposed to achieve large-scale (Haihe River Basin, 320,000 km) calibration under limited computing resources. 2 )High spatial resolution (1 arc min / 1 km) simulation can achieve effective simulation of the interaction process between surface water and groundwater in areas with strong human activities and improve the coupling accuracy of surface water and groundwater.

[0157] In one embodiment, as shown in FIG3 , FIG3 is a comparison diagram of the simulation and measured results of the groundwater reserve change anomaly in a certain area in XX-XX in one embodiment. Among them, the solid line represents the measured results, the dotted line represents the simulation results, the horizontal axis is the year, and the vertical axis is the groundwater reserve anomaly. In the rate period of XX-XX, the linear correlation coefficient between the groundwater reserves in the target area simulated by the model and the groundwater reserves in the target area based on the actual measurement can reach 0.90, and the difference in the groundwater reserve change rate obtained by linearly fitting the groundwater reserve change sequence is about 0.4 mm / yr. The model simulation results are good, and it has the potential to be promoted in other areas with strong human activities in northern China, and has high practical value.

[0158] It should be understood that, although the various steps in the flowcharts involved in the above-mentioned embodiments are displayed in sequence according to the indication of the arrows, these steps are not necessarily executed in sequence according to the order indicated by the arrows. Unless there is a clear explanation in this article, the execution of these steps does not have a strict order restriction, and these steps can be executed in other orders. Moreover, at least a part of the steps in the flowcharts involved in the above-mentioned embodiments can include multiple steps or multiple stages, and these steps or stages are not necessarily executed at the same time, but can be executed at different times, and the execution order of these steps or stages is not necessarily to be carried out in sequence, but can be executed in turn or alternately with other steps or at least a part of the steps or stages in other steps.

[0159] Based on the same inventive concept, the embodiment of the present application also provides a surface water and groundwater simulation coupling device for implementing the surface water and groundwater simulation coupling method involved above. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme recorded in the above method, so the specific limitations in one or more surface water and groundwater simulation coupling device embodiments provided below can refer to the above limitations on the surface water and groundwater simulation coupling method, which will not be repeated here.

[0160] In an exemplary embodiment, as shown in FIG4 , a surface water and groundwater simulation coupling device 40 is provided, comprising: a first determination module 41, a second determination module 42 and a third determination module 43, wherein:

[0161] A first determination module 41 is used to determine target water allocation parameters based on environmental data of a target area in a surface water simulation model;

[0162] A second determination module 42 is used to determine a groundwater flow function in a groundwater simulation model based on a finite difference simulation algorithm and a target water volume allocation parameter, and to determine a groundwater head based on the groundwater flow function; and to determine a target groundwater discharge based on the groundwater head;

[0163] The third determination module 43 is used to determine the confluence parameters corresponding to the target area in the surface water simulation model based on a preset fluid dynamics algorithm and a target groundwater discharge, where the confluence parameters represent the degree of mutual influence between surface water and groundwater.

[0164] In one of the embodiments, the environmental data includes at least: total water demand, soil unsaturated hydraulic conductivity, soil moisture content and soil residual moisture content, and the first determination module 41 is specifically used to divide the total water demand according to the population density distribution strategy in the surface water simulation model to obtain the water demand;

[0165] The product of soil unsaturated hydraulic conductivity and unit time is determined as the water transfer amount in the soil; the difference between soil water content and soil residual water content is determined as the effective use water content of the soil;

[0166] Determine the minimum value of water transport and effective soil water content, and determine the deep seepage as the product of the minimum value and the area ratio of the unsaturated degree of the underground aquifer;

[0167] Based on the product of the amount of water that reaches the soil surface after being intercepted by the vegetation canopy and the relative saturation of the soil, the preferential flow is determined, and the deep infiltration and preferential flow are determined as the infiltration amount;

[0168] The water demand and the infiltration amount are determined as initial water allocation parameters, and the target water allocation parameters are determined based on the initial water allocation parameters.

[0169] In one embodiment, the second determination module 42 is specifically configured to determine the difference between the groundwater head and the spatial range of the lower boundary of the soil if the groundwater head is higher than the spatial range of the lower boundary of the soil corresponding to the target area, and determine the product value of the difference and the hydraulic conductivity of the groundwater aquifer as the total discharge of groundwater;

[0170] The product of the total amount of discharged water and the proportion of the water area corresponding to the target area is determined as the base flow generated by the recharge from groundwater to surface water, and the difference between the total amount of discharged water and the base flow is determined as the capillary rise;

[0171] The base flow and capillary rise are determined as the initial groundwater discharge, and the target groundwater discharge is determined based on the initial groundwater discharge.

[0172] In one embodiment, the target area includes a plurality of first grids obtained by division based on a first grid division scale, and a plurality of second grids obtained by division based on a second grid division scale; a first determination module 41 is specifically used to determine, for each second grid in a groundwater simulation model, each first grid in a surface water simulation model corresponding to the second grid; determine a first overlapping area between the second grid and each corresponding first grid, perform weighted summation on the product value of the first overlapping area and the initial water allocation parameter in the first grid to obtain a first weighted summation value; and determine the ratio of the first weighted summation value to the area of ​​the second grid as the target water allocation value of the second grid.

[0173] In one embodiment, the second determination module 42 is specifically used to determine, for each first grid in the surface water simulation model, each second grid in the groundwater simulation model corresponding to the first grid; determine the second overlapping area of ​​the first grid and the corresponding second grids, and weighted sum the product value of the second overlapping area and the initial groundwater discharge in the second grid to obtain a second weighted sum value; and determine the ratio of the second weighted sum value to the area of ​​the first network as the target groundwater discharge of the first grid.

[0174] In one of the embodiments, in the surface water simulation model, the third determination module 43 is used to determine the ratio of the flow rate of surface water to the coordinate change value of the surface water as the flow rate change rate of surface water, and determine the ratio of the cross-sectional area of ​​surface water to the unit time length as the cross-sectional change rate of surface water in the surface water simulation model;

[0175] Determine the sink parameters based on the rate of change of flow, the rate of change of cross-section, as well as the capillary rise and base flow.

[0176] In one embodiment, the third determination module 43 is further used to determine an initial surface water simulation model and an initial groundwater simulation model based on the environmental data and the initial model parameters;

[0177] In an initial surface water simulation model and an initial groundwater simulation model corresponding to a first spatial resolution, initial model parameters are calibrated based on a preset calibration algorithm to obtain first model parameters;

[0178] In the initial surface water simulation model and the initial groundwater simulation model corresponding to the second spatial resolution, calibrate the first model parameters based on a preset calibration algorithm to obtain second model parameters;

[0179] The initial surface water simulation model and the initial groundwater simulation model are optimized based on the second model parameters to obtain a surface water simulation model and a groundwater simulation model.

[0180] Each module in the surface water and groundwater simulation coupling device can be implemented in whole or in part by software, hardware, or a combination thereof. Each module can be embedded in or independent of a processor in a computer device in the form of hardware, or can be stored in a memory in a computer device in the form of software, so that the processor can call and execute operations corresponding to each module.

[0181] In an exemplary embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be shown in FIG5. The computer device includes a processor, a memory, an input / output interface, a communication interface, a display unit, and an input device. The processor, the memory, and the input / output interface are connected via a system bus, and the communication interface, the display unit, and the input device are connected to the system bus via the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal in a wired or wireless manner, and the wireless manner may be implemented through WIFI, a mobile cellular network, near field communication (NFC) or other technologies. When the computer program is executed by the processor, a surface water and groundwater simulation coupling device method is implemented. The display unit of the computer device is used to form a visually visible picture, which may be a display screen, a projection device, or a virtual reality imaging device. The display screen can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer device can be a touch layer covering the display screen, or a button, trackball or touchpad set on the computer device shell, or an external keyboard, touchpad or mouse.

[0182] Those skilled in the art will understand that the structure shown in FIG. 5 is merely a block diagram of a partial structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have a different arrangement of components.

[0183] In one embodiment, a computer device is provided, including a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps in the above-mentioned method embodiments when executing the computer program.

[0184] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the steps in the above-mentioned method embodiments are implemented.

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

[0186] A person of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment method can be completed by instructing the relevant hardware through a computer program, and 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 above-mentioned methods. Among them, any reference to the memory, database or other medium used in the embodiments provided in the present application can include at least one of non-volatile memory and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. As an illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM). The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. Non-relational databases may include distributed databases based on blockchains, etc., but are not limited to this. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, an artificial intelligence (AI) processor, etc., but are not limited to this.

[0187] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, 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.

[0188] The above-described embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the present application. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the present application shall be subject to the attached claims.

Claims

1. A surface water and groundwater simulation coupling method, characterized in that: The method comprises: In the surface water simulation model, the target water allocation parameters are determined based on the environmental data of the target area; In a groundwater simulation model, a groundwater flow function is determined based on a finite difference simulation algorithm and the target water volume allocation parameter, and a groundwater head is determined based on the groundwater flow function; and a target groundwater discharge is determined based on the groundwater head; In the surface water simulation model, based on a preset fluid dynamics algorithm and the target groundwater discharge, a confluence parameter corresponding to the target area is determined, wherein the confluence parameter represents the degree of mutual influence between the surface water and the groundwater; Wherein, the method further comprises: determining an initial surface water simulation model and an initial groundwater simulation model based on the environmental data and initial model parameters; In the initial surface water simulation model and the initial groundwater simulation model corresponding to the first spatial resolution, calibrating the initial model parameters based on a preset calibration algorithm to obtain first model parameters; In the initial surface water simulation model and the initial groundwater simulation model corresponding to the second spatial resolution, calibrating the first model parameters based on the preset calibration algorithm to obtain second model parameters; The initial surface water simulation model and the initial groundwater simulation model are optimized based on the second model parameters to obtain the surface water simulation model and the groundwater simulation model.

2. The method according to claim 1, characterized in that The environmental data at least includes: total water demand, soil unsaturated hydraulic conductivity, soil moisture content and soil residual moisture content. In the surface water simulation model, the target water allocation parameters are determined based on the environmental data, including: In the surface water simulation model, the total water demand is divided according to the population density distribution strategy to obtain the water demand; The product of the unsaturated hydraulic conductivity of the soil and the unit time length is determined as the water transfer amount in the soil; the difference between the soil water content and the residual water content of the soil is determined as the effective use water content of the soil; Determine the minimum value of the water transfer amount and the effective utilization water content of the soil, and determine the product of the minimum value and the area ratio of the unsaturation of the underground aquifer as the deep leakage amount; Determine the preferential flow based on the product value of the amount of water reaching the soil surface after being intercepted by the vegetation canopy and the relative saturation of the soil, and determine the deep infiltration amount and the preferential flow as the infiltration amount; The water demand and the infiltration amount are determined as initial water allocation parameters, and the target water allocation parameters are determined based on the initial water allocation parameters.

3. The method according to claim 2, characterized in that The determining of a target groundwater discharge based on the groundwater head comprises: If the groundwater head is higher than the spatial range of the lower boundary of the soil corresponding to the target area, the difference between the groundwater head and the spatial range of the lower boundary of the soil is determined, and the product of the difference and the hydraulic conductivity of the groundwater aquifer is determined as the total discharge volume of groundwater; The product of the total amount of discharged water and the proportion of the water area corresponding to the target area is determined as the base flow generated by the recharge from groundwater to surface water, and the difference between the total amount of discharged water and the base flow is determined as the capillary rise; The base flow and the capillary rise are determined as an initial groundwater discharge, and a target groundwater discharge is determined based on the initial groundwater discharge.

4. The method according to claim 3, characterized in that The target area includes a plurality of first grids obtained by division based on a first grid division scale, and a plurality of second grids obtained by division based on a second grid division scale; The determining of the target water allocation parameter based on the initial water allocation parameter comprises: For each of the second grids in the groundwater simulation model, determining each of the first grids in the surface water simulation model corresponding to the second grid; Determine a first overlapping area between the second grid and the corresponding first grids, perform weighted summation on the product value of the first overlapping area and the initial water allocation parameter in the first grid to obtain a first weighted summation value; and determine a ratio of the first weighted summation value to the area of ​​the second grid as a target water allocation value for the second grid.

5. The method according to claim 4, characterized in that The determining of a target groundwater discharge based on the initial groundwater discharge comprises: For each of the first grids in the surface water simulation model, determine each of the second grids in the groundwater simulation model corresponding to the first grid; determine a second overlapping area between the first grid and each of the corresponding second grids, and perform weighted summation on the product of the second overlapping area and the initial groundwater discharge in the second grid to obtain a second weighted summation value; and determine the ratio of the second weighted summation value to the area of ​​the first grid as the target groundwater discharge of the first grid.

6. The method according to claim 3, characterized in that In the surface water simulation model, based on a preset fluid dynamics algorithm and the target groundwater discharge, determining the confluence parameters includes: In the surface water simulation model, the ratio of the flow rate of surface water to the coordinate change value of surface water is determined as the flow rate change rate of surface water, and the ratio of the cross-sectional area of ​​surface water to the unit time length is determined as the cross-sectional change rate of surface water; A sink flow parameter is determined based on the flow rate change rate, the cross-sectional change rate, the capillary rise, and the base flow.

7. A surface water and groundwater simulation coupling device, characterized in that: The device comprises: A first determination module is used to determine target water allocation parameters based on environmental data of a target area in a surface water simulation model; A second determination module is used to determine a groundwater flow function in a groundwater simulation model based on a finite difference simulation algorithm and the target water volume allocation parameter, and to determine a groundwater head based on the groundwater flow function; and to determine a target groundwater discharge based on the groundwater head; A third determination module is used to determine, in the surface water simulation model, a confluence parameter corresponding to the target area based on a preset fluid dynamics algorithm and the target groundwater discharge, wherein the confluence parameter represents the degree of mutual influence between the surface water and the groundwater; Wherein, the third determination module is further used to determine an initial surface water simulation model and an initial groundwater simulation model based on the environmental data and the initial model parameters; In the initial surface water simulation model and the initial groundwater simulation model corresponding to the first spatial resolution, calibrating the initial model parameters based on a preset calibration algorithm to obtain first model parameters; In the initial surface water simulation model and the initial groundwater simulation model corresponding to the second spatial resolution, calibrating the first model parameters based on the preset calibration algorithm to obtain second model parameters; The initial surface water simulation model and the initial groundwater simulation model are optimized based on the second model parameters to obtain the surface water simulation model and the groundwater simulation model.

8. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 6 are implemented.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

Images