CN107766671A - A kind of distributed water componental movement Scenario Simulating method - Google Patents
A kind of distributed water componental movement Scenario Simulating method Download PDFInfo
- Publication number
- CN107766671A CN107766671A CN201711086900.2A CN201711086900A CN107766671A CN 107766671 A CN107766671 A CN 107766671A CN 201711086900 A CN201711086900 A CN 201711086900A CN 107766671 A CN107766671 A CN 107766671A
- Authority
- CN
- China
- Prior art keywords
- mrow
- crop planting
- distribution
- area
- mfrac
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 30
- 238000009826 distribution Methods 0.000 claims abstract description 45
- 238000004088 simulation Methods 0.000 claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims abstract description 28
- 239000002689 soil Substances 0.000 claims abstract description 23
- 230000002262 irrigation Effects 0.000 claims description 24
- 238000003973 irrigation Methods 0.000 claims description 24
- 230000004720 fertilization Effects 0.000 claims description 15
- 239000003621 irrigation water Substances 0.000 claims description 10
- 238000009827 uniform distribution Methods 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 240000008042 Zea mays Species 0.000 description 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000005822 corn Nutrition 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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)
- Cultivation Of Plants (AREA)
Abstract
本发明公开了一种分布式水分运动多情景模拟方法，该方法包括获取目标区域的土壤、作物种植结构、地表高程等数据，根据这些数据划分基本计算单元，计算既定目标值，根据基本计算单元采用既定的多种参数分布类型，进行多种参数分布的组合，以及依据区域的作物种植规划及未来适于种植作物发展占比，确定不同作物种植面积占比条件下的组合形式，最后批量进行不同作物种植面积、土壤、气象参数等数据下多情景模拟，得到输出结果，计算多情景条件下的既定目标值；该情景模拟方法运算效率更高，参数分布组合条件更有针对性；同时与已有的作物种植面积占比设置相比，作物种植面积占比更精确，且能够给出作物种植的空间分布形式。
The invention discloses a multi-scenario simulation method for distributed water movement. The method includes obtaining data such as soil, crop planting structure, and surface elevation in a target area, dividing basic calculation units according to these data, and calculating predetermined target values. Use the established multiple parameter distribution types to combine multiple parameter distributions, and determine the combination form under the conditions of different crop planting area proportions according to the crop planting plan of the region and the proportion of crops suitable for planting in the future, and finally carry out in batches Multi-scenario simulation under different crop planting area, soil, meteorological parameters and other data, the output results are obtained, and the established target value under multi-scenario conditions is calculated; the scenario simulation method has higher calculation efficiency and more targeted parameter distribution and combination conditions; Compared with the existing crop planting area proportion setting, the crop planting area proportion is more accurate, and can give the spatial distribution form of crop planting.
Description
Technical Field
The invention relates to the technical field of irrigation and water conservancy, in particular to a multi-scenario simulation method under distributed different planting structures, soils, crop parameters, irrigation fertilization amount and time conditions in an irrigation area, and specifically relates to a distributed water movement multi-scenario simulation method.
Background
The scene simulation method is the most widely applied simulation technology under different crop planting conditions in the current areas and irrigation areas, and in order to design the planting structure distribution of crops, the water delivery and water rotation and irrigation system of the irrigation areas more reasonably and realize high-efficiency water saving, a large number of scenes need to be set to simulate the water circulation process under the conditions of various soil, crops, weather and irrigation and fertilization measures, namely the distributed water movement multi-scene simulation technology; the current distributed water movement multi-scenario simulation method comprises the following steps: dividing distributed basic computing units according to the elevation of the earth surface, the planting structure of crops, the soil type, the meteorological conditions, the distribution of irrigation canal systems and the like, computing the water movement process of each basic computing unit, obtaining water balance elements under the conditions of a plurality of computing units in an area, changing the soil, the meteorological conditions, the irrigation fertilization amount, the time and the crop parameters of the basic computing units, and carrying out scene simulation on the water movement process of each computing unit.
In the process of implementing the present invention, the inventor finds that the existing scene simulation technology has at least the following problems: (1) the setting of different scenes is single, and batch operation and setting under different parameter distribution types cannot be carried out; (2) the selected parameters are relatively fixed, and various subordinate parameters such as soil, weather, crops and the like cannot be combined one by one; (3) the proportion of the crop planting area is random, and the proportion attribution units of various crop areas and the spatial distribution of the units on the area cannot be accurately set. Therefore, a distributed water movement multi-scenario simulation method is provided.
Disclosure of Invention
The invention aims to provide a distributed water movement multi-scenario simulation method which runs in batches, sets and combines a plurality of groups of parameter distribution forms and accurately calculates the planting area and distribution of each crop, so as to solve the problems in the background technology.
In order to achieve the purpose, the invention provides the following technical scheme:
a distributed water movement multi-scenario simulation method comprises the following steps:
step 1: acquiring soil type, crop planting structure, surface elevation (gradient and relative altitude), meteorological parameters, irrigation and fertilization amount and time data of a target area;
step 2: dividing basic calculation units according to the data in the step 1, calculating the moisture movement process of each basic calculation unit under the current condition to obtain an output result, and calculating a set target value based on the output result;
and step 3: determining parameters for setting the scene based on the soil type, the crop planting structure, the surface elevation, the meteorological parameters, the irrigation and fertilization amount and the time data in the basic calculation unit in the step 2, and combining the distribution of various parameters by adopting various established parameter distribution types;
and 4, step 4: based on the basic computing unit determined in the step 2, determining a combination form under the conditions of different crop planting area occupation ratios according to the crop planting plan of the region and the development occupation ratio suitable for planting crops in the future;
and 5: and (4) under the combined conditions of the steps 3 and 4, carrying out multi-scenario simulation under different crop planting areas, soil types, crop planting structures, surface elevations, meteorological parameters, irrigation fertilization amounts and time data in batches to obtain output results, and calculating set target values under the multi-scenario conditions.
Preferably, the outputting result in step 2 includes: yield, evapotranspiration, irrigation water volume and leakage.
Preferably, the step of dividing the basic computing unit in step 2 includes:
the intersection function in arcgis is used for calculation:
a: inputting land utilization type, soil type and earth surface elevation data, selecting ALL (default) for connection attributes, and leaving XY tolerance unfilled (default);
b: output type INPUT (default);
and finally, calculating to obtain a basic calculation unit.
Preferably, the step of calculating the predetermined target value in step 2 includes:
according to the formula
Obtaining the given target value, wherein WUE is irrigation area scale moisture productivity (kg/m) 3 ) Q is the amount of irrigation water (m) 3 )，Y i,j Yield (kg/km) for jth field of ith crop 2 )。A i,j The area (m) of the jth field of the ith crop 2 ) The ith crop has j field pieces.
Preferably, the parameter distribution type in step 3 includes: normal distribution, uniform distribution and gamma distribution, wherein each distribution formula is as follows:
(1) normal distribution:
wherein μ is a mean value; σ is the standard deviation;
(2) uniform distribution:
wherein a and b are respectively a minimum value and a maximum value;
(3) gamma distribution:
wherein, α: a scale parameter; beta: a shape parameter.
Preferably, the parameter distribution type is uniform distribution.
Preferably, the step 4 of determining the combination form under the conditions of different crop planting area ratios comprises the following steps:
the greedy algorithm is used for carrying out combination under different crop planting area proportion conditions, and the specific calculation flow of the greedy algorithm is as follows:
determining HRU and Si with the largest area;
II, determining a maximum demand gap KdSj;
III, distributing Si:
si < KdSj, assigning Si to j,
si > KdSj, and Si is allocated to the largest Si/KSj;
wherein, si is the area of the ith HRU (basic computing unit), KSj is the jth keyword total area demand, and KdSj is the jth keyword area gap.
Compared with the prior art, the invention has the beneficial effects that: the invention provides a distributed water movement multi-scenario simulation method, which is characterized in that a basic calculation unit based on soil types, crop planting structures, surface elevations, meteorological parameters, irrigation fertilization amount and time data is established, parameters of the basic calculation unit participating in scenario setting are determined, distribution combination of various parameters is carried out according to various parameter distribution types, combination forms under different crop planting area ratio conditions are formulated, and water movement simulation is carried out under multi-scenario conditions in batches; compared with the conventional distributed water movement scene simulation method, the method has the advantages that the operation efficiency is higher, and the parameter distribution combination conditions are more targeted; meanwhile, compared with the existing arrangement of the occupied area of the crop planting area, the occupied area of the crop planting area is more accurate, and the spatial distribution form of the crop planting can be given.
Drawings
FIG. 1 is a schematic flow diagram of the present invention;
FIG. 2 is a plot of irrigation water production rate under various scenarios of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, the present invention provides a technical solution:
example 1
A distributed water movement multi-scenario simulation method comprises the following steps:
step 1: acquiring soil type, crop planting structure, surface elevation (gradient and relative altitude), meteorological parameters, irrigation and fertilization amount and time data of a target area;
step 2: dividing the basic calculation units according to the data in the step 1, calculating the moisture movement process of each basic calculation unit under the current condition to obtain an output result, and calculating a set target value based on the output result, wherein the output result comprises: yield, evapotranspiration, irrigation water volume and leakage;
the basic calculation unit for division adopts the intersection function in the arcgis to calculate:
a: inputting land utilization type, soil type and earth surface elevation data, selecting ALL (default) as a connection attribute, and not filling XY tolerance (default);
b: output type INPUT (default);
and finally, calculating to obtain a basic calculation unit.
And step 3: determining parameters for setting the scene based on the soil type, the crop planting structure, the surface elevation, the meteorological parameters, the irrigation and fertilization amount and the time data in the basic calculation unit in the step 2, and combining the distribution of various parameters by adopting various established parameter distribution types;
the parameter distribution type in the step 3 comprises the following steps: normal distribution, uniform distribution and gamma distribution, wherein each distribution formula is as follows:
(1) normal distribution:
wherein μ is a mean value; σ is the standard deviation;
(2) uniform distribution:
wherein a and b are respectively a minimum value and a maximum value;
(3) gamma distribution:
wherein, α: a scale parameter; beta: a shape parameter.
And 4, step 4: based on the basic computing unit determined in the step 2, determining a combination form under the conditions of different crop planting area occupation ratios according to the crop planting plan of the region and the development occupation ratio suitable for planting crops in the future;
the step 4 of determining the combination form under the conditions of different crop planting area ratios comprises the following steps:
the greedy algorithm is used for carrying out combination under different crop planting area proportion conditions, and the specific calculation flow of the greedy algorithm is as follows:
determining HRU and Si with the largest area;
II, determining a maximum demand gap KdSj;
III, distributing Si:
si < KdSj, assigning Si to j,
si > KdSj, and Si is distributed to the largest Si/KSj;
wherein, si is the area of the ith HRU (basic computing unit), KSj is the jth keyword total area demand, and KdSj is the jth keyword area gap.
And 5: and (4) under the combined conditions of the steps 3 and 4, carrying out multi-scenario simulation under different crop planting areas, soil types, crop planting structures, surface elevations, meteorological parameters, irrigation fertilization amounts and time data in batches to obtain output results, and calculating set target values under the multi-scenario conditions.
Wherein the given target value is calculated according to the formula
Obtaining a given target value, wherein WUE is irrigation scale moisture production rate (kg/m) 3 ) Q is the amount of irrigation water (m) 3 )，Y i,j Yield (kg/km) for jth field of ith crop 2 )。A i,j The area (m) of the jth field of the ith crop 2 ) The ith crop has j field pieces.
For example, the soil type is Huang Litu, and the further situation simulation is carried out on corn as planting, and specific parameter data are shown in table 1:
TABLE 1 data parameters of the basic calculation Unit
Wherein the maximum temperature is 24 ℃ and the minimum temperature is 14 ℃ in 26 days in 5 months, no rainfall exists, the obtained irrigation water amount of a basic calculation unit is 786mm, the evaporation amount is 571.33mm, and the yield is 9669.23Kg/hm 2 Irrigation water production rate is 1.23Kg/m 3 The parameters of the specific scenario setup are shown in table 2:
table 2 parameters embodying the scene settings
Based on the determined basic computing unit, according to the crop planting plan of the region and the development occupation ratio of crops suitable for planting in the future, the combination forms under the conditions of different crop planting area occupation ratios are determined, and the specific combination forms under the conditions of different crop planting area occupation ratios are shown in table 3:
TABLE 3 combination of different crop planting area ratios
And finally, carrying out multi-scenario simulation under different crop planting areas, soil types, crop planting structures, surface elevations, meteorological parameters, irrigation fertilization amounts and time data in batches to obtain output results, and calculating a set target value under multi-scenario conditions, wherein the irrigation water production rate under the specific multi-scenario conditions is shown in figure 2.
The invention provides a distributed water movement multi-scenario simulation method, which is characterized in that a basic calculation unit based on soil types, crop planting structures, surface elevations, meteorological parameters, irrigation fertilization amount and time data is established, parameters of the basic calculation unit participating in scenario setting are determined, distribution combination of various parameters is carried out according to various parameter distribution types, combination forms under different crop planting area ratio conditions are formulated, and water movement simulation is carried out under multi-scenario conditions in batches; compared with the conventional distributed water movement scene simulation method, the method has the advantages that the operation efficiency is higher, and the parameter distribution combination conditions are more targeted; meanwhile, compared with the existing arrangement of the occupied area of the crop planting area, the occupied area of the crop planting area is more accurate, and the spatial distribution form of the crop planting can be given.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
Claims (7)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711086900.2A CN107766671A (en) | 2017-11-07 | 2017-11-07 | A kind of distributed water componental movement Scenario Simulating method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711086900.2A CN107766671A (en) | 2017-11-07 | 2017-11-07 | A kind of distributed water componental movement Scenario Simulating method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107766671A true CN107766671A (en) | 2018-03-06 |
Family
ID=61273628
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711086900.2A Pending CN107766671A (en) | 2017-11-07 | 2017-11-07 | A kind of distributed water componental movement Scenario Simulating method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107766671A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108734347A (en) * | 2018-05-14 | 2018-11-02 | 中国水利水电科学研究院 | Waterirrigated productivity scale-transformation method and device |
| CN115062447A (en) * | 2022-04-25 | 2022-09-16 | 中国科学院地理科学与资源研究所 | Method and system for simulating hydrological environment based on ensemble simulation strategy and random layout of planting structure driving SWAT model |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030153466A1 (en) * | 2001-12-21 | 2003-08-14 | Asi Specialities, Ltd. | Phillipsitic zeolite soil amendments |
| CN103645295A (en) * | 2013-12-03 | 2014-03-19 | 中国科学院遥感与数字地球研究所 | Multilayer soil moisture simulation method and multilayer soil moisture simulation system |
| CN106874605A (en) * | 2017-02-22 | 2017-06-20 | 中国水利水电科学研究院 | A kind of gneiss soil Mountainous Area hillside yardstick hydrologic process analogy method |
-
2017
- 2017-11-07 CN CN201711086900.2A patent/CN107766671A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030153466A1 (en) * | 2001-12-21 | 2003-08-14 | Asi Specialities, Ltd. | Phillipsitic zeolite soil amendments |
| CN103645295A (en) * | 2013-12-03 | 2014-03-19 | 中国科学院遥感与数字地球研究所 | Multilayer soil moisture simulation method and multilayer soil moisture simulation system |
| CN106874605A (en) * | 2017-02-22 | 2017-06-20 | 中国水利水电科学研究院 | A kind of gneiss soil Mountainous Area hillside yardstick hydrologic process analogy method |
Non-Patent Citations (2)
| Title |
|---|
| 刘少博: "基于生态水文模型的参数估计算法和灵敏度分析算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 姜瑶: "黑河中游绿洲多尺度农业水文过程及用水效率的模拟分析与优化调控研究", 《中国博士学位论文全文数据库 农业科技辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108734347A (en) * | 2018-05-14 | 2018-11-02 | 中国水利水电科学研究院 | Waterirrigated productivity scale-transformation method and device |
| CN108734347B (en) * | 2018-05-14 | 2020-11-10 | 中国水利水电科学研究院 | Irrigation water productivity scaling method and device |
| CN115062447A (en) * | 2022-04-25 | 2022-09-16 | 中国科学院地理科学与资源研究所 | Method and system for simulating hydrological environment based on ensemble simulation strategy and random layout of planting structure driving SWAT model |
| CN115062447B (en) * | 2022-04-25 | 2023-03-10 | 中国科学院地理科学与资源研究所 | A method and system for simulating hydrological environment based on SWAT model |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7149348B2 (en) | Agricultural irrigation water demand forecast method | |
| Lv et al. | Climate change impacts on regional rice production in China | |
| CN110059362A (en) | Consider the construction method of the regional scale double source evapotranspiration model of vegetation dynamic changes | |
| CN113641946B (en) | Method and device for optimizing planting layout of water-consuming crops | |
| CN114190264B (en) | A method, system and terminal device for determining a precise irrigation scheme | |
| Leite et al. | Distribution of limited irrigation water based on optimized regulated deficit irrigation and typical metheorological year concepts | |
| CN114037351B (en) | Ecosystem service value evaluation model, establishing method and application | |
| CN106717405A (en) | Nutrient Plentiful-lack index regulation coefficient method fertilizing method | |
| CN109544047A (en) | The method for building up of corn soil nutrient Plentiful-lack index system based on height above sea level subregion | |
| CN109522655A (en) | Based on the local ground watering increment calculation method for becoming saturation moisture kinematic system | |
| CN107766671A (en) | A kind of distributed water componental movement Scenario Simulating method | |
| CN109819745A (en) | A method of phosphate fertilizer usage amount is determined based on soil available phosphorus | |
| Zhang et al. | The grain potential of cultivated lands in Mainland China in 2004 | |
| CN108269199A (en) | A kind of small hydro project spatial and temporal distributions formula of object-oriented goes out flow calculation methodologies | |
| CN116167193A (en) | A method for analyzing the impact of land use change on runoff process based on SWAT model | |
| CN108256015A (en) | A kind of Chinese population spatial grid method based on nighttime light data | |
| Liu et al. | Impacts of LUCC processes on potential land productivity in China in the 1990s | |
| CN104831776B (en) | Water network monitoring method | |
| Wu et al. | Land engineering and its role for sustainable agriculture in the agro-pastoral ecotone: A case study of Yulin, Shaanxi Province, China | |
| CN106651011A (en) | Particle swarm algorithm-based canal system optimization water distribution method | |
| CN112016211A (en) | Radish nitrogen fertilizer recommendation method capable of coordinating agriculture and environment | |
| WANG et al. | Comparative study on land use efficiency before and after rural land consolidation in different modes | |
| David et al. | The effect of downward electron heat flow and electron cooling processes in the high-latitude ionosphere | |
| Du et al. | Temporal and spatial changes of blue water and green water in the Taihang Mountain Region, China, in the past 60 years | |
| CN115953064A (en) | Comprehensive treatment and optimized regulation and control method for farmland quality |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180306 |