Полный текст

Введение Согласно инструкции ВСН 189-78, действующей на территории Российской Федерации, сложные инженерно-геологические и гидрогеологические условия строительства подземных сооружений требуют применения специальных способов производства работ, в частности искусственного замораживания грунтов. Применительно к строительству шахтных стволов искусственное замораживание грунтов и горных пород используется для создания ледопородного ограждения (ледопородной завесы, ледопородного цилиндра) вокруг запроектированной горной выработки (рис. 1). Для этого по контуру запроектированной к проходке выработки бурятся скважины, в которые опускаются замораживающие колонки. За счет работы холодильных станций по колонкам осуществляется циркуляция хладоносителя (рассола) с поддерживаемой температурой до -50 °С [1]. В результате циркуляции охлаждающего рассола в замораживающей колонке окружающий ее породный массив постепенно охлаждается, а содержащаяся в массиве вода кристаллизуется. По истечении некоторого времени вокруг замораживающих колонок образуются одиночные ледопородные цилиндры, которые в дальнейшем смыкаются, образуя сплошное ледопородное ограждение, предотвращающее попадание грунтовых вод в строящуюся горную выработку [2]. Расчеты ледопородного ограждения сводятся к статическим, определяющим необходимую толщину ледопородного ограждения из условий прочности, устойчивости или по предельно допустимым деформациям, и к теплотехническим, на основании которых определяются холодопроизводительность замораживающей установки, режим и время замораживания, температурное поле слоев грунтов и горных пород [3]. Вопросы, относящиеся к теплотехническим расчетам, освещены в классических трудах советских ученых [4, 5] и более новых трудах российских [6-8] и зарубежных [9-11] исследователей. В данной работе рассмотрены вопросы статических расчетов ледопородного ограждения. Обзор статических расчетов ледопородного ограждения Статический расчет защитных завес из замороженных грунтов должен производиться по двум предельным состояниям - по предельному напряженному состоянию (расчет на прочность) и предельным деформациям (расчет на ползучесть) [3]. На сегодняшний день в литературе представлены в основном методы расчета по предельному напряженному состоянию [3, 12-15]. Они сводятся к определению такой толщины ледопородного цилиндра, при которой действующие в заданный момент времени напряжения не превосходили бы предельного сопротивления сдвигу замороженных пород в этот же момент времени. Как правило, рассматривается стандартный набор гипотез, позволяющих получить аналитическое выражение для толщины ледопородного ограждения: 1. Замороженные породы цилиндра принимаются изотропным материалом с однородным распределением температуры и однородно распределенными прочностными характеристиками. 2. Объемные деформирования замороженной породы предельны. 3. При расчете по предельным деформациям влияние фактора времени на механические свойства замороженных пород учитывается в параметрической форме. В литературе, как правило, рассматривается случай ледопородного ограждения бесконечной протяженности в вертикальном направлении, что позволяет рассматривать двухмерную задачу в рамках гипотезы о плоском деформированном состоянии (рис. 1). Рис. 1. Горизонтальный разрез ледопородного ограждения (расчетная область): 1 - замороженный грунт (ледопородное ограждение); 2 - талый грунт; 3 - внутренний объем горной выработки Fig. 1. Horizontal section of the frozen wall: 1 - frozen wall; 2 - unfrozen soil; 3 - interior of the excavation На внешней границе ледопородного ограждения задается равномерно распределенная нагрузка P - боковое давление грунта, зависящее от глубины залегания и гидростатического давления грунтовых вод. На внутренней стенке ледопородного ограждения задается свободная поверхность с нулевым давлением. В действующей инструкции ВСН 189-78 для определения толщины E стенки ледопородного ограждения применяется известная формула сопротивления материалов для расчета толстостенных труб (получаемая из формулы Ляме) (1) или эмпирическая формула Домке (2) где - внутренний радиус ледопородного цилиндра, м; - максимальное давление на внешнюю поверхность цилиндра, Па; - предел прочности замороженного грунта на сжатие, Па; - допускаемое напряжение на сжатие для замороженного грунта, Па, В целом расчеты по приведенным формулам дают излишний запас прочности и с большим приближением могут быть применены для неглубоких выработок, когда горное давление невелико [3]. В монографии [3] предложена формула для расчета толщины ледопородного ограждения неограниченной высоты для предельного напряженного состояния замороженного грунта, описываемого линейным законом Кулона - Мора (3) где , (4) . (5) где , - коэффициенты в законе Кулона - Мора, записанном в главных напряжениях; - сцепление замороженного грунта, Па; - угол внутреннего трения замороженного грунта, рад. В работе [15], опубликованной позднее, получена формула, совпадающая с (3). В работе [12] формула (3) уточнена на предмет учета упругопластических деформаций талого грунта, окружающего ледопородное ограждение: (6) где , (7) , (8) где , - коэффициенты в законе Кулона - Мора, записанного в главных напряжениях для талого грунта. В работе [13] формула (6) уточнена с учетом изменения толщины ледопородного ограждения вследствие его деформирования после проходки горной выработки. Рассчитана толщина ледопородного ограждения до его деформирования: , (9) где , (10) где - модуль сдвига талого грунта, Па; - толщина ледопородного ограждения, рассчитанная из (6), м. Методы расчета толщины стенки ледопородного ограждения по предельным деформациям в литературе представлены слабо. Единственным известным и широко применяемым на практике является метод, предложенный в монографии [3]. В данном методе в качестве критерия разрушения выступает предельное радиальное перемещение внутренней поверхности ледопородного ограждения , а сам критерий формулируется следующим образом: . (11) С механической точки соотношение (11) определяет такой размер ледопородного ограждения, при котором развивающаяся за заданный период времени деформация ползучести не приводила бы к запредельному радиальному перемещению внутренней поверхности ледопородного ограждения. , (12) где - напряжение, Па; - деформация; - коэффициент пропорциональности, Па. Толщина E стенки ледопородного ограждения конечной протяженности [3] , (13) где - коэффициент упрочнения замороженного грунта; - высота ледопородного цилиндра, м; - параметр, характеризующий защемление торцов цилиндра. При выводе формулы (13) использовалась гипотеза о том, что касательное напряжение вблизи торцов ледопородного цилиндра ограниченной толщины зависит от вертикальной координаты Z линейно. В [3] также получена толщина E стенки ледопородного ограждения неограниченной высоты: . (14) При выводе (14) принималось, что касательное напряжение во всем объеме цилиндра равно нулю. В ГОСТ 12248-2010 «Грунты. Методы лабораторного определения характеристик прочности и деформируемости» представлена методика определения нелинейных характеристик деформируемости мерзлого грунта по результатам испытания методом одноосного сжатия. Временная зависимость коэффициента нелинейной деформации замороженного грунта при этом имеет степенной вид: , (15) где - коэффициент нелинейной деформации замороженного грунта, рассчитанный для времени нагружения , равного 1 ч, Па. Сравнительный анализ статических расчетов ледопородного ограждения Для сравнения различных формул расчета толщины ледопородного ограждения рассмотрим слой глины с механическими свойствами, представленными в табл. 1. Использование формулы для расчета по предельным деформациям наряду с дополнительными материальными константами, описывающими процесс ползучести материала, требует задания условий строительства ледопородного ограждения (таких как допустимая высота неподкрепленного участка ледопородного ограждения, предельные деформации), условий взаимодействия незакрепленного участка ледопородного ограждения с подстилающим грунтом и конструкциями шахтного ствола. В силу этого является логичным отдельно провести сравнение формул, полученных из условия предельного напряженного состояния, и затем дополнительно сравнить результаты расчета с оценкой по предельным деформациям в конкретном технологическом процессе. На рис. 2 представлены зависимости проектных толщин ледопородного ограждения для рассматриваемого слоя глины от глубины его залегания. Проектные толщины рассчитывались по формулам (1)-(3), (6), (9). Из рис. 2 видно, что формула Ламе (1) дает наиболее завышенные значения и применима в данном случае только для глубин до 100 м. Формула Домке (2) дает хорошие результаты для глубин до 200 м, а после 300 м дает существенно завышенные результаты. Формула Вялова для условия предельного напряженного состояния (3) дает более высокие значения толщины ледопородного ограждения, чем формула Домке для глубин до 170 м, в то время как при больших глубинах имеет более пологий вид и позволяет получить более адекватные значения для толщины ледопородного ограждения. Формула Янга (6), учитывающая деформируемость талого грунта, позволяет уменьшить проектную толщину ледопородного ограждения, рассчитанную по формуле Вялова, в среднем на 30-40 %. Формула Чжана (9), дополнительно учитывающая изменение толщины ледопородного ограждения вследствие его деформирования после проходки горной выработки, позволяет снизить проектную толщину ледопородного ограждения еще на 3-5 %. Таблица 1 Прочностные свойства модельного слоя глины Table 1 Strength properties of the model clay layer Свойства Значения Сцепление замороженного грунта, МПа 4,0 Угол внутреннего трения замороженного грунта, рад 8 Сцепление талого грунта, МПа 0,3 Угол внутреннего трения талого грунта, рад 15 Модуль сдвига талого грунта, МПа 38,5 Предел прочности замороженного грунта на сжатие, Па 8 Радиус ствола, м 5 Повышение давления бокового распора с увеличением глубины на 1 м, МПа 0,013 Рис. 2. Зависимость различных формул для расчета толщины ледопородного ограждения по предельному напряженному состоянию Fig. 2. Various formulas for calculating the thickness of the frozen wall for the boundary state of stress Сравнительный анализ формул (13) и (14) по предельным деформациям и формул (1) и (2) по предельным напряжениям проводился для условий строящегося рудника ОАО «Беларуськалий». На рис. 3 представлены расчетные толщины ледопородного ограждения для 7 обводненных слоев горных пород на глубинах от 0 до 275 м. Формулы (3), (6) и (9) в данном случае не анализировались по причине отсутствия надежных экспериментальных значений для сцепления и угла внутреннего трения замороженных пород. Рис. 3. Расчет проектной толщины ледопородного ограждения для строящихся стволов рудника ОАО «Беларуськалий» Fig. 3. Calculation of the design thickness of the frozen wall for the constructed shafts of the "Belaruskali" Company Использованные при расчете прочностные свойства обводненных слоев горных пород представлены в табл. 2. Таблица 2 Прочностные свойства обводненных слоев пород Table. 2 Strength properties of water-saturated rock layers Слой Глубина кровли, м Длительная прочность на сжатие (24 ч), МПа Коэффициент нелинейной деформации, МПа Коэффициент упрочнения Глина 23,2 1,4 2,2 0,37 Алевриты 58,0 2,3 3,9 0,33 Песок 85,0 4,8 5,5 0,29 Мел 136,0 6,5 6,9 0,18 Песчаник 146,0 7,3 7,7 0,27 Глинистый песок 255,8 5,1 8,2 0,25 Аргиллитовая глина 275,0 8,6 9,8 0,14 Допустимое перемещение внутренней стенки ледопородного ограждения принималось равным 10 см. Внутренний радиус ледопородного ограждения равен 5,25 м. Как видно из рис. 3, проектная толщина, рассчитанная по формуле Ламе, как и ранее, дает существенно завышенные значения для слоев горных пород начиная со 2-го. Проектная толщина, рассчитанная по формулам (13) и (14) по предельным деформациям, выше толщины (2), рассчитанной по предельным напряжениям, на 25-220 %. Наибольшие расхождения наблюдаются для слоев горных пород с более выраженной ползучестью - глина и влагонасыщенный песок. При этом формулы расчета проектной толщины ледопородного ограждения по предельному напряженному состоянию не учитывают реологические свойства горных пород, которые заложены в уравнении состояния. Заключение В работе проведен анализ существующих методов расчета проектной толщины ледопородного ограждения для шахтных стволов. Проведено количественное сравнение формул для расчета толщины ледопородного ограждения по критерию предельного напряжения. Показано, что: 1. Формула Ламе в рассматриваемом случае дает завышенные значения толщин ледопородного ограждения и применима только для глубин до 100 м. 2. Формула Домке применима при глубинах до 200 м. 3. Формула Вялова по предельному напряженному состоянию имеет более пологий вид для глубин более 200 м по сравнению с формулами Домке и Ламе и применима для больших глубин. 4. Учет дополнительных факторов в формулах Янга и Чжана (действительных деформационных характеристик талого грунта, деформированного состояния системы до проходки ствола) позволяет снизить проектную толщину ледопородного ограждения еще на 32-43 %. 5. В случае ярко выраженной ползучести породы расчетная толщина ледопородного ограждения по предельным деформациям может существенно превышать расчетную толщину по предельным напряжениям. Это подтверждает вывод, сформулированный ранее в [4], что статический расчет защитных завес из замороженных грунтов должен производиться по двум предельным состояниям - предельному напряженному состоянию и предельным деформациям, а в процессе проектирования принимается во внимание наибольшее полученное значение. Расчет по предельным деформациям может быть произведен по формуле Вялова (12).

Об авторах

Л. Ю Левин

Горный институт Уральского отделения РАН

М. А Семин

Горный институт Уральского отделения РАН

О. А Плехов

Институт механики сплошных сред Уральского отделения РАН

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

Статистика

Просмотры

Аннотация - 223

Ссылки

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