CC BY
15
- © A.M. Гальперин, B.B. Ческидов, Ю.В. Бородина,
A.B. Демидов, 2014
УДК 504.064.2
A.M. Гальперин, В.В. Ческидов, Ю.В. Бородина, A.B. Демидов
ИЗУЧЕНИЕ НЕСТАЦИОНАРНЫХ ГИДРОГЕОМЕХАНИЧЕСКИХ ПРОЦЕССОВ В ГЛИНИСТЫХ ПОРОДАХ
Освещаются вопросы оценки деформируемости, прочности и ползучести глинистых горных пород на примере объектов КМА и района строительства олимпийских объектов близ г.Сочи. Рассмотрены возможности применения дистанционного контроля геомеханических процессов при ведении горных и строительных работ, в том числе с использованием аэро- и космосъемки и GPS технологий. Ключевые слова: горные работы, гидрогеомеханические процессы, массивы горных пород, связные (глинистые) грунты, дистанционный контроль.
Изучению реологических свойств глинистых пород посвящены работы H.H. Цытовича, С.С. Вялова, М.Н. Гольдштейна, H.H. Маслова, Г.И. Тер-Степаняна, С.Р. Месчана, Г.Л. Фисенко, Ю.К. Зарецкого, З.Г. Тер-Мартиросяна, Тан Тьонг Ки, Л. Шук-ле, А. Скемтона и др. [1, 2, 3, 4].
Ползучесть горных пород рассматривалась этими исследователями с позиций механики сплошных сред. Подход к породе как к квазиоднофазной системе с использованием классической теории ползучести допустим при выполнении условия о независимости характера изменения напряженно-деформированного состояния от геометрических размеров (образцов или зон в пределах массива).
Для многофазных глинистых пород процесс деформирования определяется физико-механическими свойствами отдельных фаз и характером их взаимодействия. В этом случае интенсивность изменения соотношения фаз в единице объема зависит от геометрических размеров деформируемой области. Виды уравнений состояния определяются для каждой фазы отдельно (в условиях их невзаимодействия). Расчетные модели многофазных глинистых пород реа-
лизуются в решениях задач консолидации грунтов (К. Терцаги, Н.М. Гер-севанов, В.А. Флорин, Н.А. Цыто-вич, Д.Е. Польшин, В. Г. Короткин, М.Н. Гольдштейн, Н.Н. Веригин, Ю.К. Зарецкий, А.И. Ксенофонтов, З.Г. Тер-Мартиросян, М. Био, Н. Ка-рилло, Р. Гибсон, Д. Тейлор, М.В. Малышев, М.Ю. Абелев и др.).
Исследования устойчивости бортов карьеров с учетом фактора времени проводились Г.Л. Фисенко, В. И. Веселковым, А.М. Мочаловым, Н.Н. Куваевым, А.Н. Могилко и др. Оценке изменения во времени состояния пород насыпных и гидравлических отвалов посвящены работы МГГУ, ВНИМИ, НИИКМА, УкрНИИ Проекта, ИГД им. А.А. Скочинско-го. Фактор времени учитывается при оценке устойчивости и уплотнения пород бортовых и отвальных массивов с помощью решений задач структурной механики грунтов, фильтрационной консолидации и теории наследственной ползучести. Использование указанных теоретических решений позволяет осуществлять геомеханическое обоснование порядка выполнения вскрышных и отвальных работ. Общим для всех методов прогноза состояния карьерных откосов во време-
ни является нестационарным характер протекающих в массивах горных пород механических процессов.
Разработанными под руководством В. В. Ржевского научными основами управления карьерными откосами предусматривается гибкое реагирование на изменение состояния массивов пород. Усовершенствование методов прогноза нестационарных геомеханических процессов на карьерах позволяет, исходя из масштаба и характера изменений прочности бортовых и отвальных массивов во времени, регламентировать инженерные мероприятия по управлению состоянием пород [3, 4].
К группе пород, характеризующихся высокой степенью ползучести и пластичности, Г. Л. Фисенко относит все горные породы, содержащие глинистые минералы, каменные и калийные соли, а также все другие породы в массиве, имеющие кососекущую трещиноватость [2].
Кафедрой геологии МГГУ изучались деформируемость и прочность (сопротивление сдвигу) горных пород во времени для глинистых пород бортовых массивов карьеров КМ А, техногенных отложений различных горнопромышленных регионов и пород потенциально оползневых склонов по трассе комбинированной (автомобильной и железной) дороги до олимпийских объектов горного кластера близ Красной Поляны.
Согласно современным представлениям деформацию глинистой породы целесообразно рассматривать как следствие изменения структуры, т.е. одновременного развития дефектов и переориентации минеральных частиц, что подтверждается экспериментальными исследованиями. Таким образом, индикатором напряженно-деформированного состояния можно считать структуру породы и ее перестройку. С учетом необходимости определения
критических (разрушающих) деформаций пород бортовых массивов целесообразно использование закономерности деформирования, описывающие все три стадии процесса ползучести.
Физические основы структурно-механического подхода сформулированы С.С. Вяловым и сводятся к следующим положениям [5]:
1. Пластическое деформирование породы вызывается смещением частиц и их переориентацией. Длительное разрушение связано с развитием дефектов и микротрещинообразованием.
2. Разрушение наступает, когда плотность дефектов достигает определенной величины; плотность ориентированных частиц пропорциональна длительности процесса т>тк.
При математической формулировке этих положений процесс образования дефектов рассматривался как случайный (стохастический) по схеме мгновенных повреждений, причем изменение (активация) необходимой для смещения частиц энергии предполагается пропорциональным вероятности повреждения. Согласно кинетической теории прочности и ползучести, деформирование и разрушение горных пород связываются с преодолением минеральными частицами и их агрегатами энергетического барьера и переходом в новое положение равновесия под воздействием активирующей эти частицы силы.
Зависимость для определения скорости ползучести при использовании предложенных Ю.К. Зарецким реологических уравнений состояния принимает следующий вид [5]:
и = -
г_
По
ехр
к
(1)
где п0 = сопб1 = ехр
кв
началь-
юшем подборе параметров переходит в известные эмпирические уравнения. Поэтому это уравнение может рассматриваться в качестве закономерности деформирования как для глинистых, так и для твердых пород.
С учетом условия достижения критической плотности микротре-шин юр = const в момент разрушения горной породы и предложенного Ю.К. Зарецким и С.С. Вяловым выражения для описания процесса треши-нообразования уравнение длительной прочности принимаем в виде [5]:
Т = Т0 - (то ) T Р,
T + p (3)
где ip - момент разрушения породы, когда плотность трешин достигает критического значения; T - параметр, имеюший размерность времени.
Для определения параметров ползучести горных пород (пр, Пт, или скоростей деформации следует использовать равенства, вытекаюшие из уравнений (1) и (3):
Таблица 1
Характеристики длительной прочности и сдвиговой ползучести пород бортового массива
возраст пород ЛИТОЛО- гическая колонка мощность СЛОЯ, м краткое описание пород 7 i тнЛг W % д тн/м3 кг/см2 Сг Ф кг/см2 град тР лет Ч sm 6
Q Ж 15 суглинок лессовидный 1.93 18 2.68 0.32 0.20 29 0.142 0.267 0.85
Kf llllll 6 песок разнозернистый 2.00 22 2.65 35
^tpl'W — — — — 35 глнна песчаная, тугоп ластичн ая 2.11 20 2.69 0.43 0.10 20 0.369 0.284 0.94
jvlg-cl — 28 глина тугоштастичная 2.08 16 2.72 0.75 0.24 16 0.207 0.162 1.00
5 песок мелкозернистый 1.90 24 2.76 35 -
I-—~ 10 глина плотная, тугопластичная 2.00 26 2.71 1,30 0.80 22 0.186 0.115 1.00
5 песок разнозернистый - 1.90 - 24 2.76 35
D3 50 глина песчаная, тугоп ластич ная 2.21 18 2.77 3.00 1.70 12 0.412 0.062 1.00
ная вязкость, зависящая от начальной структуры, плотности, влажности и температуры 9; и0* - начальная энергия активации; к - постоянная Боль-цмана; 5 - безразмерный параметр; ТТ - параметр, имеющий размерность
времени; т-т« = т - уровень напря-
т0 -т
жений; т, тк, т0 - соответственно текущее, длительное и условно-мгновенное значения сопротивления сдвигу; t - время начала деформирования.
Из формулы (1) получаем зависимость
для расчета деформаций сдвига в виде:
t
и^) = и0 + а | ехр( xt ^ (2)
о
где и0 - начальная деформация; 5т
х = т~ ; а = т/п0 1Т
Универсальность выражения (2) подтверждается тем, что оно обобщается на случай сложного напряженного состояния, а также при соответству-
ит
Пт ;
-8
а2 _ т
(4)
где ^ - момент начала пластического течения; ^ - момент разрушения; а1, а2, 8 - параметры сдвиговой ползучести; ит - скорость установившейся ползучести; 'ир - скорость разрушения.
Структурно-механические уравнения ползучести и длительной прочности могут также использоваться для описания процессов сдвиговой ползучести в отвальных массивах и их основаниях, сложенных глинистыми, мерзлыми или заснеженными породами.
Для определения параметров длительной прочности (сдвиговой ползучести) выполнялись обратные оползневые расчеты по материалам выполненной институтом ВИОГЕМ
паспортизации оползней карьера Михайловского ГОКа КМА [5], а также (для девонских глин) длительных трехосных испытаний в стабилометрах УСВ-2.
В табл. 1 дана инженерно-геологическая колонка вскрышной толщи Восточного борта Михайловского карьера КМА, на рис. 1 графики зависимостей критических (разрушающих) скоростей деформаций от сдвигающих (т) и нормальных напряжений. Эти графики являются основой для контроля по деформациям состояния бортовых массивов карьеров КМА. Для отвальных насыпей, дамб гидроотвалов и хво-стохранилищ эффективный контроль за состоянием откосных сооружений осуществляется путем комплексного зондирования приоткосных зон и использованием стационарных датчиков-пьезодинамометров, заложенных по
р
Рис. 1. Критические скорости деформаций глинистых пород бортовых карьеров КМА: а - суглинки четвертичные; б - глины апт-неокомские; в - глины верхнеюрские (волж-ско-келловейские); г - глины среднеюрские (батские); д - глины девонские; параметры а , МПа: 1 - 0,05; 2 - 0,1; 3 - 0,15; 4 - 0,2; 5 - 0,25; 6 - 0,3; 7 - 0,35; 8 - 0,4 п
Рис. 2. Зависимости коэффициента запаса устойчивости от степени обводненности потенциально оползневых массивов для трех расчетных профилей в районе железнодорожного вокзала Альпика-Сервис
расчетным профилям в теле и основании дамбы на различных этапах формирования намывного массива.
В МГГУ разработана программа оперативного определения коэффициента запаса устойчивости п в зависимости от измеренного пьезо-динамометрами давления воды, приведенного к вероятным поверхностям скольжения. Расчет устойчивости выполняется методами алгебраического суммирования и многоугольника сил. Определение текущего коэффициента запаса устойчивости производится в зависимости от площади эпюры давления воды, определяемой путем снятия замеров величин Р или И по
1 ш ш
вероятной поверхности скольжения.
Комплекс исследований МГГУ, выполненных применительно к откосным сооружениям в окрестности вокзала «Альпика - Сервис» имел главной целью установление прочностных характеристик (сопротивления сдвигу) глинистых сланцев, подстилающих толщу раздельно-обломочных отложений при различном положении де-прессионной кривой. Установлено, что физико-механические параметры глинистых сланцев в среднем составляют: плотность у = 2,6 г/см3, удель-
ное сцепление С = 4,95 т/м2, угол внутреннего трения ф = 21°.
Проведены расчеты устойчивости склонов (по согласованным профилям) при различном обводнении при-откосных зон. На рис. 2 представлены зависимости коэффициента запаса устойчивости склонов от степени их обводнения, а также были установлены величины некомпенсированного оползневого давления.
В работе В.И. Осипова [8] отмечается, что при определении эффективных напряжений в тонкодисперсных глинистых образованиях с помощью теории К. Терцаги возникают расхождения экспериментальных и расчетных данных из-за неучета распределения напряжений на контактах структурных элементов. В частности в теории К. Терцаги не рассматриваются:
1. особенности распределения напряжений на контактах различных геометрических и энергетических типов;
2. существование на контактах тонких пленок связанной воды, обладающих расклинивающим давлением;
3. развитие на контактах различных процессов физико-химической природы, обусловливающих существование внутренних напряжений;
4. зависимость порового давления от характера порового пространства глин.
Для решения инженерных задач представляется целесообразным использовать упрошенный подход, предусматривающий испытания водонасы-шенных тонкодисперсных отложений в условиях трехосного сжатия (например, в универсальных стабилометрах компрессионного типа УСВ-2 конструкции ВИОГЕМ, автор - Н.П. Верещагин). В условиях закрытой системы (неконсолидированно - недре-нированные испытания) к образцу прикладывается внешняя нагрузка и определяется коэффициент порового давления, т.е. доля нагрузки, воспринимаемая поровой водой. В дальнейшем расчеты уплотнения слабоструктурных водонасыщенных глинистых грунтов производятся с использованием соотношений для расчета избыточного порового давления, учитывающих также начальный градиент фильтрации, ползучесть минерального скелета, темпы нагружения и др. [4].
Определение реального общего эффективного напряжения для водонасыщенных слабоструктурных отложений (илов, слабо литифици-рованных глинистых образований, намывных глинистых грунтов) требует установления исходных параметров к выражению [8]:
<Ф = X
а- P
-n(h) • а, 1 =
= (а- Ри) -хП( h) • ас
(5)
где х - число контактов в единице площади горизонтальной площадки, П(И) - расклинивающее давление ги-дратных пленок, ас - площадь контакта, а - общее внешнее напряжение, Ри - поровое давление.
Установление параметров к формуле (5) требует проведения сложных дорогостоящих исследований.
Решение задач фильтрационной консолидации позволяют удовлетворительно описывать развитие осадок тонкодисперсных намывных пород -в пределах 80-90% максимальной (стабилизированной) осадки. Для повышения точности описания процесса уплотнения может быть использовано решение одномерной задачи консолидации и ползучести двухфазных пород с учетом взаимодействия фаз, полученное Ю.К. Зарецким [4] .
С учетом требований безаварийной эксплуатации комбинированной дороги до горного кластера следует выполнить дополнительный комплекс работ по определению параметров длительной прочности и сдвиговой ползучести (см. выше) глинистых пород в левобережном склоне р. Мзым-та. В работе В.И. Осипова, Ю.А. Мамаева, A.A. Ястребова выполнено ранжирование по степени оползневой опасности участков строительства олимпийских объектов в г. Сочи [9]. Вполне правомерны рекомендации авторов этой статьи о создании комплексной системы геодинамического мониторинга, включающей также измерения сдвиговых склоновых смещений и уровней подземных вод. При этом приоритетными, по нашему мнению, являются бескабельные схемы контроля, что повысит жизнестойкость измерительной системы.
Для контроля за состоянием труднодоступных зон гидроотвалов МГГУ совместно с ВИОГЕМ предложен и запатентован способ дистанционного определения несущей способности намывного массива по данным аэрофотограмметрической съемки, в соответствии с которым допустимая нагрузка Рдоп определяется в зависимости от относительной осадки массива Ah/h.
Инструментальные наблюдения выполнялись институтом ВИОГЕМ с использованием метода аналитиче-
ской фототриангуляции, позволяющего определять координаты точек недоступных участков гидроотвала бесконтактным способом. Аэрофотосъемки производились в моменты начала и окончания перерыва tо в намыве («отдыха» гидроотвала).
Эталонные зависимости
Рдоп = /[ДИ(^)/И] устанавливаются с использованием решения задачи уплотнения намывного слоя в период его «отдыха» tо, из которого определяется степень уплотнения слоя и = Через степень уплотнения по эмпирическим зависимостям определяются характеристики сопротиивления сдвигу с и ф. С использованием формулы Прандтля-Рейснера по значениям с(и) и ф(и) определяется предельная критическая нагрузка на намывное основание Р кр, а затем Р = Р кр /п [4].
пр ' доп пр II-"
Предложенный способ позволяет расширить область применения аэрофотосъемки благодаря возможности оперативного определения несущей способности труднодоступных зон намывных массивов и обеспечивает проведение дистанционного контроля состояния гидроотвала, повышая тем самым безопасность ведения работ при одновременном снижении их стоимости и трудоемкости.
Методика дистанционного контроля уплотнения и, соответственно, несущей способности намывных массивов апробирована в условиях гидроотвала «Березовый Лог» (КМА). Результаты аэрофотосъемки показали достаточно хорошую сходимость с данными наземных съемок и расчетами.
На сегодняшний целесообразно рассматривать варианты применения
GPS технологий и космической снимков при дистанционном контроле уплотнения техногенных массивов. Основное преимущество проведения съемки ситуации с применением спутниковой технологии заключается в том, что при ее осуществлении отпадает необходимость создания геодезических сетей сгущения, создания съемочного обоснования и его сгущения, поскольку методы спутниковых определений по дальности и точности принципиально обеспечивают возможность проведения съемочных работ непосредственно на основе государственной геодезической и нивелирной сети.
Применение съемки GPS обусловлено уменьшением трудозатрат, повышением точности результатов измерений. В сравнении с аэрофотосъемкой данный метод является более экономичным и точным.
Аэрофотосъемка и космическая съемка намывных территорий обеспечивает наибольший эффект при сочетании ее с зондированием гидроотвалов комбинированными зондами МГГУ и пенетрационно-каротажными станциями.
Можно рекомендовать также метод контроля устойчивости дамб по деформациям, основой которого контроля является регистрация скоростей смещения грунта в отдельных точках (реперах), размещаемых по соответствующим профилям и сопоставление их с расчетными значениями. Для количественной оценки сдвиговых деформаций в теле дамбы и ее основании могут быть использованы зависимости структурной механики грунтов.
1. Панюков П.Н. Инженерная геология. - М.: Недра, 1978. - 206 с.
2. Фисенко Г.Л. Устойчивость бортов карьеров и отвалов. - М.: Недра, 1965. - 378 с.
- СПИСОК ЛИТЕРАТУРЫ
3. Гальперин A.M., Шафаренко Е.М. Реологические расчеты горнотехнических сооружений. - М.: Недра, 1977. -246 с.
4. Гальперин A.M. Геомеханика открытых горных работ. - М.: Горная книга, 2012. - 480 с.
5. Зареикий Ю.К., Вялов С.С. Вопросы структурной механики глинистых грунтов // Основания, фундаменты, механика грунтов. - 1971. - № 3. - С. 1-5.
6. Galperin A.M., Panfilov A.Y. Back calculations for the determination in open-pit slopes and hydrofilled structures. Proc. Of the Int. Workshop in Geoenvironment and Geotechnics (GEOENV 2005). Milos, Greece, 2005. P. 131-136.
КОРОТКО ОБ АВТОРАХ
7. Гальперин A.M., Зайцев B.C., Марченко С.М. и др. Способ контроля состояния намывных массивов. Патент РФ № 1188322 15.07.1993.
8. Осипов В.И. Физико-химическая теория эффективных напряжений в грунтах. -М.: ИФЗ РАН. - 74 с.
9. Осипов В.И., Мамаев Ю.А., Ястребов A.A. Оценка ранжирования по степени оползневой опасности участков строительства олимпийских объектов в Сочи // Геоэкология. - 2013. - № 6. - 530-537 с. ЕШЗ
Гальперин Анатолий Моисеевич - доктор технических наук, профессор, зав. кафедрой, Ческидов Василий Владимирович - кандидат технических наук, доцент, Бородина Юлия Владимировна - аспирант, Демидов Артем Владиславович - аспирант,
Московский государственный горный университет, e-mail: ud@msmu.ru.
UDC 504.064.2
STUDY NONSTATIONARY HYDROGEOMECHANICAL PROCESSES IN CLAY ROCKS
Galperin A.M., Doctor of Technical Sciences, Professor,
Cheskidov V. V., Candidate of Engineering Sciences, Assistant Professor,
Borodina J.V., Graduate Student,
Demidov A.V., Graduate Student,
Moscow State Mining University, e-mail: ud@msmu.ru.
The questions of deformation , strength and creeping of 'clay deposits assessment are lighted on the base of the objects KMA and the regions of the Olympic objects construction not far from the town Sochi. The possibilities of using the remote control of the geomechanical processes under the introduction of mining and constructing works, namely using aerial mapping and cosmo-mapping and GPS technologies are considered.
Key words: mining works, hydrogeomechanical processes, massifs of rocks, cohesive soils (clay), creeping, remote control.
REFERENCES
1. Panjukov P.N. Inzhenernaja geologija (Engineering geology), Moscow, Nedra, 1978, 206 p.
2. Fisenko G.L. Ustojchivost' bortov karerov i otvalov (Stability of pitwalls and dumps), Moscow, Nedra, 1965, 378 p.
3. Gal'perin A.M., Shafarenko E.M. Reologicheskie raschety gornotehnicheskih sooruzhenij (Rheological calculations of mine technical structures), Moscow, Nedra, 1977, 246 p.
4. Gal'perin A.M. Geomehanika otkrytyh gornyh rabot (Geomechanics in open pit mining), Moscow, Gornaja kniga, 2012, 480 p.
5. Zareckij Ju.K., Vjalov S.S. Osnovanija, fundamenty, mehanika gruntov, 1971, no 3, pp. 1-5.
6. Galperin A.M., Panfilov A.Y. Back calculations for the determination in open-pit slopes and hydrofilled structures. Proc. Of the Int. Workshop in Geoenvironment and Geotechnics (GEOENV 2005). Milos, Greece, 2005, pp. 131-136.
7. Gal'perin A.M., Zajcev V.S., Marchenko S.M. Sposob kontrolja sostojanija namyvnyh massivov (Aggra-dational deposit control procedure). Patent RU № 1188322 15.07.1993.
8. Osipov V.I. Fiziko-himicheskaja teorija jeffektivnyh naprjazhenij v gruntah (Physicochemical theory of actual stresses in soil), Moscow, IFZ RAN, 74 p.
9. Osipov V.I., Mamaev Ju.A., Jastrebov A.A. Geojekologija, 2013, no 6, 530-537 p.
