CN109063306A - A kind of soil bacterial diversity ability spatial spreading method of gridding Hebei Model - Google Patents
A kind of soil bacterial diversity ability spatial spreading method of gridding Hebei Model Download PDFInfo
- Publication number
- CN109063306A CN109063306A CN201810829963.0A CN201810829963A CN109063306A CN 109063306 A CN109063306 A CN 109063306A CN 201810829963 A CN201810829963 A CN 201810829963A CN 109063306 A CN109063306 A CN 109063306A
- Authority
- CN
- China
- Prior art keywords
- grid
- topographic index
- gridding
- topographic
- ability
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to a kind of soil bacterial diversity ability spatial spreading methods of gridding Hebei Model, comprise determining that size of mesh opening, obtain grid water (flow) direction, calculate catchment area, the gradient of each grid is calculated, topographic index is calculated, calculates each grid topographic index difference, establish topographic index cumulative distribution curve, the fitting of same area ratio corresponding lower infiltration ability and topographic index difference.The grid discretization of lower infiltration ability may be implemented in coding method through the invention, provides technical support to establish grid type Hebei Model, with topographic index as conversion intermediary, can react influence of the underlying surface inhomogeneities to stream is produced.
Description
Technical field
The present invention relates to a kind of methods of hydrological model building, more particularly, under a kind of soil of gridding Hebei Model
Infiltration ability spatial spreading method.
Background technique
Hebei Model is to produce on the basis of using for reference many hydrological models production confluence mechanism both at home and abroad in conjunction with Hebei province area
The conceptual lumped hydrological model in semi-dry and semi-humid area that confluence feature is developed.Currently, Hebei Model has become China
One of Main Prediction model in flood forecast system, especially northern China is semiarid, Semi-humid area, since its is good
Forecast precision and be widely applied.Hebei Model has both two kinds of runoff mechanisms of runoff yield under saturated storage and runoff yield excess, and model is by day
Right runoff is divided into two kinds of water sources of rainwash and interflow subsurface drainage, and the condition that rainwash generates is that rainfall intensity is greater than lower seep by force
Degree, the condition that interflow subsurface drainage generates is that the water of lower infiltration meets the water deficit of soil, and extra water can generate interflow subsurface drainage, model
Runoff formation belong to and first surpass infiltration, it is rear store it is full.
The entirety of a homogenization is regarded as in basin by lumped hydrological model, does not consider each parameter in watershed discretization, each water
The space nonuniformity of literary process.For a lumped model, the distribution of rainfall is spatially uniform, perimeter strip
For the geometrical characteristic in part and basin there is no being mutated, variability is not present in the hydrologic process in simulation process.Lumped hydrological model
Since the parameter of model itself is less, reduce the data-requirements for parameter calibration, but due to not considering that hydrologic process is flowing
Variation in domain, thus can not really portray the runoff process of basin each point.
It is independent from each other between each computing unit in gridding Hebei Model, can be calculated separately not in calculating process
With the hydrologic process in computing unit, without interaction between each unit, Runoff calculation is carried out respectively on each grid,
Complete to the processing of computing unit after, then carry out between grid confluence calculation, model can with the size of flexibly adjustment grid,
Realize the discretization of land surface condition in basin.
Calculating of the gridding Hebei Model in each unit still follows the production confluence principle of Hebei Model, and it is full to have both storage
Stream and two kinds of runoff mechanisms of runoff yield excess are produced, therefore how to determine the soil bacterial diversity ability in each unit lattice, directly affect super infiltration
Produce the calculating of stream.The present invention provides a kind of method established by soil relief index with soil bacterial diversity ability relationship, solves net
Lower infiltration ability in each grid of Hebei Model of formatting, solves the technical problem in gridding Hebei Model use process.
Summary of the invention
The present invention devises a kind of soil bacterial diversity ability spatial spreading method of gridding Hebei Model, the technology solved
Problem is to solve for the lower infiltration ability in each grid of gridding Hebei Model, solves in gridding Hebei Model use process
Technical problem.
In order to solve above-mentioned technical problem, present invention employs following scheme:
A kind of soil bacterial diversity ability spatial spreading method of gridding Hebei Model, comprising: step 1 determines size of mesh opening;
Step 2 obtains grid water (flow) direction;Step 3 calculates catchment area;Step 4, the gradient for calculating each grid;Step 5 calculates ground
Shape index;Step 6 calculates each grid topographic index difference;Step 7 establishes topographic index cumulative distribution curve;It is step 8, identical
The fitting of area ratio corresponding lower infiltration ability and topographic index difference.
Further, the catchment area that each grid flows through slope surface any point unit contour length is calculated in the step 3, is used
α expression, calculation formula are as follows:
α=(N+1) A;
In formula: N is the confluence grid number that cell upstream imports the cell;A is single grid area.
Further, in the step 4: the gradient of each grid is calculated, is indicated with β, calculation formula are as follows:
In formula: depth displacement of the Δ h between adjacent mesh;Δ l is the distance at adjacent mesh center, if the water flow side of the grid
To for due east, due south due west or due north, then Δ h is equal to unit grid length, if the water (flow) direction of the grid is northeast, southeast west
North or southwest, then Δ h is equal to 2 times of unit grid length.
Further, the step 5: calculating the topographic index γ of each grid according to the calculated result in step and step 4, meter
Calculate formula are as follows:
Further, the step 6: calculating each grid topographic index difference, i.e., each grid topographic index and the research area are most
Difference between ninor feature index, calculation formula are as follows:
γ′i=γi-γmin
In formula, γ 'iFor the difference of each grid topographic index and the research area minimum topographic index, γiFor the ground of each grid
Shape index, γminFor the research area minimum topographic index.
Further, the step 7: using ordinate as topographic index difference, abscissa is less than a certain topographic index difference
Grid account for the ratio of total grid number, establish the correlativity of topographic index cumulative distribution curve, it may be assumed that
In formula,The ratio of total grid number, γ ' are accounted for for the grid less than a certain topographic index differenceiFor net each in step 6
The difference of lattice topographic index and the research area minimum topographic index.
Further, the step 8: there are the distribution curves of a soil bacterial diversity ability in lump type Hebei Model, i.e.,
Wherein, a is the ratio that drainage area is accounted for less than or equal to the area of a certain infiltration rate, fmFor under the maximum point of basin
Infiltration ability, fiFor the lower infiltration ability of basin certain point, n is coefficient, is obtained by rainfall runoff data calibration lump type Hebei Model
To n and fm;
It selects corresponding on the corresponding lower infiltration capability distribution curve of same area ratio and topographic index cumulative distribution curve
Lower infiltration ability value and topographic index difference, wherein abscissa is topographic index difference, and ordinate is both lower infiltration ability value, fitting
Relationship, may be expressed as:
fj=F (γ 'i,θ);
In formula, θ is parameter when being fitted the relationship of the two, γ 'iMost for grid topographic index each in step 6 and the research area
The difference of ninor feature index, fiIt is γ ' to correspond to difference in basiniGrid lower infiltration ability.
Further, the lower infiltration ability of each grid can be calculated by the grid topographic index difference of each grid.
Further, the step 1: according to the range of survey region and research it needs to be determined that established gridding Hebei Model
Sizing grid, Nm × Nm, N > 0,30m × 30m, 500m × 500m or 1km × 1km, each grid is big in a modeling process
It is small to be consistent.
Further, by the hydrological analysis tool in Arcgis 10.2, basin the step 2: is obtained by single flow direction algorithm
Confluence add up matrix, the number for the accumulative matrix identification that converges is the confluence grid number that each unit lattice upstream imports the cell,
The flow direction of each grid is obtained simultaneously.
The soil bacterial diversity ability spatial spreading method of the gridding Hebei Model has the advantages that
The relationship between topographic index and lower infiltration ability is established in present invention research, realizes lower infiltration ability on the space of basin
It is discrete, provide technical support to establish grid type Hebei Model, with topographic index as conversion intermediary, can certain journey
Influence of the reaction underlying surface inhomogeneities of degree to stream is produced.
Detailed description of the invention
Fig. 1 is grid dividing schematic diagram in the present invention;
Fig. 2 is mesorelief index cumulative distribution curve schematic diagram of the present invention;
Fig. 3 is mesorelief index of the present invention and lower infiltration ability correspondence diagram.
Specific embodiment
Below with reference to Fig. 1 to Fig. 3, the present invention will be further described:
The present embodiment is a kind of soil bacterial diversity ability spatial spreading method of gridding Hebei Model, and process is as shown in Figure 1.
The step of the present embodiment the method, is as follows:
Step 1: according to the range of survey region and the sizing grid it needs to be determined that established gridding Hebei Model is studied,
Such as 30m × 30m, 500m × 500m or 1km × 1km, the size of each grid should be consistent in a modeling process;
Step 2: by the hydrological analysis tool in Arcgis 10.2, the confluence for obtaining basin by single flow direction algorithm is accumulative
Matrix, the number for the accumulative matrix identification that converges are that each unit lattice upstream imports the confluence grid number of the cell, while obtaining each
The flow direction of grid;
The catchment area that each grid flows through slope surface any point unit contour length is calculated in step 3, is indicated with α, is calculated
Formula are as follows:
α=(N+1) A;
In formula: N is the confluence grid number that cell upstream imports the cell;A is single grid area.
Further, in the step 4: the gradient of each grid is calculated, is indicated with β, calculation formula are as follows:
In formula: depth displacement of the Δ h between adjacent mesh;Δ l is the distance at adjacent mesh center, if the water flow side of the grid
To for due east, due south due west or due north, then Δ h is equal to unit grid length, if the water (flow) direction of the grid is northeast, southeast west
North or southwest, then Δ h is equal to 2 times of unit grid length.
Step 5: the topographic index γ of each grid, calculation formula are calculated according to the calculated result in step and step 4 are as follows:
Step 6: calculate each grid topographic index difference, i.e., each grid topographic index and the research area minimum topographic index it
Between difference, calculation formula are as follows:
γ′i=γi-γmin
In formula, γ 'iFor the difference of each grid topographic index and the research area minimum topographic index, γiFor the ground of each grid
Shape index, γminFor the research area minimum topographic index.
Step 7: using ordinate as topographic index difference, abscissa accounts for total net for the grid less than a certain topographic index difference
The ratio of lattice number establishes the correlativity of topographic index cumulative distribution curve, it may be assumed that
In formula,The ratio of total grid number, γ ' are accounted for for the grid less than a certain topographic index differenceiFor each net in 5 steps 6
The difference of lattice topographic index and the research area minimum topographic index.
Step 8: there are the distribution curves of a soil bacterial diversity ability in lump type Hebei Model, i.e.,
Wherein, a is the ratio that drainage area is accounted for less than or equal to the area of a certain infiltration rate, fmFor under the maximum point of basin
Infiltration ability, fiFor the lower infiltration ability of basin certain point, n is coefficient, is obtained by rainfall runoff data calibration lump type Hebei Model
To n and fm;
It selects corresponding on the corresponding lower infiltration capability distribution curve of same area ratio and topographic index cumulative distribution curve
Lower infiltration ability value and topographic index difference, wherein abscissa is topographic index difference, and ordinate is both lower infiltration ability value, fitting
Relationship, may be expressed as:
fj=F (γ 'i,θ)
In formula, θ is parameter when being fitted the relationship of the two, γ 'iMost for grid topographic index each in step 6 and the research area
The difference of ninor feature index, fiIt is γ ' to correspond to difference in basiniGrid lower infiltration ability.
The lower infiltration ability of each grid can be calculated by the grid topographic index difference of each grid.
Fig. 2 is the relationship between the Fuping watershed unit index difference value counted in example and area ratio, is with ordinate
Topographic index difference, abscissa are the ratio that total grid number is accounted for less than the grid of a certain topographic index difference, and established landform refers to
The correlativity of number cumulative distribution curve.
Fig. 3 is that the corresponding lower infiltration capability distribution curve of selection same area ratio counted in example and topographic index are accumulative
The corresponding lower relationship seeped between ability value and topographic index difference, point are painted in reference axis on distribution curve, wherein abscissa
For topographic index difference, ordinate is lower infiltration ability value.
Above in conjunction with attached drawing, an exemplary description of the invention, it is clear that realization of the invention is not by aforesaid way
Limitation, as long as use the inventive concept and technical scheme of the present invention carry out various improvement, or it is not improved will be of the invention
Conception and technical scheme directly apply to other occasions, be within the scope of the invention.
Claims (10)
1. a kind of soil bacterial diversity ability spatial spreading method of gridding Hebei Model, comprising: step 1 determines size of mesh opening;Step
Rapid 2, grid water (flow) direction is obtained;Step 3 calculates catchment area;Step 4, the gradient for calculating each grid;Step 5 calculates landform
Index;Step 6 calculates each grid topographic index difference;Step 7 establishes topographic index cumulative distribution curve;Step 8, identical faces
The fitting of product ratio corresponding lower infiltration ability and topographic index difference.
2. the soil bacterial diversity ability spatial spreading method of gridding Hebei Model according to claim 1, it is characterised in that: institute
It states and calculates the catchment area that each grid flows through slope surface any point unit contour length in step 3, indicated with α, calculation formula are as follows:
α=(N+1) A;
In formula: N is the confluence grid number that cell upstream imports the cell;A is single grid area.
3. the soil bacterial diversity ability spatial spreading method of gridding Hebei Model according to claim 1 or claim 2, feature exist
In: in the step 4: the gradient of each grid is calculated, is indicated with β, calculation formula are as follows:
In formula: depth displacement of the Δ h between adjacent mesh;Δ l is the distance at adjacent mesh center, if the water (flow) direction of the grid is
Due east, due south due west or due north, then Δ h be equal to unit grid length, if the water (flow) direction of the grid be northeast, all directions or
Southwest, then Δ h is equal to 2 times of unit grid length.
4. the soil bacterial diversity ability spatial spreading method of gridding Hebei Model according to claim 3, it is characterised in that: institute
It states step 5: calculating the topographic index γ of each grid, calculation formula according to the calculated result in step and step 4 are as follows:
5. the soil bacterial diversity ability spatial spreading method of gridding Hebei Model according to claim 4, it is characterised in that: institute
It states step 6: calculating each grid topographic index difference, i.e., between each grid topographic index and the research area minimum topographic index
Difference, calculation formula are as follows:
γ'i=γi-γmin
In formula, γ 'iFor the difference of each grid topographic index and the research area minimum topographic index, γiLandform for each grid refers to
Number, γminFor the research area minimum topographic index.
6. the soil bacterial diversity ability spatial spreading method of gridding Hebei Model according to claim 5, it is characterised in that: institute
State step 7: using ordinate as topographic index difference, abscissa accounts for total grid number for the grid less than a certain topographic index difference
Ratio establishes the correlativity of topographic index cumulative distribution curve, it may be assumed that
In formula,The ratio of total grid number, γ ' are accounted for for the grid less than a certain topographic index differenceiFor grid each in step 6
The difference of shape index and the research area minimum topographic index.
7. the soil bacterial diversity ability spatial spreading method of gridding Hebei Model described in any one of -6 according to claim 1,
It is characterized by: the step 8: there are the distribution curves of a soil bacterial diversity ability in lump type Hebei Model, i.e.,
Wherein, a is the ratio that drainage area is accounted for less than or equal to the area of a certain infiltration rate, fmFor the lower infiltration energy of basin maximum point
Power, fiFor the lower infiltration ability of basin certain point, n is coefficient, by rainfall runoff data calibration lump type Hebei Model obtain n and
fm;
Select corresponding lower infiltration on the corresponding lower infiltration capability distribution curve of same area ratio and topographic index cumulative distribution curve
Ability value and topographic index difference, wherein abscissa is topographic index difference, and ordinate is lower infiltration ability value, is fitted the pass of the two
System, may be expressed as:
fj=F (γ 'i,θ);
In formula, θ is parameter when being fitted the relationship of the two, γ 'iFor grid topographic index each in step 6 and the research area minimum landform
The difference of index, fiIt is γ ' to correspond to difference in basiniGrid lower infiltration ability.
8. the soil bacterial diversity ability spatial spreading method of gridding Hebei Model according to claim 7, it is characterised in that: logical
The lower infiltration ability of each grid can be calculated in the grid topographic index difference for crossing each grid.
9. the soil bacterial diversity ability spatial spreading method of gridding Hebei Model described in any one of -8 according to claim 1,
It is characterized by: the step 1: according to the range of survey region and research it needs to be determined that the net of established gridding Hebei Model
Lattice size: Nm × Nm, N > 0, the size of each grid should be consistent in a modeling process.
10. the soil bacterial diversity ability spatial spreading method of gridding Hebei Model described in any one of -9 according to claim 1,
It is characterized by: the step 2: by the hydrological analysis tool in Arcgis 10.2, obtaining the remittance in basin by single flow direction algorithm
Accumulative matrix is flowed, the number for the accumulative matrix identification that converges is the confluence grid number that each unit lattice upstream imports the cell, simultaneously
Obtain the flow direction of each grid.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810829963.0A CN109063306B (en) | 2018-07-25 | 2018-07-25 | Soil infiltration capacity space dispersion method of gridding Hebei model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810829963.0A CN109063306B (en) | 2018-07-25 | 2018-07-25 | Soil infiltration capacity space dispersion method of gridding Hebei model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109063306A true CN109063306A (en) | 2018-12-21 |
CN109063306B CN109063306B (en) | 2020-07-07 |
Family
ID=64836446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810829963.0A Active CN109063306B (en) | 2018-07-25 | 2018-07-25 | Soil infiltration capacity space dispersion method of gridding Hebei model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109063306B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112733344A (en) * | 2020-12-30 | 2021-04-30 | 北京大学 | Method and system for determining artificial drainage basin production flow |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103729483A (en) * | 2012-10-11 | 2014-04-16 | 张驰 | Device and method for simulating and predicting mountain torrent disasters |
CN106033480A (en) * | 2015-03-13 | 2016-10-19 | 西安山脉科技发展有限公司 | Method for forecasting mountain torrents based on improved Xin'anjiang model |
US20160328498A1 (en) * | 2015-05-06 | 2016-11-10 | International Business Machines Corporation | Topology-Based Method to Assess the Environmental Impact of Deicing Sodium Chloride |
CN106202790A (en) * | 2016-07-20 | 2016-12-07 | 中国水利水电科学研究院 | A kind of novel distributed Hebei Model construction method and application thereof |
CN106446394A (en) * | 2016-09-19 | 2017-02-22 | 河海大学 | Method for extracting drainage basin free water storage capacity space distribution by utilizing terrain indexes |
-
2018
- 2018-07-25 CN CN201810829963.0A patent/CN109063306B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103729483A (en) * | 2012-10-11 | 2014-04-16 | 张驰 | Device and method for simulating and predicting mountain torrent disasters |
CN106033480A (en) * | 2015-03-13 | 2016-10-19 | 西安山脉科技发展有限公司 | Method for forecasting mountain torrents based on improved Xin'anjiang model |
US20160328498A1 (en) * | 2015-05-06 | 2016-11-10 | International Business Machines Corporation | Topology-Based Method to Assess the Environmental Impact of Deicing Sodium Chloride |
CN106202790A (en) * | 2016-07-20 | 2016-12-07 | 中国水利水电科学研究院 | A kind of novel distributed Hebei Model construction method and application thereof |
CN106446394A (en) * | 2016-09-19 | 2017-02-22 | 河海大学 | Method for extracting drainage basin free water storage capacity space distribution by utilizing terrain indexes |
Non-Patent Citations (2)
Title |
---|
石朋: "网格型松散结构分布式水文模型及地貌瞬时单位线研究", 《中国优秀博硕士学位论文全文数据库 (博士) 基础科学辑》 * |
魏洁 等: "基于VIC 模型的黄河上游未来径流变化分析", 《水力发电学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112733344A (en) * | 2020-12-30 | 2021-04-30 | 北京大学 | Method and system for determining artificial drainage basin production flow |
Also Published As
Publication number | Publication date |
---|---|
CN109063306B (en) | 2020-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106884405B (en) | Inrush type mountain flood assay method for a kind of Cross Some Region Without Data | |
CN111369059B (en) | Urban waterlogging rapid prediction method and system based on rain and flood simulation coupling model | |
CN113610264B (en) | Refined power grid typhoon flood disaster prediction system | |
CN111914432B (en) | Hydrological forecasting method based on big data | |
CN102289570B (en) | Flood forecast method based on rainfall-runoff-flood routing calculation | |
CN104392111B (en) | Flood Forecasting Method based on water level sample | |
CN106202790B (en) | A kind of distribution Hebei Model construction method and its application | |
Wu et al. | Simulation of soil loss processes based on rainfall runoff and the time factor of governance in the Jialing River Watershed, China | |
CN108182543A (en) | One kind becomes more meticulous grid waterlogging water logging forecasting procedure | |
Zhang et al. | Estimating spatiotemporal variability and sustainability of shallow groundwater in a well-irrigated plain of the Haihe River basin using SWAT model | |
Liu et al. | Parameter sensitivity analysis of the WRF-hydro modeling system for streamflow simulation: A case study in semi-humid and semi-arid Catchments of Northern China | |
CN110459036B (en) | Mountain torrent early warning method based on deep learning | |
CN105160121A (en) | Finite element control based modeling method for distributed hydrological model | |
CN108897940A (en) | The unidirectional couplings method of hydrological distribution model and two-dimentional hydrodynamic model based on rectangular mesh | |
Abrishamchi et al. | Water resources management scenario analysis in Karkheh River Basin, Iran, using WEAP model | |
CN114048432A (en) | Ecological flow calculation method for area without actual measurement runoff data | |
CN109063306A (en) | A kind of soil bacterial diversity ability spatial spreading method of gridding Hebei Model | |
Li et al. | Study of groundwater using visual MODFLOW in the Manas River Basin, China | |
CN116383551A (en) | Vertical and transverse coupling mixed flow production calculation method | |
CN105976103B (en) | A kind of Flood Forecasting Method based on dynamic water storage capacity | |
Teng et al. | Early warning index of flash flood disaster: a case study of Shuyuan watershed in Qufu City | |
Mathi et al. | Spatial Technology for Mapping Suitable Sites for Run-of-River Hydro Power Plants | |
CN113435631B (en) | Flood forecasting method, flood forecasting system, readable storage medium and computing device | |
Malekani et al. | Application of GIS in modeling Zilberchai Basin runoff | |
Kazama et al. | Evaluation of flood control and inundation conservation in Cambodia using flood and economic growth models |
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 |