CN110705811A - Method and device for predicting yield based on water level burial depth value - Google Patents
Method and device for predicting yield based on water level burial depth value Download PDFInfo
- Publication number
- CN110705811A CN110705811A CN201810752216.1A CN201810752216A CN110705811A CN 110705811 A CN110705811 A CN 110705811A CN 201810752216 A CN201810752216 A CN 201810752216A CN 110705811 A CN110705811 A CN 110705811A
- Authority
- CN
- China
- Prior art keywords
- value
- stress
- burial depth
- root zone
- root
- 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
- 238000009933 burial Methods 0.000 title claims abstract description 137
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000035882 stress Effects 0.000 claims abstract description 119
- 230000008641 drought stress Effects 0.000 claims abstract description 74
- 230000006870 function Effects 0.000 claims description 64
- 230000008859 change Effects 0.000 claims description 20
- 238000003306 harvesting Methods 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 4
- 239000002131 composite material Substances 0.000 claims 1
- 238000012271 agricultural production Methods 0.000 abstract description 3
- 241000209140 Triticum Species 0.000 description 8
- 235000021307 Triticum Nutrition 0.000 description 8
- 238000004519 manufacturing process Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000006735 deficit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000009331 sowing Methods 0.000 description 2
- 208000005156 Dehydration Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 235000015816 nutrient absorption Nutrition 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000005068 transpiration Effects 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
- 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/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the invention discloses a method and a device for predicting yield based on a water level burial depth value and a computer-storable medium, and belongs to the technical field of agricultural production. The method comprises the following steps: obtaining an effective root zone buried depth value of the crop at the current stage; acquiring a waterlogging stress value and a drought stress value based on the effective root zone buried depth value; acquiring a water level buried depth value, a waterlogging stress coefficient and a drought stress coefficient of a crop growth area aiming at crops; acquiring a comprehensive stress index based on the water level burial depth value, the drought stress value, the water logging stress coefficient and the drought stress coefficient; and acquiring the relative yield of the crops based on the comprehensive stress index. By adopting the method and the device, the relative yield of the crops can be more accurately predicted based on the water level buried depth value of the area where the crops are located and considering the influence of drought and waterlogging on the yield.
Description
Technical Field
The invention relates to the technical field of agricultural production, in particular to a method and a device for predicting yield based on a water level burial depth value and a computer-storable medium.
Background
Agricultural waterlogging disasters and drought hazards are important factors for restricting food production in China. With the aggravation of extreme climate change and the development of farmland drainage engineering technology, the new forms of drought and waterlogging rush to turn and drought and waterlogging coexist further aggravate the yield loss of crops. The moisture production function is an effective tool to describe the impact of water deficit on crop yield.
At present, a water production function for describing water profit and loss is to establish a functional relation between waterlogging and yield through underground water level and yield, and establish a functional relation between drought and yield through transpiration evaporation capacity and yield so as to judge the drought and waterlogging degree of soil and further directly predict the crop yield.
In the process of implementing the invention, the inventor finds that the prior art has at least the following problems:
the traditional water production function separately establishes a functional relation representing water surplus and deficit stress respectively aiming at waterlogging and drought disasters, and cannot well represent the crop yield response relation under the condition of alternate occurrence of drought and waterlogging stress in the whole crop growth period, so that the traditional crop water production function cannot well predict crop yield and evaluate the influence of the waterlogging and drought disasters on grain production.
Disclosure of Invention
In order to solve the problem of overlarge prediction error of crop yield, the embodiment of the invention provides a method, a device and a computer storage medium for predicting yield based on a water level burial depth value, and the relative yield of crops can be directly predicted by measuring the water level burial depth value. The technical scheme is as follows:
in a first aspect, a method for predicting yield based on a water level burial depth value is provided, and the method comprises the following steps:
obtaining an effective root zone buried depth value of the crop at the current stage;
acquiring a waterlogging stress value and a drought stress value based on the effective root zone buried depth value, wherein the waterlogging stress value is the maximum water level buried depth value of crops which are not subjected to waterlogging stress, and the drought stress value is the minimum water level buried depth value of crops which are not subjected to drought stress in the growing period;
acquiring a water level buried depth value of a crop growth area aiming at crops;
acquiring a comprehensive stress index based on the water level burial depth value, the drought stress value and the waterlogging stress value;
and acquiring the relative yield of the crops based on the comprehensive stress index.
Optionally, obtaining the effective root zone burial depth value of the crop at the current stage includes:
the growth time of the crops is obtained,
inputting the growth time of the crops into an effective root zone time change function to obtain an effective root zone burial depth value, wherein the effective root zone time change function is
Wherein, REDtThe effective root zone burial depth value is represented by a unit cm, H is the maximum root system burial depth of the crops, T is the growth period duration, d (day) is represented by a unit cm, T is the crop growth time, d (day) is represented by a unit e is a natural logarithm, and k is the effective root zone burial depth and the actual root zone burial depth coefficient.
Optionally, obtaining the effective root zone burial depth value of the crop at the current stage includes:
inputting the growth time into a root region burial depth time variation function to obtain the root region burial depth, wherein the root region burial depth time variation function is
Wherein RDtIs the root zone burial depth, t is the crop growth time, unit d; t is the growth period duration, unit d, and H is the maximum root system burial depth, unit cm; e is a natural logarithm;
obtaining an effective root region burial depth coefficient;
and obtaining an effective root zone burial depth value based on the effective root zone burial depth coefficient and the root zone burial depth.
Optionally, inputting the growth time into a root region burial depth time variation function to obtain the root region burial depth, including:
inputting the crop growth time into a change function of the relative root depth and the time to obtain the relative root depth, wherein the change function of the relative root depth and the time is
Wherein RDxThe relative root depth represents the ratio of the actual root length to the maximum root length, and x is relative time, namely the ratio of the growth time of the crops to the growth period;
acquiring the maximum root system burial depth and the growth period duration;
the buried depth of the root zone is obtained based on the relative root depth, the maximum root buried depth and the growth period duration.
Optionally, acquiring a water logging stress value and a drought stress value based on the effective root zone buried depth value includes:
acquiring the effective root zone submergence rate of crops;
substituting the obtained effective root zone burial depth value into a waterlogging stress function to obtain the waterlogging stress value, wherein the waterlogging stress function is
WTWL(t)=ωwlREDt;
Wherein WTWL (t) is a waterlogging stress function, REDtAs a function of the effective root zone burial depth, omegawlIndicating the effective root zone flooding.
Optionally, acquiring a water logging stress value and a drought stress value based on the effective root zone buried depth value includes:
acquiring drought stress parameters;
substituting the obtained effective root zone burial depth value into a drought stress function to obtain the drought stress value, wherein the drought stress function is
WTDR(t)=REDt+ωdr;
Wherein WTDR (t) is a drought stress function, REDtIs buried in the effective root zoneFunction of depth change, ωdrRepresents the drought stress parameter in cm.
Optionally, obtaining a comprehensive stress index based on the water level burial depth value, the drought stress value and the water logging stress value includes:
substituting the water level buried depth value, the drought stress value and the waterlogging stress value into a comprehensive stress index function to obtain a comprehensive stress index, wherein the comprehensive stress index function is
Wherein SDI is comprehensive stress index, h (t) is actual water level buried depth value, unit cm belongs towlBelongs to the waterlogging stress coefficientdrIs drought stress coefficient, time t is crop growth time, unit d, n is crop growth period duration, unit d.
Optionally, obtaining relative crop yield based on the integrated stress index comprises:
substituting the comprehensive stress index into a relative yield relational expression to obtain the relative yield, wherein the relative yield relational expression is
Where RY is the relative yield, SDI0Is the comprehensive stress index value, SDI, of crops under the condition of no harvest100Is the value of the comprehensive stress index when the relative yield of the crops is the highest yield of the crops.
In a second aspect, there is provided an apparatus for predicting yield based on a water level burial depth value, the apparatus comprising:
a processor;
a memory for storing at least one computer-executable instruction;
wherein the processor is configured to:
the method of the first aspect is implemented when the at least one computer-executable instruction is executed.
In a third aspect, there is provided a computer readable storage medium having stored therein computer instructions which, when executed by a processor, implement the method of any one of claims 1-8.
The technical scheme provided by the embodiment of the invention has the following beneficial effects:
in the embodiment of the invention, the effective root zone burial depth value of the crop at the current stage is obtained; acquiring a waterlogging stress value and a drought stress value based on the effective root zone buried depth value; acquiring a water level buried depth value, a waterlogging stress coefficient and a drought stress coefficient of a crop growth area aiming at crops; acquiring a comprehensive stress index based on the water level burial depth value, the drought stress value, the water logging stress coefficient and the drought stress coefficient; and acquiring the relative yield of the crops based on the comprehensive stress index. Therefore, the root zone growth and the root zone effective growth change in the crop growth period are introduced into the function as variables, the fitting degree of drought and waterlogging on the crop growth is improved, the influence of drought and waterlogging on the crop yield is effectively monitored, and the yield prediction accuracy of crops is improved. By adopting the method and the device, the relative yield of the crops under the influence of the water level buried depth value can be more accurately predicted based on the water level buried depth value of the area where the crops are located and considering drought and waterlogging.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a flowchart of a method for predicting yield based on a water level burial depth value according to an embodiment of the invention;
FIG. 2 is a schematic diagram of drought and flood stress lines during the crop growth period provided by an embodiment of the present invention;
FIG. 3 is a schematic diagram of a relation between drought and flood stress lines and effective root zone burial depth provided by an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
Example one
An embodiment of the present invention provides a method for predicting yield based on a water level burial depth value, as shown in fig. 1, a processing flow of the method may include the following steps:
In the agricultural production process, the root zone of crops is an important part of nutrient absorption, and the burial depth of the root zone of the crops is different at different periods in the growth period. The root zone buried depth value refers to the height from the tail end of the root zone to the ground in the growth process of crops, and the water level buried depth value of the zone where the crops are located has great influence on the root zone of the crops, namely, the difference of the water level buried depth value can influence the growth environment of the crops, for example, the water level buried depth is shallow, water logging is caused, and the water level buried depth is deep, drought is caused, and the like. As the water stress of crops (namely the influence of drought and waterlogging on the crops) is closely related to the range of root zones, namely the influence of water on different root zones is different, the root zone with the most obvious drought and waterlogging stress is taken as an effective root zone. The growth data of the crops such as the root zone lengths and the root zone burial depth values of different environments at different stages of the growth period of the crops are relatively fixed, so the effective root zone burial depth value can be obtained from historical data, for example, the growth data of the crops such as the root lengths and the stem lengths of the crops at different periods can be consulted, and the effective root zone data can be obtained from the root zone data. The effective root zone buried depth value of the crop at the current stage can be obtained by sampling statistics in the field or other modes. The buried depth value of the effective root zone is the distance value between the earth surface and the lowest end of the effective root zone. Where the units of measure mentioned herein may be transformed in actual practice based on the habits or needs of the skilled person.
Specifically, the change relation between the relative root zone burial depth and the relative time can be constructed according to the change rule of the root zone burial depth along with the growth of the crops, and the change function of the relative root zone burial depth and the relative time is obtained
Wherein RDxThe relative root depth represents the ratio of the actual root length to the maximum root length, and the value range is 0-1; and x is relative time which represents the ratio of the actual time after crop sowing to the growth period, and the value range is 0-1. The actual time length after the crop is sown is the crop growth time, and the relative root depth of the current crop can be obtained by inputting the ratio of the crop growth time to the crop growth period time length into the formula (1).
The maximum root system burial depth and the growth period of the crops can be obtained from related research data or field monitoring, the maximum root system burial depth and the growth period are introduced into the formula (1), and the change function of the burial depth value and the time of the root zone is obtained
Wherein RDtIs the root zone buried depth, t is the crop growth time, unit d (day); t is the growth period duration, unit d (day), H is the maximum root system burial depth, unit cm; e is the natural logarithm. The length of the growing period may be counted starting on the sowing day and ending on the maturation or harvest date. The relative root depth obtained in the formula (1) and the obtained maximum root length are substituted into the formula (2) to obtain the root region burial depth of the crop.
The regions sensitive to the influence of water on the crop root zone are constantly changed and mostly concentrated on the upper part of the root zone, so that an effective region capable of reflecting the influence of drought and waterlogging stress on the crop growth is not equal to the whole root zone, the effective root zone burial depth and the actual root zone burial depth coefficient can be introduced to describe the relation between the effective root zone burial depth and the actual root zone burial depth, and a change function of the effective root zone burial depth value and time is obtained
Wherein, REDt is the buried depth value of the effective root zone and is unit cm; k is the effective root zone burial depth and the actual root zone burial depth coefficient, and the value range is 0-1; t is the crop growth time in d (days); t is the length of the birth period and the unit d (day); h is the maximum root system burial depth in cm; e is a natural logarithm and can take the value of 2.71828, and the value can be taken according to actual requirements in actual operation and is not limited here.
And (3) substituting the root zone burial depth obtained in the formula (2), the introduced effective root zone burial depth and the actual root zone burial depth coefficient into the formula (3) to obtain the effective root zone burial depth value of the crop. Or the effective root zone burial depth and actual root zone burial depth coefficients, the crop growth time, the maximum root system burial depth and the crop growth period duration can be directly input into the formula (3) to directly obtain the effective root zone burial depth value of the crop.
And 102, acquiring a water logging stress value and a drought stress value based on the effective root zone buried depth value.
The drought stress value is the minimum water level buried depth value which is not stressed by drought in the crop growth period; the waterlogging stress value is the maximum water level buried depth value of crops which are not stressed by waterlogging.
As shown in FIGS. 2 and 3, the corresponding relationship between the buried depth of the effective root zone and the water-logging stress value and the drought stress value, i.e., the relationship between the buried depth of the effective root zone and the water-logging stress value and the drought stress value, can be obtained according to the historical data of the crop growth, i.e., the relationship between the buried depth of the
WTWL(t)=ωwlREDt(4),
WTDR(t)=REDt+ωdr(5),
Wherein WTWL (t) is a waterlogging stress value; WTDR (t) is drought stress value; REDtIs the effective root zone buried depth value; omegawlThe effective root zone submergence rate is represented, and the value range is 0-1; omegadrRepresents the drought stress parameter in cm. The effective root zone inundation rate may be a ratio of an effective root zone range inundated with groundwater to an effective root zone range.
Optionally, the method for obtaining the submergence rate of the effective root zone may be to sample and measure the submergence rate of the effective root zone obtained by the crop in the field, or to measure the water level burial depth value, and then to reversely calculate the submergence rate of the effective root zone based on the burial depth value of the effective root zone obtained in the formula (3).
The drought stress parameter measuring method can be obtained by reversely calculating the obtained effective root zone buried depth value after measuring the water level buried depth value.
And (4) respectively inputting the effective root zone buried depth value obtained in the formula (3) into the formulas (4) and (5) to obtain a waterlogging stress value and a drought stress value of the crop.
And 103, acquiring a water level buried depth value, a water logging stress coefficient and a drought stress coefficient of the crop growth area aiming at the crop.
In the determined area, the water level burial depth value is constant and can be obtained through measurement. For a given crop, the waterlogging stress coefficient and the drought stress coefficient are determined and can be obtained by calibration according to the growth data of the crop in the area over the years, namely, the growth data of the over the years are counted and analyzed, the growth data are substituted into the growth model of the crop, the calculation result is compared with the measured value, and the calculation result meeting the error with the measured value is used as the calibration coefficient.
And 104, acquiring a comprehensive stress index based on the water level burial depth value, the drought stress value and the water logging stress value.
The comprehensive stress index considering both waterlogging and drought stress can be used to obtain a calculation formula based on the comprehensive stress index, i.e.
Wherein SDI is a comprehensive stress index, h (t) is the actual water level burial depth which takes the ground surface as a zero point and is negative downwards, the unit is cm, SEW is a waterlogging stress index, and belongs towlIs waterlogging stress coefficient, SDW is drought stress index, belongs todrIs drought stress coefficient, time t is time step, unit d is crop growth time, and n is crop growth time, i.e. growth period duration of crop, unit d.
Therefore, the obtained water level burial depth value, drought stress value and waterlogging stress value are substituted into the formula (6), and a comprehensive stress index can be calculated.
And 105, acquiring the relative yield of the crops based on the comprehensive stress index.
The corresponding relation between the relative yield and the comprehensive stress index can be obtained by counting and fitting the historical data of the crop growth
Where RY is the relative yield, SDI0The comprehensive stress index value of the crop under the condition of no harvest completely, namely the relative yield of the crop is 0 percent, and the SDI100The value is the value of the integrated stress index at which the relative yield of the crop is the highest, i.e., the relative yield of the crop is 100%.
The relative yield of the crop is 0% and 100% is a known parameter for the crop, and the predicted value of the relative yield of the crop can be obtained by inputting the comprehensive stress index obtained in the formula (6) into the formula (7).
Optionally, taking wheat as an example to illustrate the construction process of each function in the method.
Determining an effective root region burial depth variation function:
the maximum root system burial depth H of the wheat can be determined to be-100 cm according to relevant documents and field monitoring data of the wheat, the crop growth period T is 210 days, and the coefficient k of the effective root area burial depth and the actual root area burial depth is 0.5. Substituting the effective root zone burial depth value into a time change function, namely formula (3), to obtain the time change function of the effective root zone of wheat
The buried depth value of the effective root zone of the wheat can be obtained by inputting the growth time of the crops into a formula (8).
Wherein REDt is the buried depth value of the effective root zone, t is the growth time of the crops, and e is the natural logarithm.
Calculating a comprehensive stress index based on the water level burial depth measured data:
assuming that the waterlogging stress coefficient and the drought stress coefficient are in the same order aswl=∈drUnder the condition of 1.0, the stress index SDI is calculated by using the water level burial depth h (t) data monitored in the field and taking day as a time step in the whole growth period by adopting the formulas (4), (5) and (6). Wherein the waterlogging stress coefficient and the drought stress coefficient can also be based on wheat in the region over the yearsThe production data are obtained by calculation.
Optionally, the unknown parameter ω can be obtained by calibration based on the production data of wheat in the region over the yearswl=0.3、ωdr=30cm、SDI0727cm d and SDI10060cm d, etc. After introducing calibration parameters, the function is constructed as
WTWL(t)=0.3REDt(9),
WTDR(t)=REDt-30 (10),
Namely, finally, only the water level burial depth value is needed to be input, and the relative yield of wheat under the influence of the water level burial depth can be obtained.
In the embodiment of the invention, the effective root zone burial depth value of the crop at the current stage is obtained; acquiring a waterlogging stress value and a drought stress value based on the effective root zone buried depth value; acquiring a water level buried depth value, a waterlogging stress coefficient and a drought stress coefficient of a crop growth area aiming at crops; acquiring a comprehensive stress index based on the water level burial depth value, the drought stress value, the water logging stress coefficient and the drought stress coefficient; and acquiring the relative yield of the crops based on the comprehensive stress index. The root zone growth in the crop growth period and the change of the effective root zone based on the growth of the root zone are introduced into the function as variables, so that the fitting degree of drought and waterlogging on the crop growth is improved, the influence of drought and waterlogging on the crop yield is effectively monitored, and the yield prediction accuracy of crops is effectively improved.
Example two
The embodiment of the invention provides a device for predicting yield based on a water level burial depth value, which comprises:
a processor;
a memory for storing at least one computer-executable instruction;
wherein the processor is configured to:
upon execution of the at least one computer-executable instruction, implementing:
obtaining an effective root zone buried depth value of the crop at the current stage;
acquiring a waterlogging stress value and a drought stress value based on the effective root zone buried depth value, wherein the waterlogging stress value is the maximum water level buried depth value of crops which are not subjected to waterlogging stress, and the drought stress value is the minimum water level buried depth value of crops which are not subjected to drought stress in the growing period;
acquiring a water level buried depth value of a crop growth area aiming at crops;
acquiring a comprehensive stress index based on the water level burial depth value, the drought stress value and the waterlogging stress value;
and acquiring the relative yield of the crops based on the comprehensive stress index.
Optionally, obtaining the effective root zone burial depth value of the crop at the current stage includes:
the growth time of the crops is obtained,
inputting the crop growth time into an effective root zone time change function to obtain an effective root zone burial depth value;
effective root zone time variation function of
Wherein, REDtThe effective root zone burial depth value is expressed in the unit of cm, H is the maximum root system burial depth of the crops, T is the growth period duration in the unit of cm, T is the growth time of the crops, d and e are natural logarithm, and k is the effective root zone burial depth and the actual root zone burial depth coefficient.
Optionally, obtaining the effective root zone burial depth value of the crop at the current stage includes:
inputting the growth time into a root region burial depth time variation function to obtain the root region burial depth;
the root region burial depth time variation function is
Wherein RDtIs the root zone burial depth, t is the crop growth time, unit d; t is the growth period duration, unit d, and H is the maximum root system burial depth, unit cm; e is a natural logarithm;
obtaining an effective root region burial depth coefficient;
and obtaining an effective root zone burial depth value based on the effective root zone burial depth coefficient and the root zone burial depth.
Optionally, inputting the growth time into a root region burial depth time variation function to obtain the root region burial depth, including:
inputting the crop growth time into a change function of the relative root depth and the time to obtain the relative root depth, wherein the change function of the relative root depth and the time is
Wherein RDxThe relative root depth represents the ratio of the actual root length to the maximum root length, and x is relative time, namely the ratio of the growth time of the crops to the growth period;
acquiring the maximum root system burial depth and the growth period duration;
the buried depth of the root zone is obtained based on the relative root depth, the maximum root buried depth and the growth period duration.
Optionally, acquiring a water logging stress value and a drought stress value based on the effective root zone buried depth value includes:
acquiring the effective root zone submergence rate of crops;
substituting the obtained effective root zone burial depth value into a waterlogging stress function to obtain the waterlogging stress value, wherein the waterlogging stress function is
WTWL(t)=ωwlREDt;
Wherein WTWL (t) is a waterlogging stress function, REDtAs a function of the effective root zone burial depth, omegawlIndicating the effective root zone flooding.
Optionally, acquiring a water logging stress value and a drought stress value based on the effective root zone buried depth value includes:
acquiring drought stress parameters;
substituting the obtained effective root zone burial depth value into a drought stress function to obtain the drought stress value, wherein the drought stress function is
WTDR(t)=REDt+ωdr;
Wherein WTDR (t) is a drought stress function, REDtAs a function of the effective root zone burial depth, omegadrRepresents the drought stress parameter in cm.
Optionally, obtaining a comprehensive stress index based on the water level burial depth value, the drought stress value and the water logging stress value includes:
substituting the water level buried depth value, the drought stress value and the waterlogging stress value into a comprehensive stress index function to obtain a comprehensive stress index, wherein the comprehensive stress index function is
Wherein SDI is comprehensive stress index, h (t) is actual water level buried depth value, unit cm belongs towlBelongs to the waterlogging stress coefficientdrIs drought stress coefficient, time t is crop growth time, unit d, n is crop growth period duration, unit d.
Optionally, obtaining relative crop yield based on the integrated stress index comprises:
substituting the comprehensive stress index into a relative yield relational expression to obtain the relative yield, wherein the relative yield relational expression is
Where RY is the relative yield, SDI0Is the comprehensive stress index value, SDI, of crops under the condition of no harvest100Is the value of the comprehensive stress index when the relative yield of the crops is the highest yield of the crops.
EXAMPLE III
A computer-readable storage medium having stored thereon computer instructions, which when executed by a processor, implement the method of embodiment one.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (10)
1. A method for predicting yield based on a water level burial depth value, which is characterized by comprising the following steps:
obtaining an effective root zone buried depth value of the crop at the current stage;
acquiring a waterlogging stress value and a drought stress value based on the effective root zone buried depth value, wherein the waterlogging stress value is the maximum water level buried depth value of the crops which are not subjected to waterlogging stress, and the drought stress value is the minimum water level buried depth value of the crops which are not subjected to drought stress in the growth period;
acquiring a water level buried depth value of the crop growth area aiming at the crop;
acquiring a comprehensive stress index based on the water level burial depth value, the drought stress value and the water logging stress value;
obtaining the relative crop yield based on the integrated stress index.
2. The method of claim 1, wherein the obtaining the effective root zone burial depth value of the current stage crop comprises:
the growth time of the crops is obtained,
inputting the crop growth time into the effective root zone time variation function to obtain the buried depth value of the effective root zone, wherein the effective root zone time variation function is
Wherein, REDtThe effective root zone burial depth value is represented by the unit cm, H is the maximum root system burial depth of the crops, T is the growth period duration, d and T are the crop growth time, d and e are the natural logarithm, and k is the effective root zone burial depth and the actual root zone burial depth coefficient.
3. The method of claim 1, wherein the obtaining the effective root zone burial depth value of the current stage crop comprises:
inputting the growth time into a root region burial depth time variation function to obtain the root region burial depth, wherein the root region burial depth time variation function is
Wherein RDtIs the root zone burial depth, t is the crop growth time, unit d; t is the growth period duration, unit d, and H is the maximum root system burial depth, unit cm; e is a natural logarithm;
obtaining an effective root region burial depth coefficient;
and acquiring the effective root zone burial depth value based on the effective root zone burial depth coefficient and the root zone burial depth.
4. The method of claim 3, wherein inputting the growth time into a root region burial depth time variation function to obtain a root region burial depth comprises:
inputting the crop growth time into a change function of the relative root depth and the time to obtain the relative root depth, wherein the change function of the relative root depth and the time is
Wherein RDxThe relative root depth represents the ratio of the actual root length to the maximum root length, and x is relative time, namely the ratio of the growth time of the crops to the growth period;
acquiring the maximum root system burial depth and the growth period duration;
and acquiring the buried depth of the root zone based on the relative root depth, the maximum root system buried depth and the growth period duration.
5. The method according to claim 2, wherein the obtaining of water-logging stress values and drought stress values based on the effective root zone burial depth value comprises:
obtaining the effective root zone submergence rate of the crops;
substituting the obtained effective root zone burial depth value into the water logging stress function to obtain the water logging stress value, wherein the water logging stress function is
WTWL(t)=ωwlREDt;
Wherein WTWL (t) is a waterlogging stress function, REDtAs a function of the effective root zone burial depth, omegawlIndicating the effective root zone flooding.
6. The method according to claim 5, wherein the obtaining of water-logging stress values and drought stress values based on the effective root zone burial depth value comprises:
acquiring drought stress parameters;
substituting the obtained effective root zone burial depth value into the drought stress function to obtain the drought stress value, wherein the drought stress function is
WTDR(t)=REDt+ωdr;
Wherein WTDR (t) is a drought stress function, REDtAs a function of the effective root zone burial depth, omegadrRepresents the drought stress parameter in cm.
7. The method according to claim 6, wherein the obtaining a composite stress index based on the water level burial depth value, the drought stress value and the water logging stress value comprises:
substituting the water level burial depth value, the drought stress value and the water logging stress value into a comprehensive stress index function to obtain the comprehensive stress index, wherein the comprehensive stress index function is
Wherein SDI is comprehensive stress index, h (t) is actual water level buried depth value in cm,∈wlBelongs to the waterlogging stress coefficientdrIs drought stress coefficient, time t is crop growth time, unit d, n is crop growth period duration, unit d.
8. The method of claim 7, wherein said obtaining relative crop yield based on said combined stress index comprises:
substituting the comprehensive stress index into a relative yield relational expression to obtain the relative yield, wherein the relative yield relational expression is
Where RY is the relative yield, SDI0Is the comprehensive stress index value, SDI, of crops under the condition of no harvest100Is the value of the comprehensive stress index when the relative yield of the crops is the highest yield of the crops.
9. An apparatus for predicting yield based on a water level burial depth value, the apparatus comprising:
a processor;
a memory for storing at least one computer-executable instruction;
wherein the processor is configured to:
the method of any of claims 1 to 8 when executed by the at least one computer-executable instruction.
10. A computer-readable storage medium having stored therein computer instructions which, when executed by a processor, implement the method of any one of claims 1-8.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810752216.1A CN110705811B (en) | 2018-07-10 | 2018-07-10 | Method and device for predicting yield based on water level burial depth value |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810752216.1A CN110705811B (en) | 2018-07-10 | 2018-07-10 | Method and device for predicting yield based on water level burial depth value |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110705811A true CN110705811A (en) | 2020-01-17 |
CN110705811B CN110705811B (en) | 2022-05-13 |
Family
ID=69192825
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810752216.1A Active CN110705811B (en) | 2018-07-10 | 2018-07-10 | Method and device for predicting yield based on water level burial depth value |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110705811B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113449984A (en) * | 2021-06-25 | 2021-09-28 | 中国水利水电科学研究院 | Reed resource quantity evaluation method |
CN116227683A (en) * | 2023-01-09 | 2023-06-06 | 武汉大学 | Evaluation method for drought and waterlogging tight interaction effect in easy drought and waterlogging area |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000283800A (en) * | 1999-03-31 | 2000-10-13 | Mitsui Bussan Plant Kk | Physical geographic displacement detector and physical geographic displacement monitoring system employing it |
CN107122901A (en) * | 2017-04-26 | 2017-09-01 | 中国水利水电科学研究院 | A kind of determination method of the crucial growing stage of region of no relief field-crop draining |
CN107480888A (en) * | 2017-08-17 | 2017-12-15 | 中国水利水电科学研究院 | A kind of agricultural drought disaster methods of risk assessment based on APSIM models |
CN107918135A (en) * | 2017-11-13 | 2018-04-17 | 中国科学院寒区旱区环境与工程研究所 | Water stress state monitoring method, device and electronic equipment |
-
2018
- 2018-07-10 CN CN201810752216.1A patent/CN110705811B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000283800A (en) * | 1999-03-31 | 2000-10-13 | Mitsui Bussan Plant Kk | Physical geographic displacement detector and physical geographic displacement monitoring system employing it |
CN107122901A (en) * | 2017-04-26 | 2017-09-01 | 中国水利水电科学研究院 | A kind of determination method of the crucial growing stage of region of no relief field-crop draining |
CN107480888A (en) * | 2017-08-17 | 2017-12-15 | 中国水利水电科学研究院 | A kind of agricultural drought disaster methods of risk assessment based on APSIM models |
CN107918135A (en) * | 2017-11-13 | 2018-04-17 | 中国科学院寒区旱区环境与工程研究所 | Water stress state monitoring method, device and electronic equipment |
Non-Patent Citations (1)
Title |
---|
于颖多、焦平金: ""排水循环灌溉驱动的稻区水循环模型与评价"", 《农业工程学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113449984A (en) * | 2021-06-25 | 2021-09-28 | 中国水利水电科学研究院 | Reed resource quantity evaluation method |
CN113449984B (en) * | 2021-06-25 | 2023-12-29 | 中国水利水电科学研究院 | Reed resource quantity assessment method |
CN116227683A (en) * | 2023-01-09 | 2023-06-06 | 武汉大学 | Evaluation method for drought and waterlogging tight interaction effect in easy drought and waterlogging area |
CN116227683B (en) * | 2023-01-09 | 2023-11-17 | 武汉大学 | Evaluation method for drought and waterlogging tight interaction effect in easy drought and waterlogging area |
Also Published As
Publication number | Publication date |
---|---|
CN110705811B (en) | 2022-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3179319B1 (en) | Method for irrigation planning and system for its implementation | |
CN106028792B (en) | Method and system for crop management | |
Marek et al. | Post-processing techniques for reducing errors in weighing lysimeter evapotranspiration (ET) datasets | |
Dutta et al. | Model development in DNDC for the prediction of evapotranspiration and water use in temperate field cropping systems | |
US20150370935A1 (en) | Agronomic systems, methods and apparatuses | |
Buitink et al. | Anatomy of the 2018 agricultural drought in the Netherlands using in situ soil moisture and satellite vegetation indices | |
US20150272017A1 (en) | Method and system for automated differential irrigation | |
Cichota et al. | Ensemble pedotransfer functions to derive hydraulic properties for New Zealand soils | |
Shafiei et al. | Uncertainty assessment of the agro-hydrological SWAP model application at field scale: A case study in a dry region | |
US20190171990A1 (en) | Using historical plant-available water metrics to forecast crop yield | |
Martí et al. | An artificial neural network approach to the estimation of stem water potential from frequency domain reflectometry soil moisture measurements and meteorological data | |
CN110705811B (en) | Method and device for predicting yield based on water level burial depth value | |
y Garcia et al. | ENSO-based climate variability affects water use efficiency of rainfed cotton grown in the southeastern USA | |
CN115344815B (en) | Natural runoff change attribution analysis method and system considering vegetation space change | |
Guest et al. | Comparing the performance of the DNDC, Holos, and VSMB models for predicting the water partitioning of various crops and sites across Canada | |
CN101413935A (en) | Method for (in-situ) determining steam scatter amount of crops | |
Fry et al. | The effect of up-scaling soil properties and model parameters on predictive accuracy of DSSAT crop simulation model under variable weather conditions | |
Liakos et al. | A decision support tool for managing precision irrigation with center pivots | |
Bouksila et al. | Electromagnetic induction prediction of soil salinity and groundwater properties in a Tunisian Saharan oasis | |
Togbévi et al. | A multi-model approach for analysing water balance and water-related ecosystem services in the Ouriyori catchment (Benin) | |
Xie et al. | A new method to estimate the cover and management factor for soil loss prediction on the Loess Plateau in China: A case‐study using a soybean field | |
Yuan et al. | Evaluation of soil erosion and sediment yield from ridge watersheds leading to Guánica Bay, Puerto Rico, using the soil and water assessment tool model | |
Reddy et al. | Crop coefficients of some selected crops of Andhra Pradesh | |
Arslan | Spatial and temporal distribution of areas with drainage problems as estimated by different interpolation techniques | |
Shi et al. | Global dryland aridity changes indicated by atmospheric, hydrological, and vegetation observations at meteorological stations |
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 |