06.01.00Агрономия
УДК 502/504:631.671.1:631/635 DOI 10.26897/1997-6011/2018-5-68-74
Я.А. АКСЁНОВ
Федеральное государственное бюджетное образовательное учреждение высшего образования «Брянский государственный аграрный университет», Брянская область, Российская Федерация
определение минимально допустимого предела влажности почвы при расчётах режима орошения
Цель исследований - определение предела испарения нижнего порога влажности почвы на основе экспериментов и расчётов с использованием математической модели. Предложена математическая модель, где нижний допустимый порог влажности характеризуется параметром Zaг, что соответствует замедлению испарения, которое может наступать в разные моменты времени суток. Почвенные колонки предварительно насыщались до предельной полевой влагоёмкости. Процесс испарения влаги в течение суток фиксировался при помощи измерения веса почвенных колонок, результаты наносились на графики. Параметр Z при Ь Л = +17°С в лаборатории составил: в почвенной колонке
1 ^ 1 1 да 1 воздуха 1 1
с толщиной слоя почвы 5см - 39,85мм, 10 см - 50,08мм, 15см - 59,81 мм, 20 см - 62,77мм. Затем колонки устанавливались в поле таким образом, чтобы уровень почвы в колонке соответствовал уровню земли, дневная температура воздуха составляла +38°С. Параметр Zда достиг: в колонке с толщиной слоя почвы 5 см - 18,69 мм, 10 см - 23,95 мм, 15 см - 25,67мм, 20 см - 27,05мм. Эксперименты показали, что при низкой температуре воздуха иссушение почвы происходит сильнее, чем при высокой температуре, то есть, процесс наступления снижения тургора растягивается во времени. Обработка данных по испарению с использованием математической модели позволяет получить допустимые пределы влажности почвы в течение суток и производить корректировку поливов.
Минимально допустимый предел влажности почвы, капилляры, почвенная колонка, испарение, математическая модель, поливы.
Введение. Оптимизация водного режима почвы достигается, прежде всего, совершенствованием технологий полива, играющим решающую роль в повышении урожайности, в эффективном использовании воды, удобрений и других ресурсов. Считается, что в интервале НВ — ВРК (наименьшая влагоёмкость — влажность разрыва капилляров) влага для растений легкодоступна и высокоэффективна. Тем не менее, ещё Костяков А.Н., потом Роде А.А. и др. считали, что величина этих пределов исследована недостаточно. [1, 2, 3, 9, 10]
Методы исследования. В лабораторных опытах Брянского ГАУ четыре почвенных колонки с разной толщиной слоя суглинка (5, 10, 15 и 20 см) насыщались водой до состояния предельной полевой влагоёмко-сти (ППВ), затем они испарялись в течение месяца с ежедневным контролем изменения веса при средней температуре окружающего
воздуха +17°С и постоянной влажности воздуха, отсутствии ветра и прямого солнечного излучения. В других опытах такие же почвенные колонки насыщались до ППВ, затем испарялись шесть суток с ежедневным контролем при средней температуре окружающего воздуха +30°С и постоянных влажности воздуха, отсутствии ветра и прямого солнечного излучения. В третьей серии эксперименты с испарением с почвенных колонок проводились на открытом воздухе при максимальной дневной температуре + 38°С. В четвёртой серии почва в колонках с толщиной слоя почвы 20 см нагревалась до температуры +33°С, поливная вода имела температуру +22°. С помощью рефлектора в лаборатории имитировалось прямое солнечное излучение через каждые два часа, затем производился контроль изменения веса. В опытах исследовался процесс физического испарения влаги, влияние растений не учитывалось.
Результаты исследования. Данные опытов приведены в таблице 1. Как видно из таблицы 1, максимальное количество испарившейся влаги Z<ю, соответствующей ВРК, достигается при температуре воздуха +17°С. К примеру, в слое почвы 20 см Z<ю в лаборатории при низкой температу-
ре составила 62,77 мм, а в реальных условиях при температуре воздуха +38°С всего 17,8 мм. Воздух, находящийся в капиллярах почвы, нагревается, и его вязкость, как известно из физики, увеличивается, его поступление в освободившиеся от воды в результате испарения капилляры, затрудняется.
Таблица 1
Испарения с почвенных колонок в различных условиях
Опыты Z0, мм Zm, мм 1/^1Нвв, 1/сУт
Слой почвы 5 см, S = 70,882 см2
Испарения с колонок при ^ = +17°С в лаборатории 22,14 39,85 0,111
Испарения с колонок при ^ = +30°С в лаборатории 8,46 35,55 0,308
Испарения с колонок при t = +38°С в реальных условиях А А возд А ^ 4,11 18,69 0,318
Слой почвы 10 см, S = 70,882 см2
Испарения с колонок при ^ = +17°С в лаборатории 14,1 50,08 0,105
Испарения с колонок при ^ = +30°С в лаборатории 11,29 39,93 0,302
Испарения с колонок при t = +38°С в реальных условиях возд 4,9 23,95 0,308
Слой почвы 15 см, S = 70,882 см2
Испарения с колонок при ^ = +17°С в лаборатории 15,51 59,81 0,097
Испарения с колонок при ^ = +30°С в лаборатории 28,22 41,2 0,297
Испарения с колонок при t = +38°С в реальных условиях возд 18,34 25,67 0,318
Слой почвы 20 см, S = 70,882 см2
Испарения с колонок при ^ = +17°С в лаборатории 14,44 62,77 0,077
Испарения с колонок при ^ = +30°С в лаборатории 26,81 43,88 0,308
Испарения с колонок при t = +38°С в реальных условиях возд 19,25 27,05 0,318
Испарения под рефлектором t = +32°С Гд = +22°С Гды = +33°С почвы 10,59 19,64 0,308
Вязкость газов определяется по формуле: п = 1 mnEX, где m - масса молекул,
и - средняя скорость молекул, А - длина свободного пробега между двумя соударениями молекул, n — число молекул в единице объёма.
Так как средняя скорость молекул возрастает с ростом температуры (также, как и «А»), то вязкость газов увеличивается при нагревании пропорционально л/Т, а вязкость воды - наоборот. Таким образом, образуется вакуум, тормозящий испарение.
Способствовать образованию вакуума может избыточное насыщение влагой верхнего слоя почвы вследствие атмосферных осадков. Это также подтверждается опытами А.А. Черкасова, А.М. Алпатьева, в которых при насыщении почвы водой свыше 80% скважности, аэрация почвы серьёзно затрудняется. [8]
№ 5' 2018
Температура воздуха в течение суток постоянно меняется — возрастает днём, снижается ночью. Соответственно, меняется величина вакуума, тормозящего испарение воды. Допустимый предел минимального содержания влаги в почве (ВРК по Роде А.А.) для дневных и ночных периодов в течение суток должен меняться. Этот предел не может быть константой.
В учебнике [4] приводится график изменений запасов влаги в расчётном слое почвы в течение вегетации, на котором отображается схема поливов за сезон. Первый полив производится до ППВ в конце второй декады мая, затем влажность почвы опускается практически до ВРК, после чего производится новый полив. Начиная с третьей декады июня уровень влажности доводится практически до ППВ, далее идёт испарение до ВРК.
Ю.А. Мажайский и др.[11] приводят графики динамики внутрисуточного
1б9)
водопотребления и метеофакторов, из которых видно, что в период роста температуры (в период с 6:00 утра до 14:00) растёт и водопотребление. Но после 14:00 и вплоть до 18:00 температура продолжает расти, а водопотребление начинает снижаться и к 19:00 достигает ВРК.
Движение почвенной влаги характеризуется схемой, представленной на рисунке 1.
Рис. 1. Схема процесса влагопереноса в поверхностном слое почвенной колонки: 1-1 - граница водосбора; 2-2 - рассматриваемое горизонтальное сечение; 3-3 - поверхность почвы в почвенной колонке; W - зона почвенного профиля, занятого влагой; Z - осушенный в результате испарения влаги слой
почвенной толщи в сечении 2-2; Нв - зона почвенной толщи, занятой влагой при насыщении до ППВ
Напишем уравнение баланса в дифференциальной форме:
lq>KcdZ =
KPWdt L
(1)
где К — коэффициент влагопроводности почвы; К, — коэффициент содержания влаги в почве в процентах от сухого веса при влажности, близкой к ППВ; Р — давление менисков, пропорциональное Z: Р = aZ; а — коэффициент пропорциональности.
Так как W = H - Z, то
в '
dz = Ка(Нв -Z)Z dt LlpKc
(2)
где ф — коэффициент, характеризующий форму кривой, ограничивающей капиллярные мениски.
Введём комплексный параметр, конКа
станту скорости: п = -
LlpKc
тогда
— = п (H -Z)Z
dt hy в '
(3)
Интегрируя уравнение (3) при t = 0 и Z = Z0, получаем:
Я
Z =
i+H-Z-
z0
-ViH„t
Полученное выражение описывает закономерность изменения запасов влаги в почве по времени в процессе испарения. Но в процессе испарения, одновременно с увеличением иссушенной зоны, происходит уменьшение мощности насыщенной зоны и уменьшение скорости снижения запасов влаги под действием вакуума до dZ/ ^ = 0.При этом бывшая часть иссушенной зоны под влиянием вакуума присоединяется к насыщенной зоне и в дальнейшем ничем не отличается от насыщенной зоны.
Скорость снижения запасов влаги описывается уравнением:
dZ = п (Hв -Z)Z-nZ2
(4)
Из условия равновесия dZ/dt = 0, когда Z ^ Z<ю, можно найти соотношения между равновесной мощностью иссушенной зоны Z<ю и мощностью насыщенной зоны, разделив обе части уравнения на n2Z<ю:
Z = H
п
(5)
Уравнение (4) с учётом(5) при начальных значениях t = 0 и Z = Z0, получим аналитическое выражение:
Z =
Z
1 +
(6)
Z
о
Z^ — количество влаги, испарившееся из колонки за бесконечный период времени, или количество влаги, которое должно покинуть слой почвы для достижения им минимально допустимого порога влажности.
С учётом вышеизложенного, можно увидеть, что данная математическая модель предлагает получать допустимый предел снижения влажности, характеризуемый параметром Z^, а не привязывать ВРК по процентному отношению к предельной полевой влагоёмкости.
По словам А.А. Роде, «В интервале влажности от НВ и ВРК достаточная подвижность влаги создаёт, по-видимому,
полную обеспеченность растений влагой. Мы говорим «по-видимому», так как прямых опытных данных для этого интервала влажности не имеется.» [10]
В книге А.А. Роде, [10] приведены кривые испарения при различных условиях окружающей среды. При высокой температуре воздуха, в жаркий день, транспирация является максимальной, точка снижения турго-ра наступает при 34% от НВ, интенсивность транспирации составляет 6,5 мм/сут. в частично облачный влажный день интенсивность транспирации составляет 4,2 мм/сут, точка снижения тургора наступает при 28% от НВ. В пасмурный влажный день точка снижения тургора наступает при 23% от НВ, интенсивность транспирации составляет 1,4 мм/сут. Как видно из анализа, точки снижения турго-ра при разных погодных условиях наступают в разные моменты. В наших опытах по испарению влаги из почвенных колонок при различных условиях окружающей среды, были получены аналогичные результаты.
Определим допустимый предел снижения влажности на примере расчёта испарений с почвенной колонки с толщиной слоя почвы 20 см по предложенной математической модели. Почвенные колонки были установлены таким образом, чтобы поверхность почвы в колонке была на одном уровне с поверхностью земли. Температура воздуха составляла +38°С. Контроль веса производился каждые 2 часа, в том числе и вночные часы — с 23:00 до 08:00 следующего дня. Следовательно, М принимаем равной 2 часам.В первые два часа испарений потеря веса колонки составила 0,03 кг, в следующие два часа - 0,049 кг, через 10 часов испарения составили 0,074 кг. Нанесём данные на график рисунка 3.
Рассчитаем параметры кинетической кривой испарения: скоростной коэффициент цхН и Z<ю, принимая недостаток до ППВ 0,01 кг.в(\ ="0,01 кг).
^ = 003^001 = 0,75;
^ =
0,033
0,06 - 0,033 0,06
0,07 - 0,06
= 0,45;
^ =
0,07
0,073 - 0,07 0,073
0,074 - 0,073 0,074
= 0,142;
= 0,041;
= 0,013
Наносим полученные точки на график зависимости Ф = А^) и определяем цхН и Z".
Рис. 2. Зависимость = f(Z) для дневного периода
По оси абсцисс прямая отсекает отрезок длиной 0,075, что соответствует Z<ю = 0,075 кг. По оси ординат прямая попадает на отметку 0,79, то есть = 0,79.
Далее рассчитаем скоростной коэффициент цхНвв по формуле:
в =-
1п (1 )= 1п (1 - 0,79)
М
2
= 0,78 1/сут
Для определения теоретических точек на графике рассчитаем Z в зависимости от времени 1.
£ =
0,075
1 +
0,075
0,01
0,01 в_ ^,
Наносим теоретические точки на график рисунка 3 и видим, что они совпадают с экспериментальными. Это говорит о высокой точности расчёта параметров с использованием данной модели, а небольшие отклонения теоретических точек на 1 гр можно объяснить погрешностью измерительного прибора.
Проведём расчёт допустимого предела снижения влажности для ночных испарений, полученныхв почвенной колонке с толщиной слоя почвы 20 см при температуре ночного воздуха +18°С и отсутствии осадков.
Принимаем М = 2 часа. В первые два часа испарений потеря веса колонки составила 1,3 гр (0,0013 кг), в следующие два часа - 2,9 гр, через 10 часов испарения составили 10,1 гр (0,01 кг). Нанесём данные на график рис. 4.
Вес грунта, кг 3,430 --
13:00 15:00 17:00 19:00 21:00 23:00
Рис. 3. Испарение при максимальной температуре воздуха +38°С О — теоретические точки + — экспериментальные точки
Рассчитаем параметры прямой для определения скоростного коэффициента Ц,Нвв и Z", принимая недостаток до ППВ 0,005 кг. (£0 = 5 гр).
6 3 — 5 = 0,0 0 = 0,206;
1 6 ^ = 7,92 - 6,3 = 0,204;
^ = т 3
¥ = т 4
^ 5 =
7,92 9,9 - 7,92 9,9 12,3 - 9,9
= 0,2;
12,3 15,1 -12,3 15,1
= 0,195; = 0,185.
Наносим полученные точки на график зависимости Ф = и определяем ц,Н и Z<ю.
По оси абсцисс прямая отсекает отрезок длиной 0,084, что соответствует Z" = 84 гр. По оси ординат прямая попадает на отметку 0,21, поэтому принимаем ф = 0,21.
Далее рассчитаем скоростной коэффициент ц,Нвв по формуле:
1п(1 -Ум ) 1п(1 -0,21) ¡иНв =--^ м =--Ц^—= 0,117 1/сут
At
2
Для определения теоретических точек на графике рассчитаем Z в зависимости от времени 1.
Рассчитываем Z, подставляя уже известные пераметры — Z0, Z<ю, е
Z = 84
- М1Н„
1 +
84 - 5
-1
5
Наносим теоретические точки на график рисунка 4 и видим, что они совпадают с экспериментальными.
Вес грунта, кг 3,400 3,390 3,380 3,370 3,360 3,350 3,340 3 330
00:00 02:00 04:00 06:00 08:00 10:00 Рис. 4. Испарение в ночное время О — теоретические точки + — экспериментальные точки
Аналогичные расчёты выполняем для дневного испарения при температуре воздуха +38°С, ночного испарения при температуре + 18°С и следующего за ним дневного испарения при температуре +24°С. Теоретические точки наносим на график (рис. 5).
Время суток, час
—I-1-1-
Вес грунта, кг 3,430
13:00 15:00 17:00 19:00 21:00 23:00 00:00 02:00 04:00 06:00 08:00 10:00 12:00 14:00 16:00 18:00 20:00
Рис. 5. Обобщённый график испарения за 1,5 суток в полевых условиях О — теоретические точки + — экспериментальные точки
Первый дневной период характеризуется предельным порогом снижения влажности - 75 гр или 7,22 мм (9,17% от ППВ), ночной период — 159 гр или 15,3 мм (19,4% от ППВ), второй дневной период - 121 гр или 11,64 мм (14,8% от ППВ).
Выводы
Исходя из полученных данных можно сделать вывод о том, что при различных метеорологических условиях нижний порог влажности почвы, характеризуемый параметром Z<ю, постоянно изменяется в течение суток. К концу дня испарение почти прекращается, а в ночное время снова восстанавливается.
В производственных условиях необходимо по датчикам влажности каждый час фиксировать, как идёт испарение, и осуществлять поливы в соответствии с почасовым ходом водопотребления.
Библиографический список
1. Система капельного орошения на землях Брянского ГАУ. / Н.М. Белоус, В.Е. То-риков, В.Ф. Василенков, С.В. Василенков, Е.В. Байдакова, Я.А. Аксёнов. // Вестник Брянской ГСХА. — 2017. — № 4. — С. 16-24.
2. Экологическая и экономическая оптимизация эксплуатационного режима орошения современными дождевальными машинами. / В.Ф. Василенков, С.В. Василенков, Ю.А. Мажайский, О.Н. Демина, Е.А. Мельникова. // Вестник Рязанского государственного агро-технологического университета имени П.А. Костычева. — 2015. — № 4 (28). — С. 85-92.
3. Василенков В.Ф., Демина О.Н., Мельникова Е.А. Разработка модели вла-гопереноса с целью планирования водопользования при орошении дождевальны-
ми установками. // Инновации в АПК: проблемы и перспективы. — 2016. — № 1. — С. 3-9.
4. Природообустройство. Учебник. /Подред.
A.И. Голованова. — СПб.: Лань, 2015. — 558 с.
5. Константинов А.Р. Испарение в природе. — Л.: Гидрометеорологическое изд-во, 1968. — 529 с.
6. Разработка методики, алгоритмов и модели информационно-аналитической системы оптимизации планирования водопользования на гидромелиоративных системах с использованием методов системного анализа: отчёт о НИР. / О.Н. Дёмина,
B.Ф. Василенков, С.В. Василенков и др. — Брянск: Брянский ГАУ, 2015. — 94 с.
7. Судницын И.И. Движение почвенной влаги и водопотребление растений. — М.: Изд-во МГУ, 1979. — 255 с.
8. Харченко С.И. Гидрология орошаемых земель. — Л.: Гидрометеорологическое изд-во, 1968. — 246 с.
9. Шуравилин А.В., Кибек А.И. Мелиорация. — М.: ЭКСМОС, 2006. — 944 с.
10. Роде А.А. Основы учения о почвенной влаге: монография, т. 1 — Л.: Гидроме-теоиздат, 1965. — 663 с.
11. Экологические аспекты орошения земель в условиях техногенного загрязнения. / Ю.А. Мажайский, В.И. Желязко, Н.Н. Дубенок и др. — М.: Изд-во МГУ, 2003. — 319 с.
Материал поступил в редакцию 25.04.2018 г.
Сведения об авторе
Аксёнов Яков Андреевич, ассистент кафедры природообустройства и водопользования, ФГБОУ ВО Брянский ГАУ; 243365, Брянская обл., Выгоничский р-н, с. Кокино, ул. Советская 2а; тел.: +7(919)2921742, е-mail: poivp@bgsha.com
YA.A AKSENOV
Federal state budgetary educational institution of higher education «Bryansk state agrarian university», Bryansk region, Russian Federation
DETERMINATION OF THE MINIMALLY PERMiSSiBLE SOiL MOiSTURE LIMIT AT CALCULATIONS OF THE IRRIGATION REGIME
The purpose of investigations is determination of the evaporation limit of the lower threshold of soil moisture on the basis of experiments and calculations using a mathematical model. There is proposed a mathematical model where the lower permissible humidity threshold is characterized by the parameter Zm corresponding to the slowdown in evaporation which can occur at different moments of the day. Soil columns were pre-saturated to the limit field water capacity. The process of moisture evaporation during the day was registered by measuring the weight of soil columns, the results were plotted. Parameter Zm at t air = + 17°C in the laboratory was: in the soil column
with the soil layer thickness 5cm - 39.85 mm, 10 cm - 50.08 mm, 15cm - 59.81 mm, 20 cm - 62.77 mm. Then the columns were installed in the field so that the soil level in the column corresponded to the ground level, the daily air temperature was + 38°C. Parameter Zx reached: in the column with the thickness of the soil layer 5 cm -18.69 mm, 10 cm - 23.95 mm, 15 cm - 25.67 mm, 20 cm - 27.05 mm. It turned out that at a low air temperature the soil drying occurs more strongly than at a high temperature, that is, the process of decrease occurrence in turgor stretches in time. Processing of data on evaporation using a mathematical model allows obtaining permissible limits of soil moisture during the day and making adjustments of watering.
Minimal acceptable limit of moisture content, capillaries, soil column, evaporation,
mathematical model, irrigation.
References
1. Sistema kapeljnogo orosheniya na zem-lyah Bryanskogo GAU. / N.M. Belous, V.E. To-rikov, V.F. Vasilenkov, S.V. Vasilenkov, E.V. Baidakova, Ya.A. Aksenov. // Vestnik Bryanskoj GSHA. - 2017. - № 4. - S. 16-24.
2. Ekologicheskaya i ekonomicheskaya optimizatsiya expluatatsionnogo rezhima orosheniya sovremennymi dozhdevaljnymi mashinami / V.F. Vasilenkov, S.V. Vasilenkov, Yu.A. Mazhaisky, E.A. Melnikova. // Vestnik Ryazanskogo gosudarstvennogo agro-tehnologicheskogo universiteta imeni P.A. Ko-stycheva. - 2015. - № 4 (28). - S. 85-92.
3. Vasilenkov V.F., Demina O.N., Melnikova E.A. Razrabotka modeli vlagopereno-sa s tseljyu planirovaniya vodopoljzovaniya pri oroshenii dozhdevaljnymi ustanovka-mi. // Innovatsii v APK: problem i perspek-tivy. - 2016. - № 1. - S. 3-9.
4. Prirodoobustrojstvo. Uchebnik. / Pod red. A.I. Golovanova. - SPb.: Lan, 2015. -558 s.
5. Konstantinov A.R. Isparenie v pri-rode. - L.: Gidrometeorologicheskoe izd-vo, 1968. - 529 s.
6. Razrabotka metodiki, algoritmov i mod-eli informatsionno-analiticheskoj sistemy op-timizatsii planirovaniya vodopoljzovaniya na gidromeliorativnyh sistemah s ispoljzo-vaniem metodov sistemnogo analiza: otch-
et o NIR. / O.N. Demina, V.F. Vasilenkov, S.V. Vasilenkov i dr. - Bryanskaya oblast: Bryansky gosudarstvenny agrarny universi-tet, 2015. - 94 s.
7. Sudnitsyn I.I. Dvizhenie pochvennoj vlagi i vodopotreblenie rastenij. - M.: Izd-vo MGU, 1979. - 255 s.
8. Kharchenko S.I. Gidrlogiya oroshae-myh zemel. - L.: Gidrometeorologicheskoe izd-vo, 1968. - 246 s.
9. Shuravilin A.V., Kibek A.I. Meliora-tsiya. - M.: EKSMOS, 2006. - 944 s.
10. Rode A.A. Osnovy ucheniya o pochvennoj vlage: monografiya, t. 1 -L.: Gidrome-teoizdat, 1965. - 663 s.
11. Ekologicheskie aspekty orosheniya zemel v usloviyah tehnogennogo zagryaz-neniya. / Yu.A. Mazhaisky, V.I. Zhelyazko, N.N. Dubenok i dr. - M.: Izd-vo Mosk. Un-ta, 2003. - 319 s.
The material was received at the editorial office
25.04.2018 g.
Information about the author
Aksenov Yakov Andreevich, assistant of the chair of environmental engineering and water use, FSBEI HE SAU; 243365, Bryanskaya region, Vygonichsky district, s. Kokino, ul. Sovetskaya, 2a, tel.: +7(919)2921742, e-mail: poivp@bgsha.com