Полный текст

Введение Очистные сооружения промышленного предприятия расположены на грунтовом массиве со сложным рельефом, который представляет собой коренной склон речной долины. Площадка частично застроена зданиями и сооружениями производственного назначения с большим количеством подземных и надземных коммуникаций. Физико-геологические и техногенные процессы, опасные для эксплуатации сооружений, в том числе, связанные с процессами развития карста (поверхностные формы карстопроявления: воронки, локальные оседания) и развитием склоновых, процессов (оползни, обвалы, осыпи) на площадке не обнаружены. Перепад высот по всей площадке составляет 32 м. Непосредственно по территории обследуемого объекта перепад высот составил 12 м. В геологическом строении района территории объекта до глубины 15,0 м по данным бурения инженерно-геологических скважин принимают участие отложения двух систем - четвертичной (Q) и каменноугольной (С). Четвертичные отложения представлены техногенными (tQiv) насыпными и аллювиально-делювиальными (adQ) суглинками мягкопластичными. Каменноугольная система представлена элювиальными (еС1) суглинками дресвяными полутвердыми, суглинками полутвердыми и тугопластичными и коренными отложениями визейского яруса нижнего отдела (С1), представленными алевролитом низкой прочности и известняком средней прочности, породы размягчаемые. Поверхность территории объекта участками покрыта почвенно-растительным слоем мощностью 0,2 м. Степень размокаемости грунтов - медленная, очень медленная и неразмокающий грунт. Насыпной грунт суглинок дресвяный полутвердый ИГЭ-1, насыпной грунт суглинок дресвяный тугопластичный ИГЭ-1б, обладают средней степенью влажности Sr > 0,9. Изыскательская организация рекомендовала принять их как сильнопучинистые грунты. Насыпной грунт песок средней крупности ИГЭ-1в и насыпной гравийный грунт с песчаным заполнителем ИГЭ-1г относятся к непучинистым грунтам при любом положении уровня подземных вод. Показатели физико-механических свойств грунтов, которые определяют состав инженерно-геологических элементов, определенные в результате стандартных испытаний, полученных при бурении скважин образцов, а также нормативные и расчетные значения прочностных характеристик грунтов, приведены в табл. 1. Инженерно-геологические и гидрогеологические условия территории объекта по совокупности факторов, определенных в ходе изысканий характеризуются 3-й категорией сложности Объектом численного исследования является армогрунтовая насыпь высотой до 12 м с внешней стороны зданий и сооружений, выполненная для упорядочивания рельефа территории, размещения инженерных сетей и автопроезда. Инженерно-геологический разрез склона с обозначением инженерно-геологических элементов, свойства которых представлены в табл. 1, изображен ниже. Таблица 1 Показатели физико-механических свойств грунтов по результатам изысканий Table 1 Indicators of physical and mechanical properties of soils on the results of surveys Наименование ИГЭ Нормативные значения характеристик грунтов Расчетные значения характеристик грунтов Предел прочности на одноосное сжатие, Rc, МПа плотность грунта, г/см3 удельное сцепление, МПа угол внутреннего трения, град для расчетов по несущей способности коэффициент надежности по грунту плотность грунта, г/см3 удельное сцепление, МПа угол внутреннего трения, град Насыпной грунт: суглинок дресвяный полутвердый (tQiv), ИГЭ-1 2,05 0,028 14 gg (с) = 1,271 gg (j) = 1,043 2,02 0,022 13 - Насыпной грунт: валуны, глыбы, щебень (tQiv), ИГЭ-1а - - - - - - - - Насыпной грунт: суглинок дресвяный тугопластичный (tQiv), ИГЭ-1б 1,95 0,022 20 gg (с) = 1,198 gg (j) = 1,095 1,91 0,018 18 - Насыпной грунт: песок средней крупности (tQiv), ИГЭ-1в 1,93 0,002 38 gg (с) = 1,144 gg (j) = 1,012 1,91 0,002 38 - Насыпной гравийный грунт с песчаным заполнителем (tQiv), ИГЭ-1г 2,09 0,003 37 gg (с) = 1,300 gg (j) = 1,059 2,07 0,002 35 - Суглинок мягкопластичный (adQ), ИГЭ-2б 1,87 0,016 15 gg (с) = 1,170 gg (j) = 1,127 1,85 0,013 13 - Суглинок дресвяный полутвердый (еС1), ИГЭ-3 1,92 0,022 21 gg (с) = 1,074 gg (j) = 1,060 1,89 0,020 20 - Суглинок полутвердый (еС1), ИГЭ-4а 2,04 0,028 22 gg (с) = 1,089 gg (j) = 1,048 2,02 0,026 21 - Суглинок тугопластичный (еС1), ИГЭ-4б 1,97 0,020 19 gg (с) = 1,159 gg (j) = 1,097 1,94 0,018 17 - Алевролит низкой прочности, размягчаемый (С1), ИГЭ-5а 2,12 - - - 2,09 - - 2,1* 1,7 Известняк средней прочности, размягчаемый (С1), ИГЭ-8 2,58 - - - 2,56 - - 22,5* 17,0 Примечание: * - предел прочности на одноосное сжатие (Rс) образцов приведен в виде дроби: над чертой - в естественном состоянии, под чертой - в водонасыщенном состоянии. Рис. 1. Инженерно-геологических разрез исследуемого объекта Fig. 1. Engineering-geological section of the object under study Метод вычисление величины коэффициента запаса устойчивости грунтового массива Для вычисления величины коэффициента запаса устойчивости откоса использованы компьютерные программы [1, 2], в которых формализован метод расчета устойчивости откосов В.К. Цветкова - А.Н. Богомолова [3, 4]. Данный метод основан на анализе напряженного состояния приоткосной области либо методами теории функций комплексного переменного [5] для однородных, либо методом конечных элементов [6-8] для неоднородных откосов. Он включает в себя формулы для определения величины коэффициента запаса устойчивости в точке грунтового массива Kт и глобального коэффициента запаса устойчивости K и методику построения наиболее вероятной линии скольжения. Запишем условие прочности Кулона - Мора [9, 10] в виде, предложенном А. Како [11], введя предварительно некую функцию напряженного состояния и физико-механических свойств грунта Kт в рассматриваемой точке грунтового массива (1) где - приведенное нормаль ное напряжение [12]; - приведенное давление связности, где: С; φ; γ и H - соответственно удельное сцепление, угол внутреннего трения, удельный вес грунта и высота откоса. При Km = 1 выражение (1) принимает вид условия прочности Кулона. Рассмотрим равновесие бесконечно малой призмы единичной высоты (рис. 2), расположенной в некоторой точке грунтового массива, и составим уравнения равновесия проекций всех сил на направления (2) Рис. 2. Расчетные схемы для составления уравнений равновесия бесконечно малой треугольной призмы Fig. 2. Calculation schemes for composing the equilibrium equations of an infinitesimal triangular prism Затем выразим нормальную sn и касательную tn компоненты напряжения, действующего по наклонной площадке, ориентированной под углом a, через составляющие sz; sx; txz полного напряжения в рассматриваемой точке (3) Подставляя выражения (3) в формулу (1) и выделяя величину Km, получим выражение для вычисления величины коэффициента запаса устойчивости в точке грунтового массива . (4) Числитель и знаменатель полученной дроби равны соответственно удерживающей Fуд и сдвигающей Fсд силам, действующим в рассматриваемой точке по площадке сдвига, ориентированной под углом a (см. рис. 1). Задача об определении угла ориентации наиболее вероятной площадки сдвига решается путем исследования функции (4) на экстремум согласно тривиальной процедуре, заключающейся в выполнении условий (5) Выполнение первого из условий (5) сводится к решению квадратного уравнения, корни которого определяются выражением , (6) где ; (7) (8) Выбирать из полученных двух углов следует тот, при подстановке значения которого во второе условие (5) будет выполняться неравенство . (9) В результате вычислений установлено, что величина угла a, при котором K = Kmin соответствует знаку «+» перед радикалом в формуле (6). Построение наиболее вероятной линии скольжения (точнее говоря, следа наиболее вероятной поверхности скольжения) проводится согласно методике, основанной на предложениях, приведенных в работе [12], и усовершенствованной в работах [3-5]. Согласно этой методике, которая подробно изложена в работе [3], при построении НВЛС в каждой ее точке выполняется условие . Если S - дуговая координата на наиболее вероятной линии скольжения длиной l, то величина коэффициента запаса устойчивости грунтового массива (естественно при неизменном положении НВЛС) определяется выражением , (10) где Fуд и Fсд - удерживающие и сдвигающие силы, действующие в каждой точке НВЛС, численные значения которых определяются соответственно числителем и знаменателем формулы (4). Результаты вычислений Анализ полей напряжений Анализ напряженного состояния исследуемого склона выполнен, как отмечено выше, при помощи компьютерных программ [1, 2], разработанных в Волгоградском государственном архитектурно-строительном университете. На рис. 3 представлена расчетная схема метода конечных элементов для исследуемого геологического разреза до подрезки и устройства автодороги: дискретизация расчетной области на конечные элементы (а) и геологическая структура склона до и после подрезки и устройства автодороги - (б и в соответственно). Отметим, что расчетная схема составлена с учетом результатов топографической съемки склона левее скважины С7 (см. рис. 1). Рис. 3. Расчетная схема метода конечных элементов для исследуемого геологического разреза до подрезки и устройства автодороги: дискретизация расчетной области на конечные элементы (а) и геологическая структура склона до и после подрезки и устройства автодороги - (б) и (в) соответственно Fig. 3. Calculation scheme of the finite element method for the geological section under study before pruning and road construction: discretization of the calculated area into finite elements (a) and the geological structure of the slope before and after pruning and road construction - (б) and (в) respectively Первая расчетная схема (рис. 2, а) состоит из 4443 конечных элементов, сопряжённых в 2305 узлах. Ширина матрицы жесткости системы линейных уравнений равна 98. Вторая расчетная схема отличается от первой углом наклона участка AB к горизонту, большей степенью дискретизации (частота сетки конечных элементов (КЭ)) на участке, покрываемом отрезком AB (визуально эти отличия при данном масштабе трудно заметить), и наличием внешней нагрузки от автомобильной дороги (см. рис. 2, в). Эта расчетная схема состоит из 4772 конечных элементов, сопряженных в 2479 узлах. Ширина матрицы жесткости системы линейных уравнений равна 102. Граничные условия отрабатывались аналогично тому, как это сделано в работах [13, 14]. Численные значения физико-механических свойств слагающих горных пород приняты согласно результатам инженерно-геологических изысканий. Величина коэффициента бокового давления в обводненных породах принята равной ξо = 0,75, как это рекомендовано в [15-17], а в остальных случаях вычислялась по формуле Н.М. Герсеванова ξо = tq2 (45 - φ / 2). (11) На рис. 4 и 5 приведены картины изолиний вертикальных (а), горизонтальных (б) и касательных (в) напряжений в рассматриваемом геологическом разрезе до и после подрезки склона и устройства автодороги соответственно. Их численные значения даны в технической системе единиц с размерностью [т/м2]. Рис. 4. Изолинии вертикальных σz (а), горизоньальных σx (б) и касательных τxz (в) напряжений в исследуемом геологическом разрезе до подрезки склона и устройства автодороги Fig. 4. Isolines of vertical σz (a), horizontal σx (б) and tangential τxz (в) stresses in the studied geological section before slope trimming and road construction Рис. 5. Изолинии вертикальных (а), горизонтальных (б) и касательных (в) напряжений в исследуемом геологическом разрезе после подрезки склона и устройства автодороги Fig. 5. Isolines of vertical σz (a), horizontal σx (б) and tangential τxz (в) stresses in the studied geological section after cutting the slope and road construction Анализ картин изолиний показывает, что во втором случае при наличии нагрузки от автодороги последняя «подтянула» к себе изолинии вертикальных напряжений σz - непосредственно под нагрузкой возникла их концентрация. В остальных локациях в обоих случаях поле вертикальных напряжений однородно, форма изолиний практически повторяет форму контура откоса, что объясняется небольшим отличием величин удельного веса пород, составляющих ИГЭ. Лишь вблизи некоторых контактов слоев отмечается неровность изолиний, что связано с достаточным отличием численных значений модулей деформации Ео и коэффициентов бокового давления ξо соприкасающихся породных напластований. Изолинии горизонтальных напряжений σx практически не изменили своей конфигурации, однако зона растягивающих горизонтальных напряжений, находящаяся непосредственно под автодорогой, уменьшилась в размере примерно на 35 %, при этом положение и форма зоны растягивающих горизонтальных напряжений, расположенной справа от аэротенка, практически не изменились. Касательные напряжения τzx в разных локациях имеют различные знаки, их изолинии в обоих случаях подобны, однако, наличие нагрузки от автомобильной дороги повлекло за собой существенное повышение численных значений этих напряжений в соответствующих точках расчетной области. Из вышесказанного видно, что численные значения всех трех компонент напряжения в соответствующих точках расчетной отличаются друг от друга, что обусловит различные значения коэффициентов запаса устойчивости. Коэффициенты запаса устойчивости На рис. 6 и 7 изображены по три наиболее вероятные линии скольжения, построенные из различных точек приоткосной области. Первая из них (рис. 6, а и 7, а) определяет величину коэффициента запаса устойчивости верхнего уступа, вторая (рис. 6, б и 7, б) - двух верхних уступов, а третья (рис. 6, в и 7, в) - величину глобального коэффициента запаса. Рис. 6. Наиболее вероятные линии скольжения в исследуемом сечении при отсутствии подрезки и нагрузки от автодороги Fig. 6. The most probable sliding lines in the section under study in the absence of trimming and loading from the highway Рис. 7. Наиболее вероятные линии скольжения в исследуемом сечении после подрезки и приложения нагрузки от автодороги Fig. 7. The most probable sliding lines in the section under study after trimming and applying the load from the highway Подчеркнем, что первые, вторые и третьи НВЛС имеют одинаковые координаты точек выхода. Перемещения и осадки Из приведенных выше рис. 6 и 7 видно, что под левым краем аэротенка справа и слева от него имеются достаточно развитые области пластических деформаций (окрашены синим цветом), которые могут спровоцировать неравномерные осадки его основания и, как следствие, возникновение трещин в его железобетонной конструкции. Небольшая по размеру ОПД, расположенная в основании откоса, не оказывает влияние на его глобальную устойчивость. Однако для предотвращения местных вывалов и осыпания породы необходимо провести поверхностное локальное укрепление откоса геосинтетическим материалом. Рис. 8. Изолинии вертикальных (а) и горизонтальных (б) перемещений в теле грунтового масива Fig. 8. Isolines of vertical (a) and horizontal (б) displacements in the body of the soil massif На рис. 8 приведены изображения изолиний вертикальных и горизонтальных перемещений в теле исследуемого грунтового массива (шаг изолиний 0,05 м). «Луковица» вертикальных напряжений практически симметрична относительно вертикальной оси. Картина изолиний горизонтальных перемещений более сложная: приоткосная область практически по вертикали (см. красную «нулевую» линию горизонтальных напряжений на рис. 8, б) разделена на «сжатую» и «растянутую» зоны, последняя из которых расположена слева от «нулевой» изолинии непосредственно на склоне. Данное обстоятельство повышает степень риска разрушения поверхности склона. Выводы 1. Глобальная устойчивость склона обеспечена. 2. Устойчивость грунтовой насыпи (верхнего уступа склона) недостаточна. Необходимо проведение инженерных мероприятий по обеспечению его устойчивости. 3. Учитывая 3-ю категорию сложности инженерно-геологических условий и то, что разрушение грунтового массива может привести к нанесению значительного ущерба окружающей среде, необходимо организовать постоянный комплексный мониторинг грунтового массива и расположеннных на нем сооружений. По результатам мониторинга рекомендуется, в случае необходимости, принять решение об укреплении грунта основания аэротенка в местах локации областей пластических деформаций (слева и справа от него, под ним) во избежание неравномерных осадок основания и кренов сооружения. 4. По всей площади поверхности склона в местах локализации пластических деформаций (поверхности обоих уступов) необходимо провести инженерные (компенсирующие) мероприятия по предотвращению разрушения поверхностных слоев.

Об авторах

А. Н Богомолов

Научно-проектное экспертно-консультационное предприятие «ОиФ»

В. Г Офрихтер

Пермский национальный исследовательский политехнический университет

А. В Редин

Черномортранснефть

О. А Богомолова

Волгоградский государственный технический университет

С. А Богомолов

Российский государственный аграрный университет - Московская сельскохозяйственная академия имени К.А. Тимирязева

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

Статистика

Просмотры

Аннотация - 150

PDF (Russian) - 118

Ссылки

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