CC BY
7
УДК 622.83 DOI 10.46689/2218-5194-2022-1-1-303-316
ПЕРЕРАСПРЕДЕЛЕНИЕ МАКСИМАЛЬНЫХ ВЕРТИКАЛЬНЫХ ДЕФОРМАЦИЙ ГОРНОГО МАССИВА В ПРОЦЕССЕ ЕГО ПОДРАБОТКИ
С.Б. Кулибаба, Е.В. Федоров
Приведены результаты инструментальных натурных наблюдений за процессом сдвижения массива горных пород, полученные при подработке вертикального шахтного ствола очистными работами угольного пласта. На основе их анализа установлен характер перераспределения максимальных деформаций вертикального растяжения горного массива на участке между пластом и земной поверхностью в зависимости от местоположения движущегося очистного забоя.
Ключевые слова: угольный пласт, подработка горного массива, перераспределение деформаций, установление зависимости.
На начальных этапах изучения сдвижения горных пород основное внимание уделялось величинам деформаций земной поверхности и распределению их в границах мульды сдвижения. При существовавших тогда малых скоростях подвигания очистных забоев и глубинах разработки пластовых месторождений вопросы динамики сдвижения породной толщи были мало актуальны, и критерием степени влияния горных разработок на охраняемые объекты являлись в основном конечные значения деформаций горных пород, земной поверхности и подрабатываемых объектов.
Увеличение глубины и скорости разработки угольных пластов вызвало необходимость в смещении акцента в изучении сдвижения с земной поверхности на горный массив [1 - 3], а также в расширении области и направлений исследований вследствие ряда важных факторов. Во-первых, на фоне общего уменьшения деформаций земной поверхности, более рельефно обозначилась разница между их конечными значениями, образующимися после окончания процесса сдвижения, и временными, возникающими в динамической полумульде, перемещающейся вместе с движущимся очистным забоем. Во-вторых, в область влияния очистных работ стало попадать все большее количество действующих горных выработок и подземных сооружений, расположенных на вышележащих горизонтах, а также вертикальные выработки (стволы, технические скважины), для охраны которых знания закономерностей развития сдвижения только лишь земной поверхности оказались недостаточными, а динамика развития этого процесса в породной толще - весьма актуальной. В-третьих, трансформация напряженно-деформированного состояния, подрабатываемого угле-породного массива в новых условиях повлияла на изменение его диффузионных свойств и характеристик массопереноса, что, в свою очередь, необходимо учитывать при оценке процессов метановыделения из
подрабатываемых угольных пластов и вмещающих пород, а также при расчетах параметров дегазации и вентиляции на газовых шахтах. Принимая во внимание тот факт, что метановыделение с площади подработанной территории происходит в течение весьма длительного периода, в т. ч., и после ликвидации шахты, и зависит от газовой проницаемости подработанных пород [4-7], последний из перечисленных выше факторов заслуживает особого внимания.
Таким образом, исследования, направленные на повышение точности прогноза развития деформирования породного массива в процессе ведения очистных работ, является актуальными, поскольку позволяют обеспечивать эффективные меры охраны подрабатываемых объектов, а также обоснованно оценивать параметры зон повышенной проницаемости и интенсивность процессов десорбции метана из подрабатываемых сближенных угольных пластов.
Известно, что в динамической мульде сдвижения величины деформаций земной поверхности всегда меньше, чем в окончательно сформировавшейся, а граничные углы впереди движущегося забоя - круче их окончательных значений, сформировавшихся после затухания процесса сдвижения. Этот установленный экспериментально факт [8-11] объясняется тем, что над движущимся забоем полноценная мульда на земной поверхности не успевает сформироваться, вследствие чего момент достижения ее деформациями критических значений, характеризующих границу мульды впереди забоя, запаздывает. При этом динамические деформации подрабатываемых горных пород при определенных условиях превышают конечные значения, а скорость распространения процесса сдвижения в массиве на разных стадиях подработки различна [12].
Одним из эффективных методов исследования процесса деформирования породного массива над движущимся очистным забоем является анализ результатов экспериментов, проводимых в натурных условиях при подработке наблюдательных станций, оборудованных в вертикальных горных выработках, в частности, в шахтных стволах. Следует отметить, что вследствие большой значимости указанных горных выработок, осуществляющих все основные функции связи подземных горизонтов с земной поверхностью, подобные эксперименты проводились крайне редко, и лишь в единичных случаях их результаты оказались пригодными для всестороннего анализа. Рассмотрим результаты эксперимента, проводимого в условиях шахты "Глубокая" производственного объединения "Донецкуголь", и описанного нами в ряде публикаций [12 и др.], где в 1988-89 гг. была осуществлена подработка вентиляционного ствола № 1 под контролем инструментальных наблюдений.
Подработка проводилась двумя спаренными лавами пласта И6 шахты "Глубокая" в следующих условиях: средняя глубина разработки 541 м, вынимаемая мощность пласта 1,15 м, угол падения 5-6°, общая длина двух
подрабатывающих лав 400 м, управление кровлей - полное обрушение, месячное подвигание очистных забоев в пределах целика - до 20 м. Глубина подрабатываемого ствола - 426 м. Инструментальные наблюдения за процессом сдвижения горного массива проводились по 17-ти реперам (Яр 1 -Яр 17), заложенным на разных уровнях в стенке ствола. В результате инструментальных наблюдений определялись оседания реперов на различные моменты подработки (рис. 1). Всего за период с 14.04.1988 г. по 24.04.1991 г. на наблюдательной станции было проведено 12 серий инструментальных наблюдений.
0 0,2 0,4 0,6 Л,м
(5) - номер серии наблюдений
Рис. 1. Вертикальный разрез по простиранию пласта Не и графики оседания реперов ствола на различные даты наблюдений
Из графиков нарастающих оседаний реперов, зафиксированных на отдельные этапы наблюдений, можно видеть, что в начальный период подработки их сдвижение происходило синхронно и равномерно - величина оседаний практически не изменялась на всем интервале глубин. В последующие периоды оседания нижних реперов стали опережать оседания верхних, что явилось следствием образования в породные толще вертикальных деформаций (далее - деформаций) растяжения. Так, разница в оседаниях верхнего и нижнего реперов в 7-й серии составила более 200
мм. В дальнейшем после удаления груди очистного забоя от оси вертикального ствола, эта разница уменьшалась, и в последних сериях составила около 100 мм, хотя качественно картина не изменилась - после затухания процесса сдвижения нижние реперы осели на большую величину, чем верхние.
Рассмотрим изменение относительных деформаций в анализируемого участка массива горных пород, испытанные им при подработке, после пересечения лавой вертикальной оси подрабатываемого ствола и определяемых по формуле
в = (1)
ЛЯ
где Лц - разность оседаний двух точек массива, м; ЛН - вертикальное расстояние между данными точками, равное разности глубин Н залегания этих точек, м.
Выше было показано, что оседания точек массива распределяются неравномерно, как по глубине ствола, так и по времени протекания процесса сдвижения. На рис. 2 а показаны графики развития во времени относительных деформаций растяжения трех укрупненных участков массива, каждый из которых составляет приблизительно треть мощности наблюдаемой толщи, - верхнего, в интервалах глубин от земной поверхности 35...150 м, среднего (150...275 м) и нижнего (275...415 м), а также средних деформаций по стволу. Здесь t - время в месяцах, прошедшее после пересечения лавой вертикальной оси ствола (^ = 0 - момент пересечения).
Примечание. Здесь и далее: деформации растяжения имеют положительный знак, что является традиционным в науке о сдвижении горных пород и земной поверхности [13, 14].
Графики показывают, что нижние слои толщи при подработке подверглись существенно большим деформациям, чем верхние, а при удалении вверх от разрабатываемого пласта их величины интенсивно затухали. Так, если максимальные деформаций растяжения участка массива, расположенного в нижней части ствола, превысили 1,0 %о, то в средней части они составили 0,45 %о, а в верхней - всего 0,17 %о. Для сравнения на графике пунктиром показана кривая изменения средней относительной деформации всего участка толщи, расположенной в пределах наблюдательной станции. Кроме того, очевидно несовпадение во времени момента образования максимума деформаций участков массива, расположенных на разных глубинах. Если в нижней трети ствола этот максимум приходился на четвертый месяц после момента пересечения очистным забоем его оси, то в верхних участках он перемещался на более поздние периоды, и для средней части ствола составил 7, а для верхней - 10 - 12 месяцев. На графике траектория смещения максимума деформаций во времени показана точечной пунктирной линией.
Рис. 2. График развития относительных (е) и приведенных (Е) деформаций различных участков массива в процессе подработки
Как упоминалось выше, опережение оседаний нижних участков массива по отношению к верхним, вызывало деформации растяжения толщи. Максимум этих деформаций по мере продвижения лавы перемещался от нижних участков массива к верхним, а значения их уменьшались. После завершения процесса сдвижения весь исследуемый участок толщи горных пород оказался в области вертикального растяжения относительно своего начального состояния, но величины относительных деформаций в пределах этого участка распределились неравномерно, с уменьшением их значений по направлению снизу верх от подрабатывающего пласта к земной поверхности.
Следует подчеркнуть, что графики на рис. 2, а характеризуют развитие оседаний массива лишь в конкретных условиях описанного эксперимента. В то же время известно, что, например, величина оседания той или иной точки подрабатываемого массива при движении лавы зависит отнюдь не от времени, прошедшего после начала процесса сдвижения мас-
сива в области данной точки, а от расположения ее относительно груди очистного забоя. Теоретически процесс сдвижения земной поверхности в точке Яр 1, расположенной в устье ствола (см. рис. 1), начинается в момент нахождения лавы в точке А пласта, определяемой граничным углом 50, а заканчивается при нахождении лавы в точке В, определяемой углом полных сдвижений у3. Очевидно, что при различных скоростях подвигания очистного забоя в интервале АВ время подработки точки Яр 1 будет разным. Поэтому логично исследовать характер развития процесса сдвижения в зависимости не от параметра времени, а от местоположения лавы относительно подрабатываемых точек, например, от величины отхода лавы Ь от проекции вертикальной оси ствола на подрабатывающий пласт, обозначенной на разрезе точкой О.
Поэтому для унификации условий логично рассматривать относительные значения также и других исследуемых параметров, а именно: относительной высоты И расположения экспериментальной точки над разрабатываемым пластом, относительного отхода I груди очистного забоя от вертикальной проекции оси ствола на пласт и относительного оседания q массива в экспериментальной точке:
н = Яп; (2)
Я 0
I=Ь; (3)
Ь0
Ч = -, (4)
т
где Яп - вертикальное расстояние от разрабатываемого пласта в главном сечении мульды сдвижения по простиранию до рассматриваемой точки массива, м; Яо - средняя глубина залегания пласта, м; Ь - величина отхода лавы от проекции рассматриваемой вертикальной оси на подрабатывающий пласт, м; ц - оседание массива в рассматриваемой точке массива, м; т - вынимаемая мощность пласта, м; Ьо - отход лавы ОВ от вертикальной оси ствола, необходимый для полной подработки верхней точки исследуемого участка массива (точка Яр 1 на рис. 1), м.
Ьо = Яо ■ , (5)
где - угол полных сдвижений по простиранию пласта, градус.
С целью получения достоверных результатов последующего анализа нами были исключены из выборки нерепрезентативные исходные данные, включающие либо грубые ошибки, возникшие при проведении эксперимента, либо данные, полученные в заведомо неоднородных условиях.
Так, из 17 реперов наблюдательной станции при последующем анализе нами были исключены Яр 1 и Яр 17 по следующим причинам. Оседание Яр 1 практически не отличалось от оседания Яр 2, поскольку, во-
первых, интервал толщи, расположенной между ними, представлял собой приповерхностный слой, включающий наносы, где характер деформирования существенно отличается от остального подрабатываемого массива [15], а во-вторых, эти два репера были заложены в крепь устья ствола, отличающейся увеличенной толщиной и наличием стального армирования, вследствие чего их оседание происходило консолидировано. Характер оседания нижнего репера Яр 17 также существенно отличался от основной массы реперов, поскольку он был заложен в районе сопряжения ствола с выработками рабочего горизонта, где горный массив, как известно, подвергается значительно большим деформациям по сравнению с протяженным его участком. Кроме того, из 12-ти серий инструментальных наблюдений нами были исключены результаты серий 10, 11 и 12, поскольку после 9-й серии одна из двух спаренных лав, подрабатывающих ствол, была остановлена, что, очевидно, изменило характер дальнейшего развития процесса сдвижения. Таким образом, нами исследовался массив, состоящий из 136 исходных данных (матрица 15x9), при моделировании свойств которого использовался метод регрессионного анализа [16, 17].
На рис. 3 показаны графики распределения параметра д, рассчитанного по формуле (4), в зависимости от относительной высоты И в экспериментальных точках на моменты проведения серий 2 - 9 наблюдений относительно первой серии наблюдений (рис. 3, а), и графики распределения этого же параметра в зависимости от относительного отхода I для некоторых реперов (рис. 3, б).
Для каждой серии наблюдений фактические значения д в зависимости от И, показанные точками на графике, были аппроксимированы с высокой степенью точности функцией
q - а0 + , (6)
п + а2
где а - эмпирические коэффициенты.
Найдем приведенные деформации Е массива, определяемые первой производной параметра д (6) по И:
Е = § = -£<«> + =--^ (7)
М М п + а2 (^ + а2)2
По структуре величина Е аналогична рассмотренной выше относительной деформации массива е, определяемой формулой (1), что видно из выражения (8), однако, в отличие от нее, Е показывает изменение величины д (вместо п, м) при изменении И (вместо Н, м) на ¡-м интервале массива.
^ 1
^ т ) _ АЩ Н0_ Но Е ~Т7~ / лет л " ТГГ---- 8/-. (8)
АП
АНп Н 0
АНп т т
а
б
Рис. 3. Графики распределения параметра # в зависимости от относительной высоты Н расположения экспериментальной точки над пластом и относительного отхода I лавы
На рис. 2, б показаны графики изменения Е (7) для реперов, расположенных на различной относительной высоте И над пластом, в зависимости от относительного отхода лавы /, описываемые формулой
ВД + Ь2)
Е = Ьо +
N+Ь4
(9)
где Ь - эмпирические коэффициенты.
Из графиков видно, что с увеличением /, т. е., с удалением лавы от рассматриваемой вертикальной оси, точки максимумов величины Е для различных значений И смещаются (пунктиром показана траектория этого смещения). Так, если этот максимум для значения И=0,273 (Яр 16), образовался на относительном удалении /, равном 0,13 - 0,15, то при И=0,487 (Яр 10) - на удалении /, равном 0,29 - 0,31. При этом величины деформаций в активной стадии процесса значительно (в 1,3 - 2,4 раза) превышают конечные свои значения при его затухании.
Как можно видеть, характер этих графиков повторяет показанную выше зависимость развития относительных деформаций укрупненных участков массива в процессе подработки (см. рис. 2, а).
Более строгое определение значений /, при которых приведенные деформации Е достигают своих максимумов для точек массива, расположенных на разных относительных высотах И над пластом, можно осуществить путем приравнивания к нулю первой производной Е (9) по /:
dE Ь + Ьз Ь + Ь2 ) ЩЬц/ + Ь2 )
сИ
(Ьх/ + Ь2) еЬз/+Ь4
0.
(10)
Сопоставляя полученные в результате решения уравнения (10) значения / для каждой экспериментальной точки с ее относительной высотой Итах, на которой в этот момент будет находиться максимум приведенной
деформации Е, получим кривую, показанную на рис. 4, а, которая описывается формулой
/ Л
*тах - ±
С0
1п
С1
I - С2
-1
С 3
(11)
где с0, с1, с2, с3, - безразмерные эмпирические коэффициенты, равные для рассматриваемых условий соответственно -7,82; 0,4...0,5; -0,03...-0,04 и 2,6...2,7.
Для перехода от относительных величин в выражении (11) к абсолютным воспользуемся формулами (2), (3) и (5), в результате чего оно примет вид
Н
п тах
С0
1п
С1
ь
-1
0
- С2
с 3
(12)
где Нптах - вертикальное расстояние от разрабатываемого пласта до точки
массива с максимальной деформацией растяжения в, м.
Подставляя в формулу (12) соответствующие исходные параметры, можно рассчитать в абсолютных величинах перемещения точки максимальных деформаций растяжения породной толщи вдоль вертикальной оси при удалении от нее лавы. Так, на рис. 4, б показаны графики функции Нптах - /(Ь) для условий Донбасса (описанная выше шахта "Глубокая" объединения "Донецкуголь") и Кузбасса (лава 2462 пласта Болдыревский шахты им. С. М. Кирова, АО "СУЭК-Кузбасс").
Анализируя представленные выше графики, можно сделать следующие обобщающие выводы о характере перераспределения деформаций массива горных пород в процессе его подработки.
После пересечения лавой рассматриваемой вертикальной оси (т.е., при I — 0) в подрабатываемом массиве начался процесс развития деформаций растяжения, интенсивность которого уменьшалась при удалении от разрабатываемого пласта к земной поверхности (см. рис. 2). При возрастании параметра I до определенного значения эти деформации в точках вертикальной оси, расположенных на различном относительном удалении И от пласта, достигали максимальных величин. Причем своих максимумов деформации на разных удалениях от пласта достигали не одновременно -сначала они локализировались на более глубоких участках массива, постепенно перемещаясь вверх к земной поверхности в процессе дальнейшего отхода лавы от рассматриваемой вертикальной оси. При этом величины деформации массива снижались - в активной стадии процесса они значительно (в 1,3 - 2,4 раза) превышали свои конечные значения при его затухании.
тах, М
Донбасс 1 1
/
\
V Кузбасс
50
100
150
м
а
б
Рис. 4. Графики перемещения точки максимальных деформаций растяжения массива вдоль вертикальной оси при отходе от нее лавы: а - в относительных величинах (Нтах); б - в абсолютных (Нп тах)
Перемещение точки максимума деформаций при приближении к земной поверхности происходило с ускорением, что следует из формы графиков на рис. 4. Неравномерность скорости распространения процесса сдвижения массива явилась следствием условий зависания породных слоев над выработанным пространством - с удалением вверх от выемочного участка область сдвижения массива расширялась, в связи с чем в верхних породных слоях опорные части оказывались отнесенными в неподрабо-танную толщу на более значительные расстояния, чем в нижних, а, следовательно, зависающие породные консоли имели большую длину и обру-шались активней [12, 18, 19]. Интенсивное возрастание значений функции (11) при приближении к своей вертикальной асимптоте, определяемой выражением I = с + С2 (на рис. 4, а она показана штрих-пунктиром), говорит о том, что максимумы деформаций растяжения в верхних слоях подрабатываемой породной толщи проявлялись практически одновременно.
С удалением вверх от выработанного пространства к земной поверхности степень деформирования горных пород уменьшается на любом этапе подработки, что следует из графика распределения приведенных деформаций Е массива в системе координат I - И (рис. 5). Однако поскольку в отличие от описанных выше графиков перемещения точки максимальных деформаций при отходе лавы (см. рис. 4), значения самих деформаций Е зависят также и от других параметров подработки, а именно: от размера очистной выработки вкрест простирания, угла падения пласта и др., которые в данной статье не рассматриваются, значения Е на рис. 5 соответствуют лишь конкретным условиям рассматриваемого эксперимента.
0,2
0 0,1 0,2 0,3 0,4 0,5 0,6 /
Рис. 5. Распределение приведенных деформаций E горного массива
при его подработке
Таким образом, в результате проведенных исследований выполнен анализ изменения деформационного состояния горных пород в процессе их подработки и установлен характер перераспределения максимальных вертикальных деформаций растяжения горного массива на участке между пластом и земной поверхностью в зависимости от местоположения движущегося очистного забоя.
Список литературы
1. Коряков А. Е., Копылов А. Б., Савин И. И. Применение метода конечных элементов при моделировании сдвижения элементов горного массива // Известия Тульского государственного университета. Науки о Земле. 2017. Вып. 4. С. 344-354.
2. Kratzsch H. Bergschadenkunde. Bochum: Deutscher MarkscheiderVereine. V. 1997. 844 s.
3. Рожко М. Д. Расчет ожидаемых максимальных оседаний земной поверхности в условиях больших глубин разработки пластов // Науковi пращ УкрНДМ1 НАН Украши. 2014. № 14. С. 241-246.
4. Метановыделение на земную поверхность для территорий горных отводов ликвидированных шахт Кузбасса / Н. М. Качурин, Д. Н. Шку-
ратский, Л. Л. Рыбак, Р. В. Сидоров // Известия Тульского государственного университета. Науки о Земле. 2015. Вып. 2. С. 42-48.
5. Оценка газообмена подработанных территорий угольных бассейнов России с атмосферой / Н. М. Качурин, Т. В. Корчагина, А. Н. Качу-рин, Р. В. Сидоров // Известия Тульского государственного университета. Науки о Земле. 2021. Вып. 1. С. 290-302.
6. Анализ экологического мониторинга в районах ликвидации угольных шахт России / Т. В. Корчагина, Н. П. Иватанова, И. А. Басова, Ю. А. Воронкова // Известия Тульского государственного университета. Науки о Земле. 2017. Вып. 4. С. 39-47.
7. Исследование эмиссии угольного метана на поверхность из ликвидируемых шахт / Т. В. Корчагина, Н. В. Ефимова, А. Б. Жабин, С. А. Ишутина // Известия Тульского государственного университета. Науки о Земле. 2017. Вып. 4. С. 48-60.
8. Батугин С. А. Сдвижения и деформации земной поверхности и горных пород над движущимся забоем // Труды по вопросам горного давления, сдвижения горных пород и методики маркшейдерских работ: сб. науч. тр. Л.: ВНИМИ, 1962. Сб.47. С. 159 - 199.
9. Батугин С.А., Бежецкий А. Е. Сдвижение земной поверхности при большой скорости подвигания очистного забоя // Горное давление, сдвижение горных пород и методика маркшейдерских работ: сб. науч. тр. Л.: ВНИМИ. 1964. С. 252 - 265.
10. Охрана подрабатываемых подготовительных выработок / Н.П. Бажин, О.И. Мельников, В.С. Пиховкин, В.В. Райский. М.: Недра, 1978. 253 с.
11. Костенич В. С., Зелепукин А. М. Результаты наблюдений за сдвижением и деформацией земной поверхности при большой скорости подвигания очистного забоя и прерывной рабочей неделе // Горное давление, сдвижение горных пород и методика маркшейдерских работ. Л.: ВНИМИ, 1966. С. 137 - 146.
12. Кулибаба С. Б. Характер развития процесса сдвижения горного массива над движущимся очистным забоем // Сб. науч. тр. Междунар. науч.-практич. конф. «50 лет Российской научной школе комплексного освоения недр Земли». М.: ИПКОН РАН, 2017. С. 7-72.
13. Земисев В. Н. Расчеты деформаций горного массива. М.: Недра, 1973. 145 с.
14. Правила охраны сооружений и природных объектов от вредного влияния подземных горных разработок на угольных месторождениях. С.Петербург: ВНИМИ, 1998. 291 с.
15. Кулибаба С. Б. О характере сдвижения верхних слоев подрабатываемого массива горных пород // ГИАБ (научно-технический журнал). 2007. № 3. С. 84-89.
16. Дрейпер Н., Смит Г. Прикладной регрессионный анализ. М.: Статистика, 1973. 392 с.
17. Смирнов Н.В., Дунин-Барковский И.В. Курс теории вероятностей и математической статистики для технических приложений. М.: Наука, 1965. 512 с.
18. Шалагинов Н.Ф. Деформации различных слоев подработанной толщи и зависимость этих деформаций от расстояния до пласта, мощности пласта и ширины выработки // Сб. наут тр. по вопросам исследования горного давления и сдвижения горных пород. Л.: ВНИМИ. 1961. Вып. 43. С. 133-137.
19. Кратч Г. Сдвижение горных пород и защита подрабатываемых сооружений. М.: Недра, 1978. 494 с.
Кулибаба Сергей Борисович, д-р техн. наук, проф., вед. науч. сотр., kulihaha saipkonran. ru , Россия, Москва, Институт проблем комплексного освоения недр им. акад. Н. В. Мельникова РАН,
Федоров Евгений Вячеславович, канд. техн. наук, зав. отделом, fedo-rov eaipkonran.ru , Россия, Москва, Институт проблем комплексного освоения недр им. акад. Н. В. Мельникова РАН
REDISTRIBUTION OF THE MAXIMUM VERTICAL DEFORMATIONS OF THE ROCK MASSIF IN THE PROCESS OF ITS UNDERMINING
S. B. Kulihaha, E. V. Fedorov
The results of instrumental field observations of the process of the rock massif movement, obtained during the undermining of a vertical mine shaft hy the coal-seam face, are presented. Based on their analysis, the nature of the redistrihution of the maximum vertical tensile deformations of the rock mass in the area hetween the seam and the earth's surface, depending on the location of the moving working face, has heen estahlished.
Key words: coal seam, the undermining of rock massif, redistrihution of deformations, estahlishment of dependence.
Kulihaha Sergey Borisovich, doctor of technical sciences, professor, leading researcher, kulihaha sa ipkonran.ru , Russia, Moscow, Institute of Comprehensive Exploitation of Mineral Resources named after N. V. Melnikov Russian Academy of Sciences,
Fedorov Evgeniy Vyacheslavovich, candidate of technical sciences, department head, fedorov ea ipkonran.ru , Russia, Moscow, Institute of Comprehensive Exploitation of Mineral Resources named after N. V. Melnikov Russian Academy of Sciences
Reference
1. Koryakov A. E., Kopylov A. B., Savin I. I. Application of the finite element method in modeling the movement of elements of a mountain massif // Proceedings of Tula State University. Earth sciences. 2017. Issue 4. pp. 344-354.
2. Kratzsch H. Bergschadenkunde. Bochum: Deutscher Markscheider-Vereine. V. 1997. 844 s.
3. Rozhko M. D. Calculation of the expected maximum subsidence of the earth's surface in conditions of large depths of reservoir development // Naukovi prati UkrNDMI NAS of Ukraine. 2014. No. 14. pp. 241-246.
4. Methane release to the Earth's surface for the territories of mining branches of liquidated Kuzbass mines / N. M. Kachurin, D. N. Shkuratsky, L. L. Rybak, R. V. Sidorov // Izvestiya Tula State University. Earth sciences. 2015. Issue. 2. pp. 42-48.
5. Assessment of gas exchange of the under-worked territories of coal basins of Russia with the atmosphere / N. M. Kachurin, T. V. Korchagina, A. N. Kachurin, R. V. Sidorov // Izvestiya Tula State University. Earth sciences. 2021. Issue 1. pp. 290-302.
6. Analysis of environmental monitoring in the areas of liquidation of coal mines in Russia / T. V. Korchagina, N. P. Ivatanova, I. A. Basova, Yu. A. Voronkova // Proceedings of Tula State University. Earth sciences. 2017. Issue 4. pp. 39-47.
7. Investigation of coal methane emission to the surface from liquidated mines / T. V. Korchagina, N. V. Efimova, A. B. Zhabin, S. A. Ishutina // Proceedings of Tula State University. Earth sciences. 2017. Issue 4. pp. 48-60.
8. Batugin S. A. Displacement and deformation of the earth's surface and rocks above the moving bottom // Works on issues of rock pressure, the displacement of rocks and methods of surveys: Sat. nauch. Tr. Leningrad: VNIMI, 1962. Sat.47. S. 159 - 199.
9. Batugin S. A., bezhetskiy A. E. Displacement of the earth's surface at a high speed podvigina stope // confining pressure, the rock movement and the method of surveying work. Sat. nauch. Tr. Leningrad: VNIMI. 1964. Pp. 252 - 265.
10. Protection of part-time preparatory workings / N.P. Bazhin, O.I. Melnikov, V.S. Pihovkin, V.V. Raisky. M.: Nedra, 1978. 253 p.
11. Kostenich V. S., Zelepukin A.M. The results of observations of the displacement and deformation of the Earth's surface at a high rate of movement of the treatment face and an intermittent working week // Sb. nauch. tr. Mountain pressure, displacement of rocks and methods of surveying. L.: VNIMI, 1966. pp. 137 - 146.
12. Kulibaba S. B. nature of the development process of rock massif above the moving purification slaughter // Proc. nauch. Tr. Intern. nauch.-practical. Conf. 50 years of the Russian scientific school of comprehensive exploitation of mineral resources of the earth. M.: - IPKON RAS. 2017. Pp. 7-72.
13. Semisi V. N. The calculation of rock mass deformation. Moscow: Nedra, 1973.
145 p.
14. Rules for the protection of structures and natural objects from the harmful effects of underground mining in coal deposits. St. Petersburg, VNIMI. 1998. 291 p.
15. Kulibaba S. B. On the nature of the displacement of the upper layers of the mined rock mass // GIAB (scientific and technical journal). 2007. No. 3. pp. 84-89.
16. Draper N., Smith G. Applied regression analysis. M.: Statistics, 1973. 392 p.
17. Smirnov N.V., Dunin-Barkovsky I.V. Course of probability theory and mathematical statistics for technical applications. Moscow: Nauka, 1965. 512 p.
18. Shalaginov N. F. Deformation of the various layers of undermined strata and the dependence of these deformations of the distance to the reservoir, bed thickness and width of the production // Proc. naut Tr. the research questions of the mining pressure and strata movement. Leningrad: VNIMI. 1961. Vol. 43. S. 133-137.
19. Crutch G. rock Movement and protection of the mining structures. M.: Nedra, 1978. 494c.
