CN111461909A - Short-term prediction method for farmland evapotranspiration - Google Patents
Short-term prediction method for farmland evapotranspiration Download PDFInfo
- Publication number
- CN111461909A CN111461909A CN202010256456.XA CN202010256456A CN111461909A CN 111461909 A CN111461909 A CN 111461909A CN 202010256456 A CN202010256456 A CN 202010256456A CN 111461909 A CN111461909 A CN 111461909A
- Authority
- CN
- China
- Prior art keywords
- evapotranspiration
- short
- day
- prediction
- crop
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 66
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 29
- 238000012937 correction Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 12
- 239000002689 soil Substances 0.000 claims description 9
- 230000004907 flux Effects 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 4
- 238000000611 regression analysis Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 241000209140 Triticum Species 0.000 description 4
- 235000021307 Triticum Nutrition 0.000 description 4
- 240000008042 Zea mays Species 0.000 description 4
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 4
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 4
- 235000005822 corn Nutrition 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000002262 irrigation Effects 0.000 description 3
- 238000003973 irrigation Methods 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 239000003621 irrigation water Substances 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000005068 transpiration Effects 0.000 description 2
- 101001116926 Homo sapiens Protocadherin alpha-7 Proteins 0.000 description 1
- 102100024275 Protocadherin alpha-7 Human genes 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 239000003337 fertilizer Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000013468 resource allocation Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
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
- 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
-
- 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Animal Husbandry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Mining & Mineral Resources (AREA)
- General Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a short-term prediction method for farmland evapotranspiration, which relates to the technical field of crop monitoring and comprises the steps of determining a prediction reference day and a short-term prediction day; obtaining reference crop evapotranspiration ET of predicted reference day by PM method0(ii) a Obtaining actual crop evapotranspiration ET of target farmland prediction reference day by using vorticity correlation systemc(ii) a Calculating crop coefficient K for short-term prediction of target farmlandc(ii) a Prediction of reference crop Water demand ET 'for each short-term forecast day using the Hargreaves-Samani model'0(ii) a Obtaining ET 'of evapotranspiration amount of each short-term prediction day'cCalculating historical experience correction parameters α, and calculating predicted evapotranspiration ET of each short-term predicted day of the target farmlandcAnd completing short-term prediction of the evapotranspiration of the farmland. The method has higher application value for short-term prediction of the evapotranspiration of the farmland, high prediction precision and good stability, and provides scientific basis for farmland water management in future plain areas.
Description
Technical Field
The invention relates to the technical field of crop monitoring, in particular to a farmland evapotranspiration short-term prediction method considering crop coefficient dynamic change and historical experience.
Background
The Evapotranspiration process is an important component of land hydrologic cycle, the Evapotranspiration volume (ET)c) Has important guiding significance for irrigation planning and regional water resource allocation. The accurate prediction of the evapotranspiration amount can save the irrigation water consumption to a certain extent, and particularly in arid and semiarid regions with larger agricultural water consumption and higher dependence of crop growth on irrigation. Therefore, in order to better manage the irrigation water consumption of crops and improve the water utilization efficiency of crops, the crop ET is urgently needed to be treatedcAnd (6) carrying out accurate prediction.
The prediction of evapotranspiration is generally carried out by methods classified into 4 types: time series methods, grey model methods, intelligent algorithms, and empirical formula methods. The time-series method has uncertainty in prediction accuracy because the data used is single (only historical data of the amount of evapotranspiration is needed), and the latest data contains extremely important prediction information. The gray prediction method is essentially an exponential model, and when the target function increases by zero, the systematic error is serious, and the error is serious when the prediction period is more. The intelligent algorithm needs a dynamic data training system, needs a large amount of data in the early stage and is not suitable for short-term prediction. Compared with other more complex mathematical models, the single crop coefficient method only needs less meteorological data and basic parameters of crops and soil, can carry out parameter real-time optimization according to the types of the crops, can reflect the actual situation better and has higher calculation precision.
ET based on single crop coefficient methodcWhen predicting, the method mainly depends on the amount of the reference crop transpiration (ET)0) And the crop coefficient KcAnd (4) predicting the precision. With the improvement of the precision of weather forecast, public weather forecast data is gradually applied to ET as an important data source0Forecasting, commonly used methods are the Penman-Monteith (PM) method, Fourier analysis model andHargreaves-Samani (H-S) model, and the like. Wherein, although the PM method has higher calculation precision, the PM method needs more meteorological data; the Fourier analysis model has high forecasting precision and can be conveniently embedded into an irrigation system, but the model has great instability; the H-S model is widely applied due to the advantages of low complexity, few input factors (the highest temperature and the lowest temperature), high prediction precision and the like. Furthermore, K varies with the crop itself and the external conditionscAlso constantly changing, but the overall change law is stable during the crop growth period, KcK for crops of the same species, except in relation to the particular crop speciescAlso has obvious regional and time sequence differences which comprehensively reflect the types of crops, soil water and fertilizer conditions and field management level to ETcThe influence of (c). At present, ETcPrediction method KcMainly recommended by FAO-56, obtained by constant or linear difference between two constants, although the recommended value of FAO-56 can be adjusted according to local environmental climate conditions, the method can be used for calculating the terrestrial evapotranspiration process including vegetation transpiration and soil evaporation, is suitable for calculating the process with large time step, but cannot reflect the daily dynamic change condition of crops, and ET is carried out when ET is carried outcDuring short-term prediction, a predicted value has a slightly larger deviation from an actual value. Therefore, for KcThe acquisition of (2) needs to take into account the influence of the growth stage of the crop on its dynamic changes.
In summary, the short-term prediction of field evapotranspiration in the prior art has the following drawbacks:
(1) FAO-56 recommended fixed or simple difference KcIs suitable for long period ET0For short term ET0The prediction application value of the method is low, and the basic crop coefficient curve in the crop coefficient method is only formed by the determined 3 nodes and then linear difference values, so that the crop growth process is simplified, and large deviation can be caused;
(2)ETcthe historical experience rule of the method is not considered enough, and extreme deviation is easy to cause if the model estimation is adopted, so that the stability is poor.
Disclosure of Invention
The invention aims to provide a farmland evapotranspiration short-term prediction method considering crop coefficient dynamic change and historical experience so as to alleviate the problems.
In order to alleviate the problems, the technical scheme adopted by the invention is as follows:
a short-term prediction method for farmland evapotranspiration comprises the following steps:
s1, determining a prediction reference day and a short-term prediction day, wherein the prediction reference day is the current day for performing prediction operation, and the short-term prediction day is 1-15 days from the next day of the prediction reference day;
s2, obtaining reference crop evapotranspiration ET of the predicted reference day by utilizing a PM method based on meteorological data of the target farmland0;
S3, obtaining actual crop evapotranspiration ET of the target farmland prediction reference day by using a vorticity correlation methodc;
S4, predicting reference crop evapotranspiration ET of reference day based on single crop coefficient method0And the actual measured evapotranspiration ET of the cropscAnd calculating to obtain a crop coefficient K used for short-term prediction of the target farmlandc;
S5, predicting reference crop water demand ET 'of each short-term prediction day by adopting Hargreaves-Samani model based on historical weather forecast data of target farmland'0;
S6 according to crop coefficient KcAnd ET 'of crop Water demand for each short-term forecast day'0Obtaining ET 'of evapotranspiration amount of each short-term predicted day'c;
S7, according to the relation function f (L AI) of the target farmland and the mean value of different growth stages of the actually measured evapotranspiration amount of the vorticity correlation systemCalculating a historical experience correction parameter α, and fitting a relation function f (L AI) according to the historical crops L AI of the target farmland and the evapotranspiration;
s8, correcting parameter α and short-term prediction day evapotranspiration amount ET 'according to historical experience'cCalculating the predicted evapotranspiration ET' of each short-term predicted day of the target farmlandcAnd completing short-term prediction of the evapotranspiration of the farmland.
Further, in the step S2, the crop evapotranspiration amount ET is referred to0The calculation formula of (a) is as follows:
wherein, Delta is the slope of the saturated water-vapor pressure curve, Rn is the net surface radiation, G is the soil heat flux, gamma is the dry-wet constant, u is the2At a height of 2 m, esSaturated water vapor pressure, eaThe actual water vapor pressure.
Further, in the step S5, reference crop water demand ET'0The calculation formula of (a) is as follows:
ET′0=0.408K(Tmax-Tmin)n(Tmean+Toff)Ra
wherein K is a conversion coefficient; t ismax、TminRespectively the highest daily temperature and the lowest daily temperature; n is an exponential coefficient; t ismeanThe daily average air temperature; t isoffIs the temperature offset; raRadiating the top layer of the atmosphere.
Further, in the step S6, ET'c=ET′0×Kc。
Furthermore, in step S7, the obtaining method of the relationship function f (L AI) includes taking the historical evapotranspiration of the target farmland as a dependent variable and the historical crop L AI as an independent variable, and performing regression analysis on different growth stages of the crop to obtain the relationship function f (L AI).
Further, in the step S8, the evapotranspiration ET ″c=α×Kc×ET′c。
Compared with the prior art, the invention has the beneficial effects that: on the basis of considering the dynamic change of the crop coefficient and the historical empirical rule of the crop growth, a crop coefficient dynamic prediction method is constructed based on a step-by-step deduction algorithm, and ET is corrected by combining the historical empirical threshold of the crop growthcThe prediction model can obtain the crop evapotranspiration which is more accordant with the actual growth condition of crops in the future by improvement, has higher application value for short-term prediction of the farmland evapotranspiration, high prediction precision and good stability, and provides scientific basis for farmland water management in the future plain area.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained according to the drawings without inventive efforts.
FIG. 1 is a flow chart of a short term method of predicting field evapotranspiration in accordance with the invention;
FIG. 2 is a comparison of single crop coefficient method models with predicted results considering dynamic changes in crop coefficients during the growth period of winter wheat and summer corn;
FIG. 3 is a graph comparing predicted values for different forecast evapotranspiration values based on consideration of dynamic changes in crop coefficients.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Examples
Referring to fig. 1, the present embodiment provides a short-term farmland evapotranspiration prediction method considering dynamic changes of crop coefficients and historical experiences, which takes winter wheat and summer corn of a certain target farmland as research objects, and specifically includes the following steps:
and S1, determining a prediction reference day and a short-term prediction day, wherein the prediction reference day is the current day for performing the prediction operation, and the short-term prediction day is 1-15 days in the future from the next day of the prediction reference day.
In this embodiment, the short-term prediction day is set according to the time specification of the weather forecast, and in the weather forecast, the short-term prediction day is set within 15 days.
S2 meteorological data (pressure P) based on target farmlandaWind speed U and maximum air temperature TmaxMinimum air temperature TminAverage relative humidity HRSunshine duration n, rainfall P, etc.), and the reference crop evapotranspiration amount ET of the predicted reference day is obtained by the PM method0The calculation formula is as follows:
wherein, Delta is the slope of the saturated water-vapor pressure curve, Rn is the net surface radiation, G is the soil heat flux, gamma is the dry-wet constant, u is the2At a height of 2 m, esSaturated water vapor pressure, eaThe actual water vapor pressure.
In the present embodiment, the meteorological data is measured data for predicting a reference day for a target farmland.
S3, utilizing vorticity phaseClosing the method, obtaining the actual crop evapotranspiration ET of the target farmland prediction reference dayc。
In this example, the measured evapotranspiration ET of the cropcMeasured using the vorticity correlation system (Campbell scientific inc., USA) and calculated by the formula:
in the formula: w' is the vertical wind speed pulsation quantity, m/s; q' is the water vapor density pulsation value, g/m3Vorticity correlation system was manufactured by CSAT3 model three-dimensional ultrasonic anemometer (Campbell Scientific inc., USA), L I7500 CO2/H2O open-circuit gas analyzer (L I-COR, USA), HMP45C air temperature and humidity sensor (Campbell Scientific Inc., USA) and CR5000 type data collector (Campbell Scientific Inc., USA). The system is installed in the central and south of the test area according to the climatic features of the test area for northeast and northwest wind, ensuring large wind wave zone length.the installation height of each probe is 3.1m above the ground surface, and the net radiation R is 3.1mnNet radiation sensor (Kipp) by CNR4&Zonen, Netherlands) at an installation height of 4.0m above the surface, soil heat flux G was measured from two HFP01 soil heat flux panels (HuksefluxUSA, inc., USA) located 2cm below the surface, and the average of all measurements taken over 30min was taken as the recorded value for each time, day ETcAnd accumulating the 24h data, eliminating abnormal data according to the following principles that ① precipitation time periods and data of 1h before and after precipitation time periods, ② data obviously exceeding physical meanings, ③ data with abnormal sensor states and eliminating errors caused by energy non-closure by calculating a daily wave-text ratio and correcting latent heat flux in the process of processing the actual data related to the vorticity.
S4, predicting reference crop evapotranspiration ET of reference day based on single crop coefficient method0And the actual measured evapotranspiration ET of the cropscAnd calculating to obtain a crop coefficient K used for short-term prediction of the target farmlandcThe calculation formula is as follows:
in this embodiment, the single crop coefficient method (directly recommended by FAO-56) is a method that requires less meteorological data and basic parameters of crops and soil, and can perform real-time parameter optimization according to the types of crops, compared to other more complex mathematical models, and this method can reflect the actual situation better and has higher calculation accuracy.
S5, predicting reference crop water demand ET 'of each short-term prediction day by adopting Hargreaves-Samani model based on historical weather forecast data of target farmland'0The calculation formula is as follows:
ET′0=0.408K(Tmax-Tmin)n(Tmean+Toff)Ra
wherein, ET'0The unit is mm/d; k is a conversion coefficient; t ismax、TminRespectively the highest daily temperature and the lowest daily temperature; n is an exponential coefficient; t ismeanThe daily average air temperature; t isoffIs the temperature offset; raRadiating the top layer of the atmosphere.
In this embodiment, the maximum temperature of the meteorological station 2000-year old 2017 belonging to the target farmland day by day, the minimum temperature and the atmospheric top radiation are used as independent variables, the day-by-day reference crop water demand calculated by the PM formula is used as a dependent variable, and the HS model original parameter K is 0.0023, n is 0.5, T is 0.0023offCarrying out nonlinear regression analysis on the HS model by taking 17.8 as an initial value, carrying out parameter fitting optimization solution by using EXCE L software programming solution and search method alternate iteration, and obtaining model parameters of a target farmland through optimization, wherein the model parameters are respectively K0.001138, n 0.4925 and T0.001138off43.33, using the PM model calculation result as a control, the MAE and RMSE of the optimized HS model in the validation period (2016-.
S6 according to crop coefficient KcAnd ET 'of crop Water demand for each short-term forecast day'0Obtaining ET 'of evapotranspiration amount of each short-term predicted day'cThe calculation formula is as follows:
ET′c=ET′0×Kc
i.e. the crop coefficient KcMultiplying by the crop Water demand ET 'of the first short-term forecast day'0Obtaining an amount of evapotranspiration ET 'of the first short-term predicted day'cThe coefficient of crop KcMultiplying by the crop Water demand ET 'of the second short term forecast day'0Obtaining a second short-term predicted day evapotranspiration amount ET'c… … are provided. The crop coefficient used on each short-term prediction day is the same and is the crop coefficient Kc。
S7, according to the relation function f (L AI) of the target farmland and the mean value of different growth stages of the actually measured evapotranspiration amount of the vorticity correlation systemCalculating the historical empirical correction parameter α according to the following formula:
the relation function f (L AI) is obtained according to the historical crop L AI (leaf area index) and the evapotranspiration of the target farmland, regression analysis is carried out on different growth stages of crops by taking the historical evapotranspiration of the target farmland as a dependent variable and the historical crop L AI as an independent variable, and the relation function f (L AI) is obtained.
S8, correcting parameter α and short-term prediction day evapotranspiration amount ET 'according to historical experience'cCalculating the predicted evapotranspiration ET' of each short-term predicted day of the target farmlandcThe calculation formula is as follows:
ET″c=α×Kc×ET′0。
FIG. 2 is a comparison graph of the prediction results of the single crop coefficient method model and the prediction results considering the dynamic changes of the crop coefficients in the growing period of winter wheat and summer corn, and it can be seen from the graph that when the evapotranspiration amount is small, the goodness of fit of the prediction values of the single crop coefficient method and the prediction method considering the dynamic changes of the crop coefficients and the measured vorticity value is high, and with the increase of the evapotranspiration value, the goodness of fit is gradually reduced, but the single crop coefficient method has a larger deviation from the measured EC value than the prediction method considering the dynamic changes of the crop coefficients, and it can be found that the prediction method considering the dynamic changes of the crop coefficients has a higher accuracy and smaller errors, and when the short-term prediction of the evapotranspiration amount is performed, the prediction method considering the dynamic changes of the crop coefficients is superior to.
Fig. 3 is a comparison graph of predicted evapotranspiration values for different prediction periods based on consideration of dynamic changes of crop coefficients, and it can be seen from the graph that errors between the evapotranspiration values and actual measurement values for different prediction periods (1d, 3d, 6d, 9d, 14d) of winter wheat and summer corn predicted by a dynamic crop coefficient prediction model are small, and the prediction accuracy for the prediction periods 1-7 d is more than 87%.
Table 1 shows that the prediction accuracy comparison of the pre-and post-model evapotranspiration is improved, the corrected model can accurately track the dynamic change of the evapotranspiration, compared with the model before correction, the prediction error of the mean value of the model after α correction is reduced by 1.60%, the RMSE and the MAE are respectively reduced by 0.221mm/d, 0.134mm/d and R2And ACC are raised by 0.221 and 10.24% respectively compared to the model before correction. The correction model takes account of the historical empirical laws and regulations of crop growth, reduces the evapotranspiration prediction error possibly caused by the extreme state of crops in the crop growth process, and further improves the prediction precision and the applicability of the model.
TABLE 1
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (8)
1. A short-term prediction method for farmland evapotranspiration is characterized by comprising the following steps:
s1, determining a prediction reference day and a short-term prediction day, wherein the prediction reference day is the current day for performing prediction operation, and the short-term prediction day is 1-15 days from the next day of the prediction reference day;
s2, obtaining reference crop evapotranspiration ET of the predicted reference day by utilizing a PM method based on meteorological data of the target farmland0;
S3, obtaining actual crop evapotranspiration ET of the target farmland prediction reference day by using a vorticity correlation methodc;
S4, predicting reference crop evapotranspiration ET of reference day based on single crop coefficient method0And the actual measured evapotranspiration ET of the cropscAnd calculating to obtain a crop coefficient K used for short-term prediction of the target farmlandc;
S5, predicting reference crop water demand ET 'of each short-term prediction day by adopting Hargreaves-Samani model based on historical weather forecast data of target farmland'0;
S6 according to crop coefficient KcAnd ET 'of crop Water demand for each short-term forecast day'0Obtaining ET 'of evapotranspiration amount of each short-term predicted day'c;
S7, according to the relation function f (L AI) of the target farmland and the mean value of different growth stages of the actually measured evapotranspiration amount of the vorticity correlation systemCalculating a historical experience correction parameter α, and fitting a relation function f (L AI) according to the historical crops L AI of the target farmland and the evapotranspiration;
s8, correcting parameter α and short-term prediction day evapotranspiration amount ET 'according to historical experience'cCalculating the predicted evapotranspiration ET' of each short-term predicted day of the target farmlandcAnd completing short-term prediction of the evapotranspiration of the farmland.
2. The method for short-term prediction of agricultural evapotranspiration as claimed in claim 1, wherein in step S2, the crop evapotranspiration ET is referred to0The calculation formula of (a) is as follows:
whereinDelta is the slope of the saturated water-vapor pressure curve, Rn is the net surface radiation, G is the soil heat flux, gamma is the dry-wet constant, u2At a height of 2 m, esSaturated water vapor pressure, eaThe actual water vapor pressure.
4. The method for short-term prediction of agricultural evapotranspiration according to claim 3, wherein in step S5, reference is made to the crop water demand ET'0The calculation formula of (a) is as follows:
ET′0=0.408K(Tmax-Tmin)n(Tmean+Toff)Ra
wherein K is a conversion coefficient; t ismax、TminRespectively the highest daily temperature and the lowest daily temperature; n is an exponential coefficient; t ismeanThe daily average air temperature; t isoffIs the temperature offset; raRadiating the top layer of the atmosphere.
5. The method for short-term prediction of farmland evapotranspiration as claimed in claim 4, wherein in step S6, ET'c=ET′0×Kc。
6. The method for short-term prediction of farmland evapotranspiration as claimed in claim 5, wherein in step S7, the obtaining method of the relation function f (L AI) comprises the step of carrying out regression analysis on different growth stages of crops by taking the historical evapotranspiration of the target farmland as a dependent variable and the historical crop L AI as an independent variable to obtain the relation function f (L AI).
8. The method for short-term prediction of farmland evapotranspiration as claimed in claim 7, wherein in step S8, the evapotranspiration ET ″c=α×Kc×ET′c。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010256456.XA CN111461909B (en) | 2020-04-02 | 2020-04-02 | Short-term prediction method for farmland evapotranspiration |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010256456.XA CN111461909B (en) | 2020-04-02 | 2020-04-02 | Short-term prediction method for farmland evapotranspiration |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111461909A true CN111461909A (en) | 2020-07-28 |
CN111461909B CN111461909B (en) | 2023-02-28 |
Family
ID=71684383
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010256456.XA Active CN111461909B (en) | 2020-04-02 | 2020-04-02 | Short-term prediction method for farmland evapotranspiration |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111461909B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112035089A (en) * | 2020-08-20 | 2020-12-04 | 中国林业科学研究院林业研究所 | Data processing method for vorticity-related observation flux data |
CN112136667A (en) * | 2020-11-26 | 2020-12-29 | 江苏久智环境科技服务有限公司 | Intelligent sprinkling irrigation method and system based on edge machine learning |
CN112507558A (en) * | 2020-12-15 | 2021-03-16 | 中国水利水电科学研究院 | Reference crop evapotranspiration amount prediction method based on main control factor correction |
CN114021486A (en) * | 2022-01-05 | 2022-02-08 | 中化现代农业有限公司 | Crop evapotranspiration prediction method and device |
CN114943361A (en) * | 2022-03-15 | 2022-08-26 | 水利部交通运输部国家能源局南京水利科学研究院 | Method for estimating evapotranspiration of reference crops in data-lacking areas |
CN115049126A (en) * | 2022-06-10 | 2022-09-13 | 中国水利水电科学研究院 | Evapotranspiration prediction method based on temperature effect and historical threshold |
CN115125903A (en) * | 2022-07-08 | 2022-09-30 | 中水三立数据技术股份有限公司 | Automatic irrigation and drainage management method for field integrated gate based on water demand prediction |
CN117678506A (en) * | 2024-02-02 | 2024-03-12 | 水利部牧区水利科学研究所 | Visual irrigation decision-making auxiliary system based on weather forecast and control method |
CN118366059A (en) * | 2024-06-20 | 2024-07-19 | 山东锋士信息技术有限公司 | Crop water demand calculating method based on optical and SAR data fusion |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104521699A (en) * | 2014-11-18 | 2015-04-22 | 华北水利水电大学 | Field intelligent irrigation on-line control management method |
CN105260940A (en) * | 2015-10-22 | 2016-01-20 | 南京信息工程大学 | Crop coefficient correction method based on farmland evapotranspiration observation |
CN106355264A (en) * | 2016-08-11 | 2017-01-25 | 河海大学 | Combined prediction method of reference crop evapotranspiration |
CN107950324A (en) * | 2017-12-15 | 2018-04-24 | 上海应用技术大学 | Based on corn irrigation requirement calculates stage by stage irrigation management system and irrigation method |
CN108491963A (en) * | 2018-03-13 | 2018-09-04 | 河海大学 | A kind of evapotranspiration computational methods considering vegetation growth season succession |
WO2018173045A1 (en) * | 2017-03-20 | 2018-09-27 | Supplant Ltd. | Systems and methods for planning crop irrigation |
CN110501761A (en) * | 2019-08-23 | 2019-11-26 | 中国水利水电科学研究院 | A kind of difference leading time area crops ETc prediction methods |
-
2020
- 2020-04-02 CN CN202010256456.XA patent/CN111461909B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104521699A (en) * | 2014-11-18 | 2015-04-22 | 华北水利水电大学 | Field intelligent irrigation on-line control management method |
CN105260940A (en) * | 2015-10-22 | 2016-01-20 | 南京信息工程大学 | Crop coefficient correction method based on farmland evapotranspiration observation |
CN106355264A (en) * | 2016-08-11 | 2017-01-25 | 河海大学 | Combined prediction method of reference crop evapotranspiration |
WO2018173045A1 (en) * | 2017-03-20 | 2018-09-27 | Supplant Ltd. | Systems and methods for planning crop irrigation |
CN107950324A (en) * | 2017-12-15 | 2018-04-24 | 上海应用技术大学 | Based on corn irrigation requirement calculates stage by stage irrigation management system and irrigation method |
CN108491963A (en) * | 2018-03-13 | 2018-09-04 | 河海大学 | A kind of evapotranspiration computational methods considering vegetation growth season succession |
CN110501761A (en) * | 2019-08-23 | 2019-11-26 | 中国水利水电科学研究院 | A kind of difference leading time area crops ETc prediction methods |
Non-Patent Citations (1)
Title |
---|
段晨斐等: "基于无人机热红外遥感的夏玉米蒸散量估算及其影响因子", 《节水灌溉》, no. 12, 5 December 2019 (2019-12-05), pages 110 - 116 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112035089A (en) * | 2020-08-20 | 2020-12-04 | 中国林业科学研究院林业研究所 | Data processing method for vorticity-related observation flux data |
CN112035089B (en) * | 2020-08-20 | 2024-03-29 | 中国林业科学研究院林业研究所 | Data processing method for vorticity related observation flux data |
CN112136667A (en) * | 2020-11-26 | 2020-12-29 | 江苏久智环境科技服务有限公司 | Intelligent sprinkling irrigation method and system based on edge machine learning |
CN112136667B (en) * | 2020-11-26 | 2021-02-12 | 江苏久智环境科技服务有限公司 | Intelligent sprinkling irrigation method and system based on edge machine learning |
CN112507558B (en) * | 2020-12-15 | 2024-02-27 | 中国水利水电科学研究院 | Reference crop rising amount prediction method based on main control factor correction |
CN112507558A (en) * | 2020-12-15 | 2021-03-16 | 中国水利水电科学研究院 | Reference crop evapotranspiration amount prediction method based on main control factor correction |
CN114021486A (en) * | 2022-01-05 | 2022-02-08 | 中化现代农业有限公司 | Crop evapotranspiration prediction method and device |
CN114943361A (en) * | 2022-03-15 | 2022-08-26 | 水利部交通运输部国家能源局南京水利科学研究院 | Method for estimating evapotranspiration of reference crops in data-lacking areas |
CN115049126A (en) * | 2022-06-10 | 2022-09-13 | 中国水利水电科学研究院 | Evapotranspiration prediction method based on temperature effect and historical threshold |
CN115049126B (en) * | 2022-06-10 | 2024-06-21 | 中国水利水电科学研究院 | Evapotranspiration prediction method based on temperature effect and historical threshold |
CN115125903A (en) * | 2022-07-08 | 2022-09-30 | 中水三立数据技术股份有限公司 | Automatic irrigation and drainage management method for field integrated gate based on water demand prediction |
CN115125903B (en) * | 2022-07-08 | 2023-09-08 | 中水三立数据技术股份有限公司 | Automatic irrigation and drainage method for field integrated gate based on water demand prediction |
CN117678506A (en) * | 2024-02-02 | 2024-03-12 | 水利部牧区水利科学研究所 | Visual irrigation decision-making auxiliary system based on weather forecast and control method |
CN117678506B (en) * | 2024-02-02 | 2024-04-26 | 水利部牧区水利科学研究所 | Visual irrigation decision-making auxiliary system based on weather forecast and control method |
CN118366059A (en) * | 2024-06-20 | 2024-07-19 | 山东锋士信息技术有限公司 | Crop water demand calculating method based on optical and SAR data fusion |
Also Published As
Publication number | Publication date |
---|---|
CN111461909B (en) | 2023-02-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111461909B (en) | Short-term prediction method for farmland evapotranspiration | |
Talebmorad et al. | Evaluation of the impact of climate change on reference crop evapotranspiration in Hamedan-Bahar plain | |
CN110580657B (en) | Agricultural irrigation water demand prediction method | |
Gu et al. | Development of an irrigation scheduling software based on model predicted crop water stress | |
Choudhary et al. | Autonomous crop irrigation system using artificial intelligence | |
CN111833202B (en) | Farmland evapotranspiration short-term prediction method considering crop coefficient dynamic change and rainfall | |
CN111191193A (en) | Long-term soil temperature and humidity high-precision prediction method based on autoregressive moving average model | |
CN109447426B (en) | Response analysis method of irrigation water demand to changing environment based on crop water demand mechanism | |
CN109783774A (en) | A kind of temperature DATA PROCESSING IN ENSEMBLE PREDICTION SYSTEM method and system | |
Trajkovic | Comparison of radial basis function networks and empirical equations for converting from pan evaporation to reference evapotranspiration | |
CN111708050A (en) | Adaptive troposphere delay correction method and system for GNSS data processing | |
CN106355264B (en) | Reference crop evapotranspiration combined prediction method | |
Cáceres et al. | Economic model predictive control for smart and sustainable farm irrigation | |
CN116762539A (en) | Intelligent crop fertilizing method | |
Garcia-Mañas et al. | Leaf area index soft sensor for tomato crops in greenhouses | |
CN115912334A (en) | Method for establishing prediction model of output guarantee rate of wind power plant and prediction method | |
CN115049126B (en) | Evapotranspiration prediction method based on temperature effect and historical threshold | |
Khatri et al. | Water and salt balances at farmer fields | |
Arca et al. | Use of numerical weather forecast and time series models for predicting reference evapotranspiration | |
Lehner et al. | An improved statistical bias correction method that also corrects dry climate models | |
Nishitha et al. | Irrigation monitoring and controlling system | |
Baylon et al. | IoT-based Smart Irrigation System for Solanum Melongena or Eggplant using Evapotranspiration with Penman-Monteith Equation | |
Carey-Smith et al. | Changes in precipitation extremes for New Zealand: climate model predictions | |
US20230270060A1 (en) | System and method for managing irrigation of a field plot for a growing season | |
US12106231B2 (en) | Method and apparatus for predictive calculation of plant water need |
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 |