CN105447317B - The analysis method of crop climate yield potentiality - Google Patents

The analysis method of crop climate yield potentiality Download PDF

Info

Publication number
CN105447317B
CN105447317B CN201510867342.8A CN201510867342A CN105447317B CN 105447317 B CN105447317 B CN 105447317B CN 201510867342 A CN201510867342 A CN 201510867342A CN 105447317 B CN105447317 B CN 105447317B
Authority
CN
China
Prior art keywords
day
msub
mrow
potentiality
climate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201510867342.8A
Other languages
Chinese (zh)
Other versions
CN105447317A (en
Inventor
魏晓文
雷亚平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Cotton Research of Chinese Academy of Agricultural Sciences
Original Assignee
Institute of Cotton Research of Chinese Academy of Agricultural Sciences
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Cotton Research of Chinese Academy of Agricultural Sciences filed Critical Institute of Cotton Research of Chinese Academy of Agricultural Sciences
Priority to CN201510867342.8A priority Critical patent/CN105447317B/en
Publication of CN105447317A publication Critical patent/CN105447317A/en
Application granted granted Critical
Publication of CN105447317B publication Critical patent/CN105447317B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Breeding Of Plants And Reproduction By Means Of Culturing (AREA)

Abstract

The present invention provides a kind of analysis method of crop climate yield potentiality, including:Build crop climate Potential Production Analysis model;When needing to analyze in cotton development stage total climate yields potentiality corresponding to specified time length, start time data and deadline data are inputted in the climate yields Potentials model obtained to step 1 structure;Total climate yields potentiality in the output analyzed period.Advantage is:(1) be calculated using function approximation day water requirement, day moisture supply and demand unspent amount and daily precipitation correction coefficient, more meet the actual conditions of the Sino-Japan water requirement of process of crop growth;(2) total climate yields potentiality are obtained by the accumulation of day climate yields potentiality, can reflects the forming process of crop climate yield potentiality, accurately reflect crop dry-matter accumulation state, improve the accuracy and objectivity of climate yields potentiality.

Description

The analysis method of crop climate yield potentiality
Technical field
The invention belongs to agricultural resource studying technological domain, and in particular to a kind of analysis side of crop climate yield potentiality Method.
Background technology
Crop climate yield potentiality refers to abundant and reasonable using local light, heat, water climate resources, and other The highest biological yield or farm output obtained when part, such as soil, nutrient are in optimum situation on unit area soil. Therefore, solar radiation energy, crop photosynthesis characteristic, temperature condition become the Fundamentals that crop yield is formed.Weather produces Amount potentiality reflect the basic potentiality that somewhere agricultural production under specific Climatic possesses, can quantitatively characterizing regional climate The configuration of resource situation and its climatic elements, is the weight that science weighs region grain-production power, agricultural development and population bearing capacity Want one of index.
The computational methods of traditional climate yields potentiality are:
Yw=YT*f(W);
Wherein, YwFor climate yields potentiality total in Crop growing stage;YTDive for light temperature yield total in Crop growing stage Power;F (W) is moisture correction factor.
As it can be seen that when calculating climate yields potentiality, the definite of moisture correction factor belongs to crucial.And conventional aqueous amendment system Number f (W) is determined by the following method:
Wherein, R is gross precipitation in Crop growing stage, unit mm, E0For total evaporation in Crop growing stage, unit is mm。
It can be seen from the above that in the calculating of conventional aqueous correction factor, since the parameter of consideration is total precipitation in Crop growing stage Total evaporation in amount and Crop growing stage, therefore, the moisture correction factor being calculated is actually that Crop growing stage is always corrected Coefficient, can not reflect due to the differentiation of daily service position and to the Affected By Gradual Change process of climate yields potentiality, with crop dry Matter productive accumulation process is not inconsistent, so, the accuracy of the climate yields potentiality obtained using conventional method is very limited, usually with Actual conditions deviation is serious.
The content of the invention
In view of the defects existing in the prior art, the present invention provides a kind of analysis method of crop climate yield potentiality, can have Effect improves the accuracy and objectivity for the climate yields potentiality that analysis obtains.
The technical solution adopted by the present invention is as follows:
The present invention provides a kind of analysis method of crop climate yield potentiality, comprises the following steps:
Step 1, crop climate Potential Production Analysis model is built, the crop climate Potential Production Analysis model represents For:
Yw=∑ Ywi (1)
Wherein:
Ywi=ri·yi (2)
pi=0.6pi-1-wi+di (4)
Wherein:YwTotal climate yields potentiality in-analysis Period Length;
YwiDay climate yields potentiality in-analysis Period Length;
ri- daily precipitation correction coefficient;
yi- sunlight yield potentiality;
pi- day moisture supply and demand unspent amount;
pi-1Moisture supply and demand-yesterday unspent amount;
wi- day water requirement;
wm- this month first day water requirement;
wm+1- next month first day water requirement;
Day sequence in the n- months, its value is maximum to be no more than 31;
di- the same day actual precipitation;
Step 2, when needing to analyze in cotton development stage total climate yields potentiality corresponding to specified time length, Xiang Bu Start time data and deadline data are inputted in the climate yields Potentials model that rapid 1 structure obtains;
Step 3, after climate yields Potentials model receives start time data and deadline data, by initial time Any one day is denoted as the jth moon i-th day in the period of deadline, then, reads pre-stored of that month first day water requirement wm With first day next month water requirement wm+1Value, and substitute into formula (5), a day water requirement w be calculatedi;Wherein, n=i;
Step 4, the same day actual precipitation d is readiWith moisture supply and demand unspent amount p yesterdayi-1, by day water requirement wi, the same day Actual precipitation diWith moisture supply and demand unspent amount p yesterdayi-1Formula (4) is substituted into, a day moisture supply and demand unspent amount p is calculatedi;This Place, the 1st day for being analyzed the period, it was 0 to make moisture supply and demand yesterday unspent amount;
Step 5, by day moisture supply and demand unspent amount piFormula (3) is substituted into, daily precipitation correction coefficient r is calculatedi
Step 6, pre-stored sunlight yield potentiality y is readi, by sunlight yield potentiality yiWith daily precipitation correction coefficient riFormula (2) is substituted into, a day climate yields potentiality Y is calculatedwi
Step 7, repeat step 3- steps 6, are calculated day climate yields potentiality daily in the analyzed period, so Afterwards, based on formula (1), the day climate yields potentiality in each day being calculated are done into summation operation, that is, obtain the analyzed period Interior total climate yields potentiality;
Step 8, total climate yields potentiality in the analyzed period that output step 7 obtains.
Preferably, in step 3, the value of pre-stored each first day moon water requirement is shown in Table 1:
Each first day moon water requirement of 1 cotton of table
Wherein, since the breeding time of cotton is the 5-9 months, so each first day moon water requirement in 1 storage 5-9 month of table.
Preferably, in step 6, the sunlight yield potentiality y that is readiObtain by the following method:
Step 6.1, pre-stored three base point temperature of the analyzed moon is read, is respectively preference temperature of the cotton in current period T1, the minimum temperature T that can bear2And the maximum temperature T that can be born3;A day actual average temperature T is read againi
Step 6.2, by T1、T2、T3And TiFormula (6) is substituted into, degree/day correction factor F is calculatedi
Wherein, B=(T3-T1)/(T1-T2);
Step 6.3, pre-stored day highest real air temperature T is readu, by day highest real air temperature TuSubstitute into formula (7), meter Calculation obtains daily maximum temperature correction factor Vi
Vi=(1-li2)25-29li 2
Wherein:
Step 6.4, by daily maximum temperature correction factor ViWith degree/day correction factor FiFormula (8) is substituted into, light is calculated Warm yield potentiality day correction factor ti
ti=0.3Vi+0.7Fi (8)
Step 6.5, photosynthetic yield potentiality h of pre-stored day is readi, by tiAnd hiFormula (9) is substituted into, daylight is calculated Warm yield potentiality yi
yi=ti·hi (9)。
Preferably, in step 6.5, the day photosynthetic yield potentiality h that is readiObtain by the following method:
Step 6.5.1, establishes each first day moon weight coefficient table in breeding time;
Step 6.5.2, reads each first day moon weight coefficient table, obtains of that month first day weight vmWith first day next month weight vm+1, then by of that month first day weight vmWith first day next month weight vm+1Formula (10) is substituted into, a day weight v is calculatedi
Wherein, n is day sequence, n=i in the moon;
Step 6.5.3, reads pre-stored day radiation G, day is radiated G and day weight viFormula (11) is substituted into, is calculated To day photosynthetic yield potentiality hi
hi=vi·G (11)。
Preferably, step 6.5.1, each first day moon weight coefficient table established are shown in Table 2:
The 2 cotton photosynthetic potentiality moon of table first day weight
Wherein, v4、v5、v6、v7、v8、v9And v10The respectively first day in April weight, the first day in May weight, the first day in June power Weight, the first day in July weight, August part first day weight, September part first day weight and the first day in October weight.
The analysis method of crop climate yield potentiality provided by the invention has the following advantages:
(1) be calculated using function approximation day water requirement, day moisture supply and demand unspent amount and daily precipitation correction coefficient, more Meet the actual conditions of the Sino-Japan water requirement of process of crop growth;
(2) total climate yields potentiality are obtained by the accumulation of day climate yields potentiality, can reflects crop climate yield potentiality Forming process, accurately reflect crop dry-matter accumulation state, improve climate yields potentiality accuracy and objectivity.
Brief description of the drawings
Fig. 1 is the contrast curve for the moisture update the system that conventional method and the method for the present invention are formed;
Fig. 2 is Cotton climatic output potentiality zoning map provided by the invention;
Fig. 3 is cotton light temperature yield potentiality zoning map provided by the invention;
Fig. 4 is cotton photosynthetic yield potentiality zoning map provided by the invention;
Fig. 5 is the variation diagram of 1971-1975 Cotton climatic outputs potentiality provided by the invention;
Fig. 6 is the variation diagram of 2002-2006 Cotton climatic outputs potentiality provided by the invention.
Embodiment
In order to which technical problem, technical solution and beneficial effect solved by the invention is more clearly understood, below in conjunction with Accompanying drawings and embodiments, the present invention will be described in further detail.It should be appreciated that specific embodiment described herein only to Explain the present invention, be not intended to limit the present invention.
The present invention provides a kind of analysis method of crop climate yield potentiality, comprises the following steps:
Step 1, crop climate Potential Production Analysis model is built, the crop climate Potential Production Analysis model represents For:
Yw=∑ Ywi (1)
Wherein:
Ywi=ri·yi (2)
pi=0.6pi-1-wi+di (4)
Wherein:YwTotal climate yields potentiality in-analysis Period Length;
YwiDay climate yields potentiality in-analysis Period Length;
ri- daily precipitation correction coefficient;
yi- sunlight yield potentiality;
pi- day moisture supply and demand unspent amount;
pi-1Moisture supply and demand-yesterday unspent amount;
wi- day water requirement;
wm- this month first day water requirement;
wm+1- next month first day water requirement;
Day sequence in the n- months, its value is maximum to be no more than 31;
di- the same day actual precipitation;
Step 2, when needing to analyze in cotton development stage total climate yields potentiality corresponding to specified time length, Xiang Bu Start time data and deadline data are inputted in the climate yields Potentials model that rapid 1 structure obtains;
Step 3, after climate yields Potentials model receives start time data and deadline data, by initial time Any one day is denoted as the jth moon i-th day in the period of deadline, then, reads pre-stored of that month first day water requirement wm With first day next month water requirement wm+1Value, and substitute into formula (5), a day water requirement w be calculatedi;Wherein, n=i;
In this step, the value of pre-stored each first day moon water requirement is shown in Table 1:
Each first day moon water requirement of 1 cotton of table
Wherein, since the breeding time of cotton is the 5-9 months, so each first day moon water requirement in 1 storage 5-9 month of table.
Step 4, the same day actual precipitation d is readiWith moisture supply and demand unspent amount p yesterdayi-1, by day water requirement wi, the same day Actual precipitation diWith moisture supply and demand unspent amount p yesterdayi-1Formula (4) is substituted into, a day moisture supply and demand unspent amount p is calculatedi;This Place, the 1st day for being analyzed the period, it was 0 to make moisture supply and demand yesterday unspent amount;
Step 5, by day moisture supply and demand unspent amount piFormula (3) is substituted into, daily precipitation correction coefficient r is calculatedi
Step 6, pre-stored sunlight yield potentiality y is readi, by sunlight yield potentiality yiWith daily precipitation correction coefficient riFormula (2) is substituted into, a day climate yields potentiality Y is calculatedwi
Step 7, repeat step 3- steps 6, are calculated day climate yields potentiality daily in the analyzed period, so Afterwards, based on formula (1), the day climate yields potentiality in each day being calculated are done into summation operation, that is, obtain the analyzed period Interior total climate yields potentiality;
Step 8, total climate yields potentiality in the analyzed period that output step 7 obtains.
For example, working as needs to count a certain regional total climate yields potentiality in this period of 2015.5.1 to 2015.9.1 When, it is only necessary to crop climate Potential Production Analysis mode input initial time it is 2015.5.1, deadline 2015.9.1 The two data, crop climate Potential Production Analysis model can be calculated automatically from total climate yields potentiality in this section. Circular is:
(1) for any one day, such as 2015.8.3, day water requirement is calculated by formula 5 first.
The day water requirement=2015.8.1 water requirements+(2015.9.1 water requirement -2015.8.1 water requirements) of 8.3 this days
From above-mentioned formula as can be seen that when calculating day water requirement using formula 5, for the water day by day of the same moon each day Amount, since the of that month first day, day water requirement gradually increased with function parabolic, and constantly approach the water day by day of first day next month Amount.Using this computational methods, more meet actual demand of the cotton growth stage to water requirement.
(2) then, the day moisture supply and demand unspent amount of this days of 2015.8.3 is calculated, i.e.,:
2015.8.3 day moisture supply and demand unspent amount=0.6 (the day moisture supply and demand surplus of this days of 2015.8.2 of this day Amount) -2015.8.3 day water requirement+2015.8.3 day actual precipitation.
(3) the daily precipitation correction coefficient of this days of 2015.8.3 is calculated using formula 3.
(4) the sunlight yield potentiality of this days of 2015.8.3 and daily precipitation correction coefficient are substituted into formula 2, calculated To the day climate yields potentiality of this days of 2015.8.3.
(5) for this period of 2015.5.1 to 2015.9.1, daily day weather is calculated using the above method Yield potentiality, then sum, that is, obtain total climate yields potentiality of this period.
In the case of same water supply, if using the computational methods of traditional moisture correction factor, it can obtain in Fig. 1 Broken line curve;And if method using the present invention, daily precipitation correction coefficient is parabolic in Fig. 1.Fig. 1 is observed, can be with Find out, conventional method there are 3 shortcomings:Firstth, using megastage precipitation, such as monthly per ten days precipitation, rather than daily precipitation, from And the possibility of real-time tracking calculating yield potentiality is blocked;Secondth, megastage precipitation adds up, equivalent to using mean precipitation Amount, does not reflect the differentiation of service position completely, is not inconsistent with crop Dry Matter Production accumulation environment;3rd, using threshold zone Between divide fixed correction factor, be unified for 0 outside a value section as fixed in section, also do not meet crop Dry Matter Production and accumulate It is actual.And the improved method of the present invention overcomes above three defect completely, it can truly reflect that the change of moisture supply environment is brought Yield potentiality change.
The analysis method of crop climate yield potentiality provided by the invention, innovative point are as follows:
(1) the day water requirement being calculated using function approximation, more meets the reality of the Sino-Japan water requirement of process of crop growth Situation;
(2) approximation computation day moisture supply and demand unspent amount, result of calculation more meet the Sino-Japan moisture supply and demand knot of process of crop growth Remaining actual conditions;
(3) injury caused by approximation computation moisture off-target state, i.e., truly calculate daily precipitation correction coefficient;
(4) total climate yields potentiality are obtained by the accumulation of day climate yields potentiality, can reflects crop climate yield potentiality Forming process, accurately reflect crop dry-matter accumulation state, improve climate yields potentiality accuracy and objectivity.
(5) since daily day climate yields potentiality can be calculated, easy to calculate the automation applications such as monitoring in real time daily;
(6) it is suitable for the fields such as crop yield monitoring, agricultural production resources assessment, Resources Evolution and early warning.
In addition, the present invention has following two big innovations:
Innovation one:Sunlight yield potentiality
In above-mentioned steps 6, the sunlight yield potentiality y that is readiObtain by the following method:
Step 6.1, pre-stored three base point temperature of the analyzed moon is read, is respectively preference temperature of the cotton in current period T1, the minimum temperature T that can bear2And the maximum temperature T that can be born3;A day actual average temperature T is read againi
Step 6.2, by T1、T2、T3And TiFormula (6) is substituted into, degree/day correction factor F is calculatedi
Wherein, B=(T3-T1)/(T1-T2);
Step 6.3, pre-stored day highest real air temperature T is readu, by day highest real air temperature TuSubstitute into formula (7), meter Calculation obtains daily maximum temperature correction factor Vi
Vi=(1-li2)25-29li 2
Wherein:
Step 6.4, by daily maximum temperature correction factor ViWith degree/day correction factor FiFormula (8) is substituted into, light is calculated Warm yield potentiality day correction factor ti
ti=0.3Vi+0.7Fi (8)
Step 6.5, photosynthetic yield potentiality h of pre-stored day is readi, by tiAnd hiFormula (9) is substituted into, daylight is calculated Warm yield potentiality yi
yi=ti·hi (9)。
Innovation two:Day photosynthetic yield potentiality
In step 6.5, the day photosynthetic yield potentiality h that is readiObtain by the following method:
Step 6.5.1, establishes each first day moon weight coefficient table in breeding time;
Step 6.5.1, each first day moon weight coefficient table established are shown in Table 2:
The 2 cotton photosynthetic potentiality moon of table first day weight
Wherein, v4、v5、v6、v7、v8、v9And v10The respectively first day in April weight, the first day in May weight, the first day in June power Weight, the first day in July weight, August part first day weight, September part first day weight and the first day in October weight.
Step 6.5.2, reads each first day moon weight coefficient table, obtains of that month first day weight vmWith first day next month weight vm+1, then by of that month first day weight vmWith first day next month weight vm+1Formula (10) is substituted into, a day weight v is calculatedi
Wherein, n is day sequence, n=i in the moon;
Step 6.5.3, reads pre-stored day radiation G, day is radiated G and day weight viFormula (11) is substituted into, is calculated To day photosynthetic yield potentiality hi
hi=vi·G (11)。
Test example 1
For Huang-Huai-Hai, the climatological data day by day of 51 website 1971-2006 is counted, calculates each website respectively Light temperature yield potentiality, photosynthetic yield potentiality and weather yield potentiality, then extended according to inverse distance and neighbouring principle of similarity Calculate light temperature yield potentiality, photosynthetic yield potentiality and the weather yield potentiality in remaining 383 counties and cities level administrative area.Result of calculation is used ArcMap maps, and forms corresponding zoning map.The result is shown in Fig. 2, Fig. 3 and Fig. 4, is respectively Cotton climatic output potentiality zoning map, cotton Use up warm yield potentiality zoning map, cotton photosynthetic yield potentiality zoning map.
Observe Fig. 2, Fig. 3 and Fig. 4, it can be seen that climate yields potentiality minimum value 1506kghm in area-2, maximum 1874kg·hm-2, difference is little.The general trend of climate yields potential value be in addition to Jiaodong Peninsula from southwest northeastward direction by Step increase, it is close with the variation tendency of light temperature yield potentiality zoning map and photosynthetic yield potentiality zoning map, and with actual gas It is consistent to wait change of production trend, so as to demonstrate the accuracy of climate yields Potentials method provided by the invention.
Test example 2
According to 5 years before and after in Huang-Huai-Hai 1971-2006 terms meteorological datas day by day, 1971-1975 was calculated respectively With climate yields potentiality in two periods of 2002-2006, the differentiation of the climate yields potentiality of the two periods of comparative analysis Situation, as shown in figure 5, the variation diagram for 1971-1975 Cotton climatic output potentiality;As shown in fig. 6, it is 2002-2006 The variation diagram of Cotton climatic output potentiality.
Comparison diagram 5 and Fig. 6 can see, the northern direction transfer eastwards of Huang-Huai-Hai highest climate yields REGION OF WATER INJECTION OILFIELD.Most 2042kghm of the good climate yields potential value from 1971-1975-2Drop to the 1930kghm of 2002-2006-2, together When, 1479kghm of the minimum climate yields potentiality from 1971-1975-2Drop to the 1328kghm of 2002-2006-2, Decline 5.5% and 10.2% respectively.Wherein terrestrial surface radiation maximum declines 8.6%, and minimum value declines 12.8%;Photosynthetic yield is dived Power maximum declines 7.6%, and minimum value declines 17.2%;Light temperature yield potentiality maximum declines 10.4%, and minimum value declines 13.4%.In terms of the region change conditions of yield potentiality, it is basically identical with actual Cotton Production layout evolving trend.
The research of Huang-Huai-Hai Cotton climatic output potentiality, can not only effectively comment regional Climate-Ecology Suitability Valency, and can deeply disclose from different places climatic factor to the influence degree of production so that development to region Cotton Production with And prospective evaluation is made in the differentiation of production distribution, National agricultural and Cotton Production are acted on important guiding.
The above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications also should Depending on protection scope of the present invention.

Claims (5)

1. a kind of analysis method of crop climate yield potentiality, it is characterised in that comprise the following steps:
Step 1, crop climate Potential Production Analysis model is built, the crop climate Potential Production Analysis model is expressed as:
Yw=∑ Ywi (1)
Wherein:
Ywi=ri·yi (2)
<mrow> <msub> <mi>r</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>1</mn> <mo>-</mo> <msubsup> <mi>log</mi> <mn>25</mn> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mo>|</mo> <msub> <mi>p</mi> <mi>i</mi> </msub> <mo>|</mo> <mo>)</mo> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
pi=0.6pi-1-wi+di (4)
<mrow> <msub> <mi>w</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>w</mi> <mi>m</mi> </msub> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>w</mi> <mrow> <mi>m</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>w</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msubsup> <mi>log</mi> <mn>31</mn> <mi>n</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
Wherein:YwTotal climate yields potentiality in-analysis Period Length;
YwiDay climate yields potentiality in-analysis Period Length;
ri- daily precipitation correction coefficient;
yi- sunlight yield potentiality;
pi- day moisture supply and demand unspent amount;
pi-1Moisture supply and demand-yesterday unspent amount;
wi- day water requirement;
wm- this month first day water requirement;
wm+1- next month first day water requirement;
Day sequence in the n- months, its value is maximum to be no more than 31;
di- the same day actual precipitation;
Step 2, when needing to analyze in cotton development stage total climate yields potentiality corresponding to specified time length, to step 1 Build and start time data and deadline data are inputted in obtained climate yields Potentials model;
Step 3, after climate yields Potentials model receives start time data and deadline data, by initial time to section Only any one day is denoted as the jth moon i-th day in the period of time, then, reads pre-stored of that month first day water requirement wmWith under The first day water requirement w moonm+1Value, and substitute into formula (5), a day water requirement w be calculatedi;Wherein, n=i;
Step 4, the same day actual precipitation d is readiWith moisture supply and demand unspent amount p yesterdayi-1, by day water requirement wi, the same day it is actual Precipitation diWith moisture supply and demand unspent amount p yesterdayi-1Formula (4) is substituted into, a day moisture supply and demand unspent amount p is calculatedi;Herein, it is right In the 1st day of the analyzed period, it was 0 to make moisture supply and demand yesterday unspent amount;
Step 5, by day moisture supply and demand unspent amount piFormula (3) is substituted into, daily precipitation correction coefficient r is calculatedi
Step 6, pre-stored sunlight yield potentiality y is readi, by sunlight yield potentiality yiWith daily precipitation correction coefficient riGeneration Enter formula (2), a day climate yields potentiality Y is calculatedwi
Step 7, repeat step 3- steps 6, are calculated day climate yields potentiality daily in the analyzed period, then, base In formula (1), the day climate yields potentiality in each day being calculated are done into summation operation, that is, are obtained total in the analyzed period Climate yields potentiality;
Step 8, total climate yields potentiality in the analyzed period that output step 7 obtains.
2. the analysis method of crop climate yield potentiality according to claim 1, it is characterised in that in step 3, be pre-stored The value of each first day moon water requirement be shown in Table 1:
Each first day moon water requirement of 1 cotton of table
Wherein, since the breeding time of cotton is the 5-9 months, so each first day moon water requirement in 1 storage 5-9 month of table.
3. the analysis method of crop climate yield potentiality according to claim 1, it is characterised in that in step 6, read The sunlight yield potentiality y arrivediObtain by the following method:
Step 6.1, pre-stored three base point temperature of the analyzed moon is read, is respectively preference temperature T of the cotton in current period1, energy The minimum temperature T enough born2And the maximum temperature T that can be born3;A day actual average temperature T is read againi
Step 6.2, by T1、T2、T3And TiFormula (6) is substituted into, degree/day correction factor F is calculatedi
<mrow> <msub> <mi>F</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>)</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>B</mi> </msup> </mrow> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>)</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mi>B</mi> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
Wherein, B=(T3-T1)/(T1-T2);
Step 6.3, pre-stored day highest real air temperature T is readu, by day highest real air temperature TuFormula (7) is substituted into, is calculated To daily maximum temperature correction factor Vi
Vi=(1-li 2)25-29li 2
Wherein:
<mrow> <msub> <mi>l</mi> <mi>i</mi> </msub> <mo>=</mo> <msubsup> <mi>log</mi> <mn>26</mn> <msub> <mi>T</mi> <mi>u</mi> </msub> </msubsup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msubsup> <mi>log</mi> <mn>26</mn> <msub> <mi>T</mi> <mi>u</mi> </msub> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
Step 6.4, by daily maximum temperature correction factor ViWith degree/day correction factor FiFormula (8) is substituted into, the production of light temperature is calculated Measure potentiality day correction factor ti
ti=0.3Vi+0.7Fi (8)
Step 6.5, photosynthetic yield potentiality h of pre-stored day is readi, by tiAnd hiFormula (9) is substituted into, sunlight production is calculated Measure potentiality yi
yi=ti·hi (9)。
4. the analysis method of crop climate yield potentiality according to claim 3, it is characterised in that in step 6.5, read The day photosynthetic yield potentiality h gotiObtain by the following method:
Step 6.5.1, establishes each first day moon weight coefficient table in breeding time;
Step 6.5.2, reads each first day moon weight coefficient table, obtains of that month first day weight vmWith first day next month weight vm+1, Then by of that month first day weight vmWith first day next month weight vm+1Formula (10) is substituted into, a day weight v is calculatedi
<mrow> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>v</mi> <mi>m</mi> </msub> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mrow> <mi>m</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msubsup> <mi>log</mi> <mn>31</mn> <mi>n</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
Wherein, n is day sequence, n=i in the moon;
Step 6.5.3, reads pre-stored day radiation G, day is radiated G and day weight viFormula (11) is substituted into, daylight is calculated Close yield potentiality hi
hi=vi·G (11)。
5. the analysis method of crop climate yield potentiality according to claim 4, it is characterised in that step 6.5.1, is built Vertical each first day moon weight coefficient table is shown in Table 2:
The 2 cotton photosynthetic potentiality moon of table first day weight
Wherein, v4、v5、v6、v7、v8、v9And v10The respectively first day in April weight, the first day in May weight, the first day in June weight, 7 The first day in month weight, August part first day weight, September part first day weight and the first day in October weight.
CN201510867342.8A 2015-12-01 2015-12-01 The analysis method of crop climate yield potentiality Expired - Fee Related CN105447317B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510867342.8A CN105447317B (en) 2015-12-01 2015-12-01 The analysis method of crop climate yield potentiality

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510867342.8A CN105447317B (en) 2015-12-01 2015-12-01 The analysis method of crop climate yield potentiality

Publications (2)

Publication Number Publication Date
CN105447317A CN105447317A (en) 2016-03-30
CN105447317B true CN105447317B (en) 2018-04-20

Family

ID=55557487

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510867342.8A Expired - Fee Related CN105447317B (en) 2015-12-01 2015-12-01 The analysis method of crop climate yield potentiality

Country Status (1)

Country Link
CN (1) CN105447317B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106651619A (en) * 2016-12-30 2017-05-10 贵州大学 Division method of Chinese yamgrowth period
CN106651618A (en) * 2016-12-30 2017-05-10 贵州大学 Division method of polygonummultiflorumThunb.growthperiod
CN113112165A (en) * 2021-04-20 2021-07-13 武荣盛 Analysis method for crop climate yield potential

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1631098A (en) * 2003-12-25 2005-06-29 中国农业大学 Insufficient irrigation forecast and control method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8650069B2 (en) * 2009-12-03 2014-02-11 Telsco Industries, Inc. Irrigation water management system

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1631098A (en) * 2003-12-25 2005-06-29 中国农业大学 Insufficient irrigation forecast and control method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
不同土壤条件下节水灌溉对棉花产量和灌溉水生产力的影响;苏永中;《土壤学报》;20141130;第51卷(第6期);第1192-1201页 *
不同时期非充分滴灌对棉花光合生产及产量、品质影响的研究;范志超;《中国优秀硕士学位论文全文数据库 农业科技辑》;20140115(第1期);第D047-199页 *

Also Published As

Publication number Publication date
CN105447317A (en) 2016-03-30

Similar Documents

Publication Publication Date Title
Vazifedoust et al. Increasing water productivity of irrigated crops under limited water supply at field scale
Rezvani et al. IoT-based sensor data fusion for determining optimality degrees of microclimate parameters in commercial greenhouse production of tomato
Vera et al. Irrigation factor approach based on soil water content: A nectarine orchard case study
CN105447317B (en) The analysis method of crop climate yield potentiality
Lopez-Diaz et al. The effect of different levels of shading in a photovoltaic greenhouse with a north–south orientation
Mirás-Avalos et al. Irrigation-Advisor—A decision support system for irrigation of vegetable crops
Vujadinović Mandić et al. Observed changes in climate conditions and weather-related risks in fruit and grape production in Serbia
Zhai et al. Rice irrigation schedule optimization based on the AquaCrop model: study of the Longtouqiao irrigation district
Liu et al. Microclimate characteristics and evapotranspiration estimates of cucumber plants in a newly developed sunken solar greenhouse
Blanco et al. Research trends on greenhouse engineering using a science mapping approach
Roxani et al. Multidimensional role of Agrovoltaics in era of EU green Deal: Current status and analysis of water–energy–food–land dependencies
Begum et al. Environmental and social dynamics of urban rooftop agriculture (Urta) and their impacts on microclimate change
Li et al. Simulation of soil water and heat flow under plastic mulching and different ridge patterns
Yohanani et al. Estimating evapotranspiration of screenhouse banana plantations using artificial neural network and multiple linear regression models
Alonso-Marroquin et al. Synergy between photovoltaic panels and green roofs
Carreño-Ortega et al. Lettuce production under mini-PV modules arranged in patterned designs
Shahrokhnia et al. Evaluation of wheat and maize evapotranspiration determination by direct use of the Penman–Monteith equation in a semi-arid region
CN105512947B (en) Make the analysis method of object light temperature yield potentiality
Sun et al. Effects of different irrigation modes on the growth, physiology, farmland microclimate characteristics, and yield of cotton in an oasis
CN105373670B (en) The analysis method of crop photosynthesis yield potentiality
Yang et al. Tomato evapotranspiration, crop coefficient and irrigation water use efficiency in the winter period in a sunken chinese solar greenhouse
CN116578031A (en) Intelligent agricultural greenhouse remote monitoring regulation and control system
Gao et al. Effect of different water treatments in soil-plant-atmosphere continuum based on intelligent weighing systems
Fierro-Chacón et al. Fuzzy logic that determines sky conditions as a key component of a smart irrigation system
Liu et al. Photosynthesis Characteristics of Tomato Plants and Its’ Responses to Microclimate in New Solar Greenhouse in North China

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180420

Termination date: 20211201