Способ прогнозирования урожайности яровой пшеницы

Изобретение относится к области сельского хозяйства. В способе используют множественную регрессионную зависимость Y=a+b1X+b2Z, ц/га, где Х - температуры, Z - осадки, a, b1 и b2 - коэффициенты, при этом измеряют средние активные температуры, средние осадки и соответствующие урожайности М предшествующих вегетационных периодов, на основании которых вычисляют коэффициенты регрессионного уравнения: b1, b2 и а с помощью уравнений:

, где Y, ц/га - урожай яровой пшеницы i-го вегетационного периода, , ц/га - усредненный урожай яровой пшеницы за М лет предшествующих прогнозируемому вегетационному периоду, X, °С - суммарная активная температура i-го вегетационного периода, , °С - усредненная суммарная активная температура за М лет предшествующих прогнозируемому вегетационному периоду, Z, мм - суммарные осадки i-го вегетационного периода, , мм - усредненные суммарные осадки за М лет предшествующих прогнозируемому вегетационному периоду, измеряют среднесуточную активную температуру и соответствующие им осадки в прогнозируемом вегетационном периоде с начала вегетации до даты прогноза, на основании которых вычисляют эквивалентные активные температуры - Хэ и эквивалентные осадки - Zэ для множественного регрессионного уравнения по формулам:

где , дн - усредненная длительность вегетации, измеренная и вычисленная из предшествующих вегетационных периодов; ta , °C - суточная активная температура прогнозируемого вегетационного периода; oа, мм - соответствующие суточные осадки прогнозируемого вегетационного периода; и n - текущее количество дней в прогнозируемом вегетационном периоде на момент прогноза, которые подставляют в множественное регрессионное уравнение вместо средних температур и осадков. Способ позволяет уменьшить трудоемкость прогнозирования урожайности яровой пшеницы и получать оперативную и достоверную информацию о количестве урожая. 5 табл.

Изобретение относится к сельскому хозяйству и может быть использовано при прогнозировании урожая мягкой яровой пшеницы.

Известен способ оценки потенциальной продуктивности сельскохозяйственных растений, преимущественно озимых зерновых колосовых культур, при возделывании в условиях резко континентального климата. Прогнозируемую урожайность озимой пшеницы определяют по формуле

Y=-2,39+0,112N+0,07Ом-0,0008N2 -0,0005 Oм2,

R=0,82,

где Y - урожайность озимой пшеницы, т/га; N - доза азотных удобрений от 31 до 93 кг/га действующего вещества; Ом - количество осадков за май от 18 до 104 мм; R - коэффициент множественной корреляции (Рогачев А.Ф. (RU); Саддаев A.M. (RU); Рогачев Д.A. (RU). - Заявка 2002126981/12; заявлено 09.10. 2002; опубликовано: 20.07.2006, бюл. 2006.20).

К недостаткам описанного способа прогнозирования урожайности зерновых колосовых относится его весьма узкая применимость из-за резкого изменения числовых коэффициентов уравнения, которые сильно зависят: от изменения хозяйства возделывания растений, изменения культуры или даже сорта сельскохозяйственных растений, и наконец, при изменении тенденций в климате. Кроме того, способ применим лишь в условиях острозасушливого климата и имеет недостаточную репрезентативность при использовании в почвенно-климатических условиях Западной Сибири.

Известен способ прогнозирования урожайности озимой пшеницы с использованием математической зависимости

Y=75,6-3,14X+12,52D,

где Y - прогнозируемая урожайность, ц/га; X - среднесуточная температура в мае, °C; D - доза минеральных удобрений от 0 до 1 (0 - контроль, 1 - N100P100K 50). (Акулов П.Г., Понедельченко М.Н., Сокорева Н.Н., Сокорев Н.С. RU 2158498 C2, 7 A01G 7/00, 10.11.2000).

Как и в первом аналоге, к основному недостатку описанного способа прогнозирования урожайности зерновых колосовых относится его весьма узкая применимость из-за резкого изменения числовых коэффициентов уравнения, которые сильно меняют свое значение: при изменении района возделывания растений, изменении культуры или даже сорта сельскохозяйственных растений, и наконец, при изменении тенденций в климате. Недостаток этого способа состоит в том, что он применим только для прогнозирования урожайности пшеницы в полевых опытах (при высокой культуре земледелия) с использованием удобрений при использовании дозы N100P100K50 . Кроме того, к недостаткам описанного способа прогнозирования урожайности зерновых колосовых относятся низкая достоверность прогнозируемых данных и малый их срок службы.

Наиболее близким аналогом прогнозирования урожайности мягкой яровой пшеницы является комплексный способ прогноза урожайности яровой пшеницы в Западной Сибири (Костюков В.В., Костюкова Н.И., Черникова М.И.; опубл. в журнале «География и природные ресурсы». - 2003. - 1. - С.153-155).

Способ состоит в использовании комбинации из детерминированной и физико-статистической моделей. В основу детерминированной модели положено представление динамики урожайности в виде суммы гармоник, выявленных заранее в результате статистических расчетов, а также линейного тренда, описывающего многолетнюю культуру земледелия (сорт, технология возделывания, удобрения и т.п.):

где ai - амплитуда; i - сдвиг по фазе; Li - период; m - число гармоник; (b n+d) - тренд; n - номер года, начиная с 1956, исходного в рабочем архиве.

Общий вид единой физико-статистической модели для прогноза урожайности основных зерновых культур (в амбарном весе) может быть представлен в виде суммы этих моделей и выражена, как указано в таблице:

Примечание:

У - урожайность, ц/га (амбарный вес);

TVI1 - средняя суточная температура воздуха за первую декаду июня;

TVI2 - средняя суточная температура воздуха за вторую декаду июня;

TVII1,2 - средняя суточная температура воздуха за период с первой декады по вторую декаду июля;

OCVI1 - сумма осадков за первую декаду июня;

OCVI2 - сумма осадков за вторую декаду июня;

dVI1,2 - среднесуточный дефицит насыщения воздуха за период с первой декады по вторую декаду июня;

n - порядковый номер года, начиная с 1956 (1956 принят за 1).

ГТК VI1,2 - гидротермический коэффициент Селянинова с первой по вторую декады июня;

ОС VII1,2 - сумма осадков за период с первой по вторую декады июля;

dVII1,2 - среднесуточный дефицит насыщения воздуха за период с первой декады по вторую декаду июля;

ГТК VII1,2 - гидротермический коэффициент Селянинова с первой по вторую декады июля;

ГТК (VI2-VI1,2) - гидротермический коэффициент Селянинова со второй декады июня по вторую декаду июля, которые вычисляются с помощью уравнений:

где OCVI1,2, OCVII1,2 ; OC(VI2-VII2); ; ; - сумма осадков и температур за указанный в формуле период.

Указанный способ позволяет прогнозировать урожайность в течение вегетационного периода.

К недостаткам описанного способа, принятого нами в качестве наиближайшего аналога, относятся, несмотря на то, что коэффициенты уравнений прогноза урожайности определяются для каждой зоны возделывания сельскохозяйственных культур, культуры и даже сорта, а также учитываются тенденции в изменении климата, высокая сложность способа, потребность в дорогостоящей аппаратуре, а также большие усилия в организационной работе для выполнения прогноза.

Задачей настоящего изобретения является уменьшение трудоемкости прогнозирования урожайности яровой пшеницы и повышение точности прогноза.

Эффект и технический результат, которые могут быть получены при реализации заявленного способа, оцениваются получением оперативной и достоверной информации о количестве урожая, планировании уборочной компании, а также созданием точной стратегии по созданию товарного зерна и запасов.

Это достигается тем, что прогнозирование урожайности яровой пшеницы в течение вегетационного периода, ведут с помощью множественной линейной регрессионной зависимости, полученной на основании массива данных, состоящего из данных об урожайности - Y, ц/га хозяйства за М лет; данных от гидрометеостанции в виде суммы активных температур - X, °C и суммы осадков - Z, мм, прогнозируемой природно-климатической зоны также за эти же M лет в виде:

Y=a+b1 X+b2Z, ц/га,

где a - начало отсчета уравнения урожайности, b1 и b2 - коэффициенты частной регрессии. Параметры a, b1, b2 вычисляются методом, наименьших квадратов из которого получают соотношения:

Для прогноза также вычисляют эквиваленты сумм активных температур - Xэ, °С и эквиваленты сумм осадков - Zэ, мм с помощью выражений:

где инструментально измеряют: ta , °C - активные суточные температуры; oa, мм - суточные осадки; n - количество дней, прошедших с начала вегетационного периода до даты прогноза и вычисляют - усредненную длительность вегетации, взятую из массива данных предшествующих урожаев для конкретного сорта яровой пшеницы и хозяйства.

Вычисляют будущий урожай яровой пшеницы, прогнозируемого вегетационного периода, путем введения измеренных и вычисленных данных в регрессионное уравнение: a, b1 , b2, Xэ, Zэ.

Y=a+b1Xэ+b2Zэ

Способ прогноза урожайности реализуется следующим образом.

Для построения множественной линейной регрессионной зависимости урожайности от суммы активных температур и суммы осадков заказывают данные о суммах активных температур и суммах осадков в гидрометеостанции за десятилетний или пятилетний периоды, предшествующие прогнозируемому вегетационному периоду. Длительность периода (M=10÷5 лет) определяется из тренда многолетнего изменения урожайности. Также получают у хозяйства за эти же М лет данные об урожайности прогнозируемого сорта яровой пшеницы. Эти данные группируют в таблицу (приведены ниже в примерах). Проведя анализ таблиц вычисляют параметры a, b1, b 2, которые подставляем в множественную линейную регрессионную зависимость и тем самым получают рабочее уравнение для прогноза урожайности выбранного сорта яровой пшеницы в выбранном хозяйстве с учетом климатических тенденций.

Далее с момента начала сева измеряют активные суточные температуры и суточные осадки, которые суммируют и делят на количество прошедших дней с момента сева. Полученные значения умножают на среднее значение вегетационного периода для данного сорта яровой пшеницы. Эти величины подставляют в регрессионное уравнение и вычисляют будущую урожайность на конец вегетационного периода. По прогнозным данным определяют: количество уборочной техники; потребность в хранилищах; потребность во вспомогательных механизмах; количество товарного зерна и соответствующую выручку; количество зерна, засыпаемого в запасы; и другие показатели хозяйства, позволяющие оптимизировать прибыли и убытки.

Реализацию заявленного способа прогнозирования урожайности, проиллюстрируем на четырех примерах возделывания яровой пшеницы в разных природно-климатических зонах юго-востока Западной Сибири. Для сопоставимости приводим урожаи на полях государственных сортоиспытательных участков (ГСУ) и производственных предприятий. Эти предприятия и ГСУ расположены в степной и лесостепной зонах (Кемеровской области). Они имеют разную культуру земледелия. Прогнозирование урожайности осуществлялось по данным многолетней урожайности хозяйств в период посева и уборки в 1989 по 1998 годах, а также по данным об активных температурах и осадках гидрометеостанции приведенных районов.

Пример 1. Работа выполнялась для Ленинского ГСУ, в котором измерялся многолетний массив данных (урожайность яровой пшеницы и данные гидрометеостанции о многолетних температурах и осадках в этой местности). Причем данные о температурах и осадках гидрометеостанции трансформировали в суммарные активные температуры и суммарные осадки рассматриваемого вегетационного периода, а урожайность измерялась как амбарная урожайность. Эти данные после измерений сбора и обработки были сгруппированы в таблицу 1.

Множественная линейная регрессия для данного хозяйства, вычисленная на основании таблицы 1 данных, имеет вид:

Y=-1765,99+1,48447*X э-2,2978*Zэ,

это уравнение позволило сделать прогноз урожайности для Ленинского ГСУ на конец вегетационного периода 1998 года. Прогноз выполняли на основании инструментально измеренных суточных активных температур и суточных осадков до 25 июня и 25 июля 1998 года этого вегетационного периода. Результаты прогноза размещены в таблице 1 в соответствующих строках.

Пример 2. Работа выполнялась в степной зоне для производственного предприятия, расположенного в зоне действия Ленинского ГСУ. Поэтому эквивалентные активные температуры и эквивалентные осадки подсчитывались на одном исходном материале. Измерялся многолетний массив данных (урожайность яровой пшеницы для этого хозяйства и данные гидрометеостанции о многолетних температурах и осадках в этой местности). Причем данные о температурах и осадках гидрометеостанции трансформировали в суммарные активные температуры и суммарные осадки рассматриваемого вегетационного периода, а урожайность измерялась как амбарная урожайность. Эти данные после измерений сбора и обработки были сгруппированы в таблицу 2.

Множественная линейная регрессия для этого хозяйства, вычисленная на основании таблицы 2 данных, имеет вид:

Y=725,414-0,57386*Xэ+0,79589*Z э.

Это уравнение позволило сделать прогноз урожайности на конец вегетационного периода 1998 года. Прогноз выполняли на основании инструментально измеренных суточных активных температур и суточных осадков до 25 июня и 25 июля 1998 года этого вегетационного периода. Результаты прогноза размещены в таблице 2 в соответствующей строке.

Пример 3. Работа выполнялась для Мариинского ГСУ, в котором измерялся многолетний массив данных (урожайность яровой пшеницы и данные гидрометеостанции о многолетних температурах и осадках в этой местности). Причем данные о температурах и осадках гидрометеостанции трансформировали в суммарные активные температуры и суммарные осадки рассматриваемого вегетационного периода, а урожайность измерялась как амбарная урожайность ГСУ. Эти данные после измерений сбора и обработки были сгруппированы в таблицу 3.

Множественная линейная регрессия для Мариинского ГСУ, вычисленная на основании таблицы 3, имеет вид:

Y=1413,85-1,2399*Xэ+2,10869*Z э.

Это уравнение позволило сделать прогноз урожайности для Мариинского ГСУ на конец вегетационного периода 1998 года. Прогноз выполняли на основании инструментально измеренных суточных активных температур и суточных осадков до 25 июня и 25 июля 1998 года этого вегетационного периода. Результаты прогноза размещены в таблице 3 в соответствующих строках.

Пример 4. Работа выполнялась для предприятия, расположенного в зоне действия Мариинского ГСУ, в котором измерялся многолетний массив данных (урожайность яровой пшеницы и данные гидрометеостанции о многолетних температурах и осадках в этой местности). Причем данные о температурах и осадках гидрометеостанции трансформировали в суммарные активные температуры и суммарные осадки рассматриваемого вегетационного периода, а урожайность измерялась как амбарная урожайность. Эти данные после измерений сбора и обработки были сгруппированы в таблицу 4.

Множественная линейная регрессия для данного хозяйства, вычисленная на основании таблицы 4 данных, имеет вид:

Y=1265,41-1,1184*Xэ+1,8869*Zэ.

Это уравнение позволило сделать прогноз урожайности на конец вегетационного периода 1998 года. Прогноз выполняли на основании инструментально измеренных суточных активных температур и суточных осадков до 25 июня и 25 июля 1998 года этого вегетационного периода. Результаты прогноза размещены в таблице 4 в соответствующей строке.

Анализ прогнозов показывает, что при приближении к окончанию прогнозируемого вегетационного периода возрастает точность прогноза. Неточность прогноза может быть обусловлена потерями урожая при его уборке. Причем предложенный способ прогноза позволяет точно планировать все хозяйственные мероприятия без привлечения избыточных ресурсов и оптимизирует прибыль хозяйства.

БИБЛИОГРАФИЧЕСКИЕ ДАННЫЕ

1. Рогачев А.Ф. (RU); Салдаев A.M. (RU); Рогачев Д.A. (RU). - Заявка 2002126981/12; заявлено 09.10.2002; опубликовано: 20.07.2006, бюл. 200620.

2. Акулов П.Г., Понедельченко М.Н., Сокорева Н.Н., Сокорев Н.С. RU 2158498 C2, 7 A01G 7/00, 10.11.2000.

3. Костюков В.В., Костюкова Н.И., Черникова М.И. Комплексный способ прогноза урожайности яровой пшеницы в Западной Сибири. // Журнал «География и природные ресурсы», 2003. 1. - С.153-155.

Формула изобретения

Способ прогнозирования урожая яровой пшеницы на основании множественной регрессионной зависимости Y=a+b1X+b 2Z, ц/га, где X - температуры, Z - осадки, a, b1 и b2 - коэффициенты, отличающийся тем, что измеряют средние активные температуры, средние осадки и соответствующие урожайности М предшествующих вегетационных периодов, на основании которых вычисляют коэффициенты регрессионного уравнения: b 1, b2 и а с помощью уравнений: ,где Y, ц/га - урожай яровой пшеницы i-го вегетационного периода, , ц/га - усредненный урожай яровой пшеницы за М лет, предшествующих прогнозируемому вегетационному периоду, X, °С - суммарная активная температура i-го вегетационного периода, , °С - усредненная суммарная активная температура за М лет, предшествующих прогнозируемому вегетационному периоду, Z, мм - суммарные осадки i-го вегетационного периода, , мм - усредненные суммарные осадки за М лет, предшествующих прогнозируемому вегетационному периоду, измеряют среднесуточную активную температуру и соответствующие им осадки в прогнозируемом вегетационном периоде с начала вегетации до даты прогноза, на основании которых вычисляют эквивалентные активные температуры - Хэ и эквивалентные осадки - Zэ для множественного регрессионного уравнения по формулам: ; ,где, , дн. - усредненная длительность вегетации, измеренная и вычисленная из предшествующих вегетационных периодов; t a, °C - суточная активная температура прогнозируемого вегетационного периода; oа, мм - соответствующие суточные осадки прогнозируемого вегетационного периода; и n - текущее количество дней в прогнозируемом вегетационном периоде на момент прогноза, которые подставляют в множественное регрессионное уравнение вместо средних температур и осадков.

MM4A Досрочное прекращение действия патента из-за неуплаты в установленный срок пошлины за поддержание патента в силе

Дата прекращения действия патента: 24.03.2012

Дата публикации: 20.01.2013