CC BY
21
DOI: 10.26730/1999-4125-2017-6-76-83 УДК622.1:622.834
ОСНОВНЫЕ ПОЛОЖЕНИЯ МЕТОДИКИ ПРОГНОЗА ОЖИДАЕМЫХ ДИНАМИЧЕСКИХ СДВИЖЕНИЙ ТОЧЕК ЗЕМНОЙ ПОВЕРХНОСТИ ПО ПРОСТИРАНИЮ ЛАВЫ ПРИ ОТРАБОТКЕ ПОЛОГИХ И НАКЛОННЫХ УГОЛЬНЫХ ПЛАСТОВ КУЗБАССА ВЫСОКОМЕХАНИЗИРОВАННЫМИ ЗАБОЯМИ
THE PRINCIPLES FOR CALCULATING THE ANTICIPATED DYNAMIC DISPLACEMENTS OF THE POINTS OF THE UNDERMINED EARTH'S SURFACE ALONG THE STRIKE OF THE LONGWALL DURING MINING OF THE SLOPING AND INCLINED COAL SEAMS OF KUZBASS BY HIGHLY MECHANIZED FACES
Свирко Сергей Владимирович1,
аспирант, e-mail: svirko@ngs.ru Sergey V. Svirko1, Postgraduate Ренев Алексей Агафангелович1, докт. техн.наук, зав.каф., профессор, e-mail: raa@kuzstu.ru Aleksey A. Renev1, D. Sc. in Engineering, Professor, Head of the Department
Быкадоров Алексей Иванович2, канд.техн. наук, ген. директор, e-mail: aibykadorov@mail.ru Aleksey I. Bykadorov2, С. Sc. in Engineering, Director General
кузбасский государственный технический университет имени Т.Ф. Горбачева, 650000, Россия, г. Кемерово, ул. Весенняя, 28
1T.F. Gorbachev Kuzbass State Technical University, 28 streetVesennyaya, Kemerovo, 650000, Russian Federation
2 Сибирский институт геотехнических исследований, 653000, Россия, Кемеровская область, г. Прокопьевск, ул. К.Либкнехта, 4
2 Siberian Institute of Geotechnical Investigations, 4 street Karl Liebknecht, Prokopyevsk, 653000, Russian Federation
Аннотация: Системно излагаются основы расчета ожидаемых динамических сдвиженийточек земной поверхности по простиранию лавы при отработке пологих и наклонных угольных пластов Кузбасса высокомеханизированными забоями. Используемые аналитические зависимости установлены по данным частотных инструментальных наблюдений на профильных линиях.
Приведен пример расчетов из научно-практического опыта Сибирского института геотехнических исследований в области решения геотехнических задач при разработке угольных месторождений.
Abstract: The principles for calculating the anticipated dynamic displacements of the points of the undermined earth's surface along the longwall strike during mining of the sloping and inclined coal seams of Kuzbass by highly mechanized faces are systematically presented. The applied analytical dependencies are established based on the data of instrumental frequency observations on the profile lines.
An example of calculations from the scientific and practical experience of the Siberian Institute of Geotech-nical Research in the field of solving geotechnical problems in the development of coal deposits is given.
Ключевые слова: Земная поверхность, подработка, динамическая мульда, оседание, горизонтальное сдвижение, скорость подвигания забоя, расчет.
Key words: Daylight surface, undermining, dynamic trough, subsidence, horizontal displacements, rate of face advance, calculation.
При отработке угольных пластов системами разработки с обрушением кровли на земной поверхности образуется мульда сдвижения- участок земной поверхности, подвергшийся сдвижению под влиянием очистной выработки. Размеры
мульды, величины ожидаемых сдвижений и деформаций поверхности определяются по основному нормативному документу -«Правила охраны сооружений и природных объектов от вредного влияния подземных горных разработок на уголь-
ных месторождениях» [1]. Значения углов сдвижения, граничных углов и углов полных сдвижений, а также типовые распределения сдвижений и деформаций в мульде, полученные как средние результаты многочисленных наблюдений за сдвижением земной поверхности в определенном угольном бассейне или его части. При этом основная часть наблюдений производилась в 60-70-х годах прошлого векапри скоростях подвигания очистных забоев не более 5м/сутки и размерах выемки вкрест простирания до 80-120м. Хотя большинство ученых сходилось во мнении, что углы сдвижения имеют только локальное значение, то есть определяют точку сдвижения на земной поверхности для определенного очистного забоя и глубины разработки, в нормативные документы по прогнозному расчету сдвижений вошли четко фиксированные значения углов сдвижения вне зависимости от глубины разработки и скорости подвигания очистных забоев. Также обстоит дело и с типовыми распределениями.Кроме того, используемый метод типовых кривых предназначен для условий закончившегося процесса сдвижения подрабатываемой земной поверхности и не учитывает его развитие в динамике.
В настоящее время средняя скорость подвига-ния высокомеханизированных очистных забоев достигает 15м/сутки и более. Выросли и глубины разработки и размеры выработок вкрест простирания. На современном этапе развития науки о сдвижении земной поверхности внимание ученых и практиков сосредоточено на вопросах именно развития процесса сдвижения во времени, так как всякий объект, любая точка земной поверхности, попадая в зону влияния очистных работ, подвергается, прежде всего, воздействию изменяющихся сдвижений идеформаций, которые вызваны прохождением под ней очистного забоя.
Для условий Кузбассабыли проведены инструментальные наблюдения за сдвижением и деформациями земной поверхности над движущимся очистным забоем при разработке пологих и наклонных угольных пластов в новых условиях -с повышенными скоростями подвигания и увеличенными размерами лав вкрест простирания. Создана база инструментальных наблюдений в виде экселевских таблиц ведомостей сдвижений и деформаций по профильным линиям наблюдательных станций на различных шахтах Кузбасса, на основании которой А.С. Ягуновымрассмотрена динамика формирования профиля мульды сдвижения в зависимости от подвигания очистного забоя [2]. При этом введены термины:
- динамическая мульда сдвижения - состояния мульды сдвижения в процессе ее формирования в зависимости от положения очистного забоя лавы относительно монтажной камеры;
- динамическаямикромульда сдвижения - состояния динамической мульды сдвижения при
отходе очистного забоя лавы от монтажной камеры до образования «плоского дна»;
- динамическое оседание земной поверхности цд(мм; м) - вертикальная составляющая векторов сдвижения точек в динамической мульде сдвижения в зависимости от положения очистного забоя лавы относительно монтажной камеры;
- динамическое горизонтальное сдвижение земной поверхности ^д(мм; м) - горизонтальная составляющая векторов сдвижения точек в динамической мульде сдвижения в зависимости от положения очистного забоя лавы относительно монтажной камеры;
- полумульда«позади забоя»ЬПЗ - полумульда со стороны монтажной камеры;
- полумульда«впереди забоя»ЬВЗ - соответственно со стороны демонтажной камеры.
На основе проведенных инструментальных наблюдений намиисследованы закономерности оседаний и горизонтальных сдвижений отдельных точек поверхности при прохождении под ними очистного забоя [3 - 8]. Введенные в работе [2] термины сохранены. Установленные закономерности позволили разработать методику прогноза ожидаемых динамических сдвижений отдельных точек земной поверхности при отработке пологих угольных пластов Кузбасса высокомеханизированными одиночными лавами, основные положения которой приведены в данной статье.
Исходные данные для прогноза
Разработанная методика прогноза опирается на следующие исходные данные: т - мощность пласта, м;а - угол падения пласта, град.;Нср -средняя глубина разработки, м;П2- длина лавы по простиранию пласта (Б2х - ее текущее значение), м; длина лавы вкрест простирания пласта, м;с - скорость подвигания забоя, м/сут.;Л н - мощность наносов, м.
Кроме того, методикой используются угловые параметры Кузнецкого бассейна, установленные «Правилами охраны ...» [1]:
- граничные углы для пологих и наклонных пластов, градусы: Д = 750 - 0.9а , у0 = 700 , 80 = 700 , р0 = 450 (в обводнённых наносах и плывунах < = 30°);
- углы сдвижения, градусы: Д = 820 — а , у = 800 , 8 = 800, р = 550 (в плывунах и обводнённых галечниках р = 450);
- углы максимального оседания и пол-ныхсдвижений, градусы: в = 900 — 0.5а щ= 500 — 0.25а, = 500 + 0.38а, щ= 500.
Также применяются параметры, установленные по данным инструментальных наблюдений:
- граничный динамический угол сдвижения в
коренных породах впереди забоя — 77" (установлен А.С. Ягуновым [2]);
- граничный угол ¿>0Af со стороны монтажной
камеры при формировании динамической микромульды, зависящий от положения очистного забоя:
D,
при
2 х
H,
< 0.3 snu = 80 ■ e_1-2 + 70
CP
D
при 0.3 < 1.6
HCP
D,
S0M = 80■ exp(—4 ■ —^) + 70
D,
при
H,
H
> 1.6 ¿0M =^0
CP
-динамический угол дпм , определяющийпо-ложение максимального оседания микромульды вдоль лавы при ее формировании:
А
при
H
2х < 0.3 епм = 96.5
CP
при 0.3 < < 1.6
- CP
D„
D
при
2 х
H
Опм =—15 +101
> 1-6 впм = 77
CP
- средние динамические углы полныхсдвиже-ний у/ш = 60.30 и = 44.20 , достигнутые при формировании микромульды.
Параметры динамической мульды и её характерные зоны.
- развивающаяся динамическая мульда по простиранию dзх при различных положениях очистного забоя, оцениваемых расстоянием D2x от монтажной камеры, м
^ = Ах + (нср — К) ■ (+ с^о)+;
+ 2КН ■ оХ%Уо; Ах е [ 0. знСР; А ]
- мульда по простиранию при полной отработке лавы L3, м
ц = А + (нср - к ) ■ + сш50) + 2Кн ■
- микромульда при достижении максимальных оседаний по простиранию пласта Lзo, м
Ц30 =16 нСР + (Нср — к ) ■
(+ ^0 ) + 2Кн ■
- полумульда позади забоя Lпз, м
пРи D2 х > 1.6 ■ HCP
Ln3 = НСР ■ ctgW3M + ;
+ (НСР — Кн )ctg80 + hHCtgP0
ПРИ D2 х < 1.6 ■ HCP
Ln3M = 0 5D2х + (НСР — Кн )ctg¿0М + ;
KctgP0 — HCP ■ ctg°nM - полумульда впереди забоя Lb3, м
при D2x > 16 ■ HCP
LB3 = НСР ■ ctgW3Д + (НСР — hH )ctgS'o + huctg^0 пРи D2X < 1.6 ■ НCP
LB3M =
0.5D2 х + (Нср — К )ctgö0+.
+ hHCtgP0 + HCP ■ ctg°>
'nM
при
■ плоское дно по простиранию Lп.днo, м
т = т — т — т
ЦП. ДНО т3 тПЗ тВЗ;
■ полумульды по падению LП и восстанию Ьв
А
Н
< 1.6.
м;
CP
Ln = hH • ctgp0 + (H СР + 0,5Dj sin a- hH) •. ctgfa + 0,5 A cos a - Hcp • ctgO LB = hH • ctgp0 + (Hcp - 0,5 A sin a- hH) •. ctgy0 + 0,5 A cos a + Hcp • ctgO
- полумульды по падению ЬПд и восстанию ЬВд
Dj
при
Н.
> 1.6,
м
CP
ЬПД = hH ■ ctgp0 + (НСР + 0,5JD1 sin a- hH) ctgP0 + (0,5D1 sin a + НСР ) ■ ctg( щ +a)
1ВД = hH ■ ctg@0 + (НСР - 0,5D1 sin a- hH ) ■ ;
ctgy0 + (-0,5D1 sin a + НСР ) ■ ctg(щ2 - a) - плоское дно вкрест простирания, м Т = Т + Т - Т - Т
LBK. ДНО ТП + ТВ ТПД ТВД .
В динамической мульде вдоль подвигания очистного забоя (по простиранию лавы) выделены четыре характерные зоны, в которых различны закономерности изменения вертикальных и горизонтальных составляющих сдвижения отдельных точек поверхности (рис. 1).
Размеры зон La, Lb, Lc, Ld вычисляются по
эмпирическим граничным углам: SA = 73° ,
SB= 430, 8^= 600 и Sn= 830.
х
Рис. 1 - Характерные зоны сдвижения А, В, С и Ви их длины ЬА, Ьв, Ьс, Ьо в сформированной микромульде при положении очистного забоя Д2х = \6Нср на разрезе по простиранию лавы.
Задание положения точки М(х) в мульде сдвижения производится в системе координат XOZx, показанной на рис.1. Начало координат расположено в центре монтажной камеры, ось 0Х направлена по центру лавы в направлении подви-гания очистного забоя. От монтажной камеры в сторону подвигания забоя координата х положительна, в обратную сторону - отрицательна. По оси Zx откладываются величины оседаний точек поверхности.
Точка М:(х1) находится в зоне А при условии: 0 < х1 + Ах < ЬА . При этом все точки зоныА находятся в полумульдеЬпз.
Точки зоны В могут находиться в точка М2(х2) находится в полумульдеЬПЗ , если ЬА<х2+Ах < Ьпз;
- условие наличия плоского дна:Ь3>ЬвЗ + ЬПЗ;
- при отсутствии плоского дна точка М2(х2) находится в полумульде ЬВЗ , если ЬПЗ<х2+Ах < Ьа+Ьв;
- при наличии плоского дна точка М2(х2) находится в плоском дне, если Ьпз<хт+Ах < Ьпз+Ьп.дно;
- при наличии плоского дна точка М2(х2) находится в полумульдеЬвз если:
Ьпз+Ьп.дно<х2+Ах < Ьп.дно+Ьа +Ьв.
Точка М3(х3) находится в зонеС при условиях:
- при отсутствии плоского дна если: Ьа+ЬВ<Х3+АХ < Ьпз+Ьвз - Ьо;
- при наличии плоского дна если:
Ьа+Ьв+Ьп.дно<Х3+АХ < Ьпз+Ьвз+Ьп.дно - Ьо
Точка М4(х4) находится в зоне Б при условиях:
- при отсутствии плоского дна если: ЬА+Ьв
+Ьс<х4+Ах < Ьпз+Ьвз;
- при наличии плоского дна если: Ь3 - Ьо<х4 + Ах < Ьз .
Точка Мк(хк) находится вне зоны влияния если
0 >хк + Ах>Ьз .
В зоне А максимальное оседание произвольной точки М1(х1) рассчитываем по формуле (1): ЛтА =г,т • [0-382 •
(
arcctg
- 8.045 • X +АХ + 4.25
Л
(1)
L
•пз
- 0.089]
Л m =
0.7 + 0.25 H п
H
■ m • cosa • N • N
Динамические оседания точки в зависимости от положения очистного забоя рассчитываются по формуле:
Л,А {ь) = ЛтА • (1 — е™); 1х е[0;1] (2)
При этом положение очистного забоя йх от монтажной камеры находится по выражению:
(Щх) = (0.3 +1.3 • 1х) • НСР
Для представления динамических оседаний заданной точки в миллиметрах величина Т]х а 1)
умножается на -103.
Зона В по характеру динамических оседаний точки М2(х2) делится на подзоны:
- точки в полумульде LПЗ. Максимальное оседание "П1^пз вычисляется по формуле (1). Динамические оседания - по формуле (3):
, (1{) = ЛтВЬпз • (3)
Л
'x2BLn3
[0.362 • arcctg(-13.3 • lt + 3.5)- 0.1]
L
dx{lt ) = (0.3 +1.3 • lt ) • Hcp
- точки в плоском дне при его наличии. Максимальное оседание точек равно / Динамическое оседание вычисляется по формуле (4):
ЛХгвто {lt) = rÉ -[0.362• arcctg(-13.3 • lt + 3.5)-0.1]
(4)
dX(lt) = x2 - 0.3 • Hcp +1.3 • lt • Hcp
- точки в полумульде LB3 при отсутствии плоского дна. Максимальное оседание точки вычисляется с учетом скорости подвигания очистного забоя по формулам (5):
rmBLB3 (lt) =Лт • [0.392 • arcctg
'7.487 • Х2 +Ах - Ln3 - 11.35^
lt>n
- 0.118]
при V03 < 5 м/сут
Лmblb3 (lt) =Лт • [0 395 • arcctg
7.17 • Х2 + Ах - Ln3 - 10.4681 - 0.111 ]
LK
при V03 = 10 м/сут(5)
ЛmBLB3 (lt) = Лт • [0 43 • arcctg
5.333 • Х2 + Ах Ln3 - 7.279
- 0.138]
при Уоэ > 15 м/сут
При промежуточных значениях скоростей вычисления производятся методом интерполяции.
- 150
- 200
A B C D
50 - 25 0 25 50 75 100 125 150
Динамическая мульда оседаний в главном сечении по простиранию, м
П, мм Динамика оседания заданной точки п, мм
0г
п max = 76 мм
50
60 70 80 90
Положение очистного забоя от монтажной камеры, м
100
Ç, мм Динамические горизонтальные сдвижения в зависимости от оседаний точки
Çmax = 15 мм
20
40
60
П, мм
80
Динамика оседаний и горизонтальных сдвижений точки x. П, мм 0-
20
40
60
80
10 15
Ç, мм
Рис. 2 - Расчет к примеру. Микромульда не достигла полных сдвижений полумульдахЬПЗ, ЬвЗ или в плоском дне при его наличии. Для определения их положения необходимо произвести проверку следующих
условий:
nm мм D2x = 100 м х = 10 м Цп = 194,5 м Цщо пр. = 0
0
0
5
о
- 500 1000 1500 2000
Пт, мм D2x = 500 м х = 250 м L^ = 633 м Ьдно =180 м
A B C D
Динамика оседаний и горизонтальных сдвижений точки x.
т|, мм 0
-100 0 100 200 300 400 500 600
Динамическая мульда оседаний в главном сечении по простиранию, м
П, мм Динамика оседания заданной точки п, мм
500
1000
1500
-2000
Птах = 152 мм
150 200 250 300 350 400 450
Положение очистного забоя от монтажной камеры, м
500
Ç, мм Динамические горизонтальные сдвижения в зависимости от оседаний точки 0
- 50 -100 -150 - 200
- 250
- 300
- 350
^тах = 305 мм
- 200
400
- 600
800
1000
1200
1400
1600
0 200 400 600 800 1000 1200 1400 1600 - 350- 280- 210-140 - 70 0
1 мм
Рис.3 - Расчет к примеру. Динамическая мульда сдвижений имеет «плоское дно».
йх{11)
Динамическое оседание вычисляется по фор муле (6):
(lt) = VmBL«., ■ [0 362 • arCCtë
(6)
dx(lt ) =
(-13.3 • lt + 3.5)-0.17
— X-)
0.3 • HCP +1.3 • lt • Hcp ;
- точки в полумульде LВЗ при наличии плоского дна. В этом случае из координаты точки х2 вычитается длина плоского дна и вычисления проводятся по формулам (5) и (6).
В зонах C и D произвольные точки поверхности М3(х3), М4(х4) находятся в полумульдеLВЗ. Максимальные оседания точек определяются по формулам (5) с заменой переменного нормированного сомножителя аргумента обратной тригонометрической функции. Динамические оседания вычисляются по формулам (7):
ЛХгс (и) = ЛтС ■
[ 0.37 • агс^(-10 • И + 4.1) - 0.09 ] Ъл 1 )■ -[0.03- в3'536'" - 0.03] (7)
= х3 4 - 0.3 - Яс/) +1.3 - // - яс/)
Горизонтальные динамические сдвижения произвольной точки поверхности в главном сечении по простиранию определяются через ее динамические оседания: - в зоне A:
=Лт 'К-ва^ -а0)
(8)
- в зоне B при сдвиге в сторону монтажной камеры:
= ■ -а -в" -а0) (9)
- в зоне B при сдвиге в сторону демонтажной камеры:
= ■ - а -вауЛв - а2) (10)
- в зонах C и D вычисления производятся по формулам (9) при соответствующих оседаниях ■с ^ и своих эмпирических коэффициентах.
Эмпирические коэффициенты ас, Э1, a2 функций (8)-(10) для зон A, B, C, Dопределены методом наименьших квадратов.
Расчет динамических оседаний и горизонтальных сдвижений представляет собой многова-
0
риантные вычисления с условными и безусловными переходами на различные формулы и замены эмпирических коэффициентов. Разработан алгоритм расчета и записан на языке программирования компьютерной системы Mai^ad. Это удобно тем, что сразу получаем результат вычислений, легко изменить или добавить в алгоритм элементы расчета и посмотреть промежуточные результаты.
Аналогичные закономерности и эмпирические зависимости установлены для точек поверхности в главном поперечном сечении и в любой точке зоны влияния очистных работ.
Пример: отрабатывается лава шириной 200м, длиной 1000м на средней глубине 200м со скоростью подвигания очистного забоя 15м/сут. Угол падения пласта 10°, вынимаемая мощность 3 м.
На рис. 2 приведен расчет ожидаемых оседаний и горизонтальныхсдвижений точки поверхности в зоне А. Точка хнаходится на расстоянии 10м от монтажной камеры, при положении очистного забоя 100м от монтажной камеры. Микромульда не достигла полных сдвижений, её длина 194,5м. Оседание точки начинается при положении очистного забоя в 60м от монтажной камеры и далее точка оседает по убывающей экспоненте до максимальной величины nmax = 76мм. Горизонтальное сдвижение точки в зависимости от величины оседания растет почти линейно до максимальной величины ^max = 15мм. Точка сдвигается в
сторону демонтажной камеры.
На рис. 3 приведен расчет ожидаемых оседаний и горизонтальных сдвижений точки поверхности в зоне В. При отходе очистного забоя на 500м от монтажной камеры динамическая мульда сдвижений размером 633м имеет «плоское дно» длиной 180м. Исследуемая точка х находится на расстоянии 250м от монтажной камеры в зоне В. Оседание точки начинается при положении очистного забоя в 200м от монтажной камеры, и далее точка оседает по сигмоидедо максимальной величины Птах = 1527мм. Горизонтальное сдвижение точки в зависимости от величины оседания можно трактовать как возвратно поступательное до максимальной величины ^тах = 15мм. Точка сдвигается сначала в сторону демонтажной камеры до максимальной величины ^тах = 305мм по убывающей экспоненте, затем возвращается по возрастающей экспоненте.
Рассматриваемая методика прогноза ожидаемых динамических сдвижений отдельных точек подрабатываемой земной поверхности используется экспертами СИГИ в качестве дополнения к «Правилам охраны ...» [1] при разработке мер защиты зданий, сооружений, коммуникаций и природных объектов от вредного влияния горных работ. Это повышает научную обоснованность и надежность принимаемых решений.
СПИСОК ЛИТЕPAТУPЫ
1. Правила охраны сооружений и природных объектов от вредного влияния подземных горных разработок на угольных месторождениях. СПб, ВНИМИ, 1998, 290 с.
2. Ягунов A.C Динамика деформаций в подрабатываемом горном массиве / Минэнерго РФ, НИИ горн.геомех. и маркшейд. дела - Межотраслевой научный центр ВНИМИ, Сибирский филиал. - Кемерово: Кузбассвузиздат, 2010. - 239 с.
3. Свирко С.В., Ренев A.A. Особенности смещения точек подрабатываемой земной поверхности высокоскоростными очистными забоями при отработке полого-наклонных угольных пластов Кузбасса // Вестник КузГТУ, 2015. №5. С.32-42.
4. Свирко С.В., Быкадоров A.^, Ренев A.A. Закономерности сдвижения отдельных точек подрабатываемой поверхности в условиях Кузбасса // Геомеханика в горном деле: доклады Всероссийской научно-технической конференции с международным участием 2-4 декабря 2015г. - Екатеринбург: ИГД УрО РЛН, 2015.
5. Быкадоров AÄ, Ларичкин П.М., Свирко С.В., Ренев A.A. Динамика вертикальных и горизонтальных составляющих сдвижений земной поверхности вкрест простирания лав при отработке пологих и наклонных пластов Кузбасса // Вестник КузГТУ, - 2016, №1. - С. 25-33.
6. Свирко С.В., Ренев A.A. Закономерности развития оседаний земной поверхности вкрест простирания лав при отработке пологих и наклонных пластов Кузбасса с высокими скоростями подвигания забоев // Вестник КузГТУ, - 2014, №6. - С. 23-27.
7. Свирко С.В. О влиянии скорости подвигания очистного забоя на процессы сдвижения земной поверхности. // Вестник КузГТУ, - 2016, №3 - С.51-62.
8. Alexey Renev, Sergey Svirko, Alexey Bykadorov, Valery Fedorin «The influence of advancing speed of powered mining stope with single face on earth's surface displacing in Kuzbass» // «Environment, Energy and Earth Sciences (E3S) Web of Conferences», 2017, Volume 15, 01002.
REFERENCES
1. Pravila ohrany sooruzhenij i prirodnyh ob#ektov ot vrednogo vlijanija podzemnyh gornyh razrabotok na ugol'nyh mestorozhdenijah. SPb, VNIMI, 1998, 290 s.
2. Jagunov A.S. Dinamika deformacij v podrabatyvaemom gornom massive / Minjenergo RF, NII gorn.geomeh. i markshejd. dela - Mezhotraslevoj nauchnyj centr VNIMI, Sibirskij filial. - Kemerovo: Kuzbassvuzizdat, 2010. - 239 s.
3. Svirko S.V., Renev A.A. Osobennosti smeshhenija tochek podrabatyvaemoj zemnoj poverhnosti vyso-koskorostnymi ochistnymi zabojami pri otrabotke pologo-naklonnyh ugol'nyh plastov Kuzbassa // Vestnik KuzGTU, 2015. №5. S.32-42.
4. Svirko S.V., Bykadorov A.I., Renev A.A. Zakonomernosti sdvizhenija otdel'nyh tochek podrabatyvaemoj poverhnosti v uslovijah Kuzbassa // Geomehanika v gornom dele: doklady Vserossijskoj nauch-no-tehnicheskoj konferencii s mezhdunarodnym uchastiem 2-4 dekabrja 2015g. - Ekaterinburg: IGD UrO RAN, 2015.
5. Bykadorov A.I., Larichkin P.M., Svirko S.V., Renev A.A. Dinamika vertikal'nyh i gorizon-tal'nyh sostavljajushhih sdvizhenij zemnoj poverhnosti vkrest prostiranija lav pri otrabotke polo-gih i naklonnyh plastov Kuzbassa // Vestnik KuzGTU, - 2016, №1. - S. 25-33.
6. Svirko S.V., Renev A.A. Zakonomernosti razvitija osedanij zemnoj poverhnosti vkrest prostiranija lav pri otrabotke pologih i naklonnyh plastov Kuzbassa s vysokimi skorostjami podviganija zaboev // Vestnik KuzGTU, - 2014, №6. - S. 23-27.
7. Svirko S.V. O vlijanii skorosti podviganija ochistnogo zaboja na processy sdvizhenija zemnoj poverhnosti. // Vestnik KuzGTU, - 2016, №3 - S.51-62.
8. Alexey Renev, Sergey Svirko, Alexey Bykadorov, Valery Fedorin «The influence of advancing speed of powered mining stope with single face on earth's surface displacing in Kuzbass» // «Environment, Energy and Earth Sciences (E3S) Web of Conferences», 2017, Volume 15, 01002.
Поступило в редакцию 12.10.2017 Received 12.10.2017
