Санкт-Петербургский государственный университеТНа правах рукописиЗАЦЕПИН Михаил АлександровичМатематическое моделированиенапряженно-деформированного состояниямассива горных породпри разработке пологих месторождений01.02.04 – Механика деформируемого твердого телаА В Т О Р Е Ф Е Р А Т диссертации на соискание ученой степени кандидата физико-математических наукСанкт-Петербург 2009 Работа выполнена в государственном образовательном учреждении высшего профессионального образования Санкт-Петербургском государственном горном институте (техническом университете).Научный руководитель: доктор технических наук, профессор Господариков Александр ПетровичОфициальные оппоненты: доктор физико-математических наук, профессор Филиппов Сергей Борисович (Санкт-Петербургский государственный университет)доктор физико-математических наук, профессор Куксенко Виктор Степанович (Физико-технический институт им. А.Ф. Иоффе)Ведущая организация: Санкт-Петербургский государственный политехнический университетЗащита диссертации состоится "02" апреля 2009 г. в 15 часов на заседании совета Д 212.232.30 по защите докторских и кандидатских диссертаций при Санкт-Петербургском государственном университете по адресу: 198504, Санкт-Петербург, Петродворец, Университетский пр., 28, математико-механический факультет, ауд. 405.С диссертацией можно ознакомиться в Научной библиотеке им. М. Горького Санкт-Петербургского государственного университета по адресу: 199034, Санкт-Петербург, Университетская наб., 7/9.Автореферат разослан " " _______ 2009 г.Ученый секретарь диссертационного совета, доктор физико-математических наук, профессор С.А.Зегжда^ ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫАктуальность темы. Многообразие горно-геологических условий залегания пологих пластов и продолжающийся рост глубин разработки месторождений полезных ископаемых приводят исследователя к необходимости анализа напряженно-деформированного состояния (НДС) массивов горных пород вокруг подземных горных выработок всевозможного назначения и различного очертания. Разработка месторождений полезных ископаемых вызывает в массиве горных пород целый комплекс различных механических процессов: перераспределение деформаций и напряжений, прорывы подземных вод и плывунов, динамическое разрушение пород в форме вывалов, внезапных выбросов, горных ударов и многие другие. Изучение указанных явлений, протекающих в массиве горных пород при отработке различных месторождений, является ответственной задачей, так как они предопределяют не только экономическую часть разработки и безопасность работ, но во многих случаях и практическую целесообразность эксплуатации месторождения или какой-нибудь его части. Поэтому решение проблемы по определению прочности и надежности подземных сооружений при воздействии на них статических и динамических нагрузок имеет важное экономическое и социальное значение.^ Целью диссертационной работы является повышение безопасности проходки и поддержания подземных горных выработок в условиях неоднородного слоистого массива горных пород при разработке пологих месторождений полезных ископаемых.^ Основные методы исследований. При решении поставленных задач применялись как аналитические методы, например, вариационный метод в форме В.З. Власова, так и численные методы. На основе метода конечных элементов разработана вычислительная программа. Ряд численных результатов получен с помощью программного продукта конечно-разностного анализа FLAC v.4.0.^ Основные результаты, выносимые на защиту: - решена задача об определении НДС массива горных пород, вмещающего подземную горную выработку кругового сечения, при этом горные породы рассматривались как однородная и изотропная среда; - при помощи программного продукта FLAC v.4.0 построена модель слоистого массива горных пород, вмещающего подземную горную выработку прямоугольного сечения в начальный момент работы лавы, и определено НДС данного массива; - разработан алгоритм численно-аналитического моделирования напряженно-деформированного состояния неоднородной слоистой среды при решении важных прикладных задач геомеханики.^ Научная новизна диссертационной работы: - установлены закономерности изменения напряженного состояния неоднородного породного массива, вмещающего выработки, в рамках физически линейного процесса деформирования горных пород; - сформулирован и разработан численно-аналитический метод моделирования процесса деформации массива горных пород, вмещающего подземные горные выработки различного назначения и очертания.Достоверность научных положений и выводов диссертационной работы достигается: - строгостью математической постановки исследуемых задач и широким сопоставлением авторских решений важных прикладных задач геомеханики с решениями подобных задач другими исследователями, полученными на основе различных методов; - получением решений с высокой степенью точности за счет внутренней сходимости результатов; - сходимостью результатов, полученных на основе применения вычислительной программы FLAC v.4.0 и созданного программного продукта, с результатами натурных (шахтных) исследований, проведенных на Воркутском и Старобинском месторождениях.^ Практическое значение результатов исследований: результаты исследований по оценке НДС неоднородного слоистого породного массива с вмещающими выработками могут быть использованы на шахтах ОАО «Воркутауголь» и на рудниках РУП ПО «Беларуськалий» при разработке схем и способов управления горным давлением, в том числе и на глубоких горизонтах.^ Апробация результатов работы. Основные положения диссертационной работы были доложены на III и V научной конференции «Нелинейные проблемы механики и физики деформируемого твердого тела» (СПбГУ, 2004 г., 2006 г.), на ежегодном симпозиуме факультетов геотехники и горного дела Фрайбергской горной академии (Германия, 2005 г.), на ежегодной сессии студенческих докладов, посвященных горному делу, в Краковской горной академии (Польша, 2005 г.), на всероссийской научной конференции-конкурсе студентов выпускного курса (СПГГИ, 2006 г.), на международном форуме молодых ученых (СПб, 2006 г.), на международной научной конференции по механике (СПбГУ, 2006 г.), на ежегодной научной конференции молодых ученых «Полезные ископаемые России и их освоение» (СПГГИ, 2003 г., 2006 г., 2007 г.), на всероссийской молодежной научно-практической конференции по проблемам недропользования (Екатеринбург, 2008 г.), на 6-й межрегиональной научно-практической конференции (Воркута, 2008 г.). Диссертация в целом была доложена на научных семинарах кафедры вычислительных методов механики деформируемого твердого тела Санкт-Петербургского государственного университета и кафедры разработки месторождений полезных ископаемых Санкт-Петербургского государственного горного института.^ Публикация результатов. По теме диссертации опубликовано 6 печатных работ в вузовских и межвузовских сборниках научных трудов, из них статья [6] опубликована в журнале, рекомендованном ВАК. В совместной публикации [1] диссертанту принадлежит анализ динамических обрушений основной кровли применительно к Старобинскому месторождению калийных солей, а также предложены мероприятия по различному размещению закладки в выработанном пространстве. В публикации [2] Зацепину М.А. принадлежит определение закона изменения опорного давления в предположении, что зона разрушения пласта перед забоем отсутствует. В публикации [3] автору диссертации принадлежит разработка вычислительной программы, реализующей численный алгоритм метода конечных элементов, на основе которой определяются параметры управления кровлей и напряжения в краевой части пласта. В статье [4] Зацепину М.А. принадлежит идея учета явления зональной дезинтеграции горных пород, которая может быть использована в схеме «целик-пролет-целик» для определения НДС массива горных пород.^ Объем и структура работы. Диссертационная работа состоит из введения, трех глав, заключения, списка литературы, включающего 82 наименования, и изложена на 101 странице машинописного текста, содержит 24 рисунка и 4 таблицы.^ ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫВо введении обосновывается актуальность темы диссертации, формулируются основная идея, научная новизна и практическая значимость диссертационной работы, определяются цели и задачи исследования, указываются методы исследований, перечисляются основные научные положения, выносимые на защиту.^ В первой главе сделан краткий обзор основных исследований напряженно-деформированного состояния слоистого пологозалегающего породного массива. Рассмотрены основные математические модели, описывающие физическое состояние массива горных пород, и приводится обзор основных численных методов определения параметров полей напряжений, перемещений и деформаций для пологого слоистого массива горных пород. В последние годы наметилась тенденция увеличения глубины отработки пластовых месторождений (свыше 1000 м), которая характеризуется необходимостью применения дополнительных мер по охране выработок. Для предотвращения аварийных ситуаций при эксплуатации подземных горных выработок необходим правильный выбор способов охраны, соответствующий ожидаемым смещениям пород, вмещающих данную выработку. В различных горных выработках в зависимости от положения очистного забоя и способа управления горным давлением выделяют несколько характерных периодов их поддержания, в каждом из которых горное давление проявляется с различной интенсивностью. Практика ведения горных работ показывает, что очистные работы являются одним из главных факторов, определяющих интенсивность деформации породного контура горных выработок. Различными научно-исследовательскими институтами и Санкт-Петербургским горным институтом постоянно ведутся исследования по анализу состояния выработок и разработке предложений по улучшению их поддержания и охраны в условиях угольных шахт Печорского бассейна и калийных рудников Старобинского месторождения. Теоретическим исследованиям в области геомеханики посвящено большое количество публикаций и научных работ, среди которых следует особо отметить работы: И.В. Баклашова, А.А. Борисова, А.П. Господарикова, В.П. Зубова, Н.Н. Кайдалова, О.В. Ковалева, В.Г. Лабазина, Н.М. Проскурякова, К.В. Руппенейта, Г.Н. Савина, А.Н. Ставрогина, А.К. Черникова, И.Л. Черняка, Е.И. Шемякина и др. Очевидно, что подавляющее число задач геомеханики по оценке НДС массива решить с использованием классических методов механики сплошных сред не представляется возможным. Для исследуемого объекта – массива, вмещающего горные выработки, – характерны неоднородность, трещиноватость, анизотропность. Аналитическое описание контуров полостей, постановка граничных условий на них достаточно сложны, соответственно, получение решений в замкнутой форме, удобной для инженерной практики, весьма затруднительно. Поэтому, в таких случаях целесообразно применять различные приближенные методы, позволяющие получать решения за счет некоторого упрощения геометрии формы исследуемого объекта, понижение порядка соответствующей краевой задачи геомеханики, на основе применения, например, вариационных методов. Исходные модели, схематически отражающие условия формирования в массиве полей (полей напряжений, деформаций и перемещений) и используемые в работе, представлены на рис.1-3. Для выработок различной геометрии определяются компоненты указанных полей, и выполняется оценка необходимых параметров крепи или оценка устойчивости незакрепленных сечений. Например, задачей такого типа является задача об определении НДС неподкрепленной выработки круглого сечения в однородном и изотропном массиве (рис.1) Рис.1. Однородный массив горных пород, вмещающий подземную горную выработку круглого сечения Широкий круг проблем управления механическими процессами в массиве горных пород, связан с очистными выработками и, как следствие, выемкой полезного ископаемого из недр, то есть наличием в пределах выемочных участков выработанных пространств (рис.2). Результаты решения таких задач геомеханики позволяют не только обеспечивать безопасность работ в очистных забоях, но и определить выбор оптимальных технологических решений. Это достигается исследованием НДС пород непосредственной и основной кровли пластов, представленной неоднородными элементами массива. Рис.2. Неоднородный слоистый массив в зоне влияния выработанного пространства без закладкиРис. 3. Неоднородный слоистый массив в зоне влияния выработанного пространства с закладкой Приведены исследования НДС призабойных зон пласта, определяющих условия управления состоянием массива непосредственно в лаве, и эффективности выемки полезного ископаемого добычными комплексами (рис.3).^ Во второй главе раскрывается основные принципы применения метода конечных элементов для получения численных решений важных прикладных задач геомеханики. Приводятся расчеты некоторых важных прикладных задач механики горных пород на основе разработанной вычислительной программы. В этой главе также проводится широкое сопоставление результатов, полученных расчетов по численным и аналитическим методам для случая одиночной выработки круглого сечения. Многие прикладные задачи геомеханики тесно связаны с определением НДС массива горных пород, вмещающего выработки круглого сечения. Для решения подобных задач механики горных пород при помощи разработанной вычислительной программы необходимо знание не только конструктивных особенностей выработки, механических свойств массива горных пород (E, ), но и вид потенциала . Так, для упругой среды (грунтового массива) в работе в качестве потенциала принято выражение вида:, где - первый и второй инвариант тензора деформации, соответственно; - константа Ламе; - коэффициент Пуассона. Формулы для сжимающих подземную горную выработку напряжений определяются по формулам:, , ,где - компоненты тензора напряжений; - объемный вес массива пород; - коэффициент бокового распора; - глубина заложения выработки. Отметим, что для простоты расчетов иногда принимается равной единице, что равносильно равномерности распределения напряжений вокруг круглой выработки. На практике рассматриваются и другие значения , например, . Для подтверждения эффективности разработанной вычислительной программы была решена классическая задача по определению НДС неподкрепленной выработки круглого сечения в однородном и изотропном массиве (рис.1), физико-механические свойства которого заранее определены. Поскольку указанная задача является осесимметричной, то рассматривалась четверть выработки. После задания входных параметров вычислительная программа автоматически разбивает область на конечные элементы и выводит разбиение на экран. Выходными параметрами программы являются компоненты вектора смещений, тензоров деформаций и напряжений слоев, прилегающих к выработке массива горных пород. Последние, в конечном итоге, и определяют как безопасность ведения горных работ, так и необходимые мероприятия по закреплению подземных горных выработок. Следует отметить, что полученные результаты достаточно точно коррелируют с известными решениями подобных задач, выполненных Г. Киршем и Г.Н. Савиным. Решение задачи о деформации массива горных пород, включающего одиночную горную выработку круглого сечения, по численно-аналитическому методу было также проведено в работах А.П. Господарикова. Таким образом, разработан эффективный алгоритм численного решения прикладных задач геомеханики на основе метода конечных элементов, реализованный в вычислительной программе.^ Третья глава посвящена вопросам, связанным с приложением вариационного метода В.З. Власова в рамках упругой модели для исследования НДС пологозалегающего массива горных пород в двумерной постановке. Приводятся расчеты напряженного состояния основной кровли отрабатываемого пласта в начальный период работы лавы и в режиме установившегося движения. В этой же главе приведены примеры численного решения краевых задач горной геомеханики методом конечных разностей при помощи вычислительной программы FLAC v.4.0. Анализ произошедших обрушений основной кровли на Старобинском месторождении калийных солей за последние 20 лет показал, что при селективной выемке с закладкой подавляющее большинство (95 %) случаев происходит в полулаве со стороны массива. Так как разрушенного галита с учетом коэффициента разрыхления солей Кр=1,37 недостаточно для полной закладки выработанного пространства (заполняется около 40 %), то возникает необходимость усовершенствования схемы расположения галита в выработанном пространстве. В данной работе предлагается следующий способ управления труднообрушаемой кровлей (патент на изобретение №2312572). Рис. 4. 1 - панельный конвейерный штрек; 2 - разгружающий штрек; 3 - конвейерный штрек лавы; 4 - транспортный штрек лавы; 5 - центральный вентиляционный штрек; 6 -закладочный штрек; 7 - забойный конвейер; 8 - перегружатель; 9 - метатель; 10 -бутовые полосы; 11 - неотработанный участок шахтного поля, в котором подготавливается следующий выемочный столб; 12 - бутовая полоса отработанного столба; L1 - длина первой полулавы; L2 - длина второй полулавы; L1/3 - ширина бутовой полосы Способ осуществляется следующим образом: для подготовки выемочного столба проходят панельный конвейерный штрек 1, разгружающий штрек 2, конвейерный штрек лавы 3, транспортный штрек лавы 4, центральный вентиляционный штрек 5. Дополнительно по центру полулавы со стороны массива проходят закладочный штрек 6. Отбитая порода поступает с забойного конвейера 7 на перегружатель 8 и далее к метателю 9, при помощи которого ее забрасывают в выработанное пространство на ширину бутовой полосы L1/3. Перегружатель 8 и метатель 9 установлены под крепью сопряжения на конвейерном, вентиляционном и закладочном штреках лавы. Метатели применяются с параметрами, позволяющими возводить бутовые полосы 10 такой ширины, чтобы они перекрывали друг друга, и между ними не оставалось свободного пространства. В результате получается полная забутовка выработанного пространства первой полулавы, расположенной со стороны массива. Подготовительная выработка (панельный конвейерный штрек) с одной стороны ограничена массивом полезного ископаемого, а с другой стороны – бутовой полосой, параметры которой при данной ширине, восстанавливают свойства массива. Следовательно, данная выработка не попадает в выработанное пространство и сохраняется для дальнейшего использования при подготовке нового столба. Подготовка следующего выемочного столба производится под защитой возведенных бутовых полос без оставления межстолбового целика, при этом панельный конвейерный штрек отработанной панели используется в качестве разгружающего штрека новой панели. Такая подготовка шахтного поля позволяет производить более полное извлечение полезного ископаемого, увеличивая коэффициент извлечения на 10-15%. Существенной трудностью применения различных методов для определения НДС массива горных пород в указанном случае являются весьма приближенные геометрические параметры закладки и ее физико-механические характеристики. Однако, с привлечением явления зональной дезинтеграции горных пород эти трудности могут быть преодолены. Сущность этого явления, впервые высказанного Е.И.Шемякиным, состоит в том, что с течением времени закладка становится такой же жесткой, как и основной массив горных пород, а значит, для определения НДС кровли в указанных условиях может быть использована схема «целик-пролет-целик». В работе рассматривается плоское напряженно-деформированное состояние массива горных пород. В процессе деформации слои, слагающие массив, свободно, без образования внутренних полостей, проскальзывают относительно друг друга. Так как обжатие массива под действием собственного веса мало, то деформация в вертикальном направлении сводится к оседанию массива в целом, что позволяет выбрать координатные функции метода В.З. Власова в виде: (1) Обобщенные перемещения произвольной точки массива горных пород тогда представляются в виде: (2) где U,V – обобщенные угловые и вертикальные перемещения соответственно. Функции , i=1,2,…,n, и V1 являются искомыми, а координатные функции , i=1,2,…,n, и - известны и выбираются таким образом, чтобы удовлетворить принятым допущениям. Так как вертикальные оседания массива не зависят от координаты y, то и V=V1. Закон Гука в предположении несжимаемости массива горных пород в вертикальном направлении примет вид:, где модуль упругости E, модуль сдвига G и коэффициент Пуассона могут быть различными для различных слоев массива. С учетом конкретного вида координатных функций (1) система (3) обыкновенных дифференциальных уравнений примет вид: (3) где - цилиндрическая жесткость k-го слоя; - толщина слоев, перенумерованных от дневной поверхности; - модуль упругости, модуль сдвига, коэффициент Пуассона и объемный вес i-го слоя массива соответственно; k - коэффициент жесткости упругого основания породного массива. Массив предполагается свободным от нагрузок на дневной поверхности и опирающимся на упругое основание с коэффициентом жесткости k. Отметим, что для свободно висящей части массива k=0. Для постановки граничных условий введем величины (с учетом вида координатных функций (1)): где ^ Mk - изгибающий момент в произвольном вертикальном сечении k-го слоя, N - перерезывающая сила в том же сечении для всего массива. Отметим, что приближенное решение задач геомеханики, помимо аналитического метода, может быть обеспечено и за счет применения различных численных методов исследования НДС, например, в работе с помощью программы FLAC были определены вертикальные смещения в центре выработки (рис.5) и опорное давление на различных удалениях от забоя (рис.6) для подобного массива. Рис.5. Распределение опорного давления в начальный период работы лавы Расчет показал, что при численном решении опорное давление экспоненциально затухает по мере удаления от забоя. То же самое происходит и с деформациями. На центр исследуемой выработки, как и при аналитическом решении, приходится максимум вертикальных смещений. В обоих решениях были не учтены многие особенности: учет трения между слоями, обжатие по вертикали и другие факторы. Но при этом решение получилось пригодным для прикладных инженерных исследований. Таким образом, предложенные алгоритмы расчета слоистого массива для прогноза напряженного состояния кровли и соляного пласта в зоне очистных работ показал правомерность использования математической модели для определения главных напряжений на участках образований закладочным материалом при различных схемах бутовых полос. В настоящее время результаты данных исследований используются при отработке лав второго и при проектировании опытного участка в условиях четвертого рудников РУП ПО «Беларуськалий». Рис.6. Вертикальные смещения кровли над выработкой^ В заключении приводятся основные выводы по результатам исследований в соответствии с поставленными задачами и рекомендации по применению разработанного метода расчета для обеспечения безопасности разработки пологих месторождений полезных ископаемых.^ ПУБЛИКАЦИИ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ Господариков А.П., Зацепин М.А. Расчет напряженного состояния основной кровли соляного пласта при оптимальном выборе рациональных параметров технологических схем // Сборник трудов «Нелинейные проблемы механики и физики деформируемого твердого тела» под ред. проф. Черныха К.Ф., СПб, вып. 8, 2004. – с.31-38. Господариков А.П., Зацепин М.А. Расчет напряженного состояния кровли и угольного пласта в зоне очистных работ // Сб. трудов международной конференции школы симпозиум «Современные проблемы механики и прикладной математики», Воронеж, ВГУ, ч.1, т.1, 2004. - с.152-156. Господариков А.П., Зацепин М.А., Мансурова С.В., Васильев С.В. Математическое моделирование прогноза напряженного состояния кровли и угольного пласта в зоне очистных работ // Сборник научных трудов «Проблемы рудничной аэрологии и безопасности месторождений полезных ископаемых», Пермь, 2004. – с.51-55. Господариков А.П., Сиренко Ю.Г., Зацепин М.А. Алгоритм расчета напряженно-деформированного состояния кровли с учетом оптимального выбора параметров технологических схем в условиях Старобинского месторождения // Записки Горного института, СПб, т.170, ч.1, 2007. – с.106-110. Зацепин М.А. Прогноз напряженно-деформированного состояния пологозалегающего массива горных пород на основе численных методов. – В сб.: «Проблемы недропользования 2008». Материалы II Всероссийской молодежной научно-практической конференции по проблемам недропользования, УРО РАН, Екатеринбург, 2008. – с. 373-378. Зацепин М.А. Математическое моделирование прогноза напряженно-деформированного состояния пологозалегающего массива горных пород // Вестн. С.-Петерб. ун-та, сер.10: Прикладная математика, информатика, процессы управления, вып. 1, 2009. - с.68-73.