CC BY
81
Гидравлика, гидрология, водные ресурсы
УДК 502/504 : 556.18
В. Ф. ЖАБИН, Д. В. КОЗЛОВ, Л. Д. РАТКОВИЧ, Е. В. ФРИЗЕНА
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Московский государственный университет природообустройства»
ЭКОЛОГИЧЕСКИЕ ОСОБЕННОСТИ СОВМЕСТНОГО ИСПОЛЬЗОВАНИЯ ПОВЕРХНОСТНЫХ И ПОДЗЕМНЫХ ВОД
Дан анализ взаимодействия подземных и поверхностных вод. Приведены составляющие баланса и параметры потока подземных вод. Рассмотрены вопросы питания и расходования подземных вод, испарения с поверхности грунтовых вод, сопротивления в зоне контакта поверхностных и подземных вод. Определены факторы защищенности подземных вод.
Баланс подземных вод, питание и разгрузка подземных вод, коэффициент фильтрации, единичный расход инфильтрационного потока грунтовых вод, ущерб речному стоку, защищенность подземных вод.
There is considered an interaction of ground and surface water. Components of balance and parameters of ground water flow are given. Questions of ground water feeding and discharge,evaporations from the ground water surface,resistance in the contact zone of surface and ground water are considered. There are defined factors of ground water protectability.
Balance of ground water,feeding and discharge of ground water,filtration factor,unit discharge of ground water infiltration flow, damage to a river flow, ground water protectability.
Совместное нахождение поверхност- задачи на основе анализа влияния изме-ных и подземных вод в верхней части нения экологических условий водосбора земной коры определяет тесную взаимо- на режим и параметры взаимосвязи по-зависимость их балансов, отражающих верхностных и подземных вод. Анализ изменение количественных и качествен- экологической ситуации, обусловленной ных составляющих. Процессы притока П изменением режима подземных вод, пред-и оттока О подземных вод (ПВ) связывают полагает оценку подтопления территории области питания и разгрузки их в есте- их распространения, изменения их ресур-ственных и в нарушенных хозяйственной сов (объема и качества), а также условий деятельностью условиях. Приток подзем- взаимосвязи с поверхностными водами. ных вод во многом обусловлен питанием Взаимодействие подземных и по-
их в пределах водосборного бассейна (во- верхностных вод. Закономерной являет-досбора). Площади водосборов (тыс. км2) ся взаимосвязь атмосферных процессов, принято считать так: малых рек < 1; поверхностных и подземных вод в гло-средних 1...10; больших >10. бальном цикле водообмена. Баланс под-
Изучение взаимосвязи поверхност- земных вод целесообразно рассматри-ных и подземных вод позволяет решать вать в границах симметричного потока, большие экологические и хозяйственные направленного от водораздела к дрене
(рис.1). Морфологический водораздел не всегда совпадает с водоразделом подземных вод, который определяется наивысшим положением свободной поверхности в этом сечении h = h и условием
1 max ^
(dq/dx)x1= 0. Составляющие водного баланса должны определяться независимо друг от друга, как в естественных условиях, так и при нарушении последних в результате хозяйственной деятельности. Анализ стационарного состояния выполняют для среднемноголетних значений составляющих водного баланса. Измене -ние составляющей баланса приводит к его нарушению, хотя разбалансированная гидрогеохимическая природная система стремится к равновесному (устойчивому) состоянию в новых условиях. Для восстановления или улучшения прежнего состояния необходимы мероприятия по регулированию гидрогеохимических процессов с учетом динамики антропогенного воздействия.
♦ и, Тр
AL
Основное питание подземных вод происходит при инфильтрации благодаря нисходящему движению воды от поверхности земли к их зеркалу. Разгрузка осуществляется в виде родникового стока, подземного стока в гидрографическую сеть, испарения со свободной поверхности, а также в виде перетекания в нижележащий водоносный слой.
Различное соотношение составляющих водного баланса на поверхности земли и в зоне неполного водонасыщения (зоне аэрации) определяет величину питания (расходования) грунтовых вод через капиллярную кайму. Результирующее питание, выраженное в миллиметрах, принято называть инфильтрационным w, а расходование - испарением е. Интенсивность процессов ю и е измеряют в метрах в сутки - как правило, это постоянная величина в течение расчетного периода времени. Коэффициент инфильтрации п в определенных условиях может быть равен коэффициенту подземного стока. При неизменной норме испаряемости коэффициент инфильтрации оценивают по формуле [1]:
V Зона разгрузки ГВ
«И
Рис. 1. Составляющие баланса и параметры потока подземных вод: И, Тр - испарение, транспирация соответственно; Ос
- атмосферные осадки; w (д),е (-д) - соответственно питание и испарение грунтовых вод; р, -р - составляющая взаимосвязи грунтовых и напорных вод (питание и разгрузка); де - естественный приток грунтовых вод (динамические запасы); АЬ
- дополнительное фильтрационное сопротивление в зоне разгрузки грунтовых вод; К , К - горизонтальная и вертикальная составляющие водопроницаемости водоносных фильтрационно-анизотропных горных пород, коэффициенты фильтрации; К2, т2 - коэффициент фильтрации и мощность относительного водоупора между грунтовым и напорным потоками подземных вод, со свободной и пьезометрической поверхностями соответственно; I - длина потока грунтовых вод.
n = n0(i+);
N
(1)
п0 =
N
где п0 — среднемноголетний коэффициент инфильтрации при норме осадков N0 и слое инфильтра-ционного питания w0; п1 - то же в изменившихся условиях по водности при норме осадков N1.
Отношение величин атмосферных осадков и испаряемости Е0 за один и тот же промежуток времени определяет водный режим зоны неполного насыщения и условия питания грунтовых вод. Эта предпосылка и анализ экспериментального материала по различным климатическим областям позволили получить следующие соотношения, справедливые для равнинных территорий и нормального просачивания.
N0/E0 £ 0,1 0,6 0,7 1,0
п0 - 0,07 0,085 0,09 0,11
N0/E0 1,2 1,3 1,4 -
0,13 0,16 0,2 -
Испарение с поверхности грунтовых вод (далее - испарение) определяет наряду с инфильтрацией направленность влагообмена в зоне аэрации, влияет на
почвообразовательный процесс, повышает содержания солей в грунтах и грунтовых водах. Испарение может равняться доле естественных ресурсов грунтовых вод и иногда представлять собой основную расходную статью их баланса, обеспечивая стабилизацию свободной поверхности. Этот фактор значим при анализе гидрохимического режима зоны аэрации и подземных вод.
Исследование связи испарения с глубиной грунтовых вод и свойствами грунтов ненасыщенной зоны [1] привело к зависимости
_ "Е - Е,
(2)
^ , К:
Е
z
■ Е1( ), 4аГ
если -< 0,05, Я = 4,6 ^
4at z
где у - коэффициент пропорциональности (экспериментальный); е. - испарение грунтовых вод; Е. , Е0 - испарение с открытой водной поверхности за месяц и год соответственно; Е1 - интегральная экпоненциальная функция; а = 1 /(с о) - коэффициент температуропроводности, м2/сут (1 - коэффициент теплопроводности, ккал/(м . сут х ° С); с - теплоемкость, ккал/(кг • ° С); о - объемный вес грунта, кг/м3; t - общая продолжительность процесса, сут; z - средняя мощность зоны аэрации, м.
Теплофизические характеристики грунтов зоны неполного насыщения зависят от их состава, влажности, а также от морфогенетических особенностей и структуры порового пространства. Заметное влияние оказывает минерализация грунтовых вод. Для ориентировочных расчетов при песчано-глинистом составе грунтов зоны аэрации значения коэффициента температуропроводности можно выбирать в диапазоне (1...2) • 10-2м2/сут (у0 < 0,08), затем уточнять выбранные характеристики по результатам эксперимента.
Практически полное прекращение испарения наступает при достижении поверхностью грунтовых вод «критической» глубины, которую можно отождествлять с климатической характеристикой региона гк (по аналогии с глубиной промерзания). Для ее расчета при определении годовой величины испарения по формуле (2) рекомендуется использовать предложение В. А. Ковды:
=(170 см + 8ГС) = фл/о7, (3)
где t - среднегодоваятемпературавоздухав регионе, С.
Определение соотношений напоров по всей длине потока, как в горизонте грунтовых вод, так и в находящихся под ним пластах, необходимо при оценке составляющих баланса подземных вод.
Если по пути движения потока грунтовых вод не происходит их питания или расходования, то в этом случае в уравнении общего водного баланса П = О. Во всех остальных случаях величина оттока подземных (грунтовых) вод будет отличаться от притока к границе балансового участка на величину, равную разгрузке (родниковой, испарению) или питанию (инфильтрационному, напорному) в пределах участка. В дальнейших рассуждениях будем использовать интегральную величину питания грунтовых вод с учетом ее изменение по длине потока.
При анализе взаимодействия поверхностных и подземных вод большое значение имеет сопротивление в зоне их контакта АЬ, прямо пропорциональное потерям напора в пределах зоны кольма-тации («дополнительный слой») [2]. Величину сопротивления принято выражать в метрах эквивалентной дополнительной длины потока, полученной расчетным путем или графическим построением [3, 4].
Формулы для расчета этого параметра и условия их применения приведены в табл. 1.
Аналитические способы определения сопротивления дрен
Таблица 1
0
2
2
Расчетная схема Строение пласта Расчетная формула АХ
Водоем или водоток шириной 2Ьг Однородно-анизотропный АЬу~0,5атй, а = \11Г (4) V К
Двухслойный [3, 4] АЬ. = АЬ2с1Ц^) , АЬ2 =
Береговая дрена, по С. Ф. Аверьянову Однородно-анизотропный кЛ= г = (6) 1+ П<1 8Ш хА +АЬ 2ат&
В формуле (6), по которой определяется общее сопротивление береговой дрены Ьы, перед тА добавляется множитель, учитывающий фильтрационную анизотропию, как в формуле (4).
Из табл. 1 ясно, что практически все дренирующие границы потоков грунтовых вод имеют дополнительное фильтрационное сопротивление, которое необходимо учитывать путем замены действительной координаты х на расчетную: X = х + АЬ.
Заметим, что в формулу (5) входит половина ширины естественной дрены, которая может быть определена так: Ьг > 3АЬ2 для анизотропного двухслойного строения ложа водоема или водотока.
AL _. —. - -
7
Л d
1 + Lbd
при L^ = AL, (*d = 0) => L^ = 0,5 и 0,67q0.
При этом на урезе разгружается qro = 0,33q0. Иными словами, независимо от величины сопротивления на урезе
водотока (водоема) разгружается около 1/3 расхода потока ГВ, идущего с водораздела. Близкий по значению результат получается, если использовать решение В. М. Шестакова с заменой х на X [4].
В этом случае расход грунтовых вод в каждом сечении:
Д ГГ х+АЬ х+АЬ
= Т —0 в^ = дп в^; (9)
под рекой - дх
DL
Рис. 2. Береговая дрена, ряд скважин (расчетная схема): А Ь - дополнительное фильтрационное сопротивление, м; -расстояние от уреза воды до центра скважины, м; d - диаметр скважины, м; д0 -единичный установившийся расход, м2/сут
Единичный установившийся расход, т.е. удельная величина динамических запасов в условиях естественной дрениро-ванности
"о = ю1. (7)
Однако расход воды, поступающий в дрену, распределится в зависимости от ее несовершенства и условий формирования потока подземных вод. Так, единичный расход инфильтрационного потока грунтовых вод дг, поступающего в реку, следуя рекомендациям С. Ф. Аверьянова, определим так:
на урезе - х = 0 и qro = д0 е"1 « 0,368д0.
Тот же результат получается при решении дифференциального уравнения: -1±I = л_; _ ^ = х + х = 0; ^ = (10)
йХ АЬ 4 АЬ
а X х+аь
— 1п~=АЬ' ах = аов АЬ; ах=о« а36^- (11) а0 АЬ
Определим ширину зоны разгрузки в реке, считая от уреза. При х = 2А Ь и qr = 0,05q0 поток грунтовых вод разгрузится только на 95 %. При меньшей ширине реки «проскок» потока грунтовых вод под рекой будет больше, на что необходимо обращать внимание при проектировании береговых скважинных водозаборов в зонах возможного загрязнения грунтовых вод на противоположном берегу.
Аналитические расчеты дополнительного сопротивления затруднительны, поэтому наиболее пригодными для практического использования являются данные измерений положения свободной поверхности стационарного потока грунтовых вод вблизи зоны разгрузки.
По рекомендации В. М. Шестакова, когда инфильтрация отсутствует, то
Н, — Н р
АЬ = —-р • х2 — х,. (12)
Н 2 — Н, 2 1
При инфильтрации к величинам Н , Н2 добавляются слагаемые юх,2 ¡2Т юх
(8)
Рис. 3. К расчету дополнительного фильтрационного сопротивления
Экспериментальные методы определения АЬ представлены в [4].
Для горных районов наряду с вариантами питания водотоков характерны случаи потерь речного стока
№ 5' 2011
(бэ)
характерны. В долинах горных рек формирование подземного стока в верхнем и среднем течениях происходит в основном за счет потерь из русла водотока П, в то время как в нижнем наблюдается увеличение стока реки за счет разгрузки подземных вод. Величину одностороннего удельного водообмена ПВ и реки можно определить по разнице расходов реки в двух створах, расположенных на расстоянии У, км:
43,2(б2 - 01)
Ч.
У
(13)
Для замкнутого экранированного водоема И. С. Пашковский рекомендует следующую формулу:
Л (Ьг0)
0 = 2 рТАИЬг0-^-^-, (14)
Л0 (ЬГ) )
где - расход на перетекание, м3/сут; Т - коэффициент водопроводимости; АН - превышение пьезометрического уровня водоносного горизонта над поверхностью водоема, м; г0 - приведенный радиус водоема, м; 11(Ьг0), 10(Ьг0) - модифицированные функции Бесселя первого и нулевого порядков соответственно; Ь2 = ^ к1/(Т2т1).
При работе водозабора, когда общий дебит < складывается из двух составляющих - речных вод (привлекаемых ресурсов) и притока ПВ со стороны водораздела (динамических запасов), Ф. М. Бочевер полагает [5], что _
26
У
0г = — агс*§ - 2ЧоУо, % = X
1.(15)
р X ^РХд0
Половина зоны захвата естественно-
го потока подземных вод:
О X ^ -^аг^—+ У).
В = ■
У
(16)
й И о Он
ч
А
Хо '
V
Рис. 4. Скважинный водозабор у реки
Например, если дебит водозабора <г составляет 5186 м3/сут, X = 522 м, У0 = 1047 м, то из реки поступает <р = 2561 м3/сут. Со стороны водораздела при q0 = 0,525 м2/сут поступает примерно 2625 м3/сут. При размерах зоны захвата 2В ~ 5000 м.
При отсутствии естественного притока производительность водозабора обеспечивается целиком за счет фильтрации из реки. При этом 80...85 % дебита водозабора обеспечиваются на длине 2У0~ 6Х за счет притока из реки. Соотношение расходов позволяет прогнозировать качество воды, направляемой потребителю.
Ущерб речному стоку при работе водозабора подземных вод теоретически будет всегда, однако величина его зависит от строения, питания и гидрогеологических параметров системы подземных вод, расстояния от уреза реки, величины сопротивления в зоне разгрузки ДЬ и времени эксплуатации t. Соответствующие расчеты подробно представлены М. М. Черепанским [2]. Для инженерной оценки этого фактора можно использовать уравнение М. Хантуша:
О - - Ь+АЬ
- ег&Я ; Я=—=■ , (17),
где Ь - расстояние водозабора от реки, м; а - коэффициент пьезопроводности (уровнепроводности), м2/с.
Справедливо считать, что
- О аО
0 = Ог = 0Г '
где а - доля статических запасов, не связанных с поверхностным стоком.
При уменьшении испарения с поверхности грунтовых вод величина а увеличивается на расчетную величину инверсии испарения р. Тогда допустимо полагать,_что
а + р- егЖ. (18)
В данном соотношении критерием сравнения целесообразно выбирать зимний меженный сток реки обеспеченностью 50 % и более.
В табл. 2 приведены результаты расчета потерь речного стока при отборе подземных вод в Калужской области РФ с водосбора площадью 24,56 тыс. км2. Здесь при работе 700 водозаборов, эксплуатирующих напорный водоносный комплекс в известняках нижнего карбона, образовалась региональная депрессионная воронка диаметром 20 км.
В таблице использованы данные Геоцентра (Москва) и ЗАО Совинтервод. Величина водоотдачи трещиноватых известняков принята 0,05, понижение при отборе подземных вод составляет при этом 30 м. Зная значение водопроводи-мости, можно оценить «обобщенную»
(во)
№ 5' 2011
Таблица 2
Оценка статической составляющей запасов подземных вод (ПВ) и ущерба речному стоку (округлено)
Водные ресурсы, тыс. м3/сут Время отбора, сут Понижение, м Статические запасы ^ст ^отб ПЭЗ, связанные с рекой Ущерб речному стоку, %
Отбор ПВ (19982008), Я с ^ото Прогнозные эксплуатационные запасы (ПЭЗ) Меженный сток рек V, м3 б , тыс. м3/ сут
250 3612 3601 7300 30 1,6.108 22 0,09 3287 91
водоотдачу, учитывающую изменение характеристик пласта, вертикальный переток и другие факторы, проявившиеся при эксплуатации водоносной системы [5]. Обобщенный модуль подземного стока близок к 1,5 л/с . км2.
Заметим, что определение ущерба речному стоку в сопоставлении с прогнозными эксплуатационными запасами имеет в большей степени теоретическое значение. При фактическом отборе подземных вод (250 тыс. м3/сут) этот показатель составляет около 6 %. В случае полного использования эксплуатационных запасов (877 тыс. м3/сут) потери речного стока равны 22 %.
Динамика загрязнения подземных вод. Очевидно, что экологическая обстановка на земной поверхности, проникновение вредных веществ в первый водоносный горизонт и дальнейшее распространение загрязнения непосредственно связаны с условиями питания и разгрузки грунтовых вод. Защищенность подземных вод от загрязнения определяется наличием слабоводопроницаемых отложений, препятствующих проникновению загрязняющих веществ с поверхности земли в водоносные системы. Факторы, определяющие защищенность подземных вод, можно разбить на три группы - природные, техногенные, физико-химические.
К основным природным факторам относятся следующие: наличие в разрезе слабоводопроницаемых отложений; глубина залегания подземных вод; состав, фильтрационные и сорбционные свойства перекрывающих пород; соотношение уровней исследуемого горизонта и вышележащего (верховодки).
Группа факторов, контролирующих загрязнение подземных вод, включает такие свойства загрязнителей, как их миграционная способность, плотность, сорбируемость, химическая консервативность или время распада (разложения),
№ 5' 2011
взаимодействие загрязняющих веществ с горными породами и подземными водами.
К техногенным факторам следует отнести условия накопления, хранения, сброса загрязняющих веществ, а также их очистку, утилизацию и транспортировку в штатной и аварийной ситуациях. Исследования в этой области показывают, что ликвидация любого начавшегося загрязнения поверхностных вод - весьма сложная в техническом отношении, трудоемкая и дорогостоящая задача. В отдельных случаях загрязнение грунтовых вод оказывается необратимым. В этой связи выполнение прогнозных оценок по загрязнению с предварительной реализацией специальных гидрогеологических исследований, а также профилактические (предупредительные) мероприятия оказываются более выгодными в технико-экономическом отношении, чем борьба с последствиями загрязнения подземных вод.
Основным агентом переноса химических веществ и бактерий ниже поверхности земли является вода, с которой, как и с водовмещающими породами, эти вещества взаимодействуют, меняют химический состав и концентрацию. Миграция и трансформация химических соединений происходит в результате процессов молекулярной диффузии, движения воды, дисперсии, сорбции, ионного обмена, растворения и кристаллизации [1, 3, 4].
Биохимический и температурный обмен между речными и подземными водами происходит в зоне гипорейка [4], включающей частично или полностью «дополнительный слой». Мощность зоны зависит от строения и состава аллювия, скорости и направленности водообмена между рекой и подземными водами. Латеральные размеры зоны гипорейка увеличиваются при образовании стариц и меандр. Колебания размеров тесно связаны с факторами формирования и существования зоны, и антропогенное
(б1)
вмешательство в процесс взаимодействия речных и подземных вод может привести к сокращению размеров зоны или к ее исчезновению.
Процессы, происходящие в зоне сопряжения поверхностных и подземных вод, требуют постоянного мониторинга как основы планирования совместного использования водных ресурсов [1].
Выводы
Отбор подземных вод приводит к уменьшению поверхностного стока, что особенно сказывается на состоянии малых и средних рек в период межени. Ущерб речному стоку необходимо учитывать при оценке эксплуатационных запасов подземных вод.
Запасы подземных вод, не связанные с поверхностным стоком, соответствуют статическим и упругим запасам водоносных систем, а также потерям подземных вод при испарении с их поверхности.
Разгрузка подземных вод в реки зависит от сопротивления ложа водотока (водоема), являющегося его характеристикой (дополнительный слой).
Качественный водообмен поверхностных и подземных вод происходит в контактной зоне гипорейка, мониторинг которой необходим при использовании водных ресурсов.
1. Манукьян Д. А., Жабин В. Ф. Гидрогеоэкологические проблемы в задачах
природ ообустройства. - М.: ФГОУ ВПО МГУП, 2006. - 194 с.
2. Черепанский М. М. Теоретические основы гидрогеологических прогнозов влияния отбора подземных вод на речной сток. - М.: НИА-ПРИРОДА, 2005. - 260 с.
3. Шестаков В. М. Теоретические основы оценки подпора, водопонижения и дренажа. - М.: Изд-во МГУ, 1965. - 233 с.
4. Шестаков В. М., Невечеря И. К., Авилина И. В. Методика оценки ресурсов подземных вод на участках береговых водозаборов - М.: Изд-во Книжный дом «Университет», 2009. - 192 с.
5. Бочевер Ф. М. Расчеты эксплуатационных запасов подземных вод. - М.: Недра, 1968. - 325 с.
Материал поступил в редакцию 04.04.11. Жабин Виктор Федорович, кандидат геолого-минералогических наук, доцент Тел. 8 (495) 312-34-70
Козлов Дмитрий Вячеславович, доктор
технических наук, профессор, ректор
E-mail: kozlovdv@mail.ru
Раткович Лев Данилович, кандидат
технических наук, профессор
Тел. 8 (499) 976-21-56
E-mail: levkivr@mail.ru
Фризена Елена Владимировна, магистр
Тел. 8 (909) 647-08-20
E-mail: elenafrizena87@mail.ru
