со скоростью подвигания забоя v (м/сут) и радиусом площади полной подра ботки R (м) [407]:
|
dvzi |
- (max) = 0,6 |
R - v |
(364) |
|||
|
dt |
|
|||||
|
dvхДИН |
(max) = 0,25 |
R |
v |
(365} |
||
|
dt |
|
|||||
|
В СССР для условий Донецкого бассейна применяется формула |
||||||
|
|
*ин (max) = 30 |
, |
мм/сут, |
(366) |
||
где |
v — подвигание |
очистных работ, м/мес |
[158]. |
||||
|
Применяемые в ПНР формулы для расчета максимальных деформаций |
||||||
земной |
поверхности имеют следующий вид [202, 345]: |
||||||
|
для |
наклонов |
|
|
|
|
|
|
^2 max |
0,7-М ,0 |
аRМ |
, |
мм/м; |
(367) |
|
для |
радиуса кривизны |
|
|
||||
|
|
|
|
Д2 |
|
(368) |
|
|
Prmin = |
650 -f- 720 а М |
м; |
||||
для |
деформаций сжатия |
|
|
||||
|
—esmax=— “TjT- |
ИЛИ |
О.бУг max, мм/м; |
(369) |
|||
для деформаций растяжения |
|
||||||
|
e s щ а х = |
0,45 |
, |
мм/м. |
(370) |
||
В формулах (367)—(370) вынимаемая мощность указывается в миллимет рах. Получаемые по этим формулам деформации (за исключением наклонов) больше, чем наблюдающиеся в Рурской области (табл. 20), однако еще большие деформации наблюдаются в Саарском бассейне1 [75, 234].
1 Первые формулы для приближенной оценки общего горизонтального сдвижения А«общ, по профилю мульды оседания были предложены в 1919 г. К. Леманом; они были выведены, исходя из предположения о поворотах породных блоков навстречу друг другу с обеих сторон выемочного поля (Ь — половина длины очистной выработки), причем оседание кровли в се редине очистной выработки принималось равным аМ. Эти формулы имели вид [227]:
оН
—А$общ 1—2яД/ ^ ;
“Ь А$общ — 0,5 (—А^общ) •
Т А Б Л И Ц А |
20 |
|
|
Условия |
Формулы, |
Определяемый |
выведенные |
|
подработки |
аналитически |
|
параметр |
(полная или |
для условий |
сдвижения |
неполная) |
Рурской |
|
|
области ФРГ |
Оседание, мм |
Полная |
а М |
|||
|
Неполная |
а М е |
|||
Наклон, мм/м |
Полная |
. |
- |
аМ |
|
|
|
1 2 |
|
R |
|
|
Неполная |
0,7- |
аМ |
||
|
R |
||||
Сжатие, мм/м |
Полная |
|
|
аМ |
|
0-25- Н |
|||||
|
Неполная |
0,5- |
аМ |
||
|
Н |
||||
Растяжение, |
Полная |
0,25- |
аМ |
||
Н |
|||||
мм/м |
|
|
|
|
|
|
Неполная |
0,25- |
аМ |
||
|
Н |
||||
Минимальный |
Полная |
500- |
R 2 |
||
а М |
|||||
радиус кри |
|
|
|
||
визны, м |
Неполная |
|
|
/?2 |
|
|
250- а М |
||||
Формулы, |
|
полученные |
Место, для которого |
эмпирически |
|
для условий |
определяются параметры |
Саарского бас |
|
сейна ФРГ |
|
аМ Середина мульды над цент
аМ е ром очистной выработки
а А/ 1,9- Над контуром очистной вы
I работки
Середина мульды над цент ром очистной выработки
1,3- |
а А/ |
|
Н |
|
|
0,65 - |
а А/ |
Над серединой зоны опор |
Н |
||
|
|
ного давления |
12
250 |
а М |
Над серединой очистной |
|
|
выработки |
(кривизна вогну |
|
|
|
тости) или над серединой зон ел |
|
|
|
опорного |
давления [(кривизна |
|
|
выпуклости) |
|
П р и м е ч а н и е : а — коэффициент оседания (0,5 при полной закладке и 0,9 при выемке с обру шением кровли); М —вынимаемая мощность (здесь в мм); е — коэффициент влияния, определяемый
при помощи интеграционной сетки; Н —глубина |
разработки, м; я — радиус площади полной |
попра- |
ботки, м (Я = 0,7Н); г — расстояние по горизонтали от края мульды сдвижения до ее середины, |
м (для |
|
неполной подработки площадь выработки принята |
равной 50% площади полной подработки). |
|
В Голландии для оценки ожидаемых деформаций сжатия и растяжения
пользуются |
приближенными формулами: |
|
|
e smax — l f 6 yz max^ |
\ мм/м; |
(371) |
|
max = |
0»8l?z maxZ-1, |
MM/M, |
|
где I — длина области деформаций сжатия или растяжения (м) при выемке площади неполной подработки. Границей между зонами деформации сжатия и растяжения является линия, па которой оседание равно 50% максимального 1468].
|
В Остравском бассейне ЧССР для глубины разработки от 500 до 800 м |
|||||||
получены |
следующие максимальные деформации: |
) |
||||||
для |
наклонов |
|
|
|
|
|
||
|
max = |
4 , 4 ^ |
m ax, |
мм/м; |
(372) |
|||
для |
деформаций |
растяжения |
|
|||||
|
% max = 0,'45l>zmax. мм/м; |
(373) |
||||||
для |
деформации |
сжатия |
|
|||||
|
— ■e s n ia x = 0 , 6 3 y ; max, мм/м. |
(374) |
||||||
В этих формулах максимальное оседание выражается в метрах |
[480]. |
|||||||
|
Аналогичные формулы, применяемые в Великобритании, имеют следу |
|||||||
ющий вид: |
|
|
|
|
|
|
||
ДЛЯ |
наклонов |
|
|
|
|
|
||
|
-- Q/ |
vz max . |
|
(375) |
||||
|
Uzщах —^ |
|
|
j.[ |
, |
|
|
|
для |
растяжения |
|
|
|
|
|||
|
max |
о 8 |
|
max |
•* |
|
(376) |
|
Для сжатия |
|
|
|
|
|
|
||
|
fcs max “ |
О О |
vz max |
» |
(377) |
|||
|
|
|
// |
|
|
|||
где максимальное оседание vZmax, которое в зависимости от величины выемоч ного поля меньше, чем полное оседание аМ , должно переводиться в милли метры. Деформации сжатия достигают критических значений при ширине выемочного поля 0,24Н и его длине 2R. Коэффициенты в первых двУх Фор мулах приведены для выемочного поля длиной 2R и шириной 0,5# [276].
Если принять в формулах для кривизны и наклонов R = 0,7Н и выра зить коэффициенты округленными числами, то можно построить наглядную диаграмму, по которой могут быть ориентировочно определены экстремальные значения параметров сдвижения земной поверхности для М = 1 м в зависимо сти от глубины разработки и способа выемки полезного ископаемого (рпс- 185). Для других значений вынимаемой мощности и коэффициентов полуяенные‘ результаты могут быть вычислены при условии, что их изменение проИсх°Дит по линейному закону, т. е., например, при М = 1,8 м полученные значения наклонов и горизонтальных деформаций нужно умножить, а значения радиуса кривизны — разделить на 1,8. Из этой диаграммы можно видеть, что с ростом глубины разработки деформации земной поверхности уменьшаются. Есл0 изве стны размеры области влияния на земной поверхности или радиус пЛ01ДаДи полной подработки R, то можно определить сдвижения для любых точек про филя по их расстоянию от точек, где сдвижения достигают экстремальных значений. Например, наклон в точке 1 на рис. 186 составляет vz
Рис. 185.
Зависимость экстремальных значений деформации земной поверхности над очистной выра боткой от глубины разработки, мощности пласта и системы разработки:
I — при выемке с закладкой выработанного пространства; II — при выемке с обрушением кровли
Рис. 186.
Схема расположения точек с экстремальными значениями деформаций земной поверхности относительно очистной выработки:
I, II, III, IV — линейное распределение соответственно растяжений, сжатий, кривизны, наклонов; V — кривая оседаний; VI и VII — растяжение соответственно при полной и неполной подработке
1 0 .6 .
Критические замечания по расчету сдвижений земной поверхности
Наука о сдвижении горных пород, опирающаяся на работы сравнительно небольшого круга исследователей-марШнейдеров, за несколько десятилетий XX в. сделала, несомненно, заметные успехи. Если не говорить об особо слож ных горнотехнических и горно-геологических условиях, создаваемых крутым падением пластов, наличием нарушенных зон и высокой концентрацией очист ных работ, то в настоящее время имеется возможность с помощью существу ющих методов расчета определять ожидаемые сдвижения земной поверхности с достаточной для производственных целей точностью (+10% при пологом
и±20% при крутом падении). Однако нельзя не признать, что еще ждут реше ния некоторые проблемы, которые или возникали вновь или же существовали
иранее, но лишь теперь стали актуальными в связи с возросшей степенью
подработанности породных массивов, ростом |
глубины |
разработок |
до 1000 м |
||||
и возведением чувствительных к подработке |
крупных |
зданий и |
сооружений |
||||
в районах горных разработок с непрерывно растущей плотностью населения. |
|||||||
Поскольку в будущем в науке о сдвижении горных пород уже нельзя будет |
|||||||
отказываться от решения этих специальных задач, необходимо внести ясность |
|||||||
в некоторые |
вопросы, вытекающие |
из следующих |
критических |
замечаний. |
|||
1. |
До |
настоящего времени |
проведено |
еще |
слишком мало |
сравнительных |
|
наблюдений влияния геомеханических, геометрических и горнотехнических факторов на процесс сдвижения земной поверхности. В ФРГ не уделяется достаточного внимания сбору данных опыта, например о величине и местополо
жении максимальных значений параметров сдвижения, а также о развитии процесса сдвижения во времени и о соотношении оседаний и горизонтальных сдвижений; кроме того, научная обработка опубликованных результатов изме рений часто затрудняется из-за недостатка данных для привязки этих резуль татов к конкретным условиям горного предприятия.
2. Ни одним методом расчета, за исключением метода конечных элементов, до сих пор не поддавалось учесть фактическое местное распределение конвер генции в краевой зоне очистной выработки и явление вторичного оседания, вызванного влиянием ранее отработанной горной выработки, граничащей с рас четной выработкой (активизация сдвижения, вызванная влиянием старых горных работ). Правда, для методов расчета с помощью интеграционных сеток было предложено компенсировать влияние недостаточной или избыточной кон вергенции в краевой зоне введением различных значений коэффициента оседа ния для двух краевых полос (см. рис. 72) или добавлением к площади очистной выработки полосы шириной от 5 до 15 м [24, 87 J, однако из-за сложности практического осуществления этих рекомендаций их чаще всего не выполняют, -хотя локальное распределение величин конвергенции в выемочном поле как исходных данных, от которых зависит характер процесса сдвижения пород ного массива, должно вводиться в расчет с возможно большей точностью. Кроме того, часто не принимается в расчет влияние околоштрековых бутовых полос на процесс оседания (см. рис. И).
3.Чтобы иметь возможность с достаточной точностью определять ход раз вития процесса сдвижения во времени, необходимо вводить в расчет данные on изменении с течением времени таких факторов, как конвергенции в очистной выработке (см. рис. 171) или деформирования целиков (см. рис. 14), степени разрыхления и последующего уплотнения пород основной кровли в зоне обру шения (см. рис. 7), об уменьшении с ростом степени подработанности массива изгибной жесткости породных слоев, а также данные об изменении во времени деформированного состояния пород в средней и верхней части горного мас сива (см. рис. 173). Все это вероятно, можпо осуществить только теоретиче скими методами с многократным повторением расчета для различных сочетаний влияющих факторов.
4.Наблюдающаяся даже при горизонтальном залегании пласта асиммет рия профиля мульды оседания (см. рис. 165), явно связанная с направлением развития очистных работ, а также зависящая от степени подработанности и ориентировки систем трещиноватости породного массива [324] может быть учтена при расчете путем повышения «веса» ячеек интеграционной сетки в той
ееполовине, которая находится в стороне движущегося забоя, хотя более близкие к действительности результаты можно было бы получить методами, основанными на теоретических моделях горного массива с учетом упругопласти ческих и реологических свойств пород, или же использованием асимметричной функции распределения (см. рис. 43).
5.Дополнительное оседание породных слоев кровли, обусловленное возобновлением процесса сдвижения под влиянием соседних с очистной выра боткой старых (ранее отработанных) выработок, приводит к тому, что оседание земной поверхности над очистной выработкой оказывается больше, чем полу-