Full Text

Введение Подтопление подземными водами городских территорий усложняет проведение земляных работ в ходе строительства. Технологии разработки траншей и котлованов в условиях городской застройки требуют искусственного ограничения контура земляных выемок для защиты от подтопления, обрушения и оплывания. Применение для этого различных креплений в виде шпунтовых свай, стен в грунте и прочих аналогичных конструкций создает участки неоднородности в отношении проницаемости окружающего выработки грунта. Актуальным направлением исследований является совершенствование фильтрационных расчетов при проведении земляных работ в условиях подтопления в городском строительстве с учетом появления техногенных участков неоднородности по проницаемости грунта, окружающего траншеи и котлованы. В представленной работе рассмотрены наиболее распространенные расчетные фильтрационные схемы при выполнении земляных работ в городском строительстве. Произведена исходная постановка соответствующих фильтрационных задач. Поставлены граничные и начальные условия для областей фильтрации подземных вод. Решение задач найдено с применением методов конформных отображений, эквивалентных фильтрационных сопротивлений и интегральных преобразований Лапласа в нестационарной постановке движения подземных вод. Решения снабжены конкретными практическими примерами. Полученная усовершенствованная методика фильтрационных расчетов рекомендуется к применению при выполнении земляных работ в городском строительстве. Основная часть Рассмотрим наиболее распространенные расчетные фильтрационные схемы при проведении земляных работ в городском строительстве, в которых имеются характерные участки искусственных неоднородностей в виде различных креплений - шпунтовых свай, стен в грунте и прочих аналогичных конструкций. Вначале перечислим эти характерные схемы, затем произведем постановку задач и найдем их решения аналитическими методами, комбинируя конформные отображения, эквивалентные фильтрационные сопротивления и интегральные преобразования Лапласа в нестационарной постановке. Это позволит найти достаточно удобные и не очень сложные формулы для расчета реальных процессов фильтрации подземных вод, которые образуют притоки в строительные выработки. Полученные зависимости для расходов и напоров фильтрационных потоков при выполнении земляных работ будут носить оценочный характер, то есть достаточно приближенный. Тем не менее, именно с таких оценок и следует начинать предварительные фильтрационные расчеты при защите от подтопления в городском строительстве. Дальнейшее уточнение режима движения подземных вод в зонах строительных выработок можно произвести с помощью компьютерного моделирования фильтрации воды, однако такие последующие, более точные, шаги в данной работе не рассматриваются. Три распространенные расчетные фильтрационные схемы при проведении земляных работ в городском строительстве, в которых имеются характерные участки искусственных неоднородностей в виде различных креплений или конструкций изображены на рис. 1-3. Приведем некоторые пояснения к иллюстрациям. На рис. 1-3 ось симметрии относится к траншее или вытянутому котловану. Непроницаемая стенка 1 обозначает конструкцию в грунте, создающую преграду для движения подземных вод и крепление грунта от обрушения. В случае проницаемой стенки 1, показанной на рис. 3, считаем, что такая конструкция удерживает от обрушения грунт, но немного пропускает подземную воду. Водоносный пласт 3 может быть напорный или безнапорный. Грунт 4, как вариант, состоит из водоупорных пород, либо, напротив, он несколько пропускает подземную воду, то есть является проницаемым грунтом в виде относительного водоупора. УПВ - уровень подземных вод. Lвл - зона влияния земляных работ на естественный уровень подземных вод, который обозначен как «Ест. УПВ». Другие обозначения, изображенные на рисунках, оговорены ниже по ходу изложения. Рис. 1. Расчетная схема 1: 1 - непроницаемая стенка; 2 - пониженный уровень подземных вод; 3 - водоносный пласт; 4 - водоупорный грунт Fig. 1. Calculation scheme 1: 1 - impermeable wall; 2 - reduced groundwater level; 3 - aquifer; 4 - waterproof soil Рис. 2. Расчетная схема 2: 1 - непроницаемая стенка; 2 - пониженный уровень подземных вод; 3 - водоносный пласт; 4 - проницаемый грунт (относительный водоупор); 5 - движение подземных вод в сторону выработки Fig. 2. Calculation scheme 2: 1 - impermeable wall; 2 - reduced groundwater level; 3 - aquifer; 4 - permeable soil (relative aquiclude); 5 - groundwater movement towards the excavation Рис. 3. Расчетная схема 3: 1 - проницаемая стенка; 2 - пониженный уровень подземных вод; 3 - водоносный пласт; 4 - водоупорный грунт Fig. 3. Calculation scheme 3: 1 - permeable wall; 2 - lowered groundwater level; 3 - aquifer; 4 - waterproof soil По решению подобных задач вначале сделаем небольшое замечание по обзору имеющихся публикаций за последние годы. Недавно появилась работа исследователей из Китая и Австралии [1] по рассматриваемой тематике в виде аналитических решений, которые дают возможность инженеру выполнить предварительные расчеты напоров и притоков подземных вод в грунтах, окружающих строительные выемки в виде круглых в плане котлованов с участками фильтрационных неоднородностей - ограждающих непроницаемых конструкций. Сравнение результатов нашей работы и упомянутой статьи [1] обсуждено в следующем разделе. Современная всеобщая тенденция проектирования подобных задач опирается преимущественно на компьютерные технологии в виде численного моделирования фильтрационных потоков в грунте. При этом наибольшее распространение получил метод конечных элементов (МКЭ), а также, как альтернатива, - метод конечных разностей (МКР). Данные методы положены в основу программных комплексов гидрогеологической направленности. МКЭ, например, применен в продукте Seep/W [2-5], а МКР - в известном MODFLOW, с примерами которого можно ознакомиться онлайн в работах [6-9]. Метод же получения аналитических решений, позволяющий найти на основе решения краевых задач приближенные инженерные форму стал несколько менее популярен среди исследователей. Поэтому представляется актуальным направлением возобновление научного интереса к аналитике, дающей приближенные зависимости, которые позволяют сделать более общие, все охватывающие выводы на их основе в отличие от более точной, но узко направленной технологии численного моделирования геофильтрации. Приступим к решению поставленных задач. В случае расчетной схемы 1 непроницаемая стенка 1 немного не доходит своим низом до нижележащего водоупора (см. рис. 1). Поэтому движение подземных вод по водоносному пласту в сторону строительной выработки является двумерным в вертикальной плоскости. Аналитически такие задачи обычно решают методом конформных отображений, позволяющих учесть двумерную фильтрацию и получить зависимость напоров и расходов подземных вод с учётом горизонтальной и вертикальной координат. Решим эту задачу в два этапа с использованием конформных отображений, фильтрационных сопротивлений и интегральных преобразований Лапласа для учета нестационарной фильтрации подземных вод. На первом этапе расчетную схему 1 (см. рис. 1) условно заменим на другую схему, не содержащую вертикальных линий фильтрационного тока. На этой, другой схеме, вместо непроницаемой стенки, не доходящей своим нижним концом до водоупорных пород, покажем участок проницаемой неоднородности, перекрывающей полностью движение подземных вод в строительную выработку, но немного пропускающей воду. Материал этой искусственной неоднородности будем считать имеющим коэффициент фильтрации kn, который по своему значению несколько меньше коэффициента фильтрации k водоносного пласта. Толщину этой неоднородности в горизонтальном направлении обозначим d. В результате получим новый рис. 4 взамен рис. 1. Линии тока, подходящие к строительной выработке, стали параллельными горизонтали, то есть движение подземных вод условно от двумерного перешло в одномерное. Решать такие задачи существенно проще. Коэффициент фильтрации искусственной неоднородности kn надо определить такой величины, чтобы фильтрационное сопротивление реальной неоднородности, показанное на рис. 1, было эквивалентно фильтрационному сопротивлению искусственной неоднородности, показанной на рис. 4. На втором этапе решения задачи составим уравнение баланса фильтрационного потока в строительную выработку, сопоставляя рис. 1 и 4. Вначале запишем значение расхода подземных вод для случая, показанного на рис. 1 в виде где z1, z2, z3 - фильтрационные сопротивления трех участков размером соответственно B, d, Lвл; DH - перепад напоров между Hе и Hк, вызывающий движение подземных вод в строительную выработку. Рис. 4. Искусственная неоднородность с одномерной фильтрацией: 1 - проницаемая стенка толщиной d с коэффициентом фильтрации kn; 2 - пониженный уровень подземных вод; 3 - водоносный пласт; 4 - водоупорный грунт Fig. 4. Artificial inhomogeneity with one-dimensional filtration: 1 - permeable wall with thickness d with filtration coefficient kn; 2 - lowered groundwater level; 3 - aquifer; 4 - waterproof soil Значение каждого сопротивления z1, z2, z3 можно определить с применением конформных отображений. Для этого применим методику, приведенную в [10]. Эти сопротивления будут такие: Теперь запишем значение расхода подземных вод для случая, показанного на рис. 4 в виде где x1, x2 - фильтрационные сопротивления двух участков размером соответственно d, Lвл; DH - перепад напоров. Значение каждого сопротивления x1, x2 можно также определить с применением конформных отображений. Для этого опять используем рекомендации из [10]. Эти сопротивления выглядят так: Задавшись условием равенства фильтрационных расходов по вышеприведенным формулам для q1 и q4, выразим коэффициент фильтрации искусственной неоднородности kn для расчетной схемы 1 (см. рис. 1) в виде (1) Для сравнения приведем известную упрощенную формулу для kn, которая отражает фильтрационное движение подземных вод лишь параллельно разнородным по проницаемости слоям [11] (2) однако заметим, что эта формула не учитывает двумерное движение подземных вод в выработки. В этой формуле для частного случая схемы 1 ноль соответствует непроницаемой конструкции. В реальных условиях в строительных выработках при выполнении земляных работ движение подземных вод происходит в нестационарном режиме. Поэтому расчетные формулы для напоров и расходов получим на основе решения системы исходных дифференциальных уравнений фильтрации в виде Здесь t - время; m - коэффициент гравитационной или упругой водоотдачи; H1 - напор внутри искусственной неоднородности; H2 - напор в прилегающей внешней зоне водоносного пласта. Граничные и начальные условия см. на рис. 4. Hе - начальное распределение напоров при t = 0. Граничное условие для координаты x в виде бесконечности условно, так как в действительности влияние возмущений на уровень подземных вод (УПВ) имеет вполне ограниченное действие [12]. Наиболее интересна в рамках рассматриваемой работы 1-я зона фильтрации с напорами H1 вблизи контура строительной выработки, где возникают проблемы с повышенными градиентами напоров подземных вод и опасностью обрушения грунта. Приведенная система уравнений представляет разновидность гидравлической теории фильтрации [13]. Это означает, что вертикальными компонентами скоростей фильтрации можно пренебречь. По меткому замечанию в книге [14], эта предпосылка Дюпюи-Форхгеймера, не взирая на ее приближенность, дает неплохие результаты совпадения теории и практики в гидрогеологических исследованиях. Мощность (толщина) M относится к напорному водоносному пласту. В случае безнапорного пласта значение M надо заменить на величину средней мощности водоносного пласта hc = (Hе + Hк) / 2. Это соответствует 1-му способу линеаризации нелинейных дифференциальных уравнений безнапорной фильтрации, рекомендованному для прогнозных задач в [13]. Общее решение приведенной выше системы параболических дифференциальных уравнений известно [15]. Научный интерес представляет получение частного решения для мало изученного начального периода времени земляных работ в зонах грунта, примыкающих к строительным выработкам с фильтрационными неоднородностями. Применяя к системе вышеприведенных уравнений интегрирование по Лапласу, то есть прямое преобразование Лапласа, найдем в частном случае при малых величинах t (начальный период времени земляных работ) для зоны внутри искусственной неоднородности выражение в изображениях напоров (со знаком надчеркивания) в виде Здесь p - промежуточный параметр Лапласа, который нужен для эффективного решения с последующим обратным преобразованием выражения в изображениях напоров. Применим к последнему выражению в изображениях обратное преобразование Лапласа. В результате, после несложных выкладок, получим формулу для расчета в нестационарном режиме напоров подземных вод в зоне искусственной неоднородности толщиной d (см. рис. 4), которая выглядит так: (3) В этой формуле erfc (…) - математическая функция дополнительного интеграла вероятности, которая является встроенной в большинство распространенных электронных таблиц на персональных компьютерах. Умножим полученное выражение на величину knM/B (см. обозначения на рис. 1) и продифференцируем по x. Тогда получим формулу для расчета нестационарного водопритока подземных вод в строительную выработку с учетом искусственной неоднородности (см. рис. 4) в виде (4) Рассмотрим теперь расчетную схему 2, показанную на рис. 2. В этом случае непроницаемая стенка 1 полностью перекрывает водоносный пласт и заходит своим низом в относительный водоупорный слой. Движение подземных вод по водоносному пласту в сторону строительной выработки является двумерным в вертикальной плоскости с обтеканием защитной конструкции снизу через относительный водоупор. Заменим условно схему 2 (см. рис. 2) на схему 4 (см. рис. 4) с целью возможности проведения расчета нестационарной фильтрации подземных вод в траншею или вытянутый котлован. Вначале необходимо найти коэффициент фильтрации kn для искусственной неоднородности, показанной на схеме 4, эквивалентный по расходу схеме 2. Используя опять метод конформных отображений [16], запишем приток в строительную выработку при обтекании фильтрационным потоком непроницаемой стенки через относительный водоупор (см. рис. 2) в виде Эквивалентный расход через искусственную неоднородность с одномерной фильтрацией (см. рис. 4) будет выглядеть так: Приравнивая расходы qв и qn, после несложных выкладок выразим расчетный коэффициент фильтрации условной искусственной неоднородности для схемы 4 (см. рис. 4), заменяющей расчетную схему 2, в виде (5) Дальнейший фильтрационный расчет напоров и расходов при выполнении земляных работ для схемы 2 (см. рис. 2) можно производить по ранее полученным формулам (2) и (3), подставляя в них коэффициент фильтрации условной искусственной неоднородности по формуле (5). Наконец, рассмотрим расчетную схему 3, показанную на рис. 3. В этом случае стенка 1 на краю строительной выработки вовсе не является водонепроницаемой, а лишь удерживает грунт от обрушения. Движение подземных вод по водоносному пласту в сторону строительной выработки является одномерным в вертикальной плоскости с просачиванием через защитную конструкцию стенки. Фильтрационный расчет напоров и расходов при выполнении земляных работ для схемы 3 (см. рис. 3) можно производить по ранее полученным формулам (3) и (4), подставляя в них коэффициент фильтрации проницаемой стенки. Таким образом, мы рассмотрели три распространенные расчётные схемы фильтрационных водопритоков к строительным выработкам (траншеям и вытянутым котлованам), которые представили на рис. 1-3. Для всех трех случаев получена расчётная методика приближённого определения напоров и расходов подземных вод, просачивающихся к строительным выработкам при наличии в прилегающем грунте участков искусственных неоднородностей в виде ограждающих стенок. Нетрудно заметить, что все три расчетных схемы, показанные на рис. 1-3, сводятся к условной эквивалентной по расходу схеме одномерной фильтрации подземных вод, изображенной на рис. 4. Данные расчёты позволяют найти характеристики фильтрационных потоков в нестационарном режиме движения подземных вод. Это соответствует реальным условиям проведения земляных работ на городских стройках. Приведем пример сравнительного расчета по разработанной методике фильтрационного прогноза подтопления и дренирования при выполнении земляных работ в городском строительстве. Допустим, что в ходе проектирования на основе данных изысканий выяснилось, что строительная выработка соответствует расчетной схеме 1, показанной на рис. 1. Пусть будут такие исходные данные. Железобетонная стена в грунте толщиной d = 0,6 м перекрывает напорный водоносный пласт мощностью M = 11 м, грунт которого имеет коэффициент фильтрации k = 0,1 м/сут и коэффициент водоотдачи m = 0,05. Стена в грунте не доходит до водоупора на величину m = 0,5 м. Напор водоносного пласта Hе = 16 м. Напор в котловане Hк = 6 м. Разность напоров в пласте и выработке DH = 16 - 6 = 10 м. Половина ширины выработки B = 12 м. Этот случай, например, характерен для строительства котлована станции метрополитена мелкого заложения открытым способом. Вначале найдем по формуле (1), учитывающей двумерное движение фильтрационных потоков, коэффициент фильтрации искусственной неоднородности kn, м/сут: Затем, для методического сравнения, определим по упрощенной формуле (2) этот же коэффициент фильтрации искусственной неоднородности: и видим, что формула (2) в 4,2 раза завышает проницаемость искусственной неоднородности на контуре строительной выработки вследствие отсутствия учета в ней двумерного движения подземных вод вблизи котлована или траншеи. Поэтому надежнее использовать формулу (1). Таким образом, принимаем для дальнейшего расчета коэффициент фильтрации искусственной неоднородности со значением Определим теперь по формуле (3) напор подземных вод с внешней стороны непроницаемой стенки при x = d (см. рис. 1) спустя t = 10 сут от начала выполнения земляных работ с водоотливом: Опять для сравнения рассчитаем по формуле (3) напор подземных вод при x = d в случае отсутствия непроницаемой стенки (см. рис. 1) также как и в предыдущем расчете спустя t = 10 сут от начала выполнения земляных работ с водоотливом. При этом водоносный грунт всюду предполагаем однородным с коэффициентом фильтрации k = 0,1 м/сут. Получим следующее значение напора Это сравнение показывает, что недоучет наличия искусственной неоднородности при выполнении земляных работ с водоотливом приводит к слишком оптимистичному прогнозу снижения напора подземных вод вблизи строительных выработок. На самом деле напор снижается гораздо медленнее, что утяжеляет условия проведения строительства. Поэтому для повышения надежности проведения таких строительных работ рекомендуется к использованию при проектировании методика фильтрационных расчётов, предложенная в представленной статье. Обсуждение Наиболее свежей и близкой научной публикацией по затронутой теме является статья авторов из Китая и Австралии [1]. Имеет смысл сделать сравнительное обсуждение данной и нашей работ, раз они близки и перекликаются. Все обоснования, выкладки и расчёты приведены нами в предыдущем разделе. Авторы [1] пишут, что при откачке внутри котлована, огражденном водонепроницаемой завесой, направление движения подземных вод за пределами выемки изменяется с горизонтального на вертикальное за счет наличия завесы-барьера. При этом они утверждают, что не существует аналитического метода расчета распределения напора подземных вод, вызванного осушением внутри котлована. Затем они предлагают простое уравнение, основанное на аналитическом решении, для расчета напоров грунтовых вод внутри и снаружи котлована с водонепроницаемой завесой в напорном водоносном горизонте. Решение в [1] основано на комбинациях применения известной формулы Дюпюи для стационарного режима радиального течения подземных вод. Наше же решение, сравнительно с упомянутой работой [1], получено для более реалистичного неустановившегося режима движения подземных вод в земляную выработку, что обычно более свойственно для технологий городского строительства с конкретными сроками выполнения работ. Тем не менее, и та, и другая методики дополняют друг друга и могут быть рекомендованы к практическому применению при защите от подтопления в городском строительстве. Выводы Таким образом, рассмотрена нестационарная фильтрация подземных вод и поведение их уровня вблизи строительной выработки с характерными участками искусственных неоднородностей в виде различных креплений или конструкций. Исследование показало, что недоучет наличия искусственной неоднородности при выполнении земляных работ с водоотливом приводит к слишком оптимистичному прогнозу снижения напора подземных вод вблизи строительных выработок. На самом деле напор снижается гораздо медленнее, что утяжеляет условия проведения строительных работ. Полученная усовершенствованная методика фильтрационных расчетов рекомендуется к применению при выполнении земляных работ в городском строительстве.

About the authors

V. I Sologaev

Siberian State Automobile and Highway University

References

Statistics

Views

Abstract - 86

PDF (Russian) - 107

Refbacks

• There are currently no refbacks.