CN110686862A - Flow process rasterization method based on soil infiltration capacity - Google Patents
Flow process rasterization method based on soil infiltration capacity Download PDFInfo
- Publication number
- CN110686862A CN110686862A CN201911016712.1A CN201911016712A CN110686862A CN 110686862 A CN110686862 A CN 110686862A CN 201911016712 A CN201911016712 A CN 201911016712A CN 110686862 A CN110686862 A CN 110686862A
- Authority
- CN
- China
- Prior art keywords
- flow
- grid
- distribution
- grids
- soil
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 239000002689 soil Substances 0.000 title claims abstract description 45
- 230000008569 process Effects 0.000 title claims abstract description 35
- 238000001764 infiltration Methods 0.000 title claims abstract description 26
- 230000008595 infiltration Effects 0.000 title claims abstract description 26
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 28
- 230000035699 permeability Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 abstract description 16
- 230000007246 mechanism Effects 0.000 abstract description 10
- 238000004519 manufacturing process Methods 0.000 description 8
- 238000005056 compaction Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 244000025254 Cannabis sativa Species 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M10/00—Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Chemical & Material Sciences (AREA)
- Pathology (AREA)
- Remote Sensing (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Dispersion Chemistry (AREA)
- Computer Graphics (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Sewage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a flow process rasterization method based on soil infiltration capacity, which comprises the following steps: step S1, converting the flow process line of the cross section of the basin outlet into a water level process line through a water level-flow relation curve; step S2, determining flow distribution weights in different grids according to the land utilization type and the soil characteristic data; step S3, determining the distribution proportion of the radial flow of different grids, and calculating the distribution proportion of all grids in the watershed according to the flow distribution weight of each grid; and step S4, distributing the flow of the outlet section of the basin into each grid according to the distribution proportion to realize flow rasterization. The method considers the influence of the soil infiltration capacity of different grids of the drainage basin on the yield and the flow, follows the drainage basin yield and confluence mechanism and the physical mechanism of the hydrological process, reduces the problem of calculation errors caused by average distribution, and has the advantages of simple calculation process, good operation and easy popularization.
Description
Technical Field
The invention relates to the technical field of data downscaling, in particular to a flow process rasterization method based on soil infiltration capacity.
Background
The coupling of hydrologic and hydrodynamic predictions is one of the hot problems in hydrologic research. The important link involved in the hydrographic-hydrodynamic coupling method is the data connection between the hydrographic model and the hydrodynamic model, the calculation result of the hydrographic model is that the flow process of the outlet section of the drainage basin is used as the input of the two-dimensional hydrodynamic model, the two-dimensional hydrodynamic model performs calculation based on the grid data, and the calculation can be performed only by converting the hydrographic data of one-dimensional space scale into the two-dimensional space scale, namely rasterization. However, no mature flow data rasterization method exists at present, and a common method is to equally distribute the flow to each grid, and the method has the disadvantage that the spatial heterogeneity of a watershed is not considered, that is, the soil infiltration capacity of each grid is different, so that the processing directly causes a large error in the spatial distribution of the two-dimensional hydrodynamic calculation result.
Disclosure of Invention
The object of the present invention is to solve at least one of the technical drawbacks mentioned.
Therefore, the invention aims to provide a flow process rasterization method based on soil infiltration capacity.
In order to achieve the above object, an embodiment of the present invention provides a flow process rasterization method based on soil infiltration capacity, including the following steps:
step S1, converting the flow process line of the cross section of the basin outlet into a water level process line through a water level-flow relation curve;
step S2, determining flow distribution weights in different grids according to the land utilization type and the soil characteristic data;
step S3, determining the distribution proportion of the radial flow of different grids, and calculating the distribution proportion of all grids in the watershed according to the flow distribution weight of each grid;
and step S4, distributing the flow of the outlet section of the basin into each grid according to the distribution proportion to realize flow rasterization.
Further, in the step S2, the land use type is a water permeability factor; the soil characteristic data includes: soil structure and initial water content.
Further, in step S2, determining the flow distribution weights in different grids includes determining weights representing the soil infiltration capacity of different grids.
Further, in step S3, the calculating an allocation ratio of all the grids in the drainage basin according to the flow rate allocation weight of each grid includes:
wherein, Pi，jCharacterizing a water permeability factor of the grid (i, j); si，jCharacterizing the soil structure of grid (i, j); wi，jIs the initial moisture content of the grid (i, j); i and j respectively represent the serial numbers of the grids in the x and y directions; m and n are the total number of the grids in the x and y directions of the two-dimensional grid; PRi，jThe distribution proportion of net rainfall for different grids.
Further, in the step S4, the distributing the flow rate of the outlet section of the flow field into each grid according to the distribution ratio includes:
EPi，j＝PRi，j×NR，
among them, EPi，jA grid value representing a flow rate; NR is the flow value of the outlet section of the basin at a certain moment. 6
According to the flow process rasterization method based on the soil infiltration capacity, provided by the embodiment of the invention, soil characteristic data (water permeability factors, soil structures and initial water content) are adopted to represent the soil infiltration capacity, the influence of the three factors on the production flow of a drainage basin is fully considered, the influence of the soil infiltration capacity of different grids of the drainage basin on the production flow is considered, the physical mechanism of a drainage basin production convergence mechanism and a hydrological process is followed, the problem of calculation errors caused by average distribution is reduced, and the calculation process is simple, easy to operate and easy to popularize. The method follows a basin production convergence mechanism and a physical mechanism of a hydrological process, reduces the calculation error caused by simple averaging in the traditional method, and the result can be directly applied to the simulation calculation of the two-dimensional hydrodynamic model.
Additional aspects and advantages of the invention will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the invention.
Drawings
The above and/or additional aspects and advantages of the present invention will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
FIG. 1 is a flow chart of a flow process rasterization method based on soil infiltration capacity in accordance with an embodiment of the present invention;
fig. 2 is a schematic diagram of a flow process rasterization method based on soil infiltration capacity according to an embodiment of the present invention.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the drawings are illustrative and intended to be illustrative of the invention and are not to be construed as limiting the invention.
As shown in fig. 1 and 2, a flow process rasterization method based on soil infiltration capacity of an embodiment of the invention includes the following steps:
and step S1, converting the flow process line of the cross section of the outlet of the basin into a water level process line (NR, unit: meter) through a water level-flow relation curve. This is because the flow rate must be converted to the water level for subsequent calculation.
And step S2, determining flow distribution weights in different grids according to the land use type and the soil characteristic data.
In the step, the water permeability factor P is determined with reference to land utilization types, wherein the weights of common land utilization types are ranked in order of forest land < sparse land < shrubbery land < grassland < cultivated land < bare land < water area < construction land, and specific values are determined with reference to roughness ratios of different land utilization types.
The soil property data includes: soil structure S and initial water content W. The soil structure is generally divided into three conditions of compaction, undisturbed state and scarification, specific values are determined by referring to infiltration rates of different soil structures, and the reference values are 1 for compaction, 0.7 for undisturbed state and 0.3 for scarification. The initial water content is determined according to actual conditions, and the percentage can be directly adopted for calculation in practical application.
It should be noted that, determining the flow distribution weights in different grids, i.e. determining the weights representing the soil infiltration capacity of different grids.
Specifically, in this step, reference values of different grid water permeability factors are determined according to land use types, as shown in table 1. The land use types include: woodland, sparsely, brush, grass, cultivated land, bare land, waters and construction land, table 1 was obtained based on application experience and physical experiments. Reference values for the soil structure (compaction value of 1, undisturbed value of 0.7, scarification value of 0.3).
TABLE 1 reference values for water permeability factors of different grids
| Type of land use | Forest land | Dredging field | Shrub forest | Grass land | Cultivation of land | Bare land | Water area | Land for construction |
| Factor of water permeability | 1 | 2 | 3 | 6.5 | 7 | 8.5 | 10 | 10 |
Step S3 is to determine the distribution ratio of the radial flows of different grids, and calculate the distribution ratio of all the grids in the watershed according to the flow distribution weight of each grid.
Specifically, calculating the allocation proportion of all the grids in the basin according to the flow allocation weight of each grid includes:
wherein, Pi，jCharacterizing a water permeability factor of the grid (i, j); si，jCharacterizing the soil structure of grid (i, j); wi，jIs the initial moisture content of the grid (i, j); i and j respectively represent the serial numbers of the grids in the x and y directions; m and n are the total number of the grids in the x and y directions of the two-dimensional grid; PRi，jThe distribution proportion of net rainfall for different grids.
And step S4, distributing the flow of the outlet section of the basin into each grid according to the distribution proportion to realize flow rasterization.
Specifically, the following formula is adopted to distribute the surface runoff process, namely the drainage basin yield (net rainfall) into each grid according to the distribution proportion, and the method comprises the following steps:
EPi，j＝PRi，j×NR，
among them, EPi，jA grid value representing a flow rate; NR is the surface runoff process of the outlet section of the drainage basin at a certain moment, namely drainage basin yield (net rainfall).
According to the flow process rasterization method based on the soil infiltration capacity, provided by the embodiment of the invention, soil characteristic data (water permeability factors, soil structures and initial water content) are adopted to represent the soil infiltration capacity, the influence of the three factors on the production flow of a drainage basin is fully considered, the influence of the soil infiltration capacity of different grids of the drainage basin on the production flow is considered, the physical mechanism of a drainage basin production convergence mechanism and a hydrological process is followed, the problem of calculation errors caused by average distribution is reduced, and the calculation process is simple, easy to operate and easy to popularize. The method follows a basin production convergence mechanism and a physical mechanism of a hydrological process, reduces the calculation error caused by simple averaging in the traditional method, and the result can be directly applied to the simulation calculation of the two-dimensional hydrodynamic model.
In the description herein, references to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention, and that variations, modifications, substitutions and alterations can be made in the above embodiments by those of ordinary skill in the art without departing from the principle and spirit of the present invention. The scope of the invention is defined by the appended claims and equivalents thereof.
Claims (5)
1. A flow process rasterization method based on soil infiltration capacity is characterized by comprising the following steps:
step S1, converting the flow process line of the cross section of the basin outlet into a water level process line through a water level-flow relation curve;
step S2, determining flow distribution weights in different grids according to the land utilization type and the soil characteristic data;
step S3, determining the distribution proportion of the radial flow of different grids, and calculating the distribution proportion of all grids in the watershed according to the flow distribution weight of each grid;
and step S4, distributing the flow of the outlet section of the basin into each grid according to the distribution proportion to realize flow rasterization.
2. The soil infiltration capacity-based flow process rasterization method of claim 1 wherein in said step S2, said land use type is a water permeability factor; the soil characteristic data includes: soil structure and initial water content.
3. The soil infiltration capacity-based flow process rasterization method of claim 1 wherein in said step S2, determining flow distribution weights within different grids includes determining weights characterizing the soil infiltration capacity of different grids.
4. The soil infiltration capacity-based flow process rasterization method of claim 1, wherein in the step S3, the calculating the distribution proportion of all grids in the drainage basin according to the flow distribution weight of each grid includes:
wherein, Pi，jCharacterizing a water permeability factor of the grid (i, j); si，jCharacterizing the soil structure of grid (i, j); wi，jIs the initial moisture content of the grid (i, j); i and j respectively represent the serial numbers of the grids in the x and y directions; m and n are the total number of the grids in the x and y directions of the two-dimensional grid; PRi，jThe distribution proportion of net rainfall for different grids.
5. The soil infiltration capacity-based flow process rasterization method of claim 1 wherein in said step S4, said distributing the flow of watershed outlet sections into each grid according to distribution proportions comprises:
EPi，j＝PRi，j×NR，
among them, EPi，jA grid value representing a flow rate; NR is the flow value of the outlet section of the basin at a certain moment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911016712.1A CN110686862B (en) | 2019-10-24 | 2019-10-24 | Flow process rasterization method based on soil infiltration capacity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911016712.1A CN110686862B (en) | 2019-10-24 | 2019-10-24 | Flow process rasterization method based on soil infiltration capacity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110686862A true CN110686862A (en) | 2020-01-14 |
| CN110686862B CN110686862B (en) | 2021-04-16 |
Family
ID=69114249
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911016712.1A Expired - Fee Related CN110686862B (en) | 2019-10-24 | 2019-10-24 | Flow process rasterization method based on soil infiltration capacity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110686862B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH06123671A (en) * | 1992-10-09 | 1994-05-06 | Nippon Steel Corp | Method of indicating flow trace from flow field data |
| CN102034001A (en) * | 2010-12-16 | 2011-04-27 | 南京大学 | Design method for distributed hydrological model by using grid as analog unit |
| CN103093114A (en) * | 2013-02-05 | 2013-05-08 | 河海大学 | Distributed-type river basin water deficit calculating method based on terrain and soil characteristics |
| CN107092796A (en) * | 2017-04-26 | 2017-08-25 | 沈阳市水利建筑勘测设计院 | A kind of evaluation method of soil watershed rainwash regulating power |
| CN107290129A (en) * | 2017-05-04 | 2017-10-24 | 中国水利水电科学研究院 | A kind of domatic hydraulics model test flow field observation system and method |
| CN108154270A (en) * | 2017-12-25 | 2018-06-12 | 广州地理研究所 | Middle Flood of small drainage area feature is to the response model construction method of changing environment |
| CN109101706A (en) * | 2018-07-25 | 2018-12-28 | 北京大学 | The coupling process of lumped hydrological model and two-dimentional hydrodynamic model |
| CN109190213A (en) * | 2018-08-20 | 2019-01-11 | 中国科学院地理科学与资源研究所 | Rainwash basin Spatial discretization methods based on land use pattern |
-
2019
- 2019-10-24 CN CN201911016712.1A patent/CN110686862B/en not_active Expired - Fee Related
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH06123671A (en) * | 1992-10-09 | 1994-05-06 | Nippon Steel Corp | Method of indicating flow trace from flow field data |
| CN102034001A (en) * | 2010-12-16 | 2011-04-27 | 南京大学 | Design method for distributed hydrological model by using grid as analog unit |
| CN102419788A (en) * | 2010-12-16 | 2012-04-18 | 南京大学 | Method for designing distributed-type hydrographical model based on penetration-storage integrated dynamic runoff yield mechanism |
| CN103093114A (en) * | 2013-02-05 | 2013-05-08 | 河海大学 | Distributed-type river basin water deficit calculating method based on terrain and soil characteristics |
| CN107092796A (en) * | 2017-04-26 | 2017-08-25 | 沈阳市水利建筑勘测设计院 | A kind of evaluation method of soil watershed rainwash regulating power |
| CN107290129A (en) * | 2017-05-04 | 2017-10-24 | 中国水利水电科学研究院 | A kind of domatic hydraulics model test flow field observation system and method |
| CN108154270A (en) * | 2017-12-25 | 2018-06-12 | 广州地理研究所 | Middle Flood of small drainage area feature is to the response model construction method of changing environment |
| CN109101706A (en) * | 2018-07-25 | 2018-12-28 | 北京大学 | The coupling process of lumped hydrological model and two-dimentional hydrodynamic model |
| CN109190213A (en) * | 2018-08-20 | 2019-01-11 | 中国科学院地理科学与资源研究所 | Rainwash basin Spatial discretization methods based on land use pattern |
Non-Patent Citations (3)
| Title |
|---|
| YUELING WANG 等: "Sensitivity Analysis of the Surface Runoff Coefficient of HiPIMS in Simulating Flood Processes in a Large Basin", 《WATER》 * |
| 王莉莉 等: "基于栅格的分布式超渗产流水文模型构建及比较", 《河海大学学报(自然科学版)》 * |
| 都金康 等: "基于栅格的分布式降雨径流模拟系统及应用", 《地理科学》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110686862B (en) | 2021-04-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chen et al. | Land-use change from arable lands to orchards reduced soil erosion and increased nutrient loss in a small catchment | |
| CN108664647B (en) | Basin fine management system of integrated water environment model | |
| Benli et al. | A non-linear model for farm optimization with adequate and limited water supplies: application to the South-east Anatolian Project (GAP) Region | |
| Fan et al. | Comparative evaluation of crop water use efficiency, economic analysis and net household profit simulation in arid Northwest China | |
| Malveira et al. | Hydrological impact of a high-density reservoir network in semiarid northeastern Brazil | |
| Liu et al. | Evaluating potential impacts of land use changes on water supply–demand under multiple development scenarios in dryland region | |
| CN112001010B (en) | Design method of rainwater regulation and storage facility for controlling runoff pollution of flow distribution system | |
| Slatni et al. | Assessing alternate furrow strategies for potato at the Cherfech irrigation district of Tunisia | |
| Guo et al. | Optimizing irrigation schedule in a large agricultural region under different hydrologic scenarios | |
| KhazaiPoul et al. | Multi-objective optimization for interactive reservoir-irrigation planning considering environmental issues by using parallel processes technique | |
| Min et al. | Quantitative simulation of the relationships between cultivated land-use patterns and non-point source pollutant loads at a township scale in Chaohu Lake Basin, China | |
| Li et al. | Simulation-based optimization for spatiotemporal allocation of irrigation water in arid region | |
| Pei et al. | Study on the optimization of staple crops spatial distribution in China under the influence of natural disasters | |
| Wu et al. | Particulate nutrient loss from drylands to grasslands/forestlands in a large-scale highly erodible watershed | |
| Bluemling et al. | Making water productivity operational—A concept of agricultural water productivity exemplified at a wheat–maize cropping pattern in the North China plain | |
| Deng et al. | Spatial matching and flow in supply and demand of water provision services: A case study in Xiangjiang River Basin | |
| CN110686862B (en) | Flow process rasterization method based on soil infiltration capacity | |
| Li et al. | A daily water balance modelling approach for simulating performance of tank-based irrigation systems | |
| CN109190213B (en) | Surface runoff watershed space discretization method based on land utilization type | |
| Zhang et al. | Runoff prediction under different precipitation scenarios based on SWAT model and stochastic simulation of precipitation | |
| Huang et al. | Agricultural water optimization coupling with a distributed ecohydrological model in a mountain-plain basin | |
| Gudissa et al. | Evaluation of surge and cutback flow furrow irrigation systems for pepper (Capsicum annuum) production | |
| Berihun et al. | Predicting runoff and sediment responses to climate-resilient land use and management scenarios | |
| Yigzaw et al. | Land use and land cover impact on probable maximum flood and sedimentation for artificial reservoirs: case study in the Western United States | |
| Xue et al. | Assessing sustainability of agricultural water saving in an arid area with shallow groundwater |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210416 |