Работа в PLAXIS 3D. Трехмерное численное исследование эффективности подземного дренажа при крупномасштабных оползнях — ГеоИнфо — метапортал для инженеров
Реклама
  • Реклама, 0+. ООО «ИнжПроектСтрой» ИНН 5902163884
  • erid: 2vfnxvifrnd
Баннер MalininSoft правая колонка Баннер MalininSoft правая колонка
Реклама
  • Реклама, 0+. АО «Мостдоргеотрест» ИНН 7716750744
  • erid: 2vfnxwa1cem
Баннер МОСТДОРГЕОТРЕСТ правая колонка Баннер МОСТДОРГЕОТРЕСТ правая колонка

Работа в PLAXIS 3D. Трехмерное численное исследование эффективности подземного дренажа при крупномасштабных оползнях

Работа в PLAXIS 3D. Трехмерное численное исследование эффективности подземного дренажа при крупномасштабных оползнях
Линь Дэж-Гуэй (Lin Der-Guey)
Линь Дэж-Гуэй (Lin Der-Guey)
Факультет охраны грунтов и вод Национального университета Чун Син, г. Тайчжун, Тайвань

Предлагаем вниманию читателей адаптированный перевод статьи «Трехмерное численное исследование эффективности подземного дренажа при крупномасштабных оползнях». Ее авторы – тайваньские специалисты Дэж-Гуэй Линь, Ко-Чин Чан, Чэн-Юй Ку и Цзюй-Чин Чоу. Они являются сотрудниками Национального университета Чун Син (г. Тайчжун) и Тайваньского национального университет океана (г. Килунг). Указанная работа была опубликована в 2020 году в рецензируемом журнале Applied Sciences («Прикладные науки»), выпускаемом интернет-издательством MDPI. Статьи на сайте MDPI находятся в открытом доступе и распространяются по лицензии Creative Commons Attribution (CC BY), которая разрешает их неограниченное использование, воспроизведение, распространение, перевод, редактирование, использование в своей работе (даже в коммерческих целях) при условии правильного указания ссылки на первоисточник.

В указанной статье представлены результаты исследований, основанных на мониторинге оползневого склона и направленных на проектирование противооползневого мероприятия в виде глубинного дренажа с использованием конечноэлементного программного комплекса PLAXIS 3D для анализа фильтрации грунтовых вод и устойчивости склона. Система подземного дренирования состояла из двух дренажных колодцев диаметром по 4 м с многоуровневыми горизонтальными дренами. Она была установлена для понижения уровня грунтовых вод и стабилизации неустойчивого оползневого склона. Полученные результаты показали, что такая система вполне функциональна и способна ускорять отвод грунтовых вод, накапливающихся при инфильтрации воды во время проливных дождей в сезон тайфунов. Это может защитить склон от дальнейшего разрушения и поддержать его устойчивость на приемлемом уровне.

Адаптированный перевод выполнен при поддержке партнера журнала «ГеоИнфо» – компании «НИП-Информатика».

Чан Ко-Чин (Chang Kuo-Ching)

Факультет охраны грунтов и вод Национального университета Чун Син, г. Тайчжун, Тайвань

Ку Чэн-Юй (Ku Cheng-Yu)

Центр передового опыта в области морской инженерии Тайваньского национального университет океана, г. Килунг, Тайвань

Чоу Цзюй-Чин (Chou Jui-Ching)

Факультет гражданского строительства Национального университета Чун Син, г. Тайчжун, Тайвань

1. ВВЕДЕНИЕ

Ю-Е-Линь (You-Ye-Lin) – крупный оползень в горном районе Тайваня с долгой историей периодических крупных смещений грунта вниз по склону во время сезонов дождей начиная с тайфуна «Моракот» 2009 года. Чтобы предотвратить дальнейшее развитие этого оползня и стабилизировать его, государственные службы начали проводить (и проводят до сих пор) систематические полевые исследования, прорабатывают проектные решения и реализуют стабилизирующие мероприятия для предотвращения негативных последствий. Основные работы заключаются в понижении уровня грунтовых вод (УГВ) с помощью системы глубинного подземного дренажа, состоящей из двух дренажных колодцев с многоуровневыми горизонтальными дренами.

При подземном дренаже эффективность работы горизонтальных дрен для понижения УГВ и повышения устойчивости склона тесно связана с их длиной, расстоянием между ними, их количеством и местами установки [1–5]. Функции и влияние горизонтальных дрен на понижение УГВ и на стабилизацию склона также были изучены несколькими исследователями с использованием численных методов [6–12]. Цай с соавторами [6] пришли к выводу, что удлинение группы горизонтальных дрен более эффективно, чем сокращение расстояний между ними и увеличение их количества. Рахарджо с коллегами [7–9] установили, что основным преимуществом использования горизонтальных дрен является понижение уровня грунтовых вод, причем наибольший вклад вносят дрены, расположенные в нижней части склона. Эберхардт с соавторами с уверенностью предположили [10], что система глубинного дренажа (дренажная галерея с субвертикальными дренами) является ключевой мерой для успешной стабилизации крупных оползней. Кроме того, с помощью трехмерных численных моделей для крупных оползней также было изучено потенциальное влияние на стабилизацию склона со стороны субвертикальных дрен, пробуренных из дренажных галерей [11, 12]. Линь с коллегами [13, 14] продемонстрировали, что подземная дренажная система в крупном оползне способна ускорить отвод грунтовых вод, накопившихся в результате инфильтрации дождевой воды при ливнях большой интенсивности и продолжительности во время сезона тайфунов, и защитить склон от быстрого разрушения.

Однако системы горизонтальных дрен в параллельных конфигурациях в предыдущих исследованиях были неглубоко расположенными и простыми по сравнению с системой горизонтальных дрен в оползне Ю-Е-Линь, где они были соединены с дренажными колодцами в многоуровневой веерообразной конфигурации. В частности, несмотря на то что на Тайване для повышения устойчивости крупных оползней было построено большое количество дренажных колодцев, ни один из достигнутых эффектов по снижению уровней грунтовых вод и стабилизации склонов до сих пор не подвергся тщательной проверке. В этом же исследовании оценивается эффективность системы подземного дренажа, установленной в оползне Ю-Е-Линь, на основе данных мониторинга и численных методов. Оно нацелено на выработку практического оперативного трехмерного метода для количественной оценки эффективности работы горизонтальных дрен при проливных дождях во время сезона тайфунов. В использованном методе оценка изменений уровня грунтовых вод с помощью сложного анализа инфильтрации дождевой воды в неводонасыщенный грунт заменяется данными мониторинга УГВ.

2. ОПОЛЗЕНЬ Ю-Е-ЛИНЬ

2.1. Расположение оползня

Как показано на рисунке 1, оползень Ю-Е-Линь расположен в округе Чай-И и занимает площадь около 40 га. Для движения транспорта в этом округе в основном используется автодорога № 166. Основными видами деятельности здесь являются производство чая, обработка побегов бамбука и туризм. Проливные дожди во время тайфуна «Моракот» с 6 по 9 августа 2009 года вызвали сход массивного оползня, который поставил под угрозу безопасность жилых домов, сельскохозяйственных культур и транспортной системы. Согласно данным долгосрочного мониторинга и распределению жилых районов был выбран участок в пределах оползня Ю-Е-Линь для дальнейшего изучения и расширенного анализа (он обведен на рисунке 1 красной штриховой линией).

Рис. 1. Спутниковые снимки острова Тайвань и территории расположения оползня Ю-Е-Линь (обведен желтой штриховой линией). Исследуемый участок обведен красной штриховой линией

2.2. Климат, топография и геология территории оползня

Средняя температура воздуха в районе оползня Ю-Е-Линь составляет около 23 °C. Среднегодовое количество атмосферных осадков (дождей) составляет около 3 050 мм, а среднегодовое количество дождливых дней составляет 166 на основании данных об осадках с 2003 по 2016 год. Наибольшая часть осадков выпадает в летний период (с мая по сентябрь). С октября по следующий март наблюдается сухой сезон, а сезон сильных дождей начинается в конце мая – начале июня.

Рельеф поверхности оползня Ю-Е-Линь волнистый, с небольшим ровным участком. Ее абсолютная высота уменьшается с севера на юг с 1 110 м до 820 м, в среднем составляя 1 000 м. По уклону S (%) или углу наклона β (град.) поверхность оползневой зоны варьирует от покатой (15% < S&lE; 30%; 8,53 град. < β&lE; 16,70 град.) до очень крутой и обрывистой (55% < S&lE; 100%; 28,81 град. < β&lE; 45 град.), но в основном является крутой (40% < S&lE; 55%; 21,80 град. < β&lE; 28,81 град.). Зона оползня может быть очерчена в водосборном бассейне реки Шэн-Мао-Шу, которая протекает через долину с южной стороны (см. рис. 1). Склон этой долины крутой, с хорошо развитыми эрозионными траншееобразными оврагами, борта которых склонны к обрушению. При этом из-за сильной продольной и поперечной эрозии во время сезона тайфунов русла реки Шэн-Мао-Шу и ее притоков становятся очень нестабильными. Кроме того, ручей Хай-Ту-Лунь, протекающий через оползневую зону с севера на юг, разделяет ее на восточный и западный топографические блоки (см. рис. 1).

Из-за водоразделов, образованных горными хребтами, тянущимися с восточной и с западной стороны от ручья Хай-Ту-Лунь, сток с водосборного бассейна собирается в этот ручей и стекает на юг в реку Шэн-Мао-Шу.

Удаление нижней части склона у ручья Хай-Ту-Лунь из-за сильной эрозии берегов может привести к обрушению этого склона. Такой механизм разрушения подтверждается широким распространением уступов и трещин растяжения на берегах данного ручья. Кроме того, нарушения, вызванные строительством главной дороги округа Чай-И (автодороги № 166), которая идет через весь оползневый склон Ю-Е-Линь, часто вызывают обвалы участков склона в сезон тайфунов.

Можно предположить, что сползающая грунтовая масса оползня Ю-Е-Линь представляет собой древний коллювий, подстилаемый коренной породой. Эти коллювиальные отложения скользят вниз по границе раздела с нижележащим песчаником, потому что они поглощают дождевую воду и увеличивают собственный вес. Обследование показало, что, похоже, нет никаких признаков того, что плоскость скольжения где-то проходит через коренную породу. Как следствие, движение грунта происходит только в пределах залегающих сверху коллювиальных отложений.

Простирание плоскостей напластования слоев грунта в оползне Ю-Е-Линь направлено на северо-восток. Поверхность склона параллельна падению слоев. Разлом Лу-Ку проходит через южный край оползня с простиранием на северо-восток и превращается во взброс с уклоном на восток и углом наклона плоскости разлома 50 град. Обнажения оползня Ю-Е-Линь, расположенные на западной стороне разлома Лу-Ку, могут быть отнесены к плиоценовой формации Чо-Лань, состоящей из мощного песчаника (Ss по американской единой системе классификации грунтов USCS), глинистого сланца (Sh) и чередующихся песчаника и глинистого сланца (Ss/Sh), как показано на рисунках 2, а–в.

Рис. 2. Геологические составляющие оползня Ю-Е-Линь: a – песчаник; б – глинистый сланец; в – чередующиеся песчаник и глинистый сланец; г – коллювий

Полевые исследования и бурение инженерно-геологических скважин показывают, что поверхностный коллювий оползня Ю-Е-Линь имеет мощность 8,65–35,8 м и состоит из неуплотненных песчаных, пылеватых или глинистых грунтов, смешанных с выветрелыми скальными породами или их обломками, как показано на рисунке 2, г. В выветрелых скальных породах или их обломках плоскости напластования отсутствуют. В результате в общем механическом поведении коллювия основную роль играют мелкозернистые дисперсные грунты.

Коллювиальные отложения можно разделить на глинистые или пылеватые с низкой пластичностью (CL или ML) и на пылеватый песок (SM), смешанный с обломками скальных горных пород. В средней и нижней части склона коллювий подстилается тонким слоем серо-черного глинистого сланца с прослоями мелкозернистого пылеватого песчаника (Sh/Ss), ниже которого залегает коренная порода (у авторов здесь, по-видимому, была ошибка в изложении на английском языке, которая была исправлена переводчиком в соответствии, например, с рисунком 11, a. – Ред.). Считается, что в районе оползня Ю-Е-Линь сошел крупный древний оползень и потенциальная поверхность скольжения (ППС) развивалась вдоль вышеуказанного тонкого слоя глинистого сланца. Как доказательство, в образцах керна этого слоя, взятых при бурении в зоне оползня, часто наблюдаются следы расслаивания и фрагментации. Основная коренная порода здесь представляет собой слоистый песчаник (Ss) серовато-белого цвета мощностью более 25 м.

2.3. Оползень Ю-Е-Линь

В пределах оползня Ю-Е-Линь часто происходило обрушение участков склона из-за неблагоприятных топографических, геологических и гидрологических факторов влияния и плохих дренажных условий. Во время тайфуна «Моракот» с 6 по 9 августа 2009 года сильные ливни (с максимальным суточным количеством осадков 1 080 мм) нанесли серьезный ущерб транспортным линиям, жилым зданиям и общественным объектам, как показано на рисунке 3. Несколько лет спустя после трехдневного дождя с суммарным количеством осадков 629 мм во время тайфуна «Суделор» (с 6 по 9 августа 2015 года) последовала еще одна серия осадков, продолжавшаяся с 9 августа по 1 сентября 2015 года и вызвавшая многочисленные трещины растяжения и оседание на главной дороге округа (автодороге № 166). Ливень во время тайфуна «Мэги» с 25 по 28 сентября 2016 года (с суммарным количеством осадков 776 мм и максимальным суточным количеством осадков 481 мм 26 сентября 2106 г.) снова увеличил существовавшие трещины и вызвал множество мелких обрушений на склоне.

Рис. 3. Ущерб, причиненный тайфунами «Моракот» (8–9  августа 2009 г.) и «Мэги» (27–28 сентября 2016 г.) в результате смещений оползня Ю-Е-Линь:  а – трещина растяжения на главной дороге округа; б – трещина растяжения в здании; в – обрушение склона на главной дороге; г – оседание дороги из-за сдвигового разрушения; д – оседание дороги из-за повреждения подпорной стенки; е – повреждение подпорной стенки

Несмотря на то что движение оползня Ю-Е-Линь относительно медленное и незначительное по сравнению с другими серьезными оползнями, сошедшими на Тайване, длительное смещение грунта на склоне все же вызывает повреждения сооружений, поэтому по-прежнему необходимы работы по его техническому обслуживанию и стабилизации.

2.3.1. Полевой мониторинг

После схода крупного оползня в 2009 году был инициирован план восстановительных работ, который был одобрен Агентством по охране грунтов и вод Совета по сельскому хозяйству Тайваня (SWCB). Эти работы касались грунтовых террас, подпорных стен, дорожных одежд, креплений склона, защитных дамб, неглубоких горизонтальных дрен и поверхностного дренажа. Кроме того, чтобы спроектировать высокоэффективную систему подземного дренажа (два дренажных колодца с горизонтальными дренами), агентство SWCB наладило систематический полевой мониторинг для выяснения механизма, вызывающего смещения оползня. Как показано на рисунке 4, на исследуемой территории были установлены инклинометры и станции мониторинга уровня грунтовых вод (B04-W, B08-W, B09-W, B10-W и B12-W). Жилые зоны (ЖЗ-1, ЖЗ-2 и ЖЗ-3) являлись основными объектами защиты. В этом исследовании для подробного обсуждения были рассмотрены результаты численного анализа дренажных конструкций и их виляние на устойчивость склона вдоль профилей A-A’ и B-B’.

Рис. 4. Точки мониторинга и дренажные системы в исследуемой зоне оползня Ю-Е-Линь

2.3.2. Уровни грунтовых вод

На рисунках 5 а–в показаны результаты измерений количества осадков с помощью самозаписывающего дождемера (рис. 5, г), показывающего подъем/снижение УГВ в ответ на проливные дожди во время тайфунов. Как видно из рисунка 5, a, в точке мониторинга B09-W уровень грунтовых вод постепенно поднялся до самого высокого через 4–6 дней после прихода тайфуна «Дудзюань» («Дуджуан»). Затем УГВ за 10–15 дней постепенно снизился до такого же уровня, какой был перед приходом тайфуна. Запаздывающая реакция уровня грунтовых вод на дожди означает, что общая проницаемость слоя грунта в зоне оползня сравнительно низка и может быть даже ниже, чем средняя интенсивность осадков. Как следствие, просочившаяся в коллювиальные отложения (CL-ML) дождевая вода сохраняется там из-за низкой проницаемости очень мелкозернистого грунта. Эта вода становится неблагоприятным фактором для устойчивости склона и вызывает смещения грунта.

Рис. 6. Профили смещений инклинометров и потенциальная поверхность скольжения по профилю B-B' после тайфуна «Мэги» 25–28 сентября 2016 г. (ЖЗ – жилая зона)

Во время тайфуна «Мэги» (25–28 сентября 2016 г.), как показано на рисунке 5, б, уровень грунтовых вод в точке мониторинга B09-W не сразу среагировал на основные дождевые осадки (вплоть до 1 октября 2016 г.). Это отставание подъема УГВ во времени может быть связано с установкой нескольких групп параллельных неглубоких горизонтальных дрен (см. рис. 4) для подземного дренажа. Сравнивая данные для точек мониторинга B12-W (рис. 5, в) и B09-W (рис. 5, б), можно обнаружить, что реакции УГВ в них во время тайфуна «Мэги» коррелировали между собой и в обоих случаях ответы были с задержкой во времени.

В соответствии с топографическими и геологическими особенностями оползня Ю-Е-Линь, по оценкам, последующие дожди после основных ливней с суммарным количеством осадков более 800 мм могут вызвать локальные смещения грунта.

Результаты мониторинга УГВ во время трех тайфунов сведены в таблицу 1.

Таблица 1. Подъем уровней грунтовых вод (УГВ) во время тайфунов

Рис. 8. Неглубокие подземные дрены в оползне Ю-Е-Линь

2.5. Дренажный колодец для долгосрочной стабилизации оползня

2.5.1. Проектирование дренажного колодца

На Тайване подземный дренаж на сегодняшний день является наиболее часто используемым методом стабилизации склонов – как сам по себе, так и в сочетании с другими методами. Были предприняты попытки оптимизировать количество горизонтальных дрен и расстояний между ними [1, 2, 15]. Расчетная длина горизонтальной дрены должна не только доходить до тела оползня, но и проходить через потенциальную поверхность скольжения на дополнительные 3–5 м [4–6, 16]. Из-за несоответствия стандартам проектирования расчетные расстояния между горизонтальными дренами в литературе [1, 17–20] в целом сильно различаются (от 3 м до 30 м). Креншоу и Санти [21] предложили определять расстояние между горизонтальными дренами в первую очередь в соответствии с условиями дренажа in situ и корректировать его в процессе установки.

Ксантакос с соавторами [22] указали, что естественные склоны редко бывают достаточно однородными, чтобы обеспечить надежное проектирование подземного дренажа в соответствии с простыми принципами осушения. Кроме того, Хаусманн [23] предположил, что для успешной системы осушения проектировщик должен выбрать схему системы дренажа, которая увеличивает вероятность пересечения основного водоносного горизонта. Рахарджо с коллегами [7, 8] доказали, что горизонтальная дрена наиболее эффективна, когда она установлена в нижней части склона. Корнфорт [20] также показал, что горизонтальная дрена, пробуренная в сторону верха склона, может максимизировать площадь пересечения с поверхностью скольжения и эффективно понизить уровень грунтовых вод. При этом Кук с соавторами [24] указали, что горизонтальная дрена может эффективно функционировать, если она может улавливать поток грунтовых вод до того, как он просочится в нестабильную зону оползня. По мнению Линя [13, 14], чтобы рационально спроектировать место строительства дренажного колодца, необходимо заранее понять расположение поверхности скольжения, а также то, как проходят основные пути потоков грунтовых вод и каковы виды их фильтрации в теле оползня. Таким образом, при проектировании системы подземного дренажа в оползне Ю-Е-Линь расположение и высотные отметки дренажных колодцев вместе с подходящей конфигурацией дренажных скважин (или горизонтальных дрен), имеют решающее значение для эффективности дренажа.

Итак, большое количество грунтовых вод может быть отведено из подземной части склона посредством дренажных колодцев, понижающих УГВ. Для удаления грунтовых вод из глубоких слоев грунта в оползне Ю-Е-Линь были спроектированы два ствола колодцев (1 и 2) с дренажными скважинами в соответствии с конфигурацией, показанной на рисунке 9. Их строительство было завершено в 2018 году.

Рис. 9. Типичная конфигурация ствола дренажного колодца с многоуровневыми дренажными скважинами (или горизонтальными дренами) для крупного оползня [25]: а – схематический вертикальный разрез; б – схематический горизонтальный разрез

Построенные дренажные колодцы имеют диаметры по 4 м, а их стволы обсажены гофрированными кольцами из оцинкованной стали. Колодцы 1 и 2 проходят вертикально сквозь коллювий до глубины 20 и 16 м соответственно и находятся близко к поверхности скольжения. Из ствола каждого колодца на четырех различных высотных отметках (с вертикальными шагами 4 и 5 м) было пробурено по группе из 5–6 радиальных необсаженных дренажных скважин диаметрами по 76,2 мм, длиной 50 м, с углами подъема 8–10 град. Эти многоуровневые горизонтальные дрены расположены в веерообразной конфигурации с центральными углами между ними 30 град. и горизонтальными шагами 1,0–2,5 м.

В таблице 2 приведены проектные параметры систем дренажных колодцев 1 и 2.

Таблица 2. Проектные параметры систем дренажных колодцев 1 и 2 в оползне Ю-Е-Линь

Рис. 10. Типичная последовательность строительства ствола дренажного колодца и многоуровневых горизонтальных дрен на Тайване: а – первоначальная выемка грунта с дневной поверхности; б – установка первого сверху обсадного кольца; в – закрепление входа в колодец бетоном после установки первого обсадного кольца; г – поэтапная выемка грунта и установка очередных обсадных колец; д – установка стальных обручей для закрепления стальных колец; е – бурение скважины для установки горизонтальной дрены; ж – обматывание горизонтальной дрены фильтром; з – трубы из ПВХ для отвода собранных грунтовых вод в водобойный колодец на дне ствола дренажного колодца, лестница для технического обслуживания и ремонта; и – крышка дренажного колодца и защитное ограждение вокруг

На Тайване обычно проектируются дренажные колодцы диаметром 3,5 или 4,0 м и глубиной 15–30 м или более. Дренажный колодец не должен проходить через поверхность скольжения, чтобы избежать сдвигового разрушения его ствола из-за движения грунтовой массы оползня вниз по склону. Согласно опыту прошлых работ, если дренирующая способность труб из ПВХ в нижней части ствола колодца недостаточна для приема воды, поступающей из горизонтальных дрен, включаются насосы для откачки воды в водобойный колодец.

3. МЕТОДОЛОГИЯ

С помощью трехмерного численного моделирования в этом исследовании рассматривались эффективность подземных дренажных систем и их влияние на снижение уровня грунтовых вод (hw,max) и на соответствующие коэффициенты устойчивости (Ку) потенциальной поверхности скольжения оползня Ю-Е-Линь.

Расчеты проводились поэтапно: сначала были определены природные напряжения (in situ), смоделировано строительство двух дренажных колодцев, а после этого – выполнена серия расчетов нестационарной фильтрации совместно с определением напряженно-деформированного состояния (fully coupled flow-deformation), на основе чего был определен коэффициент устойчивости с использованием метода снижения прочности (итерационного снижения прочностных характеристик грунта – угла внутреннего трения и удельного сцепления). Последнее выполнялось для расчета значения Ку в ответ на понижение уровня грунтовых вод из-за функционирования горизонтальных дрен. Следует подчеркнуть, что более реалистичное и более высокое значение Ку будет получено именно при полностью совмещенном анализе фильтрационных потоков и деформаций при учете поглощения влаги водоненасыщенным грунтом (в PLAXIS выполнение такого расчета возможно в том числе с учетом инфильтрации и испарения. – Ред.).

3.1. Численная модель с граничными и начальными условиями

На основе проектных параметров дренажных колодцев (см. таблицу 2), топографических, геологических и гидрологических данных длительных полевых исследований и испытаний дисперсных и скальных грунтов, отобранных из скважин на исследуемой территории, была получена трехмерная численная модель оползня Ю-Е-Линь, показанная на рисунке 11, a. Слои грунта моделируются 10-узловыми тетраэдрическими элементами (или 10-узловыми объемными элементами в программном комплексе PLAXIS 3D [26]) с четырьмя точками численного интегрирования. Эти элементы имеют три степени свободы на узел (или три компонента смещения на узел) и обеспечивают интерполяцию смещений второго порядка (рис. 11, б). Дренажный колодец состоит из ствола и многоуровневых горизонтальных дрен. Ствол колодца, непрерывно обсаженный гофрированными кольцами из оцинкованной стали, моделируется 6-узловыми треугольными плитными элементами с тремя точками численного интегрирования, обеспечивающими интерполяцию смещений второго порядка (рис. 11, в). Поскольку плитные элементы не могут выдерживать крутящие моменты, они имеют только пять степеней свободы на узел (три поступательных и две вращательных). Плитные элементы – это специальные конечные элементы, используемые для моделирования тонкой двумерной конструкции стенки ствола в грунте со значительной жесткостью на изгиб и основанные на теории пластин Миндлина [26], которая допускает смещения и изгибы пластины из-за сдвига. Структурные силы оцениваются в точках интеграции и экстраполируются на узлы элементов. Наконец, горизонтальные дрены моделируются серией 3-узловых линейных дренажных элементов с гидравлическими условиями (см. рис. 11, в), которые совместимы со сторонами 6-узловых треугольных плитных элементов и 10-узловых тетраэдрических элементов грунта, поскольку эти элементы также имеют три узла на стороне. Линейные дренажные элементы с 3-мя узлами обеспечивают интерполяцию второго порядка (квадратичную) расхода подземных вод в горизонтальных дренах.

Рис. 11. Численная модель оползня Ю-Е-Линь с системой подземного дренажа (a); вид сетки конечных элементов и системы подземного дренажа внутри численной модели (б); смоделированные стволы дренажных колодцев с горизонтальными дренами (в)

В идеале границы модели располагаются вдоль естественных гидрологических границ, таких как водоразделы (границы водосборных бассейнов), водоемы (ручьи) и непроницаемая коренная порода. Расчетная схема (на рисунке 12) выбрана так, чтобы охватить границу водораздела для уменьшения количества элементов и сокращения времени расчета. Для расчета фильтрационных потоков грунтовых вод требуются местоположение, мощность и гидравлические свойства каждого слоя грунта. Начальный УГВ для оползня определяется и устанавливается путем расчета установившегося потока грунтовых вод, который калибруется с помощью измеренных величин УГВ в наблюдательных скважинах B08-W, B09-W, B10-W и B12-W (см. рис. 4). Для расчета установившегося потока грунтовых вод и совместного анализа фильтрации и деформаций, за исключением нижней границы (плоскости X-Y), которая задана как граница нулевого потока, все окружающие границы (плоскости X-Z и Y-Z) определяются как открытые границы с постоянным напором. Кроме того, для мониторинга результатов расчетов выбирается профиль A-A’, вдоль которого сооружены два дренажных колодца, как показано на рисунке 12.

Рис. 12. Размеры численной модели и расположение защищаемых жилых зон (ЖЗ-1, ЖЗ-2), дренажных колодцев 1 и 2 и профилей A-A', C1-C1', C2-C2' для мониторинга результатов численных расчетов

3.2. Моделирование дренажа для дренажного колодца и горизонтальных дрен

При анализе устойчивости склона на контактах между стволом дренажного колодца (плитными элементами) и окружающим грунтом добавляются интерфейсные элементы, чтобы обеспечить корректное моделирование взаимодействий «грунт – конструкция». По умолчанию в программе PLAXIS 3D плитные элементы полностью водопроницаемы, поэтому ствол колодца также становится проницаемым. Однако на Тайване, чтобы увеличить сопротивление ствола колодца сдвигу и моменту при действии направленной вниз склона силы, строительный зазор между стволом колодца и окружающим грунтом часто заполняется чистым (без примесей) бетоном, что делает ствол водонепроницаемым. Следовательно, при численном моделировании интерфейсные элементы, относящиеся к стволу дренажного колодца, устанавливаются полностью непроницаемыми, чтобы блокировать потоки влаги из окружающих слоев грунта.

Как уже упоминалось, горизонтальный дренаж в численной модели имитируется линейными дренажными элементами (линейными дренами). В местах расположения горизонтальных дрен накладывается серия линейных дрен для управления распределением порового давления при совместном анализе фильтрационных потоков и деформаций. Линейные дрены используются для задания линий внутри геометрической модели в местах, где снижается поровое давление. При расчете поровое давление во всех узлах дренажной линии уменьшается до заданного напора, эквивалентного нулю для гравитационного (самотечного) режима дренажа.

3.3. Параметры материалов в модели

Для определения прочностных параметров была проведена серия испытаний образцов дисперсных и скальных грунтов на прямой сдвиг. Образцы грунтов, отобранные из коллювиальных отложений, классифицируются как глинистые или пылеватые грунты с низкой пластичностью (CL, ML, CL-ML). Остаточные удельные сцепления cres и остаточные углы внутреннего трения φres находятся в диапазоне 15–35 кПа и 26–30 град. соответственно. Кроме того, образцы, отобранные из скальной коренной породы, залегающей ниже коллювия, классифицируются как глинистые сланцы (Sh), а параметры их остаточной прочности (cres, φres) составляют 0 кПа и 23–32 град. соответственно (прочность на одноосное сжатие qu = 3 060 кПа). Для песчаника (Ss) максимальные параметры прочности (cpeak, φpeak) равны 0 кПа и 32–35 град. соответственно, а остаточные (cres, φres) составляют 0 кПа и 25–33 град. соответственно (прочность на одноосное сжатие qu = 4 000 кПа). Прочность на сдвиг границы раздела между песчаником и сланцем в этом исследовании не проверялась. При анализе используется остаточная прочность коллювия и подстилающего его глинистого сланца, чтобы включить снижение прочности из-за большого смещения грунта в полевых условиях (в РФ остаточную прочность определяют методом «плашка по плашке» и также используют ее в расчетах оползней. – Ред.). Напротив, песчаник рассматривался как коренная порода, претерпевающая лишь незначительные сдвиги, и его свойства определялись максимальной прочностью на сдвиг. Линейно-упругая идеально-пластическая модель Мора – Кулона с критериями разрушения Мора – Кулона используется для моделирования механического поведения слоев грунта в оползне Ю-Е-Линь. Входные параметры модели Мора – Кулона приведены в таблице 3. Параметры линейной изотропной модели такого материала, как обсадное кольцо ствола дренажного колодца из гофрированной оцинкованной стали, использованные для численного моделирования, перечислены в таблице 4.

Таблица 3. Входные параметры модели Мора  Кулона для слоев грунта при анализе устойчивости склона

Рис. 14. Сравнение величин снижения УГВ по данным для наблюдательной скважины B12-W с подземным дренажем и без него для разных коэффициентов фильтрации коллювиальных отложений

В таблице 6 приведены изменения глубины залегания подземных вод от дневной поверхности HW(t) в наблюдательной скважине B12-W и коэффициента устойчивости КУ(t) оползня во времени с учетом работы подземного дренажа. Для трех значений коэффициента фильтрации коллювия KS тенденции изменений HW(t) и КУ(t) аналогичны, поэтому для краткости представлены только результаты расчетов для KS = 1 × 10–4 см/с. Видно, что коэффициент устойчивости КУ(t) неуклонно растет при понижении уровня грунтовых вод.

Таблица 6. Изменения снижения УГВ и коэффициента устойчивости во времени при подземном дренаже и коэффициенте фильтрации коллювия KS = 1 × 10–4 см/с

Рис. 15. Механизм развития оползня в исследуемой зоне без подземного дренажа в сухой сезон: а – поле полных приращений смещений; б – направление скольжения на основе данных обследования поверхности; в – потенциальная оползневая зона; г – потенциальная поверхность скольжения (ППС) по профилю В-В'

4.3. Сравнение потенциальных поверхностей скольжения при проведении подземного дренажа и без него

Рисунок 16 демонстрирует эффективность подземного дренажа для улучшения устойчивости склона. Потенциальная оползневая зона (рис. 16, a-1) без дренажа распространяется на большую площадь вниз по склону, и соответствующая потенциальная поверхность скольжения вдоль профиля A-A’ (см. рис. 4) также находится глубоко – вблизи уровня дна дренажных колодцев (рис. 16, a-2). Во время сезона тайфунов исследуемая территория находится в состоянии нестабильности с более низким коэффициентом устойчивости КУ(t=0 сут) = 1,383. Коэффициент устойчивости КУ < 1 означает, что склон находится в состоянии разрушения, и было решено, что для обеспечения долгосрочной стабильности в результате работ по ее сохранению и повышению требуется КУ = 1,5 [28]. На Тайване в качестве технических критериев для инженерного проектирования стабилизации склонов были приняты три значения КУ:

  • для обычного времени КУ 1,50;
  • в случае землетрясения КУ 1,2;
  • в случае проливных дождей КУ 1,1.

Более того, Попеску [29] предложил трехпериодный непрерывный спектр КУ для определения состояния устойчивости склонов:

  • КУ > 1,3 (стабильный);
  • 1,0 < КУ <1,3 (минимально стабильный);
  • КУ < 1,0 (активно нестабильный).

Напротив, после подземного дренажа и снижения уровня грунтовых вод потенциальная оползневая зона (см. рис. 16, б-1) ограничивается участком ниже по склону, а потенциальная поверхность скольжения располагается на сравнительно небольшой глубине (см. рис. 16, б-2). В конечном итоге склон сохраняется в более стабильном состоянии с более высоким коэффициентом устойчивости КУ(t=20 сут) = 1,521.

Рис. 16. Потенциальная оползневая зона и потенциальная поверхность скольжения вдоль профиля A-A' в сезон тайфунов без подземного дренажа (а-1, а-2) и с подземным дренажем (б-1, б-2)

4.4. Эффективность расходов воды в горизонтальных дренах

На рисунке 17 показаны расходы для систем двух дренажных колодцев в сезон тайфунов. Индексы w1, w2 и d1–d4 относятся к дренажным колодцам 1; 2 и к уровням горизонтальных дрен 1–4 (сверху вниз) соответственно. В сезон тайфунов средний начальный уровень грунтовых вод перед подземным дренажем составляет приблизительно 11 и 5 м ниже поверхности земли для дренажных колодцев 1 и 2 соответственно. Для колодца 1 горизонтальные дрены 1-го и 2-го уровней размещаются соответственно на 5 и 10 м ниже дневной поверхности, что выше среднего начального уровня грунтовых вод (11 м ниже поверхности земли). При таких обстоятельствах расход в горизонтальных дренах 1-го и 2-го уровней (Qw1d1 и Qw1d2) становится нулевым и неэффективным для снижения УГВ, как показано на рисунке 17, a. В горизонтальных дренах 1-го и 2-го уровней системы дренажного колодца 2 (Qw2d1 и Qw2d2) наблюдается аналогичная ситуация (рис. 17, б).

Рис. 18. Корректировка местоположения дренажных колодцев (места контакта горизонтальных дрен показаны кружками, скорректированные местоположения колодцев показаны заштрихованными прямоугольниками)

Если бы положение дренажного колодца 1 было изменено на скорректированное, то зона эффективного снижения УГВ выше по склону была бы расположена в пределах основного смещающегося тела оползня и скомпенсировала бы зону неэффективного снижения УГВ ниже по склону. С другой стороны, если бы положение дренажного колодца 2 было изменено на скорректированное, то зона эффективного снижения УГВ выше по склону могла бы охватить более нестабильную область в пределах основного смещающегося тела оползня. Как следствие, общая эффективность работы системы подземного дренажа в оползне могла бы быть значительно увеличена.

5. ВЫВОДЫ

На основе полевого мониторинга и совместного фильтрационно-деформационного трехмерного численного анализа проверялась пригодность двух глубоких дренажных колодцев с системами многоуровневых горизонтальных дрен для снижения уровня грунтовых вод и стабилизации склона. Устойчивость склона оценивалась по развитию потенциальной зоны скольжения и соответствующего коэффициента устойчивости с использованием метода итерационного снижения прочности. Два глубоких дренажных колодца с многоуровневыми горизонтальными дренами оказались эффективными для повышения устойчивости склона и снижения вероятности схода оползня. По результатам численных расчетов можно сделать следующие выводы.

1. Эффективность численной модели, учитывающей взаимное влияние фильтрационных потоков и деформаций, для расчета установившегося и переходного потоков грунтовых вод и анализа устойчивости склона с использованием метода итерационного снижения прочности была подтверждена:

  • калибровкой начальной глубины залегания подземных вод от поверхности земли в точке мониторинга B12-W (минус 13 м по результатам моделирования и минус 12,97 м по итогам полевых наблюдений);
  • направлением движения грунта оползня (к юго-востоку по результатам и моделирования, и наблюдений);
  • глубиной потенциальной поверхности скольжения вдоль профиля B-B’ (минус 11,4 м для точки B12-W и минус 35,6 м для B09-W по итогам моделирования; минус 12 м для B12-W и минус 36 м для B09-W по результатам наблюдений).

2. В целом, проектная дина дрены Ld должна проходить через потенциальную поверхность скольжения, а переходная длина Lt1 должна располагаться в основной части ППС. Однако в этом исследовании длина Lt1 выходит далеко за пределы ППС, что уменьшает эффективность снижения уровня грунтовых вод для обеспечения стабильности потенциальной поверхности скольжения.

3. Для повышения эффективности подземного дренажа дренажные колодцы 1 и 2 можно переместить вниз по склону на горизонтальные расстояния соответственно 24 и 10,5 м. После такой корректировки снижение УГВ и стабилизация склона становятся более эффективными (по отношению к ППС) из-за того, что длина Lt1 оказывается внутри основного смещающегося тела оползня. Кроме того, уровень грунтовых вод ниже по склону от дренажного колодца 1 может быть эффективно понижен колодцем 2 из-за уменьшения расстояния между этими двумя колодцами.

4. В конечном итоге дренажный колодец может дать более хорошую эффективность дренажа и более высокую скорость понижения УГВ, если соответствующие горизонтальные дрены в основном расположены в средней зоне потенциальной поверхности скольжения и проходят через большую область массива грунта внутри ППС, чтобы обеспечить фильтрационные потоки грунтовых вод в дрены.

5. Для горизонтальных дрен, расположенных выше уровня грунтовых вод, подземный дренаж становится неэффективным в отношении снижения УГВ. Для систем обоих дренажных колодцев самые нижние горизонтальные дрены имеют наибольшие расходы. Кроме того, скорость снижения УГВ (эффективность дренажа) подземной дренажной системой может быть повышена за счет увеличения коэффициента фильтрации KS (водопроницаемости) коллювия, однако это может стать неэкономичным, если коллювий имеет более низкую проницаемость

(KS < 1 × 10-5 см / с).

6. И наконец, следует отметить, что подземный дренаж с высокими инженерными затратами на рассматриваемой территории остается необходимым из-за невозможности перемещения трех существующих жилых зон.

Это исследование было финансово поддержано Агентством по охране грунтов и вод Совета по сельскому хозяйству Тайваня (SWCB).

24 Июнь 2021