Относительно мощности изгибаемого слоя и расстояния до нейтральной линии можно сказать следующее: на рис. 141 показаны три элементарные мульды оседания, образовавшиеся на земной поверхности под влиянием выемки элементов очистной выработки А , В и С. Исходя из форм этих элементарных мульд, можно определить долю участия элементов А, В и С в деформации кривизны у центральной точки Р. Теоретически деформирование породного массива должно начаться уже при выемке весьма малой части площади очист ной выработки и привести к возникновению на земной поверхности элементар ной мульды сдвижения, так что даже при весьма малом развитии очистных работ в породном массиве должен образоваться нейтральный слой. Если пред ставить взаимосвязь переменной мощности изгибаемого слоя d с радиусом площади полной подработки Л, в свою очередь зависящим от глубины раз работки и деформационных свойств горных пород, в виде линейной зависи мости
d = С (R — r), |
(264) |
где С — показатель деформации, меньший единицы, то доля участия отдель ных элементов очистной выработки в возникновении кривизны земной поверх ности у точки Р может быть охарактеризована глубиной деформирования,
равной CR для центрального элемента А |
(г = 0), С (R — гв) для среднего |
||
элемента В жО для краевого элемента (г = |
R). |
|
|
Если подставить выражения (263) и (264) в уравнение (260), то искомая |
|||
функция |
распределения горизонтальных |
сдвижений по |
оси х , выведенная |
из кривизны, будет иметь вид |
|
|
|
kx = ^ C ( R - r ) = - “g f ( R - r ) ( R * - r % |
(265) |
||
где к — постоянная; |
|
|
|
к - |
ЯД* |
|
(266) |
|
|
||
Построить интеграционную сетку для расчета горизонтальных сдвиже ний можно, если в каждой точке площади полной подработки отложить верти кальные отрезки, равные составляющим горизонтального сдвижения и по лученное тело влияния для горизонтальных сдвижений разбить вертикаль ными цилиндрическими и плоскими секущими плоскостями на кольцевые зоны и секторы равного объема раздельно для двух половин площади полной подработки, лежащих слева и справа от линии DD1. Такое тело влияния имеет только одну плоскость симметрии по оси я, в направлении которой произво дится расчет горизонтальных движений точки Р\ составляющие горизонталь ных сдвижений кху по другим направлениям, получающиеся при вращении кривой функции распределения [уравнение (265)] относительно точки Р\
1 Еще проще можно получить значения зональных радиусов, если площадь, ограни ченную кривой гкх и осью х%разбить на пять равновеликих полос, подобно тому, как это сделано на рис. 100 [284].
должны в соответствии с уравнением (253) проектироваться на ось х путем умножения на cos б, так что на оси DD они уменьшатся до нуля (см. рис. 141). Площадь любого полуцилиндрического сечения с зональным радиусом г опре деляется по формуле
+Л/2 |
(267) |
|
A = rkx l |
cos6d6 = rfcx [sin6]iJ5/2 |
|
- Л / |
2 |
|
А = 2гкх. |
|
(268) |
Объем половины тела влияния для горизонтальных сдвижений, соответ ствующий полному горизонтальному сдвижению, при концентрических се чениях возрастает от центра к периферии, в соответствии с уравнением (151), следующим образом:
|
R |
R |
R |
|
К ху = § A dr = 2^ rkxdr = ^ r ^ r2( R - г) (Д2 - г2) dr. |
(269) |
|||
|
0 |
0 |
о |
|
Решение |
интеграла дает |
окончательно |
|
|
К ху — |
8кС |
|
(270) |
|
Д4 |
|
|||
|
|
|
||
Это уравнение служит основой для разбивки кольцевых зон интегрирова ния сетки (см. табл. 12), которые затем только нужно разделить на секторы равного влияния пропорционально функции синуса. Если пренебречь малыми значениями горизонтальных сдвижений
ихуп= — 0лЪ82Сияп |
(271) |
(при С = 0,25 получается, что vxyn = 0,1у2п)» то зональное деление, в соответ ствии с уравнением (259), оказывается вполне пригодным для многих производственных расчетов горизонтальных сдвижений. Для получения близ ких к действительности горизонтальных сдвижений в условиях Рурской об ласти ФРГ нужно принимать С = 1,1 [416].
Выводам уравнений (257)—(270) здесь было уделено такое большое место потому, что в них входит мощность изгибаемого слоя, как новая характери стика свойств горных пород. Вместо перемещения точек к центру тяжести очистной выработки, связанного с представлением о породном массиве, как о несвязной зернистой среде, здесь рассматривалась наглядная модель про гиба породного слоя. В связи с этим горизонтальное сдвижение оказалось зависящим не только от конвергенции аМ в очистной выработке, но и от мощ ности ближайшего к земной поверхности изгибаемого породного слоя, дефор мирующегося как единое целое, а также от линии прогиба, подстилающей этот
слой толщи коренных пород. Однако, вследствие деформаций сдвига по поверх ностям скольжения, образующимся на некоторых межслоевых контактах, действительная мощность поверхностного изгибаемого слоя в слоистом породном массиве несколько меньше, чем в сплошном массиве, где нейтральный слой находится примерно в нижней трети глубины разработки (см. рис. 23). Не будет большой ошибки, если мы для сложенного рыхлыми породами верхнего слоя, покоящегося на слое крепких пород, для всей области влияния очистных работ примем в уравнении (265) постоянную мощность изгибаемого слоя d, так что объем половины тела влияния будет возрастать в направлении радиуса в соответствии с выражением
(272)
о
Эксперименты на моделях, в которых породный пласт имитировался пла стиной пенопласта, показали также, что на горизонтальные сдвижения оказы вают влияние мощность породного слоя и трение на межслоевом контакте. Пока еще остается невыясненным, в какой мере уменьшение пористости и взаим ные перемещения частиц грунтового слоя с пластическими свойствами способны компенсировать горизонтальные деформации, возникающие в верхнем пород ном слое без возникновения предельных напряжений.
9.3.3.
Интеграционные сетки для расчета горизонтальных сдвижений при наклонном залегании пласта
Так же как и при расчете оседаний, в виде опыта можно выполнить расчет горизонтальных сдвижений для наклонно залегающего пласта при помощи овальной интеграционной сетки, трансформированной из обычной (круговой) сетки путем центрального проецирования [207].
Задачу построения сетки можно решить и иначе, если, по аналогии с урав нением (160), принять функцию распределения372
(273)
для составляющих горизонтального сдвижения в направлении падения пласта, обусловленных влиянием элемента очистной выработки, и проинтегрировать ее по половине площади полной подработки.
На рис. 142 показаны интегральные кривые ку в плоскости падения пласта, построенные в соответствии с уравнением (273) для углов падения, равных 50 и 70°. По этим кривым могут быть определены зональные углы отдельных значений kL и отложены на разрезе от вертикальной оси РР\ как показано для примера на рис. 143, для заданных граничных углов. В плоскости прости рания, проходящей через точку Р\ для построения зон сетки использованы Зо-
Рис. 142.
Интегральная кривая, описываемая уравнением (273), для разбивки на зоны равного влия ния в плоскости падения пласта [370]:
I я II — участок сетки соответственно по падению и восстанию пласта
Рис. 143.
К определению горизонтального сдвижения при наклонном залегании пласта по методу центра тяжести при полной отработке и зо нальных углов для интеграционной сетки
нальные углы, взятые с рис. 143 для угла а = 0°. В диагональных плоско стях контуры зон могут быть вычерчены или по данным, взятым из специальных таблиц [370], или по точкам эллипсов. Лучи, делящие сетку на секторы в соот ветствии с функцией синуса, должны строиться не из центра сетки Z, как это делается при построении сетки для расчета оседаний, а из точки Р'
Максимальные горизонтальные сдвижения, имеющие место над верхней и нижней границами очистной выработки, при наклонном залегании пласта имеют различную величину (см. рис. 89), причем большее сдвижение возни кает над нижней границей выработки. Это объясняется как тем, что в этом месте мощность изгибаемого слоя пород больше, так и тем, что здесь имеют место нормальные к напластованию сдвижения горных пород, направленные к очистной выработке и создающие в покрывающей толще пород сдвигающие усилия в направлении лежачего бока пласта. Вопрос о порядке величин этих
т а б л и ц а 17 |
|
|
|
Полное горизонтальное сдвижение в плоскости падения, |
|
|
% от полного оседания |
|
Угол падения, градус |
над нижней границей |
над верхней границей |
|
||
|
очистной выработки |
очистной выработки |
0 |
45 |
45 |
30 |
40 |
25 |
50 |
35 |
5—10 |
70 |
30 |
2—5 |
максимальных и полных горизонтальных сдвижений в литературе пока освещен очень мало. Одна из возможностей определить их по известным значениям оседаний заключается в том, что величину оседания в плоскости падения пласта, найденную для точки, расположенной над границей очистной выработки, умножают, в соответствии с уравнением (237), на тангенс угла центра тяжести, отсекающего верхнюю и нижнюю половины площади полной подработки (см. рис. 143). Угол центра тяжести для половины площади полной подработки определяется для соответствующих исходных данных (угол падения пласта, граничный угол) по приведенной на рис. 109 номограмме для определения
зональных углов для |
1/2 |
K zi и 1/2 К 22. Ориентировочное представление |
||
о соотношении обоих |
полных |
горизонтальных сдвижений может |
дать табл. |
|
17, в которой приведены |
их |
значения, вычисленные описанным |
способом. |
|
Кроме описанных выше методов были предложены и другие функции распределения для расчета горизонтальных сдвижений, например основанные
на преобразовании функции оседания Гаусса |
(357) по формуле |
К = К ± - |
(274) |
или эмпирические методы разбивки площади полной подработки на зоны влия ния на основе результатов натурных наблюдений за сдвижением точек земной поверхности [196, 483].
9.4.
Расчет горизонтальных сдвижений на основе величин наклонов земной поверхности
Достаточно взглянуть на кривые параметров сдвижения земной поверхности, показанные на рис. 84, чтобы убедиться, что кривая горизонтальных сдвиже ний по форме сходна с кривой наклонов. У края мульды и в ее середине обе кривые дают значения, равные нулю, а над границей очистной выработки имеют максимумы, причем между экстремальными значениями обе кривые имеют форму, приближающуюся к синусоиде. Это сходство кривых vx и v'z исполь зуется при вычислении сдвижений земной поверхности при помощи ЭВМ —