05.23.07 Гидротехническое строительство
УДК 502/504 : 627.8
А. И. ГОЛОВАНОВ, Ю. И. СУХАРЕВ
Федеральное государственное бюджетное образовательное учреждение высшего образования «Российский государственный аграрный университет — МСХА имени К. А. Тимирязева» Институт природообустройства имени А. Н. Костякова
С. В. СОЛЬСКИЙ
Открытое акционерное общество «Всероссийский научно-исследовательский институт имени Б. Е. Веденеева»
МЕТОД РАСЧЕТА ФИЛЬТРАЦИОННЫХ КУПОЛОВ
ПОД ДНИЩЕМ БАССЕЙНОВ СУТОЧНОГО РЕГУЛИРОВАНИЯ ГЭС
Предложен способ расчета фильтрационных куполов под днищем бассейнов суточного регулирования (БСР) деривационных гидроэлектростанций. Рассмотрены вопросы схематизации, граничные условия, алгоритм решения уравнения движения грунтовых вод, результаты расчетов с различными параметрами пласта и граничными условиями. Оценено вакуумирование дренажа в основании БСР. Дан анализ результатов расчетов.
Бассейны суточного регулирования деривационных ГЭС, фильтрационные купола, трехмерная математическая модель.
There is proposed a method of calculation of filtration cupolas under the basins bottom of daily regulation (BDR) of diversion hydropower plants. There are considered questions of schematization,boundary conditions,algorithm of equation solution of ground water movement, results of calculations with different parameters of layer and boundary conditions. Drainage vacuumizing is assessed in the basis of BDR. There is given an analysis of calculations results.
Basins of daily regulation of diversion HPP, filtration cupolas, three-dimensional mathematical model.
Для разработки диагностических устройств с целью локализации трещин в днищах бассейнов суточного регулирования (БСР) на деривационных гидроэлектростанциях необходимо знать параметры фильтрационных куполов под днищем БСР. Эта фильтрационная задача затруднительна из-за сложной конфигурации области фильтрации, явно выраженной многомерности потока и наличия зоны неполного водонасыщения. Аналитического решения такая задача не имеет, но современные методы описания фильтрационных потоков позволяют это сделать, применив математическое моделирование.
Условия моделирования были следующие: размеры бассейна суточного регулирования в плане - 20...50 м2, днище мощностью 0,3.0,5 м из монолитного
железобетона и уложенного по слою геотекстиля на песчаной подушке. Нижний слой - щебень, примерно такой же мощности, выполняющий роль пластового дренажа для снятия фильтрационного напора и взвешивающих сил, которые при переменном и быстро изменяющемся уровне воды в бассейне могут вызвать нежелательные знакопеременные нагрузки, приводящие к образованию и развитию трещин в теле днища. Для ликвидации этого инфильтрационные воды должны отводиться дренажем, который с целью повышения его эффективности требует вакуумирования. Размеры этих трещин могут составлять несколько миллиметров.
Помимо этого возможно нарушение температурно-деформационных швов и порывы гидроизоляции, образующие свищи
и инфильтрацию в дренажную подушку. Размеры таких отверстий могут измеряться сантиметрами. Авторы считают, что таких источников инфильтрации, а следовательно, фильтрационных куполов, может быть несколько на один бассейн, они могут интерферировать, взаимодействуя друг с другом.
В связи с этим область фильтрации можно представить в виде пласта-бочки, приходящейся на одно отверстие, расположенное в центре. Пласт имеет изолированную боковую поверхность, что позволяет учесть интерференцию куполов. Верхней поверхностью пласта является днище бассейна суточного регулирования, нижней - непроницаемый или слабопроницаемый слой грунта основания, на котором лежат песчаная и щебеночная прослойки. Движение в таком пласте будет явно трехмерным, но задача упрощается, если движение представить радиально осевым. Фильтрационное тело принято многослойным с существенным отличием водно-физических свойств слоев.
На верхней поверхности расчетной схемы необходимо задать расположение отверстия и рассмотреть для общности два случая: когда в отверстие поступает поток воды с заданной интенсивностью; когда над отверстием формируется постоянный напор воды. В модели предусмотрена изменчивость во времени этого напора. На остальной площади, перекрытой фундаментной плитой, отсутствуют как атмосферные осадки, так и испарение (конденсацией влаги пренебрегаем).
Нижняя граница может быть представлена в двух вариантах: уровень грунтовых вод (УГВ) в дренируемой толще постоянный или переменный, когда происходит заполнение водой дренируемой толщи (в этом случае нет оттока фильтрата). Первое условие реализуется, когда в дренажной системе поддерживается постоянный напор, который с помощью сифона может образовывать некоторый регулируемый вакуум.
В значительной части фильтрационного тела движение воды будет происходить при неполном влагонасыщении, что требует учета капиллярных сил. Действующий напор надо представлять в следующем виде [1]: Н = -х + у ,
где Н - действующий напор (отсчет - от
fee)
поверхности днища и оси x, направленной вниз), м; x - гравитационная составляющая напора, м; у - напор, эквивалентный каркасно-капиллярному давлению в зоне неполного насыщения (у <= 0) и эквивалентный гидростатическому давлению в зоне полного насыщения, м.
Такая обобщенная трактовка у позволяет описывать передвижение влаги в слое, включающем в общем случае изменяющийся в нем уровень грунтовых вод, а также зоны полного и неполного влагонасыщения пористой среды. Для зоны неполного влагонасыщения нужно иметь зависимость между влажностью ю и у, которую называют основной гидрофизической характеристикой пористой среды (ОГХ). Она может быть определена по данным полевых экспериментов, например, при одновременных измерениях влажности и напора пористой среды с помощью тензиометров. Авторами эта зависимость принята в следующем виде [2]:
Г L..I \п~
— СО
со = —
со.
тп- ю.
= ехр
А
где йк - максимальная высота капиллярного поднятия, м; <вм - максимальная гигроскопичность, мв3/м3; т - влажность при полном заполнении пор, численно равная пористости, мв3/м3; V и п - эмпирические безразмерные коэффициенты, подбираемые так, чтобы наилучшим видом описать экспериментально определенную ОГХ (в наших расчетах приняты V = 1 и п = 1).
Вместо коэффициента фильтрации используется величина коэффициента влагопроводности. При изменении влажности в максимально широком диапазоне - от полного насыщения до максимальной гигроскопичности - эта зависимость принята в виде, несколько отличающемся от известной зависимости С. Ф. Аверьянова из-за расширения диапазона влажности грунта:
(
к — к*
ю-CO.,
Л
7П-00
м у
где к{ - коэффициент фильтрации, м/сут.
Вместо обычного коэффициента водоотдачи при неполном влагонасыщении применяют коэффициент влагоемкости а = ды1 дН , характеризующий изменение влажности грунта при небольшом изменении каркасно-капиллярного напора. Коэффициент влагоемкости определяется дифференцированием аппроксимации ОГХ. Для принятой зависимости
5' 2014
О = 71-
т-ш, А
1п
а)
п-1
п
а при п = 1 и V = 1
к
(например, для песка при влажности ю = 0,15; юм = 0,05 и \ = 0,5 м о = 0,2 м-1).
Естественно, что при у > 0 с =0, т. е. под уровнем грунтовых вод (при полном влагонасыщении) коэффициент влагоемкости равен нулю, а водоотдача проявляется при изменении влажности над уровнем грунтовых вод при их колебании. При влажности близкой к максимальной гигроскопичности коэффициент влагоемкости также стремится к нулю.
Запишем уравнение движения грунтовых вод в обеих зонах в случае радиально-осевой двумерной фильтрации, выраженное в полных напорах влаги, в конечно-разностной форме (неявная схема) для блока i, j:
ЯП+1 ГГП ГГП+1 П-П ттп+1 ТТЛ
1,3 ~ДЦ _-НЦ-1~-НЦ -"ц
А* АД"Н НЖ;
гп+1 _ ттп 1,3 "1,3 1+1,3
_п+1 а1,3
■
3 1,3
+
КЙ-Щ
_е1,3'
(1)
где - вертикальное сопротивление потоку влаги между центрами г,у и г,у +1 блоков, сут.
На моменты времени п и п + 1:
где Дг - горизонтальное сопротивление потоку влаги между центрами г, у и г, у +1 блоков, сут/м:
1 Яа>1+и >1+0,5&/
пласта-бочки. В примере площадь на одно отверстие принята равной 30 м2, а радиус 3,09 м.
Определение напоров почвенной влаги Н^1 с помощью системы алгебраических уравнений представляет собой громоздкую вычислительную задачу, так как сводится к нахождению порядка 250 неизвестных (при принятой разбивке на блоки) с малым шагом ввиду сильного возмущения потока. Следует также отметить существенную нелинейность этой системы уравнений, в которой емкостной коэффициент и проводимость существенно зависят от напоров почвенной влаги, следовательно, и от влажности почвы, что требует 3.7 итераций на каждом временном шаге. Поэтому алгоритм решения системы (1) должен быть наиболее эффективным. В настоящее время таковым является метод матричной прогонки, который введением вектора напоров по всем г-м столбцам для каждого слоя у позволяет понизить размерность задачи до одномерной [3]: и = 1; 1; Н+1; ... при
этом у = 0, 1, 2, 3 ...
С помощью этого вектора систему (1) уравнений запишем в матричном виде:
ллрУ1 - . + BBjUj+l = F,
где ЛЛ.. и ВВ. - квадратные диагональные матрицы размером (Мг -1) • (N -1), учитывающие вертикальные потоки влаги между г, у -1 и г, у блоками и между г, у и г,у +1 блоками; СС. -квадратная трехдиагональная матрица размером (Ыг -1) • (Ыг -1), учитывающая потоки влаги между г -1, у' и г, у блоками и между г, у и г + 1у блоками, а также емкостной член; вектор ^ объединяет все свободные члены:
Сш^1 Сш9п+1 Си).
где Дt - расчетный шаг по времени, сут.
Вся исследуемая толща была разбита на элементарные слои И. (1 < у < ) с толщиной 0,05 м, при этом Н0 = Н = 0 . Число слоев - 21, включая и нулевой, т. е.
= 21. Рассматриваемый «пласт-бочка» разбивался вертикальными цилиндрами для образования столбцов и расчетных блоков. Ширина этих блоков (по радиусу) Ь (1 < Ь < Nг - 1) принималась переменной, а Ь0 = Ьг = 0 . Число блоков в одном слое было равно 12, включая и нулевой, т. е. Nг = 21. Центр блока находился на расстоянии г = г-1 + (Ь-1 + Ь{) / 2 от оси «пласта-бочки».
Задавшись площадью днища, приходящейся на одно отверстие, и заменив ее на равновеликий круг, получаем радиус
Aí
|
Дí
Aí
Источники или стоки (если они имеются) на вертикальных границах или внутри области фильтрации учитываем при вычислении этого вектора.
Решение системы матричных уравнений находим в виде рекуррентной формулы:
= рР • . +
Для этого при прямой прогонке вычисляем матрицы прогоночных коэффициентов РР. и прогоночные векторы-столбцы QQ.:
РЦ =11ССj - ^ • рры II1 ВВ.} ;
=|| СС-АА^- 1Г1 • || АА, • ^ - Щ ||.
В этих формулах || а1к ||-1 обозначает
обращенную матрицу.
Матрицы РР. квадратные, размером (Ыг -1) • (Ыг -1) , их общее количество равно при принятой разбивке области фильтрации они содержат 12х21 = 252 элемента. Длина вектора-столбца ЯЯ. 21 элемент, всего таких векторов N .
При обратной прогонке вычисляем искомый вектор напоров почвенной влаги и. на конец временного шага Д£ .
Верхнее граничное условие учитываем, особым образом вычисляя первую матрицу прогоночных коэффициентов РР0 и первый вектор-столбец QQ0. Так, если через верхнюю границу (х = 0) нет потока влаги, то Н1+1 = Я7+1 при заданном напоре —о"1 = Н0 . Элементы диагонали матрицы РР0 равны Ри = 1, а остальные
- нулевые. Все элементы вектора-столбца QQ0 равны нулю. При заданном потоке воды через поверхность грунта интенсивностью q м/сут диагональные элементы матрицы РР0 равны Рц = 1 , а остальные -нулевые. Элементы вектора-столбца QQ0 равны: Qi = -qR^0).
Условия на нижней границе (х = хм ) реализуются при особом вычислении последнего вектора напоров и^:
1) при отсутствии потока влаги через нижнюю границу (водоупор) Н1,,? = Н, поэтому напоры вычисляют по формуле
и^ =1ПЕЕ-РР^ II"1
где ЕЕ - квадратная диагональная единичная матрица, т. е. у которой элементы Е.. = 1;
2) при наличии гидравлического взаимодействия с пластами более глубокого залегания (напорное подпитывание) интенсивностью
а1 - К0 т » о
где к0 - коэффициент фильтрации раздельного слабопроницаемого пласта; Т0 - его мощность; Ър
- отметка пьезометрического уровня в напорном пласте; йп - уровень грунтовых вод в ¿-м столбце, считая от днища БСР,
напоры на нижней границе вычисляем по формуле
и^ =| Н^ |=|| ЕЕ - РР^-11|-1 •(ЯЯ^-1 - ddi),
где dd¡ - вектор-столбец с элементами d¡;
3) при заданной постоянной во времени глубине грунтовых вод Ъ
-К;
4) при очень глубоких грунтовых водах, не участвующих в круговороте почвенных вод,
и^ =|Я£ |=||1 - ррКх
-1 .(QQn-! -10, 5Hn - |),
где |0,5ftNx1| - вектор-столбец с одинаковыми элементами, равными 0,5^Nx1.
Компьютерная программа фильтрационных расчетов написана на алгоритмическом языке Visual Basic, версия 6. Блок-схема обычная: ввод исходных данных, подготовка к расчету, решение системы уравнений, их анализ и формирование файлов результатов расчетов.
По разработанной математической модели выполнена серия расчетов с различными параметрами пласта и граничными условиями. В основании БСР исследовалась песчаная подушка мощностью 0,5 м со свойствами: пористость 0,35; максимальная гигроскопичность 0,05; коэффициент фильтрации 4 м/сут; максимальная высота капиллярного поднятия 0,6 м. Уложенная ниже пластовая дрена толщиной 0,5 м из щебня имела следующие свойства: пористость 0,30; максимальная гигроскопичность 0,04; коэффициент фильтрации 40 м/сут; максимальная высота капиллярного поднятия 0,5 м.
Было принято, что днище бассейна имеет отверстие диаметром 20 мм. Площадь днища, приходящаяся на одно отверстие, равна 30 м2. Поступление воды в отверстие определено по двум расчетным схемам: под действием напора в бассейне или при заданном расходе утечки. Слой щебня дренируется вакуумным дренажем, разряжение в дренажной трубе принято равным 0,5 м водного столба, т. е. напор в нем равен 1,5 м, считая от дна бассейна. Такая схема объясняется требованием эффективно дренировать основание бассейна для снятия взвешивающего гидростатического давления. Следует отметить, что если дренажная труба соединяется с атмосферой, то вода из не полностью влагонасыщенной песчано-щебенчатой подушки вытекать не будет, пока давление в ней не превысит или не станет равным атмосферному, а пористое пространство не будет полностью насыщено водой. При этом под днищем создается противодавление, особенно большое около отверстия, т. е. создаются условия для дальнейшего раскрытия температурно-деформа-ционных швов.
Было важно сравнить результаты расчетов при разных схемах поступления воды в отверстия: под заданным напором или при заданном расходе утечки. Оказалось, что они практически мало отличаются.
Были оценены размеры куполов
растекания, а точнее, увлажнения, под отверстием, для чего использовалась конечная установившаяся матрица влажностей. Эпюры влажности в куполе увлажнения при напоре в отверстии 20 м показаны на рис. 1. На рисунке выделены песчаный и щебеночный слои, отражено действие вакуумного дренажа (малая влажность в слое щебня). Под отверстием формируется зона промачивания, которая не всегда достигает щебеночного слоя. О дальности промачивания можно судить по изоплете влажности 0,10 объема на подошве песчаного слоя (в данном примере около 1,8 м).
0,35^
К 0,30 л
S 0,25
0,001
0,0 0,4 0,8 1,2 1,6 2,0 2,4 2,8 Расстояние от отверстия, м
Рис. 3. Объемная влажность песчаного слоя на его подошве при разных напорах в отверстии: 1 - напор 1м; 2 - напор 5 м; 3 -напор 20 м; 4 - напор 50 м; 5 - напор 100 м
Рис. 1. Установившиеся эпюры влажности при напоре 20 м
Дальность промачивания, или радиус фильтрационного бугра существенно зависит от напора в отверстии (рис. 2). При редко расположенных отверстиях и при больших напорах радиус может доходить до двух метров, а при частом размещении они могут смыкаться. Распределение влажности песчаного слоя на его подошве показано на рис. 3.
(X Напор в отверстии, м
Рис. 2. Дальность увлажнения под отверстием в зависимости от напора в нем
5' 2014
Суточный слой впитывания, отнесенный к площади на одно отверстие, пропорционален напору в нем и при наших условиях колеблется от 1,5 до 150 мм/сут при напорах 1...100 м. Он больше зависит от водопроницаемости песчаной подушки, чем от гидравлического сопротивления самого отверстия.
Качественный анализ результатов расчета показывает, что предложенный способ расчета влияния отверстий в днище бассейнов суточного регулирования деривационных ГЭС достаточно подробно описывает размеры куполов увлажнения и может быть использован для разработки схемы размещения датчиков утечки воды через трещины.
1. Голованов А. И., Новиков О. С. Математическая модель переноса влаги и растворов солей в почвогрунтах на орошаемых землях: Труды МГМИ. - М.: МГМИ, 1974. - Т. 36. - С. 87-95.
2. Плюснин И. И., Голованов А. И. Мелиоративное почвоведение. - М.: Колос, 1983. - 318 с.
3. Голованов А. И., Сухарев Ю. И. Математическая модель влагопереноса в ландшафтных катенах: Природообустройство и рациональное природопользование -необходимые условия социально-экономического развития России: сб. науч. трудов.- М.: МГУП, 2005. - Ч. 2. - С. 3-11.
Материал поступил в редакцию 11.09.14. Голованов Александр Иванович, доктор технических наук, профессор Тел. 8-916-341-35-51 E-mail: a.i.golovanov@mail.ru Сухарев Юрий Иванович, доктор
(4i)
технических наук, профессор кафедры «Мелиорация и рекультивация земель» Тел. 8-916-611-80-23 E-mail: vodoem@mail.ru Сольский Станислав Викторович,
доктор технических наук, заведующий отделом «Основания, грунтовые и подземные сооружения» Тел. 8-921-964-36-05 E-mail: solskiysv@vniig.ru
УДК 502/504 : 626/627 В. Н. АВЕРЬЯНОВ
ЗАО «Проектно-изыскательское научно-исследовательское бюро «ГИТЕСТ»
ИСПОЛЬЗОВАНИЕ НАБУХАЮЩЕЙ ГЛИНЫ ДЛЯ ВОЗВЕДЕНИЯ ПРОТИВОФИЛЬТРАЦИОННЫХ ЭКРАНОВ ГАЭС
Дана оценка возможности использования набухающих глин для возведения противофильтрационных экранов ГАЭС. Выполненными исследованиями показано, что нарушение естественной структуры глины при строительстве экрана и устройство его пригрузки в виде защитного покрытия практически исключают набухание глины и увеличение ее водопроницаемости. Экран, возведенный из этих глин, характеризуется коэффициентом фильтрации Кф = А-10-9 см/с, где А = 1...10.
Противофильтрационный экран, возведение, набухающая глина,
водопроницаемость, структура глины, коэффициент фильтрации, плотность грунта, влажность грунта.
The assessment of the possibility of using swelling clays for construction of watertight screens of PSP (pump storage plant) is given. The carried out researches have shown that breakage of the clay natural structure during construction of the screen and the way of its cantledging as a protective coating practically eliminate clay swelling and increase of its water tightness. The screen built of these clays is characterized by filtration coefficient Kf=A-10-9 cm/s, where A= 1...10.
Watertight screen, construction, swelling clay, permeability, clay structure, filtration coefficient, soil density, soil humidity.
В состав основных сооружений ги- эксплуатации за счет набухания глины в
дроаккумулирующих электростанций процессе ее водонасыщения. (ГАЭС) обязательно входит верхний Влияние набухания глинистого
аккумулирующий водоем, в котором грунта экрана на водопроницаемость
накапливается вода для выработки электро- изучали на примере Днестровской
энергии. Чтобы минимизировать филь- ГАЭС. Исследовали набухание об-
трационные потери воды, закаченной в разцов неогеновой глины ненарушенной
водоем снизу, и исключить возможность структуры при естественной влажности
ухудшения гидрогеологической обстанов- и образцов нарушенной структуры, сфор-
ки в районе расположения ГАЭС, верхний мированных уплотнением валовых проб
аккумулирующий водоем экранируют. глины до состояния, характеризующегося
При проектировании и строитель- полутвердой консистенцией и коэффи-
стве противофильтрационных экранов из циентом водонасыщения 0,89...0,94, соот-
глинистых грунтов для водоемов ГАЭС ветствующим требованиям технических
большое внимание уделяется оценке условий на возведение экрана. возможности изменения их водопро- В естественной залежи тонкоди-
ницаемости во время строительства и сперсные неогеновые глины характери-