CC BY
62
ИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ
УДК 622.142.5:004.42
С.В. Шаклеин, В.В. Гетман
АЛГОРИТМ ОЦЕНКИ АНИЗОТРОПИИ МОЩНОСТИ УГОЛЬНОГО ПЛАСТА И ЕГО ПРОГРАММНАЯ РЕАЛИЗАЦИЯ
В связи с вступлением в силу поправки в Налоговый кодекс, предусматривающей взимание налога за добычу угля с каждой извлеченной их недр тонны по фиксированной ставке, коренным образом изменился правовой статус маркшейдерского замера добычи по чистым угольным пачкам. Ранее результаты замера использовались при осуществлении государственного учета движения запасов и решения ряда корпоративных задач.
Новый статуса замера требует использования жесткой методики маркшейдерских измерений мощности пласта в контуре отрабатываемого выемочного столба, обеспечивающей достоверность данных, необходимых для расчета платежей за
тикальная) анизотропия для них не имеет принципиального значения
Собственно анизотропию мощности в контуре выемочного столба возможно оценить только геометрическими методами, поскольку использование геостатистических подходов невозможно, т.к. предполагает необходимость исследования характера вариограмм не менее чем по четырем различным направлениям [2]. Вместе с тем, следует обратить внимание на геостатистическую идею описания анизотропии с помощью эллипса анизотропии [2], которая, в целом, согласуется с отечественными горно-геометрическими представления-
ми.
пользование недрами при добыче полезных ископаемых. Характер сети измерений мощности, определяемый системой горных выработок, предполагает прямоугольную ее конструкцию. Соотношение длин ее сторон явно должно соответствовать параметрам анизотропии геополя мощности пласта.
В связи с этим разработка методики оценки анизотропии мощности пласта в контуре выемочного столба является актуальной задачей, не имеющей в настоящее время общепризнанного решения.
На практике используют несколько простых моделей анизотропии: геометрическая и зональная. Для условий угольных месторождений вполне достаточно ограничится учетом только геометрической анизотропии, поскольку зональная (вер-
К числу наиболее разработанных методов геометрического подхода следует отнести предложение Л. И. Четверикова [1], который исходит из того, что всесторонняя характеристика анизотропии может быть осуществлена при помощи построения предлагаемой им индикатрисы анизотропии - воображаемой вспомогательной поверхности, выражающей изменчивость параметра по различным направлениям внутри объекта. Каждый радиус-вектор такой индикатрисы геометрически в некотором масштабе пропорционален значению координированной изменчивости признака в направлении этого радиуса-вектора. Собственно показателем изменчивости полагается средняя величина первых разностей между значениями признака по избранному направлению радиуса-вектора. В отличие от традиционных разностных
методов, Л. И. Четвериков предлагает использовать не данные в точках производства замеров, а полученные в результате геометризации признака его топографические поверхности представленные в виде системы изолиний.
Собственно методика построения индикатрисы анизотропии Л. И. Четверикова сводится к тому, что на анализируемый фрагмент топографической поверхности (рис. 1) «набрасывается» палетка - сетка параллельных линий (пунктирные линии на рис. 1), центр которой совпадает с положением «точки» измерения анизотропии (точкой О).
Для каждой линии палетки определяется ее длина іі, находящаяся в пределах контура оценки и подсчитывается количество интервалов между изолиниями на плане, пересекаемое всеми линиями палетки при данной ее ориентировке - величина п.
Полагая, что величина сечения изолиний К, при котором выполнено их построение постоянным, рассчитывается величина показателя изменчивости
I = К —. (1)
І іі
Величина этого показателя полагается соответствующей направлению ориентации линий палетки.
Затем все линии палетки разворачиваются на угол ф (рис. 1б) и определяется показатель изменчивости I, соответствующий данной ориентации палетки. После этого палетка вновь поворачивается на угол ф и вычисляется новое значение показателя изменчивости. Разворот палетки осуществляется до тех пор, пока общий угол ее поворота не составит 180°.
По завершению описанных действий, осуществляется построение индикатрисы анизотропии.
Основой для построения индикатрисы является серия равновеликих направлений, формирующих сектора, угол между которыми равен ф (рис. 2). Выбрав необходимый масштаб, вдоль каждого из них от центра индикатрисы О откладывается соответствующее направлению величина показателя изменчивости.
В результате этих действий на каждом направлении будут нанесены две точки. На рис. 2 точки А и А' соответствуют показателю изменчивости по направлению «юг - север», В и В' - направлению под углом ф. Собственно индикатриса представляет собой многоугольник, соединяющие все нанесенные на нее точки.
Для упрощения процедуры оценки показателя анизотропии I предлагается использовать квадратную палетку с длиной диагонали, равной наименьшему расстоянию между подготовительными штреками. При вращении такой палетки вокруг центра О все ее линии будут постоянно находится внутри контура выемочного столба.
Поэтому, при условии постоянства сечения изолиний, при ее применении показатель анизотропии Л. И. Четверикова, вычисляемый по реко-
мендуемой формуле (1) при любом угле ее поворота ф будет равен:
I = tn, (2)
K
где t =-----= const - отношение сечения изоли-
I h
ний к суммарной длине линий палетки; n - суммарное количество изолиний, пересекаемых линиями палетки.
Поскольку, при оценке анизотропии необходимо ориентироваться на соотношения значений ее показателей по различным направлениям, то, учитывая постоянство коэффициента пропорциональности t, в качестве эквивалентного показателя может быть использовано только число n пересечений линий палетки изолиниями.
Апробированные геостатистической практикой представления о геометрической (эллиптической) анизотропии наталкивают на очевидную мысль о том, что построенную индикатрису анизотропии следует аппроксимировать эллипсом, сгладив, тем самым, отдельные частные колебания оценок по направлениям. Такую аппроксимацию следует рассматривать в качестве метода математической обработки индикатрис.
Собственно аппроксимацию предлагается осуществлять с помощью метода наименьших квадратов.
На рис. 3 представлена индикатриса анизотропии и аппроксимирующий ее эллипс (рис. 3.2а), параметры которого характеризуются ориентировкой его главной оси а и длинами большой и малой полуосей - а и b (рис. 36).
В качестве условия определения трех характе-
ризующих аппроксимирующий эллипс параметров (угла а и полуосей а и Ь) предлагается условие: к 2 $ = X (П - ш 1 ) = шт , (3)
г=1
где и,- и ш{ - соответственно, длина г-го луча индикатрисы анизотропии (количество пересечений изолиний линиями палетки) и расстояние между центром и аппроксимирующим эллипсом по направлению г-го луча; к - количество лучей.
Отыскание значений параметров а, а и Ь путем аналитического решения уравнения (3) никем ранее не производилось и, безусловно, является достаточно сложной задачей. Учитывая строго дискретный характер используемых данных (постоянное число лучей индикатрисы и углов между ними) и возможности, представляемые современной вычислительной техникой, поставленная задача была решена численным методом.
Разработанный алгоритм ее решения состоит в следующем.
Во-первых, выделяются все возможные направления расположения большой оси аппроксимирующего эллипса. При угле между лучами индикатрисы в 30° выделяется шесть (рис. 4), а при 15° - 12 таких направлений.
Во-вторых, каждый из выделенных вариантов положения главной оси рассматривается индивидуально. Рассмотрение имеет целью определение конкретных значений длин полуосей а и Ь, для которых выполняется условие (3).
Поиск осуществляется методом перебора возможных значений полуосей. В качестве их максимально возможного значения принимается макси-
мальная наблюдаемая на индикатрисе величина п,. Минимально возможная длина полуосей принимается равной 1 (линии палетки пересекли только одну изолинию). Шаг перебора принимается равным 1.
Для каждого задаваемого в указанном диапазоне значения полуоси а рассматриваются все возможные величины полуоси Ь. Для каждой комбинации значений а и Ь по формуле (3) рассчитывается сумма квадратов отклонений $.
В ходе формирования вариантов возможных значений а и Ь возникают и такие, для которых Ь>а. Понятно, что они некорректны по своей сути и для них расчет суммы $ просто не производится.
Для каждого рассматриваемой величины а, зависимость получаемой суммы $ от Ь имеет характер функции, обладающей только одним экстремумом. Это создает возможность ограничения рассматриваемой совокупности полуосей Ь, т. е. сокращения объема перебора. Однако скорость современных рядовых средств вычислительной техники такова, что достигаемое сокращение времени расчетов совершенно не заметно, а, следовательно, и малозначимо.
Расчет величин шi предлагается выполнять по формуле:
ш =
Ь
2
где е - эксцентриситет эллипса.
В результате расчетов для каждого /-го варианта сочетаний возможных значений полуосей а и Ь для рассматриваемого варианта ориентации главной оси эллипса определяются значения суммы $. В качестве полуосей ориентированного под углом а эллипса выбираются те, при использовании которых сумма $ была минимальна.
И, наконец, проводя указанные вычисления
для всех возможных вариантов ориентации эллипса, из всей совокупности значений полуосей выбирается та, для которой сумма $ минимальна. Полученные таким образом полуоси и ориентация эллипса являются параметрами эллипса, аппроксимирующего индикатрису анизотропии.
В качестве показателей описанной соответствующим эллипсом геометрической анизотропии предлагается использовать следующие основные характеристики.
Во-первых, ориентацию главной оси эллипса а, которая характеризует направления наибольшей изменчивости признака. А, во-вторых, отношение длины малой полуоси эллипса Ь к большой полуоси а.
Это отношение является характеристикой различий интенсивности изменения значений признака по направлению наибольшей и наименьшей изменчивости.
Геометрически, это отношения равно косинусу угла, на который необходимо повернуть плоскость проекции, в которой изображен эллипс, вокруг его большой оси так чтобы в новой проекции он превратился в окружность. Т. е. преобразовать координат так, чтобы анизотропное поле «превратилось» в изотропное.
Поскольку именно эта операция используется в геостатистике, то по аналогии с применяемой в ней терминологии отношение полуосей будем именовать отношением анизотропии и обозначать как
к=Ь.
а
Применительно к решаемой задаче возникает необходимость введения еще одного показателя анизотропии, а именно соотношение анизотропии признака по заданным ортогональным направлениям, а именно по направлению расположения
подготовительных и очистной выработок.
Необходимость введения такого показателя можно проиллюстрировать рис. 5. На нем изображено два совершенно одинаковых эллипса анизотропии, при различных направлениях ведения горных работ. Проводя через центр эллипсов направления, параллельные направлениям подготовительных штреков и линии очистного забоя, можно определить расстояния характеризующие степень анизотропии вдоль этих направлений.
Из рис. 5а следует, что скорость изменения признака вдоль направления очистного забоя почти в два раза выше, чем по направлению штреков. Отсюда следует, что для получения равноточных результатов по двум указанным направлениям измерений, сеть замеров по лаве должна быть примерно в два раза более густой, чем по штрекам.
При взаимной ориентации эллипса и горных выработок, имеющей место в условиях рис. 56 наблюдается равновеликая скорость изменения признака, как в направлении очистного забоя, так и в направлении штреков. Следовательно, плотность сети замеров в обоих направлениях может быть одинаковой.
Исходя из изложенного, требуется ввести еще один относительный показатель, который, в конце концов, и должен определять отношение интерва-
лов измерений между замерами в двух взаимно перпендикулярных направлениях.
Величину этого показателя, который, в связи с его зависимостью от технологических решений может быть условно назван технологическим отношением анизотропии, предлагается определять по формуле:
V = LШ
^ ’
где Ьш и Ьл - проходящие через центр эллипса анизотропии расстояния между его токами, ориентированными, соответственно, вдоль подготовительных штреков и по линии очистного забоя.
Нетрудно заметить, что построение эллипса анизотропии является трудоемкой задачей, решение которой без применения средств вычислительной техники практически невозможно.
Разработанное программное обеспечение,
оформленное в виде прикладной программы «А№2», реализует описанный выше алгоритм и снабжено рядом сервисных элементов, повышающих удобство работы. Вид панели программы представлен на рис. 6.
В левой части панели расположено графическое окно с координатной системой индикатрисы анизотропии. В правой части - окна ввода тексто-
Рис. 6. Панель программы «АNIZ»
вой и числовой информации, а также управляющие кнопки.
Подготовка программы к работе состоит в заполнении общей, не имеющей обязательного характера информации, такой как наименовании шахты (окно «Шахта»), пласта («Пласт») и диаметра оценки (длины диагонали палетки Л. И. Четверикова). С помощью переключателя «Данные» указывается, по каким данным была построена горно-геометрическая модель изменения мощности пласта - по данным разведки или по данным горных работ.
С помощью окна «Шаблон имени файла» пользователь может занести любую фразу, с которой будет начинаться автоматически генерируемое имя файла, в который будут записываться все введенные данные.
Все перечисленные данные могут быть записаны на диск для дальнейшего использования путем «нажатия» кнопки «Сохранить установки». При этом автоматически будет создан файл «а^.орс», в котором общие данные будут сохранены. Файл формируется в том каталоге, в котором находится файл программы.
По завершению подготовки программы к работе в нее водятся исходные данные. Как это понятно из предыдущего текста, эти данные подготавливаются вручную с использованием плана горных работ, на котором показаны контуры выемочного столба и топографическая поверхность мощности пласта в изолиниях.
На план горных работ накладывается палетка Л. И. Четверикова таким образом, чтобы ее вертикальная ось совпадала с направлением на север. После этого подсчитывается суммарное количество пересечений изолиний всеми линиями палетки. Это целое число вводится в таблицу данных в колонку «Число изолиний» в строку соответствующую углу поворота палетки «0». Затем палетка разворачивается на 15° по направлению движения часовой стрелки и находится новое суммарное число пересечений, которое записывается в строку для угла 15° и т. д.
По заполнению данных о числе пересеченных
палеткой изолиний в окне «Простирание» указывается дирекционный угол (в принятой системе координат) направления подготовительных штреков в районе точки наблюдений.
Это направление может быть впоследствии показано на индикатрисе красной линией, заканчивающейся небольшой окружностью. Но ее построение производится только в случае, если на необходимость этого действия указывает наличие «птички» в окне «показать», размещенном под окном ввода дирекционного угла.
Помимо ввода данных по индикатрисе непосредственно в программе, возможен их импорт из табличного редактора «Excel». Для этого, сформированная в «Excel» колонка данных должна быть выделена и сохранена в буфере обмена.
После ввода данных пользователь получает возможность сохранить их на диске. Для этого служит кнопка «Записать в файл». Построение индикатрисы анизотропии производится программой после «нажатия» кнопки «Построить индикатрису анизотропии». После ее «нажатия» в графической части экрана появляется масштабное изображение индикатрисы с нанесением на нее аппроксимирующего эллипса и всех полученных результатов.
Экспорт полученного изображения в другие программные продукты осуществляется путем «нажатия» кнопки «Рисунок в буфер». Это обеспечивает помещение изображения в буфер обмена Microsoft Office, из которого его можно внедрить в документы текстовых, табличных и графических редакторов.
Реализованный программой «ANIZ» алгоритм позволяет выполнять все расчеты по оценке анизотропии мощности, необходимые для обоснования конструкции маркшейдерской сети измерения мощностей чистых угольных пачек пласта в контуре отрабатываемого выемочного столба представляемой органам Росприроднадзора, на которые возложены контрольно-надзорные функции за достоверностью данных для расчета платежей за пользование недрами при добыче полезных ископаемых.
СПИСОК ЛИТЕРАТУРЫ
1. Четвериков, Л. И. Оценка анизотропии наблюдаемой изменчивости параметров тел полезных ископаемых / Л. И. Четвериков // Известия вузов. Горный журнал. - 1972. - № 4. - С. 36-40.
2. Armstrong, M. Basic Linear Geostatistics / М. Armstrong. - Berlin: Springer Verlag, 1998. - 153 p.
□ Авторы статьи
Шаклеин Сергей Васильевич, докт. техн. наук, проф. каф. маркшейдерского дела, кадастра и геодезии КузГТУ, вед. научн. сотр. Кемеровского филиала ИВТ СО РАН. E-mail: svs1959@mail.ru
Гетман Валерий Валериевич, зам. технического директора по землепользованию и лицензированию ОАО «СУЭК-Кузбасс» E-mail: getmanvv@suek.ru
