Прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы

За допомогою комп'ютерного моделювання зроблено прогноз динаміки обвалень при відробітку примикаючої до раніше виробленого простору лави. Through the computer design prognosis of dynamics of bringing down at working off lava joining to the before produced space is done....

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Datum:2007
Hauptverfasser: Назимко, Е.В., Канин, В.А., Демченко, А.И., Ефремов, И.А.
Format: Artikel
Sprache:Russian
Veröffentlicht: Український науково-дослідницький і проектно-конструкторський інститут гірничої геології, геомеханіки і маркшейдерської справи НАН України 2007
Online Zugang:https://nasplib.isofts.kiev.ua/handle/123456789/13968
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы / Е.В. Назимко, В.А. Канин,А.И. Демченко,И.А. Ефремов // Наукові праці УкрНДМІ НАН України. — 2007. — № 1. — С. 173-189. — Бібліогр.: 15 назв. — рос.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id nasplib_isofts_kiev_ua-123456789-13968
record_format dspace
spelling Назимко, Е.В.
Канин, В.А.
Демченко, А.И.
Ефремов, И.А.
2010-12-07T13:32:06Z
2010-12-07T13:32:06Z
2007
Прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы / Е.В. Назимко, В.А. Канин,А.И. Демченко,И.А. Ефремов // Наукові праці УкрНДМІ НАН України. — 2007. — № 1. — С. 173-189. — Бібліогр.: 15 назв. — рос.
1996-885X
https://nasplib.isofts.kiev.ua/handle/123456789/13968
622.287
За допомогою комп'ютерного моделювання зроблено прогноз динаміки обвалень при відробітку примикаючої до раніше виробленого простору лави.
Through the computer design prognosis of dynamics of bringing down at working off lava joining to the before produced space is done.
ru
Український науково-дослідницький і проектно-конструкторський інститут гірничої геології, геомеханіки і маркшейдерської справи НАН України
Прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы
Article
published earlier
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
title Прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы
spellingShingle Прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы
Назимко, Е.В.
Канин, В.А.
Демченко, А.И.
Ефремов, И.А.
title_short Прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы
title_full Прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы
title_fullStr Прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы
title_full_unstemmed Прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы
title_sort прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы
author Назимко, Е.В.
Канин, В.А.
Демченко, А.И.
Ефремов, И.А.
author_facet Назимко, Е.В.
Канин, В.А.
Демченко, А.И.
Ефремов, И.А.
publishDate 2007
language Russian
publisher Український науково-дослідницький і проектно-конструкторський інститут гірничої геології, геомеханіки і маркшейдерської справи НАН України
format Article
description За допомогою комп'ютерного моделювання зроблено прогноз динаміки обвалень при відробітку примикаючої до раніше виробленого простору лави. Through the computer design prognosis of dynamics of bringing down at working off lava joining to the before produced space is done.
issn 1996-885X
url https://nasplib.isofts.kiev.ua/handle/123456789/13968
citation_txt Прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству лавы / Е.В. Назимко, В.А. Канин,А.И. Демченко,И.А. Ефремов // Наукові праці УкрНДМІ НАН України. — 2007. — № 1. — С. 173-189. — Бібліогр.: 15 назв. — рос.
work_keys_str_mv AT nazimkoev prognozdinamikiobrušeniipriotrabotkeprimykaûŝeikraneevyrabotannomuprostranstvulavy
AT kaninva prognozdinamikiobrušeniipriotrabotkeprimykaûŝeikraneevyrabotannomuprostranstvulavy
AT demčenkoai prognozdinamikiobrušeniipriotrabotkeprimykaûŝeikraneevyrabotannomuprostranstvulavy
AT efremovia prognozdinamikiobrušeniipriotrabotkeprimykaûŝeikraneevyrabotannomuprostranstvulavy
first_indexed 2025-11-26T21:32:34Z
last_indexed 2025-11-26T21:32:34Z
_version_ 1850777317091049472
fulltext Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 УДК 622.287 ПРОГНОЗ ДИНАМИКИ ОБРУШЕНИЙ ПРИ ОТРАБОТКЕ ПРИМЫКАЮЩЕЙ К РАНЕЕ ВЫРАБОТАННОМУ ПРОСТРАНСТВУ ЛАВЫ Назимко Е.В., Канин В.А. (УкрНИМИ НАНУ, г. Донецк), Демченко А.И. (шх. «Красноармейская-Западная №1» г. Красноармейск, Украина), Ефремов И.А. (шх. им. А.Ф. Засядько, г. Донецк, Украина) За допомогою комп'ютерного моделювання зроблено прогноз динаміки обвалень при відробітку примикаючої до раніше виробленого простору лави. Through the computer design prognosis of dynamics of bringing down at working off lava joining to the before produced space is done. Эффективность работы угольных шахт в условиях современной рыночной экономики определяется, как правило, интенсивностью производственного процесса (нагрузкой на очистной забой и скоростью проходческих работ). Исходя из реальных возможностей добычной техники, нагрузка на очистной забой для пластов мощностью порядка 1 м может составлять 2,5- 3,0 тыс. т в сутки, т.е. в 3-4 раза выше существующей. Однако, как показал опыт работы некоторых шахт, интенсификация технологии добычи угля провоцирует ряд негативных явлений в поведении геомеханической системы, включающей углепородный массив, горные выработки, крепь и охранные конструкции. В частности, увеличивается частота газодинамических и геодинамических явлений, пучение пород и ухудшаются условия управления горным давлением в лаве, возрастают нагрузки на крепь выемочных штреков, содержание метана в рудничной атмосфере. Научно-практические наработки, 173 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 имеющиеся в физике массива горных пород, не дают однозначного ответа, а часто и не объясняют причин происходящих явлений. Из этого следует необходимость расширения и углубления знаний о геомеханических процессах, протекающих в массиве при больших скоростях обнажения массива. В настоящее время особую актуальность приобрела проблема совершенствования теории и практики познания геомеханических процессов, протекающих в подземной угледобывающей геотехнической системе с целью создания основ управления ее состоянием. Решение этой проблемы состоит, прежде всего, в осмыслении современного состояния геомеханики углепородного массива, оценке новых горно- геологических и горно-технологических условий, выделении ключевых звеньев в функционировании пространственно- временной системы «массив-выработки-охранные конструкции», выявлении закономерностей протекания геомеханических процессов в этой системе, а также создании концептуальных основ контроля и управления горным давлением в подземных геотехнологических системах. Объем публикаций, посвященных динамическим процессам, сопровождающим активную стадию сдвижений в отечественной литературе, вряд ли составит несколько процентов. Это свидетельствует о недостаточной проработанности данного направления. В качестве типичных примеров можно назвать публикации УкрНИМИ [1], ДонНТУ [2, 3], МакНИИ [4] и других институтов. Сложность проблемы динамических проявлений сдвижения горных пород при их подработке быстро движущимися очистными забоями иллюстрируется тем, что даже сегодня влияние изменения скорости подвигания очистного забоя на геомеханические процессы, протекающие в массиве при больших скоростях обнажения массива, оценивается неоднозначно. Ряд исследователей, основываясь на наблюдениях за сдвижением горных пород и земной поверхности под влиянием горных разработок, указывают, что увеличение скорости подвигания 174 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 очистных работ во многих случаях уменьшает деформации земной поверхности. Целый ряд специалистов придерживаются противоположного взгляда на эффект скорости подвигания очистного забоя. Так Р.А. Муллер [5], А.Н. Медянцев [6], М.А. Иофис [7], В.Н. Земисев [8], В.С. Костенич [9], Е.Г. Петрук [10], В.И. Мякенький [11], А. Ковальски [12] и др. установили, что скорость подвигания очистного забоя не оказывает существенного влияния на величины деформаций земной поверхности, а увеличивает или уменьшает продолжительность и интенсивность процесса сдвижения. Указанные работы основаны на прикладном подходе к изучению параметров мульды сдвижений. Однако если рассматривать динамические проявления процесса сдвижения не только с точки зрения устойчивости наземных сооружений, но и как важные индикаторы параметров этого процесса в динамике, следует сделать вывод о преждевременном заключении о несущественной роли подобных проявлений. Таким образом, анализ существующих представлений о динамических проявлениях сдвижений подрабатываемого горного массива на их активной стадии показывает, что скорость подвигания очистного забоя существенно влияет на геомеханические процессы, протекающие в массиве. Механизм этого влияния действует через удлинение зависающих консольных слоев горных пород, в результате чего перемещаются точки перегиба кривой оседания, уменьшаются их наклоны и текущие величины самих оседаний, изменяется траектория точек земной поверхности в процессе оседаний. Зависания прочных породных слоев, не склонных к ползучести, приводят к их неизбежному обрушению. Полости зависания являются потенциальными емкостями для накопления взрывоопасного газа метана, что представляет опасность для рабочих очистного забоя. Поэтому является целесообразным сделать прогноз динамики обрушений при отработке примыкающей к ранее выработанному пространству 13-й восточной лавы шахты им. А.Ф. Засядько по пласту l1, так как шахта им. А.Ф. Засядько работает на большой глубине (1150 м) и при высоких скоростях подвигания. 175 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 В данном случае прогноз динамики обрушений сделан при помощи математического моделирования. Для решения поставленной задачи использовался метод конечных разностей для прогноза напряженного состояния подрабатываемого породного массива, который позволил получить распределение высоты зоны полных обрушений основной кровли, а также распределение эквивалентных напряжений. Выкопировка с плана горных выработок представлена на рис. 1. 0 50 Рис. 1. Выкопировка с плана горных выработок по пласту l1 в районе 13-й восточной лавы шахты м. А.Ф. Засядко Пласт l1 отрабатывается 13-й восточной лавой на глубине 1150 м. Лава длиной 240 м отходила от монтажного ходка шириной 6 м, который пройден в целике и не примыкает к ранее выработанному пространству 12 восточной лавы, а расположен 176 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 на расстоянии 25 м от угловой части у нижней границы выработанного пространства. Управление кровлей осуществляется путем ее полного обрушения с помощью механизированной крепи 3КД-90. Вынимаемая мощность пласта составляет 1,9-2,0 м. Пласт сложного строения залегает в устойчивых вмещающих породах (рис. 2). Рис. 2. Стратиграфическая колонка толщи в окрестности 13-й восточной лавы пласта l1 шахты им. А.Ф. За- сядько Вблизи монтажного ходка 13-й восточной лавы расположен замок флексурой антиклинальной складки, ось которой указана на левой части рис. 3 штрихпунктирной кривой. 177 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 Ось антиклинали Рис. 3. Расположение замка флексурой антиклинальной складки в районе 13-й восточной лавы пласта l1 шахты им. А.Ф. Засядько Закономерности перераспределения напряженно - деформированного состояния вокруг очистного забоя 13-й восточной лавы исследовалось с помощью компьютерного моделирования. С его помощью исследовалось пространственное распределение горного давления, а также оседаний основной кровли и ее предельного состояния. Подрабатываемая толща горных пород рассматривается как толстая плита, опирающаяся на неоднородное основание 3 и прогибающаяся под собственным весом q (рис. 4). 178 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 4 5 1 2 q 3 Рис. 4. Расчетная схема для определения напряженно-дефор- мированного состояния подрабатываемой толщи горных пород: 1 – ранее выработанное пространство; 2 – действующая лава; 3 – основание; 4 – земная поверхность; 5 – мульда Жесткость основания 3 зависит от наличия выработанных пространств и вынимаемой мощности m пластов. Позициями 1 и 2 указаны ранее выработанное пространство и то, которое формируется примыкающей действующей лавой. В результате отработки лав на земной поверхности 4 формируется мульда сдвижений 5. Распределение напряжений и деформаций в толстой плите мощностью 150 m можно найти с помощью трехмерных моделей. Однако, если заменить толстую 1 2 q 150m 179 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 плиту тонкой, можно подобрать такие эквивалентные величины жесткости основания и деформационных параметров плиты, при которых некоторые параметры напряженно-деформированного состояния будут соответствовать фактическим. Например, с помощью такой модели можно с высокой достоверностью прогнозировать параметры мульды сдвижений на земной поверхности, а также прогибы и напряжения в породных слоях в пределах высоты подрабатываемой толщи, равной 150 m. Задача о прогибе плиты-слоя над выработанными пространствами решается путем численного дифференцирования уравнения в частных производных методом конечных разностей. , D K D q yyx 2 x 4 4 22 4 4 4 ω −= ∂ ω∂ + ∂∂ ω∂ + ∂ ω∂ (1) где ω – прогибы точек плиты; q – нормальная к срединной плоскости и распределенная по поверхности плиты пригрузка; K – коэффициент жесткости основания, изменяющийся в зависимости от наличия или отсутствия выработанного пространства; D – цилиндрическая жесткость плиты; ( )2 3 112 µ− = EhD , (2) где h – мощность плиты; E – приведенный модуль деформации слоистого массива; Μ – коэффициент Пуассона вмещающих пород. В итоге получают опускания подработанной толщи на уровне верхней границы основной кровли пласта. Для учета механических свойств всех слоев, составляющих толщу, их величины приводились к срединной плоскости плиты- слоя по следующей зависимости: ∑ ∑= i ii p Ep E , (3) 180 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 где Ei – модуль деформации i-го слоя; pi – вес влияния i-го слоя, который находится по формуле i i i y hp = , (4) где hi – мощность i-го слоя; yi – расстояние от центра тяжести i-го слоя до срединной плоскости. Такая зависимость позволила учесть свойства всех слоев подработанной толщи. Указанная модель учитывает также релаксацию напряжений и ползучесть пород, слагающих подрабатываемый массив. Это является существенным преимуществом выбранной модели, поскольку динамические процессы проявлений сдвижения горных пород являются выраженными неравновесными физическими процессами и существенно зависят от времени. Сначала определяется начальное напряженное состояние выработанного пространства. Это позволит учесть неоднородность напряженного состояния до начала работы 13-й восточной лавы, а также даст возможность учесть влияние ранее выработанных пространств. На рис. 5 показано распределение вертикальной компоненты горного давления с учетом отработки 11-й и 12-й восточных лав до начала работы 13-й восточной лавы. Вокруг 12-й восточной лавы присутствует зона остаточного повышенного горного давления. Это означает, что 13-я восточная лава (ее верхняя часть) будет работать в критических условиях этой зоны, что будет представлять определенные сложности. Значения прочности вводились по скважинам, расположенным по всей площади выемочного пространства, что позволяет учесть ее неоднородность. В пределах выемочного столба приведенная прочность подработанной толщи меняется (в западной части больше, в восточной меньше), это окажет свое влияние на динамику обрушения кровли (рис. 6). 181 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 200 400 600 800 1000 1200 1400 1600 1800 2000 2200 2400 200 400 600 800 1000 1200 1400 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 Рис. 5. Распределение горного давления вокруг выработан- ного пространства с учетом отработки 11-й и 12-й восточных лав до начала работы 13-й восточной лавы пласта l1 шахты им. А.Ф. Засядько 182 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 5 10 15 20 25 30 35 40 45 5 10 15 20 25 0.145 0.155 0.165 0.175 0.185 0.195 0.205 0.215 0.225 0.235 0.245 0.255 0.265 0.275 0.285 Рис. 6. Распределение приведенной прочности подработан- ной толщи в пределах выемочного столба пласта l1 шахты им. А.Ф. Засядько В районе замка флексурной складки нет скважин, из-за чего прочность подработанной толщи определить невозможно. Поэтому на первом этапе расчеты велись без учета влияния замка флексурной складки (рис. 7, а, б, в). Единица соответствует нетронутому массиву либо оседанию подработанной толщи без разрушения. а) Скорость подвигания лавы 160 м/мес. 183 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 б) Скорость подвигания лавы 120 м/мес. Скорость подвигания лавы 80 м/мес. в) Рис. 7. Распределение высоты зоны полных обрушений ос- новной кровли 13-й восточной лавы пласта l1 шахты им. А.Ф. Засядько Чем больше впадина на распределении, тем на большую высоту разрушилась кровля в пределах 30 вынимаемых мощностей пласта. Из разрушенной кровли в этих зонах начнется выделение газа. 184 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 Расчет производился для трех скоростей подвигания лавы: 160 м/мес., 120 м/мес. и 80 м/мес. для того, чтобы учесть всевозможные неблагоприятные ситуации по мере подвигания лавы. Участки на распределении, указанные стрелкой, – это зоны зависания кровли, которые стремятся к границе ранее выработанного пространства. В зонах зависания газовыделение задерживается, но происходит позже, при отходе лавы. При скорости подвигания лавы 160 м/мес. шаг глобальной посадки около 150 м, всего 13 глобальных зависаний, при скорости подвигания лавы 120 м/мес. шаг посадки уменьшается, количество зависаний кровли увеличивается, около 16. При скорости подвигания лавы 80 м/мес. глобальных зависаний нет, общее количество зависаний приближается к 20, газовыделение будет происходить стабильно. Для учета влияния замка флексурной складки принимаем в этой зоне приведенную прочность подработанной толщи в 1,3 раза меньше (рис. 8, а, б, в). Скорость подвигания лавы 160 м/мес. а) 185 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 б) Скорость подвигания лавы 120 м/мес. Скорость подвигания лавы 80 м/мес. в) Рис. 8. Распределение высоты зоны полных обрушений ос- новной кровли 13-й восточной лавы пласта l1 шахты им. А.Ф. Засядько с учетом флексурной складки 186 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 В районе складки толща максимально разбита трещинами. После прохода лавой зоны складки общая картина зависаний кардинально меняется. Шаг зависаний неравномерный, высота обрушений увеличивается, зависаний становится меньше. При скорости подвигания лавы 120 м/мес. зависания и обрушения кровли будут происходить равномернее, чем при скорости подвигания лавы 160 м/мес. и 80 м/мес. Глобальные зависания могут возникнуть в начале отработки лавы (160 м/мес.) и в конце (80 м/мес.). ВЫВОДЫ Скорость подвигания лавы существенно влияет на характер динамики сдвижений подрабатываемой толщи: на длину зависаний и их высоту. При больших скоростях подвигания длина зависаний примерно 150 м, при малых – около 50 м. Наблюдается общая неравномерность зависаний по всей длине лавы. При уменьшении скорости подвигания уменьшается шаг зависаний. При малых скоростях появляются локальные зависания в некоторых участках лавы. Динамика газовыделения определяется скоростью подвигания лавы, чем меньше скорость, тем равномернее будет происходить газовыделение. СПИСОК ССЫЛОК 1. Земисев В.Н. Расчеты деформаций горного массива. – М.: Недра, 1973. – 144 с. 2. Зборщик М.П., Грищенков Н.Н. Компьютерное моделирование геомеханических параметров обработки свит угольных пластов // Днепропетровск: НГА Украины. – 1999. – Т. 2, №7. - С. 3 – 7. 3. Геомеханические основы предотвращения провалов земной поверхности при ликвидации шахт / Звягильский Е.Л., Минаев А.А., Назимко В.В. и др. - Донецк: ООО "Норд Компьютер", 2001. – 332 c. 187 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 4. Инструкция по безопасному ведению горных работ на пластах, опасных по внезапным выбросам угля, пород и газа. - М: Минуглепром СССР, 1989. - 192 с. 5. Муллер Р.А. Определение мульды сдвижения и скорости деформаций земной поверхности при движущемся забое // Сдвижение горных пород. – Л.: Изд. ВНИМИ, 1975. – Сб. 96. – С. 13-27. 6. Медянцев А.Н., Черняев В.И. Сдвижение и деформации горных пород в подрабатываемой толще // Уголь Украины. – 1961. – № 11. – С. 25–27. 7. Иофис М.А., Фастов Г.А. Характер развития деформаций в полумульде над движущимся забоем // Горное давление, сдвижение горных пород и методика маркшейдерских работ / ВНИМИ. – 1965. – Сб. 55. – С. 143-149. 8. Земисев В.Н. Современные проблемы сдвижения горных пород. // Уголь. - 1999.– № 10.– С. 23. 9. Костенич В.С., Зелепукин А.М. Результаты наблюдений за сдвижением и деформацией поверхности при большой скорости подвигания забоя и прерывной рабочей неделе // Сдвижения и деформации массива при разработке месторождений с учетом структуры и механических свойств горных пород / ВНИМИ. – 1966. – Сб. 58. – С. 137-146. 10. Петрук Е.Г. Исследование деформаций земной поверхности в мульде сдвижения по времени // Изв. вузов. горный журнал. – 1969. – № 1. – С. 40-43. 11. Мякенький В.И. Сдвижение и дегазация пород и угольных пластов при очистных работах. – Киев: Наукова думка, 1975. – 100 с. 12. Ковальски А. Деформация поверхности над быстро двигающимся фронтом горных работ // Proceedings IX Congress ISM, Praha 18.4-22.4. 1994. – С. 320-329. 13. Гавриленко Ю.Н., Папазов Н.М., Морозова Т.В. Динамика оседаний земной поверхности при большой глубине разработки и высокой скорости подвигания забоя // Проблеми гірського тиску.-Донецьк, ДонДТУ, 2000. - № 4. - С. 108-119. 14. Сдвижение массива горных пород и влияние скорости подвигания очистного забоя на выделение метана / Четверик 188 Наукові праці УкрНДМІ НАН України, № 1, 2007 Transactions of UkrNDMI NAN Ukraine, № 1, 2007 М.С., Андрощук Е.В., Бубнова Е.А., Гаврюк Г.Ф. // Сб. научн. тр. НГУ. – Днепропетровск:- 2003. – вып. 16. – С. 42-53. 15. Кузяра С.В., Назимко В.В. Геомеханические основы модели обрушения пород в окрестности движущегося очистного забоя // Проблеми гірського тиску. - Донецьк, ДонНТУ. – 2003 - №10. - С. 158-171. 189