CN106570627A - Crop irrigation water requirement calculation method on future climatic conditions - Google Patents
Crop irrigation water requirement calculation method on future climatic conditions Download PDFInfo
- Publication number
- CN106570627A CN106570627A CN201610951658.XA CN201610951658A CN106570627A CN 106570627 A CN106570627 A CN 106570627A CN 201610951658 A CN201610951658 A CN 201610951658A CN 106570627 A CN106570627 A CN 106570627A
- Authority
- CN
- China
- Prior art keywords
- crop
- data
- history
- future
- future climate
- 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
- 238000004364 calculation method Methods 0.000 title abstract description 7
- 239000003621 irrigation water Substances 0.000 title abstract description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 53
- 230000002262 irrigation Effects 0.000 claims abstract description 48
- 238000003973 irrigation Methods 0.000 claims abstract description 46
- 238000000034 method Methods 0.000 claims abstract description 23
- 239000002689 soil Substances 0.000 claims abstract description 11
- 238000012360 testing method Methods 0.000 claims abstract description 9
- 230000012010 growth Effects 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims description 23
- 238000012937 correction Methods 0.000 claims description 21
- 230000001186 cumulative effect Effects 0.000 claims description 19
- 238000009395 breeding Methods 0.000 claims description 16
- 230000001488 breeding effect Effects 0.000 claims description 16
- 238000000205 computational method Methods 0.000 claims description 14
- 238000009826 distribution Methods 0.000 claims description 14
- 230000005855 radiation Effects 0.000 claims description 7
- 238000001556 precipitation Methods 0.000 claims description 6
- 230000006735 deficit Effects 0.000 claims description 5
- 238000005303 weighing Methods 0.000 claims description 2
- 238000009825 accumulation Methods 0.000 claims 2
- 230000003698 anagen phase Effects 0.000 claims 1
- 238000013439 planning Methods 0.000 abstract description 3
- 238000011156 evaluation Methods 0.000 abstract 1
- 238000007726 management method Methods 0.000 abstract 1
- 238000010792 warming Methods 0.000 abstract 1
- 241000209094 Oryza Species 0.000 description 16
- 235000007164 Oryza sativa Nutrition 0.000 description 16
- 235000009566 rice Nutrition 0.000 description 16
- 238000012913 prioritisation Methods 0.000 description 6
- 230000010354 integration Effects 0.000 description 5
- 208000005156 Dehydration Diseases 0.000 description 4
- 241000196324 Embryophyta Species 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or 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
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
- Y02A40/28—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture specially adapted for farming
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Primary Health Care (AREA)
- Animal Husbandry (AREA)
- General Health & Medical Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Marine Sciences & Fisheries (AREA)
- Agronomy & Crop Science (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention discloses a crop irrigation water requirement calculation method on future climatic conditions. The method comprises the steps of collecting the future climatic mode data, and correcting the future climatic mode data based on the historical actually measured meteorological data to enable the future climatic mode data to be suitable for the climatic change influence evaluation of a region or station scale; according to the growth period data of a crop field test and an accumulated temperature formula, constructing a response model of the crop planting date and the crop growth period length to the temperature; utilizing a penman formula to combine a single crop coefficient method and a soil moisture stress coefficient to calculate the crop daily water requirement; based on a crop irrigation system and a water balance principle to calculate the crop daily irrigation water requirement. Aiming at the influence of the climatic change on the agricultural water resource security, the method of the present invention considers the change of the crop planting date and the growth period caused by the global warming, so that a water resource planning management department can forecast the regional future agricultural water resource utilization amount more accurately, and accordingly, a water resource planning scheme is proposed more reasonably.
Description
Technical field
The present invention relates to crop irrigation requirement computational methods under the conditions of a kind of Future Climate, belong to field of agricultural irrigation.
Background technology
Climate change with temperature rising, precipitation fluctuation as principal character, needs water to generate significantly agricultural irrigation
Impact.Crop irrigation requirement is the important component part of agricultural water, and agricultural water is the important step of hydrologic cycle
One of, temperature rising, Abnormal Precipitation, disaster take place frequently etc., and climate change directly affects hydrologic cycle process, so as to agricultural water
General layout produces significant impact.So estimate crop irrigation requirement under the conditions of Future Climate providing to ensureing agricultural production and saving water
Source is significant.Utilize for China's agricultural safety under the influence of Future Climate Change and water resource simultaneously and theories integration be provided,
Also the more accurate discreet area the future of agriculture water resource utilization of large-scale system optimum department can be facilitated, so as to propose more
Rational water resources scheme.
At present, the computational methods of crop irrigation requirement are mainly with history meteorological data and general circulation model(GCMs)Knot
Cooperation thing Y-factor method Y and principle of water balance combine crop modeling, and crop irrigation requirement is inquired into day by day.But, different fillings
The mode of irrigating correspond to different irrigation programs, have different bounds of pouring water in different breeding times.Additionally, the liter of temperature
Height, not only affects the evapotranspiration of crop, while can also change initial plantation date and the breeding time length of crop.
Therefore, the water balanced calculation based on water requirement, effective rainfall and average seepage, can not completely react paddy rice filling
Irrigate process.How the shadow of climate change to crop initial plantation date, breeding time and irrigation requirement is more accurately estimated
Ring, be the task of top priority of agricultural planting planning and following water resources, be also the key issue that this technology is solved.
The content of the invention
The technical problem to be solved is to provide crop irrigation requirement calculating side under the conditions of a kind of Future Climate
Method, builds the cumulative probability bearing calibration of a kind of history actual measurement meteorological data and GCM climatic model period of history data, will correct
Future Climate mode data afterwards, according to accumulated temperature principle, determines crop starting plantation date and breeding time length.Filled based on crop
System is irrigate with reference to single crop coefficient and principle of water balance, crop irrigation requirement under the conditions of Future Climate, and this is calculated
Method is inquired into the water demand of crop, crop irrigation requirement and breeding time, and the method has reasonability and operability.
The present invention is employed the following technical solutions to solve above-mentioned technical problem:
The present invention provides crop irrigation requirement computational methods under the conditions of a kind of Future Climate, comprises the following steps that:
Step 1, collects the GCM climatic models data and history actual measurement meteorological data of agricultural weather website;
Step 2, based on history the cumulative probability distribution of meteorological data and period of history data in GCM climatic models are surveyed
Cumulative probability is distributed, and to Future Climate mode data offset correction process is carried out;
Step 3, according to crop field test data and accumulated temperature principle, determines that crop initially plants date and cycle breeding time
Required equivalent accumulated temperature;
Step 4, it is flat based on deficit irrigation schedule, Penman formula and the water yield with reference to the Future Climate mode data after correction
Weighing apparatus principle, calculates crop irrigation requirement.
Used as the further prioritization scheme of the present invention, history actual measurement meteorological data includes the drop no less than 30 years in step 1
Water, temperature, radiation, wind speed and vapour pressure day by day data;Period of history data and history actual measurement gas in GCM climatic model data
Image data is corresponded to completely.
It is raw by extracting history actual measurement meteorological data month and sorting in step 2 as the further prioritization scheme of the present invention
Into the cumulative probability distribution of the history actual measurement meteorological data in correspondence month.
As the further prioritization scheme of the present invention, period of history data in GCM climatic models are extracted by month in step 2
And sort, generate the cumulative probability distribution of period of history data in the GCM climatic models in correspondence month.
As the further prioritization scheme of the present invention, offset correction process is carried out to Future Climate mode data in step 2
Concrete grammar be:The probable value of meteorological factor in Future Climate mode data is calculated, history reality is distributed in based on this probable value
Survey in the cumulative probability distribution curve and GCM climatic models of meteorological data and inserted in the cumulative probability distribution curve of period of history data
Value find respective value, using the difference or ratio of the two respective values as addition or multiplicative correction coefficient to Future Climate pattern count
According to being corrected, the Future Climate mode data after correction is generated with this.
Used as the further prioritization scheme of the present invention, step 3 is further included:The temperature in time is observed based on crop test
Data and field test data, it is determined that actual measurement crop initially plants date and breeding time length;The work is determined using accumulated temperature formula
Thing grows required equivalent accumulated temperature;Again following crop is determined based on the temperature data in the Future Climate mode data after correction
Breeding time from date and Crop growing stage length.
Used as the further prioritization scheme of the present invention, step 4 is further included:First, crop is calculated using Penman formula
Evapotranspiration amount in growth and development stage;Then, with reference to single crop coefficient and soil moisture stress coefficient, the water demand of crop is calculated;
Finally, using deficit irrigation schedule, with reference to field water balance principle, crop irrigation requirement is calculated.
The present invention adopts above technical scheme compared with prior art, with following technique effect:The present invention adopts crop
Irrigation program.With reference to single crop coefficient and accumulated temperature computing formula and principle of water balance.Calculate crop irrigation under future condition
Water requirement, utilizes for China's agricultural safety under the influence of Future Climate Change and water resource and provides theories integration.
Description of the drawings
Fig. 1 is following crop irrigation requirement calculation flow chart.
Fig. 2 is based on the offset correction schematic diagram of cumulative probability.
Fig. 3 is the change of rice transplanting date.
Fig. 4 is growth period duration of rice change.
Fig. 5 is the change of Rice irrigation water requirement.
Specific embodiment
Below in conjunction with the accompanying drawings and specific embodiment is described in further detail to technical scheme:
The invention provides crop irrigation requirement computational methods under the conditions of a kind of Future Climate, as shown in figure 1, concrete side
Method flow process is as follows:
(1)Collect the GCM climatic models data and history actual measurement meteorological data of agricultural weather website;Wherein, history actual measurement
Meteorological data includes being no less than precipitation, temperature, radiation, wind speed and the vapour pressure day by day data of 30 years;GCM climatic model data
Middle period of history data and history actual measurement meteorological data are corresponded to completely.
(2)The cumulative probability of the cumulative probability distribution and GCM period of history data of surveying meteorological data based on history is distributed,
Offset correction process is carried out to following GCM data.
First, sequence and sort for many years by extracting history actual measurement meteorological data month, generate the cumulative probability in correspondence month
Distribution function fobs(x)=Pobs(x), wherein, x is probable value, fobsX () is history actual measurement meteorological data correspondence number under the probability
Value;
Then, to the period of history data of GCM climatic models also according to extracting data different months, and the correspondence moon is generated
Part cumulative distribution function fm-o(x)=Pm-o(x), wherein, x is probable value, fm-oX () is climatic model history under the probability
Period meteorological data correspondence numerical value;
Finally, the probable value of meteorological factor in Future Climate mode data is calculated, it is general in history actual measurement based on this probable value
Interpolation search respective value in rate curve and weather mode history period probability curve, by the difference of the two respective values(Or ratio)
Future Climate mode data is corrected as correction coefficient, the Future Climate mode data after correction is generated with this.
xm-p.adjst=xm-p+fobs(x)-fm-o(x) (a)
Wherein, xm-p.adjstFor the meteorological factor after correction, xm-pTo correct front meteorological factor, for temperature, radiation, aqueous vapor
Pressure adopts formula a, for precipitation and wind speed adopt formula b.
(3)According to field test data and accumulated temperature principle, determine that crop initially plants product needed for date and cycle breeding time
Temperature.
Described accumulated temperature computing formula, specially:It was divided into three phases by one day, the first stage rises the moment from the sun
(Hn)To the highest temperature correspondence moment(Hx);Second stage is from the highest temperature correspondence moment to carving at sunset(Ho);Phase III is from sunset
Moment was to second day lowest temperature correspondence moment (Hp).One two-stage temperature is fitted with two sections of the tracks of line voltage, the phase III
It is fitted with square root function.
H0And HnDetermined according to local longitude and latitude, Hx=Ho- 4, Hp=Hn+24.Four moment Hn、Hx、Ho、HpCorresponding temperature
Respectively:Work as Daily minimum temperature(Tn), work as max. daily temperature(Tx), temperature at sunset(To), second day minimum temperature(Tp), its
Middle To=Tx-0.39(Tx-Tp)
The temperature funtion computing formula of each moment t is in one day:
In formula:α=Tx-Tn, R=Tx-To,
To formula(3)Integration, by integration after three results summation, you can obtain the accumulated temperature value of a day.
(4)Crop irrigation requirement is calculated based on Penman formula and principle of water balance, concrete calculation procedure is as follows:
First, Reference Evapotranspiration is calculated using Penman formula:
Wherein, ET0For Reference Evapotranspiration(mmd-1);RnFor net radiation(MJ m-2d-1);G is soil heat flux(MJ
m-2d-1);esFor average saturation vapour pressure(kPa));T is mean temperature(℃);u2For the wind speed of two meters of eminences(ms-1);Δ is
Saturation vapour pressure slope(kPa℃-1);γ is dry and wet constant(kPa℃-1);eaFor actual water vapor pressure(kPa).
Then, water stress factor is combined using single crop coefficient, the water demand of crop is calculated day by day(ETci):
ETci=KcKsET0i(5)
Wherein, ET0iFor Reference Evapotranspiration on the i-thth, mm;KcFor crop coefficient;KsFor water stress factor.
Finally, according to deficit irrigation schedule, with reference to field water balance principle, crop irrigation requirement is determined.
Wherein, water balance formula is:
hi-hi-1=Ri+Ii-Di-Si-ETci(6)
In formula:hi、hi-1Respectively i-th day, i-th -1 water depth, mm;Ri、Ii、Di、SiThe drop of i-th day is represented respectively
Rainfall, irrigation quantity, displacement, leakage, mm.
Technical scheme is further elaborated below by specific embodiment:
By taking the Rice irrigation of Kunshan station as an example, crop irrigation requirement computational methods under the conditions of a kind of Future Climate of the invention,
Specific implementation step is as follows:
(1)Collect 3 periods of future of Kunshan website(2011-2040(2020s), 2041-2070(2050s)With
2071-2100(2080s))Climatic model data and history actual measurement meteorological data(1961-2010);Wherein the present embodiment is not
Carry out climatic model data source and compare the 5th stage of plan in CGCM(CMIP5)Climatic model data, historical data source
China Meteorological data network, the present embodiment adopts four kinds of Climate Scenarios of BCC-CSM1.1 (m) climatic models.
(2)Based on history actual measurement meteorological data cumulative probability distribution and the distribution of GCM periods of history cumulative probability, to future
GCM data carry out offset correction process, as shown in Fig. 2 comprising the following steps that:
1)Sequence and sort for many years by actual measurement meteorological data is extracted month, generate month cumulative probability function fobs(x)=
Pobs(x), wherein, x is probable value, fobsX () is that meteorological data correspondence numerical value is surveyed under the probability;
2)To climatic model period of history data also according to extraction of different months data, and it is tired to generate correspondence month probability
Product function fm-o(x)=Pm-o(x), wherein, x is probable value, fm-oX () is climatic model period of history meteorological data under the probability
Correspondence numerical value;
3)The probable value of meteorological factor in following Meteorological series is calculated, based on this probable value in actual measurement probability curve and weather
Interpolation search respective value in mode history period probability curve, by the difference of the two respective values(Or ratio)As correction coefficient
Following meteorological factor is corrected, the Meteorological series after correction are generated with this.
xm-p.adjst=xm-p+fobs(x)-fm-o(x) (a)
Wherein, xm-p.adjstFor the meteorological factor after correction, xm-pTo correct front meteorological factor, for temperature, radiation, aqueous vapor
Pressure adopts formula a, for precipitation and wind speed adopt formula b.
(3)According to field test data and accumulated temperature principle, determine that crop initially plants date and breeding time length;Wherein,
The accumulated temperature formula that accumulated temperature is adopted is calculated, is comprised the following steps that:It was divided into three phases by one day, the first stage rises the moment from the sun
(Hn)To the highest temperature correspondence moment(Hx);Second stage is from the highest temperature correspondence moment to carving at sunset(Ho);Phase III is from sunset
Moment was to second day lowest temperature correspondence moment (Hp).One and two-stage temperature be to be fitted with two sections of the tracks of line voltage, the 3rd rank
Section is fitted with square root function.H0And HnDetermined according to local longitude and latitude, Hx=Ho- 4, Hp=Hn+24.Four moment Hn、Hx、Ho、
HpCorresponding temperature is respectively:Work as Daily minimum temperature(Tn), work as max. daily temperature(Tx), temperature at sunset(To), second day lowest temperature
Degree(Tp), wherein To=Tx-0.39(Tx-Tp).Each moment temperature funtion computing formula in one day:
In formula:α=Tx-Tn, R=Tx-To,
To above formula integrate, by integration after three results summation, obtain the accumulated temperature value of a day.
On the basis of according to Kunshan station field rice breeding time data in 2011, obtaining the survey region rice transplanting date ought
Amount accumulated temperature is 53000 DEG C, and breeding time equivalent accumulated temperature is 73000 DEG C, and so as to calculate Future Climate Rice under Condition the date is transplanted
With growth period duration of rice length.Result of calculation is as shown in Figures 3 and 4.
(4)Crop irrigation requirement is calculated based on Penman formula and Irrigation Water Requirement for Paddy Rice combined water equilibrium principle, specifically
Calculation procedure is as follows:
First, Reference Evapotranspiration is calculated using Penman formula:
Wherein, ET0For Reference Evapotranspiration(mmd-1);RnFor net radiation(MJ m-2d-1);G is soil heat flux(MJ
m-2d-1);esFor average saturation vapour pressure(kPa));T is mean temperature(℃);u2For the wind speed of two meters of eminences(ms-1);Δ is
Saturation vapour pressure slope(kPa℃-1);γ is dry and wet constant(kPa℃-1);eaFor actual water vapor pressure(kPa).
Then, adopt but crop coefficient method and water stress factor, calculate crop water requirement day by day(ETci):
ETci=KcKsET0i
In formula:ET0iFor Reference Evapotranspiration on the i-thth, mm;KcFor crop coefficient;KsFor water stress factor.KcCoefficient
Using Kunshan mono-season medium rice crop coefficient amendment, breeding time is initial, mid-term and latter stage are taken as respectively 1.05,1.2,1.0.
Different Irrigation pattern correspondence different irrigation systems, the present embodiment Rice irrigation is filled using control irrigation program, control
The difference referred to according to paddy rice different times to moisture-sensitive degree is irrigate, reasonable soil moisture supply is set in each breeding time.Remove
Keep outside 5-25mm water layers in period of seedling establishment, remaining each growing stage does not set up water layer, only holding soil moisture upper control limit is
Saturation moisture content, lower limit takes respectively the 60%-80% of saturation moisture content in different growing,(Wherein, the present invention is in calculating process
Water depth negative value is converted into into water depth for soil moisture content).
Finally, crop irrigation requirement is determined according to field water balance principle, water balance formula is as follows:
hi-hi-1=Ri+Ii-Di-Si-ETci(6)
In formula:hi、hi-1Respectively i-th day, i-th -1 water depth, mm;Ri、Ii、Di、SiThe drop of i-th day is represented respectively
Rainfall, irrigation quantity, displacement, leakage, mm.
Seepage computational methods are as follows:
I. when there is water layer in field,
Si=Ki
In formula:KiFor the per day leakage in rice field in the case of normal water supply, take when there is water layer in field and survey leakage average,
mm。
Ii. when field no water layer, estimate as the following formula:
In formula:SiFor the rice field leakage of i-th day, mm;K0It is mainly relevant with the soil texture for saturation hydraulic conductivity, one
As take 0.1~1.0m/d;β is empirical, generally 50~250, the more glutinous weight of soil, its value is bigger;tiFor soil moisture content
The time that i-th sky and water are experienced at ordinary times, d are reached from saturation state;H be the main root layer depth of paddy rice, m.Calculating Rice irrigation needs
The water yield change such as Fig. 5.
The above, the only specific embodiment in the present invention, but protection scope of the present invention is not limited thereto, and appoints
What be familiar with the people of the technology disclosed herein technical scope in, it will be appreciated that the conversion expected or replacement, all should cover
The present invention include within the scope of, therefore, protection scope of the present invention should be defined by the protection domain of claims.
Claims (7)
1. crop irrigation requirement computational methods under the conditions of a kind of Future Climate, it is characterised in that comprise the following steps that:
Step 1, collects the GCM climatic models data and history actual measurement meteorological data of agricultural weather website;
Step 2, based on history the cumulative probability distribution and the accumulation of period of history data in GCM climatic models of meteorological data are surveyed
Probability distribution, to Future Climate mode data offset correction process is carried out;
Step 3, according to crop field test data and accumulated temperature principle, needed for determining that crop initially plants date and cycle breeding time
Equivalent accumulated temperature;
Step 4, it is former based on deficit irrigation schedule, Penman formula and water balance with reference to the Future Climate mode data after correction
Reason, calculates crop irrigation requirement.
2. crop irrigation requirement computational methods under the conditions of a kind of Future Climate according to claim 1, it is characterised in that
History actual measurement meteorological data includes being no less than precipitation, temperature, radiation, wind speed and the vapour pressure day by day data of 30 years in step 1;
Period of history data and history actual measurement meteorological data are corresponded to completely in GCM climatic model data.
3. crop irrigation requirement computational methods under the conditions of a kind of Future Climate according to claim 1, it is characterised in that
By extracting history actual measurement meteorological data month and sorting in step 2, the history for generating correspondence month surveys the accumulation of meteorological data
Probability distribution.
4. crop irrigation requirement computational methods under the conditions of a kind of Future Climate according to claim 1, it is characterised in that
Extract by month in step 2 and period of history data and sorted in GCM climatic models, in generating the GCM climatic models in correspondence month
The cumulative probability distribution of period of history data.
5. crop irrigation requirement computational methods under the conditions of a kind of Future Climate according to claim 1, it is characterised in that
It is to the concrete grammar that Future Climate mode data carries out offset correction process in step 2:In calculating Future Climate mode data
The probable value of meteorological factor, based on this probable value the cumulative probability distribution curve and GCM gas of history actual measurement meteorological data are distributed in
Interpolation search respective value in the cumulative probability distribution curve of period of history data in time pattern, by the difference of the two respective values or
Ratio is corrected as addition or multiplicative correction coefficient to Future Climate mode data, and with this Future Climate after correction is generated
Mode data.
6. crop irrigation requirement computational methods under the conditions of a kind of Future Climate according to claim 1, it is characterised in that
Step 3 is further included:The temperature record and field test data in time are observed based on crop test, it is determined that actual measurement crop is initial
Plantation date and breeding time length;Equivalent accumulated temperature needed for the crop growth is determined using accumulated temperature formula;It is based on again after correction
Future Climate mode data in temperature data determine following Crop growing stage from date and Crop growing stage length.
7. crop irrigation requirement computational methods under the conditions of a kind of Future Climate according to claim 1, it is characterised in that
Step 4 is further included:First, evapotranspiration amount in the crop growth phase is calculated using Penman formula;Then, with reference to nonoculture thing
Y-factor method Y and soil moisture stress coefficient, calculate the water demand of crop;Finally, it is flat with reference to the field water yield using deficit irrigation schedule
Weighing apparatus principle, calculates crop irrigation requirement.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610951658.XA CN106570627B (en) | 2016-11-02 | 2016-11-02 | Method for calculating crop irrigation water demand under future climate conditions |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610951658.XA CN106570627B (en) | 2016-11-02 | 2016-11-02 | Method for calculating crop irrigation water demand under future climate conditions |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106570627A true CN106570627A (en) | 2017-04-19 |
CN106570627B CN106570627B (en) | 2020-09-29 |
Family
ID=58535243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610951658.XA Active CN106570627B (en) | 2016-11-02 | 2016-11-02 | Method for calculating crop irrigation water demand under future climate conditions |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106570627B (en) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107274040A (en) * | 2017-08-10 | 2017-10-20 | 中国水利水电科学研究院 | Crop level water requirement estimation method based on soil moisture content transformation |
CN107292461A (en) * | 2017-08-14 | 2017-10-24 | 中国水利水电科学研究院 | Regional industry level water requirement estimation method based on Efficiency Statistics |
CN107466816A (en) * | 2017-07-24 | 2017-12-15 | 东华大学 | A kind of irrigation method based on dynamic multilayer extreme learning machine |
CN107950324A (en) * | 2017-12-15 | 2018-04-24 | 上海应用技术大学 | Based on corn irrigation requirement calculates stage by stage irrigation management system and irrigation method |
CN107990936A (en) * | 2017-11-27 | 2018-05-04 | 武汉大学 | A kind of water requirement of rice combined geometry |
CN109258417A (en) * | 2018-12-06 | 2019-01-25 | 中国水利水电科学研究院 | A kind of automatic irrigation method |
CN109258330A (en) * | 2018-09-21 | 2019-01-25 | 安徽省农业科学院土壤肥料研究所 | A kind of cultural method reducing rice field winter green manure production risk |
CN109447426A (en) * | 2018-10-12 | 2019-03-08 | 河海大学 | Response analysis method of the crop structure based on crop water mechanism to changing environment |
CN109657854A (en) * | 2018-12-13 | 2019-04-19 | 中国农业科学院农田灌溉研究所 | A kind of method, system and the equipment of water demand of crop prediction |
CN109858678A (en) * | 2018-12-29 | 2019-06-07 | 航天信息股份有限公司 | A kind of method and system of determining sunflower Meteorological Output |
CN110580657A (en) * | 2019-10-12 | 2019-12-17 | 中国水利水电科学研究院 | agricultural irrigation water demand prediction method |
CN111401649A (en) * | 2020-03-23 | 2020-07-10 | 成都信息工程大学 | Method and system for predicting irrigation flow according to climate change |
CN112163703A (en) * | 2020-09-25 | 2021-01-01 | 中国水利水电科学研究院 | Farmland reference crop evapotranspiration prediction method considering meteorological factor uncertainty |
CN115088595A (en) * | 2022-05-16 | 2022-09-23 | 中联智慧农业股份有限公司 | Automatic irrigation method, device and system for crops and processor |
CN115640956A (en) * | 2022-09-27 | 2023-01-24 | 中国长江三峡集团有限公司 | Future water resource supply and demand balance analysis method |
CN116664340A (en) * | 2023-07-25 | 2023-08-29 | 长江信达软件技术(武汉)有限责任公司 | Irrigation area level water quantity scheduling method suitable for arid regions |
CN117114450A (en) * | 2023-10-19 | 2023-11-24 | 江西农业大学 | Irrigation water decision-making method and system |
CN117435849A (en) * | 2023-12-15 | 2024-01-23 | 华中科技大学 | Agricultural drought monitoring method, system, equipment and storage medium |
CN117455041A (en) * | 2023-10-19 | 2024-01-26 | 上海应用技术大学 | Multi-target management and control utilization method for future water resources obtained through prediction |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120290140A1 (en) * | 2011-05-13 | 2012-11-15 | Groeneveld David P | Method and System to Control Irrigation Across Large Geographic Areas Using Remote Sensing, Weather and Field Level Data |
CN103413035A (en) * | 2013-07-30 | 2013-11-27 | 中国科学院遥感与数字地球研究所 | Model for farmland net irrigation water and irrigation water capacity estimation method |
CN103886392A (en) * | 2014-03-26 | 2014-06-25 | 中国农业科学院农田灌溉研究所 | Winter wheat water consumption predicting method based on weather forecast information |
-
2016
- 2016-11-02 CN CN201610951658.XA patent/CN106570627B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120290140A1 (en) * | 2011-05-13 | 2012-11-15 | Groeneveld David P | Method and System to Control Irrigation Across Large Geographic Areas Using Remote Sensing, Weather and Field Level Data |
CN103413035A (en) * | 2013-07-30 | 2013-11-27 | 中国科学院遥感与数字地球研究所 | Model for farmland net irrigation water and irrigation water capacity estimation method |
CN103886392A (en) * | 2014-03-26 | 2014-06-25 | 中国农业科学院农田灌溉研究所 | Winter wheat water consumption predicting method based on weather forecast information |
Non-Patent Citations (3)
Title |
---|
王卫光 等: "多模式集合模拟未来气候变化对水稻需水量及水分利用效率的影响", 《CNKI:水利学报》 * |
王卫光 等: "气候变化下长江中下游水稻灌溉需水量时空变化特征", 《CNKI:水科学进展》 * |
鲍金丽 等: "控制灌溉条件下水稻灌溉需水量对气候变化的响应", 《CNKI:中国农村水利水电》 * |
Cited By (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107466816A (en) * | 2017-07-24 | 2017-12-15 | 东华大学 | A kind of irrigation method based on dynamic multilayer extreme learning machine |
CN107466816B (en) * | 2017-07-24 | 2020-04-07 | 东华大学 | Irrigation method based on dynamic multilayer extreme learning machine |
CN107274040B (en) * | 2017-08-10 | 2018-09-11 | 中国水利水电科学研究院 | Crop level water requirement estimation method based on soil moisture content transformation |
CN107274040A (en) * | 2017-08-10 | 2017-10-20 | 中国水利水电科学研究院 | Crop level water requirement estimation method based on soil moisture content transformation |
CN107292461A (en) * | 2017-08-14 | 2017-10-24 | 中国水利水电科学研究院 | Regional industry level water requirement estimation method based on Efficiency Statistics |
CN107990936B (en) * | 2017-11-27 | 2019-07-23 | 武汉大学 | A kind of water requirement of rice combined geometry |
CN107990936A (en) * | 2017-11-27 | 2018-05-04 | 武汉大学 | A kind of water requirement of rice combined geometry |
CN107950324A (en) * | 2017-12-15 | 2018-04-24 | 上海应用技术大学 | Based on corn irrigation requirement calculates stage by stage irrigation management system and irrigation method |
CN109258330B (en) * | 2018-09-21 | 2021-02-19 | 安徽省农业科学院土壤肥料研究所 | Cultivation method for reducing production risk of green manure in rice field in winter |
CN109258330A (en) * | 2018-09-21 | 2019-01-25 | 安徽省农业科学院土壤肥料研究所 | A kind of cultural method reducing rice field winter green manure production risk |
CN109447426A (en) * | 2018-10-12 | 2019-03-08 | 河海大学 | Response analysis method of the crop structure based on crop water mechanism to changing environment |
CN109258417B (en) * | 2018-12-06 | 2021-04-06 | 中国水利水电科学研究院 | Automatic irrigation method |
CN109258417A (en) * | 2018-12-06 | 2019-01-25 | 中国水利水电科学研究院 | A kind of automatic irrigation method |
CN109657854B (en) * | 2018-12-13 | 2023-07-04 | 中国农业科学院农田灌溉研究所 | Crop water demand prediction method, system and equipment |
CN109657854A (en) * | 2018-12-13 | 2019-04-19 | 中国农业科学院农田灌溉研究所 | A kind of method, system and the equipment of water demand of crop prediction |
CN109858678A (en) * | 2018-12-29 | 2019-06-07 | 航天信息股份有限公司 | A kind of method and system of determining sunflower Meteorological Output |
CN109858678B (en) * | 2018-12-29 | 2023-04-25 | 航天信息股份有限公司 | Method and system for determining meteorological yield of sunflowers |
WO2021068540A1 (en) * | 2019-10-12 | 2021-04-15 | 中国水利水电科学研究院 | Agricultural irrigation water demand prediction method |
CN110580657A (en) * | 2019-10-12 | 2019-12-17 | 中国水利水电科学研究院 | agricultural irrigation water demand prediction method |
US11209574B2 (en) | 2019-10-12 | 2021-12-28 | China Institute Of Water Resources And Hydropower Research | Method for forecasting an agricultural irrigation water requirement |
JP2022512047A (en) * | 2019-10-12 | 2022-02-02 | チャイナ インスティチュート オブ ウォータ リソースィズ アンド ハイドロパワー リサーチ | Agricultural irrigation water demand forecast method |
JP7149348B2 (en) | 2019-10-12 | 2022-10-06 | チャイナ インスティチュート オブ ウォータ リソースィズ アンド ハイドロパワー リサーチ | Agricultural irrigation water demand forecast method |
CN111401649A (en) * | 2020-03-23 | 2020-07-10 | 成都信息工程大学 | Method and system for predicting irrigation flow according to climate change |
US20220100997A1 (en) * | 2020-09-25 | 2022-03-31 | China Institute Of Water Resources And Hydropower Research | Farmland reference crop evapotranspiration prediction method considering uncertainty of meteorological factors |
CN112163703A (en) * | 2020-09-25 | 2021-01-01 | 中国水利水电科学研究院 | Farmland reference crop evapotranspiration prediction method considering meteorological factor uncertainty |
CN112163703B (en) * | 2020-09-25 | 2024-02-02 | 中国水利水电科学研究院 | Farmland reference crop evapotranspiration prediction method considering weather factor uncertainty |
CN115088595A (en) * | 2022-05-16 | 2022-09-23 | 中联智慧农业股份有限公司 | Automatic irrigation method, device and system for crops and processor |
CN115640956A (en) * | 2022-09-27 | 2023-01-24 | 中国长江三峡集团有限公司 | Future water resource supply and demand balance analysis method |
CN116664340A (en) * | 2023-07-25 | 2023-08-29 | 长江信达软件技术(武汉)有限责任公司 | Irrigation area level water quantity scheduling method suitable for arid regions |
CN116664340B (en) * | 2023-07-25 | 2023-12-12 | 长江信达软件技术(武汉)有限责任公司 | Irrigation area level water quantity scheduling method suitable for arid regions |
CN117114450A (en) * | 2023-10-19 | 2023-11-24 | 江西农业大学 | Irrigation water decision-making method and system |
CN117455041A (en) * | 2023-10-19 | 2024-01-26 | 上海应用技术大学 | Multi-target management and control utilization method for future water resources obtained through prediction |
CN117114450B (en) * | 2023-10-19 | 2024-02-02 | 江西农业大学 | Irrigation water decision-making method and system |
CN117455041B (en) * | 2023-10-19 | 2024-05-28 | 上海应用技术大学 | Multi-target management and control utilization method for future water resources obtained through prediction |
CN117435849A (en) * | 2023-12-15 | 2024-01-23 | 华中科技大学 | Agricultural drought monitoring method, system, equipment and storage medium |
CN117435849B (en) * | 2023-12-15 | 2024-03-29 | 华中科技大学 | Agricultural drought monitoring method, system, equipment and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN106570627B (en) | 2020-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106570627A (en) | Crop irrigation water requirement calculation method on future climatic conditions | |
US11617313B2 (en) | Controlling agricultural production areas | |
CN103645295B (en) | Multilayer soil moisture simulation method and multilayer soil moisture simulation system | |
CN107103040B (en) | Irrigated area basic data acquisition system | |
CN107945042A (en) | A kind of plant growth irrigation decision control system | |
Minda et al. | The combined effect of elevation and meteorology on potato crop dynamics: a 10-year study in the Gamo Highlands, Ethiopia | |
CN115376016B (en) | Actual rice field irrigation area identification method based on combination of vegetation water index and evapotranspiration | |
Yang et al. | Estimation of groundwater use by crop production simulated by DSSAT‐wheat and DSSAT‐maize models in the piedmont region of the North China Plain | |
CN111280019A (en) | Soil moisture digital prediction and irrigation early warning method | |
EP3986107A1 (en) | System for optimizing the use of water in irrigation based on the predictive calculation of the soil's water potential | |
CN106934534A (en) | Rattan knot melon Irrigation Project Design water effective utilization coefficients computational methods long | |
CN112837169A (en) | Gridding drought comprehensive monitoring, early warning and evaluating method in drought process | |
CN110599360A (en) | High-resolution remote sensing estimation method for evapotranspiration of crops in arid region | |
CN109829134A (en) | Long sequence field irrigation quota measuring method | |
Schneider et al. | Potential and versatility of WEAP model (Water Evaluation and Planning System) for hydrological assessments of AWD (Alternate Wetting and Drying) in irrigated rice | |
CN108617355A (en) | A kind of cluster greenhouse irrigation decision-making technique and system | |
AU2023203307A1 (en) | Controlling Agricultural Production Areas | |
Adhikari et al. | Calibration and validation of CSM‐CROPGRO‐Cotton model using lysimeter data in the Texas High Plains | |
Igbadun | Estimation of crop water use of rain-fed maize and groundnut using mini-lysimeters | |
CN101647336A (en) | Topsoil salt-control method of planting salt-tolerant plant in mudflat kaline soil | |
Western et al. | An evaluation of a methodology for seasonal soil water forecasting for Australian dry land cropping systems | |
CN112715322B (en) | Method and device for obtaining agricultural irrigation water | |
Jovanovic et al. | Total evaporation estimates from a Renosterveld and dryland wheat/fallow surface at the Voëlvlei Nature Reserve (South Africa) | |
Abdul Karim et al. | FAO 56 model and remote sensing for the estimation of crop-water requirement in main branch canal of the Bhadra Command area, Karnataka State | |
CN111192315A (en) | Actual irrigation area extraction method based on multi-source information |
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 |