2.2. Фильтрация в зоне аэрации

Специфическим признаком зоны аэрации является неполное водонасыщение, тесно связанное с проявлением капиллярных сил.

При гравитационном влагопереносе в зоне неполного насыщения скорость фильтрации определяется законом, идентичным закону Дарси, в котором коэффициент пропорциональности существенно зависит от влажности.

Рисунок 41

где Рисунок 42 - коэффициент влагопереноса в зоне неполного насыщения (коэффициент фильтрации), м/сут.

Из теоретических соображений и анализа экспериментальных данных следует, что зависимость коэффициента влагопереноса от относительной влажности имеет степенной характер вида:

Рисунок 43

где Рисунок 44 - коэффициент влагопереноса при полном водонасыщении (коэффициент фильтрации), м/сут.

Рисунок 45

где: Рисунок 46 - относительная влажность; Рисунок 47 - текущее влагосодержание свободной (несвязанной) воды; Рисунок 48 - влажность, при которой движение влаги практически отсутствует; Рисунок 49 - полная влагоемкость.

Для показателя степени рекомендуется использовать значения n = 3 Рисунок 50 4.

В определенном диапазоне изменения влажности можно пользоваться экспоненциальной зависимостью коэффициента фильтрации от влажности:

Рисунок 51

где по литературным данным принимают Рисунок 52

В теории влагопереноса обычно вместо высоты давления используется величина давления всасывания. Зависимость коэффициента фильтрации (влагопереноса) от давления всасывания рекомендуется выражение:

Рисунок 53

где: Рисунок 54: Рисунок 55 - давление всасывания, м; Рисунок 56 - высота давления, м; параметр m для песков принимают равным 2, для тяжелых суглинков - 4.

Уравнение вертикального влагопереноса:

Рисунок 57

где Рисунок 58 - капиллярная диффузивность, или диффузивность почвенной влаги, м2/сут.

Скорость инфильтрационного просачивания рассчитывается:

Рисунок 59

где: Рисунок 60 - скорость инфильтрационного просачивания, м/сут.; w - интенсивность дождевания, м/сут.; Рисунок 61 - начальное значение влажности.

При n = 3 получено решение:

Рисунок 62

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