CN110490473A - A kind of crop production water footprints measuring method based on soil moisture dynamic balance - Google Patents
A kind of crop production water footprints measuring method based on soil moisture dynamic balance Download PDFInfo
- Publication number
- CN110490473A CN110490473A CN201910784579.8A CN201910784579A CN110490473A CN 110490473 A CN110490473 A CN 110490473A CN 201910784579 A CN201910784579 A CN 201910784579A CN 110490473 A CN110490473 A CN 110490473A
- Authority
- CN
- China
- Prior art keywords
- water
- crop
- soil
- irrigation
- unit
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 199
- 239000002689 soil Substances 0.000 title claims abstract description 132
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000012272 crop production Methods 0.000 title claims abstract description 55
- 230000002262 irrigation Effects 0.000 claims abstract description 80
- 238000003973 irrigation Methods 0.000 claims abstract description 80
- 241001062009 Indigofera Species 0.000 claims abstract description 40
- 230000012010 growth Effects 0.000 claims abstract description 23
- 238000009826 distribution Methods 0.000 claims abstract description 14
- 230000005540 biological transmission Effects 0.000 claims abstract description 13
- 239000003621 irrigation water Substances 0.000 claims abstract description 9
- 238000004519 manufacturing process Methods 0.000 claims abstract description 9
- 239000008400 supply water Substances 0.000 claims abstract description 8
- 230000008020 evaporation Effects 0.000 claims description 27
- 238000001704 evaporation Methods 0.000 claims description 27
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 8
- 238000001556 precipitation Methods 0.000 claims description 8
- 230000005068 transpiration Effects 0.000 claims description 8
- 244000037666 field crops Species 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 208000005156 Dehydration Diseases 0.000 claims description 3
- 230000003698 anagen phase Effects 0.000 claims description 3
- 238000012512 characterization method Methods 0.000 claims description 3
- 230000007812 deficiency Effects 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 230000008595 infiltration Effects 0.000 claims description 3
- 238000001764 infiltration Methods 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 238000009395 breeding Methods 0.000 claims 1
- 230000001488 breeding effect Effects 0.000 claims 1
- 230000000149 penetrating effect Effects 0.000 claims 1
- 230000000630 rising effect Effects 0.000 claims 1
- 230000008569 process Effects 0.000 description 7
- 230000004048 modification Effects 0.000 description 5
- 238000012986 modification Methods 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 4
- 230000001580 bacterial effect Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000012271 agricultural production Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
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/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
-
- 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/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Business, Economics & Management (AREA)
- Pure & Applied Mathematics (AREA)
- Operations Research (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Animal Husbandry (AREA)
- Probability & Statistics with Applications (AREA)
- General Health & Medical Sciences (AREA)
- Algebra (AREA)
- Agronomy & Crop Science (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A kind of crop production water footprints measuring method based on soil moisture dynamic balance, comprising the following steps: obtain measuring and calculating regional soil type and soil moisture content S；Obtain crop-planting condition data in the measuring and calculating period；Obtain the production characteristic attribute of measuring and calculating object crop；Calculate the soil evaporativeness and crop transpirstion amount in crop growth period；Calculate field flow path surface；Calculate milliosmolarity under deep subsoil；Soil moisture dynamic is characterized, blue, green soil water balance is distinguished；Calculate corresponding supply water and crop indigo plant water resources consumption amount, clear water consumed resource under irrigation method；Crop is calculated accordingly to supply water and crop production indigo plant water footprints, crop production clear water footprint under irrigation method；Calculate crop production water conveyance loss indigo plant water footprints；Crop production water conveyance loss indigo plant water footprints are equal to the product of regional irrigation transmission & distribution water and the region transmission & distribution water loss coefficient, and regional irrigation transmission & distribution water is total irrigation water capacity.The present invention effectively increases the estimation precision of area crops production water footprints.
Description
Technical field
Measuring and calculating field is used and consumed the present invention relates to crop production process water resource, and in particular to one kind is based on soil
The crop production water footprints measuring method of moisture dynamic equilibrium considers irrigation method and water conveyance loss influence factor.
Background technique
Crop production water footprints are that statement crop production process water consumption, water consumption efficiency and the comprehensive and quantitative of water consumption type are commented
Valence index.By distinguishing blue water (irrigation water) and clear water (rainwater), crop production water footprints may be implemented different in agricultural production
The organic unity of type water resource utilization efficiency is evaluated.Accurately, thoroughly evaluating region different time and space scales crop production indigo plant, clear water
Footprint is to establish agricultural water resources regulating the market mechanism, realize the sustainable important base efficiently managed of regional agriculture water resources
This premise.
In existing crop production water footprints computing technique, the fact that soil moisture dynamic balances often is had ignored, only
The total water footprints size of crop growth period can be expressed and constituted；Soil evaporation and two water consumption processes of Crop evapotranspiration cannot be distinguished, still
Different water supply and influence of the irrigation method to crop water consumption intensity are not paid close attention to, can only be calculated under more satisfactory situation in this way and be made material consumption
Water efficiency, calculated result can not accurately reflect the actual consumption of regional agriculture production water resource, limit to water footprints measuring and calculating in agriculture
Applicability in industry water resources assessment, it is difficult to be calculated and regional agriculture water-saving strategy system for different time and space scales crops water footprints
It is fixed that reference is provided.
Summary of the invention
It is an object of the invention to the problems not high for the production water footprints estimation precision of area crops in the prior art, mention
For a kind of crop production water footprints measuring method based on soil moisture dynamic balance, the shadow of different water supply and irrigation method is introduced
The mechanism of sound.
To achieve the goals above, the present invention has the following technical solution:
A kind of crop production water footprints measuring method based on soil moisture dynamic balance, comprising the following steps:
Step 1: obtaining measuring and calculating regional soil type and soil moisture content S；Obtain crop-planting conditional number in the measuring and calculating period
According to the crop-planting condition data includes that the measuring and calculating region day highest temperature, day lowest temperature and intra day ward PR [t] and reference are steamed
Tatol evapotranspiration ET0[t]；Obtain the production characteristic attribute of measuring and calculating object crop, the production characteristic attribute of the measuring and calculating object crop
Crop topsoil moistens rate f under irrigation method, corresponding irrigation method including measuring and calculating objectw, irrigating water quota IRR [t], crop
Growing area water efficiency of irrigation, Crop growing stage gp, crop yield per unit area Y, growth phase vegetation coverage CC*[t]；
Step 2: calculating the soil evaporativeness and crop transpirstion amount in crop growth period；
Crop growth period soil evaporativeness is multiplied by evapotranspiration rate of referential crops, the soil evaporation coefficient with soil coefficient of water pressure
It obtains:
E [t]=Kr[t]×Ke[t]×ET0[t] (1)
In formula (1), E [t] is the t days field soil evaporation capacity, unit mm；Kr[t] is soil moisture pressure coefficient, immeasurable
Guiding principle indicates the degree that maximum evaporation ability is not achieved because of soil moisture deficiency；Ke[t] be the soil evaporation coefficient, dimensionless, by
Vegetation coverage CC*Crop topsoil moistens rate f under [t] and corresponding irrigation methodwCodetermine its size:
Ke[t]=(1-CC*)×fw×Kex (2)
In formula (2), KexFor soil maximum evaporation coefficient, dimensionless refers to evaporation intensity under the complete wet condition of soil；
Crop growth period transpiration quantity is codetermined by vegetation coverage, crop coefficient and soil moisture stress coefficient；
Tr [t]=Ks[t]×CC*[t]×KC,Tr[t]×ET0[t] (3)
In formula (3), Tr [t] is the t days crop evapotranspirations, unit mm；Ks[t] is to influence stomata closing or soil infiltration energy
The water stress factor of power, dimensionless；KC,Tr[t] is Crop transpirstion coefficient, dimensionless；
ET [t]=E [t]+Tr [t] (4)
In formula (4), ET [t] is the t days evaporation in the field transpiration rates, unit mm；
Step 3: field flow path surface is calculated according to rainfall intensity；
Step 4: milliosmolarity under water content measuring and calculating deep subsoil when according to rainfall, irrigation norm and soil saturation；
Step 5: characterization soil moisture dynamic distinguishes blue, green soil water balance；
Step 6: calculating corresponding supply water and crop indigo plant water resources consumption amount CWU under irrigation methodb, clear water consumed resource
CWUg；
Step 7: calculating crop and accordingly supply water and crop production indigo plant water footprints WF under irrigation methodb, crop production clear water foot
Mark WFg；
Step 8: calculating crop production water conveyance loss indigo plant water footprints；Crop production water conveyance loss indigo plant water footprints are filled equal to region
The product of transmission & distribution water Yu the region transmission & distribution water loss coefficient is irrigate, regional irrigation transmission & distribution water is total irrigation water capacity.
The measuring and calculating regional soil type and soil moisture content S of the step 1 are to delimit measuring and calculating pair according to the setting of user
As including crop, region, period, being transferred from system database or the attribute data of field survey and measuring and calculating object.
The expression formula of the step 3 measuring and calculating field flow path surface is as follows:
RO [t]=f (PR [t]) (5)
In formula (5), RO [t] is the t days field flow path surfaces, unit mm；PR [t] is rainfall, unit mm.
The expression formula of milliosmolarity is as follows under the step 4 measuring and calculating deep subsoil:
DP [t]=f (PR [t], IRR [t], Sm) (6)
In formula (6), DP [t] is milliosmolarity under the t days deep subsoils, unit mm；IRR [t] is irrigation norm, and unit is
mm；SmVolumetric(al) moisture content when for soil saturation, unit mm/m.
In crop growth period, the daily soil moisture dynamic in field is characterized the step 5 by formula (7):
S [t]=S [t-1]+PR [t]+IRR [t]-ET [t]-RO [t]-DP [t] (7)
In formula (7): S [t] is the t days soil moisture contents of crop growth period, unit mm；
The fact that balanced according to crop root zone soil moisture dynamic, it is assumed that crop growth period initial soil water is clear water, raw
Educating irrigation and precipitation in the phase is blue water footprints and clear water footprint source respectively, tracks day by day irrigation quantity and precipitation to soil moisture
The contribution proportion of each element is balanced, Crop growing stage indigo plant, green soil moisture dynamic equilibrium are expressed respectively are as follows:
In above formula, Sb[t] is the t days soil indigo plant water contents, unit mm；Sg[t] is the t days soil clear water contents, unit
mm。
The step 6 crop indigo plant water resources consumption amount CWUb, clear water consumed resource CWUgCalculation it is as follows:
In formula (10) (11): t is different water supply and irrigation method, dimensionless；CWUbFor crop indigo plant water resource under t mode
Consumption, unit m3/hm2；CWUgFor crop clear water consumed resource, unit m under t mode3/hm2；Gp is Crop growing stage day
Number；" 10 " are mm and m3/hm2Between Units conversion factor.
The step 7 crop production indigo plant water footprints WFb, crop production clear water footprint WFgCalculation it is as follows:
In formula (12) (13): WFbBlue water footprints, unit m are produced for field crops3/kg；WFgIt is produced for field crops green
Water footprints, unit m3/kg；Y is crop yield per unit area, units/kg/hm2。
The step 8 obtains on the basis of known crop field indigo plant water footprints according to Irrigation Project Design Water application rate:
In formula, WFb_eRefer to and irrigates transmission & distribution water loss indigo plant water footprints, unit m3/kg；η is crop growing spots irrigation water system
Number.
The step 1 carries out the calculating of evapotranspiration rate of referential crops by following formula:
In formula (15): Δ be saturation vapour pressure and temperature curve slope, kPa/ DEG C of unit；RnTo input canopy
Net radiation amount, unit MJ/ (m2·d)；G is soil heat flux, unit MJ/ (m2·d)；γ is psychrometer constant, unit
kPa/℃；
T is 2m eminence mean daily temperature, unit DEG C；μ2For 2m eminence wind speed, unit m/s；esTo be saturated vapour pressure, unit
kPa；eaFor actual water vapor pressure, unit kPa.
Crop topsoil moistens rate f under the corresponding irrigation method of the step 2wRecommendation range it is as follows:
Flood irrigation: 1.0；
Border irrigation: 1.0；
Furrow irrigation, wide bottom: 0.4-0.6；
Furrow irrigation, narrow bottom: 0.6-1.0；
Furrow irrigation is spaced ditch: 0.3-0.5；
Sprinkling irrigation: 1.0；
Trickle irrigation: 0.3-0.4.
Compared to the prior art, the present invention is with following the utility model has the advantages that combining survey based on soil moisture dynamic balance
Calculate regional soil type and soil moisture content S, crop-planting condition and the production feature of measuring and calculating object crop in the measuring and calculating period
Attribute calculates separately soil evaporativeness and crop transpirstion amount in crop growth period, determines crop growth period according to rainfall intensity
Field rainwash, rainfall form rainwash by the processes such as Crop transpirstion, soil evaporation, rest part.Soil is calculated
Milliosmolarity under deep layer, by the processes such as precipitation and irrigation, Crop transpirstion, soil evaporation, rainwash, residual moisture is in soil saturation
Water percolating capacity under the deep soil formed afterwards, the recharge of ground water as corresponding location.Soil moisture dynamic is finally characterized, is distinguished
Blue, green soil water balance calculates crop indigo plant water resource under corresponding water supply and irrigation method according to blue, clear water tatol evapotranspiration and disappears
Consumption and clear water consumed resource supply water and crop indigo plant water resources consumption amount, clear water resource consumption under irrigation method according to corresponding
It is green to acquire crop production indigo plant water footprints, crop production under corresponding water supply and irrigation method respectively for amount and crop yield per unit area
Water footprints finally acquire crop production water conveyance loss indigo plant water footprints.The present invention is based on the crop productions of soil moisture dynamic balance
Water footprints measuring method introduces different water supply and irrigation method to the influencing mechanism of crop water consumption, by different water supply and irrigates
Soil water balance feature under mode can distinguish crop production water footprints under each working condition and form in detail, improve region
The estimation precision of crop production water footprints, when being different, empty scale crop production water footprints are calculated and regional agriculture water-saving strategy
Formulate providing method reference.
Detailed description of the invention
To describe the technical solutions in the embodiments of the present invention more clearly, make required in being described below to embodiment
Attached drawing is briefly described, it should be apparent that, drawings in the following description are some embodiments of the invention, for ability
For the those of ordinary skill of domain, without creative efforts, it can also be obtained according to these attached drawings other attached
Figure.
Crop production water footprints measuring method flow diagram Fig. 1 of the invention.
Specific embodiment
Below with reference to examples and drawings, technical solution of the present invention is clearly and completely further described, is shown
So, described embodiment is a part of the embodiments of the present invention, instead of all the embodiments.
Based on the embodiments of the present invention, those of ordinary skill in the art may be used also in the premise for not making creative work
To carry out several simple modifications and retouching, every other embodiment obtained be shall fall within the protection scope of the present invention.
Referring to Fig. 1, the present invention is based on the crop production water footprints measuring methods of soil moisture dynamic balance, including following step
It is rapid:
Step 1: according to the setting of user delimit measuring and calculating object (crop, region, period), transferred from system database or
The attribute data of field survey and measuring and calculating object, including measuring and calculating regional soil type and soil moisture content S；It is surveyed within the measuring and calculating period
Calculate the region day highest temperature, day lowest temperature and intra day ward PR [t], reference evapotranspiration ET0The crop-plantings conditional number such as [t]
According to；Calculate crop topsoil under the irrigation method of object, corresponding irrigation method and moistens rate fw, irrigating water quota IRR [t], crop
Growing area water efficiency of irrigation, Crop growing stage gp, crop yield per unit area Y, growth phase vegetation coverage CC*[t]
Deng measuring and calculating object crop production characteristic attribute.When evapotranspiration rate of referential crops can not be transferred from database, referred to by following formula
Crop Evapotranspiration calculates:
In formula (15): Δ be saturation vapour pressure and temperature curve slope, kPa/ DEG C of unit；
RnTo input canopy net radiation amount, unit MJ/ (m2·d)；
G is soil heat flux, unit MJ/ (m2·d)；
γ be psychrometer constant, kPa/ DEG C of unit；
T is 2m eminence mean daily temperature, unit DEG C；
μ2For 2m eminence wind speed, unit m/s；
esTo be saturated vapour pressure, unit kPa；
eaFor actual water vapor pressure, unit kPa.
Step 2: calculating soil evaporativeness and crop transpirstion amount in crop growth period.Soil evaporativeness and crop transpirstion amount structure
At evaporation in the field transpiration rate, i.e. water consumption total amount.Crop growth period soil evaporativeness refers to that field does not have the soil of vegetative coverage to steam
Hair, is multiplied to obtain by evapotranspiration rate of referential crops, the soil evaporation coefficient with soil coefficient of water pressure.
E [t]=Kr[t]×Ke[t]×ET0[t] (1)
In formula (1), E [t] is the t days field soil evaporation capacity, unit mm；Kr[t] is soil moisture pressure coefficient, immeasurable
Guiding principle indicates the degree that maximum evaporation ability is not achieved because of soil moisture deficiency；Ke[t] be the soil evaporation coefficient, dimensionless, by
Vegetation coverage CC*Crop topsoil moistens rate f under [t] and corresponding irrigation methodwCodetermine its size:
Ke[t]=(1-CC*)×fw×Kex (2)
In formula (2), KexFor soil maximum evaporation coefficient, dimensionless refers to evaporation intensity under the complete wet condition of soil；
Crop growth period transpiration quantity is codetermined by vegetation coverage, crop coefficient and soil moisture stress coefficient.
Tr [t]=Ks[t]×CC*[t]×KC,Tr[t]×ET0[t] (3)
In formula (3), Tr [t] is the t days crop evapotranspirations, unit mm；Ks[t] is to influence stomata closing or soil infiltration energy
The water stress factor of power, dimensionless；KC,Tr[t] is Crop transpirstion coefficient, dimensionless；
ET [t]=E [t]+Tr [t] (4)
In formula (4), ET [t] is the t days evaporation in the field transpiration rates, unit mm.
Crop topsoil moistens rate f under corresponding irrigation methodwRecommendation range it is as follows:
Flood irrigation: 1.0；
Border irrigation: 1.0；
Furrow irrigation, wide bottom: 0.4-0.6；
Furrow irrigation, narrow bottom: 0.6-1.0；
Furrow irrigation is spaced ditch: 0.3-0.5；
Sprinkling irrigation: 1.0；
Trickle irrigation: 0.3-0.4.
Step 3: measuring and calculating field flow path surface.Crop growth period field rainwash, is determined by rainfall intensity, by making
The rest part of the processes such as object transpiration, soil evaporation, rainfall forms rainwash.
RO [t]=f (PR [t]) (5)
In formula (5), RO [t] is the t days field flow path surfaces, unit mm；PR [t] is rainfall, unit mm.
Step 4: milliosmolarity under measuring and calculating deep subsoil.Soil bacterial diversity amount is steamed by precipitation and irrigation, Crop transpirstion, soil
The processes such as hair, rainwash, water percolating capacity under the deep soil that residual moisture is formed after soil saturation, the ground as corresponding location
It is lauched supply.Water content when soil bacterial diversity amount is by rainfall, irrigation norm and soil saturation codetermines.
DP [t]=f (PR [t], IRR [t], Sm) (6)
In formula (6), DP [t] is milliosmolarity under the t days deep subsoils, unit mm；IRR [t] is irrigation norm, unit mm；Sm
Volumetric(al) moisture content when for soil saturation, unit mm/m.
Step 5: characterization soil moisture dynamic distinguishes blue, green soil water balance.
In crop growth period, the daily soil moisture dynamic in field is characterized by formula (7):
S [t]=S [t-1]+PR [t]+IRR [t]-ET [t]-RO [t]-DP [t] (7)
In formula (7): S [t] is the t days soil moisture contents of crop growth period, unit mm.
The fact that balanced according to crop root zone soil moisture dynamic, it is assumed that crop growth period initial soil water is clear water, raw
Educating irrigation and precipitation in the phase is blue water footprints and clear water footprint source respectively, tracks day by day irrigation quantity and precipitation to soil moisture
The contribution proportion of each element is balanced, Crop growing stage indigo plant, green soil moisture dynamic equilibrium are expressed respectively are as follows:
Above formula: Sb[t] is the t days soil indigo plant water contents, unit mm；Sg[t] is the t days soil clear water contents, unit mm.
Step 6: according to blue, clear water tatol evapotranspiration in step 5, calculating corresponding water supply and irrigation side according to formula (10) (11)
Crop indigo plant water resources consumption amount CWU under formulab, clear water consumed resource CWUg:
In formula (10) (11): i is different water supply and irrigation method, dimensionless；CWUbDisappear for crop indigo plant water resource under t mode
Consumption, unit m3/hm2；CWUgFor crop clear water consumed resource, unit m under t mode3/hm2；Gp is Crop growing stage number of days；
" 10 " are mm and m3/hm2Between Units conversion factor；
Step 7: calculating field crops and produce water footprints.It is supplied water and crop indigo plant water resources consumption under irrigation method according to corresponding
Measure CWUb, clear water consumed resource CWUgWith crop yield per unit area Y, corresponding supply water and crop under irrigation method is acquired respectively
Produce blue water footprints WFb, crop production clear water footprint WFg:
In formula (12) (13): WFbBlue water footprints, unit m are produced for field crops3/kg；WFgClear water is produced for field crops
Footprint, unit m3/kg；Y is crop yield per unit area, units/kg/hm2；
Step 8: calculating crop production water conveyance loss indigo plant water footprints.Crop production water conveyance loss indigo plant water footprints are filled equal to region
Irrigate the product of transmission & distribution water (total irrigation water capacity) Yu the region transmission & distribution water loss coefficient.It can be in known crop field indigo plant water footprints
On the basis of, (including field irrigation water consumption and irrigation water loss amount) is obtained according to Irrigation Project Design Water application rate.
In formula (14), WFb_eRefer to and irrigates transmission & distribution water loss indigo plant water footprints, unit m3/ kg, η are crop growing spots irrigation water
Usage factor, dimensionless.
Step 2,3,4 can quickly be calculated by modeling tool.
Proposed by the present invention is the crop production that irrigation method and water conveyance loss are balanced and considered based on soil moisture dynamic
Water footprints measuring method introduces different water supply and irrigation method to the influencing mechanism of crop water consumption, by different water supply and irrigates
Soil water balance feature under mode can distinguish crop production water footprints under each working condition and form in detail, improve region
The estimation precision of crop production water footprints is calculated and regional agriculture water-saving strategy system for different time and space scales crop production water footprints
Determine providing method reference.
Combining specific features and embodiment above, invention has been described, it is clear that, do not departing from this hair
In the case where bright spirit and scope, it can also be carry out various modifications and is combined.Correspondingly, the specification and drawings are only
The exemplary illustration of the invention that appended claims are defined, and be considered as covered it is any and all in the scope of the invention
Modification, variation, combination or equivalent.Obviously, various changes and modifications can be made to the invention by those skilled in the art, these
Do not depart from the spirit and scope of the present invention modifications and variations also belong to the claims in the present invention and its equivalent technologies range it
It is interior.
Claims (10)
1. a kind of crop production water footprints measuring method based on soil moisture dynamic balance, which is characterized in that including following step
It is rapid:
Step 1: obtaining measuring and calculating regional soil type and soil moisture content S；Obtain crop-planting condition data in the measuring and calculating period, institute
Crop-planting condition data is stated to include the measuring and calculating region day highest temperature, day lowest temperature and intra day ward PR [t] and steam with reference to evaporation
The amount of rising ET0[t]；The production characteristic attribute of measuring and calculating object crop is obtained, the production characteristic attribute of the measuring and calculating object crop includes
Calculate crop topsoil under the irrigation method of object, corresponding irrigation method and moistens rate fw, irrigating water quota IRR [t], crop-planting
Area water efficiency of irrigation, Crop growing stage gp, crop yield per unit area Y, growth phase vegetation coverage CC*[t]；
Step 2: calculating the soil evaporativeness and crop transpirstion amount in crop growth period；
Crop growth period soil evaporativeness is multiplied by evapotranspiration rate of referential crops, the soil evaporation coefficient with soil coefficient of water pressure
It arrives:
E [t]=Kr[t]×Ke[t]×ET0[t] (1)
In formula (1), E [t] is the t days field soil evaporation capacity, unit mm；Kr[t] is soil moisture pressure coefficient, and dimensionless refers to
Show the degree that maximum evaporation ability is not achieved because of soil moisture deficiency；Ke[t] is the soil evaporation coefficient, dimensionless, by vegetation
Coverage CC*Crop topsoil moistens rate f under [t] and corresponding irrigation methodwCodetermine its size:
Ke[t]=(1-CC*)×fw×Kex (2)
In formula (2), KexFor soil maximum evaporation coefficient, dimensionless refers to evaporation intensity under the complete wet condition of soil；
Crop growth period transpiration quantity is codetermined by vegetation coverage, crop coefficient and soil moisture stress coefficient；
Tr [t]=Ks[t]×CC*[t]×KC, Tr[t]×ET0[t] (3)
In formula (3), Tr [t] is the t days crop evapotranspirations, unit mm；Ks[t] is to influence stomata closing or soil infiltration ability
Water stress factor, dimensionless；KC, Tr[t] is Crop transpirstion coefficient, dimensionless；
ET [t]=E [t]+Tr [t] (4)
In formula (4), ET [t] is the t days evaporation in the field transpiration rates, unit mm；
Step 3: field flow path surface is calculated according to rainfall intensity；
Step 4: milliosmolarity under water content measuring and calculating deep subsoil when according to rainfall, irrigation norm and soil saturation；
Step 5: characterization soil moisture dynamic distinguishes blue, green soil water balance；
Step 6: calculating corresponding supply water and crop indigo plant water resources consumption amount CWU under irrigation methodb, clear water consumed resource CWUg；
Step 7: calculating crop and accordingly supply water and crop production indigo plant water footprints WF under irrigation methodb, crop production clear water footprint WFg；
Step 8: calculating crop production water conveyance loss indigo plant water footprints；It is defeated that crop production water conveyance loss indigo plant water footprints are equal to regional irrigation
The product of water distribution quantity and the region transmission & distribution water loss coefficient, regional irrigation transmission & distribution water are total irrigation water capacity.
2. the crop production water footprints measuring method according to claim 1 based on soil moisture dynamic balance, feature exist
In the measuring and calculating regional soil type and soil moisture content S of: the step 1, be measuring and calculating object delimited according to the setting of user, including
Crop, region, period, transferred from system database or field survey and measuring and calculating object attribute data.
3. the crop production water footprints measuring method according to claim 1 based on soil moisture dynamic balance, feature exist
In the expression formula of the step 3 measuring and calculating field flow path surface is as follows:
RO [t]=f (PR [t]) (5)
In formula (5), RO [t] is the t days field flow path surfaces, unit mm；PR [t] is rainfall, unit mm.
4. the crop production water footprints measuring method according to claim 1 based on soil moisture dynamic balance, feature exist
In the expression formula of milliosmolarity is as follows under the step 4 measuring and calculating deep subsoil:
DP [t]=f (PR [t], IRR [t], Sm) (6)
In formula (6), DP [t] is milliosmolarity under the t days deep subsoils, unit mm；IRR [t] is irrigation norm, unit mm；SmFor
Volumetric(al) moisture content when soil saturation, unit mm/m.
5. the crop production water footprints measuring method according to claim 1 based on soil moisture dynamic balance, feature exist
In in crop growth period, the daily soil moisture dynamic in field is characterized the step 5 by formula (7):
S [t]=S [t-1]+PR [t]+IRR [t]-ET [t]-RO [t]-DP [t] (7)
In formula (7): S [t] is the t days soil moisture contents of crop growth period, unit mm；
The fact that balanced according to crop root zone soil moisture dynamic, it is assumed that crop growth period initial soil water is clear water, breeding time
Interior irrigation is blue water footprints and clear water footprint source respectively with precipitation, tracks day by day irrigation quantity and precipitation to soil water balance
The contribution proportion of each element, Crop growing stage indigo plant, green soil moisture dynamic equilibrium are expressed respectively are as follows:
In above formula, Sb[t] is the t days soil indigo plant water contents, unit mm；Sg[t] is the t days soil clear water contents, unit mm.
6. the crop production water footprints measuring method according to claim 1 based on soil moisture dynamic balance, feature exist
In step 6 crop indigo plant water resources consumption amount CWUb, the clear water consumed resource CWUgCalculation it is as follows:
In formula (10) (11): t is different water supply and irrigation method, dimensionless；CWUbFor crop indigo plant water resources consumption under t mode
Amount, unit m3/hm2；CWUgFor crop clear water consumed resource, unit m under t mode3/hm2；Gp is Crop growing stage number of days；
" 10 " are mm and m3/hm2Between Units conversion factor.
7. the crop production water footprints measuring method according to claim 1 based on soil moisture dynamic balance, feature exist
In the step 7 crop production indigo plant water footprints WFb, crop production clear water footprint WFgCalculation it is as follows:
In formula (12) (13): WFbBlue water footprints, unit m are produced for field crops3/kg；WFgClear water foot is produced for field crops
Mark, unit m3/kg；Y is crop yield per unit area, units/kg/hm2。
8. the crop production water footprints measuring method according to claim 1 based on soil moisture dynamic balance, feature exist
In the step 8 obtains on the basis of known crop field indigo plant water footprints according to Irrigation Project Design Water application rate:
In formula, WFb_eRefer to and irrigates transmission & distribution water loss indigo plant water footprints, unit m3/kg；η is crop growing spots water efficiency of irrigation.
9. the crop production water footprints measuring method according to claim 1 based on soil moisture dynamic balance, feature exist
In the step 1 carries out the calculating of evapotranspiration rate of referential crops by following formula:
In formula (15): Δ be saturation vapour pressure and temperature curve slope, kPa/ DEG C of unit；RnTo input the net spoke of canopy
The amount of penetrating, unit MJ/ (m2·d)；G is soil heat flux, unit MJ/ (m2·d)；γ is psychrometer constant, unit kPa/
℃；
T is 2m eminence mean daily temperature, unit DEG C；μ2For 2m eminence wind speed, unit m/s；esTo be saturated vapour pressure, unit kPa；
eaFor actual water vapor pressure, unit kPa.
10. the crop production water footprints measuring method according to claim 1 based on soil moisture dynamic balance, feature exist
In crop topsoil moistens rate f under the corresponding irrigation method of the step 2wRecommendation range it is as follows:
Flood irrigation: 1.0；
Border irrigation: 1.0；
Furrow irrigation, wide bottom: 0.4-0.6；
Furrow irrigation, narrow bottom: 0.6-1.0；
Furrow irrigation is spaced ditch: 0.3-0.5；
Sprinkling irrigation: 1.0；
Trickle irrigation: 0.3-0.4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910784579.8A CN110490473B (en) | 2019-08-23 | 2019-08-23 | Crop production water footprint measuring and calculating method based on soil moisture dynamic balance |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910784579.8A CN110490473B (en) | 2019-08-23 | 2019-08-23 | Crop production water footprint measuring and calculating method based on soil moisture dynamic balance |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110490473A true CN110490473A (en) | 2019-11-22 |
| CN110490473B CN110490473B (en) | 2022-11-01 |
Family
ID=68553312
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910784579.8A Active CN110490473B (en) | 2019-08-23 | 2019-08-23 | Crop production water footprint measuring and calculating method based on soil moisture dynamic balance |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110490473B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112052428A (en) * | 2020-08-31 | 2020-12-08 | 西北农林科技大学 | Method for measuring and calculating production water footprints of live pigs and pork with differentiated production scales |
| CN113408826A (en) * | 2021-07-16 | 2021-09-17 | 西北农林科技大学 | Life cycle-based crop production water footprint measuring and calculating method and system |
| CN113516392A (en) * | 2021-07-16 | 2021-10-19 | 西北农林科技大学 | Virtual water flow resource effect and cooperation degree evaluation method, system and storage medium |
| CN114637353A (en) * | 2022-03-24 | 2022-06-17 | 四川省水利科学研究院 | Agricultural irrigation control method, system and terminal based on multi-environment factor analysis |
| CN115362922A (en) * | 2022-09-30 | 2022-11-22 | 中慧高芯技术(山东)有限公司 | Irrigation system and irrigation method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103226791A (en) * | 2013-04-12 | 2013-07-31 | 西北农林科技大学 | Measuring and calculating method of grain production water footprint of region |
| CN105230451A (en) * | 2015-10-16 | 2016-01-13 | 西北农林科技大学 | Automatic irrigation forecasting device for water shortage of field crops |
| CN107133881A (en) * | 2017-04-11 | 2017-09-05 | 河海大学 | A kind of method for calculating production estimation water footprints using process based on field liquid manure |
| CN108077042A (en) * | 2017-12-04 | 2018-05-29 | 北京农业智能装备技术研究中心 | A kind of winter wheat time of infertility irrigates early warning decision method |
-
2019
- 2019-08-23 CN CN201910784579.8A patent/CN110490473B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103226791A (en) * | 2013-04-12 | 2013-07-31 | 西北农林科技大学 | Measuring and calculating method of grain production water footprint of region |
| CN105230451A (en) * | 2015-10-16 | 2016-01-13 | 西北农林科技大学 | Automatic irrigation forecasting device for water shortage of field crops |
| CN107133881A (en) * | 2017-04-11 | 2017-09-05 | 河海大学 | A kind of method for calculating production estimation water footprints using process based on field liquid manure |
| CN108077042A (en) * | 2017-12-04 | 2018-05-29 | 北京农业智能装备技术研究中心 | A kind of winter wheat time of infertility irrigates early warning decision method |
Non-Patent Citations (2)
| Title |
|---|
| 单纯宇等: "海河流域作物水足迹研究", 《灌溉排水学报》, 31 May 2016 (2016-05-31), pages 50 - 55 * |
| 吴普特等: "区域主要作物生产实体水-虚拟水耦合流动过程", 《科学通报》, 30 June 2019 (2019-06-30), pages 1953 - 1965 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112052428A (en) * | 2020-08-31 | 2020-12-08 | 西北农林科技大学 | Method for measuring and calculating production water footprints of live pigs and pork with differentiated production scales |
| CN112052428B (en) * | 2020-08-31 | 2024-05-07 | 西北农林科技大学 | Measuring and calculating method for water footprint of live pigs and pork in different production scales |
| CN113408826A (en) * | 2021-07-16 | 2021-09-17 | 西北农林科技大学 | Life cycle-based crop production water footprint measuring and calculating method and system |
| CN113516392A (en) * | 2021-07-16 | 2021-10-19 | 西北农林科技大学 | Virtual water flow resource effect and cooperation degree evaluation method, system and storage medium |
| CN113516392B (en) * | 2021-07-16 | 2023-11-21 | 西北农林科技大学 | Virtual water flow resource effect and synergy evaluation method, system and storage medium |
| CN113408826B (en) * | 2021-07-16 | 2024-04-09 | 西北农林科技大学 | Crop production water footprint measuring and calculating method and system based on life cycle |
| CN114637353A (en) * | 2022-03-24 | 2022-06-17 | 四川省水利科学研究院 | Agricultural irrigation control method, system and terminal based on multi-environment factor analysis |
| CN114637353B (en) * | 2022-03-24 | 2022-11-15 | 四川省水利科学研究院 | Agricultural irrigation control method, system and terminal based on multi-environment factor analysis |
| CN115362922A (en) * | 2022-09-30 | 2022-11-22 | 中慧高芯技术(山东)有限公司 | Irrigation system and irrigation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110490473B (en) | 2022-11-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110490473A (en) | A kind of crop production water footprints measuring method based on soil moisture dynamic balance | |
| Fengrui et al. | Water use patterns and agronomic performance for some cropping systems with and without fallow crops in a semi-arid environment of northwest China | |
| Mao et al. | Effects of deficit irrigation on yield and water use of greenhouse grown cucumber in the North China Plain | |
| Adeniran et al. | Water requirements of some selected crops in Kampe dam irrigation project | |
| Vozhehova et al. | Assessment of the CROPWAT 8.0 software reliability for evapotranspiration and crop water requirements calculations | |
| De Wit et al. | Modelling production of field crops and its requirements | |
| Ben-Gal et al. | Temporal robustness of linear relationships between production and transpiration | |
| Zhang et al. | Integrating rainwater harvesting and drip irrigation for water use efficiency improvements in apple orchards of northwest China | |
| Liu et al. | Drip irrigation enhances water use efficiency without losses in cucumber yield and economic benefits in greenhouses in North China | |
| Memon et al. | Crop water requirement and irrigation scheduling of soybean and tomato crop using CROPWAT 8.0 | |
| CN103340132A (en) | Irrigation method for improving water use efficiency of peanuts in dry land | |
| Laphatphakkhanut et al. | IoT-based smart crop-field monitoring of rice cultivation system for irrigation control and its effect on water footprint mitigation | |
| López‐López et al. | Water productivity of rice genotypes with irrigation and drainage | |
| Friedman | Is the crop evapotranspiration rate a good surrogate for the recommended irrigation rate? | |
| CN103918448A (en) | Maize cultivation method for group structure optimization regulation | |
| CN105010057A (en) | Machine-transplanted seedling sposh seedling raising method | |
| CN105075565A (en) | A narrow mulching film side transplanting method for corns | |
| Soussa | Effects of drip irrigation water amount on crop yield, productivity and efficiency of water use in desert regions in Egypt | |
| Zhang et al. | Yield and quality response of greenhouse tomato to different micropores group spacing and irrigation amount. | |
| Mahata et al. | Effect of drip irrigation with mulch cover on yield and water use efficiency in pot culture of brinjal | |
| Balan et al. | Modeling the crop water requirement and irrigation scheduling of banana using CROPWAT 8 model: A case study of Manali watershed, Thrissur, Kerala | |
| Mushab | Forecasting crop coefficient values for cucumber plant (Cucumis sativus) | |
| Hommadi | Evaluation crop evapotranspiration and crop coefficient of squash (zucchini) with subsurface water retention technology (SWRT) and without SWRT | |
| Nasrabad et al. | Soil water distribution and simulation under subsurface drip irrigation in cotton (Gossypium hirsutum) | |
| Koontanakulvong et al. | Improved Irrigation Demand Estimation via Soil Moisture Data from Satellite Images |
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 |