CC BY
43
Reference
1. Belyak, V.B. Metodicheskie rekomendacii po vozdelyvaniyu i semenovodstvu ryzhika [Tekst] I V. B. Belyak. - M.: Rossel'hozakademiya, 2004. - 40 s.
2. Belyak, V. B. Introdukciya ozimogo ryzhika [Tekst] I V. B. Belyak, T. Ya. PrahovaII Novye metody selekcii ozimyh kolosovyh kul'tur: sb. trudov. - Ufa, 2001. - S. 158-163.
3. Buyankin, V. I. Maslosemena ryzhika raznyh sortov kak istochnik nezamenimyh zhirnyh kislot v nashem pitanii [Tekst] IV. I. Buyankin, V. M. Fedorova II Nauchno-agronomicheskij zhurnal. -2013. - № 2. - S. 17-20.
4. Buyankin, V. I. Poseesh' ryzhik, pozhnesh' vygodu [Tekst] I V. I. Buyankin II Pole deyatel'nosti. - 2013. - № 5. - S. 50-53.
5. Buyankin, V. I. Osobennosti vyraschivaniya ryzhika i znachenie ego v sevooborotah Yu-zhnogo i Povolzhskogo federal'nyh okrugov strany [Tekst] I V. I. Buyankin II Nauchnoe obespechenie agropromyshlennogo kompleksa v sovremennyh usloviyah: materialy Mezhdunarodnoj nauchno-prakticheskoj konferencii i zasedaniya Soveta po vedeniyu zemledeliya v zasushlivyh usloviyah 9-10 iyunya 2014 g. - Volgograd, 2014. - S. 15-22.
6. Buyankin, V. I. Istoriya osvoeniya, tehnologiya, ]konomicheskaya ]ffektivnost' maslichnoj kul'tury ryzhik [Tekst] I V. I. Buyankin, G. I. Nesterenko. - Ural'sk: ZKF AO "NCNTI", 2014. - 88 s.
7. Voskresenskaya, G. S. Ryzhik [Tekst] I G. S. Voskresenskaya. - M.: Sel'hozgiz, 1952. - 47 s.
8. Kshnikatkina, A. N. Kak regulyatory rosta rastenij vliyayut na semena ryzhika [Tekst] IA. N. Kshnikatkina, T. Ya. Prahova, A. E. Safronkin II Fermer. Povolzh'e. - 2016. - № 8(50). - S. 15-18.
9. Rekomendacii po vozdelyvaniyu perspektivnoj maslichnoj kul'tury ryzhika posevnogo [Tekst] I V. V. Zubkov, O. V. Terent'ev, V. I. Buyankin. - Samara: Samara-ARIS, 2014. - 16 s.
10. Chekmarev, P. A. Introdukciya netradicionnyh maslichnyh kul'tur [Tekst] I P. A. Che-kmarev, A. A. Smirnov, T. Ya. Prahova II Dostizheniya nauki i tehniki APK. - 2013. - №7. - S. 3-5.
E-mail: agro@volgau.com
УДК 631.432.3
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ВЕРТИКАЛЬНОГО ДВИЖЕНИЯ ВЛАГИ В ЗОНЕ АЭРАЦИИ
NUMERICAL SIMULATION OF VERTICAL MOISTURE MOVEMENT
IN UNSATURATED ZONE
А.Н. Салугин, доктор сельскохозяйственных наук А.К. Кулик, кандидат сельскохозяйственных наук М.В. Власенко, кандидат сельскохозяйственных наук
A. N. Salugin, A.K. Kulik, M.V. Vlasenko
ФГБНУ «Федеральный научный центр агроэкологии, комплексных мелиораций и защитного лесоразведения Российской академии наук», г. Волгоград
Federal Research Centre ofAgroecology, amelioration and protective afforestation RAS, Volgograd
Численное решение модифицированного уравнения Дарси реализовано в среде имитационного моделирования Anylogic с использованием библиотеки системной динамики. В квазистационарном приближении рассмотрен перенос влаги в верхних слоях почвы после полива. Процесс вертикального движения воды в ненасыщенной почвенной толще представлен в виде уравнения Дарси, адаптированного к пористой среде. Капиллярный потенциал в приближении квазистационарности влагопереноса связан с влажностью линейной функции, что позволяет перейти к дифференциальному уравнению параболического типа (уравнение диффузии или теплопроводности). Это обстоятельство позволило применить традиционные методы для описания вертикального движения влаги после дождевания. Решение дифференциального уравне-
ния осуществлялось по явной схеме метода конечных разностей в имитационной модели системной динамики. Сеточная аппроксимация представлялась в виде разбиения зоны аэрации на почвенные слои, совпадающие с шагом сетки по вертикальной координате. Дискретизация по времени задавалась алгоритмом системной динамики Anylogic. Дальнейшее развитие метода предусматривает нестационарное моделирование, когда коэффициент влагопроводности зависит от влажности или почвенного потенциала.
The numerical solution of the modified Darcy equation is realized in the Anylogic simulation environment using the system dynamics library. In the quasi-steady-state approximation, the transfer of moisture in the upper layers of the soil after irrigation is considered. The process of vertical movement of water in an unsaturated soil stratum is represented as the Darcy equation, adapted to a porous medium. The capillary potential in the approximation of the quasistationary nature of moisture transfer is related to moisture by a linear function, which allows us to proceed to a differential equation of the parabolic type (the diffusion equation or the heat conduction equation). This circumstance allowed the use of traditional methods for describing the vertical movement of moisture after sprinkling. The solution of the differential equation was carried out according to the explicit scheme of the method of finite differences in the simulation model of system dynamics. The grid approximation was represented in the form of a division of the aeration zone into soil layers, which coincide with the grid spacing along the vertical coordinate. The time sampling was specified by the Anylogic system dynamics algorithm. Further development of the method involves non-stationary modeling, when the coefficient of moisture conductivity depends on humidity or soil potential.
Ключевые слова: влажность, влагопроводность, зона аэрации, численные методы, системная динамика, матричное давление.
Key words: humidity, moisture conductivity, aeration zone, numerical methods, system dynamics, matrix pressure.
Введение. Перемещение воды в ненасыщенной почве отражается на эффективности ирригации и имеет практическое значение для всех регионов России. Теоретически процессы влагопереноса в пористых средах изучены достаточно полно и математически описываются дифференциальными уравнениями в частных производных параболического типа [4, 5, 12, 9]. Особенно важны они при прогнозировании влагообеспече-ния зоны аэрации, уровня грунтовых вод, оценки норм полива и т.д. Следует отметить, что полученные к настоящему времени результаты в этой области малодоступны для рядового фермера, как и программное сопровождение этой проблемы в виде готовых интегрированных пакетов [2]. Это регламентируется не столько стоимостью программного продукта, сколько трудозатратами на освоение теоретических положений и алгоритма самой программы.
В связи с этим в данной статье мы предлагаем доступный и вместе с тем адекватный подход к решению задач вертикального водотока в зоне аэрации почв с легким гранулометрическим составом. Заметим, что большинство задач гидрологии в агроэкологии связано с параметрами объектов, которые изменяются во временных и пространственных координатах. В таких случаях создаются модели с распределенными параметрами, значения которых изменяются в пространстве решений. Они строятся, как правило, с использованием систем дифференциальных уравнений в частных производных [12].
Материалы и методы. В работе используется метод математического моделирования как метод исследования процессов влагопереноса. В качестве исходной использовалась упрощенная модель в виде уравнения диффузии. При численном решении дифференциальных уравнений методом сеток для пространственных (zi) и временных (ti) переменных задаются значения шагов их изменения. Затем строится конечно-разностная схема, связывающая коэффициенты и искомые влажности системой алгеб-
ИЗВЕСТИЯ"
№ 2 (50), 2018
раических уравнений. Далее задача сводится к решению полученной этой системы, которое определят значения влажностей в точках сетки. Каждому временному уровню ? соответствует свое значение м>(г,?{) для z-координаты в пространстве уравнений (горизонтов) зоны аэрации.
Результаты и обсуждение. Рассмотрим пример передвижения влаги в зоне аэрации с разными значениями влажностей на границах (рисунок 1). Требуется построить модель распределения влаги для внутренних точек разреза (профиля).
а а а
сс
&
я; Щ
О Г"
ш,
гг.
щ
Ж
Ь I
V
Щ
Рисунок 1 - Зона аэрации, разделенная на слои для численного расчета передвижения влаги: И - шаг интегрирования по координате z, ~м0 и м8 - начальные значения влажности на границах решения
В теории передвижения влаги в ненасыщенной зоне [1, 3, 7, 8] (в нашем случае -зоне аэрации) используется модифицированный закон Дарси в виде (1):
q(w) = -К (*м)
\
dP л — +1
V dz у
(1)
где q(w) - количество воды, прошедшей через единичную площадку за одни сутки (мм/сут.); К(м) - коэффициент влагопроводности (мм/сут.); Р - капиллярный почвенный потенциал (м); г - координата уровня с осью, поправленной вниз с началом на дневной поверхности (м).
Для скорости расхода влаги, передвигающейся по вертикали на основе классической термодинамики необратимых процессов, можно записать [6]:
Ж (м)
(2)
дм = К (м) Р2 М
д? 4 У дг2 дг где К(м) - коэффициент влагопроводности почвы в зоне аэрации.
Уравнение (2) получено с учетом принципа неразрывности потока для стационарного режима передвижения влаги. Зависимость влагопроводности от влажности К(м) различна для разных типов почв и зависит в основном от гранулометрического состава [13, 15, 16]. Для небольших временных интервалов уравнение (2) можно переписать в виде:
д(Дм>) ид 2(Дм) Л (ДК (м))
= В
+ д-
(3)
д? дг2 дг
где Дм,=м,-ме; АК(м)-К(м)-К(ме) (здесь м, К(м) - влажность и коэффициент влагопроводности в момент времени ?, а ме и К(ме) - их значения в начальный момент ?=0; В=2К(м)) [6].
Таким образом, в случае правомерности подобного допущения можно перейти к уравнению теплопроводности (диффузии), решая его методом конечных разностей. Для стационного случая, при постоянном коэффициенте влагопроводности (^(ч)=сом1), из уравнения (2) получим:
дч д2 ч
щ - Х •&*-■ (4)
где Х0 - коэффициент влагопроводности почвы, усредненный по всей зоне аэрации; ( - время.
Данное уравнение является дифференциальным уравнением в частных производных параболического типа и описывает распределение влажности во времени и пространстве в виде функции двух переменных: времени и координаты - (). Фактически это и есть уравнение, решение которого определяет динамику движения влаги в зоне аэрации в приближении стационарности инфильтрации. Для построения разностной схемы наложим сетку на область изменения независимых переменных в виде заданного шаблона. Множество узлов сетки будет использовано для аппроксимации частных производных.
Чтобы представить уравнение (4) в некоторой точке (.^,4), введем шаблон, изображенный на рисунке 2 и состоящий из четырех узлов (.ш, (п), . (п), . (п+Д ^и, (п), среди которых временной слой (п связан с (п+1 -м слоем дугой, ведущей к узлу . (п+1). Производную дw/дt заменим в точке разностным отношением:
дчг ^ +1 -
д( т
а производную д2ч/д.2 - второй разностной производной
д2 ч ^ чП+1 - чП + чП-1
И2
Рисунок 2 - Разностная схема для модели системной динамики (явная схема) [4]: - координата /-го слоя профиля, (п - п-й временной уровень наблюдения
Дополнительные (начальные и граничные) условия зададим на краевых узлах сетки. Пусть вначале (перед поливом) влажность по всей зоне аэрации была неизменной, то есть: ч (., (= 0) и равнялась МГ (максимальная гигроскопическая влажность). Таким образом, мы определили для нашей проблемы начальные условия (НУ). Граничные условия (ГН) зададим в виде: ч(.0,() =ч0; = Это означает, что величина на поверхности зоны аэрации z=z0 и на глубине z=zk заранее задана и не изменяется со временем. Такая упрощенная модель была взята с целью удобства при оценке адекватности результатов имитационного моделирования на начальном этапе его освоения.
60
Вместе с тем НУ и ГН такой задачи соответствуют случаю в практике орошения, когда реализуется стационарный во времени режим полива - краевые и начальные условия не меняются (условие 1-го рода).
Переходя к численной реализации модели, отметим, что метод конечных разностей используется с большим успехом при решении подобных дифференциальных уравнений. При этом применяются различные способы решения систем алгебраических уравнений с использованием линейной алгебры [12, 9, 6, 15, 16]. Как отмечено ранее, такие задачи решаются исследователями математического направления и недоступны для мелиоратора в практическом плане. В связи с этим мы сделаем попытку преодолеть это обстоятельство с помощью программы имитационного моделирования Апу^ю.
Пакет AnyLogic был разработан в России и к настоящему времени нашел широкое применение среди исследователей разных научных отраслей. Отметим, что, наряду с парадигмой системной динамики, AnyLogic поддерживает и другие концепции моделирования (дискретно-событийную, динамических систем и агентную) [2, 14, 10, 11]. Этот инструмент содержит средства аналитического решения уравнений в замкнутом виде, описывает изменение переменных во времени, учитывая модельное время, и содержит средства его масштабирования (имеются средства выражения логики и описания реакции систем на любые внутренние и внешние события). Одним из преимуществ AnyLogic является возможность наглядного представления поведения модели в виде анимации поведения системы и изучения ее характера во времени для различных значений параметров модели.
Непрерывные процессы описываются в AnyLogic заданием вещественных переменных формул и уравнений (алгебраических и дифференциальных) в привычной аналитической записи. Переменные определяют состояние динамического объекта и изменяются с течением времени по законам, определяемым дифференциальными уравнениями. Значения динамических переменных можно изменять и контролировать в ходе вычислительного эксперимента. У каждого экземпляра активного объекта существует свой набор переменных и параметров, поэтому поведение различных экземпляров одного и того же объекта могжет быть различным.
Для построения модели в AnyLogic в нашей модели необходимо задать выражения для накопителей влажности м в разных сечениях зоны аэрации и времени фильтрации (рисунки 1, 3).
Рисунок 3 - Слой искомой влажности с нормой i в разные моменты наблюдения
Частная производная второго порядка - правая часть уравнения (4) - приближенно выражается через значения влажности м в точке г и в двух соседних с ней: г-1 и г+1:
д2 м>
-1- [м,+1(?) - 2 м, (?) + м^)]
дг2 И2
Число таких уравнений будет равно числу пространственных сечений. Построим модель для шести сечений. Начальное значение для ~м0 возьмем 45 % - оно соответствует весовой влажности насыщения. Для нижнего слоя: м=5 % (максимальная гигроскопическая влага). Аналогично сформируем остальные накопители, начальное значе-
ИЗВЕСТИЯ"
№ 2 (50), 2018
ние для которых в виде постоянной влажности, равной нулю (сухая почва). Коэффициент влагопроводности зададим равным 0,1 мм/сут. Для объединения всех элементов в систему уравнений установим связи (рисунок 4).
Рисунок 4 - Математическая модель влагопроводности
С помощью временного графика и диаграммы результаты модели отображаются на экране: значения влажности в соответствующих накопителях в процессе анимации работы имитационной модели (рисунок 5). Из рисунка видно, что процесс продвижения влаги вниз носит нелинейный характер, замедляясь при приближении к границе раздела. В данном модельном варианте можно только качественно представить картину вертикального перемещения влаги. Сначала впитывание довольно интенсивно: на первых двух шагах от 45 % до 27 % , затем - до 15 %; внизу замедляется: 51 % ^ 4,5 % ^3,8 %. Если время перевести из единиц модельного в единицы реального времени, то динамика вла-гопереноса будет описана в том виде, который можно использовать на практике.
Жеес. %
27
Рисунок 5 - Результаты моделирования вертикального влагопереноса. Начальные условия: ^о=45 %, w7=3%. Коэффициент влагопроводности ^0,1 мм/ед. мод. врем.
Таким образом, если за единицу модельного времени взять 1 сутки, а коэффициент К0=0,1 м/сут., то полученный в нашей модели результат будет соответствовать именно этому случаю.
Заключение. Развитие данного подхода для нестационарных задач с распределенными параметрами позволит расширить область приложений имитационного моделирования в почвенной гидрологии. Конкретными задачами, для которых мы находим практическое применение при решении уравнения влагопереноса в системе AnyLogic, являются:
- экстраполяция экспериментальных данных по влажности от заданного уровня вниз по глубине до грунтовых вод;
- оценка интенсивности влагопереноса на различных глубинах зоны аэрации;
- расчет изменения питания грунтовых вод при изменении краевых условий на поверхности (выпадение осадков).
Библиографический список
1. Глобус, А. М. Почвенно-гидравлическое обеспечение агроэкологических математических моделей [Текст]/ А.М. Глобус. - Л.: Гидрометиоиздат, 1987. - 427 с.
2. Карпов, Ю.Г. Имитационное моделирование систем. Введение в моделирование с AnyLogic 5 [Текст]/ Ю.Г. Карпов. - СПб.: БХВ-Петербург, 2006. - 400 с.
3. Качинский, Н. А. Физика почвы [Текст] / Н.А. Качинский. - М.: Высшая школа, 1979. - 357 с.
4. Кулик, А.К. Водный режим и баланс влаги песчаных земель Нижнего Дона [Текст]: диссер... канд. с.-х. наук: 06.03.03 / А.К. Кулик. - Волгоград: ВНИАЛМИ, 2005. - 143 с.
5. Кулик, Н.Ф. Водный режим песков аридной зоны [Текст]/ Н.Ф. Кулик. - Л.: Гидро-метеоиздат, 1979. - 277 с.
6. Лебедев, А.В. Методы изучения баланса грунтовых вод [Текст]/ А.В. Лебедев. - М.: Недра, 1976. - 223 с.
7. Нерпин, С.В. Физика почв [Текст]/ С.В. Нерпин, А.Ф. Чудновский. - М.: Наука, 1967. - 584 с.
8. Роде, А.А. Вопросы водного режима почв [Текст]/ А.А. Роде. - Л.: Гидрометеоиз-дат, 1978. - 215 с.
9. Салугин, А.Н. Динамическое моделирование деградационных процессов в агроэкологии [Текст] : диссер. докт. с.-х. наук / А.Н. Салугин. - Волгоград, 2006. - 313 с.
10. Салугин, А.Н. Влагопроницаемость ненасыщенных почвогрунтов аридной зоны [Текст]/ А.Н. Салугин, А.К. Кулик, М.В. Власенко // Российская сельскохозяйственная наука. -2017. - №1. - С. 21-24.
11. Салугин, А.Н. Восстановление гидрофизических характеристик почв с помощью математического моделирования [Текст]/ А.Н. Салугин // Пути повышения эффективности орошаемого земледелия. - 2017. - №66(2). - С. 205-209.
12. Тихонов, А.Н. Уравнения математической физики [Текст]/ А.Н. Тихонов, А.А. Самарский. - М.: Наука, 1977. - 736 с.
13.Шеин, Е. В. Курс физики почв [Текст]/ Е В. Шеин. - М.: МГУ, 2005. - 432 с.
14. AnyLogic. Учебное пособие по Enterprise Library. URL: http: // www.xjtek.com / products / anylogic5 / enterpriselibrarytutorial.pdf (дата обращения: 21.03.2018).
15. Gardner W. R., D. I. Hiller, J. Geophys. Res., 1962. - v. 67. - P. 4319.
16. Mualem, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media/ Y. Mualem / Wafer Resour. Res., 1976. - V. 12. - P. 513-522.
Reference
1. Globus, A. M. Pochvenno-gidravlicheskoe obespechenie agrojekologicheskih matematich-eskih modelej [Tekst]/ A. M. Globus. - L.: Gidrometioizdat, 1987. - 427 s.
2. Karpov, Yu. G. Imitacionnoe modelirovanie sistem. Vedenie v modelirovanie s AnyLogic 5 [Tekst]/ Yu. G. Karpov. - SPb.: BHV-Peterburg, 2006. - 400 s.
3. Kachinskij, N. A. Fizika pochvy [Tekst] / N. A. Kachinskij. - M.: Vysshaya shkola, 1979. -
357 s.
4. Kulik, A. K. Vodnyj rezhim i balans vlagi peschanyh zemel' Nizhnego Dona [Tekst]: disser... kand. s. - h. nauk: 06.03.03 / A. K. Kulik. - Volgograd: VNIALMI, 2005. - 143 s.
5. Kulik, N. F. Vodnyj rezhim peskov aridnoj zony [Tekst]/ N. F. Kulik. - L.: Gidrometeoiz-dat, 1979. - 277 s.
6. Lebedev, A. V. Metody izucheniya balansa gruntovyh vod [Tekst]/ A. V. Lebedev. - M.: Nedra, 1976. - 223 s.
7. Nerpin, S. V. Fizika pochv [Tekst]/ S. V. Nerpin, A. F. Chudnovskij. - M.: Nauka, 1967. - 584 s.
8. Rode, A. A. Voprosy vodnogo rezhima pochv [Tekst]/ A. A. Rode. - L.: Gidrometeoizdat, 1978. - 215 s.
9. Salugin, A. N. Dinamicheskoe modelirovanie degradacionnyh processov v agrojekologii [Tekst] : disser. dokt. s. - h. nauk / A. N. Salugin. - Volgograd, 2006. - 313 s.
10. Salugin, A. N. Vlagopronicaemost' nenasyschennyh pochvogruntov aridnoj zony [Tekst]/ A. N. Salugin, A. K. Kulik, M. V. Vlasenko // Rossijskaya sel'skohozyajstvennaya nauka. - 2017. -№1. - S. 21-24.
11. Salugin, A. N. Vosstanovlenie gidrofizicheskih harakteristik pochv s pomosch'yu ma-tematicheskogo modelirovaniya [Tekst]/ A. N. Salugin // Puti povysheniya jeffektivnosti oroshaemogo zemledeliya. - 2017. - №66(2). - S. 205-209.
12. Tihonov, A. N. Uravneniya matematicheskoj fiziki [Tekst]/ A. N. Tihonov, A. A. Samar-skij. - M.: Nauka, 1977. - 736 s.
13. Shein, E. V. Kurs fiziki pochv [Tekst]/ E. V. Shein. - M.: MGU, 2005. - 432 s.
14. AnyLogic. Uchebnoe posobie po Enterprise Library. - URL: http: // www.xjtek.com / products / anylogic5 / enterpriselibrarytutorial.pdf (data obrascheniya: 21.03.2018).
15. Gardner W. R., D. I. Hiller, J. Geophys. Res., 1962. - V. 67. - P. 4319.
16. Mualem, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media/ Y. Mualem / Wafer Resour. Res., 1976. - V. 12. - P. 513-522.
E-mail: saluginan@mail.ru
УДК 631.67.03
СОЗДАНИЕ ЭФФЕКТИВНЫХ ЭКОЛОГИЧЕСКИ БЕЗОПАСНЫХ ТЕХНОЛОГИЙ И ТЕХНИЧЕСКИХ СРЕДСТВ ОРОШЕНИЯ НОВОГО ПОКОЛЕНИЯ
THE CREATION OF EFFECTIVE AND ENVIRONMENTALLY-FRIENDLY TECHNOLOGIES AND TECHNICAL MEANS FOR IRRIGATION OF THE NEW GENERATION
С.Я. Семененко, доктор сельскохозяйственных наук М.Н. Лытов, кандидат сельскохозяйственных наук Е.И. Чушкина, кандидат сельскохозяйственных наук А.Н. Чушкин, кандидат технических наук
S.Y. Semenenko, M.N. Litow, E.I. Chushkina, A.N. Chushkin
ФГБНУ «Федеральный научный центр агроэкологии, комплексных мелиораций и защитного лесоразведения Российской академии наук», г. Волгоград
Federal scientific center for Agroecology, land reclamation and protective afforestation RAS,
Volgograd
В статье рассматриваются способы полива, раскрывается возможность применения активированной воды для полива, а также оснащение систем капельного орошения модулем электрохимической подготовки оросительной воды. В исследованиях рассматривается эффективность нового направления сельскохозяйственного использования электрохимически активированной воды для
