Полный текст

Введение Рассмотрим фильтрационные задачи строительной геотехники при движении воды в тонких проницаемых слоях грунтов и искусственных материалов. Геологические разрезы в условиях городской застройки имеют существенно техногенный характер [1-6]. Центры городов постоянно ремонтируют и перестраивают. В ходе такой строительной деятельности оказываются погребенными слои искусственных материалов, например из асфальта, бетона, щебня, песка. Многие строительные конструкции в грунте имеют слоистое строение с чередованием проницаемых и водоупорных слоев. Типичным примером такой конструкции является городская автомобильная дорога, где есть слои как гидроизоляционные, так и дренирующие. В качестве гидроизоляции обычно выступает асфальтобетон, а в качестве дренирующих слоев традиционно применяют щебень, песок и в последнее время геотекстильный материал в виде сеток и нетканых полотен. Вода в городских условиях в тонких проницаемых слоях грунтов и искусственных материалов может иметь разное происхождение. Вода может возникать из атмосферы, просачиваясь через трещины и поры в виде инфильтрации или конденсируясь из влажного воздуха. Часть воды в эти тонкие слои попадает из поверхностных источников, водоемов. Кроме того, другой причиной появления фильтрационной влаги в рассматриваемых слоях может быть подток подземных вод. При повышении уровня подземных вод такой процесс называют подтоплением. Вода, появившаяся в рассматриваемых слоях, в редком случае может иметь горизонтальный пьезометрический уровень, при котором ее фильтрационное движение отсутствует. Исключением из этого случая может быть наличие либо вертикального оттока подземной влаги в нижележащий пласт, либо вертикальное испарение воды в атмосферу под влиянием температурных перепадов. Такие случаи в нашей работе не рассматриваются. Обычно пьезометрические уровни воды в проницаемых слоях находятся на разных высотных отметках. Это создает разность напоров, что приводит в фильтрационное движение воду в слоях. Вода начинает двигаться от большего напора к меньшему. Будем считать, что такие процессы происходят при соблюдении известного закона фильтрации Дарси, принимая ламинарный режим движения воды. Исходная постановка дифференциальных уравнений, описывающих нестационарную фильтрацию воды в изучаемых слоях, будет соответствовать гидравлической теории фильтрации без линеаризации. Основная часть Фильтрационное движение воды в тонких проницаемых слоях всегда происходит под влиянием разности напоров. Аналитические методы решения таких задач фильтрации наиболее удобно сводить к условно одномерным постановкам. Формулы, которые получают в результате решения этих задач, широко используют в проектно-изыскательском деле. При этом условно одномерные постановки в рамках гидравлической теории фильтрации могут быть рассмотрены для двух основных расчетных случаев. Первый случай связан с плоскопараллельными задачами, которые применяются для линейных источников подтопления или дренажных систем. Например, разрабатывают траншею, к которой может возникнуть с двух сторон приток подземных вод по техногенным тонким проницаемым слоям. Или, наоборот, произошел внезапный аварийный прорыв воды из инженерной сети в траншее и образовался линейно вытянутый источник подтопления окружающей территории. Требуется рассчитать процесс нестационарной фильтрации воды по техногенным тонким проницаемым слоям с целью дальнейшей разработки защитных мероприятий от подтопления. Второй случай относится к радиальной в плане фильтрации воды, которая течет по направлению к центру или, наоборот, от центра источника подтопления, вернее, от его границы с пластом к наружным границам области фильтрации. Такой случай типичен в изыскательской практике при определении фильтрационных параметров водоносных пластов с помощью откачек воды из скважин. Однако в рассматриваемой теме фильтрации в техногенных тонких проницаемых слоях более актуально рассмотреть приток воды к небольшим приямкам или же, наоборот, решить задачу о подтоплении окружающей территории при постепенной утечке воды из небольшого котлована или выемки с растеканием купола подтопления. Такие две характерные для рассматриваемой темы концептуальные задачи рассмотрим далее по порядку. Пример для первого случая рассмотрим классический, который ранее был исследован в работе [7]. Однако мы решим его другим, новым, способом. Сравнение позволит нам критически оценить качество новой методики. Исходная постановка фильтрационной задачи о подтоплении первоначально не обводненного тонкого проницаемого слоя грунта или искусственного материала при подпоре воды из линейного источника подтопления показана на рис. 1. Для этого случая постановка дифференциального уравнения, описывающего нестационарную фильтрацию воды в проницаемом тонком слое, будет соответствовать гидравлической теории фильтрации без линеаризации. Исходное уравнение такого процесса фильтрации можно записать в виде (1) где левая часть уравнения соответствует изменению расхода воды по закону Дарси, причем мощность (глубина, напор) слоя воды H стоит под знаком дифференциала, указывая на отсутствие линеаризации. Правая часть уравнения отражает нестационарное изменение насыщения пористой среды водой. Кроме того, обозначены: k - коэффициент фильтрации грунта или искусственного материала, м/сут; m - коэффициент недостатка насыщения; t - время процесса подтопления, сут; x - горизонтальная координата, м. Рис. 1. Зона подтопления Lп из линейно вытянутого источника Fig. 1. Underflooding zone Lп from a linearly elongated source Краевое условие на левой границе области фильтрации (см. рис. 1): (2) Краевое условие на правой границе области фильтрации (фронт подтопления): (3) Величина зоны подтопления Lп изменяется с течением времени, увеличиваясь в размере (см. рис. 1). Начальное распределение напоров отсутствует, и область фильтрации увеличивается в размере. Вначале, решая краевую задачу (1) с условиями (2) и (3), будем следовать методологии работы [7]. Дойдя до определенного шага, отметим его и покажем, как дальше можно решить эту задачу по-новому. Преобразуем размерные переменные x и H к двум безразмерным величинам. Первая переменная отражает безразмерную пространственную координату: (4) Вторая переменная является безразмерным относительным напором: (5) Тем самым вместо размерных переменных были введены автомодельные безразмерные переменные. Основы данного метода решения с помощью автомодельных безразмерных переменных более подробно изложены в работах [8-12]. Заменим в уравнении (1) размерные переменные на безразмерные. После выкладок приходим к дифференциальному уравнению в следующем виде: (6) Обозначим в (6) первую и вторую производные U по S соответственно одним и двумя штрихами и перепишем это уравнение в более компактном виде: (7) Приравнивая в (7) U = 0 и учитывая, что U'' ≠ ∞, получим откуда при U' ≠ 0 вытекает, что (8) Теперь переосмыслим задачу от рис. 1 к рис. 2 для лучшего пояснения метода. Рис. 2. Линейная зона подтопления в автомодельных переменных Fig. 2. Linear underflood zone in self-similar variables Длина зоны подтопления из формулы (4) может быть записана как (9) В работе [7] аналитическим путем найдено, что B = 1,143. Поэтому окончательно в [1] была получена длина зоны подтопления в виде (9а) Начиная с этого момента покажем, как найти другим способом, по-новому, коэффициент B, отражающий линейную зону подтопления в автомодельных переменных (см. рис. 2). Для этого используем комбинированный аналитико-численный способ с применением метода конечных разностей (МКР). Запишем уравнение (6) в конечно-разностном виде так: (10) Для МКР-уравнения (10) приводим поясняющий рис. 3. Рис. 3. Конечно-разностная сетка для автомодельных переменных Fig. 3. Finite difference mesh for self-similar variables Из уравнения (10) выразим значение безразмерного напора в i-той центральной ячейке МКР-сетки в виде (11) Дальнейшее нахождение искомой величины B можно выполнить с помощью компьютерной программы на любом языке программирования интерпретирующего или компилирующего типа. Главным требованием к языку программирования является поддержка рекурсии. Запишем алгоритм действий по отысканию B в виде псевдокода, который можно переписать на любом из упомянутых языков программирования. Пример псевдокода по нахождению величины B может быть изложен так: - пусть первое приближение B = 1,000; - назначаем шаг МКР-сетки DS = 0,005; - принимаем количество шагов МКР-сетки равным 1000; - дано напор на левой границе области фильтрации U = 1; - дано напор на правой границе области фильтрации U = 0; - записываем напор в узел сетки по формуле (11); - рассчитываем всю матрицу напоров МКР-модели; - проверяем, что счет устойчив и процесс расчета сходится; - если «да», то увеличиваем на один шаг DS величину B и рассчитываем матрицу; - если «нет», то предыдущее B найдено и записано в файл, расчет закончен. Программировать вовсе не обязательно по данному псевдокоду. Можно воспользоваться любой электронной таблицей, например Excel, LibreCalc, Gnumeric и т.д. [13]. Подбор величины B следует производить в табличной МКР-модели, вводя в соседние ячейки таблицы формулу (11). Перед такой подготовкой модели надо в настройках таблицы разрешить выполнение циклических ссылок на соседние ячейки таблицы. Таким путем, вручную, постепенно, можно подобрать величину B. В результате получим ответ B = 1,143. Это полностью совпадает со значением B, найденном аналитически в работе [7]. Следовательно, наш новый метод позволяет решать нелинейные задачи фильтрации смешанным аналитико-численным способом. Для исследователей, больше имеющих склонность к программированию, нежели к аналитическим выкладкам, такой метод может быть полезен для получения поисковых результатов. Для исследователей, пока еще неопытных в программировании, можно воспользоваться МКР-моделированием вручную с помощью любого известного редактора электронных таблиц. Рис. 4. Радиус подтопления Rп из компактного источника Fig. 4. Radius of underflooding Rп from a compact source Рассмотрим теперь второй случай, относящийся к радиальной в плане фильтрации воды, которая течет по направлению к центру или, наоборот, от центра источника подтопления, вернее, от его границы с пластом к наружным границам области фильтрации. В рассматриваемой теме фильтрации в техногенных тонких проницаемых слоях актуально рассмотреть приток воды к небольшим приямкам или же, наоборот, решить задачу о подтоплении окружающей территории при постепенной утечке воды из небольшого котлована или выемки с растеканием купола подтопления. Далее рассмотрим радиальную задачу о подтоплении окружающей территории от небольшого котлована, приямка или компактного углубления (рис. 4). Для нахождения решения в радиальных координатах будем применять новый способ, обоснованный в предыдущей задаче плоскопараллельной фильтрации воды. Предположим, требуется найти радиус подтопления Rп (см. рис. 4) при фильтрации воды из компактного источника, в качестве которого может быть котлован, приямок, компактное углубление, заполненные водой, например, после аварии на какой-либо инженерной системе водоснабжения, водоотведения или теплоснабжения. Радиус подтопления надо найти в виде (12) Здесь BR - неизвестный пока коэффициент, который надо подобрать; r0 - радиус источника подтопления; все остальные обозначения аналогичны предыдущей задаче. Исходное дифференциальное уравнение, описывающее радиальную нестационарную фильтрацию воды в проницаемом тонком слое, будет соответствовать гидравлической теории фильтрации без линеаризации в следующем виде: (13) где левая часть уравнения соответствует изменению расхода воды по закону Дарси, причем мощность (глубина, напор) слоя воды H стоит под знаком дифференциала, указывая на отсутствие линеаризации. Правая часть уравнения отражает нестационарное изменение насыщения пористой среды водой. Вода радиально растекается в плане. Кроме того, обозначены: k - коэффициент фильтрации грунта или искусственного материала, м/сут; m - коэффициент недостатка насыщения; t - время процесса подтопления, сут; r - радиальная координата, м. Краевое условие на внутренней границе области фильтрации (см. рис. 4): (14) Краевое условие на внешней границе области фильтрации (фронт подтопления): (15) Величина радиуса подтопления Rп изменяется с течением времени, увеличиваясь в размере (см. рис. 4), т.е. сама область фильтрации расширяется, причем за ее пределами, снаружи радиуса подтопления Rп, воды нет. Преобразуем размерные переменные r и H к двум безразмерным величинам n и U. Первая переменная отражает безразмерную радиальную пространственную координату: (16) Вторая переменная U является безразмерным относительным напором по формуле (5). Тем самым вместо размерных переменных были введены автомодельные безразмерные переменные для радиальной задачи. Подставляя переменные n и U в уравнение (13) и переписывая его, соответственно, будем иметь другое уравнение в виде Последнее уравнение запишем в МКР-представлении в виде (17) Из уравнения (17) выразим безразмерный напор в i-й центральной ячейке МКР-сетки для радиальной задачи в виде (18) Действуя по вышеизложенному алгоритму псевдокода аналогично плоскопараллельной задаче, применяя смешанный аналитико-численный метод, получим в результате формулу для нахождения без линеаризации радиуса подтопления при фильтрации воды из компактного источника, в качестве которого могут быть небольшой котлован, приямок, компактное углубление и т.д., заполненные водой, например, после аварии на какой-либо инженерной системе водоснабжения, водоотведения или теплоснабжения в виде (см. рис. 4): (19) Таким образом, в представленной работе была сделана попытка показать на примерах, как решить часть задач нестационарных фильтрационных течений в тонких водоносных пластах, свойственных строительной геотехнике, с получением аналитических решений конечного вида без введения искусственных упрощений в виде линеаризации исходных дифференциальных уравнений фильтрации. Предложена методология сочетания аналитических и численных методов в процессе отыскания решений фильтрационных задач строительной геотехники с получением автомодельных зависимостей. Такой подход позволит избежать ошибок линеаризации и проводить фильтрационные расчеты в тонких проницаемых слоях грунтов и искусственных материалов более точно. Обсуждение Одной из первых работ, близкой к примерам, представленным в настоящей статье, является публикация [14]. Постановка задачи соответствовала плоскопараллельной фильтрации воды по горизонтальному водоупорному слою (см. рис. 1). Исходное дифференциальное уравнение фильтрации было нелинейным по типу (1). Начальные и граничные условия задачи почти соответствовали (2) и (3), но дополнительно на фронте растекания воды было записано еще одно условие в виде (20) где все обозначения такие же, как для (1), (2) и (3), при этом имеется в виду, что Lп - величина переменная. По-видимому, условие (20) было избыточным и слишком усложнило задачу, поэтому решение в [14] выписано лишь в общем виде, но не доведено до расчетной формулы для величины Lп (см. рис. 1). Окончательное решение в замкнутом виде для плоскопараллельной фильтрации воды было получено в статье [7]. Это решение (9). В практике гидрогеологии оно считается классическим. В дальнейшем оно было опубликовано в монографии [15]. Наше решение по радиальной фильтрации получено в виде формулы (19) для радиуса подтопления при растекании воды в грунте или искусственном материале из компактного источника, в качестве которого могут быть небольшой котлован, приямок, компактное углубление и т.д., заполненные водой, например, после аварии на какой-либо инженерной системе водоснабжения, водоотведения, теплоснабжения. Это решение прошло независимую проверку в лабораторных условиях и показало достаточную надежность с погрешностями результатов не более 5 %. Обзор иностранных публикаций по рассматриваемой тематике показал некоторый застой в этом направлении. Тем не менее следует отметить интересную работу [16], где сделан достаточно категоричный вывод о неопределенности проверок теоретических методов и методик по определению эффективных параметров гидрогеологических моделей в условиях природных пластов. Переосмысляя данную работу применительно к нашему исследованию, следует отметить гораздо большую эффективность при проверке полученных расчетных формул теории фильтрации именно лабораторных условий, нежели природных. Это можно объяснить гарантированной определенностью тщательной подборки фильтрационных параметров лабораторных моделей. Подобный подход ранее встречался в работе [17]. Выводы Таким образом, рассмотрены некоторые фильтрационные задачи строительной геотехники при движении воды в тонких проницаемых слоях грунтов и искусственных материалов. В таких условиях традиционные подходы, заимствованные из классической гидрогеологии, дают погрешности. Это связано с тем, что в обычных условиях подземных водоносных пластов достаточно большой глубины в гидрогеологической практике используют гидравлическую теорию фильтрации, при которой усредняют скорости фильтрационных течений по вертикальной координате. Исходные нелинейные дифференциальные уравнения фильтрационных потоков подвергают принудительной линеаризации тем или иным известным способом. При этом широко используют способ линеаризации введением средней мощности (глубины) фильтрационного течения. Этот способ изначально разработан для водоносных пластов большой глубины. В представленной работе сделана попытка отказаться от подобной линеаризации и решить часть задач нестационарных фильтрационных течений в тонких водоносных пластах, свойственных строительной геотехнике, с получением аналитических решений конечного вида без введения искусственных упрощений. Предложена методология сочетания аналитических и численных методов в процессе отыскания решений фильтрационных задач строительной геотехники с получением автомодельных зависимостей. Такой подход позволит избежать ошибок линеаризации и проводить фильтрационные расчеты в тонких проницаемых слоях грунтов и искусственных материалов более точно. Обзор отечественных и зарубежных научных публикаций по затронутой теме исследования показал некоторый застой в этом направлении. Сопоставление выявило, что в тонких проницаемых слоях грунтов и искусственных материалов предпочтение можно отдать фильтрационным расчетам по представленной методике, так как она более точно описывает нестационарное движение воды в строительной геотехнике.

Об авторах

В. И. Сологаев

Сибирский государственный автомобильно-дорожный университет

Список литературы

Статистика

Просмотры

Аннотация - 114

PDF (Russian) - 53

Ссылки

• Ссылки не определены.