Научная статья на тему 'МОДУЛЬ «ПЛАНИРОВАНИЕ В УСЛОВИЯХ РИСКОВ» ПРОГРАММНОГО КОМПЛЕКСА «ЭКОЛОГО-МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ АГРАРНОГО ПРОИЗВОДСТВА»'

МОДУЛЬ «ПЛАНИРОВАНИЕ В УСЛОВИЯХ РИСКОВ» ПРОГРАММНОГО КОМПЛЕКСА «ЭКОЛОГО-МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ АГРАРНОГО ПРОИЗВОДСТВА» Текст научной статьи по специальности «Математика»

CC BY
16
10
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПРОГРАММНЫЙ КОМПЛЕКС / ЭКОЛОГО-МАТЕМАТИЧЕСКАЯ МОДЕЛЬ / ОПТИМИЗАЦИЯ / АГРАРНАЯ ПРОДУКЦИЯ / РИСКИ

Аннотация научной статьи по математике, автор научной работы — Ковалева Евгения Александровна, Иваньо Ярослав Михайлович

Предложено специальное информационное, математическое и алгоритмическое обеспечение для модуля «Планирование в условиях рисков» в составе программного комплекса «Эколого-математическое моделирование аграрного производства» - АГРЭКО. Описана параметрическая эколого-математическая модель для оптимизации производства аграрной продукции с учетом рисков, связанных с климатическими событиями. С помощью алгоритма выделения аномальных уровней путем построения многоуровневых трендов и стохастического анализа оценены производственные риски в будущем с определением ущербов от потери урожая. Реализованы тестовые задачи на примере сельскохозяйственной организации Иркутской области СХЗАО «Приморский» Нукутского района, деятельность которой подвержена значительному влиянию рисков, связанных с экстремальными климатическими явлениями, прежде всего, засухами. Многоуровневые тренды построены на примере временного ряда урожайности зерновых культур.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Ковалева Евгения Александровна, Иваньо Ярослав Михайлович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

MODULE "PLANNING UNDER RISKS" OF THE SOFTWARE PACKAGE "ECOLOGICAL AND MATHEMATICAL MODELING OF AGRICULTURAL PRODUCTION"

On the basis of the created software package, special information, mathematical and algorithmic support is proposed for the module "Planning under risk conditions" of the software package "Ecological and mathematical modeling of agricultural production". A parametric ecological and mathematical model for optimizing the production of agricultural products, taking into account the risks associated with climatic events, is described. With the help of an algorithm for identifying abnormal levels by constructing multilevel trends and stochastic analysis, production risks in the future are estimated with the determination of damages from crop loss. The test tasks were implemented on the example of the agricultural organization of the Irkutsk region of the Primorsky Agricultural District of the Nukutsky district, whose activities are significantly affected by the risks associated with extreme climatic events, primarily droughts. Multilevel trends are based on the example of a time series of grain yields.

Текст научной работы на тему «МОДУЛЬ «ПЛАНИРОВАНИЕ В УСЛОВИЯХ РИСКОВ» ПРОГРАММНОГО КОМПЛЕКСА «ЭКОЛОГО-МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ АГРАРНОГО ПРОИЗВОДСТВА»»

УДК 004.94: 504:631.1:045 DOI:10.38028/ESI.2022.27.3.013

Модуль «Планирование в условиях рисков» программного комплекса «Эколого-математическое моделирование аграрного производства»

Ковалева Евгения Александровна, Иваньо Ярослав Михайлович

Иркутский государственный аграрный университет имени А.А. Ежевского, Россия, Иркутский район, п. Молодежный, zhenia-93com@ya.ru

Аннотация. Предложено специальное информационное, математическое и алгоритмическое обеспечение для модуля «Планирование в условиях рисков» в составе программного комплекса «Эколого-математическое моделирование аграрного производства» - АГРЭКО. Описана параметрическая эколого-математическая модель для оптимизации производства аграрной продукции с учетом рисков, связанных с климатическими событиями. С помощью алгоритма выделения аномальных уровней путем построения многоуровневых трендов и стохастического анализа оценены производственные риски в будущем с определением ущербов от потери урожая. Реализованы тестовые задачи на примере сельскохозяйственной организации Иркутской области СХЗАО «Приморский» Нукутского района, деятельность которой подвержена значительному влиянию рисков, связанных с экстремальными климатическими явлениями, прежде всего, засухами. Многоуровневые тренды построены на примере временного ряда урожайности зерновых культур.

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

Цитирование: Ковалева Е.А. Модуль «Планирование в условиях рисков» программного комплекса «Эколого-математическое моделирование аграрного производства» / Е.А. Ковалева, Я.М. Иваньо // Информационные и математические технологии в науке и управлении. - 2022. - № 3(27). - С. 135-147. -DOI: 10.38028/ESI.2022.27.3.013.

Введение. Эколого-математические модели с критерием оптимальности в виде минимизации ущерба окружающей среде или максимизации доходов с учетом минимизации отрицательного воздействия на природную систему представляют собой научно-практический интерес в решении задач планирования производства аграрной продукции. На практике апробированы подобные модели с интервальными и вероятностными коэффициентами при неизвестных целевой функции и ограничениях. Вместе с тем сельское хозяйство в значительной степени подвержено воздействию экстремальных климатических явлений, которые увеличивают риски производства сельскохозяйственной продукции. Поэтому, помимо эколого-математических моделей, позволяющих оптимизировать деятельность сельскохозяйственного производителя для некоторых усредненных природно-климатических условий, в работе [1] предлагается модель, учитывающая производство продукции при формировании экстремальных событий.

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

Целью статьи является описание модуля «Планирование в условиях рисков» программного комплекса АГРЭКО, который позволяет оптимизировать производство аграрной продукции в условиях климатических рисков с минимальными ущербами водным и земельным ресурсам. Для достижения цели решались следующие задачи:

1) рассмотрение математического и алгоритмического обеспечения для минимизации ущербов окружающей среде при производстве аграрной продукции в условиях климатических рисков;

2) реализация модели минимизации ущербов окружающей среде при производстве продукции с учетом климатических рисков на реальном объекте с помощью разработанного программного комплекса АГРЭКО.

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

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

Для оценки свойств изменчивости показателей, входящих в эколого-математические модели оптимизации производства аграрной продукции, применены методы теории вероятностей и математической статистики [2-6]. В частности, к ним относятся способы построения законов распределения вероятностей, корреляционно-регрессионный, автокорреляционный анализ, метод статистических испытаний. При определении потерь плодородного слоя почвы использован метод оценки потенциальной опасности эрозии почв под воздействием дождей и расчетная формула интенсивности выноса почвы от ветровой эрозии [7-8]. Решение прикладных экстремальных задач основано на методах математического программирования в условиях неопределенности [9-11]. Кроме того, для построения модуля «Планирование в условиях рисков» использованы методы проектирования программных комплексов [12-14].

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

Результаты работы и их обсуждение. Авторами статьи в соавторстве создан программный комплекс АГРЭКО - с разработанным специальным математическим, алгоритмическим и информационным обеспечением [16, 17]. Программный комплекс АГРЭКО (рис. 1) решает прикладные задачи на основе информационного, математического, алгоритмического и программного обеспечения.

Рис. 1. Схема функционирования программного комплекса «Эколого-математическое моделирование аграрного производства» - АГРЭКО

Реализация интерфейса программного комплекса АГРЭКО и методов математической статистики осуществлялась при помощи интегрированной среды разработки Embarcadero Berlin 10.1. Алгоритмическое обеспечение представляет собой разработанные алгоритмы решения эколого-математических задач с использованием имитационного моделирования. Для реализации моделей планирования производства сельскохозяйственной продукции применен программный продукт LP Solve 5.5. Результаты статистической обработки данных, решений задач математического программирования и сведения из базы данных можно сохранять в виде таблиц в приложении MS Excel. Кроме того, с помощью программной среды реализованы различные методы оценки статистических параметров, критерии соответствия аналитических законов распределения вероятностей эмпирическим функциям и другие способы статистической обработки производственно-экономических, природно-климатических и экологических показателей.

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

На рисунке 2 приведен пример реализации оптимизационной эколого-математической задачи. В верхней части «Целевая функция» значения Х1, Х2,..., Х15 - отображают прибыль, получаемую от продажи 1 ц товарной продукции отрасли растениеводства или животноводства. Ниже приведены расчетные значения коэффициентов негативного воздействия техногенных и природных процессов на почву на богарных и орошаемых землях при возделывании сельскохозяйственных культур - 0,2093, 0,1938, ..., 0,2015. В нижней

части окна «Результат» показан определенный в процессе решения задачи оптимальный план со значениями целевой функции в виде прибыли и ущербов предприятия.

Оптимизация Файл Е!ыкод

Целевая функция Ограничения Размерность задачи Количество решений

U

X

Итераций для поиска вероятностного значения Район

15

Ш

Нукутский район

Поиск решения ® Максимум О Минимум

Целевая функция

Задача Детерменированнная

XI 1Я49 X2 912 X3 587 XI 1« X5 41 X6 238 X7 232 X8 183 X9 1125

< >

1 Зерновые 0,2093 2 Рапс 0,1938 3 Однолетние,, 0,1938 4 Однолетние,, 0,1938 5 Многолетни,, 0,0248 6 Многолетни,. 0,0248 7 Многолетни.. 0,0248 8 Силосные к.. 0,1938 9 Кукуруза на 0,2015

< >

Найти решение

Сжать результаты

Результат График Решение

Х14 Х15 Ущербы Значение целевой функции

91 580 3471205,671 22136686,94926

Рис. 2. Окно определения оптимального решения задачи эколого-математического моделирования продовольственной продукции с учетом ущербов от негативного влияния на окружающую среду для СХЗАО «Приморский» Нукутского района С помощью программного комплекса (ПК) АГРЭКО апробированы эколого-математические модели с детерминированными, интервальными и вероятностными коэффициентами при неизвестных целевой функции и ограничений с критерием оптимальности в виде минимизации ущерба окружающей среде или максимизации доходов с учетом минимизации отрицательного воздействия на природную систему. При этом ПК позволяет получать оптимальные решения для разных отраслей сельского хозяйства и их сочетания в условиях орошаемого и неорошаемого земледелия.

Отдельно следует выделить задачи математического программирования, позволяющие оптимизировать производство аграрной продукции в условиях проявления экстремальных гидрометеорологических и биологических событий [ 1].

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

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

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

Целевая функция:

I=I (1 - К )ср «)X - X (1 -1, +X (1 - ЪШ)х - X (1 - ООГ^ +

tel tel tel tel

+ Z(1 -ak )ckrk(t)yk - Z(1 -ak )dkrkœk ^ max,

keK keK (1)

c c

где t, ^ - прибыль, получаемая от продажи 1 ц товарной продукции культуры вида i при

c

богарном и орошаемом земледелии; k - прибыль, получаемая от продажи 1 ц товарной продукции животных вида k ; df, d-F - потери прибыли от продажи 1 ц товарной продукции культуры вида i при богарном и орошаемом земледелии, связанные с влиянием климатического события; - потери прибыли от продажи 1 ц товарной продукции

животных вида k, связанные с влиянием климатического события; Pt (t),Pt (t) -урожайность при богарном и орошаемом земледелии i-вида сельскохозяйственной культуры зависимая от параметра t; Xj, х' - площади богарных и орошаемых земель; ш' - площади

F 'F

богарных и орошаемых земель, подвергшиеся влиянию климатического события; Pi , Pi -потери урожайности на богарных и орошаемых землях i-вида сельскохозяйственной культуры под влиянием климатического события; rk(t) - продуктивность животных k, связанная с параметром t; rk - продуктивность животных k с учетом влияния климатического события; /¿, 1[ - коэффициенты негативного воздействия техногенных и природных процессов на почву на богарных и орошаемых землях при возделывании культуры i; ak - коэффициент негативного воздействия на окружающую среду при выращивании животных к; ук - количество животных; шк - количество животных, подвергшихся влиянию климатического события; F - вероятность события; I, K - число видов сельскохозяйственных культур и животных. Параметр t изменяется в интервале

[ Т1,Т2 ].

Ограничения по минимальному объему производства товарной растениеводческой продукции на предприятии выглядят так:

Z (1 - h )Pi (t)xt+Z (1 - Qrnt-Z (1 -h )pF ч -Z (1-Vp? > s - ZF - Z'F,

tF

iel iel iel iel (2)

где - минимальный объем производства растениеводческой продукции; Zf, Z'f -слагаемые уменьшения минимального объема производства растениеводческой продукции в условиях влияния климатического события на богарных и орошаемых землях.

Ограничения по минимальному объему производства животноводческой продукции имеют вид:

2(1-ак) rk(t)yk (1 ) rkœk *Sk -ZFk,

keK keK (3)

где Sk - минимальный объем производства животноводческой продукции; Zf - слагаемое уменьшения минимального объема производства животноводческой продукции в условиях влияния климатического события.

Ограничения по увязке производства растениеводческой продукции и потребности в ней животноводства записываются так:

X Кук > Б, (I е I),

кеК (4)

где ^ - потребность животных вида к в корме с использованием культуры /. Ограничения по наличию трудовых ресурсов имеют вид:

XЬх + ХЬХ + +X+ XЫч + * в,

1е1 1е1 кеК 1е1 1е1 кеК (5)

где Ь, Ь - затраты трудовых ресурсов на обработку 1 га богарных и орошаемых земель соответственно; Ь[ и Ь[р - дополнительные затраты трудовых ресурсов на обработку 1 га богарных и орошаемых земель с учетом влияния климатического события; Ьк, - затраты трудовых ресурсов для ухода за животными без учета и с учетом влияния климатического события; В - имеющиеся трудовые ресурсы.

Условия по предельно допустимой концентрации (ПДК) некоторых вредных веществ в почве выглядят так:

Х%(Х ) + х^(х'+ ч) + Хуи(х + ч) + ХК(х'+ ч) * П (/ е 1\

1е1 1е1 ¡е1 ¡е1 (6)

где , - исходная концентрация вредного вещества у на богарных и орошаемых землях,

соответственно; Уч' - концентрация вредного вещества у, попавшего на богарные и П .

орошаемые земли; J - величина предельно допустимой концентрации вредного вещества у

в почве; 3 - количество вредных веществ. Показатели ^ч, ^ч, и являются приведенными к единице площади сельскохозяйственных угодий. Ограничения на водозабор в реке имеют вид

X * Т£

1е1 (7)

где - оросительная норма культуры /; Т - вегетационный период; £ - расход воды реки.

Условие по предельно допустимой концентрации некоторых вредных веществ в реке записывается в виде:

■£ + ^ XV(ql + Л)(х\ +щ ) + Л(х, +щ )] * WJ (, е I)

е , (8)

где ^ ■ - исходная концентрация вредного вещества у в реке; ^ ■ - концентрация вредного вещества у в единице объема возвратных вод хозяйства; X - осадки, выпадающие за

Ж, п

вегетационный период; ■ 0 заданные величины предельно допустимой концентрации вредного вещества у в реке.

Ограничение на почвенные потери от водной и ветровой эрозии имеет вид

X яи руср (х + ч)+X яи руср (X+ч)+X М Т * п,

1е1 Ш I е1 (9)

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

C P

фактор крутизны склона; • - фактор растительности и севооборота; • - фактор

эффективности противоэрозионных мероприятий; Mi - интенсивность выноса почвы, т/га в 1 час, Т - время, в течение которого почва разрушается (в часах). Для расчета параметров, входящих в ограничение (9,) используются нормативные методики. Условие неотрицательности переменных:

X > х» Щ, Щ> Ук > Щ ^ 0. (10)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Функции (t),(t) могут быть описаны в виде линейных и нелинейных выражений. Согласно [16], наиболее приемлемыми для описания урожайности сельскохозяйственных культур для Иркутской области являются логистическая, асимптотическая и степенная функции:

p max

Р• (t) = '

(1+eA) (11)

pi (t) = РГ - (РГ - РГ )е, (12)

Р' (t) = P0tЛ , (13)

max nmin p

где Pi Pi - верхняя и нижняя оценки урожайности культуры i; pi - коэффициенты,

.0

характеризующие степень роста показателя; р - значение функции при 1=\.

Особенностью предложенной эколого-математической модели оптимизации производства аграрной продукции (1)-(13) является оценка ущербов от негативного влияния техногенных и природных факторов в условиях природных рисков, а также применение задачи параметрического программирования в условиях неопределенности. Ущербы окружающей среде включены в качестве коэффициентов негативного влияния в критерий оптимальности и ограничения по установленным объемам производства аграрной продукции. В статье [18] приведена схема последовательных операций по определению коэффициента негативного влияния природных и техногенных факторов на почву и водную среду.

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

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

\) выделять события разными способами в зависимости от особенностей временного ряда;

2) оценивать их вероятности появления;

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

В работе [15] описан алгоритм выделения аномальных уровней временных рядов, характеризуемых значимыми трендами. К таким рядам относятся урожайности сельскохозяйственных культур, для которых приведены результаты моделирования. Идея разделения значений многолетних рядов на последовательности с разными уровнями, принадлежащая И.П. Дружинину с соавторами [19], использована для оценки динамики трех последовательностей - всех уровней, нижних и верхних значений.

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

Рис. 3. Выделение аномальных значений (событий) по многоуровневым трендам логистической функции всего ряда (1), последовательностей нижних (2) и верхних уровней (3) для урожайности зерновых культур в Нукутском районе по данным 2000-2020 гг.: х - аномальные значения, соответствующие высоким рискам; + - аномальное значение, связанное с благоприятными условиями; • - планируемый уровень насыщения

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

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

По эмпирическим данным, вероятность появления аномальных уровней колеблется от 0.043 до 0.261, а в соответствии с аналитической функцией распределения Пирсона III типа -в пределах 0.044 - 0.233. Отметим, что в приведенном примере распределение Пирсона III типа лучше согласуется с эмпирическими значениями по сравнению с нормальным законом. Число аномальных уровней, соответствующих высоким рискам, равно 7. Следует иметь в виду, что вид функции влияет на количество выделяемых аномальных уровней (табл. 1).

Рис. 4. Функция распределения Пирсона III типа разностей между фактическими значениями ряда и уровнями, полученными по тренду логистической функции последовательности нижних членов для урожайности зерновых культур в Нукутском районе по данным 2000-2020 гг.

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

зерновых по данным Нукутского района за 2000-2020 гг.

Показатель Год Логистическая функция Вероятность

Потери по трендам, ц/га 2024 5,5

Потери по аномальному уровню, ц/га 20\6 8,3 0,044

20\5 7,\ 0,068

20\8 5,9 0,\20

20\7 5,4 0,\38

2020 5,2 0,\67

20\0 3,8 0,\77

20\9 3,6 0,233

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

Пирсона III типа. Самые большие потери приходятся на 2015 и 2016 годы. Предложенный алгоритм позволяет прогнозировать возможные потери.

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

При реализации параметрической эколого-математической модели для СХЗАО «Приморский» Нукутского района в экстремальных условиях многоуровневые тренды построены по временному ряду урожайности зерновых культур.

В качестве неизвестных в задаче оптимизации производства аграрной продукции использованы площади десяти сельскохозяйственных культур. К ним относятся: зерновые (xi), рапс (x2), однолетние травы на сено (xj) и зеленый корм (хД многолетние травы на сено (хД зеленый корм (хб) и семена (xz), силосные культуры (xs), кормовая кукуруза (хр), кукуруза орошаемая (хю). Кроме того, в моделях использованы пять переменных отрасли животноводства: поголовье коров основного стада (хц), поголовье коров на откорме (xi2), количество лошадей основного стада (xij), количество лошадей на откорме (xi4) и количество пчелосемей (хи).

С помощью программного комплекса АГРЭКО в условиях риска были получены следующие результаты (табл. 2).

Таблица 2. Значения целевой функции f) в виде прибыли при оптимизации

производства аграрной продукции при сочетании орошаемого и богарного земледелия с минимизацией ущербов окружающей среде для СХЗАО «Приморский»

в условиях риска, тыс. руб.

Тренд Год Значение целевой функции f Экологические ущербы

Средний уровень тренда 2024 24244,6 3814,6

Нижний уровень тренда 2024 21559,9 3350,8

Средний уровень тренда 2020 23767,5 3812,8

Нижний уровень тренда 2020 21409,2 3319,2

Нижний уровень тренда 2016 21147,2 3264,4

Аномальное значение 19304,0 3755,5

Согласно полученным результатам, в 2024 году прибыль по нижнему уровню составит 21,6, по среднему уровню - 24,2 млн руб. В 2020 году эти значения соответствовали 21,4 и 23,8 млн руб.. Для аномального уровня урожайности зерновых (ситуация 2016 года) прибыль составит 19,3 млн руб. В этом случае потери прибыли с учетом экологических потерь относительно нижнего уровня составили 2,3 млн руб.

Оптимальный план с учетом влияния засухи на производство продукции отличается от аналога, полученного по тренду среднего уровня. Для 2020 года: в неблагоприятной ситуации уменьшение площади зерновых культур составило 3589 га (67%), однолетних трав на зеленый корм - 4439 га (50%), а площади силосных культур сократились до 0 (100%). Количество коров основного стада уменьшилось до 3744 голов (сокращение на 24%). При этом увеличились площади однолетних трав на сено до 13180 га (в 4 раза), кукурузы - на силос - до 1130 га (на 125%), а также поголовье коров на откорме - до 1846 голов (на 85%).

Задача решена с помощью модуля «Планирование в условиях рисков» программного комплекса АГРЭКО - «Эколого-математическое моделирование аграрного производства». С использованием пункта меню «Регрессия» построены многоуровневые тренды. Тренды можно получить в виде значений функций и графиков. Аномальные уровни или события определены благодаря пункту меню «Статистические показатели». Задача математического программирования решена с помощью пункта меню «Оптимизация». При этом расширены возможности для решения задачи параметрического программирования вводом функции выбора в ограничениях вида тренда и его значений.

Заключение. Предложено математическое, алгоритмическое и информационное обеспечение для модуля «Планирование в условиях рисков» программного комплекса АГРЭКО - «Эколого-математическое моделирование аграрного производства».

Рассмотрена модель оптимизации производства аграрной продукции на богарных и орошаемых землях при сочетании отраслей растениеводства и животноводства в условиях формирования сильных засух, реализованная в тестовом режиме на примере сельскохозяйственной организации СХЗАО «Приморский» Иркутской области.

Список литературы

1. Ivanyo Y.M., Petrova S.A., Kovaleva E.A. Ecological-Mathematical modeling in planning production of agricultural products in conditions of risks. IOP Conference series:earth and environmental science, virtual, online, 10-12 января 2022 года, Virtual, Online, 2022, p. 022083, DOI 10.1088/1755-1315/988/2/022083, EDN GUYDCD.

2. Андерсон Т. Статистический анализ временных рядов / Т. Андерсон: пер. с англ. - М.: Мир, 1976. -356 с.

3. Гриневич Г.А., Композитное моделирование гидрографов / Г.А. Гриневич, Н.А. Петелина, А.Г. Гриневич. - М.: Наука, 1972. - 182 с.

4. Статистическое моделирование и прогнозирование / под. ред. А.Г. Гранберга. -М.: Финансы и статистика, 1990. -383 с.

5. Хастингс Н. Справочник по статистическим распределениям / Н. Хастингс, Дж. Пикок. -М.: Статистика, 1980. -95 с.

6. Юзбасиев М.М. Статистический анализ тенденций и колеблемости / М.М. Юзбасиев, А.М. Манелл. -М.: Финансы и статистика, 1998. -207 с.

7. ГОСТ 17.4.4.03-86. Метод определения потенциальной опасности эрозии под воздействием дождей : издание официальное: утвержден и введен в действие Постановлением Государственного комитета. СССР по стандартам от 10 ноября 1986 г. № 3401 : дата введения 1987-07-01 / разработан Государственным агропромышленным комитетом СССР. - Москва: Стандартинформ, 1986. - 8 c.

8. Волков С.Н. Землеустройство. - Т. 2. Землеустроительное проектирование. Внутрихозяйственное землеустройство / С.Н. Волков. -М.: Колос, 2001. -648 с.

9. Кирилюк, В.С. Меры риска в задачах стохастической оптимизации для получения робастных решений / В.С. Кирилюк // Стохастическое программирование и его приложения. - 2012. - С. 104-124.

10. Кнопов П.С. Метод эмпирических средних в задачах стохастической оптимизации и оценивания / П.С. Кнопов // Стохастическое программирование и его приложения. - 2012. - С.12 -149.

11. Кузьменко В.Н. Решение двухэтапных задач стохастического программирования большой размерности PNK-методом / В.Н. Кузьменко // Стохастическое программирование и его приложения. - 2012. -С. 150159.

12. Гагарина Л. Г. Разработка и эксплуатация автоматизированных информационных систем / Л.Г. Гагарина. - М.: ФОРУМ: ИНФРА-М, 2013. -384 с.

13. Липаев В.В. Программная инженерия. Методологические основы / В.В. Липаев. -М.: ТЕИС, 2006. -608 с.

14. Проектирование экономических информационных систем / под. ред. Ю.Ф. Тельнова. -М.: Финансы и статистика, 2003. -512с.

15. Иваньо Я.М. Об одном алгоритме выделения аномальных уровней временного ряда для оценки рисков / Я.М. Иваньо, С.А. Петрова // Актуальные вопросы аграрной науки, 2022. - №42 - С. 48-57.

16. Ковалева Е.А., Программный комплекс для эколого-математического моделирования производства сельскохозяйственной продукции / Е.А. Ковалева, Я.М. Иваньо // Информационные и математические технологии в науке и управлении, 2020. - № 2(18). - С. 94-106. - DOI 10.38028/ESI.2020.18.2.008

17. Иваньо Я.М. Программный комплекс моделирования природных и техногенных рисков / Я.М. Иваньо, С.А. Петрова // Известия Иркутской государственной экономической академии, 2015. - Т. 25 - №3 - С. 533-541. - Б01: 10.17150/1993-3541.2015.25(3).533-541.

18. Ковалева Е.А. Коэффициенты негативного влияния на окружающую среду в эколого-математическом моделировании аграрного производства / Я.М. Иваньо, Е.А. Ковалева // Научные исследования и разработки к внедрению в АПК: материалы науч.-практ. конф. молодых ученых (20 мая 2020 г.). -Иркутск: Изд-во ИрГАУ, 2020. - С. 153-160.

19. Дружинин И.П. Динамика многолетних колебаний речного стока. / И.П. Дружинин., В.Р. Смага, А.Н. Шевнин. -М.: Наука, 1991. -176 с.

Ковалева Евгения Александровна. Аспирант заочного обучения, SPIN-код: 2446-8670, zhenia-93com@ya.ru, Россия, Иркутский район, п. Молодежный.

Иваньо Ярослав Михайлович. Д.т.н., профессор, SPIN-код: 9654-8057, AuthorID: 106394, iymex@rambler.ru, Россия, Иркутский район, п. Молодежный.

UDC 004.94: 504:631.1:045 DOI:10.38028/ESI.2022.27.3.013

Module "Planning under risks" of the software package "Ecological and mathematical modeling of agricultural production"

Eugenia A. Kovaleva, Yaroslav M. Ivanyo

Irkutsk State Agrarian University named after A.A. Ezhevsky, Irkutsk district, Molodezhny, zhenia-93com@ya.ru

Abstract. On the basis of the created software package, special information, mathematical and algorithmic support is proposed for the module "Planning under risk conditions" of the software package "Ecological and mathematical modeling of agricultural production". A parametric ecological and mathematical model for optimizing the production of agricultural products, taking into account the risks associated with climatic events, is described. With the help of an algorithm for identifying abnormal levels by constructing multilevel trends and stochastic analysis, production risks in the future are estimated with the determination of damages from crop loss. The test tasks were implemented on the example of the agricultural organization of the Irkutsk region of the Primorsky Agricultural District of the Nukutsky district, whose activities are significantly affected by the risks associated with extreme climatic events, primarily droughts. Multilevel trends are based on the example of a time series of grain yields.

Keywords: software package, ecological and mathematical model, optimization, agricultural products, risks

References

1. Ivanyo Y.M., Petrova S.A., Kovaleva E.A. Ecological-Mathematical modeling in planning production of agricultural products in conditions of risks. IOP Conference series:earth and environmental science, virtual, online, 10-12 января 2022 года, Virtual, Online, 2022, p. 022083, DOI 10.1088/1755-1315/988/2/022083, EDN GUYDCD.

2. Anderson T. Statisticheskij analiz vremennyh ryadov [Statistical analysis of time series]: per. s angl. M., Mir [World], 1976, p.356.

3. Grinevich G.A., Petelina N.A., Grinevich A.G. Kompozitnoe modelirovanie gidrografov [Composite modeling of hydrographs]. M., Nauka [The science], 1972, p.182.

4. Statisticheskoe modelirovanie i prognozirovanie [Statistical modeling and forecasting]. pod. red. A.G. Granberga, M., Finansy i statistika [Finance and statistics], 1990, p.383.

5. Hastings N., Pikok Dzh. Spravochnik po statisticheskim raspredeleniyam [Handbook of statistical distributions]. M., Statistika [Statistics], 1980, p. 95.

6. Yuzbasiev M.M., Manell A.M. Statisticheskij analiz tendencij i koleblemosti [Statistical analysis of trends and fluctuations]. M., Finansy i statistika [Finance and statistics], 1998, p. 207.

7. GOST 17.4.4.03-86. Metod opredeleniya potencial'noj opasnosti erozii pod vozdejstviem dozhdej [Method for determining the potential hazard of erosion from rainfall]. Moscow, Standartinform Publ., 1986, p. 8.

8. Volkov S.N. Zemleustrojstvo. Vol. 2. Zemleustroitel'noe proektirovanie. Vnutrihozyajstvennoe zemleustrojstvo [Land management. Vol. 2. Land management design. On-farm land management]. M., Kolos [Ear], 2001. p. 648.

9. Kirilyuk V.S. Mery riska v zadachah stohasticheskoj optimizacii dlya polucheniya robastnyh reshenij [Risk measures in stochastic optimization problems for obtaining robust solutions]. Stokhasticheskoye programmirovaniye i yego prilozheniya [Stochastic programming and its applications], 2012, pp. 104-124.

10. Knopov P.S. Metod empiricheskih srednih v zadachah stohasticheskoj optimizacii i ocenivaniya [The method of empirical averages in stochastic optimization and estimation problems]. Stokhasticheskoye programmirovaniye i yego prilozheniya [Stochastic programming and its applications], 2012, pp. 125-149.

11. Kuz'menko V.N. Reshenie dvuhetapnyh zadach stohasticheskogo programmirovaniya bol'shoj razmernosti PNK-metodom [Solving two-stage large-dimensional stochastic programming problems using the PNK method]. Stokhasticheskoye programmirovaniye i yego prilozheniya [Stochastic programming and its applications], 2012, pp. 150-159.

12. Gagarina L. G. Razrabotka i ekspluataciya avtomatizirovannyh informacionnyh system [Development and operation of automated information systems]. M., FORUM: INFRA-M, 2013, p. 384.

13. Lipaev V. V. Programmnaya inzheneriya. Metodologicheskie osnovy [Software engineering. Methodological basis]. M., TEIS, 2006, p. 608.

14. Proektirovanie ekonomicheskih informacionnyh system [Design of economic information systems]. pod. red. Yu. F. Tel'nova. M., Finansy i statistika [Finance and statistics], 2003, p. 512.

15. Ivan'o Ya.M., Petrova S.A. Ob odnom algoritme vydeleniya anomal'nyh urovnej vremennogo ryada dlya ocenki riskov [On one algorithm for the allocation of anomalous levels of a time series for risk assessment]. Aktual'nyye voprosy agrarnoy nauki [Topical issues of agricultural science], 2022, no. 42, pp. 48-57.

16. Kovaleva E.A., Ivan'o Ya.M. Programmnyj kompleks dlya ekologo-matematicheskogo modelirovaniya proizvodstva sel'skohozyajstvennoj produkcii [Software package for ecological and mathematical modeling of agricultural production]. Informatsionnyye i matematicheskiye tekhnologii v nauke i upravlenii [Information and mathematical technologies in science and management], 2020, no. 2(18), pp. 94-106., DOI: 10.38028/ESI.2020.18.2.008.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

17. Ivan'o Ya. M., Petrova S.A. Programmnyj kompleks modelirovaniya prirodnyh i tekhnogennyh riskov [Software package for modeling natural and man-made risks]. Izvestiya Irkutskoy gosudarstvennoy ekonomicheskoy akademii [Proceedings of the Irkutsk State Academy of Economics], 2015, no.25-3, pp. 533-541,DOI: 10.17150/1993 -3541.2015.25(3).533-541.

18. Ivanyo Ya., Kovaleva E. Koefficienty negativnogo vliyaniya na okruzhayushchuyu sredu v ekologo-matematicheskom modelirovanii agrarnogo proizvodstva [Coefficients of negative influence on the environment in ecological and mathematical modeling of agricultural production]. Nauchnyye issledovaniya i razrabotki k vnedreniyu v APK: materialy nauch.-prakt. konf. molodykh uchenykh (20 maya 2020 g.) [Scientific research and development for implementation in the agro-industrial complex: materials of the scientific and practical conference of young scientists (May 20, 2020)], 2020, pp. 153-160.

19. Druzhinin I.P., Smaga V.R., Shevnin A.N. Dinamika mnogoletnih kolebanij rechnogo stoka [Dynamics of long-term fluctuations of river flow]. M., Nauka [The science], 1991, p. 176.

Kovaleva Evgeniya Aleksandrovna. Postgraduate student of correspondence studies, SPIN: 2446-8670, zhenia-93com@ya.ru, Irkutsk district, Molodezhny.

Ivanyo Yaroslav Mikhailovich. Doctor of Technical Sciences, Professor, SPIN: 9654-8057, AuthorID: 106394, iymex@rambler.ru, Irkutsk district, Molodezhny.

Статья поступила в редакцию 31.08.2022; одобрена после рецензирования 16.09.2022; принята к публикации 19.09.2022.

The article was submitted 08/31/2022; approved after reviewing 09/16/2022; accepted for publication 09/19/2022.

i Надоели баннеры? Вы всегда можете отключить рекламу.