CC BY
31
УДК:622.272 М.В. Писаренко
ОПРЕДЕЛЕНИЕ ОЖИДАЕМЫХ ОСЕДАНИЙ ЗЕМНОЙ ПОВЕРХНОСТИ С ИСПОЛЬЗОВАНИЕМ ГИС ТЕХНОЛОГИЙ
Разработан алгоритм, позволяющий быстро и наглядно оценить степень влияния подземных горных разработок на земную поверхность.
Ключевые слова: подработка земной поверхности, параметров мульды, грид-тема, оседание земной поверхности.
Семинар № 1
M. V. Pisarenko
THE ESTIMATING OF EXPECTED LAND SUBSIDENCE BY HYBRID TECHNOLOGIES
The algorithm for fast and visually estimation of the level of underground mining impact on the soil surface is developed.
Key words: undermining, mould parameters, parameter grid, land subsidence.
Wb процессе добычи угля и проведении горных выработок изменяется напряженное состояние горных пород, это приводит к их деформации и сдвижению. Сдвижение горных пород при определенном соотношении размеров выработанного пространства и глубины достигает земной поверхности, и она так же претерпевает деформацию. Часть породного массива, подвергшегося сдвижению под влиянием горных разработок, принято называть областью сдвижения горных пород, о соответствующую часть земной поверхности -мульдой сдвижения.
Границы мульды сдвижения определяются граничными углами. Это внешние относительно выработанного пространства углы, образованные на вертикальных разрезах по главным сечениям мульды горизонтальной линией и линиями, соединяющими границы вырабо-
танного пространства с граничными точками сдвижения (в качестве которых принимают обычно точки, получившие оседание 10-15 мм или наклоны и горизонтальные деформации 0,5х10-3. [1].
При определенном соотношении размеров выработанного пространства и глубины разработки в мульде сдвижения образуется плоское дно, т.е. участок с максимально возможными оседаниями при данной мощности, глубине и угле падения пласта. На этом участке сдвижения носят в основном равномерный характер. Условия, при которых в мульде сдвижения образуется плоское дно принято называть полной подработкой земной поверхности.
Границы плоского дна определяются углами полных сдвижений у/х,^2.^3. При неполной подработке положение точки, имеющей максимальное оседание определяется углом максимального оседания 0.
Расчет параметров мульды сдвижения и ожидаемых оседаний земной поверхности выполняется согласно [1]. Однако процесс вычисления и особенно построение мульды сдвижения в плане с учетом рельефа является весьма трудоемкий процесс, и поэтому ограничива-
ется расчетами в главных сечениях (по простиранию и падению пласта).
Использование интегрированных ге-информационных систем, которые имеют возможности пространственного представления, манипулирования, анализа и визуализации пространственных данных, позволяют решать задачи такого плана достаточно быстро и эффективно.
Для создания модели расчета ожидаемых сдвижений земной поверхности использовались следующие программные продукты:
- MapInfo 6.5 (можно AutoCAD или любую другую ГИС) - для визуализации, привязки к планам горных работ;
- Surfer 8 - расчетов гридов и построения карт;
- Microsoft Office Excel - проведения расчетов.
В качестве исходных данных использовались совмещенные с поверхностью планы горных работ (в представленных материалах использовались данные по ш. Березовская), горно-гео-логическая информация необходимая для выполнения расчетов. Для автоматизации процессов расчета был создан макрос на языке программирования средствами Visual Basic в среде Microsoft Office Excel. Оседание земной поверхности в любой точке мульды сдвижения определяется по формуле [1]:
%у = (1)
где S(zx) и S(zy) — функция типовой кривой оседания, определяемая в зависимости от коэффициентов N1 (для точек главного сечения вкрест простирания пластов) и N2 (для точек главного сечения по простиранию пла-стоАлгоритм вычисления, заложенный в макросе следующий:
- вводятся исходные горнотехнологические данные: глубина отработки, длина лавы по простиранию и по падению, угол падения и мощность пласта, мощность наносов, наличие подработки;
- вычисляются углы сдвижения, максимального оседания и полных сдвижений, коэффициенты N1 и N2, максимальное оседание, длины мульд по простиранию и падению;
- определяются координаты точек мульды сдвижения и ожидаемые оседания согласно формуле 1. При этом вначале координаты определяются в условных системах координат xy, yy, когда за начала координат принимается точка, соответствующая середине выработки, при этом ось Х проходит в главном сечении мульды по простиранию и направлена по простиранию, ось Y проходит в главном сечении мульды вкрест простирания и направлена в строну восстания пласта. Затем координаты пересчитываются в исходную систему координат Хи, YH (в которой представлен план горных работ) по известным в аналитической геометрии формулам:
Xu= xo+xy*cosy - yyxsmp,
Yu= yo+xyxcosy + yyXsiny,
где Хо, уо- координаты центра начала отсчета условных систем координат в истинных, ф - угол поворота координат ( угол образуемый между направлениями оси Х в условных и истинных координат).
В результате расчетов создается табличный файл, в котором хранятся пространственные координаты и значения ожидаемых оседаний земной поверхности в точках мульды.
Этот файл передаются Surfer где строится грид-темы. Грид тема - это растровый эквивалент темы объектов.
Рис. 1. Мульда сдвижения представленная в двух и трех мерном пространствах (для лава16-16 ш. Березовская)
Она представляет собой графический слой, где пространство Вида разделено на квадратные ячейки. Каждая ячейка хранит числовое значение данных, передающее информацию о географическом слое, который она представляет. На основе грид-тем можно проводить следующие аналитические операции: сложение двух гридов, определение уклонов склона, создание рельефа с отмывкой, вычисления кривизны поверхности, математические и статистические расчеты. Перечисленные выше операции использовались при определении ожидаемых оседаний земной поверхности от влияния горных работ. Полученная грид-тема может быть представлена в виде трехмерной поверхности или двух мерной в виде изолиний (рис.
Рис. 2. Оседание поверхности от 3х лав при средних геометрических размерах
1). Нормативными документами рекомендуется [1] при расчете параметров ожидаемых оседаний земной поверхности от влияния нескольких лав при мощности целика меньше 0,1 глубины, выполнять как от одной лавы, принимая ос-редненные значения длины лавы и выемочного столба, глубины.
Нормативными документами рекомендуется [1] при расчете параметров ожидаемых оседаний земной поверхности от влияния нескольких лав при
Среднегодо-
вая
Рис. 3. Оседание поверхности представленные в изолиниях от влияния 2-х лав, одна из которых имеет переменную длину представленная в MapInfo
мощности целика меньше 0,1 глубины, выполнять как от одной лавы, принимая осред-ненные значения длины лавы и выемочного столба, глубины. Однако лавы могут иметь сложную форму, поэтому расчет по средним значениям вносит ошибки и неточности. Демонстрация этого показана на рис. 2, где представлен фрагмент построенной мульды сдвижения в изолиниях, от влияния трех лав по средним значениям выемочного столба, при этом граница мульды сдвижения имеет меньший размер, нежели площадь, которая была отработана или больше.
Для преодоления этих неточностей расчет ожидаемых оседаний производим для каждой отдельной лавы с построением гридов, далее складывая их, получаем результирующий. Последний через обменный формат передается в МарМо (А^оСАБ или любой ГИС пакет), где полученный результат накладывается на планы горных работ.
1. Правила охраны сооружений и природных объектов от вредного влияния подземных горных
Результат описанного алгоритма выше представлен на рис. 3.
Путем вычитания гридов поверхности и ожидаемых оседаний можно получить новую поверхность, которая отразит результат влияния ведения горных работ на земную поверхность.
Таким образом, описанный выше алгоритм с использованием информационных систем позволяет довольно быстро и наглядно оценить степень влияния подземных горных разработок на земную поверхность.
-------------- СПИСОК ЛИТЕРАТУРЫ
разработок на угольных месторождениях. - М.,
Недра, 1996,288 с. ЕШ
г Коротко об авторе
Писаренко М.В. - кандидат технических наук, доцент, старший научный сотрудник Института угля и углехимии СО РАН, iuu@kemsc.ru
