Water resource allocation method based on influence of riverside exploitation underground water

Abstract

The invention discloses a water resource allocation method based on influence of riverside exploitation of underground water, belongs to the field of water resource allocation, and can solve the problem of excessive exploitation of underground water. The method comprises the steps that S1, minimization optimization is carried out on the first total water deficit of a drainage basin in a deployment period before underground water exploitation, and the reservoir outlet water amount of each reservoir in each time period in the deployment period is obtained; s2, the total water supply guarantee rate of all the water supply objects in the allocation period is optimized to the maximum degree according to the reservoir outlet water amount and the river channel loss water amount of the drainage basin, and the water supply amount of the drainage basin for each water supply object in each time period in the allocation period is obtained; s3, performing minimization optimization on the second total water deficit of the drainage basin in the allocation period after the underground water exploitation according to the reservoir-out water volume and the water supply volume to obtain the deprived water volume of the underground water exploitation on the drainage basin, and correcting the water volume loss of the river channel according to the deprived water volume; and S4, repeating the steps S2 and S3 to obtain a surface water-underground water combined allocation scheme. The method is used for water resource allocation.

Description

Technical Field

[0001] This invention relates to a water resource allocation method based on the impact of groundwater extraction along rivers, belonging to the field of water resource allocation technology. Background Technology

[0002] Surface water resources in the basin provide crucial water security for agricultural production and oasis economy. However, with rapid population growth, expanding arable land, and increasing artificial oases, the supply of surface water is gradually becoming insufficient to meet the water demands of these activities, necessitating groundwater extraction to compensate for the shortfall. However, indiscriminate and excessive groundwater extraction can undermine the sustainability of groundwater resources and lead to the deterioration of the basin's ecological environment. Therefore, it is essential to develop a scientifically sound and rational surface water-groundwater co-allocation plan to meet the basin's water needs while maintaining the sustainability of groundwater resources. Summary of the Invention

[0003] This invention provides a water resource allocation method based on the impact of groundwater extraction along rivers, which can solve the problem of watershed ecological environment deterioration caused by blind and excessive groundwater extraction.

[0004] This invention provides a water resource allocation method based on the impact of groundwater extraction along rivers, the method comprising:

[0005] S1. Based on the design parameters of all reservoirs in the basin, minimize the first total water shortage of the basin during the allocation period before groundwater extraction, and obtain the outflow of water from each reservoir in each time period during the allocation period.

[0006] S2. Based on the outflow of water and the water loss in the river channels of the basin, maximize the overall water supply guarantee rate for all water supply targets during the allocation period to obtain the water supply volume of the basin to each water supply target in each time period during the allocation period.

[0007] S3. Based on the outflow and supply of water, minimize the second total water shortage in the basin during the allocation period after groundwater extraction to obtain the water loss from groundwater extraction in the basin, and correct the water loss in the river channel based on the water loss.

[0008] S4. Repeat S2 and S3 to determine a joint surface water-groundwater allocation scheme based on the obtained multiple sets of water supply and water capture amounts, and allocate water resources based on the joint surface water-groundwater allocation scheme.

[0009] Optionally, in S1, based on the design parameters of all reservoirs within the basin, the first total water shortage in the basin during the allocation period before groundwater extraction is minimized and optimized, specifically including:

[0010] A joint scheduling model for the reservoir group is constructed based on the design parameters of all reservoirs in the basin.

[0011] Based on the joint scheduling model of the reservoir group, the first total water shortage in the basin during the allocation period before groundwater extraction is minimized and optimized.

[0012] Optionally, a joint scheduling model for the reservoir group is constructed based on the design parameters of all reservoirs within the basin, specifically including:

[0013] A first objective function is constructed based on the ecological base flow of the watershed and the water demand of all water supply objects. The first objective function is used to minimize the first total water shortage of the watershed during the allocation period before groundwater extraction.

[0014] The first constraint condition is constructed based on the design parameters of all reservoirs in the basin. The first constraint condition is used to constrain the operation process of each reservoir.

[0015] The joint scheduling model of the reservoir group is constructed based on the first objective function and the first constraint.

[0016] Optionally, in S2, based on the outflow of water from the reservoir and the water loss in the river channels of the basin, the overall water supply guarantee rate for all water supply targets during the allocation period is maximized and optimized, specifically including:

[0017] A surface water resource allocation model for the watershed is constructed based on the water demand of all water-supplying entities.

[0018] Based on the outflow of water from the reservoir and the water loss in the river channels of the basin, the total water supply guarantee rate for all water supply targets during the allocation period is maximized using the surface water resource allocation model.

[0019] Optionally, a surface water resource allocation model can be constructed based on the water demand of all water-supplying entities, specifically including:

[0020] A second objective function is constructed based on the water demand of all water supply targets. The second objective function is used to maximize the total water supply guarantee rate of all water supply targets during the allocation period.

[0021] A second constraint is constructed based on the water demand of all water-supplying entities. This second constraint is used to constrain the water resource allocation process.

[0022] The surface water resource allocation model is constructed based on the second objective function and the second constraint.

[0023] Optionally, in S3, based on the outflow and supply of water, the second total water shortage in the watershed during the allocation period after groundwater extraction is minimized and optimized, specifically including:

[0024] Construct a groundwater resource allocation model for the aforementioned watershed;

[0025] Based on the outflow and supply of water, the second total water shortage in the basin during the allocation period after groundwater extraction is minimized using the groundwater resource allocation model.

[0026] Optionally, a groundwater resource allocation model for the watershed is constructed, specifically including:

[0027] The groundwater level in different river sections of the basin after groundwater extraction is determined based on the outflow and supply volume.

[0028] The amount of water extracted from different river sections after groundwater extraction is determined based on the groundwater level.

[0029] The groundwater resource allocation model is constructed based on the groundwater level and the amount of water seized.

[0030] Optionally, the groundwater level in different river sections of the basin after groundwater extraction is determined based on the outflow and supply volumes, specifically including:

[0031] The water level of different river sections in the basin is determined based on the outflow and supply of water.

[0032] The groundwater level in different river sections after groundwater extraction is determined based on the river water level.

[0033] Optionally, the amount of water captured from different river sections after groundwater extraction is determined based on the groundwater level, specifically including:

[0034] The exchange flow rate between the river channel and the underground aquifer after groundwater extraction is determined based on the groundwater level.

[0035] The amount of water extracted from different river sections after groundwater extraction is determined based on the exchange flow rate.

[0036] The present invention also provides a water resource allocation system based on the impact of groundwater extraction along rivers, the system comprising:

[0037] The outflow water calculation module is used to minimize the first total water shortage in the basin before groundwater extraction during the allocation period based on the design parameters of all reservoirs in the basin, and to obtain the outflow water of each reservoir in each time period during the allocation period.

[0038] The water supply calculation module is used to maximize the overall water supply guarantee rate of all water supply objects during the allocation period based on the outflow of water and the water loss in the river channels of the basin, so as to obtain the water supply volume of the basin to each water supply object in each time period during the allocation period.

[0039] The water loss calculation module is used to minimize the second total water shortage in the basin during the allocation period after groundwater extraction based on the outflow and supply of water, to obtain the water loss of the basin caused by groundwater extraction, and to correct the river channel water loss based on the water loss.

[0040] The joint allocation module is used to iteratively calculate the water supply and water capture volume, determine the surface water-groundwater joint allocation scheme based on the obtained multiple sets of water supply and water capture volume, and allocate water resources based on the surface water-groundwater joint allocation scheme.

[0041] The present invention also provides a storage device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that the processor, when executing the computer program, implements the steps of any of the above-described water resource allocation methods based on the impact of groundwater extraction along a river.

[0042] The beneficial effects that this invention can produce include:

[0043] This invention optimizes the water supply by minimizing the total water shortage in the watershed during the allocation period before groundwater extraction, thus determining the outflow from each reservoir in each time period. It then optimizes the water supply guarantee rate for all water-supplying entities during the allocation period, yielding the water supply volume to each entity in each time period. Based on the outflow, supply, and river loss, the invention determines the water loss from groundwater extraction. The river loss is then adjusted based on the water loss, and the supply and loss volumes are continuously optimized through iterative calculations to ultimately determine a joint surface water-groundwater allocation scheme. This approach ensures the water needs of the target water-supplying entities while avoiding indiscriminate and excessive groundwater extraction, achieving a rational allocation of surface water and groundwater. This is beneficial for maintaining the sustainability of groundwater resources and protecting the watershed's ecological environment. Attached Figure Description

[0044] Figure 1 A flowchart of a water resource allocation method based on the impact of groundwater extraction along a river, provided as an embodiment of the present invention. Detailed Implementation

[0045] The present invention will now be described in detail with reference to the embodiments, but the present invention is not limited to these embodiments.

[0046] This invention provides a water resource allocation method based on the impact of groundwater extraction along rivers, such as... Figure 1 As shown, the method includes:

[0047] S1. Based on the design parameters of all reservoirs within the basin, minimize the first total water shortage during the allocation period before groundwater extraction, and obtain the outflow volume of each reservoir for each time period during the allocation period. Specifically, this includes:

[0048] A joint scheduling model for the reservoir group is constructed based on the design parameters of all reservoirs in the basin. Then, based on the joint scheduling model, the first total water shortage in the basin during the allocation period before groundwater extraction is minimized to obtain the outflow of water from each reservoir during each period of the allocation period when the first total water shortage is minimized.

[0049] In this embodiment, the joint scheduling model of the reservoir group consists of two parts: a first objective function and a first constraint condition. The construction process of each part is as follows:

[0050] 1. Construct the first objective function.

[0051] A first objective function is constructed based on the ecological base flow of the watershed and the water demand of all water-supplying entities. This first objective function is used to minimize the initial total water shortage of the watershed during the allocation period before groundwater extraction. The water-supplying entities include agricultural irrigation areas and oasis landscapes.

[0052] The first objective function can be expressed as:

[0053] (1)

[0054] In formula (1), This represents the first total water shortage (in 100 million m³) in the watershed during the allocation period before groundwater extraction. 3 ); Indicates the period of allocation Water shortage for all water-supply recipients during the period (100 million m³) 3 ); Indicates the period of allocation Water shortage during the period of ecological flow release (100 million m³) 3 ); Indicates the period of allocation Natural groundwater recharge into the river during a given period (100 million m³) 3 ); This indicates the total number of time periods (inclusive) in the allocation period.

[0055] in, Water shortage for all water-supply recipients during the period According to Once the water demand of all water-supplying entities is determined for a given period, the difference between the water demand of all water-supplying entities and the water supply from the basin to all water-supplying entities is the water shortage of all water-supplying entities. Water shortage during the period of ecological flow release According to The ecological base flow of the watershed is determined for a given time period. The difference between the ecological base flow and the ecological discharge from the watershed is the water shortage due to ecological flow. The total number of allocation periods and time periods is also considered. The allocation period can be flexibly set according to the allocation needs. For example, in this embodiment, the allocation period is set to 50 years, with one month as a time period, for a total of 600 time periods.

[0056] 2. Construct the first constraint condition.

[0057] The first constraint condition is constructed based on the design parameters of all reservoirs within the basin. This first constraint condition is used to constrain the operation of each reservoir. The design parameters of the reservoirs include the minimum design capacity, the maximum design capacity, and the maximum design discharge volume.

[0058] The first constraint includes:

[0059] 1) Water balance constraints:

[0060] (2)

[0061] In formula (2), express Reservoir capacity during a given period (100 million m³) 3 ); express Reservoir capacity during a given period (100 million m³) 3 ); express Reservoir inflow during a given period (100 million m³) 3 ); express Reservoir outflow during a given period (100 million m³) 3 ); Indicate Water loss due to evaporation and seepage during a given period (100 million m³) 3 ).

[0062] in, and Water can be taken from the reservoir separately. Time period and Initial storage capacity for a given period of time.

[0063] 2) Reservoir capacity constraints:

[0064] (3)

[0065] In formula (3), express Reservoir capacity during a given period (100 million m³) 3 ); This indicates the minimum design capacity of the reservoir (in billions of cubic meters).3 ); This indicates the maximum designed storage capacity of the reservoir (in billions of cubic meters). 3 ).

[0066] 3) Reservoir discharge constraints:

[0067] (4)

[0068] In formula (4), express Reservoir discharge volume during a given period (100 million m³) 3 ); This indicates the reservoir's maximum designed discharge capacity (100 million m³). 3 ).

[0069] 4) Non-negativity constraint: All variables must be non-negative.

[0070] In this embodiment, when running the joint scheduling model of the reservoir group, the final storage capacity of each reservoir in each time period is used as the decision variable. Based on the first constraint, the first objective function is used to minimize the first total water shortage of the watershed during the allocation period before groundwater extraction, and the output is the outflow of water from each reservoir in each time period during the allocation period when the first total water shortage is minimized. Among them, the upper limit of the final storage capacity is the storage capacity corresponding to the flood limit water level during the flood season, and the storage capacity corresponding to the normal storage water level during the non-flood season. The lower limit of the decision variable is the dead water level.

[0071] S2. Based on the outflow of water from the reservoir and the water loss in the river channels of the basin, the overall water supply guarantee rate for all water supply recipients during the allocation period is maximized through optimization, resulting in the water supply volume to each water supply recipient for each time period within the basin during the allocation period. Specifically, this includes:

[0072] A surface water resource allocation model for the watershed is constructed based on the water demand of all water supply recipients. Then, based on the outflow of water and the water loss in the river channels of the watershed, the total water supply guarantee rate for all water supply recipients during the allocation period is maximized using the surface water resource allocation model. This yields the water supply volume to each water supply recipient in each time period when the total water supply guarantee rate is maximized.

[0073] In this embodiment, the surface water resource allocation model consists of two parts: a second objective function and a second constraint condition. The construction process of each part is as follows:

[0074] 1. Construct the second objective function.

[0075] A second objective function is constructed based on the water demand of all water supply targets. This second objective function is used to maximize the overall water supply guarantee rate of all water supply targets during the allocation period.

[0076] The second objective function can be expressed as:

[0077] (5)

[0078] In formula (5), This indicates the total water supply guarantee rate for all water-supplying entities during the allocation period; Indicates the period of allocation Time-based water supply targets Water demand satisfaction coefficient; Indicates the watershed during the allocation period Time period for water supply recipients Water supply volume (100 million m³) 3 ); Indicates the period of allocation Time-based water supply targets Water demand (billion m³) 3 ); Indicates the total number of time periods (in terms of the allocation period); This indicates the total number of water supply recipients.

[0079] 2. Construct the second constraint.

[0080] A second constraint is constructed based on the water demand of all water-supplying entities. This second constraint is used to constrain the water resource allocation process.

[0081] The second constraint includes:

[0082] 1) Node water balance constraints:

[0083] (6)

[0084] In formula (6), Indicates the watershed Node at Water inflow during the period (100 million m³) 3 ); Indicates the watershed Node at Water inflow during the period (100 million m³) 3 ); Indicates the watershed Node at Water supply volume (in 100 million m³) to the target water user during different time periods 3 ); Indicates the watershed Node at River channel water loss during the period (100 million m³) 3 ).

[0085] The location of each node can be determined based on the location of water conservancy projects such as reservoirs and canal heads within the basin. Node located Downstream of the node.

[0086] 2) Water supply constraints:

[0087] (7)

[0088] In formula (7), Indicates the watershed during the allocation period Time period for water supply recipients Water supply volume (100 million m³) 3 ); Indicates the period of allocation Time-based water supply targets Water demand (billion m³) 3 ).

[0089] 3) Water supply satisfaction constraints:

[0090] (8)

[0091] In formula (8), Indicates the watershed during the allocation period Time period for water supply recipients Water supply volume (100 million m³) 3 ); Indicates the period of allocation Time-based water supply targets Water demand (billion m³) 3 ); This indicates the preset water supply satisfaction level, which can be flexibly set according to the actual situation of the specific water supply target. For example, the preset water supply satisfaction level for agricultural irrigation areas can be 70%.

[0092] 4) Non-negativity constraint: All variables must be non-negative.

[0093] In this embodiment, when running the surface water resource allocation model, the sum of the water allocation proportions of the watershed at each node is used as the decision variable. Based on the second constraint, the second objective function is used to maximize the total water supply guarantee rate for all water supply recipients during the allocation period, and the output is the water supply volume of the watershed to each water supply recipient in each time period when the total water supply guarantee rate is maximized. The upper limit of the sum of the water allocation proportions of the watershed at each node is 1, and the lower limit is 0.

[0094] S3. Based on the outflow from the reservoir output in S1 and the supply output in S2, minimize the second total water shortage in the basin during the allocation period after groundwater extraction to obtain the water loss from groundwater extraction, and correct the river channel loss based on the water loss. Specifically, this includes:

[0095] A groundwater resource allocation model is constructed. Based on the outflow of water from the reservoir output by S1 and the supply of water output by S2, the model is used to minimize the second total water shortage in the watershed during the allocation period after groundwater extraction. The model outputs the amount of water extracted from the watershed when the second total water shortage is minimized, and the amount of water lost from the river channel is corrected based on the amount of water extracted.

[0096] In this embodiment, the operation process of the groundwater resource allocation model is as follows:

[0097] 1. Determine the flow rate at the control section of different river sections based on the outflow, supply, and loss of water in the river channel.

[0098] Each river section can be divided according to the location of water conservancy projects such as reservoirs and canal heads within the basin. Control sections are set on each river section, and the water flow through the control section is the control section flow.

[0099] Based on the outflow, supply, and water loss of each river section, the flow rate at the control section of each river section can be obtained in this embodiment.

[0100] To accurately calculate the flow at the control section, it is necessary to determine the water loss in each river segment. In this embodiment, the calculation formula for the water loss in each river segment is obtained by fitting the measured data of each river segment.

[0101] Taking the Yarkand River Basin as an example, the formula for calculating the water loss in the river section between the head of the main canal and the head of the Minsheng canal in this basin, based on fitting, is as follows:

[0102] (9)

[0103] The formula for calculating the water loss in the river section between the Minsheng Canal Head and the Ailiktam Canal Head in this basin is as follows:

[0104] (10)

[0105] In formulas (9) and (10), This indicates the water loss in the river section between the head of the main canal and the head of the people's livelihood canal (in billions of cubic meters). 3 ); This indicates the water loss (in billions of cubic meters) in the river section between the Minsheng Canal Head and the Ailiktam Canal Head. 3 ); and These represent the monthly and last month's water inflow for the corresponding river sections, respectively (in 100 million cubic meters). 3 ).

[0106] 2. Based on the flow rate at the control section, determine the river level for different river sections using the Manning formula.

[0107] This embodiment establishes a flow-water level conversion relationship based on the geometric parameters of each control section and the Manning formula. Using this conversion relationship, the flow rate at each control section can be converted into the river water level at that section, thereby determining the river water level in different river segments.

[0108] 3. Determine the groundwater level in different river sections after groundwater extraction based on the river water level in different sections.

[0109] This embodiment first uses FloPy to construct a groundwater simulation model of the watershed. Considering the irregularity of the watershed boundary, directly using FloPy for modeling is quite complex, affecting the efficiency and accuracy of the modeling. Therefore, this embodiment utilizes the MODFLOW flow module in the Groundwater Modeling System (GMS) to initially construct the groundwater simulation model of the watershed. Then, this embodiment imports the initially constructed groundwater simulation model into FloPy and further adjusts and simulates the groundwater simulation model in the Python environment.

[0110] The steps for importing and adjusting a groundwater simulation model are as follows:

[0111] 1) Importing GMS Model: Through the API interface, FloPy can import the groundwater simulation model initially built in GMS into the Python environment to realize the simulation operation of the model.

[0112] 2) Revise the RIV module: Since there are many water conservancy projects in the basin, the water storage process of the water conservancy projects may cause the water level of each river section to be discontinuous. Therefore, in this embodiment, the water level of each river section is input into the RIV module to adjust the groundwater simulation model.

[0113] 3) Revising the Groundwater Evaporation Module: The groundwater evaporation module in GMS calculates groundwater evaporation based on a given groundwater evaporation rate. However, in reality, the groundwater evaporation rate varies with groundwater depth and land cover type. Therefore, this embodiment uses measured data to fit the correspondence between groundwater evaporation and groundwater depth under different land cover types, and improves the groundwater evaporation module in FloPy based on this correspondence to enhance the accuracy of groundwater evaporation calculation in the groundwater simulation model.

[0114] For example, the fitted relationship between groundwater evaporation and groundwater depth under different land cover types is as follows:

[0115] (11)

[0116] (12)

[0117] (13)

[0118] (14)

[0119] (15)

[0120] In formulas (11) to (15), The groundwater evaporation rate (mm) in arbor forests; The groundwater evaporation rate (mm) is for shrubland. The groundwater evaporation rate of the grassland (mm); The value represents the groundwater evaporation rate (mm) under bare soil conditions. The depth of groundwater (m); This is the evaporation limit intensity for submerged water; The observed evaporation rate (mm) of the evaporating dish; For parameters related to the burial depth of the submersible; The depth of the buoyancy depth is measured in meters (m).

[0121] Then, in this embodiment, the river water level and groundwater extraction plan are input into the adjusted groundwater simulation model for simulation operation.

[0122] When groundwater is extracted along a river, the continuous drop in groundwater level at the extraction well will form a regional groundwater level drawdown cone centered on the extraction well, which will enhance the migration of surrounding groundwater to the extraction well and increase the hydraulic gradient between the river and the underground aquifer. This process breaks the natural replenishment and discharge balance between surface water and groundwater, and induces unidirectional seepage replenishment from the river to the underground aquifer.

[0123] The groundwater simulation model constructed in this embodiment can simulate the groundwater transport process using a three-dimensional unsteady flow equation, which can be expressed as:

[0124] (16)

[0125] In formula (16), , , They represent the underground aquifers along direction, direction, Hydraulic conductivity in direction (anisotropy) (m / s); The water head (m) of an underground aquifer is the height of the free water surface of the water body in the underground aquifer relative to a certain reference surface. Represents source and sink terms (m) 3 / d), including replenishment and excretion items; Indicates water supply degree.

[0126] In the scenario of groundwater extraction along a river, the newly added pumping well discharge term in the source-sink term W will alter the regional water flow field, leading to a drop in the regional groundwater level. When the pumping well is near the river, as groundwater migrates, the regional groundwater level drawdown cone will spread to the river channel, causing the groundwater level along the river to drop further and increasing the amount of water seepage from the river channel. This phenomenon seriously encroaches on the surface water resources of the watershed.

[0127] Because the pumping volume and the hydrogeological conditions of the pumping well location are different, the impact range of groundwater extraction along the river will also be different. Therefore, this embodiment uses a groundwater simulation model to simulate the dynamic water exchange process between the groundwater aquifer and the river under different groundwater extraction schemes, so as to optimize the water resource allocation process.

[0128] This embodiment uses the above-described simulation process to obtain the groundwater level in different river sections after groundwater extraction.

[0129] 4. Determine the exchange flow between the river channel and the underground aquifer after groundwater extraction based on the groundwater level in different river sections after groundwater extraction.

[0130] When the river water level is higher than the groundwater level, the river replenishes the groundwater aquifer; conversely, when the river water level is lower than the groundwater level, the groundwater aquifer replenishes the river. The exchange flow between the river and the groundwater aquifer can be calculated using the RIV module, with the following formula:

[0131] (17)

[0132] In formula (17), Indicates the exchange flow between the river channel and the underground aquifer (m³). 3 / d); The coefficient of conductivity of the riverbed (m) 2 / d); Indicates the groundwater level (m); Indicates the river water level (m).

[0133] Formula (17) shows the exchange flow between the river channel and the underground aquifer. This is mainly determined by the hydraulic gradient between the river channel and the underground aquifer. When the groundwater level is higher than the river level, groundwater flows into the river; when the groundwater level is lower than the river level, surface water seeps into the underground aquifer; and when the river level is lower than the bottom elevation of the riverbed, water exchange ceases.

[0134] 5. Determine the amount of water extracted from different river sections after groundwater extraction based on the exchange flow between the river channel and the underground aquifer.

[0135] The amount of river water indirectly captured by pumping wells after groundwater extraction is called the captured water volume, and its calculation process is as follows:

[0136] (18)

[0137] In formula (18), Indicates the interception of traffic (m) 3 / d); This indicates the exchange flow rate (m³) between the river channel and the underground aquifer after groundwater extraction. 3 / d); This represents the exchange flow between the river channel and the underground aquifer when the watershed is in its natural state before groundwater extraction (m³). 3 / d).

[0138] Then, the amount of water seized in each river section can be obtained based on the seized flow rate.

[0139] 6. Adjust the amount of water lost from the river channel based on the amount of water seized.

[0140] The correction process for water loss in the river channel can be expressed as:

[0141] (19)

[0142] In formula (19), This indicates the corrected river channel water loss (in billions of cubic meters). 3 ); This indicates the amount of water lost from the river channel before the correction (in billions of cubic meters). 3 ); This indicates the amount of water seized (in billions of cubic meters). 3 ), based on the traffic seized Calculated.

[0143] In this embodiment, when running the groundwater resource allocation model, the groundwater extraction volume corresponding to each river section in the basin is used as a decision variable, and a groundwater depth threshold is set according to the actual situation of the basin to constrain the groundwater extraction volume. Then, the second total water shortage in the basin during the allocation period after groundwater extraction is minimized and optimized, the water loss volume of the basin due to groundwater extraction is output, and the river channel water loss volume is corrected according to the water loss volume.

[0144] S4. Repeat S2 and S3. Based on the obtained multiple sets of water supply and water diversion, determine the surface water-groundwater joint allocation scheme, and allocate water resources based on the surface water-groundwater joint allocation scheme.

[0145] This embodiment corrects for river channel water loss based on the amount of water seized, then iteratively calculates the water supply volume using the surface water resource allocation model constructed in S2, and then iteratively calculates the seized water volume based on the water supply volume and the groundwater resource allocation model constructed in S3. This iterative calculation is performed a preset number of times to finally obtain multiple combinations of water supply volume and seized water volume. Subsequently, the optimal combination can be selected from multiple combinations according to the actual situation and allocation needs, and the surface water-groundwater allocation process corresponding to the optimal combination is taken as the surface water-groundwater joint allocation scheme.

[0146] This embodiment uses a single-objective genetic algorithm to jointly iteratively solve the surface water resource allocation model and the groundwater resource allocation model.

[0147] Another embodiment of the present invention provides a water resource allocation system based on the impact of groundwater extraction along a river, the system comprising:

[0148] The outflow water calculation module is used to minimize the first total water shortage in the basin before groundwater extraction during the allocation period based on the design parameters of all reservoirs in the basin, and to obtain the outflow water of each reservoir in each time period during the allocation period.

[0149] The water supply calculation module is used to maximize the overall water supply guarantee rate of all water supply objects during the allocation period based on the outflow of water and the water loss in the river channels of the basin, so as to obtain the water supply volume of the basin to each water supply object in each time period during the allocation period.

[0150] The river loss correction module is used to minimize the second total water shortage in the basin during the allocation period after groundwater extraction based on the outflow and supply of water, to obtain the amount of water extracted from the basin by groundwater extraction, and to correct the amount of water lost in the river channel based on the amount of water extracted.

[0151] The joint allocation module is used to iteratively calculate the water supply and water capture volume, determine the surface water-groundwater joint allocation scheme based on the obtained multiple sets of water supply and water capture volume, and allocate water resources based on the surface water-groundwater joint allocation scheme.

[0152] Specifically, the river damage correction module includes:

[0153] The cross-sectional flow calculation unit is used to determine the control section flow of different river sections based on the outflow, supply, and river loss.

[0154] The river water level calculation unit is used to determine the river water level of different river sections based on the flow rate at the control section of different river sections.

[0155] The groundwater level calculation unit is used to determine the groundwater level in different river sections after groundwater extraction, based on the river water level in different sections.

[0156] The water extraction calculation unit is used to determine the water extraction volume of different river sections after groundwater extraction, based on the groundwater level of different river sections.

[0157] The river loss correction unit is used to correct the river channel loss volume of the corresponding river section based on the amount of water captured in different river sections.

[0158] Another embodiment of the present invention provides a storage device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that the processor executes the computer program to implement the steps of any of the above-described water resource allocation methods based on the impact of groundwater extraction along a river.

[0159] This invention optimizes the water supply by minimizing the total water shortage in the watershed during the allocation period before groundwater extraction, thus determining the outflow from each reservoir in each time period. It then optimizes the water supply guarantee rate for all water-supplying entities during the allocation period, yielding the water supply volume to each entity in each time period. Based on the outflow, supply, and river loss, the invention determines the water loss from groundwater extraction. The river loss is then adjusted based on the water loss, and the supply and loss volumes are continuously optimized through iterative calculations to ultimately determine a joint surface water-groundwater allocation scheme. This approach ensures the water needs of the target water-supplying entities while avoiding indiscriminate and excessive groundwater extraction, achieving a rational allocation of surface water and groundwater. This is beneficial for maintaining the sustainability of groundwater resources and protecting the watershed's ecological environment.

[0160] The above description is merely a few embodiments of this application and is not intended to limit this application in any way. Although this application discloses preferred embodiments as described above, it is not intended to limit this application. Any changes or modifications made by those skilled in the art without departing from the scope of the technical solution of this application using the disclosed technical content are equivalent to equivalent implementation cases and fall within the scope of the technical solution.

Claims

1. A water resource allocation method based on the impact of groundwater extraction along rivers, characterized in that, The method includes: S1. Based on the design parameters of all reservoirs in the basin, minimize the first total water shortage of the basin during the allocation period before groundwater extraction, and obtain the outflow of water from each reservoir in each time period during the allocation period. S2. Based on the outflow of water and the water loss in the river channels of the basin, maximize the overall water supply guarantee rate for all water supply targets during the allocation period to obtain the water supply volume of the basin to each water supply target in each time period during the allocation period. S3. Based on the outflow and supply of water, minimize the second total water shortage in the basin during the allocation period after groundwater extraction to obtain the water loss from groundwater extraction in the basin, and correct the water loss in the river channel based on the water loss. S4. Repeat S2 and S3 to determine a joint surface water-groundwater allocation scheme based on the obtained multiple sets of water supply and water capture amounts, and allocate water resources based on the joint surface water-groundwater allocation scheme.

2. The method according to claim 1, characterized in that, S1 optimizes the minimization of the first total water shortage in the basin during the allocation period before groundwater extraction, based on the design parameters of all reservoirs within the basin. Specifically, this includes: A joint scheduling model for the reservoir group is constructed based on the design parameters of all reservoirs in the basin. Based on the joint scheduling model of the reservoir group, the first total water shortage in the basin during the allocation period before groundwater extraction is minimized and optimized.

3. The method according to claim 2, characterized in that, A joint scheduling model for the reservoir group is constructed based on the design parameters of all reservoirs within the basin, specifically including: A first objective function is constructed based on the ecological base flow of the watershed and the water demand of all water supply objects. The first objective function is used to minimize the first total water shortage of the watershed during the allocation period before groundwater extraction. The first constraint condition is constructed based on the design parameters of all reservoirs in the basin. The first constraint condition is used to constrain the operation process of each reservoir. The joint scheduling model of the reservoir group is constructed based on the first objective function and the first constraint.

4. The method according to claim 1, characterized in that, S2 optimizes the overall water supply guarantee rate for all water supply targets during the allocation period based on the outflow of water from the reservoir and the water loss in the river channels of the basin. Specifically, this includes: A surface water resource allocation model for the watershed is constructed based on the water demand of all water-supplying entities. Based on the outflow of water from the reservoir and the water loss in the river channels of the basin, the total water supply guarantee rate for all water supply targets during the allocation period is maximized using the surface water resource allocation model.

5. The method according to claim 4, characterized in that, A surface water resource allocation model is constructed based on the water demand of all water-supplying entities, specifically including: A second objective function is constructed based on the water demand of all water supply targets. The second objective function is used to maximize the total water supply guarantee rate of all water supply targets during the allocation period. A second constraint is constructed based on the water demand of all water-supplying entities. This second constraint is used to constrain the water resource allocation process. The surface water resource allocation model is constructed based on the second objective function and the second constraint.

6. The method according to claim 1, characterized in that, S3 optimizes the minimization of the second total water shortage in the watershed during the allocation period after groundwater extraction, based on the outflow and supply volumes. Specifically, this includes: Construct a groundwater resource allocation model for the aforementioned watershed; Based on the outflow and supply of water, the second total water shortage in the basin during the allocation period after groundwater extraction is minimized using the groundwater resource allocation model.

7. The method according to claim 6, characterized in that, Constructing a groundwater resource allocation model for the aforementioned watershed specifically includes: The groundwater level in different river sections of the basin after groundwater extraction is determined based on the outflow and supply volume. The amount of water extracted from different river sections after groundwater extraction is determined based on the groundwater level. The groundwater resource allocation model is constructed based on the groundwater level and the amount of water seized.

8. The method according to claim 7, characterized in that, Determining the groundwater level in different river sections of the basin after groundwater extraction based on the outflow and supply volumes specifically includes: The water level of different river sections in the basin is determined based on the outflow and supply of water. The groundwater level in different river sections after groundwater extraction is determined based on the river water level.

9. The method according to claim 7, characterized in that, The amount of water captured from different river sections after groundwater extraction is determined based on the aforementioned groundwater level, specifically including: The exchange flow rate between the river channel and the underground aquifer after groundwater extraction is determined based on the groundwater level. The amount of water extracted from different river sections after groundwater extraction is determined based on the exchange flow rate.

10. A water resource allocation system based on the impact of groundwater extraction along a river, the system comprising: The outflow water calculation module is used to minimize the first total water shortage in the basin before groundwater extraction during the allocation period based on the design parameters of all reservoirs in the basin, and to obtain the outflow water of each reservoir in each time period during the allocation period. The water supply calculation module is used to maximize the overall water supply guarantee rate of all water supply objects during the allocation period based on the outflow of water and the water loss in the river channels of the basin, so as to obtain the water supply volume of the basin to each water supply object in each time period during the allocation period. The river loss correction module is used to minimize the second total water shortage in the basin during the allocation period after groundwater extraction based on the outflow and supply of water, to obtain the amount of water extracted from the basin by groundwater extraction, and to correct the amount of water lost in the river channel based on the amount of water extracted. The joint allocation module is used to iteratively calculate the water supply and water capture volume, determine the surface water-groundwater joint allocation scheme based on the obtained multiple sets of water supply and water capture volume, and allocate water resources based on the surface water-groundwater joint allocation scheme.

Images