CN110543688B - Hydrological model confluence method of multi-blocking two-dimensional diffusion wave equation under reservoir group regulation - Google Patents
Hydrological model confluence method of multi-blocking two-dimensional diffusion wave equation under reservoir group regulation Download PDFInfo
- Publication number
- CN110543688B CN110543688B CN201910716827.5A CN201910716827A CN110543688B CN 110543688 B CN110543688 B CN 110543688B CN 201910716827 A CN201910716827 A CN 201910716827A CN 110543688 B CN110543688 B CN 110543688B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- regulation
- grid
- water surface
- storage
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000009792 diffusion process Methods 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 55
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 123
- 238000003860 storage Methods 0.000 claims abstract description 96
- 238000011144 upstream manufacturing Methods 0.000 claims abstract description 74
- 230000008859 change Effects 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000011160 research Methods 0.000 claims abstract description 13
- 230000001133 acceleration Effects 0.000 claims description 12
- 230000005484 gravity Effects 0.000 claims description 12
- 230000000694 effects Effects 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 22
- 238000004088 simulation Methods 0.000 abstract description 11
- 230000006870 function Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 2
- 230000002457 bidirectional effect Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 235000019687 Lamb Nutrition 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000007427 paired t-test Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000003998 snake venom Substances 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Strategic Management (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Data Mining & Analysis (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a hydrological model confluence method of a multi-blocking two-dimensional diffusion wave equation under reservoir group regulation, which comprises the following steps: (1) Collecting related data in a research area, and determining a dispatching rule, a regulation and storage node and an upstream water collecting area range of a reservoir; (2) Calculating the change of the corresponding water surface elevation of the upstream grid of the reservoir caused by the change of the reservoir water storage capacity according to the scheduling rule and the regulation and storage node of the large reservoir, and constructing a two-dimensional diffusion wave equation set of the upstream grid of the reservoir; (3) According to the dispatching rule and the regulation and storage node of the large reservoir, a two-dimensional diffusion wave equation set of a downstream grid of the reservoir is constructed; (4) According to the regulation nodes and the dispatching rules of the small reservoir, a two-dimensional diffusion wave equation set of the grid where the small reservoir is located is constructed; (5) And performing basin confluence calculation on the equation set to obtain a final confluence result. The invention can improve the model structure of the distributed hydrological model, thereby realizing the simulation of the confluence process under the regulation and storage of the reservoir group and improving the simulation precision of the hydrological model.
Description
Technical Field
The invention relates to the technical field of hydrological model confluence, in particular to a hydrological model confluence method of a multi-blocking two-dimensional diffusion wave equation under reservoir group regulation.
Background
The reservoir is a tool for scientifically, actively and reasonably developing, utilizing and storing water resources and changing regional water circulation rules by human beings. Wherein, the number of large and medium-sized reservoirs (hereinafter collectively referred to as large reservoirs) is relatively small, but the regulation and storage capacity is strong, and the hydrological effect is obvious; the small reservoirs are widely and densely distributed, but the regulation capacity and the hydrological effect are relatively weak, and the reservoir group is the combination of the two.
According to the 'nature-society' binary water resource theory, the watershed reservoir groups can obviously influence the cyclic conversion speed between surface water and underground water, and change the rules of hydrologic cyclic evolution and updating. As an important tool for researching hydrologic cycle, the hydrologic model can simulate the space-time distribution of watershed water resources through modules such as runoff generation, confluence, underground water and the like. Limited by the limitation of computational power, early hydrological models are mostly based on lumped or semi-distributed structures, and simplified convergence calculation methods such as unit lines, linear reservoirs, ma Sijing are adopted. With the development of computer technology, fully distributed models are gradually developed, and such models usually use a river course grid as a computing unit, and perform convergence calculation directly through a discrete one-dimensional saint-venin equation set or a simplified form thereof.
In recent years, the concept of traditional river channel convergence is abandoned by a partial distributed hydrological model, and convergence calculation is carried out on a full-basin grid by adopting a two-dimensional diffusion wave equation. However, due to the regulation of the reservoir group, the hydraulic connection between the upstream and downstream of the natural river is broken, the river flow no longer has continuity, and the two-dimensional diffusion wave equation cannot be applied here. Under the background, because a current hydrological model confluence calculation method of two-dimensional diffusion waves under reservoir group regulation is lacked, the influence of the reservoir group regulation on a confluence process is ignored by the model, the hydrological effect of a hydraulic engineering group cannot be described, and the interaction between a reservoir group consisting of large, medium and small reservoirs and water circulation is difficult to reflect on a physical mechanism.
Disclosure of Invention
The purpose of the invention is as follows: in order to overcome the defects of the prior art, the invention provides a hydrological model confluence method for blocking a two-dimensional diffusion wave equation more under reservoir group regulation.
The technical scheme is as follows: the invention relates to a hydrological model confluence method of a multi-blocking two-dimensional diffusion wave equation under reservoir group regulation, which comprises the following steps:
(1) Collecting the storage capacity, longitude and latitude coordinates, vector river network files and scheduling rule data of a large reservoir in a research area, establishing grid point coordinates, and determining scheduling rules, regulation and storage nodes and upstream catchment area ranges of the large reservoir; collecting the storage capacity and longitude and latitude coordinates of the small reservoir in the research area, establishing grid point coordinates, and determining a scheduling rule and a regulation and storage node of the small reservoir;
(2) Calculating the change of the corresponding water surface elevation of the upstream grid of the reservoir caused by the change of the reservoir water storage capacity according to the scheduling rule and the regulation and storage node of the large reservoir, and constructing a two-dimensional diffusion wave equation set of the upstream grid of the reservoir according to the change of the corresponding water surface elevation;
(3) Determining the inflow rate of the reservoir according to the regulation and storage nodes of the large reservoir, and calculating the discharge rate of the reservoir by using the dispatching rule of the large reservoir, wherein the dispatching rule can adopt the collected actual dispatching rule or the conceptualized dispatching rule obtained by generalization. Constructing a two-dimensional diffusion wave equation set of the downstream grids of the reservoir according to the above;
(4) Determining the inflow rate of the reservoir according to the regulation nodes of the small reservoir, calculating the regulation and storage effect of the small reservoir on the runoff by using the dispatching rules of the large reservoir, and generalizing the regulation and storage effect into a source and sink item of a grid where the small reservoir is located, wherein the dispatching rules can adopt collected actual dispatching rules or generalized conceptual dispatching rules. Constructing a two-dimensional diffusion wave equation set of the grid where the small reservoir is located;
(5) And (5) performing basin confluence operation on the equation sets from the step (2) to the step (4) to obtain a final confluence result.
Further, comprising:
in the step (2), calculating the change of the corresponding water surface elevation of the upstream grid of the reservoir caused by the change of the reservoir water storage capacity, and expressing the change as follows:
wherein O is the outlet flow of the large reservoir, I is the inlet flow of the reservoir, A w Dt is the time step for the upstream grid surface area.
Further, comprising:
in the step (2), a two-dimensional diffusion wave equation set of the upstream grid of the reservoir is constructed, namely:
wherein H 'is the water surface elevation of the upstream grid under reservoir regulation, H' = H + delta H, u and v are the flow velocity in the x direction and the y direction respectively, H is the water surface elevation of the upstream grid under the natural condition,
the partial derivative of the water surface elevation of the upstream grid under reservoir regulation in the time direction is obtained,
the partial derivative of the water surface elevation of the upstream grid under reservoir regulation in the x direction is obtained,
the partial derivative of the elevation of the water surface of the upstream grid under reservoir regulation in the y direction is shown, D is the water depth of the upstream grid, g is the gravity acceleration, and K is the flow modulus.
Further, it includes:
in the step (3), a two-dimensional diffusion wave equation set of the downstream grid of the reservoir is constructed and expressed as:
wherein, V x 、V y The flow velocities of the large reservoir discharge flow in the x direction and the y direction respectively, u and v are the flow velocities in the x direction and the y direction respectively, A w Is the water surface area of the downstream grid, H is the water surface elevation of the downstream grid under the natural condition,
the partial derivative in time of the downstream grid water surface elevation in the natural situation,
the partial derivative in the x direction of the downstream grid water surface elevation in the natural situation,
the method is characterized in that the partial derivative of the elevation of the water surface of a downstream grid in the y direction under the natural condition is shown, D is the water depth of the downstream grid, g is the gravity acceleration, K is the flow modulus, and O is the outflow of a large reservoir.
Further, comprising:
in the step (4), the regulation and storage function of the small reservoir to the runoff is calculated according to the regulation and storage node and the dispatching rule of the small reservoir, and is generalized to the source and sink item of the grid where the small reservoir is located, and the source and sink item is expressed as:
wherein S is a source exchange item caused by reservoir regulation, O ' is the outlet flow of the small reservoir, I ' is the inlet flow of the small reservoir, A ' w The area of the water surface of the grid where the small reservoir is located.
Further, it includes:
in the step (4), a two-dimensional diffusion wave equation set of the grid where the small reservoir is located is constructed, and the method comprises the following steps:
wherein H' is the water surface elevation of the grid under the regulation and storage of the small reservoir,
the partial derivative of the grid water surface elevation in the time direction under the regulation and storage of the small reservoir,
the partial derivative of the grid water surface elevation in the x direction under the regulation and storage of the small reservoir,
the partial derivative of the grid water surface elevation in the y direction under the regulation and storage of the small reservoir is obtained; d ' is the water depth of the grid where the small reservoir is located, g is the gravity acceleration, K is the flow modulus, and u ' and v ' are the flow velocities in the x direction and the y direction respectively; v' x ，V' y The flow rates of the small reservoir discharge in the x direction and the y direction are respectively, and A is the area of the water cross section.
Has the advantages that: the model structure of the distributed hydrological model can be perfected by constructing the hydrological model convergence method of the multi-blocking two-dimensional diffusion wave equation under the regulation and storage of the reservoir group, so that the simulation of the convergence process under the regulation and storage of the reservoir group is realized, the simulation precision of the hydrological model is improved, the regional hydrological process under the influence of the reservoir group is quantized, the evolution rule of the regional hydrological process is revealed, theoretical guidance is provided for the development and utilization of regional water resources, and the method has important practical value.
Drawings
FIG. 1 is a schematic flow diagram of a hydrological model confluence method of a multi-blocking two-dimensional diffusion wave equation under reservoir group regulation;
FIG. 2 is the observed value of the monthly runoff, the simulated value of the model under the condition of no reservoir, and the simulated value of the model considering the process of regulating and storing the runoff of the reservoir group in the Ganjiang river upstream river of 1978-1986.
Detailed Description
(1) Collecting reservoir group data of a research area: collecting the storage capacity and longitude and latitude coordinates of a large reservoir and a small reservoir in a research area, and digital elevation models and vector river network data of a drainage basin from a water conservancy department respectively, and determining a regulation node of the reservoir in a hydrological model and a catchment area range of the upstream of the reservoir, wherein if the longitude and latitude coordinate data are incomplete, the longitude and latitude positions of the reservoir can be determined by using remote sensing images; collecting the scheduling rule data of the reservoir group, if the data are not complete, determining the specific position of the reservoir by using the remote sensing image, and approximating the actual scheduling rule of the reservoir by using the conceptual scheduling rule or the optimized scheduling rule.
(2) Dividing sub-basins according to a large reservoir: and converting the longitude and latitude coordinates of the large reservoir into grid point coordinates of a hydrological model according to the collected longitude and latitude coordinates of the large reservoir and the vector river network file, ensuring that the reservoir is positioned in the river network, and determining the regulation and storage nodes of the large reservoir. And processing the watershed digital elevation model, determining the boundary longitude and latitude coordinates of each reservoir water collection area (namely the upstream watershed of the reservoir regulation and storage node), converting the boundary longitude and latitude coordinates into grid point coordinates of the hydrological model, and independently dividing the boundary longitude and latitude coordinates into sub watersheds of the research area in the hydrological model and coding the sub watersheds. Keeping the upstream watershed and the downstream watershed of the reservoir relatively independent, and taking the reservoir as a coupling node of the upstream watershed and the downstream watershed;
(3) Constructing a two-dimensional diffusion wave equation of regulation and storage influence upstream of the large reservoir: calculating the change of the corresponding water surface elevation of the upstream grid of the reservoir caused by the change of the reservoir water storage capacity according to the collected large reservoir dispatching rules and the storage regulation nodes, and constructing a two-dimensional diffusion wave equation set of the upstream grid of the reservoir according to the assumption that the momentum is unchanged; the method comprises the following specific steps:
calculating the change of the corresponding water surface elevation of the upstream grid of the reservoir caused by the change of the reservoir water storage capacity, and expressing as:
wherein O is the output flow of the large reservoir and is determined by reservoir dispatching rules, I is the input flow of the reservoir, and A w Dt is the time step for the upstream grid surface area.
Further, comprising: constructing a two-dimensional diffusion wave equation set of the reservoir upstream grid, namely:
wherein H 'is the water surface elevation of the upstream grid under reservoir regulation, H' = H + delta H, u and v are the flow velocity in the x direction and the y direction respectively, H is the water surface elevation of the upstream grid under the natural condition,
the partial derivative of the water surface elevation of the upstream grid under reservoir regulation in the time direction is obtained,
the partial derivative of the water surface elevation of the upstream grid under reservoir regulation in the x direction is obtained,
and D is the partial derivative of the water surface elevation of the upstream grid under reservoir regulation in the y direction, D is the water depth of the upstream grid, g is the gravity acceleration, and K is the flow modulus.
(4) Constructing a two-dimensional diffusion wave equation of the large reservoir regulation influence downstream: determining the inlet flow of the reservoir according to the regulation and storage nodes of the large reservoir, calculating the discharge flow of the reservoir by using the dispatching rule of the large reservoir, and taking the discharge flow as a source exchange item of the downstream grid of the reservoir, wherein the dispatching rule can adopt the collected actual dispatching rule or the generalized conceptual dispatching rule. And constructing a two-dimensional diffusion wave equation set of the downstream grid of the reservoir according to the equation, wherein the equation comprises a continuity equation and a momentum equation, namely:
wherein, V x 、V y The flow rates of the large reservoir discharge in the x direction and the y direction are respectively, the free discharge is assumed, the calculation is carried out according to the free discharge formula, u and v are the flow rates in the x direction and the y direction respectively, A w Is the water surface area of the downstream grid, H is the water surface elevation of the downstream grid under the natural condition,
the partial derivative in time of the downstream grid water surface elevation in the natural situation,
the partial derivative in the x direction of the downstream grid water surface elevation in the natural situation,
the partial derivative of the elevation of the water surface of the downstream grid in the y direction under the natural condition is shown, D is the water depth of the downstream grid, g is the gravity acceleration, K is the flow modulus, and O is the outflow of the large reservoir.
(5) Constructing a two-dimensional diffusion wave equation under the influence of regulation and storage of the small reservoir group: and converting the longitude and latitude coordinates of the small reservoir into grid point coordinates of a hydrological model according to the collected longitude and latitude coordinates of the small reservoir, and determining a regulation and storage node of the reservoir in the model. Determining the inflow rate of the reservoir according to the regulation nodes of the small reservoir, calculating the regulation and storage effect of the small reservoir on the runoff by utilizing the dispatching rule of the large reservoir, and approximately generalizing the regulation and storage effect into a source collection item of a grid where the small reservoir is located, wherein the dispatching rule can adopt the collected actual dispatching rule or a conceptualized dispatching rule obtained by generalization. In the process, the continuity of water flow is maintained, and a two-dimensional diffusion wave equation set of the grid where the reservoir is located is constructed according to the continuity; the method comprises the following specific steps:
calculating the regulation and storage function of the small reservoir to the runoff, and approximating the regulation and storage function to a source and sink item of a grid where the regulation and storage function is located:
wherein S is a source sink item caused by reservoir regulation, O ' is the outlet flow of the small reservoir, I ' is the inlet flow of the small reservoir, A ' w The area of the water surface of the grid where the small reservoir is located. When the inflow of the reservoir is larger than the outflow, a part of runoff in the period is transferred out of the watershed confluence process; when the outflow of the reservoir is larger than the inflow, a part of runoff in the period is additionally added into the watershed converging process.
And maintaining the continuity of the water flow in the process, and constructing a corresponding two-dimensional diffusion wave equation system according to the continuity equation and the momentum equation, wherein the two-dimensional diffusion wave equation system comprises the following steps:
wherein H' is the elevation of the grid water surface under the regulation and storage of the small reservoir,
the partial derivative of the grid water surface elevation in the time direction under the regulation and storage of the small reservoir,
the partial derivative of the grid water surface elevation in the x direction under the regulation and storage of the small reservoir,
the partial derivative of the grid water surface elevation in the y direction under the regulation and storage of the small reservoir is obtained; d ' is the water depth of the grid where the small reservoir is located, g is the gravity acceleration, K is the flow modulus, u ' and V ' are the flow velocity in the x and y directions, V ' respectively ' x ，V' y The flow rates of the small reservoir discharge in the x direction and the y direction are respectively, and A is the area of the water cross section.
(6) Solving a multi-blocking two-dimensional diffusion wave equation under reservoir group regulation: in the hydrological model, in order to ensure the stability of the solution, the two-dimensional diffusion wave equation set can be respectively solved in each sub-flow domain determined in (2) by adopting an explicit difference with a time step of about 5 minutes or an implicit difference with a time step of half an hour. When the later equation is adopted, the equation set becomes a group of n linear equations of the next time step, wherein n is the grid number, then sparse matrix solving methods such as a bidirectional conjugate gradient iteration method and the like can be adopted for solving, basin convergence calculation is carried out, and a final convergence result is obtained. Therefore, a hydrological model confluence method of a multi-blocking two-dimensional diffusion wave equation under reservoir group regulation is established.
To judge the effectiveness of the present invention, the following experiments were made:
the invention provides a hydrological model confluence method of a multi-blocking two-dimensional diffusion wave equation under regulation and storage of a reservoir group, which comprises the steps of firstly dividing an upstream basin and a downstream basin of a large reservoir into relatively independent sub-basins according to regulation and storage nodes of the large reservoir according to the characteristics of small number and strong regulation and storage capacity of the large reservoir, and calculating the changes of upstream water surface elevation and lower discharge quantity of the reservoir by utilizing a scheduling rule to construct the multi-blocking two-dimensional diffusion wave equation under regulation and storage of the large reservoir; and secondly, aiming at the characteristics that the small reservoirs are large in number, high in density and difficult to obtain in water collecting areas, the regulation and storage effects of the reservoirs on the runoff are calculated according to the regulation and storage nodes and the regulation and storage rules, a two-dimensional diffusion wave equation under the regulation and storage of the small reservoirs is constructed, and the two-dimensional diffusion wave equation are used for the confluence calculation of hydrological models under the regulation and storage of the reservoir groups. The invention effectively makes up the defect that the existing hydrological model adopting two-dimensional diffusion wave equation to carry out confluence calculation cannot simulate the confluence process under reservoir group regulation, and can perfect the model structure of the distributed hydrological model, thereby realizing the simulation of the confluence process under reservoir group regulation, improving the simulation precision of the hydrological model, and having important values in the aspects of revealing the influence rule of the reservoir group on hydrological circulation and guiding the development and utilization of regional water resources.
In this embodiment, the drainage basin at the upstream of the Jiangxiang river is taken as a research area, and is dispersed into 1188 grids, a distributed Hydrologic Model System (HMS) of the research area in a natural state is established, and the Model performs convergence calculation through a two-dimensional diffusion wave equation. A reservoir group consisting of 8 large reservoirs and 432 small reservoirs is built in the watershed, and the runoff process of the watershed outlet section is obviously influenced by the storage regulation function of the reservoir group. Therefore, in order to quantify the influence of the reservoir group on the runoff process and improve the simulation precision, the hydrological model convergence method of the multi-blocking two-dimensional diffusion wave equation under reservoir group regulation is introduced into the HMS, and the model convergence calculation under reservoir group regulation is realized.
The embodiment comprises the following steps:
s1, collecting reservoir group data of a research area: collecting HYDO1K DEM digital elevation models of Ganjiang upstream drainage basin from USGS; the reservoir capacity of 440 reservoirs in the research area, the scheduling rule data of 8 large-scale reservoirs in the research area and the vector river network data of the river basin are collected from the Water and interest agency in Jiangxi province.
Due to the lack of longitude and latitude coordinates and scheduling rule data of the small reservoir, water surface information is extracted by utilizing Landset-8 satellite remote sensing images, longitude and latitude coordinates of 432 small reservoirs are determined, and the actual scheduling rule of the reservoir is replaced by the conceptual scheduling rule. In this embodiment, the following conceptual scheduling rule is adopted:
wherein Q is t For reservoir discharge, U is the human water demand, determined by water resources communique, V t For the current reservoir storage capacity, V d 、V c Respectively a dead storage capacity and a Xingli storage capacity, wherein 20 percent and 65 percent of the storage capacity are respectively taken, and r is a coefficient;
s2, dividing sub-basins according to the large reservoir: and (3) converting the longitude and latitude coordinates of the reservoir into grid point coordinates of an HMS model according to the collected longitude and latitude coordinates and vector river network files of the large reservoir in the step (1), ensuring that the reservoir is positioned in the river network, and determining the regulation and storage nodes of the large reservoir. The method comprises the steps of analyzing a HYDRO1K DEM digital elevation model by using an ARCGIS, determining boundary longitude and latitude coordinates of 8 large reservoir water collection areas (namely upstream drainage areas of reservoir regulation and storage nodes), converting the boundary longitude and latitude coordinates into grid point coordinates of HMS by using a lamb set subprogram, dividing 8 sub drainage areas in a hydrological model together and independently, and coding. Keeping the upstream watershed and the downstream watershed of 8 large-scale reservoirs relatively independent, and taking the reservoirs as coupling nodes of the upstream watershed and the downstream watershed;
s3, constructing a two-dimensional diffusion wave equation of the regulation and storage influence upstream of the large reservoir: calculating the change of the corresponding water surface elevation of the upstream grid of the reservoir caused by the change of the reservoir water storage capacity according to the large reservoir dispatching rule and the regulation node collected in the step (1), and assuming that the momentum is unchanged, and constructing a two-dimensional diffusion wave equation set of the upstream grid of the reservoir according to the change of the corresponding water surface elevation; the method comprises the following specific steps:
calculating the change of the corresponding water surface elevation of the upstream grid of the reservoir caused by the change of the reservoir water storage capacity, and expressing as:
wherein O is macrowaterThe reservoir outlet flow is determined by reservoir dispatching rules, I is the inlet flow of the reservoir, A w Obtaining the water surface area of the upstream grid from an HMS model database, wherein dt is a time step length, and taking half an hour in the HMS model.
Further, comprising: constructing a two-dimensional diffusion wave equation set of the upstream grid of the reservoir, namely:
wherein H 'is the water surface elevation of the upstream grid under reservoir regulation, H' = H + delta H, u and v are the flow velocity in the x direction and the y direction respectively, H is the water surface elevation of the upstream grid under the natural condition,
the partial derivative of the water surface elevation of the upstream grid under reservoir regulation in the time direction is obtained,
the partial derivative of the water surface elevation of the upstream grid under reservoir regulation in the x direction is obtained,
and D is the partial derivative of the water surface elevation of the upstream grid under reservoir regulation in the y direction, D is the upstream water depth, g is the gravity acceleration, and K is the flow modulus. The variables are all from the HMS model database or calculated by the model.
S4, constructing a two-dimensional diffusion wave equation of the regulation and storage influence downstream of the large reservoir: determining the inflow rate of the reservoir according to the large reservoir regulation and storage nodes collected in the step (1), calculating the discharge rate of the reservoir according to the scheduling rules collected in the step (1), and taking the discharge rate as a source and sink item of the downstream grids of the reservoir, thereby constructing a two-dimensional diffusion wave equation set of the downstream grids of the reservoir, wherein the two-dimensional diffusion wave equation set comprises a continuity equation and a momentum equation, namely the continuity equation and the momentum equation
Wherein, V x 、V y The flow rates of the large reservoir discharge in the x direction and the y direction are respectively, the free discharge is assumed, the calculation is carried out according to the free discharge formula, u and v are the flow rates in the x direction and the y direction respectively, A w Is the water surface area of the downstream grid, H is the water surface elevation of the upstream grid under the natural condition,
the partial derivative of the upstream grid water surface elevation in the time direction in the natural situation,
the partial derivative in the x direction of the upstream grid water surface elevation in the natural situation,
the method is characterized in that the partial derivative of the elevation of the water surface of an upstream grid in the y direction under the natural condition is shown, D is the depth of downstream water, g is the gravity acceleration, K is the flow modulus, and O is the outflow of a large reservoir.
S5, constructing a two-dimensional diffusion wave equation under the influence of regulation and storage of the small reservoir group: and (3) converting the longitude and latitude coordinates of the small reservoir into HMS grid point coordinates according to the collected longitude and latitude coordinates of the small reservoir in the step (1), and determining a storage regulation node of the reservoir in the HMS model. Determining the inflow rate of the reservoir according to the obtained regulation nodes, calculating the regulation and storage functions of the small reservoir on the runoff by using the small reservoir scheduling rules collected in the step (1), and approximating and generalizing the regulation and storage functions into the source and sink items of the grid where the small reservoir is located. In the process, the continuity of water flow is maintained, and a two-dimensional diffusion wave equation set of the grid where the reservoir is located is constructed according to the continuity; the method comprises the following specific steps:
calculating the regulation and storage function of the small reservoir to the runoff, and approximating the regulation and storage function to a source and sink item of a grid where the regulation and storage function is located:
wherein S is a source sink item caused by reservoir regulation, O ' is the outlet flow of the small reservoir, I ' is the inlet flow of the small reservoir, A ' w The area of the water surface of the grid where the small reservoir is located. When the inflow of the reservoir is larger than the outflow, a part of runoff in the period is transferred out of the watershed confluence process; when the outflow of the reservoir is larger than the inflow, a part of runoff in the period is additionally added into the watershed converging process.
And maintaining the continuity of the water flow in the process, and constructing a corresponding two-dimensional diffusion wave equation system according to the continuity equation and the momentum equation, wherein the two-dimensional diffusion wave equation system comprises the following steps:
wherein H' is the water surface elevation of the grid under the regulation and storage of the small reservoir,
deviation of grid water surface elevation in time direction under regulation and storage of small reservoirThe number of the first and second groups is counted,
the partial derivative of the grid water surface elevation in the x direction under the regulation and storage of the small reservoir,
the partial derivative of the grid water surface elevation in the y direction under the regulation and storage of the small reservoir is obtained; d ' is the water depth of the grid where the small reservoir is located, g is the gravity acceleration, K is the flow modulus, u ' and V ' are the flow velocity in the x and y directions, V ' respectively ' x ，V' y The flow rates of the small reservoir discharge flow in the x direction and the y direction are respectively, A is the water passing section area, and the rest symbols are the same as above and are obtained from an HMS model database or calculated by a model.
S6, solving a multi-blocking two-dimensional diffusion wave equation under reservoir group regulation: in the distributed hydrological model HMS, in order to ensure the solving stability, the two-dimensional diffusion wave equation set is respectively solved in each sub-flow domain determined in the step (2) by adopting the time implicit difference with the time step length of half an hour. The equation set becomes a set of 1188 linear equations of the next time step, wherein 1188 is the grid number, then a bidirectional conjugate gradient iteration method is adopted to solve, the drainage basin convergence calculation is performed, and the final convergence result of the outlet section of the Ganjiang upstream drainage basin is obtained.
The simulation method effectively makes up the defect that the distributed hydrological model HMS cannot simulate the converging process under the regulation and storage of the reservoir group, and perfects the model structure of the HMS, thereby realizing the simulation of the converging process of the Ganjiang river upstream drainage basin under the regulation and storage of the reservoir group. The method can also obviously improve the simulation precision of the lunar runoff of the model in the Ganjiang upstream drainage basin. As shown in FIG. 2, by considering the influence of reservoir group regulation on the confluence process through the application of the method, the Nash efficiency coefficient (NSE) of the runoff simulation value of the Ganjiang upstream drainage basin in 1978-1986 is improved to 0.97 from 0.93, and the significance level of the paired T test is passed by 95%. In addition, the method has important values in the aspects of revealing the influence rule of the reservoir group on the hydrological cycle of the Ganjiang upstream drainage basin and guiding the development and utilization of water resources of the Ganjiang upstream drainage basin.
It is to be noted that, in this document, relational terms such as first and second, and the like are used solely to distinguish one entity or operation from another entity or operation without necessarily requiring or implying any actual such relationship or order between such entities or operations.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely application embodiment, or an embodiment combining application and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
While preferred embodiments of the present invention have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. Therefore, it is intended that the appended claims be interpreted as including the preferred embodiment and all changes and modifications that fall within the scope of the invention.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.
Claims (3)
1. A hydrological model confluence method of a multi-blocking two-dimensional diffusion wave equation under reservoir group regulation is characterized by comprising the following steps:
(1) Collecting the storage capacity, longitude and latitude coordinates, vector river network files and scheduling rule data of a large reservoir in a research area, establishing grid point coordinates, and determining a scheduling rule, a regulation node and an upstream water collecting area range of the large reservoir; collecting the storage capacity and longitude and latitude coordinates of the small reservoir in the research area, establishing grid point coordinates, and determining a scheduling rule and a regulation and storage node of the small reservoir;
(2) Calculating the change of the corresponding water surface height of the upstream grid of the reservoir caused by the change of the reservoir water storage capacity according to the scheduling rule and the regulation and storage node of the large reservoir, and constructing a two-dimensional diffusion wave equation set of the upstream grid of the reservoir;
(3) Determining the inlet flow of the reservoir according to the regulation and storage nodes of the large reservoir, and calculating the discharge flow of the reservoir by using the scheduling rule of the large reservoir, wherein the scheduling rule adopts the collected actual scheduling rule or the conceptual scheduling rule obtained by generalization, and a two-dimensional diffusion wave equation set of the downstream grid of the reservoir is constructed;
(4) Determining the inflow rate of the reservoir according to the regulation nodes of the small reservoir, calculating the regulation and storage effect of the small reservoir on the runoff by utilizing the dispatching rule of the large reservoir, and generalizing the regulation and storage effect into a source collection item of a grid where the small reservoir is located, wherein the dispatching rule adopts the collected actual dispatching rule or the conceptualized dispatching rule obtained by generalization, and a two-dimensional diffusion wave equation set of the grid where the small reservoir is located is constructed;
(5) Performing basin confluence calculation on the equation sets in the steps (2) to (4) to obtain a final confluence result;
in the step (2), the change of the corresponding water surface elevation of the upstream grid of the reservoir caused by the change of the reservoir water storage capacity is expressed as:
wherein O is the outlet flow of the large reservoir, I is the inlet flow of the reservoir, A w The water surface area of the upstream grid is obtained, and dt is a time step;
in the step (2), a two-dimensional diffusion wave equation set of the upstream grid of the reservoir is constructed, namely:
wherein H 'is the water surface elevation of the upstream grid under reservoir regulation, H' = H + delta H, u and v are the flow velocity in the x direction and the y direction respectively, H is the water surface elevation of the upstream grid under the natural condition,
the partial derivative of the water surface elevation of the upstream grid under reservoir regulation in the time direction is obtained,
the partial derivative of the water surface elevation of the upstream grid under reservoir regulation in the x direction is obtained,
the partial derivative of the elevation of the water surface of the upstream grid under reservoir regulation in the y direction is shown, D is the water depth of the upstream grid, g is the gravity acceleration, and K is the flow modulus;
in the step (3), a two-dimensional diffusion wave equation set of the downstream grid of the reservoir is constructed and expressed as follows:
wherein, V x 、V y The flow velocities of the large reservoir discharge flow in the x direction and the y direction respectively, u and v are the flow velocities in the x direction and the y direction respectively, A w Is the water surface area of the downstream grid, h is the water surface elevation of the downstream grid under the natural condition,
the partial derivative in time of the downstream grid water surface elevation in the natural situation,
for the natural situation that the elevation of the water surface of the downstream grid is xThe partial derivative in the direction of the light,
the method is characterized in that the partial derivative of the elevation of the water surface of a downstream grid in the y direction under the natural condition is shown, D is the water depth of the downstream grid, g is the gravity acceleration, K is the flow modulus, and O is the outflow of a large reservoir.
2. The hydrological model confluence method of multi-blocking two-dimensional diffusion wave equation under reservoir group regulation according to claim 1, wherein in the step (4), the regulation and storage effect of the small reservoir on the runoff is calculated according to the regulation and storage node and the dispatching rule of the small reservoir, and is generalized to the source and sink term of the grid where the small reservoir is located, and the source and sink term is expressed as:
wherein S is a source exchange item caused by reservoir regulation, O ' is the outlet flow of the small reservoir, I ' is the inlet flow of the small reservoir, A ' w Is the water surface area of the small reservoir grid.
3. The hydrological model confluence method of multi-blocking two-dimensional diffusion wave equation under regulation and storage of reservoir group according to claim 1, wherein in the step (4), the construction of the two-dimensional diffusion wave equation set of the grid where the small reservoir is located comprises:
wherein H' is the elevation of the grid water surface under the regulation and storage of the small reservoir,
the partial derivative of the grid water surface elevation in the time direction under the regulation and storage of the small reservoir,
the partial derivative of the grid water surface elevation in the x direction under the regulation and storage of the small reservoir,
the partial derivative of the grid water surface elevation in the y direction under the regulation and storage of the small reservoir is obtained; d ' is the water depth of the grid where the small reservoir is located, g is the gravity acceleration, K is the flow modulus, u ' and V ' are the flow velocity in the x and y directions, V ' respectively ' x ，V' y The flow rates of the small reservoir discharge in the x direction and the y direction are respectively, and A is the area of the water passing section.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910716827.5A CN110543688B (en) | 2019-08-05 | 2019-08-05 | Hydrological model confluence method of multi-blocking two-dimensional diffusion wave equation under reservoir group regulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910716827.5A CN110543688B (en) | 2019-08-05 | 2019-08-05 | Hydrological model confluence method of multi-blocking two-dimensional diffusion wave equation under reservoir group regulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110543688A CN110543688A (en) | 2019-12-06 |
| CN110543688B true CN110543688B (en) | 2022-10-14 |
Family
ID=68710550
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910716827.5A Active CN110543688B (en) | 2019-08-05 | 2019-08-05 | Hydrological model confluence method of multi-blocking two-dimensional diffusion wave equation under reservoir group regulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110543688B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104281780A (en) * | 2014-10-11 | 2015-01-14 | 水利部交通运输部国家能源局南京水利科学研究院 | Linear reservoir retention confluence and nested watershed (multiple sub watersheds) confluence method |
| CN108256193A (en) * | 2018-01-10 | 2018-07-06 | 中国水利水电科学研究院 | A kind of river flood due to dam-break one, two-dimensional numerical model coupling analytical method |
| CN109101706A (en) * | 2018-07-25 | 2018-12-28 | 北京大学 | The coupling process of lumped hydrological model and two-dimentional hydrodynamic model |
| CN109190261A (en) * | 2018-09-07 | 2019-01-11 | 中国水利水电科学研究院 | A kind of flood risk analysis method that one-dimensional river network generalization is coupled with one two dimension of high-performance |
-
2019
- 2019-08-05 CN CN201910716827.5A patent/CN110543688B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104281780A (en) * | 2014-10-11 | 2015-01-14 | 水利部交通运输部国家能源局南京水利科学研究院 | Linear reservoir retention confluence and nested watershed (multiple sub watersheds) confluence method |
| CN108256193A (en) * | 2018-01-10 | 2018-07-06 | 中国水利水电科学研究院 | A kind of river flood due to dam-break one, two-dimensional numerical model coupling analytical method |
| CN109101706A (en) * | 2018-07-25 | 2018-12-28 | 北京大学 | The coupling process of lumped hydrological model and two-dimentional hydrodynamic model |
| CN109190261A (en) * | 2018-09-07 | 2019-01-11 | 中国水利水电科学研究院 | A kind of flood risk analysis method that one-dimensional river network generalization is coupled with one two dimension of high-performance |
Non-Patent Citations (2)
| Title |
|---|
| CMADS-Driven Simulation and Analysis of Reservoir Impacts on the Streamflow with a Simple Statistical Approach;Ningpeng Dong等;《water》;20190121;全文 * |
| 分布式水文模型建模过程研究;夏积德等;《西安文理学院学报：自然科学版》;20081031;第11卷(第4期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110543688A (en) | 2019-12-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | A coupled 1D–2D hydrodynamic model for flood simulation in flood detention basin | |
| Kang et al. | High-resolution numerical simulation of turbulence in natural waterways | |
| Chen et al. | Determination of input for artificial neural networks for flood forecasting using the copula entropy method | |
| Delelegn et al. | Multi-objective optimisation of cost–benefit of urban flood management using a 1D2D coupled model | |
| CN111768502A (en) | Non-structural grid two-dimensional flood simulation system based on GPU acceleration technology | |
| CN104598701B (en) | A kind of river pollution emulation deduction method of three-dimensional scenic segment iteration | |
| CN107315847B (en) | A kind of generation method and device in river and underground water coupled simulation parameter | |
| Kremer et al. | Simulating property exchange in estuarine ecosystem models at ecologically appropriate scales | |
| CN116822400B (en) | One-dimensional unsteady flow simulation method suitable for large-scale plain river network | |
| Wang et al. | Testing a new adaptive grid‐based shallow flow model for different types of flood simulations | |
| Tavelli et al. | A high‐order parallel Eulerian‐Lagrangian algorithm for advection‐diffusion problems on unstructured meshes | |
| Bittinger et al. | Mathematical simulations for better aquifer management | |
| Alarcon et al. | Nested hydrodynamic modeling of a coastal river applying dynamic-coupling | |
| Huang et al. | Modelling flow in the middle and lower Yangtze River, China | |
| Gu et al. | Optimizing operation rules of sluices in river networks based on knowledge-driven and data-driven mechanism | |
| CN110543688B (en) | Hydrological model confluence method of multi-blocking two-dimensional diffusion wave equation under reservoir group regulation | |
| CN103970610A (en) | Method for monitoring node flow of water supply network | |
| CN108256266B (en) | One-dimensional hydrodynamic model and two-dimensional hydrodynamic model coupling method and system | |
| Sun et al. | Analysis of numerical factors affecting large eddy simulation of pollutant diffusion around buildings | |
| Thakur et al. | Exploring CCHE2D and its sediment modelling capabilities | |
| Fang et al. | Diagonal Cartesian method for the numerical simulation of flow and suspended sediment transport over complex boundaries | |
| Jiang et al. | Responses of river runoff to climate change based on nonlinear mixed regression model in Chaohe River Basin of Hebei Province, China | |
| Podhoranyi et al. | A Web-Based Modelling and Monitoring System Based on Coupling Environmental Models and Hydrological-Related Data. | |
| Hu et al. | Three-dimensional simulation of scalar transport in large shallow water systems using flux-form Eulerian–Lagrangian method | |
| Maddaus et al. | A regional groundwater resource management model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |