Экология и водное хозяйство, No 2(02), 2019 г., [117—129]
В общем случае влагообмен грунтовых вод с зоной аэрации является
функцией влажности зоны аэрации: @; = /(Ф,), гдеф, — влажность зоны
аэрации. Функция зависимости потока от влажности зоны аэрации опреде-
ляется параметрами почвогрунтов [8].
Подставив в уравнение (11) выражения для потока (9) и объема вла-
ги (10), получим формулу для уровня грунтовых вод 1-го элемента:
И . 
К ‚АГ. 
КАГ 
КА. 
. 
.
а Г/ — 
1 
[ХЁі(КШ'_1 Ф [1—1 Т 
Ф21 
[1 + К...' 
Ф1А[1 Т КФ1+1%[1+1 ) +
@а — М(-о,) 
(МА 
(М; + М.
+К. 
КФі—1А1і—1 +КФ1'А11' 23_ —К . КФіА]і + КФ:'+1А[1'+1 23 
+Ч°а]-
" (МА 
О“О (МА
Взаимодействующий поток между рекой и грунтовыми водами
на }-м участке можно представить в виде:
г :No(у. - 7гм)Г.,.
/ 
А[„ 
л 
/
По изложенной выше модели были составлены алгоритм и програм-
ма расчета выклинивающегося расхода грунтовых вод для рассматривае-
мого участка реки.
Решение системы обыкновенных дифференциальных уравнений вы-
полняется методом Рунге — Кутта [9, 10]. Для численной проверки модели
были использованы данные, полученные в результате полевых исследова-
ний на участке р. Маныч в районе х. Тузлуков (таблица 1) сотрудниками
РосНИИПМ.
Таблица 1 — Данные полученные в результате полевых исследований
на участке р. Маныч в районе х. Тузлуков
Т, сут 
1 
| 
2 |3| 4 | 5 |6| 7 | $ |9| 10 | 11 | 12 | 13 | 14 | 15
О, м`/с 
6,2 | 5,8 | 6 | 5,4 | 4,9 | 5 | 3,8 | 4,5 | 4 | 4,6 | 5 | 4 | 3,8 | 4,4 | 3,4
Продолжение таблицы 1
16 | 17 | 18 | 19 | 20 | 21 
22 | 23 | 24 | 25 | 26 | 27 | 28 | 29 | 30
4,2 | 4,2 | 3,6 | 2,9 | 2,5 | 2,85 | 2,85 | 2,5 | 2,6 | 2,65 | 2 | 2,9 | 2,3 | 2,1 | 2,1