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

Изобретение относится к области сельского хозяйства. Способ включает многолетние ежегодные наблюдения за фактическими значениями урожайности и преобразование фактических значений урожайности в качественные величины. При этом фактические значения урожайности преобразуют в ряд лингвистических значений «Высокая урожайность», «Средняя урожайность», «Низкая урожайность». Далее выделяют фрагменты предпрогнозного периода преобразованного ряда продолжительностью 1, 2 и 3 года. Прогнозируемую урожайность на следующий год находят с учетом повторяемости фрагментов предпрогнозного периода в преобразованном ряду значений урожайности за все годы наблюдений и поправочного коэффициента, учитывающего относительную повторяемость фрагментов предпрогнозного периода, по формуле , где i - номер фрагмента предпрогнозного периода, i=1, 2, 3; yi - прогнозная урожайность на следующий год, определяемая с учетом повторяемости i-го фрагмента предпрогнозного периода в преобразованном ряду значений урожайности за все годы наблюдений; ki - поправочный коэффициент, учитывающий относительную повторяемость i-го фрагмента предпрогнозного периода. Изобретение позволяет повысить достоверность прогнозирования урожайности. 5 табл.

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

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

у=-21,12+0,33х1-8·10-4х1 2+0,023х2-8·10-6х 2 2,

у=-13,0+0,2х1-4·10 -4х1 2+0,02x2-7·10-6х 2 2,

где у - урожайность зерна озимой пшеницы, ц/га;

x1 - запасы продуктивной влаги в метровом слое почвы;

х2 - число стеблей на 1 м2.

(Уланова Е.С. "Методические указания по составлению долгосрочных прогнозов урожайности озимой пшеницы в основных районах ее возделывания". - Л.: Гидрометеоиздат, 1997).

Недостатком их является высокая трудоемкость (бурение почвы, определение запасов продуктивной влаги, подсчет числа стеблей на 1 м2) и низкая достоверность.

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

у=75,6-3,14х+12,52d,

где у - урожайность озимой пшеницы, ц/га;

x - среднесуточная температура воздуха в мае, °С;

d - доза минеральных удобрений от 0 до 1 (0 - без удобрений, 1 - N120P120 K60)

(патент RU 2158498 С2, МПК 7 A01G 7/00. Способ прогнозирования урожайности озимой пшеницы. / Акулов П.Г., Понедельченко М.Н., Сокорева И.Н., Сокорев Н.С. Заявл. 30.11.1998; опубл. 10.11.2000).

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

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

где max - максимальное значение в множестве С,

M(C ) - среднее число элементов для а-уровневого множества С ,

F(C) - точечное значение удовлетворительности для множества С.

(Андрейчиков А.В. Анализ, синтез, планирование решений в экономике. / А.В.Андрейчиков, О.Н.Андрейчикова. - М.: Финансы и статистика, 2000. - С.157-206, гл. «Методы принятия решений на основе теории нечетких множеств»).

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

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

,

где yi, xi - урожайность озимой пшеницы в индицирующих регионах в первом предшествующем прогнозируемому году;

max у, max x - максимальная урожайность озимой пшеницы в соответствующем индицирующем регионе;

a, b - эмпирические коэффициенты, характеризующие связи колебаний природных условий прогнозируемого региона с индицирующим (температура, осадки, гидротермические показатели и т.п.), а [0,2; 0,4], b [0,3; 0,8], прогнозируют подъем или спад урожайности, если соответственно sign( )=1 или sign( )=-1.

(заявка RU 2005114095 А. МПК A01G 7/00 (2006.01). Способ прогнозирования межгодовых колебаний урожайности озимой пшеницы. / И.Б.Загайтов (RU), Л.П.Яновский (RU). - Заявлено 11.05.2005; опубл. 20.11.2006).

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

Сущность заявленного изобретения заключается в следующем.

Задача, на решение которой направлено заявленное изобретение, - повышение достоверности прогнозирования.

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

Указанный технический результат достигается тем, что в способе прогнозирования урожайности озимой пшеницы, возделываемой в условиях засушливого климата, включающем многолетние ежегодные наблюдения за фактическими значениями урожайности и преобразование фактических значений урожайности к качественным величинам, согласно изобретению фактические значения урожайности преобразуют в ряд лингвистических значений «Высокая урожайность», «Средняя урожайность», «Низкая урожайность», выделяют фрагменты предпрогнозного периода преобразованного ряда продолжительность 1, 2 и 3 года, а прогнозируемую урожайность на следующий год находят с учетом повторяемости фрагментов предпрогнозного периода в преобразованном ряду значений урожайности за все годы наблюдений и поправочного коэффициента, учитывающего относительную повторяемость фрагментов предпрогнозного периода, по формуле

где i - номер фрагмента предпрогнозного периода, i=1, 2, 3;

yi - прогнозная урожайность на следующий год, определяемая с учетом повторяемости i-го фрагмента предпрогнозного периода в преобразованном ряду значений урожайности за все годы наблюдений;

k i - поправочный коэффициент, учитывающий относительную повторяемость i-го фрагмента предпрогнозного периода, определяемый по формуле

где i, j - номер фрагмента предпрогнозного периода, i, j=1, 2, 3;

i, j - эмпирические коэффициенты, учитывающие влияние на прогнозную урожайность озимой пшеницы фрагмента (продолжительностью 1, 2 и 3 года соответственно) предпрогнозного периода преобразованного ряда урожайностей, 1 [0,9; 2,2]; 2 [2,8; 5,1]; 3 [2,3; 3,7];

mi, mj - продолжительность фрагмента предпрогнозного периода преобразованного ряда урожайностей, mi, mj=1, 2 и 3 года;

Di, Dj - число вхождений фрагмента предпрогнозного периода преобразованного ряда урожайностей продолжительностью mi, mj лет соответственно;

n - количество лет ежегодных наблюдений.

Относительная повторяемость фрагментов предпрогнозного периода в формуле (2) вычисляют на основе выражения Di/((n-mi +1)/3m), где числитель - фактическое число вхождений фрагмента предпрогнозного периода преобразованного ряда урожайностей продолжительностью mi лет. Знаменатель - теоретически ожидаемое число вхождений фрагмента предпрогнозного периода преобразованного ряда урожайностей продолжительностью mi лет (соотношение получено на базе формул комбинаторики и математической статистики). Полученные значения нормируют с использованием эмпирических коэффициентов i.

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

Примеры реализации способа.

Пример 1.

В таблице 1 представлен многолетний временной ряд (более 27 лет) значений урожайности озимой пшеницы по Волгоградской области, отличающейся засушливым климатом. Перевод в лингвистические значения «Высокая урожайность», «Средняя урожайность», «Низкая урожайность» выполняют по формуле

где q1=miny+d - нижняя граница интервала величин урожайности (ц/га) для значения лингвистической переменной «Средняя урожайность»;

q 2=maxy-d - нижняя граница интервала величин урожайности (ц/га) для значения лингвистической переменной «Низкая урожайность»;

d=(maxy-miny)/3 - размер интервала величин урожайности (ц/га) для значений лингвистической переменной, определяемый как треть от размаха вариации величин урожайности за все годы наблюдений;

miny, maxy - наименьшее и наибольшее значения урожайности за все годы наблюдений соответственно.

Ряд лингвистических значений урожайности также представлен в таблице 1.

Далее выделяют фрагменты предпрогнозного периода преобразованного ряда продолжительность 1, 2 и 3 года (выделено заливкой в таблице 1).

После этого определяют количество появлений фрагментов предпрогнозного периода, предшествующих соответствующему лингвистическому значению «высокая», «средняя», «низкая» (таблица 2).

Прогнозную урожайность на следующий год (без учета поправочного коэффициента ki) определяют по формуле

где , - количество появлений фрагментов предпрогнозного периода, предшествующих лингвистическому значению «низкая», «средняя», «высокая» соответственно.

Например, для составления прогноза на 1997 год фрагменты предпрогнозного периода, количество появлений фрагментов предпрогнозного периода, предшествующих лингвистическому значению, и прогнозная урожайность на следующий год представлены в таблице 2. В этом случае для прогнозирования используется n=29 лет (число лет с 1968 г. по 1996 г.). В частности, количество появлений фрагмента предпрогнозного периода «низкая», предшествующих лингвистическому значению «средняя», (соответствующие сочетания выделены в таблицах 1 и 2 жирной рамкой). В используемом для прогнозирования ряду не встречаются фрагменты предпрогнозного периода продолжительностью 3 года.

Приведенные в формуле изобретения диапазоны значений эмпирических коэффициентов 1 [0,9; 2,2]; 2 [2,8; 5,1]; 3 [2,3; 3,7] получены для многолетних рядов урожайности озимой пшеницы, возделываемой в условиях засушливого климата, в результате обобщения результатов расчета на ЭВМ по критерию минимума среднеквадратического отклонения прогнозного значения от фактического. В расчетах приняты значения 1=1; 2=5; 3=3, принадлежащие рекомендуемым диапазонам.

Значение поправочного коэффициента ki, учитывающего относительную повторяемость i-го фрагмента предпрогнозного периода, определяют по формуле (2).

Численный расчет поправочных коэффициентов ki для i=1, 2, 3 по формуле (2) приведен ниже:

После этого рассчитывают прогнозируемую урожайность на следующий год (с учетом повторяемости фрагментов предпрогнозного периода в преобразованном ряду значений урожайности за все годы наблюдений и поправочного коэффициента, учитывающего относительную повторяемость фрагментов предпрогнозного периода) по формуле (1), ц/га

у=21,58·0,243+19,45·0,757+0·0=19,968.

По сравнению с фактической урожайностью за 1997 г. (18,8 ц/га) относительная погрешность прогноза составит 6,21%.

Пример 2.

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

В этом случае для прогнозирования используется уже n=36 лет (число лет с 1968 г. по 2003 г.). Также, в отличие от примера 1, в используемом для прогнозирования ряду встречается 1 раз фрагмент предпрогнозного периода продолжительность 3 года (выделено жирной рамкой в таблице 4).

Значение поправочного коэффициента ki , учитывающего относительную повторяемость i-го фрагмента предпрогнозного периода, определяют по формуле (2).

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

у=16,255·0,099+19,45·0,615+30,1·0,286=22,18.

По сравнению с фактической урожайностью за 2004 г.(19,3 ц/га) относительная погрешность прогноза составит 14,92%.

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

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

Способ прогнозирования урожайности озимой пшеницы, возделываемой в условиях засушливого климата, включающий многолетние ежегодные наблюдения за фактическими значениями урожайности и преобразование фактических значений урожайности в качественные величины, отличающийся тем, что фактические значения урожайности преобразуют в ряд лингвистических значений «Высокая урожайность», «Средняя урожайность», «Низкая урожайность», выделяют фрагменты предпрогнозного периода преобразованного ряда продолжительностью 1, 2 и 3 года, а прогнозируемую урожайность на следующий год находят с учетом повторяемости фрагментов предпрогнозного периода в преобразованном ряду значений урожайности за все годы наблюдений и поправочного коэффициента, учитывающего относительную повторяемость фрагментов предпрогнозного периода, по формуле где i - номер фрагмента предпрогнозного периода, i=1, 2, 3;yi - прогнозная урожайность на следующий год, определяемая с учетом повторяемости i-го фрагмента предпрогнозного периода в преобразованном ряду значений урожайности за все годы наблюдений;ki - поправочный коэффициент, учитывающий относительную повторяемость i-го фрагмента предпрогнозного периода, определяемый по формуле где i, j - номер фрагмента предпрогнозного периода, i, j=1, 2, 3; i, j - эмпирические коэффициенты, учитывающие влияние на прогнозную урожайность озимой пшеницы фрагмента продолжительностью 1, 2 и 3 года соответственно предпрогнозного периода преобразованного ряда урожайностей, 1 [0,9; 2,2]; 2 [2,8; 5,1]; 3 [2,3; 3,7];mi, mj - продолжительность фрагмента предпрогнозного периода преобразованного ряда урожайностей, mi, mj = 1, 2 и 3 года;Di , Dj - число вхождений фрагмента предпрогнозного периода преобразованного ряда урожайностей продолжительностью mi , mj лет соответственно;n - количество лет ежегодных наблюдений.

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

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

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