CC BY
47
УДК 556.14:556.332.2:556.136:556.164 С.О. Гриневский1
СХЕМАТИЗАЦИЯ СТРОЕНИЯ И ПАРАМЕТРОВ ЗОНЫ АЭРАЦИИ ДЛЯ МОДЕЛИРОВАНИЯ ИНФИЛЬТРАЦИОННОГО ПИТАНИЯ ПОДЗЕМНЫХ ВОД2
На основе обобщения опубликованных данных поинтервального опробования разрезов почв и зоны аэрации получены типовые зависимости всасывающего давления и коэффициента влагопереноса от влажности для основных интервалов почвенного разреза в зависимости от литологического состава пород при разных ландшафтных условиях на поверхности земли. Рассчитаны соответствующие им средние параметры аппроксимации кривых уравнениями Ван Генухтена и показано, что их использование не приводит к существенным ошибкам при моделировании инфильтрационного питания подземных вод.
Ключевые слова: зона аэрации, гидрофизические параметры, влагоперенос, моделирование, инфильтрационное питание подземных вод.
On the base of published soil profiles data the typical retention curves and hydraulic conductivity functions were obtained depending on soil depth, lithological characters and different surface landscape conditions. The corresponding average soil hydraulic parameters for M.Th. van Genuchten equations were estimated. These average soil parameters are suitable for the groundwater recharge estimations based on unsaturated flow modeling.
Key words: vadose zone, soil retention curve, unsaturated hydraulic conductivity, unsaturated flow modeling, groundwater recharge.
Введение. Инфильтрационное питание подземных вод за счет атмосферных осадков — сложный природный многофакторный процесс, непосредственная (in situ) количественная оценка которого весьма сложна и неоднозначна. В современной практике гидрогеологических исследований все большее распространение получают методы оценки инфильтрационного питания подземных вод на основе моделирования процессов трансформации атмосферной влаги на поверхности земли и последующего влагопереноса в зоне аэрации [Пашковский, 1985; Гриневский, Поздняков, 2010]. Такие модели, будучи, по сути, комплексными, поскольку рассматривают инфиль-трационное питание подземных вод как часть единого водно-балансового цикла [Шестаков, Поздняков, 2003], используются не только в гидрогеологии, но и в смежных науках — почвоведении [Глобус, 1987], гидрологии [Кучмент, Гельфан, 1993] и др., а также при разработке глобальных метеоклиматических моделей отдельных регионов [Гусев, Насонова, 2004].
Важнейшая задача при использовании таких комплексных моделей в региональном масштабе — их параметрическое обеспечение, которое может быть выполнено на основе схематизации типовых природных условий и обобщения данных фактических наблюдений. В частности, для моделирования процессов формирования инфильтрационного питания подземных вод необходимы схематизация разреза зоны аэрации и обоснование гидрофизических параметров
влагопереноса, экспериментальное определение которых весьма затруднительно [Бадов и др., 1988].
В статье на основе обобщения опубликованных данных о строении и гидрофизических характеристиках почвенных профилей рассмотрены принципы схематизации разреза зоны аэрации и обоснована возможность использовать обобщенные параметры модели влагопереноса для оценки инфильтрацион-ного питания подземных вод.
Принципы оценки инфильтрационного питания на основе моделирования влагопереноса в зоне аэрации. Принципиальная модель формирования инфильтра-ционного питания подземных вод состоит из двух взаимосвязанных блоков:
— трансформации осадков на поверхности земли;
— влагопереноса в зоне аэрации.
Первый расчетный блок определяет условие на верхней границе модели влагопереноса — расход поступления влаги в зону аэрации, формирующийся из атмосферных осадков с учетом временной изменчивости их поступления, которая определяется метеоклиматическими условиями, а также процессов задержания и испарения атмосферной влаги растительностью, снегонакопления и снеготаяния, формирования поверхностного стока. Собственно величина инфильтрации оценивается как расход влаги, поступающей на нижнюю границу модели влагопереноса в зоне аэрации, где задается либо заданный напор, отвечающий глубине залегания
1 Московский государственный университет имени М.В. Ломоносова, геологический факультет, кафедра гидрогеологии, доцент, канд. геол.-минер. н., e-mail: sogrin@geol.msu.ru
2 Работа выполнена при финансовой поддержке РФФИ (грант 08-05-00720-а).
уровня грунтовых вод (условие 1-го рода), либо условие 3-го рода, определяющее латеральный отток к ближайшей дрене или вертикальный водообмен с более глубокими водоносными горизонтами. Такая модель представлена уравнениями одномерной вертикальной насыщенно-ненасыщенной фильтрации от поверхности земли до уровня грунтовых вод и реализована в широко используемом программном коде НТОЯи8-Ш [йшйпек ег а1., 2005]:
ММН-"^ (1)
8= =Г| + (аШУГ, т = 1 - - при Л < 0,
(- Штпт Ь 4 ' \ п
(1а) (1б)
0 = 1 при й> 0;
где ю — объемная влажность пород; ramin — доля неподвижной влаги; ютах — полная влагоемкость; h — высота давления; k — коэффициент фильтрации; TRp — транспирация; a, n — эмпирические параметры.
Уравнения (1а)—(1б) в модели влагопереноса представляют собой численную аппроксимацию (в данной записи — по М.Т. Ван Генухтену [van Genuchten, 1980]) зависимостей всасывающего давления от влажности (основная гидрофизическая характеристика — ОГХ) и коэффициента влагопереноса от влажности (кривая влагопроводности). Входящие в эти уравнения переменные, собственно, и являются параметрами модели влагопереноса, а практическая сложность экспериментального определения кривых ОГХ и влагопроводности [Бадов и др., 1988], а также изменчивость гидрофизических свойств пород в разрезе зоны аэрации определяют проблемы параметрического обеспечения расчетных моделей.
Таким образом, цель исследований — обоснование расчетных параметров модели влагопереноса в зависимости от строения и свойств пород почвенного покрова и зоны аэрации, а также анализ возможности их использования для оценки инфильтрационного питания подземных вод.
Схематизация строения зоны аэрации. Мощность зоны неполного водонасыщения определяется глубиной залегания уровня грунтовых вод (УГВ). Для процессов инфильтрационного питания определяющей является верхняя часть зоны аэрации мощностью в среднем не более 5—6 м, где наиболее существенно изменяется влажность пород и происходит формирование восходящего или нисходящего потока влаги. На большей глубине зоны аэрации (при глубоком залегании УГВ) формируется «транзитный» режим влагопереноса, при котором расход нисходящего потока (инфильтрация) уже практически не меняется [Пашковский, 1985].
Наиболее динамичный режим влажности формируется в почвенной зоне до глубины ~0,6 м (в черноземах — до 1 м), где при любом литологическом составе грунтов наиболее существенно меняются их свойства — плотность, пористость и др., что обусловлено воздействием корневой системы растительности, а также деятельностью живых организмов (рис. 1). Примечательно, что на этой глубине свойства почвы существенно зависят от ландшафта (рис. 1, А, по [Федоров, 1977]) — в лесу они с поверхности гораздо менее плотные и обладают более высокой влагоем-костью, чем в условиях полевого ландшафта. Ниже по разрезу зоны неполного водонасыщения свойства пород становятся более однородными и типичными для их литологического состава и практически не зависят от растительности, развитой на поверхности. Эти же закономерности прослеживаются и по данным других авторов [Роде, 1965; Субботин, 1966; Вериго, Разумова, 1973] (рис. 1, Б) и присущи практически всем основным типам почв.
Таким образом, разрез зоны аэрации можно схематизировать тремя расчетными интервалами. Максимальная изменчивость водно-физических свойств характерна для почвенных горизонтов типа А (дернина, лесная подстилка, перегнойный и гумусовый горизонты) и Е (А2) — элювиальных почв до глубины ~0,2 м. Затем до глубины, в среднем составляющей 0,6 м, тенденция к увеличению объемного веса и уменьшению параметров влагоемкости почвы сохраняется, но менее интенсивна. Этому интервалу отвечают иллювиальные почвенные горизонты типа В—Вх. На глубине свыше 0,6 м, отвечающей горизонту B2 и подстилающим материнским породам (С), водно-физические свойства пород зоны аэрации практически не меняются и отвечают литологическому составу отложений. При этом первые два слоя отражают почвенные горизонты; их водно-физические свойства, помимо механического состава (типа почвы), могут зависеть и от характера ландшафта и растительности на поверхности земли. В разрезе черноземных почв мощность 1-го слоя (горизонта А) следует увеличить до 0,6 м, а интервал иллювиального горизонта принимать от 0,6 до 1 м (рис. 1, Б).
Для выделенных трех модельных слоев гидрофизические параметры принимаются средними по глубине, поскольку, как показало специальное модельное тестирование, учет более мелкой внутренней неоднородности зоны аэрации, особенно на большей глубине, практически не сказывается на итоговых величинах инфильтрационного питания [Пашковский, 1985].
Обоснование параметров влагопереноса для типовых разрезов зоны аэрации. Из-за сложности экспериментальных оценок параметров ОГХ их расчет возможен на основе так называемых педотрансфер-ных функций PTF (pedotransfer functions), в качестве которых наиболее распространены регрессионные уравнения, связывающие равновесные значения давления и влажности с основными физическими
29 ВМУ, геология, № 6
свойствами почв [Глобус, 1987], которые получены в результате обобщения лабораторных исследований. В зависимости от того, какой набор характеристик почвы используется в PTF, однозначность параметрической характеристики ОГХ может различаться [Chanzy et al., 2008].
Для практической оценки параметров ОГХ на основе PTF широко применяется программный код Rosetta (ver 1.1) [Schaap et al., 2001], использующий базу данных почвоведческой службы США (USDA), в этом коде расчет параметров ОГХ осуществляется по фактическому относительному содержанию песчаных, пылеватых и глинистых фракций, объемному весу почвы и характерным величинам влажности (наименьшая влагоемкость и влажность завядания), отвечающим всасывающему давлению 33 и 1500 кПа.
Правомерность использования PTF, реализованных в программном коде Rosetta, проанализирована на основе лабораторных определений зависимости всасывающего давления от влажности, выполненных под руководством автора студентом М. Ширниным. Для этого на территории Звенигородского учебного полигона геологического факультета МГУ, на первой надпойменной террасе р. Москва, сложенной верхнечетвертичными аллювиальными отложениями [Полевые методы..., 2000], на краю лесного массива отобраны 2 монолита ненарушенной структуры размером 30x30 см с глубины 0,5 и 1 м. Определение гранулометрического состава и водно-физических свойств пород образцов проведено в лаборатории грунтоведения геологического факультета МГУ (табл. 1).
При этом значение наименьшей влагоемкости (НВ) получено расчетным путем по формуле Мичурина [Роде, 1965]:
НВ = 1,35ОВ/УВ(1 — ОВ/УВ),
а влажности завядания — по эмпирическому уравнению [Campbell, 1985]:
ВЗ = 1,5396(НВ)2 + 0,0967НВ — 0,0044.
Экспериментальное определение зависимости всасывающего давления от влажности для рассматриваемых образцов, представленных тяжелой супесью (по классификации В.В. Охотина), выполнялось на капилляриметре [Бадов, Киселев, 1987] в лаборатории института ВНИИ ВОДГЕО под руководством А.А. Киселева. Хорошее соответствие расчетной кривой ОГХ экспериментальным значениям (рис. 2) достигается только при максимальном наборе водно-физических характеристик образца — показателях гранулометрического состава, плотности, наименьшей влагоемкости и влажности завядания.
В контексте исследований необходимо оценить возможность использовать типовые параметры ОГХ применительно к выделенным интервалам разреза почвы и зоны аэрации для основных видов почв, ландшафтов и растительности на поверхности.
Таблица 1
Водно-физические свойства исследуемых образцов
Процентное содержание фракций, мм Образец 1, глубина 0,5 м Образец 2, глубина 1,0 м
>0,5 6 8
0,5-0,1 17 20
0,1-0,05 50 31
0,05-0,01 13 25
0,01-0,002 7 7
<0,002 7 9
Удельный вес (УВ), г/см3 2,45 2,65
Объемный вес (ОВ), г/см3 1,34 1,77
Пористость 0,454 0,333
Наименьшая влагоемкость (НВ) 0,335 0,300
Влажность завядания (ВЗ) 0,2 0,163
Принципиально такая возможность обусловлена тем, что генетическое различие почв, как правило, не проявляется в форме их ОГХ, и именно это положение лежит в основе создающихся национальных банков почвенно-гидрофизической информации [Глобус, 1987]. Этим объясняются весьма схожие закономерности изменчивости водно-физических свойств почв с глубиной независимо от генетического типа почвы (рис. 1, А, Б), что позволяет ожидать и соответствующих относительно незначительных вариаций параметров ОГХ в пределах выделенных интервалов зоны аэрации, схожих по механическому (литологи-ческому) составу отложений.
Вместе с тем в ряде исследований отмечается значительная пространственная изменчивость почвенно-гидрофизических характеристик на незначительных по площади территориях [Глобус, 1987]. При этом, однако, возникает вопрос: насколько такая статистическая вариабельность параметров ОГХ значима с точки зрения формирования суммарного потока влаги через зону аэрации, который и определяет величину инфильтрационного питания подземных вод?
Таким образом, для обоснования параметров модели влагопереноса необходимо решить две задачи. Во-первых, обосновать типовые кривые ОГХ и оценить вариацию их параметров применительно к принятой схематизации разреза зоны аэрации, а во-вторых, проанализировать, насколько существенно дисперсия величин параметров ОГХ в модели влагопереноса проявляется в расчетных значениях инфильтрационного питания подземных вод.
Решение первой задачи проведено на основе анализа опубликованных материалов, среди которых выбраны наиболее полно охарактеризованные по механическому составу и водно-физическим свойствам почвенные разрезы, что позволило рассчитать параметры ОГХ на основе педотрансферных функций. Всего использовано около 30 почвенных разрезов, преимущественно для Нечерноземной зоны России,
Объемный вес, г/см3
ОМ 1 1.2 11 1.8 1.3
Полная влагоемкость, %
20 30 40 50 60 70
—•— суглинок-поле —о— суглинок-лес
—супесь-поле —-л— супесь-лес
—•— песок-поле — -X — песок-лес
Б
Объемный вес, г/см3
08 1 1.2 1.4 1.6 1.8 2
Полная влагоемкость, %
10 20 30 40 50 60 70 80
—)— -1-поле, залежь, дер ново-слаб о подзолиста я сред несуглинистая почва
--2-березовый лес, д ерно в о-с иль ко подзолиста я
суглинистая почва
—ь--3- пахотное поле, дерново-слабопод зол истая
сред несу глинистая почва х -4-смененный лес. дер но в о-средне подзолиста я
супесчаная почва
* —5- березовый лес, дер ново-средне подзолистая
тяжелосуглинистая почва о -6-дер но в о-подзолистый суглинок (Батищэво)
—□--7-дер но в о-подзолистый лылеватый суглинок
(Белогорка)
—■--8— пахотная дер нов о-сред не подзолиста я
тяжелосуглинистая почва
—«----9-лесная, дер ново-средне под золиста я
легкосуглинистая почва . -о- --10-целинный тучный суглинистый чернозем
- -о- - —11—пахотный южный солонцеватыйчернозем
Рис. 1. Изменение объемного веса и влагоемкости пород в зоне неполного водонасыщения с глубиной в зависимости от их лито-логического состава и ландшафта поверхности по данным: А — С.Ф. Федорова; Б — 1-5 — А.И. Субботина, 6-7 — С.А. Вериго и
Л.А. Разумовой, 8-11 — А.А. Роде
Рис. 2. Сопоставление экспериментальных (точки) и расчетных (линии) зависимостей всасывающего давления от влажности: 1 — образец 1, глубина 0,5 м; 2 — образец 2, глубина 1 м
для которых имеются поинтервальные показатели гранулометрического состава, характерные значения влажности и объемного веса почв по глубине почвенного разреза с дискретностью от 0,1 до 0,3—0,5 м [Роде, 1965; Субботин, 1966; Вериго, Разумова, 1973; Гусев и др., 2008]. Для них в соответствии с принятой схематизацией разреза зоны аэрации на три расчетных интервала по программе Яо8еИа получены параметры и рассчитаны кривые ОГХ и влагопроводности,
соответствующие их аппроксимации уравнениями (1а) и (1б).
Анализ расчетных кривых ОГХ, характеризующих свыше 200 интервалов почвенных разрезов, показал, что в пределах каждого почвенного горизонта (А, В, С) независимо от генетического типа почвы они группируются по механическому (гранулометрическому) составу. При этом выделено 4 основных типа разреза — песчаный (<6%), супесчаный (6—15%), суглинистый (15—30%) и глинистый (>30%) — по процентному содержанию глинистых частиц в почвенном интервале в соответствии с классификацией дисперсных грунтов В.В. Охотина.
Характер кривых ОГХ проанализирован также с учетом ландшафтных условий на поверхности, с выделением закрытого (с древесной, лесной растительностью) и открытого (полевого) типов ландшафта; в последнем отдельно рассматривались почвенные разрезы на многолетних травах (далее — луг) и на пахотных землях (далее — пашня).
Поскольку в соответствии с существующими представлениями почвенные гидрофизические характеристики имеют распределение, близкое к нормальному, за исключением коэффициента фильтрации (влагопереноса), распределенного логнормально [Глобус, 1987], каждый интервал разреза зоны аэрации в соответствии с литологическим типом разреза может быть охарактеризован средними значениями параметров ОГХ (табл. 2).
Таблица 2
Расчетные параметры ОГХ для выделенных интервалов разреза зоны аэрации в зависимости от литологического (гранулометрического) состава почвы и подстилающих пород в различных ландшафтах по результатам обобщения данных почвенных разрезов
Тип разреза / ландшафта Число определений Минимальная влажность, ют(Г| Максимальная влажность, ютах а, 1/м п к, м/сут
среднее КУ, % среднее КУ, % среднее КУ, % среднее КУ, % среднее* КУ, %
Горизонт А, 0—0,2 м
Суглинки / лес 3 0,037 30,4 0,571 8,4 2,84 40,3 1,317 1,2 0,56 51,7
Суглинки / пашня 7 0,060 29,2 0,472 19,5 1,51 45,0 1,305 4,4 0,053 61,5
Суглинки / луг 19 0,054 27,0 0,529 6,4 1,58 48,1 1,352 4,6 0,049 55,3
Супесь / лес 14 0,038 15,7 0,532 7,0 0,90 59,0 1,486 7,8 0,73 90,6
Супесь / пашня 7 0,044 15,4 0,514 5,6 0,80 59,3 1,645 6,9 0,62 45,6
Супесь / луг 5 0,042 38,4 0,494 8,2 0,72 85,6 1,631 13,7 0,27 27,7
Песок / лес 9 0,029 24,3 0,482 17,6 0,34 36,7 2,280 8,5 3,85 60,9
Песок / пашня 3 0,028 7,7 0,364 7,4 0,22 8,5 2,205 1,2 0,31 31,3
Песок / луг 6 0,030 31,8 0,413 11,9 0,55 40,9 1,744 3,8 1,50 21,3
Горизонт B, 0,2-0,6 м
Суглинки 30 0,055 17,6 0,427 6,8 0,60 48,8 1,352 7,4 0,07 48,5
Супеси 32 0,038 20,8 0,412 13,3 0,90 41,8 1,442 6,7 0,21 45,2
Пески 17 0,022 11,1 0,366 14,0 1,25 19,3 1,578 12,0 1,21 178,7
Горизонт С, >0,6 м
Глины - 0,057 - 0,380 - 0,74 - 1,145 - 0,007 -
Суглинки 48 0,052 19,2 0,383 7,9 0,49 29,0 1,269 6,1 0,04 83,1
Супеси 28 0,034 7,5 0,379 10,9 1,39 44,9 1,361 7,3 0,28 66,1
Пески 18 0,018 28,8 0,370 8,8 2,06 48,0 1,581 8,3 2,96 201,0
* Для коэффициента фильтрации значения среднего и вариации рассчитаны исходя из логнормального закона распределения.
ГоризонтА, 0-0,2 м
Горизонт А, 0-0,2 м
0,2 0,3 0,4 Влажность
—й— суглинок^ес
—л— сугл инок-пашня
■■■А— суглинок-лес
супесьлес
-О— супесь-пашня
---о--- супесьлес
—о— песок^аес
— 0— песок-пашня
- --п- - - песоклес
0,4 0,6 0,8 1
Относительная влажность
1.Е+00 1Б
> 1.Е-01
се 8
5 1.Е-02
о, о
I
§ 1.Е-03
Ё к
В 1.Е-04
Й 1.Е-05
Горизонт А
0,1 0,2 0,3 0,4 0,5 0,6 Влажность
Й 1.Е-05
0,2 0,4 0,6 0,8
Относительная влажность
Рис. 3. Расчетные кривые ОГХ (А, Б) и влагопроводности (В, Г) для верхнего горизонта зоны аэрации в зависимости от литологического состава пород и ландшафтных
условий на поверхности
Полученные средние значения и анализ соответствующих обобщенных кривых ОГХ и влагопроводно-сти для верхнего интервала (горизонта А) показывают (рис. 3), что здесь прослеживается зависимость параметров влагопереноса не только от литологического состава почвы, но и от типа ландшафтных условий на поверхности, что соответствует рассмотренным выше закономерностям изменчивости основных водно-физических свойств с глубиной (рис. 1, А). Из анализа кривых ОГХ видно, что роль ландшафтных условий на поверхности проявляется главным образом в характерных величинах влажности пород и значениях коэффициента фильтрации (рис. 3, А, В), в то время как сама форма кривых, определяющаяся значениями коэффициентов а и п уравнений Ван Генухтена (1а) и (1б), зависит преимущественно от литологического состава пород, что отчетливо проявляется на графиках в относительных величинах (рис. 3, Б и Г). При этом наименее значимое влияние ландшафтных условий на поверхности проявляется при суглинистом составе горизонта А.
Наиболее «консервативными» параметрами ОГХ с минимальными значениями коэффициента вариации (Ку) являются величины максимальной влажности и параметра п; максимальная вариабельность характерна для значений коэффициента фильтрации (табл. 2). Анализ вариаций параметров а и п, определяющих форму кривых ОГХ и влагопроводности, показывает, что значение а более стабильно в песчаных почвах горизонта А, а величина п, наоборот, в суглинистых (рис. 4). Также хорошо прослеживается уменьшение коэффициента а и увеличение параметра п при облегчении состава пород от суглинистого к песчаному.
Для второго рассматриваемого интервала почвенного разреза, отвечающего иллювиальному горизонту В, влияние ландшафтных условий на поверхности проявляется значительно слабее (рис. 5, А). Это вполне объяснимо с позиций генетического формирования иллювиального почвенного горизонта — за счет вымывания солей и наиболее тонких фракций из верхнего слоя, что обусловливает его относительную гидрофизическую однородность. На основании полученных кривых ОГХ для горизонта В можно заключить, что его гидрофизические характеристики следует рассматривать вне зависимости от ландшафтных условий на поверхности и в расчетах влагопереноса использовать осредненные ОГХ, отвечающие различным литологическим разностям. На рис. 5, Б для сравнения приведены также точки экспериментальной кривой ОГХ, полученной по образцу супесчаных пород, отобранному на Звенигородском полигоне с глубины 0,5 м (табл. 1), точки соответствуют усредняющей расчетной кривой для отложений этого типа.
Для горизонта В, так же как и для горизонта А, минимальным диапазоном изменчивости обладают параметр п и максимальная влагоемкость, а макси-
мальная вариация характерна для значений коэффициента фильтрации (табл. 2). Однако, в отличие от горизонта А, значения параметров а и п возрастают от суглинистых отложений к песчаным, при этом диапазон изменчивости коэффициента а существенно выше (рис. 4).
Для подстилающих почву материнских пород зоны аэрации (горизонт С), где водно-физические свойства пород гораздо менее изменчивы с глубиной (рис. 1), анализировать влияние ландшафтных условий на поверхности земли на форму кривых ОГХ и влагопроводности не имеет смысла. Для этих интервалов разреза (глубина свыше 0,6 м) также наблюдается однозначная дифференциация кривых в соответствии с литологическим составом отложений, что позволяет охарактеризовать их обобщающими параметрами (табл. 2, рис. 6). Хотя глубина анализировавшихся почвенных разрезов, как правило, не превышала 1,5 м, выявленные закономерности можно считать характерными и для более глубоких интервалов разреза зоны аэрации — до глубины залегания УГВ. Поскольку в анализируемых данных не описаны отложения, которые, согласно принятой типизации, следует отнести к глинам, параметры ОГХ для данного типа разреза были синтезированы с помощью педотрансферных функций кода Яо8ейа по относительной величине содержания глинистых частиц 0,45 и плотности 1,8 г/см3.
Расчетные кривые ОГХ и влагопроводности для нижнего интервала разреза зоны аэрации (горизонт С) наиболее контрастно (по сравнению с двумя верхними зонами) отражают различие литологического состава пород (рис. 6). Это объясняется тем, что на данной глубине породы зоны аэрации практически не изменены процессами почвообразования, корневыми системами растений и жизнедеятельностью живых организмов, в связи с чем породы обладают «нормальными» водно-физическими свойствами (плотность, общая пористость, влагоемкость и т.д.), типичными для их литологического (гранулометрического) состава. Иллюстрацией этого служит весьма схожий вид расчетных кривых для этого интервала и типовых кривых соответствующей литологии из базы данных Института почвоведения США (ШВА) ^сИаар е! а1., 2001] (рис. 6).
На графике (рис. 6) для сравнения также приведены точки экспериментальной кривой ОГХ, полученной по образцу супесчаных пород, отобранному на Звенигородском полигоне с глубины 1,0 м (табл. 1). Опытные точки отвечают диапазону возможной вариации усредняющей кривой супесчаных отложений, хотя при более высоких значениях давления они больше соответствуют кривой для суглинков.
Вариация параметров влагопереноса для этого интервала разреза зоны аэрации представлена в табл. 2 и на рис. 4. Здесь минимальным диапазоном изменчивости также обладают параметр п и максимальная влагоемкость пород.
6,0 -1
.4,0 -
Горизонт А, 0-0,2м
К 2,0
2,0-1 Горизонт В, 0,2-0,6 м
лес пашня луг лес пашня луг лес пашня луг суглинок супесь песо«
3,0 -1
лас паиня луг суглинок
пашня луг г в с пашня луг супесь песок
1.0 -
2,0-1
1,6"
« 1,2 &
0,4-
суглинок
супесь
4,0
3,0-
а &
2,0
с, я С
1,0-
Горизонте, >0,6 м
2,0 1,6-
с
¡.и
I ¿0,8-
0,4
глина суглинок супесь песок
глина суглинок супесь песок
суглинок
Рис. 4. Гистограммы средних значений параметров ОГХ (а и п) для горизонтов А, Б и С зоны аэрации для различных по литологическому составу пород
и ландшафтных условий
о
Рис. 5
Расчетные дифференцированные (вверху) и обобщенные (внизу) кривые ОГХ и влагопроводности для второго горизонта зоны аэрации в зависимости от литологического состава пород и ландшафтных условий на поверхности
Таким образом, в результате рассчитаны обобщенные кривые ОГХ и влагопроводности для выделенных литологических типов пород (рис. 7), их параметры и пределы вариации средних значений (табл. 2). При этом получено, что для верхнего почвенного интервала (горизонт А) определенное влияние на форму кривых и параметры ОГХ имеет также тип ландшафтных условий на поверхности. В целом в результате анализа обобщенных кривых ОГХ и влагопроводности (рис. 7) выявлено их существенное отличие для однотипных по литологическому составу пород в зависимости от рассматриваемых интервалов разреза зоны аэрации. Наиболее значимо это различие проявляется при высокой влажности.
Однако, как показал анализ, значения коэффициента вариации полученных средних параметров уравнения Ван Генухтена оказываются во многих случаях весьма значительными (табл. 2). При этом возникают вопросы: как влияет вариация параметров ОГХ на результаты модельной оценки значений инфильтрационного питания подземных вод и насколько правомерно использовать средние значения параметров при моделировании? Для выявления этого проведена серия модельных расчетов, в которых при одинаковых граничных условиях модели влагопере-носа (1), представленной тремя схематизированными горизонтами зоны аэрации, оценивалась величина
среднемноголетнего инфильтрационного питания (отток через нижнюю границу модели).
При этом на верхней границе заданы среднесуточные значения поступления влаги на поверхность почвы и величины потенциальной эвапотранспира-ции, полученные моделированием трансформации осадков на поверхности земли по программе 8игШа1 [Гриневский, Поздняков, 2010] с использованием многолетних суточных метеоклиматических рядов метеостанции г. Мосальск (Калужская область). На нижней границе модели поддерживался постоянный уровень грунтовых вод на глубине 6 м. В процессе моделирования, проведенного с суточным шагом по времени, последовательно изменялись расчетные параметры влагопереноса каждого горизонта разреза в пределах оцененной вариации значений (табл. 2) при средних параметрах по другим горизонтам. При этом для верхнего горизонта А вариация параметров рассматривалась не только в зависимости от литологического состава отложений, но и от типа ландшафтных условий на поверхности.
Результаты моделирования показывают, что вариация расчетных параметров влагопереноса относительно средних значений несущественно сказывается на итоговых величинах среднемноголетнего инфильтрационного питания подземных вод — в большинстве случаев до 5 мм/год (менее 5-10%)
Таблица 3
Результаты моделирования инфильтрационного питания при вариации расчетных параметров влагопереноса в зоне аэрации
Тип ландшафта Тип пород (почвы) Горизонт разреза Расчетная инфильтрация, мм/год, при значениях параметров влагопереноса
средних минимальных максимальных
Лес Суглинки А 69 76 62
В 79 60
С 70 69
Супеси А 113 119 111
В 119 111
С 116 113
Пески А 147 149 148
В 150 148
С 151 157
Луг Суглинки А 2 6 1
В 6 0
С 4 8
Супеси А 44 47 50
В 43 50
С 47 44
Пески А 102 103 105
В 104 111
С 105 105
Пашня Суглинки А 62 67 57
В 69 56
С 60 62
Супеси А 82 81 86
В 86 82
С 85 82
Пески А 86 88 84
В 89 90
С 90 93
(табл. 3). Это подтверждает возможность использовать средние параметры кривых ОГХ и влагопроводности в модели влагопереноса для расчета инфильтрационного питания подземных вод в соответствии с принятой схематизацией разреза зоны аэрации и типизацией литологического состава пород.
Выводы. 1. На основе анализа изменчивости водно-физических свойств почвенного покрова и подстилающих коренных пород с глубиной разрез зоны аэрации схематизирован тремя расчетными интервалами, отвечающими основным генетическим почвенным горизонтам А, В и С.
2. Количественный анализ гидрофизических параметров влагопереноса для почвенных разрезов, отличающихся литологическим составом отложений, показывает, что кривые ОГХ и влагопроводности практически не зависят от генетического типа почвы, а в гораздо большей степени определяются литоло-гическим составом почвообразующих пород. При этом для верхнего интервала разреза (горизонт А)
-глина -суглинки---супесь .........песок
□ Loamy Sand х Loam й Sandy Clay ■ Образец 2;
глубина I м
Рис. 6. Расчетные кривые ОГХ и влагопроводности для третьего интервала разреза зоны аэрации (горизонт С) при разном литоло-гическом составе пород (линии) в сравнении с типовыми кривыми из базы данных USDA (точки)
расчетные параметры влагопереноса зависят также от типа ландшафтных условий на поверхности земли за счет процессов разуплотнения и преобразования почвенного покрова, связанных с сельскохозяйственной обработкой земель, деятельностью корневой системы растений и живых организмов. За счет этого характер кривых ОГХ и влагопроводности для верхнего почвенного интервала отличается от характерных зависимостей для соответствующего литологического типа пород. Для более глубоких интервалов разреза зоны аэрации ландшафтные условия и характер растительности практически не сказываются на параметрах влагопереноса и форма кривых ОГХ и влагопровод-ности определяется только литологическим составом отложений.
3. На основе проведенной для различных почвенных разрезов оценки параметров влагопереноса, отвечающих аппроксимации кривых ОГХ и влагопро-водности уравнениями Ван Генухтена, обоснованы средние значения расчетных параметров модели вла-гопереноса, характеризующие основные интервалы
1.Е- 01-
1.Е- 02-
св
8
X <и 1.Е- 03-
& -
1.Е- 04-
а
и щ 1.Е- 05-
35
I 1.Е- 06-
5 1.Е- 07-
Пески
0,2 0,3 0,4 Влажность
0,2 0,3 0,4 Влажность
0,1 0,2 0,3 0,4 Влажность
1.Е+00
0,2 0,3 0,4 Влажность
0,2 0,3 0,4 Влажность
0,1 0,2 0,3 0,4 Влажность
0,5 0,6
0,5 0,6
1 А 2 г —о — ■ 4
—■—
Рис. 7. Результирующие кривые ОГХ и влагопроводности для выделенных интервалов разреза зоны аэрации в зависимости от литологического состава пород и ландшафтных
условий на поверхности: 1—3 — горизонт А (1 — луг, 2 — пашня, 3 — лес); 4 — горизонт В; 5 — горизонт С
разреза зоны аэрации в соответствии с типом лито-логического состава отложений, и оценен масштаб их вариации. Результаты тестового моделирования инфильтрационного питания показывают, что изменчивость расчетных параметров влагопереноса в диапазоне их возможной вариации незначительно сказывается на итоговых величинах среднемноголет-ней инфильтрации. Это позволяет использовать полученные параметры модели влагопереноса, характеризующие основные интервалы разреза зоны аэрации и литологический состав пород, для моделирования процессов формирования инфильтрационного питания подземных вод.
4. Схематизация строения зоны аэрации и параметров влагопереноса проведена преимущественно на основе анализа почвенных разрезов, типичных для Нечерноземной зоны России, где наиболее распространены почвы дерновые с разной степенью оподзоленности и серые лесные почвы [Почвенно-геологические..., 1984]. Можно предположить, что выявленные закономерности в целом справедливы и для других регионов, однако в связи с малой представительностью проанализированных разрезов других типов почв эта предпосылка нуждается в дальнейшем обосновании.
СПИСОК ЛИТЕРАТУРЫ
Бадов В.В., Киселев А.А. Капилляриметр: Авт. свид-во № 1354068 СССР, МКИ 01 15.08. БИ № 43. 23.11.1987.
Бадов В.В., Киселев А.А., Романова Т.П. Методы экспериментального изучения влагопереноса в ненасыщенных горных породах // Обзор ВИЭМС. Сер. Гидрогеология и инж. геология. М., 1988.
Вериго С.А., Разумова Л.А. Почвенная влага. Л.: Гидро-метеоиздат, 1973. 328 с.
Глобус А.М. Почвенно-гидрофизическое обеспечение агроэкологических математических моделей. Л.: Гидроме-теоиздат, 1987. 427 с.
Гриневский С.О., Поздняков С.П. Принципы региональной оценки инфильтрационного питания подземных вод на основе геогидрологических моделей // Водн. ресурсы. 2010. Т. 37, № 5. С. 1-15.
Гусев Е.М., Насонова О.Н. Проблемы изучения и моделирования тепло- и влагообмена в системе почва-растительный и (или) снежный покров-приземный слой атмосферы // Вод. ресурсы. 2004. Т. 31, № 2. С. 148-164.
Гусев Е.М., Штекауерова В., Стеглова К. и др. Определение гидрофизических характеристик почв водосбора Иваньковского водохранилища // Водн. ресурсы. 2008. Т. 35, № 3. С. 348-357.
Кучмент Л.С., Гельфан А.Н. Динамико-стохастические модели формирования речного стока. М.: Наука, 1993. 101 с.
Пашковский И.С. Разработка геофильтрационных моделей системы «зона аэрации—грунтовые воды» и их применение при изучении взаимосвязи подземных и поверхностных вод: Автореф. докт. дис. М., 1985.
Полевые методы гидрогеологических, инженерно-геологических, геокриологических, инженерно-геофизических и эколого-геологических исследований / Под ред. В.А. Ко-
ролева, Г.И. Гордеевой, С.О. Гриневского, В.А. Богословского. М.: Изд-во Моск. ун-та, 2000. 352 с.
Почвенно-геологические условия Нечерноземья / Под ред. Е.М. Сергеева. М.: Изд-во Моск. ун-та, 1984. 608 с.
Роде А.А. Основы учения о почвенной влаге. Т. 1. Водные свойства почв и передвижение почвенной влаги. Л.: Гидрометеоиздат, 1965. 663 с.
Субботин А.И. Сток талых и дождевых вод (по экспериментальным данным). М.: Гидрометеоиздат, 1966. 374 с.
Федоров С.Ф. Исследование элементов водного баланса в лесной зоне европейской территории СССР. Л.: Гидро-метеоиздат, 1977. 264 с.
Шестаков В.М., Поздняков С.П. Геогидрология. М.: Академкнига, 2003. 173 с.
Campbell G.S. Soil Physics with BASIC: Transport models for soil-plant systems. Amsterdam: Elsevier, 1985. 150 p.
Chanzy A., Mumen M., Richard G. Accuracy of top soil moisture simulation using a mechanistic model with limited soil characterization // Water Res. 2008. Vol. 44 (W03432, doi:10.1029/2006WR005765, 2008).
Genuchten M.Th. van. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils // Soil Sci. Soc. Amer. J. 1980. Vol. 44. P. 892-898.
Schaap M.G., Leij F.L., Genuchten M.Th. van. Rosetta: A computer program for estimating soil hydraulic parameters with hierarchical pedotransfer functions // J. of Hydrology. 2001. Vol. 251. P. 163-176.
Simünek J., Genuchten M.Th. van, SejnaM. The HYDRUS-1D software package for simulating the one-dimensional movement of water, heat, and multiple solutes in variably-saturated media. Ver. 3.0. // Prepr. Depart. of Environ. Sci. University of California Riverside. California, Riverside, 2005. 270 р.
Поступила в редакцию 15.12.2009
