Фрактальная особенность структуризации массива горных пород в изменениях давления на призабойную часть угольного пласта, отрабатываемого длинным очистным забоем

Размещено на http://www.allbest.ru/

Института угля СО РАН

Фрактальная особенность структуризации массива горных пород в изменениях давления на призабойную часть угольного пласта, отрабатываемого длинным очистным забоем

М.В. Шинкевич

Из классической геомеханики следует, что движение очистного забоя приводит к периодическому изменению размеров зависающих консолей подрабатываемых слоев пород и их числа, участвующих в формировании опорного давления. Формирование консолей сопровождается поднятием слоев пород над пластом впереди зоны максимальных напряжений. Известен гармонический ряд «волны Вебера» (1920 г.), согласно которому с удалением от очистного забоя амплитуда волны снижается, а ее период увеличивается. Натурными исследованиями Канлыбаевой Ж.М. [1] инструментально доказано, что даже в малопрочных породах (аргиллиты) при глубине залегания отрабатываемого пласта около 100 м и длине очистного забоя 100 м поднятие нижней части кровли впереди забоя достигает 3 см на расстоянии 10-15 м от его плоскости. Черняком И.Л. и Зайденваргом В.Е. [2] аналитически показано, что в процесс вовлекаются и удаленные слои подрабатываемого массива, обладающие своей периодичностью; происходит «топтание» пласта горным давлением.

Отмеченные особенности достаточно близки к изложенным ниже результатам, но эти исследования выполнены на основе иных представлений о геомеханических процессах во вмещающем массиве горных пород. В принятой модели основой развития является не вес пород, а их удельная упругая энергия. Соответственно, причиной разрушения (дезинтеграции) является избыток упругой энергии в определенном слое массива. Снижение энергетического потенциала обеспечивается отторжением соответствующей массы пород в виде сводов параболической формы [3].

Рассматривая углегазоносный массив как геосреду с вложенными в нее природой пластинами-индикаторами (углеметановые пласты), ранее была установлена возможность [4] по показаниям системы аэрогазового контроля шахт отследить динамику изучаемых процессов и тем самым оценить особенности реакции массива на изменение внешних условий при отработке угольного пласта длинными очистными забоями [5]. Подтверждено формирование иерархии нелинейных блочных структур при вертикальной разгрузке углегазоносного массива от действовавших исходных напряжений. Геометрические параметры структур подчиняются правилу суперпозиции с кратностью 2 (фрактальная размерность Ф относительно длины очистного забоя lоч). В разработанной модели геосреда, стремясь к новому устойчивому состоянию в изменившихся внешних условиях, снижает свой энергетический потенциал путем выделения соответствующей массы энергонасыщенной породы в форме параболоидов. Для образования боковой поверхности увеличивающихся сводов необходимая величина полной упругой энергии пород в соответствующем объеме формирующего своды геомеханического слоя может обеспечиваться ростом только их высоты, так как основания условно ограничены квадратами в технологически создаваемом прямоугольном вырезе на отрабатываемом пласте, границы которого безусловны.

Линейное снижение напряжений в среде с приближением к дневной поверхности обусловливает изменение высоты в иерархии сводов с нарастаниием фрактальной размерности. Для удобства алгоритмизации этой части модели представим массив в виде совокупности вложенных геомеханических слоев, отношение мощностей которых hi,в соответствует Ф=2.

Принимая линейное нарастание вертикальных напряжений с глубиной разработки, получим среднее значение упругой энергии на единицу мощности рассматриваемого геомеханического слоя:

, МДж/м,

где k - размерный коэффициент; г - плотность пород, кг/м3; Нр - глубина залегания разрабатываемого пласта, м; Нс - глубина верхней границы рассматриваемого геомеханического слоя, м; Е - средневзвешенный линейный модуль упругости пород в рассматриваемом слое, МПа.

Поскольку нас интересует относительное влияние потенциальной упругой энергии на рост высоты сводов, то условие hi = Фih0, где h0 - высота первого от отрабатываемого пласта геомеханического слоя; i - уровень иерархии, приобретает вид:

м.

Численный расчет по этой формуле при hi,у? Нр показал, что отношение hi,у / 2ih0 возрастает с ростом уровня иерархии. В диапазоне реальных глубин подземной разработки и значений h0 оно соответствует тангенсу углов 45-650 (рисунок 1), что можно принять в качестве объяснения установленного горным опытом факта: величина угла полных сдвижений подрабатываемого массива горных пород до дневной поверхности составляет: непосредственно над пластом ш3 = 450; в пределах средних глубин горных работ (Кузбасс) - 500; больших глубин (Донбасс) - 600.

Рисунок 1 - Изменение угла полных сдвижений с приближением фронта разгрузки к дневной поверхности при Нр = 1500 м и h0 = 10 м

Рассмотрим горно-технологические условия при отработке пласта с труднообрушаемой кровлей.

Отрабатывается пласт 3-3а Алардинского месторождения на глубине 520-660 м. Полная мощность пласта до 5,6 м, вынимаемая мощность 4,5 м с оставлением пачки угля у почвы пласта. Угол падения пласта 140. Длина выемочного столба 1945 м, длина очистного забоя 220 м. Коэффициент крепости угля f =1,5, следовательно, сопротивление сжатию угля 15 МПа. По мощности пласт переслаивается крепким алевролитом, мощность его слоев колеблется от 0,1 до 2,1 м крепостью f = 4-6.

В кровле пласта залегает углистый алевролит мощностью от 0,4 до 0,95 м, средней крепости f = 2-3. Выше находится слоистый зернистый алевролит средней устойчивости, мощностью 3,3 - 20,9 м, крепостью f = 6. Основная кровля представлена устойчивым крупнослоистым песчаником крепостью f = 6-9, мощностью от 1,0 до 13,5 м.

В почве пласта залегает среднезернистый слоистый алевролит крепостью f = 5.

Регистрация давлений в стойках механизированной крепи выполнялась 1 раз в сутки в ремонтную смену. Каждая секция имеет две стойки. Схема периодичности измерений показана на рисунке 2.

На интервале от 517 до 530 м давления в стойках регистрировались на каждой секции крепи из 112, т.е. через 1,8 м, а на всех прочих интервалах - через 5 секций, т.е. через 9 м. Для удобства обобщающего анализа значения давления интерпретированы через мощность слоя пород, вес которого соответствует реакции крепи.

При средней плотности пород кровли, равной 2,5 т/м3, значения высоты слоя, нагружающего перекрытия секций крепи, составляли:

, м,

где Рст - давление в гидроцилиндрах стоек крепи, МПа; Sст - площадь сечения гидроцилиндра, м2; Sв - площадь перекрытия секции крепи, м2.

Результаты обработки данных показаны на рисунке 3, где видно явное проявление фрактальной размерности Ф=2. Имеем периодичность с шагом 50 м, что меньше длины забоя в 4 раза как по его длине, так и по направлению его подвигания (lоч/Ф2). Но внутри этой периодичности имеются более мелкие периоды, которые при осреднении полиномом на большой площади сглаживаются.

Рисунок 2 - Схема регистрации давлений в стойках механизированной крепи при отработке выемочного столба на Алардинском месторождении Кузбасса

Рисунок 3 - Своды давления на секции механизированной крепи при отработке пласта с труднообрушаемой кровлей

Проявление более мелких фракталов представлено на рисунках 4 и 5.

Здесь отчетливо видна параболическая форма свода сдвижений (средняя линия замеренных значений) с длиной основания, несколько большей длины очистного забоя, и высотой 35 м на интервале подвигания от 517 до 530 м, а на интервале от 530 до 553 м - высотой 43 м. Технически возможная максимальная плотность замеров через 1,8 м позволяет заметить (рисунок 6) в средней части очистного забоя небольшое снижение значений, указывающее на присутствие свода с периодичностью 100 м (Ф=2). В условиях вертикальной разгрузки его высота должна быть равна 50 м. В данном случае давление на секции крепи регистрируется у границы зоны радиальной разгрузки. Тогда получаем 50/35 ?, что отвечает условию удвоения энергетического потенциала структурной иерархии массива при радиальной разгрузке.

На первом интервале при контроле давления по длине забоя через 1,8 м проявляется периодичность с шагом 12,5 м (Ф3), а при контроле через 9 м имеем шаг 25 м (Ф2).

Рисунок 4 - Своды давления на секции механизированной крепи при отработке пласта с труднообрушаемой кровлей на интервале подвигания от 517 до 530 м

Рисунок 5 - Своды давления на секции механизированной крепи при отработке пласта с труднообрушаемой кровлей на интервале подвигания от 530 до 553 м при выполнении замеров на каждой пятой секции крепи (9 м)

Рисунок 6 - Расчетные высоты сводов давления на секции механизированной крепи и их средние значения на интервале подвигания от 517 до 530 м при выполнении замеров на каждой секции крепи (1,8 м)

На втором интервале подвигания забоя горное давление на перекрытия секций крепи повышается и регистрируется формирование вложенных сводов только с периодичностью 50 м, равной шагу сводообразования по оси выемочного столба (квадраты). Если исключить приштрековые зоны, то отклонения замеренных значений от размеров вертикального сечения основного свода во всех случаях лежат в пределах ±10 %. По направлению подвигания очистного забоя (рисунок 3) регистрируется повышение давления через 70 м. Это соответствует периодичности 50 м, умноженной на .

Применение полученных результатов требует знания величины h0 - высоты первого геомеханического слоя в конкретных горно-технологических условиях. Это затруднение снимается, если в качестве опорной величины принять шаг первичного обрушения основной кровли пласта, который является технологически значимым показателем и фиксируется для каждого очистного забоя. Его предварительное определение обеспечено в Кузбассе многолетним опытом применения соответствующей методики, разработанной проф. В.И. Мурашевым [6]. Согласно этой методике отношение первичного шага обрушения к вторичным шагам равно 5, что достаточно близко к кратности параметров . Расхождение величин и составляет 12 %, причем в большую сторону. Это может быть объяснено тем, что в эмпирических коэффициентах [6] присутствует влияние взаимодействия сводов сдвижения при отработке длинных столбов с межлавными целиками недостаточной несущей способности. Это влияние всегда сказывается в большую сторону. Тогда для удобства анализа процессов примем величину вторичного шага обрушения основной кровли (четверти первичного шага) в качестве контролируемого размера критического элемента геосреды, который интегрально отражает влияние напряжений в породах, может быть, от их механических характеристик и длины очистного забоя. Для равномерного обрушения кровли количество критических элементов, а также сводов более высокого уровня иерархии в пределах свода с основанием, равным длине забоя, должно быть целым четным числом.

С учетом вышесказанного появляется возможность оценок параметров структурных элементов массива с позиций накопленной ими упругой энергии.

Оценим изменения механических свойств массива, влияющих на периодичность обрушения пород кровли, по длине выемочного столба, используя значения линейного модуля упругости и коэффициента крепости пород кровли междупластья пластов 1 и 3 по оси выемочного столба.

Подготовка данных заключается в снятии информации с плана горных работ и геологических разрезов по штрекам и по разведочным линиям. Величина устойчивого пролета кровли l кр= 50 м принята равной половине первичного шага первичного обрушения при отходе забоя от монтажной камеры. Соответственно, через каждые 50 м по всей длине выемочного столба с разрезов по конвейерному и вентиляционному штреку сняты мощности вмещающих пород - алевролита и песчаника. Рассчитаны средневзвешенные значения коэффициента крепости пород кровли (4) и построены графики их изменений по длине конвейерного и вентиляционного штреков и средних значений по оси выемочного столба (рисунок 7):

,

где f1 - коэффициент крепости алевролита по справочным данным; f2 - коэффициент крепости песчаника по справочным данным; m1 - мощность алевролита, м; m2 - мощность песчаника, м; mпор - суммарная мощность пород, м.

Рисунок 7 - Изменение средневзвешенного коэффициента крепости пород кровли по длине выемочного столба 3-32 на Алардинском месторождении Кузбасса

Удельная (на единицу объема) потенциальная упругая энергия по оси выемочного столба 3-32 рассчитывается как:

, Дж/м3 ; , МПа,

где у - напряжение в рассматриваемой области угольного пласта, МПа; Е - модуль упругости угля в массиве (справочные данные), МПа; Hр - глубина залегания пласта, м.

Средневзвешенная упругая энергия в слое равна:

, Дж/м3,

где - удельная потенциальная энергия алевролита, Дж/м3; - удельная потенциальная энергия песчаника, Дж/м3.

Для нахождения критической энергии слоя, обеспечивающей образование новой поверхности, используем размеры сводов для соответствующих геомеханических слоев на первичном шаге обрушения основной кровли:

, Дж/м2,

где - упругая энергия при образовании сводов, Дж/м2; - объем свода, м3; - площадь боковой поверхности свода, м2.

Изменения по длине выемочного столба критической энергии (с учетом колебаний средневзвешенного коэффициента крепости пород слоя) позволяют оценить энергетические возможности среды по образованию поверхности разрушения показателем относительной устойчивости геомеханических слоев (рисунок 8) вида:

,

где - средневзвешенный коэффициент на первичном шаге обрушения;- упругая энергия пород на первичном шаге обрушения кровли.

Рисунок 8 - Показатель относительной устойчивости геомеханических слоев

Если принять, что геомеханический слой мощностью 50 м является доминирующим в процессе сводообразования, то разрушения в двух вложенных слоях будут системными при соответствии кратности показателей Лк фрактальной размерности. Однако этого не наблюдается. Следовательно, на интервалах существенного отклонения кратности от фрактальной размерности Ф необходимо принудительное обрушение соответствующего слоя.

Библиографический список

1 Канлыбаева, Ж.М. Закономерности сдвижения горных пород в массиве / Ж.М. Канлыбаева. - М.: Наука, 1968. - 108 с.

2 Черняк, И.Л. Периодичность изменения напряженно-деформированного состояния массивов угля и пород впереди очистного забоя / И.Л.Черняк, В.Е. Зайденварг // Горный журнал. Известия высших учебных заведений. - 1993. - №3. - С. 25-28. энергия массив горный порода

3 Полевщиков, Г.Я. Нелинейные геомеханические процессы в динамике метанообильности выемочных участков / Г.Я. Полевщиков, Е.Н.Козырева // Материалы Всероссийской конференции «Геодинамика и напряженное состояние недр земли», посвященной 80-летию акад. М.В. Курлени (с участием иностранных ученых), 3-6 октября 2011 года: в II т. - Т. I. - Новосибирск: Ин-т горного дела им. Н.А. Чинакала СО РАН, 2011. - С. 449-455.

4 Полевщиков, Г.Я. Газокинетический паттерн разрабатываемого массива горных пород / Г.Я. Полевщиков, Е.Н. Козырева // Горный информационно-аналитический бюллетень. - 2002. - №11. -С. 117-120.

5 Шинкевич, М.В. Динамика геомеханических процессов в призабойной части массива при движении длинного очистного забоя / М.В. Шинкевич, Н.В. Рябков, Е.Н. Козырева // Горный информационно-аналитический бюллетень. - 2010. - № 3. - С. 356-359.

6 Мурашев, В.И. Разработка научных основ безопасного ведения горных работ в угольных шахтах на основе исследования геомеханических процессов: автореферат дис. … д-ра техн. наук (05.26.01). - М.: ИГД им. А.А. Скочинского, 1980. - 36 с.

Аннотация

На основании анализа горно-технологических данных о движении очистного забоя показана периодичность реализации упругой энергии массива горных пород, подтверждающая существенность отличий современной нелинейной геомеханики от классических представлений. Показано, что реализация упругой энергии массива приводит к формированию системы блочных структур, последствия чего выявляются при мониторинге давления в стойках механизированной крепи.

Работа выполнена при финансовой поддержке междисциплинарного интеграционного проекта СО РАН № 99, партнерского интеграционного проекта СО РАН № 100 и гранта РФФИ № 10-05-98009-р_сибирь_а.

Ключевые слова: геомеханика, давлениЕ, напряжение, энергия, шаг, обрушение

Analyzed mining technological data on coal face advances exhibited a periodical realization of rock mass elastic energy, justifying the significant difference of modern nonlinear geomechanics from the classical concepts. It is demonstrated that massif's elastic energy realization leads to the formation of block structures and is also detected during pressure monitoring in props of powered supports.

The work is performed under the financial support of Interdisciplinary Integration Project SB RAS № 99, Partner Integration Project SB RAS № 100 and RFFI grant и № 10-05-98009-р_sibir_а.

Key words: Geomechanics, pressure, stress, energy, pich, caving

Размещено на Allbest.ru