книги / Сдвижение горных пород и земной поверхности при подземных разработках
..pdfРис. 1.2. Схема расчета влияния опорного давления на точку в массиве
Рис. 1.3. Распределение напряже ний вокруг выработки:
а — изолинии напряжений в долях максимальной нагрузки; б — распреде ление напряжении в сечении А —А
1.1.2. Сдвижение горных пород при повторной подработке толщи
В Донбассе шахты разрабатывают, как правило, несколько пластов. Влияние горных работ на земную поверхность при этом может быть одновременным и разновременным. Под разновремен ным влиянием понимается такая ситуация, когда время между от работкой пластов в зоне влияния на данный участок земной по верхности превышает общую продолжительность процесса сдвиже ния, регламентированную «Правилами охраны...» [31]. Подра ботка в условиях уже нарушенной толщи имеет ряд особенностей в сдвижении пород и земной поверхности (прежде всего возникно вение активизации сдвижения). Основным проявлением активиза ции является увеличение максимального оседания и деформации земной поверхности при повторной подработке по сравнению с первичной подработкой, осуществленной в идентичных условиях.
Существующее мнение о причине активизации при отработке свиты сводится к тому, что при первичной подработке в массиве возникают устойчивые зависания, которые при подработке другим пластом по тем или иным причинам погашаются и вызывают до полнительные к расчетным оседания и деформации земной поверх ности. Зависания возникают в условиях крепких вмещающих по род, а также при оставлении в выработанном пространстве цели ков. Во многих случаях активизация наблюдается в породах сред ней крепости и даже слабых, а также тогда, когда первый пласт на данном участке отработан без целиков. Возникновение зависа ний в этих условиях невозможно, а активизация все же происхо дит. Следовательно, существующее мнение объясняет лишь част ную причину и не дает объяснение общей причины активизации.
В общем случае причину этого явления можно объяснить сле дующим. В ненарушенном массиве слои горных пород находятся в сжатом состоянии, обусловленном давлением вышележащей тол щи. При отработке угольного пласта часть веса пород над выра боткой в виде опорного давления передается на призабойный мас сив и целики. В' результате этого породы в области полных сдви жений частично разгружаются от нормального давления, что ве дет к их расширению [44]. Часть расширения обуславливается упругим восстановлением, другая часть — разуплотнением массива системами трещин, образующихся в породах кровли после отра ботки пласта. В процессе подработки точки массива проходят че рез область необратимых деформаций [26], поэтому при вторич ном уплотнении пород в области полной подработки слои пород не сжимаются до прежнего размера, т. е. они получают необрати мое расширение. Если такой массив подработать еще раз, то его повторное расширение будет меньше, чем первичное, и оседание земной поверхности в этом случае будет больше. По данным инст рументальных наблюдений, проведенных над горными выработка ми в Донбассе, среднее увеличение максимального оседания при повторной подработке составляет около 30 % [30, 41].
Рассмотрим случай полной подработки земной поверхности пластом 1 в условиях ненарушенной толщи (рис. 1.4). Здесь и да лее будем считать, что по простиранию пласта лавы имеют раз меры, превышающие 1,4Яср, а способ управления горным давле нием осуществляется полным обрушением кровли. Такое условие
подработки
св Ы 0'3,мм/м
вполне реально для Донбасса, на шахтах которого, в основном, применяются столбовая и сплошная системы разработки со значи тельными размерами выемочных полей по простиранию, а число лав, отрабатываемых с закладкой, невелико. Допустим также, что деформации в точках массива являются непрерывными функция ми координат, т. е. выполняется общий принцип механики сплош ной среды. Деформации расширения подработанного массива е* по мере удаления от земной поверхности возрастают от нуля до некоторого значения ео на уровне отработанного -пласта. С доста точной для практики точностью их изменение можно аппроксими ровать линейной зависимостью
6л=8ох/Н. (1.1)
Тогда оседание земной поверхности от выемки первого пласта после окончания процесса сдвижения
|
|
/л |
|
|
Vol = т |
1 —£oi |
J - щ |
- d x = ni\ - e01 H J 2. |
( 1.2) |
|
|
О |
|
|
Из выражения |
(1.2) |
имеем |
|
|
|
е01 = 2 |
(от, — ifoi)/#,. |
(1 .3 ) |
Величину ео выразим в зависимости от горно-геологических па раметров. Для этого используем фактические данные об оседа ниях земной поверхности при полной подработке и закончившемся процессе сдвижения по 59 наблюдательным станциям. В эти дан ные вошли только случаи, когда коэффициент подработаниости земной поверхности вкрёст простирания пласта П\ был близким к единице, т. е. свод зоны полных сдвижений замыкался на земной поверхности.
Высоту |
зоны полных сдвижений обозначим через # 0. В слу |
чаях, когда |
угол падения аф О , нормальная мощность заменяется |
еепроекцией на вертикаль т cos а.
Врезультате корреляционного анализа выяснилось, что вели чина ео находится в наиболее тесной связи с нормальной кратно стью подработки K = H 0/{m co sa). Для выбора наилучшего вида аппроксимирующей функции и определения коэффициентов в уравнении связи использовался регрессивный анализ. Были ис следованы две формы связи: гиперболическая и экспоненциаль ная. Оценка' значимости коэффициентов корреляции производи
лась по критерию Стыодента, оценка адекватности уравнений — по критерию Фишера. Сравнение эмпирических критериев с тео ретическими осуществлялось на 5 % -ном уровне значимости. Наи более адекватной формой связи является гиперболическая зависи
мость (рис. 1.5). Упростим эту зависимость, |
не снижая ее |
точ |
|||
ности в наиболее реальном диапазоне изменения переменных: |
|||||
|
ео= 640т cos alНо, |
|
(1.4) |
||
где H0— n\H— Q,7Di — высота |
зоны полных сдвижений; Dx— дли |
||||
на выработки (лавы), м. |
величина |
ео |
представляет |
собой |
|
Полученная |
таким образом |
||||
максимальную |
остаточную деформацию |
расширения подработан |
ных пород в направлении, нормальном к напластованию, после окончания процесса сдвижения.
Предположим теперь, что поверхность полностью подработана вторым пластом (см. рис. 1.4). В этом случае расширение будет происходить в ненарушенных породах междупластья и в породах, нарушенных отработкой первого пласта. Характер убывания рас ширения в обоих случаях можно принять одинаковым. Отличие заключается в том, что выше глубины залегания отработанного пласта повторное расширение пород будет меньше, чем первичное.
Величина еов (индекс в означает вторичную подработку) свя-
зала с величиной ео постоянным коэффициентом к, зависящим от физико-механических свойств пород. В слабых (пластичных) по родах вёличина к больше, чем в крепких (упругих) породах.
Для определения к поступим следующим образом. Запишем формулу (1.3) для второго и первого пластов:
802= |
2 (m2—rio2)/#2; |
|
eoi = |
2(m1—r]0i)/# i. |
|
Составим отношение |
|
|
eo2/eoi= ( т 2—1402) Нх1 [(т\—rioi) # 2]. |
(1.5) |
Для определения к необходимо, чтобы Hi = H2, т. е. чтобы пла сты сливались (в Донбассе такие условия практически отсутству ют). Используем результаты натурных наблюдений за сдвижени ем земной поверхности, полученные при слоевой отработке пласта
«Верхняя Марианна» в Карагандинском угольном бассейне |
[19]. |
|
В этом случае выражение (Ь5) |
можно записать следующим обра |
|
зом |
|
|
еог/е01= Еов/е 0 = /с= |
( т 2— 1102) / ( « i — Л01), |
( 1 .6 ) |
где п\\ и га2 — мощности отрабатываемых слоев; rjoi и т^ог— мак симальные оседания поверхности после отработки каждого слоя. Обработка данных, приведенных в работе [19], показала, что зна чение коэффициента к, вычисленное по формуле (1.6), колеблет ся в пределах 0,5—0,8. Породы Карагандинского бассейна соот ветствуют наиболее слабым породам Донецкого бассейна. Поэто му для соответствующих пород Донбасса (марки угля Д, Г) коэф фициент к равен 0,6. Исходя из сравнения результатов расчетов оседания с данными наблюдений, коэффициент к следует принять равным 0,4 для пород средней крепости (марки угля Ж, К, ОС) и 0,2 для крепких пород (марки угля Т, А). Таким образом
е0“= /с е о = (0,2— 0,6)ео.
Тогда оседание земной поверхности после отработки второго пласта
|
N.. |
//, |
|
|
Оп'Л== |
J ~~РП~ |
^*01 J |
“0‘- |
2 |
|
д, |
о |
|
|
Площадь треугольника аВс на рис. 1.4 численно равна добавке к максимальному оседанию земной поверхности после отработки второго пласта
Д = e0i J -щ - d x — kz0l j -щ -dx = (1 —Æ)e(II//,,2.
о1 о
Подставляя вместо eoi выражение (1.4) и учитывая, что при полной подработке земной поверхности HQ— H\, получаем
Д = ( 1 — /c)640mi (cos а ) /2.
Причем гп\ выражается в метрах, а Д — в миллиметрах.
Для определения условий возникновения активизации необхо димо знать высоту распространения расширения массива над вы работанным пространством лавы. Данные о высоте распростране ния расширения в подработанном массиве могут быть получены в результате натурных инструментальных наблюдений за сдвиже нием массива, проведенных в подработанных выработках и сква жинах.
Приуправлении горным давлением обрушением кровли под выработанным пространством образуется область полных сдвиже ний. В зависимости от крепости вмещающих пород высота этой области может составлять (0,54-1,0)/)!. В процессе подвигания за боя лавы породы в этой области проходят через зоны опорного давления, разгрузки и повторного сжатия. При этом слои пород испытывают деформации сжатия, изгиба, расслоения, сноса и т. д. По сравнению с другими участками массива, испытавшими влия ние подработки, породы области полных сдвижений являются наи более деформированными и нарушенными многочисленными тре щинами.
При достаточно большой мощности разрабатываемого пласта породы в области полных сдвижений отделяются по контуру этой области от остального массива (под отделением пород понимается появление открытых трещин по контуру зоны). В этом случае на грузки на почву пласта будут определяться весом пород в области полных сдвижений. Слои пород, залегающие выше области полных сдвижений, передают свой вес на целики. Такое перераспределение давления создает условия для сохранения в области полных сдви жений деформаций расширения в течение длительного периода времени после окончания процесса сдвижения. Осносительные ос таточные расширения пород, расположенных выше области пол ных сдвижений, значительно меньше, чем в пределах этой обла сти. Например, по данным наблюдений, проведенных на шахте «Красная Звезда» производственного объединения «Торезантрацит», относительное остаточное расширение пород над областью полных сдвижений в 15 раз меньше, чем внутри этой области. При разработке пологих пластов мощностью до 2,5 м и длине лавы не менее 120 м остаточное расширение подработанных пород распро страняется в вышележащем массиве на расстояние, не превышаю щее высоту области полных сдвижений. Слои пород, залегающих под этой зоной, прогибаются без существенного расширения. По этому породы, нарушенные первичной подработкой, способны ока зывать существенное влияние на процесс сдвижения при повтор ной подработке только в том случае, если они попадают в область полных сдвижений. Если же высота области полных сдвижений над действующей лавой меньше мощности междупластья, то ак тивизации сдвижения не происходит.
С принятых позиций произведен расчет добавочного оседания активизации при неполной подработке земной повёрхности вкрест простирания пласта. Вершина области полных сдвижений в даи-
ном случае не совпадает с земной поверхностью, поэтому исход ная функция (1.1) примет вид
е*= е0 [х— (Н—Но) ]/Н0.
Рассмотрим случай, когда не полностью подработанная толща подвергается вторичной неполной подработке (рис. 1.6), где ф — угол полной подработки. При нисходящем порядке отработки пла-
ул*..су,«л
Рис. 1.6. Схема к расчету активизации при нисходящем порядке отработки
пластов
Рис. 1.7. Схема, иллюстрирующая граничные условия формулы (1.7)
стов, преимущественно применяемом в Донбассе, такой вариант встречается наиболее часто. Увеличение максимального оседания, вызванное активизацией, будет численно равно площади треуголь ника abc:
д = (1 — А) в,,, j я, — (Н, — н^) d x ’ н%-н+
Учитывая, что Но2= # 2Я1(2), находим:
■Д= [(1—к)eoi/2] [# i-f-# 2(rti(2)— 1)], |
(1.7) |
где лгI(2) — коэффициент подработки земной поверхности |
действу |
ющей лавой во втором пласте в направлении вкрест простирания;
8oi — максимальное |
остаточное расширение пород на |
уровне пер |
|
вого |
(ранее отработанного) пласта, определяемое |
по формуле |
|
(1.4) |
с учетом высоты зоны полных сдвижений под лавой в этом |
||
пласте. |
Я) и Н2 имеют размерность в метрах, а Д — |
||
В формуле (1.7) |
в миллиметрах.
Рассмотрим два варианта взаимного расположения лав в раз рабатываемом и ранее отработанном пластах, которые характери зуют граничные условия формулы (1.7). Максимальная активиза ция произойдет тогда, когда пласты 2 и 1 сольются, а границы
действующей |
лавы |
совпадут с границами ранее отработанной |
|
(рис. 1.7,а). |
В этом случае H2= H i = |
H и формула (1.7) прини |
|
мает вид: |
|
|
|
|
,Д= |
0,5 ( 1—к) 8oi/7nU2) = |
320 ( 1—к) гп\. |
Если же Hoz^h |
(рис. 1.7, в), то оседание от активизации ста |
новится равным нулю, т. е. активизации не произойдет. Случаи, показанные на рис. 1.7,6, и на рис. 1.6, являются промежуточны ми [для них и выведена формула (1.7)].
Величина eoi тем больше, чем меньше давление пород в зоне полных сдвижений. Таким образом, с уменьшением длины ранее пройденной выработки величина eoi возрастает, а активизация увеличивается.
Рассмотрим схемы, показанные на рис. 1.8. Если длину лавы, вторично подрабатывающей массив, принять неизменной, а дли ну ранее отработанной лавы уменьшать, то активизация будет из меняться следующим образом.
При переходе от схемы / к схеме // активизация возрастает и достигает максимального значения, так как расширение пород при повторной подработке происходит целиком в пределах зоны полных сдвижений и величина eoi достаточно велика.
Дальнейшее уменьшение D\ (схема III) может привести к то му, что при отработке первой лавы в массиве образуется устойчи вое зависание. Величина eoi при этом будет больше расчетного значения, а функция (1.1) приобретает разрыв непрерывности. Активизация при этом возрастет скачкообразно. Поэтому необхо димо ввести нижний предел изменения величины D ь
Согласно исследованиям, зависание основной кровли над сред? ней частью выработанного пространства при глубине разработки 500—700 м возможно при длине лавы около 60 м. В ближайшие годы средняя длина лавы по бассейну составит 180 м. Следова тельно, длина D | лавы в ранее отработанном пласте, при которой может возникнуть зависание основной кровли, составляет третью часть длины D2 новой, действующей лавы.
Таким образом, расчет активизации по формуле (1.7) следует
производить при условии |
|
£ ,> 0 ,З Д 2. |
(1.8) |
Активизация сдвижения при разработке свиты пологих пластов проявляется в увеличении оседания и деформаций земной поверх ности, а также в увеличении размеров мульды сдвижения. Эти параметры оказывают противоположное воздействие на величины деформаций. В разнообразных вариантах подработки земной по верхности свитой пластов сочетание этих параметров может уве личивать деформации по сравнению с первичной подработкой или,
J8
наоборот, уменьшать их. С этих позиций неблагоприятным являет ся вариант, когда лавы в отработанном пласте на плане распола гаются в пределах границ действующей лавы. В этом случае де формации земной поверхности увеличиваются, а размеры мульды сдвижения остаются неизменными.
Ш гЩ Щ Ш Ж yAy/^V/xSV/AV/Ay/Ay//^y.
Рис. 1.8. Схема к определению влияния соотношения длин лав в отработанном и разрабатываемом пластах на проявление активизации
Максимальное оседание при повторной подработке предлагает ся рассчитывать по формуле
T)Вmах == Ц т а х + 'Д ,
где Tjmax — расчетное максимальное оседание при первичной под работке; Д — добавочное оседание, определяемое по форму
ле (1.7).
Для иллюстрации предлагаемого способа произведен расчет максимального оседания земной поверхности на шахте «Куй бышевская» производственного объединения «Донецкуголь»
(рис. 1.9).
В зону полных сдвижений над лавой I в пласте /с8 попали слои пород, надработанные лавой II в пласте 1\.
Длина лавы в пласте 1\ — 66 м, в пласте /с8— 140 м, поэтому условие (1.8) выполняется. Следовательно, сдвижение горных по род под лавой в пласте /г8 будет сопровождаться активизацией, которую'можно рассчитать по формуле (1.7).
По степени метаморфизма угля (марка Ж) вмещающие поро ды имеют среднюю крепость и коэффициент к можно принять равным 0,4. По формуле (1.4) определим максимальное остаточ ное расширение пород над ранее отработанным пластом
640 |
640 |
л |
, |
е 01 ~ HitiM |
m i c o s a — 198-0,23 |
1 0 ,9 м м /м . |
Вычислим добавочное оседание, вызванное активизацией сдви жения:
Д = — £— е01 [//t + 7/г (яца) 1)] = g—1— 10,9 [198 +
+ 224(0,46 -1)] =2 52 мм.
40 |
35 |
30 |
25 |
20 |
15 |
10 |
5 |
1 |
------- Расчетные оседания и деформации, вычис
ленные по правилам 1972г.
--------Расчетные деформации и оседания, вычис ленные по предлагаемому способу
Рис. 1.9. Графики измеренных и расчетных оседаний н де формаций земной поверхности