Наиболее сильное засоление — непосредственно над дреной {0,3% по со-
держанию иона хлора}; примерно в два раза ниже — в середине междренья
(справа и слева} и еще более низкое (до 0,04% по СГ 
) — ниже глубины
заложения. После окончания промывки произошло значительное перераспре-
деление иона хлора: в придренной области содержание его снизилось более
чем в 20 раз; на междренье (с поверхности земли) — в 5...10 раз. Наиболее
интенсивное опреснение произошло в придренной десятиметровой области
(здесь изолиния с допустимой концентрацией по хлору О,01% проходит на
глубине чуть более } м}. Минерализация дренажного стока по плотному 0с-
татку снизилась с 20 до # г/л, а по иону хлора © #,75 до 1,7 г/л.
Для выполнения теоретических расчетов по перераспределению солей в
почвогрунте на фоне промывки была использована следующия исходная_ ин-
формация: коэффициент фильтрации — 0,5 м/сут; диаметр дрены — 0,25
м, глубина заложения верха дрены — 2,85; избыточное давление - 2,6; слой
воды на поверхности земли — 0,25; глубина залегания до водоупора — 50
м; концентрация подаваемой воды го хлору — О,005%; активная пористость
— 0,20; половина междренного расстояния — 25 м,
Кроме того, был еведен в матричной форме по исходному засолению поч-
ногрунтов коэффициент конвективной диффузии. Непрерывную функцию на-
чального засоления рассчитывали с помощью линейной аппроксимации
натурных данных. Коэффициент конвективной диффузии в рассматриваемой
модели принят как функция скорости, О = Эм + & И“ где Оы — коэффици-
ент молекулярной диффузии, равный для рассматриваемых условий 8-10”5
м*/сут; И - скорость фильтрации; 4, « — параметры, зависящие от геомет-
рми пористой среды.
Для определения параметров й и « были выполнены расчеты по опреде-
лению коэффициента конвективной диффузии при разных скоростях фильт-
рации и на раазном удалении от дрены по программе ЗА!МТ 1 (Рекс, Баранова,
1974} для однопараметрической модели. В середине междренья при скорости
филыраькли 0,00465 м/сут коэффициент конвективной диффузии составил
8...8-10°% м*/сут; в 10 м от дрены при вертикальной скорости О,012 м/сут
@ составил 1...3-0”' м" /сут; при скорости 0,043 м/сут О=0,1...0,06 ме/ут,
Полученные результаты аппроксимированы уравнением вида у = 12% +
+ 0.0008, где 4 = 1,23; а = 1, а коэффициент корреляции / = +0,89, При
аппроксимации квадратичной зависимостью параметры уравнения следую-
щие: 1 = З3?, а = 2; козффициент корреляции / = +0,78.
В результате расчета перераслределения солей, выполненного для тактов
промывки на 26, 37 и 68-е сут (в узлах сетки с шагами по у = 0,5, х = 5
м), получены (в матричной форме) вертикальная и горизонтальная состав-
ляющие, модуль вектора скорости, а также концентрация иона.
Сравнение теоретических изолиний равных концентраций по хлору на конец
промывки с натурными данными (как видно из рисунка} показывает, что по
картине и характеру рассоления теоретическое распределение аналогично на--
турным данным. Однако теоретическое рассоление происходит интенсивнее,
чем фактическое, Наиболее близкая сходимость результатов — для середины
междренья, где скорости фильтрации малы и преобладают вертикальные токи
воды. Отклонение натурных данных от теоретических в этом спучае не пре-
вышает 40 %. Значительное расхождение наблюдается в пределах придрен-
ной области ширинеой Ю м. Таким образом, эту модель можно использовать
только в случае двухмерной задачи для приближенной характеристики про-
цесса солепереноса.
2. Моделирование двухмерных задач влагосолепереноса в почвогрунтах на
ЭВМ,
Рассмотренные выше математические модели обладают рядом недостат-
ков: двухмерные модели насыщенной фильтрации не позволяют прогнози-
106