Файл: Количественные методы в мелиорации засоленных почв..pdf

стых случаев, при этом на поверхности почвы было принято условие, отражающее баланс солей:

где Ci — концентрация поливной воды.

Это условие впервые предложено Г. Бреннером (Brenner,. 1962) при исследовании конвективной диффузии в пористой среде и в дальнейшем было применено для рассмотрения явлений миграции солей в почве. Оно отражает равенство расходов солей, поступающих к поверхности почвы благо­ даря диффузии и фильтрации, с одной стороны, и фильтра­ ции поливных вод — с другой. Решение получено для полуограниченного пласта со следующими дополнительно к усло­ вию (1.2.12) условиями:

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

где k — коэффициент фильтрации; АН — разность уровней: воды над поверхностью поля и в дрене; Т — глубина залега­ ния водоупора; Ъ— расстояние между дренами; I — рас­ стояние от края промываемого чека до дрены; Н0— глуби­ на залегания уровня воды в дрене от поверхности поля. Ве­

Кроме того, Н. Н. Веригиным было предложено рассматри­ вать две стадии рассоления: первую — до достижения пол­ ного рассоления грунта на поверхности почвы, вторую — после образования зоны рассоления почвы у ее поверхности. Для второй стадии им предлагается решать задачу с нели-

14

шейными краевыми условиями и с подвижными границами области (Веригин, Шержуков, Шалинекая, 1967).

В случае, когда начальное распределение солей в почвогрунтовом слое неравномерно по глубине, оно аппроксими­ руется ломаной или ступенчатой функцией. Решение урав­ нения (1.2.3) получено Л. М. Рексом (1967, 1968, 1971) в

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

ции ; k, о — проницаемость и пористость грунта; Р(х, t) — давление в растворе; t — время; у — коэффициент еолеобмена, характеризующий скорость растворения солей; С* — концентрация водного раствора полностью насыщенного солями. Причем вязкость ц(С) и плотность р(С) раствора принимаются зависящими от концентрации:

ц(С)= ро+рС, р(С)=ро+аС.

Подбирая числовые значения коэффициентов (3 и а, можно добиться более реального совпадения величин р(С) и р(С) с действительными значениями.

Для устойчивого засоления (или рассоления) почв при малом содержании солей, находящихся в твердой фазе, можно 'воспользоваться уравнением (1.2.4) в виде

(1.2.18)

15

Г. П. Шапинской (Веригин, 1969) было получено решение этого уравнения для стационарного распределения солей для зоны аэрации. При этом скорость конвективного пере­ носа инфильтрациоиных вод и принималась постоянной, в общем случае не зависящей от влажности почвы в пределах рассматриваемой зоны. Условия на границе выбирались в следующем виде:

диффузии солей ;в зоне аэрации. Первое из этих условий отражает изолированность границы (ж= 0) для потока солей в атмосферу, а второе задает постоянную во все время про­ текания процесса концентрацию на границе области (х = 1).

М. П. Чиркиным и Д. Ф. Шульгиным (1964) для описа­ ния переноса солей в области капиллярной каймы при нали­ чии испарения рассматривается случай неглубокого залега­ ния уровня грунтовых вод, когда все соли выносятся в верх­ ние слои почвы, и в процессе испарения концентрация в них будет близкой к концентрации насыщения. Все оказанное можно записать следующим равенством: С = С п, имеющим место на поверхности почвы. Сп— концентрация предельно­ го насыщения раствора солями. Это равенство можно заме­ нить более точным условием, отражающим баланс солей на этой границе:

(1.2.21)

^ t + v C = °

Здесь р — пористость грунта. Условие (1.2.21) может быть выведено также из уравнения движения солей (Веригин, 1953). Аналитическое решение уравнения солепереноса с учетом члена, описывающего растворение солей, получено А. С. Хабировым (1969). Им был также исследован случай выноса солей из почвы вместе с испаряющейся влагой при

дС/дх = 0 (Хабиров, 1963).

При рассмотрении вопроса нестационарного изменения концентрации на границах области Н. Н. Веригиным берет­ ся в первом приближении диффузионное уравнение в виде

где V s ■— среднее значение скорости фильтрации.

16

Граничные условия принимаются (Веригин, 1969) сле­ дующие:

д * дС£ г) + U SC(0, 0 = 0 ,

C(l, t) —C i .

Начальное условие выбирается в предположении того, что влага в зоне капиллярной каймы имеет концентрацию грун­ товых вод, т. е. С(х, 0) = Сь

Для практических целей мелиорации представляет ин­ терес изучение солевого режима почвогрунтов как в есте­ ственно сложившихся условиях, так и при орошении не только в зоне полного насыщения грунтов водой, но и ниже уровня грунтовых вод. М. П. Чиркиным и Д. Ф. Шульгиным была предпринята попытка учесть зону полного насыщения грунта, а также влияние на распределение солей в этой зоне солевого режима области, лежащей выше уровня грун­ товых вод. Решение полученной задачи предполагает, что дренированноетъ почвогрунта и уровень грунтовых вод яв­ ляются постоянными величинами, при этом результирую­ щая скорость фильтрации (определенный по времени дре­ нажный модуль по М. П. Чиркину) имеет одинаковую вели­ чину по глубине как в зоне полного насыщения, так и в грунтовых водах.

В. А. Бароном (1967) была рассмотрена задача опреде­ ления зависимости между минерализацией грунтовых вод и почвенного раствора, когда покровная толща мелкоземов подстилается хорошо проницаемыми грунтами с напорными подземными водами, имеющими постоянную минерализа­ цию. Вертикальная скорость движения влаги рассматрива­ лась как средняя величина:

т ’

где V1— суммарное испарение; VQ— ежегодное поступление поверхностных вод; т — средняя объемная влажность слоя (эффективная влажность толщи грунтов, полностью насы­ щенных водой).

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

''' . !

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

В. М. Деловым, Н. М. Марченко и Э. А. Соколенко (1973) было предложено при изучении миграции солей в почвогрунте рассматривать единый почвенно-гидрогеологический процесс от дневной поверхности до первого регионального водоупора. Это позволило учесть многие важные для мелио­ рации параметры (как, например, скорость потока грунто­ вых вод, мощность водоносного пласта, глубину залегания грунтовых вод и т. д,). При аналитическом решении исполь­ зовались совместно уравнение конвективной диффузия (1.2.4), имеющей место в капиллярной кайме, и уравнение баланса солей в потоке грунтовых вод. Изменение дипрессионной кривой вдоль по потоку в стационарном состоянии находилось из решения уравнения для фильтрационного потока.

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

В практике мелиоративных рассчетов сравнительно ред­ ко встречаются случаи, когда фильтрация и массообмен

18