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

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

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

00000041.wmz

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

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

00000043.wmz

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

00000045.wmz

где: 00000046.wmz - относительная влажность; 00000047.wmz - текущее влагосодержание свободной (несвязанной) воды; 00000048.wmz - влажность, при которой движение влаги практически отсутствует; 00000049.wmz - полная влагоемкость.

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

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

00000051.wmz

где по литературным данным принимают 00000052.wmz

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

00000053.wmz

где: 00000054.wmz: 00000055.wmz - давление всасывания, м; 00000056.wmz - высота давления, м; параметр m для песков принимают равным 2, для тяжелых суглинков - 4.

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

00000057.wmz

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

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

00000059.wmz

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

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

00000062.wmz

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