МЕЛЕШКО Андрей Викторович МЕТОДИКА РАСЧЕТА НАПРЯЖЕННО- ДЕФОРМИРОВАННОГО СОСТОЯНИЯ НЕОДНОРОДНОГО СЛОИСТОГО МАССИВА
advertisement
На правах рукописи МЕЛЕШКО Андрей Викторович МЕТОДИКА РАСЧЕТА НАПРЯЖЕННОДЕФОРМИРОВАННОГО СОСТОЯНИЯ НЕОДНОРОДНОГО СЛОИСТОГО МАССИВА ПРИ ОТРАБОТКЕ ПОЛОГИХ УГОЛЬНЫХ ПЛАСТОВ Специальность 25.00.20 – Геомеханика, разрушение горных пород, рудничная аэрогазодинамика и горная теплофизика Автореферат диссертации на соискание ученой степени кандидата технических наук САНКТ-ПЕТЕРБУРГ 2012 Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Санкт-Петербургском государственном горном университете. Научный руководитель – доктор технических наук, профессор Господариков Александр Петрович Официальные оппоненты: Зубков Виктор Васильевич доктор технических наук, Национальный минерально-сырьевой университет «Горный», главный научный сотрудник лаборатории геодинамики научного Центра геомеханики и проблем горного производства Никишин Даниил Юрьевич кандидат технических наук, ООО «ГорноХимический инжиниринг», начальник отдела подземных горных работ Ведущее предприятие– ООО «СПб-Гипрошахт». Защита диссертации состоится 31 мая 2012 г. в 15 ч на заседании диссертационного совета Д 212.224.06 при СанктПетербургском государственном горном университете по адресу: 199106, Санкт-Петербург, 21-линия, д. 2, (boguslEI@yandex.ru), ауд. 1166. С диссертацией можно ознакомиться в библиотеке СанктПетербургского государственного горного университета. Автореферат разослан 28 апреля 2012 г. УЧЕНЫЙ СЕКРЕТАРЬ диссертационного совета д.т.н., профессор Э.И. БОГУСЛАВСКИЙ ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность темы диссертации. Многообразие горногеологических условий залегания пластов и увеличение глубины отработки требует постоянного мониторинга и прогнозирования напряженно-деформированного состояния (НДС) породных массивов, вмещающих горные выработки различного назначения и очертания, базирующегося, как правило, на эффективных численных методах решения. Многие прикладные задачи горной геомеханики связаны с определением НДС техногенно нарушенного массива горных пород (МГП). Такие задачи решались многими исследователями различными методами математического и физического моделирования. Математическое моделирование имеет определенные преимущества перед физическим, поскольку обладает наибольшей общностью как при описании сущности геомеханических процессов, так и дает возможность исследовать и прогнозировать последние в наиболее широком спектре их определяющих параметров. В связи с развитием вычислительной техники и методов математического моделирования, наряду с традиционными (широко апробированными) аналитическими методами, все шире применяются численные методы: метод конечных разностей (МКР), метод конечных элементов (МКЭ), метод граничных элементов (МГЭ). Эффективное применение этих методов для решения важных прикладных задач горной геомеханики определяется не только возможностями используемого программного комплекса, но и наличием соответствующей методической базы решения подобных задач. Существенный вклад в теорию и практику горногеомеханического обеспечения горных работ при отработке пластовых месторождений внесли такие ученые и специалисты, как Ардашев К.А., Борисов А.А., Бич Я.А., Ковалев О.В., Комиссаров С.В., Лабазин В.Г., Линьков А.М., Петухов И.М., Протосеня А.Г., Проскуряков Н.М., Слесарев В.Д., Шик В.М. и др. В то же время механизм деформирования, разрушения и оседания горных пород при работе высокопроизводительных очистных забоев имеет особенности, требующие дополнительных исследований, поскольку существенно влияет на состояние целиков и подготовительных выработок – особенно при парной (многоштрековой) подготовке выемочных столбов. Получение надежных прогнозных оценок механических состояний массива (включая компоненты тензоров напряжения и деформации и вектора перемещений), учитывающих временной фактор, позволит своевременно предотвращать опасные проявления горного давления в подготовительных выработках, повысить безопасность ведения горных работ и минимизировать затраты на проведение и эксплуатацию горных выработок. Решение указанной задачи требует дальнейшего эффективного использования методов механики сплошных сред (МСС) с соответствующим отображением процессов обрушения горных пород и оседания подработанной толщи и, соответственно, разработки специальной методики, базирующейся на данных практики (маркшейдерские наблюдения, специальные исследования процессов деформирования и обрушения подработанной толщи, механических характеристик массивов обрушенных пород и т.п.). Все это и обуславливает актуальность темы исследования, а так же структуру и содержание работы. Цель диссертационной работы. Разработка методики расчета НДС неоднородного слоистого массива, основанной на применении эффективных численных методов и позволяющей повысить надежность прогноза «горно-геомеханической» обстановки на выемочном участке (состояние целиков и участковых выработок), а так же безопасность отработки угольных пластов. Основные задачи исследования: 1.Обосновать выбор горно-геомеханической модели МГП, адекватной по своим основным свойствам реальному массиву для условий отработки угольных пластов Ерунаковского месторождения. 2.Разработать методику численного моделирования (на основе МКЭ) процесса деформирования и обрушения горных пород в выработанном пространстве при отработке пологих угольных пластов. 3.Установить закономерности изменения НДС неоднородного слоистого МГП, вмещающего выработки выемочного участка, на различных этапах отработки угольного пласта с учетом нелинейного процесса деформирования горных пород для условий отработки угольных пластов Ерунаковского месторождения. Идея диссертационной работы. Математическое моделирование НДС неоднородного слоистого массива при управлении кровлей полным обрушением на пологих пластах необходимо проводить на основе разработанного алгоритма, учитывающего закономерности обрушения пород в выработанном пространстве и формирования больших массивов с новыми механическими свойствами. Методы исследований. Для решения поставленных задач использован комплексный метод исследований: анализ и обобщение данных, опубликованных в научной и горнотехнической литературе по проблемам устойчивости участковых подготовительных выработок; методы современного математического моделирования; шахтные исследования процессов деформирования и разрушения краевых частей целиков, форм обрушений пород непосредственной кровли в подготовительных выработках; регистрация и анализ величин смещений почвы и кровли; сравнительный анализ и сопоставление результатов моделирования с данными натурных измерений. Научная новизна диссертационной работы: установлены закономерности изменения напряженного состояния неоднородного породного массива, вмещающего комплекс выработанных пространств и участковых выработок, в рамках применимости физически линейного процесса деформирования горных пород; установлены закономерности изменения напряженного состояния целиков и краевых частей горного массива, вмещающего участковые выработки, на различных этапах отработки выемочных столбов в рамках физически нелинейного процесса деформирования горных пород; Основные защищаемые положения: 1. Моделирование НДС МГП, нарушенного отработкой длинных столбов, целесообразно осуществлять на однородных изотропных геомеханических моделях с последующей оценкой влияния полученного НДС на функции геометрических (мощность), механических (модуль деформации, коэффициент Пуассона) и контактных параметров слагающих слоев непосредственной и основной кровли до момента подбучивания вышележащих слоев обрушенными породами. 2. Моделирование геомеханических процессов обрушения МГП, подработанного длинными столбами, должно производиться с учетом оценки высоты обрушения и деформационных характеристик пород в выработанном пространстве, как функций коэффициента разрыхления обрушенных пород кровли. 3. Моделирование процессов оседания подработанных пород кровли вплоть до поверхности должно учитывать изменчивость деформационных свойств ответственных элементов слоистого МГП как функций трещиноватости последнего, а также использовать нелинейные зависимости физического закона (пошагово) в характерных зонах с сопоставлением с данными маркшейдерских наблюдений. Практическая значимость работы заключается в разработке методики расчета НДС МГП на выемочном участке, позволяющей учесть изменение механических характеристик краевых частей массива и обрушенных пород в выработанном пространстве. Результаты диссертационной работы переданы в ООО «СПб-Гипрошахт» для использования при проектировании и используются для горно-геомеханического обоснования проектных решений. Достоверность и обоснованность научных положений и рекомендаций подтверждается большим объемом проанализированной информации по математическим моделям и методам численного моделирования НДС МГП, корректностью выполненных экспериментально-аналитических исследований и удовлетворительной сходимостью их с результатами численного моделирования, широким сопоставлением результатов по разработанной методике как с натурными данными, так и с результатами, полученными другими авторами и по другим методикам. Апробация работы. Основные положения диссертационной работы докладывались на XXXVIII международной научной конференции аспирантов и студентов «Процессы управления и устойчивость» (Санкт-Петербург, 2008 г.); научном симпозиуме «Неделя горняка-2010» (Москва, 2010 г.); на ежегодных научных конференциях молодых ученых СПГГУ «Полезные ископаемые России и их освоение» (Санкт-Петербург, 2009, 2010 гг.), семинарах кафедры РМПИ СПГГУ. Личный вклад автора заключается в постановке задач исследований, разработке методики проведения исследований, анализе геологических и горнотехнологических условий отработки угольных пластов Ерунаковского месторождения, проведении численного моделирования НДС углевмещающего массива для различных участков по длине выемочного столба, выполнении анализа результатов по разработанной методике моделирования НДС МГП с пошаговым изменением геометрии и механических характеристик различных элементов последнего, отображающей процессы разрушения, обрушения и оседания горных пород, выявлении функциональной взаимосвязи между механическими параметрами обрушенных пород в выработанном пространстве и степенью их уплотнения. Публикации. По теме диссертации опубликовано 8 печатных работ, из них 4 в журналах, рекомендованных ВАК. Структура и объём работы. Диссертационная работа изложена на 145 страницах машинописного текста, состоит из введения, четырех глав, заключения, списка литературы из 98 источников, приложения, включает 104 рисунка и 12 таблиц. ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ В первой главе диссертационной работы отражены состояние и перспективы развития методов моделирования геомеханических процессов в недрах при подземной разработке угольных пластов, выполнен сравнительный анализ программных продуктов, разработанных на основе современных численных методов, проанализированы технологии отработки угольных пластов Ерунаковского месторождения и сформулированы основные требования к научному обоснованию безопасных и рациональных параметров используемых технологических схем, сформулированы цель и задачи исследований. Во второй главе проанализированы и обобщены горногеологические условия отработки угольных пластов на шахтах Ерунаковского месторождения, приведены материалы по обобщению натурных данных о механических характеристиках различных литотипов горных пород и угольных пластов для рассматриваемого района Кузбасса, приведены разработанные горно-геомеханические модели и расчетные схемы для моделирования НДС МГП в пределах выемочного участка, изложены методические рекомендации по расчету механических характеристик различных элементов МГП при моделировании НДС для случая отработки угольных пластов длинными столбами и многоштрековой подготовке. В третьей главе приведены результаты моделирования НДС МГП вмещающего угольные пласты в пределах выемочного участка для различных моментов отработки выемочного столба; проанализирована динамика изменения НДС в зависимости от длины лавы, ширины охранного угольного целика, мощностей непосредственной и основной кровли, мощности обрушения пород в выработанном пространстве, степени разрыхления пород, обрушенных в выработанное пространство, наличия или отсутствия факта «подбучивания» вышележащих пород, приведены принципиальные подходы к расчету механических характеристик ответственных элементов массива при моделировании процесса оседания подработанной толщи пород от непосредственной кровли разрабатываемого пласта до поверхности (наносов). В четвертой главе приведены результаты шахтных исследований по выявлению основных горно-геологических и горнотехнических факторов, влияющих на состояние участковых подготовительных выработок, обоснована методика шахтных исследований, используемая для сопоставления результатов математического моделирования и натурных наблюдений, изложены методические рекомендации по расчету механических характеристик различных элементов МГП при моделировании НДС для случая отработки угольных пластов длинными столбами и многоштрековой подготовке. Основные результаты исследований отражены в следующих защищаемых положениях: 1. Моделирование НДС МГП, нарушенного отработкой длинных столбов, целесообразно осуществлять на однородных изотропных геомеханических моделях с последующей оценкой влияния полученного НДС на функции геометрических (мощность), механических (модуль деформации, коэффициент Пуассона) и контактных параметров слагающих слоев непосредственной и основной кровли до момента подбучивания вышележащих слоев обрушенными породами. Проблемы, возникающие при отработке пластовых месторождений полезных ископаемых, свидетельствуют о необходимости дальнейшего развития горной геомеханики. Совершенствование уже известных методов изучения НДС МГП тесно связано с необходимостью обоснованного принятия решений в производственных вопросах, а также обеспечения безопасности на горном предприятии как с технической, так и с экологической точки зрения. Очевидно, что изучение геомеханических процессов в МГП не может производиться с использованием только натурных данных, результатов лабораторных исследований или аналитических расчетов. К тому же постановка такой задачи нецелесообразна, так как для каждого из приведенных методов исследования «характерен» только свой класс задач. С точки зрения получения достоверного прогноза горно-геомеханических условий подземной разработки месторождений полезных ископаемых наиболее значимыми являются два метода: метод «поиска месторождений аналогов» и математическое моделирование НДС МГП. Наиболее информативными и сравнительно малозатратными методами математического моделирования являются численные методы. Среди численных методов наиболее эффективными являются МКР, МКЭ, МГЭ, метод граничных интегральных уравнений и другие методы, получившие большое распространение в инженерной и научной среде при высоком развитии вычислительной техники. Несмотря на эффективность численных методов при решении задач горной геомеханики, при использовании последних, зачастую, приходится прибегать к различного рода модификациям. В частности, при исследовании геомеханических процессов в неоднородных, многосвязных, трещиноватых и анизотропных МГП, учет таких явлений как разрушение, обрушение, пригрузка и уплотнение значительных объемов массива при управлении кровлей полным обрушением, требует, как правило, такой модификации. Широкий круг проблем управления геомеханическими процессами в МГП связан с очистными выработками (лавами) и наличием в пределах выемочных участков значительных выработанных пространств. Решение таких задач обеспечивает как безопасность работ в таких забоях, так и определяет выбор рациональных технологических решений. Исследование НДС массива неоднородных пород непосредственной и основной кровли пластов, а также НДС призабойных зон пласта и вмещающих пород, определяет условия управления состоянием массива непосредственно в лаве и участковых выработках, а следовательно, влияет на эффективность выемки полезного ископаемого. Большую практическую значимость имеет оценка эффекта влияния выработанного пространства на расположенные в краевых частях массива горные выработки (систему целиков и выработок). Расчет НДС в этом случае обеспечивает возможность оценки удароопасности целиков и разработки комплекса мероприятий для исключения возникновения таких явлений на выемочных участках. В настоящее время на смену бесцеликовым схемам отработки выемочных столбов и использованию рамных крепей в участковых выработках пришли схемы отработки пластов с охраной выработок целиками угля и анкерным креплением кровли, боков, а в некоторых случаях и почвы выработок. Использование современных механизированных комплексов позволяет достигать высоких нагрузок на очистной забой, и соответственно, существенных скоростей подвигания лавы. Задача поддержания участковых выработок на границе с выработанным пространством в случае использования анкерной крепи или значительной мощности отрабатываемого пласта становится практически невыполнимой. В связи с этим на шахтах ОАО «СУЭК-Кузбасс» и получили распространение технологические схемы с охраной участковых выработок целиками угля. На границе с выработанным пространством такая выработка попадает в зону интенсивно протекающих геомеханических процессов, связанных не только с обрушением непосредственной и основной кровли, но и с довольно длительными процессами сдвижения всей налегающей толщи, и как следствие этих процессов – выработка находится в зоне интенсивного трещинообразования, перехода части пласта в «запредельное состояние», что приводит к перераспределению напряжений в пласте и налегающей толще и так до момента оседания всего массива пород до поверхности. Основными геомеханическими проблемами при применении данной технологической схемы являются выбор ширины целика, оценка его удароопасности и расчет параметров анкерного крепления и крепи усиления для различных «этапов эксплуатации» участковых выработок. Схематически «этапы эксплуатации» подготовительных выработок при подготовке выемочных столбов двумя выработками приведены на рис. 1. Условно можно выделить следующие этапы: проведение выработок, время до подхода первой лавы, зона повышенного горного давления от выработанного пространства первой лавы, время до подхода второй лавы, зона опорного давления второй лавы. Рис. 1. Общая схема для выделения характерных сечений для оценки НДС массива в течение эксплуатации участковых выработок Оценка напряжений и перемещений, возникающих в МГП на всех этапах эксплуатации выработок, является основой для оценки состояния угольного пласта и кровли выработки и, соответственно, параметров ее крепления и способов охраны. Такую оценку можно получить только с помощью численных методов. В работе оценивались результаты применения различных программных пакетов по определению задач о НДС неоднородного слоистого массива. На основе анализа горно-геологических условий отработки пластов Ерунаковского месторождения была построена горногеомеханическая модель исследуемого МГП, представленная на рис.2. Численное моделирование осуществлялось в рамках плоской деформации для «тяжелой» полуплоскости. В первой расчетной схеме моделировалось наличие свободного от обрушенных пород выработанного пространства. Рис. 2. Горно-геомеханическая модель МГП На первом этапе была решена задача о степени влияния неоднородности (мощности слоев, их деформационных характеристик и контактных условий) на распределение напряжений в краевой части пласта и на уровне основной и непосредственной кровли (при различных пролетах зависающих слоев основной кровли). На рис. 3. Представлены графические зависимости вертикальных напряжений (по кровле пласта, на уровне верхней границы основной кровли и на уровне верхней границы непосредственной кровли) в краевой части угольного пласта. Из представленных зависимостей следует, что основными факторами, определяющими напряженное состояние краевой части пласта и его кровли, являются пролет выработанного пространства и контактные условия на границе непосредственной и основной кровли. При этом уровни вертикальных напряжений остаются неизменными, а горизонтальные напряжения в характерных зонах существенно изменяются, переходя из области сжатия в область растяжения. 2. Моделирование геомеханических процессов обрушения МГП, подработанного длинными столбами, должно производиться с учетом оценки мощности обрушения и деформационных характеристик пород в выработанном пространстве, как функций коэффициента разрыхления обрушенных пород кровли. Численные методы расчета НДС МГП не позволяют в полной мере моделировать процессы трещинообразования, разрушения и обрушения подработанных пород. Поэтому, для моделирования этапов обрушения пород в выработанное пространство и процессов оседания необходимо оценить геометрию обрушенного пространства, а так же деформационные свойства обрушенных пород. Для принятого при отработке пластов Ерунаковского месторождения способа управления горным давлением (УГД) в лавах – полным обрушением, в целом, можно констатировать, что обрушенные породы (непосредственной и основной кровли) могут «подбучивать» или «не подбучивать» вышележащие породы. Условие подбучивания при разрыхлении пород непосредственной кровли можно оценить из выполнения неравенства: hнк (ксрнк 1) m mc , (1) где hнк – суммарная мощность обрушенных пород, м; m – вынимаемая мощность пласта, м; mс – величина предельного свободного опускания пород основной кровли: mс (0,10,4)m, м; ксрнк – среднее значение разрыхляемости пород непосредственной кровли: kср 1,151,3, достигая в некоторых случаях значения 1,45. Для построения расчетных схем, отражающих процесс оседания подработанных пород, после определения высоты обрушения пород кровли до момента «подбучивания», необходимо определить высоту обрушения пород и оценить механические (деформационные) характеристики формируемого в выработанном пространстве материала. Последний представлен породами непосредственной и основной кровли. Коэффициенты разрыхления обрушаемых слоев кровли, соответстующие определенным литотипам горных пород, приняты равными для аргиллитов -1,4, алевролитов: 1,2-1,4; песчаников: 1,05-1,1. На базе лаборатории ООО «ИК «КузНИУИ» были проведены исследования по определению деформационных характеристик разрыхленных пород при варьировании крупности кусков и исследовании различных литотипов пород (рис.4). Е, МПа , МПа Рис. 4. Зависимость модуля деформации пород выработанного пространства от действующих вертикальных напряжений (аргилит и алевролит различных фракций с коэффициентом разрыхления, равным 1,4) Характерным этапом моделирования является момент подбучивания обрушенными породами вышележащих пород. Исходной информацией для таких расчетов являются сведения о геологическом строении кровли рассматриваемого пласта и на этой основе принимаются значения коэффициентов разрыхления. О подбучивании вышележащих пород свидетельствует выполнение неравенства: (2) hi (m mc ) hi ki . i i В этом случае моделирование обрушенных пород в выработанном пространстве осуществляется в виде слоев с определенными механическими характеристиками, пригруженных вышележащими слоями обрушенных пород. Взаимосвязь деформационных характеристик обрушенных пород от коэффициента разрыхления принята в виде полученных эмпирических зависимостей: Е 11780 k p 19 0,1750 k p , (3) 2 , 36 На рис. 5 представлены зависимости вертикальных и горизонтальных напряжений (по пласту, по верхней границе основной кровли и по верхней границе непосредственной кровли) в краевой части угольного пласта на момент подбучивания вышележащих пород разрушенными породами кровли пласта. Как видно из представленных зависимостей, основным фактором, определяющим напряженное состояние краевой части пласта и его кровли на момент подбучивания вышележащих пород, является длина пролета выработанного пространства. Характер распределения напряжений в краевой части пласта полностью идентичен предыдущему распределению (отличие только в высоте выработанного пространства). 3. Моделирование процессов оседания подработанных пород кровли вплоть до поверхности должно учитывать изменчивость деформационных свойств ответственных элементов слоевого МГП как функций трещиноватости последнего, а также использовать нелинейные зависимости физического закона (пошагово) в характерных зонах с сопоставлением с данными маркшейдерских наблюдений. На момент подбучивания обрушенными породами вышележащих слоев кровли опорное давление достигает максимума и остается относительно выдержанным по мощности от почвы пласта до верхней границы обрушенных пород. На рис. 6 представлены вертикальные перемещения МГП на момент беспорядочного обрушения пород кровли пласта. Как видно, упругие деформации пласта и вмещающего массива незначительны, то есть процесс оседания пород налегающей толщи невозможно описать в рамках упругой модели без соответствующего изменения деформационных характеристик краевой части пласта, а также налегающих слоев кровли в краевой части мульды сдвижения. На рис. 7 представлена схема оседания земной поверхности с указанием углов полных сдвижений и граничных углов. Рис. 7. Схема элементов мульды сдвижения налегающей толщи и поверхности в приконтурных зонах массива Полученные данные о распределении напряжений и деформаций служат основанием для оценки изменчивости модуля деформаций и коэффициента Пуассона в характерных зонах массива. Анализ указанных полей показал, что характерными зонами в МГП при выемке угольных пластов являются краевая часть пласта как по вертикальным напряжениям, так и по горизонтальным напряжениям, а также линии раздела слоев, отвечающие наименьшему сцеплению на контакте – по горизонтальным деформациям. Данные о полях напряжений и деформаций были пошагово использованы для расчета деформационных параметров характерных зон массива. Ширина зон с измененными деформационными свойствами массива принималась, исходя из геологического строения слоевой толщи и составляла от 5 до 15 м. Необходимые функциональные зависимости деформационных характеристик принимаются в виде Е, = f(Е;Т); (4) Е, =f1(Е;Т); Е=f2(Е;Т), где Е,, Е - модули деформации трещиноватого МГП в вертикальном и горизонтальном направлениях, соответственно; Е- модуль деформации по результатам лабораторных испытаний; Тпараметр, отражающий трещиноватость МГП, и являющийся функцией вида Т=(n; ; ). (5) Здесь n - количество систем трещин (функция напряжений и паспорта прочности); - комплексная геометрическая характеристика для i-й системы трещин; - угол i-й системы трещин с горизонтом. На рис. 8 представлено поле вертикальных перемещений согласно расчетной схеме, отражающей сформировавшиеся зоны массива со сниженными деформационными характеристиками. Как видно из представленного рисунка оседание, земной поверхности в пределах плоского дна мульды сдвижения соответствует мощности отрабатываемого пласта. Таким образом, разработанная методика позволяет рассчитывать напряженно-деформированное состояние МГП на различных стадиях эксплуатации выработок выемочного участка и оценивать удароопасность целиков угля и условия поддержания подготовительных выработок в течение всего срока их службы. ЗАКЛЮЧЕНИЕ Диссертация представляет собой законченную научноквалификационную работу, в которой предлагается решение актуальной задачи разработки методики моделирования НДС слоистого массива горных пород при отработке угольных пластов высокопроизводительными очистными забоями с многоштрековой подготовкой выемочных столбов, имеющей существенное значение при подземной разработке пологих пластовых месторождений. Основные научные и практические результаты: 1. Разработана горно-геомеханическая модель слоистого углевмещающего МГП для условий Ерунаковского месторождения. 2. На основе моделирования горно-геомеханических процессов установлены закономерности изменения параметров НДС слоистого углевмещающего массива, в частности, компонентов тензоров деформации и напряжения, а также вектора перемещений, техногенно нарушенного отработкой длинных столбов. 3. Установлена зависимость коэффициентов концентрации напряжений в краевых частях выемочных столбов (зависимость вертикальных напряжений) от длины лавы, строения и мощности непосредственной и основной кровли, контактных условий на границе раздела ответственных элементов массива. 4. Выявлено, что основными факторами, определяющими нагрузку на краевые части массива, являются величина пролета выработанного пространства (длина лавы) и степень подбучивания в выработанном пространстве вышележащих слоев обрушенными породами. 5. Установлен качественный характер распределения опорного давления в краевых частях МГП в зависимости от ширины охранного целика. 6. Разработана методика расчета высоты зоны беспорядочного обрушения горных пород изначения коэффициента разрыхления пород в выработанном пространстве и его взаимосвязь с деформационными характеристиками обрушенных пород, которые формируют «отпор» в выработанном пространстве лав при сдвижении подработанного МГП. 7. Разработана методика пошаговой оценки НДС МГП при отработке длинных столбов для характерных моментов: «до обрушения основной кровли», «до подбучивания» налегающих пород, «до образования плоского дна мульды сдвижения», на момент полных оседаний подработанной толщи. 8. На основе анализа динамики изменения НДС МГП для характерных этапов отработки столба получена методика расчета деформационных параметров ответственных элементов МГП от факторов действующих напряжений и деформаций. 9. Разработана методика моделирования НДС МГП в течение оседания всех подработанных слоев – от момента обрушения непосредственной кровли до полных оседаний земной поверхности. Основные положения диссертации опубликованы в следующих работах: в изданиях, входящих в Перечень ВАК Минобрнауки России: 1. Мелешко А.В. Численные решения нелинейных краевых задач геомеханики/ А.П. Господариков, М.А.Зацепин, А.В.Мелешко // Записки Горного института. Т. 196. СПб., 2012. С.306-310 2. Мелешко А.В. Оценка удароопасности междуштрековых целиков при отработке запасов угля шахты «Котинская// Горный информационно-аналитический бюллетень (научно-технический журнал). Отдельные статьи (специальный выпуск).-2012.-№3.-М.: издательство «Горная книга». С. 3-9. 3. Мелешко А.В. Численное моделирование на основе метода конечных разностей некоторых прикладных задач геомеханики/ А.П. Господариков, М.А. Зацепин, А.В Мелешко //Записки Горного института. Т. 182. СПб. 2009. С. 238-240. 4. Мелешко А.В. Совершенствование селективной выемки мощных калийных пластов при камерной системе разработки/ Ю.Г. Сиренко, А.В. Мелешко// Горный информационно-аналитический бюллетень, № 1, М., 2010. С. 386-389. в прочих изданиях: 5. Мелешко А.В. Оценка удароопасности междуштрековых целиков. В сб.: «Проблемы недропользования 2010». Материалы IV Всероссийской молодежной конференции, УрО РАН, Екатеринбург, 2010. С. 496-504. 6. Мелешко А.В. Математическое моделирование прогноза напряженно-деформированного состояния пологозалегающего массива горных пород / А.П. Господариков, М.А. Зацепин, А.В Мелешко// 6-я Межрегиональная научно-практическая конференция «Освоение минеральных ресурсов Севера. Проблемы и решения», Воркута , 2008.- с. 174-178. 7. Мелешко А.В. Об одном подходе математического моделирования важных геомеханических задач. В сб.: «Процессы управления и устойчивость». Труды XXXVIII Международной научной конференции аспирантов и студентов, СПб., 2008. С.181186. 8. Meleshko A. Modeling stress and strain state roofs and pillars for development bed deposits// Materialy XLVIII Sesji Pionu Gorniczego, Krakow, 6 grudnia 2007. Р. 62. Рис. 3. Зависимость горизонтальных и вертикальных напряжений в краевой части массива при «зависании пород кровли» (sх1, sy1-по пласту, sх2, sy2-по непосредственной кровле, sх3, sy3-по основной кровле) Рис. 5. Зависимость горизонтальных и вертикальных напряжений в краевой части массива после обрушения пород основной и непосредственной кровли (sх1, sy1-по пласту, sх2, sy2-по непосредственной кровле, sх3, sy3-по основной кровле) Рис.6. Поле распределения вертикальных перемещений V (максимальная величина опускания 5 см по кровле и 3 см по поверхности) Рис. 8. Поле распределения полных смещений (максимальная величина опускания 4 м по поверхности)