CC BY
50
УДК 622.142.5:004.9
Т. Б. Рогова, С. В. Шаклеин
ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ МОНИТОРИНГА ДОСТОВЕРНОСТИ ЗАПАСОВ УГЛЕДОБЫВАЮЩИХ ПРЕДПРИЯТИЙ
Действующая отечественная классификация запасов [1], а также разрабатываемый «Российский Кодекс публичной отчетности о результатах геологоразведочных работ, ресурсах и запасах твердых полезных ископаемых», относящийся к кодексам семейства СКЖБСО, предполагают необходимость использования при оценке запасов специальных количественных методов.
В связи с этим авторами в содружестве с группой работников ФГУ «Государственная комиссия по запасам полезных ископаемых» (ГКЗ) были разработаны «Методические рекомендации по проведению количественной оценки степени соответствия геологических моделей месторождения угля его истинному состоянию» [2].
Данные рекомендации прошли
соответствующее рассмотрение и даны к практическому применению экспертнотехническим советом ФГУ «ГКЗ» (протокол от 22.05.2007) и признаны в качестве официального нормативно-методического документа,
рекомендованного «Обществом экспертов России по недропользованию» (протокол экспертнотехнического совета ОЭРН от 12.05. 2011). Здесь рассмотрена специфика количественной оценки достоверности запасов действующих
угледобывающих предприятий с использованием технологии мониторинга достоверности запасов.
Мониторинг достоверности запасов,
выполняемый на участке месторождения, предполагает следующие этапы, обеспечивающие установление значений коэффициентов пропорциональности между степенью
неоднозначности горно-геомет-ри ческой модели и ее погрешностью по уровню, которой и классифицируются запасы (рис. 1).
Основные этапы мониторинга -
I. Подготовка геологоразведочной
информации по планам подсчета запасов и геологическим разрезам.
К числу исходных данных относятся:
- номер скважины;
- плоские геодезические координаты точек пластоподсечения (Х, У, м)
- высотная отметка пласта в точке пластоподсечения (2, м);
- дирекционный угол направления линии падения пласта в точке замера (А, градус);
- угол падения пласта в точке замера (5, °);
- данные о значениях признаков, характеризующих форму, условия залегания и качественные свойства пласта (рп, например,
мощность угольных пачек ту, зольность пачек Ау, и т. п.).
II. Квадриангулирование разведочной сети (разбиение на оценочные четырехугольные блоки с вершинами - разведочными скважинами) и расчет критериев разведанности гипсометрии (X) и принятых к оценке геометризируемых признаков (Д,).
III. Накопление по мере отработки пласта информации по данным горных работ:
К числу исходных данных относятся:
- номер замера в горной выработке;
- плоские геодезические координаты точек измерения (Х, У, м);
- высотная отметка пласта в точке замера (2,
м);
- данные о значениях признаков, по которым построены горно-геометрические модели.
IV. Расчет фактических погрешностей (ЯР) моделей гипсометрии и признаков путем сравнения данных геологоразведочных и горных работ в оценочных блоках (при отработке участка пласта площадью в 20 и более оценочных блоков).
V. Корректировка коэффициентов перехода (ККр) от критериев разведанности (X, Д) к погрешностям моделей (Я).
VI. Построение картограмм разведанности горно-геометрических моделей признаков.
Для автоматизации работ по мониторингу достоверности запасов разработано специальное программное обеспечение, реализующее разработанные авторами алгоритмы, - программа «МЭ2», правомочность применения которой подтверждена протоколами экспертно-
технических советов ФГУ «ГКЗ» и ОЭРН.
Данная программа обеспечивает выполнение работ по оценке достоверности запасов как при отсутствии, так и при наличии горных работ (в режиме мониторинга) и имеет дружественный интерфейс, ориентированный на
непрофессионального пользователя.
При старте программы пользователю предлагается основное меню, состоящее из 7 режимов: «Данные»; «Расчет критериев»;
«Адаптация»; «Установки»; «Сервис»;
«Информация»; «Выход».
Работа программы начинается с режима «Установки», где указываются наименования рассматриваемых признаков (обязательный признак «гипсометрия» и 8 дополнительных признаков, горно-геометрические модели которых могут оцениваться), число отображаемых знаков в их значениях, предельные значения погрешностей
Рис. 1. Общая схема порядка ведения мониторинга достоверности запасов
для различных категорий запасов, а также коэффициенты перехода от критериев разведанности к погрешностям. Для всех дополнительных признаков, достоверность моделей которых оценивается с помощью дельта-критериев разведанности, указывается вид используемого для классификации критерия: относительный или
абсолютный. Здесь же выбирается уровень надежности оценок (0,68 или 0,80).
Режим «Данные» имеет два подрежима: «Данные геологоразведки» и «Данные горных работ», которые обеспечивают, соответственно, ввод данных по геологоразведочным скважинам и по данным замеров в горных выработках. Данные горных работ используются только при выполнении оценки достоверности в режиме мониторинга. Ввод данных соответствует общим
правилам и настройкам Windows.
Для расчета критериев разведанности используется режим «Расчет критериев» основного меню. Данный режим состоит из двух подрежимов: «Выбор признака» и «Оценка
геологоразведки».
После выбора признака запускается подрежим «Оценка геологоразведки». При старте подрежима на панели в масштабе показывается положение точек пластоподсечений. Пользователь имеет возможность подписать номера, увеличить или уменьшить изображение.
Разделение сети на систему выпуклых четырехугольников с вершинами - точками пластоподсечений выполняется интерактивно в режиме «Квадриангулирование». Для этого пользователь должен указать маркером мыши
2072
Все номера
Удаляъ номера
К&адриак-.ч»»<>овамие
3-я скважина 2072
Просмотр б попов
Удаление блоков
Г" Покаэьевгьдагныепрм просмотре блоков ____________Расчет критериев___________|
Рис. 2. Фрагмент панели подрежима «Оценка геологоразведки» при выполнении квадриангулирования
первую скважину четырехугольника и нажать левую клавишу мыши (в качестве первой скважины может приниматься любая скважина сети). Затем, обходя четырехугольник по или против направления движения часовой стрелки, пользователь нажатием левой клавиши мыши выбирает вторую, третью (как на рис. 2) и четвертую скважины четырехугольника.
Управление режимом «Квадриангулирование» осуществляется кнопками панели. Кнопка «Просмотр блоков» обеспечивает возможность удаления любого из ранее выделенных
четырехугольников. Нажатие кнопки «Удаление блоков» приводит к удалению всех ранее выделенных блоков.
После завершения квадриангулирования
осуществляется расчет критериев разведанности. Для этого нажимается кнопка «Расчет критериев», и на экране появляется панель результата расчета критериев (при оценке гипсометрии используются ламбда-критерии, при оценке остальных
признаков - дельта-критерии) в табличной форме.
Таблица расчета содержит колонки с порядковым номером оценочного
четырехугольного блока, номера скважин - его вершин, значения критериев разведанности (для гипсометрии ламбда-критериев, в вертикальном и в нормальном пласту направлении, для остальных признаков - абсолютный и относительный дельтакритерии), ожидаемая погрешность модели, площадь проекции оценочного блока в тысячах квадратных метрах и рекомендуемая категория
запасов.
Погрешность модели гипсометрии выражается только в метрах, моделей остальных признаков в абсолютных или в относительных единицах (в зависимости от размерности, указанной пользователем в режиме
При «нажатии» кнопки «Категоризация» выполняется построение картограммы
разведанности, квадриангулированная часть сети скважин окрашивается соответствующими различным категориям запасов цветами в зависимости от величины ожидаемых погрешностей (рис. 3).
Построение картограммы разведанности является итоговой процедурой оценки достоверности запасов, выполняемой как на основе геологоразведки, так и в режиме мониторинга.
Для реализации предусмотренной мониторингом процедуры корректировки
коэффициентов пересчета критериев
разведанности в погрешности признаков используется режим «Адаптация». Собственно установление коэффициентов пересчета критериев разведанности в погрешности признаков (ранее упомянутых коэффициентов К и Кр) осуществляется в подрежиме «Коэффициенты перехода».
При старте режима появляется панель, на которой показывается положение
пластоподсечений (в виде окружностей) и замеров в горных выработках (в форме «косых» крестов).
Резуъты расчета г» бпок« № 7 (2072 ■ 1620 • 656 • 33801. критерий 16.8 Признак Гюсометрия
Замер Данхл разведки | Данные разработки Разность -<1 и а [ЧС\2 ±1
341 2048 1973 7.5 7.5 56 25
342 203 3 2020 1.3 1.3 1 69
343 1950 1938 1.2 1.2 1.44
229 1683 1682 01 01 0.01
339 211.3 211 1 02 02 0 04
340 208 4 207 8 06 06 0 36
344 2138 2139 -01 01 001
286 214 2 211 7 25 25 625
195 1192 1196 ■04 ■04 016
291 245 8 _ —! у> < •> 11 с 115 13225
Игнорировать соооставле**1в 291 Сохранить сопоставление
345 1896 05 025
Су<ма 5081 9 400 296 92
Среднее 195 5 1939 1.5 1.5 11.5
Средне* маратич 34 1
3
Систематическая погрешность (С) отсутствует Передать в ЕхсН Вернуться
Не пр»*с*.»ать в обработку
Принять для обработки
Рис. 5. Фрагмент панели подрежима «Коэффициенты перехода» в режиме выбора расчета погрешности по блоку
В верхней части панели указывается наименование признака.
Первая задача работы в подрежиме «Коэффициент перехода» состоит в определении истинной погрешности горно-геометрической модели в пределах контура оценочного четырехугольного блока. Она решается путем расчета среднеквадратического отклонения данных горных и геологоразведочных работ в точках замеров, относящихся к рассматриваемому оценочному блоку.
Предусмотрено три варианта порядка рассмотрения оценочных блоков,
устанавливаемых с помощью переключателей: «Все последовательно», «Необработанные -последовательно» «Произвольный порядок».
Контур выбранного блока показывается на графической части панели, его номер и значение критерия разведанности указываются в окнах, расположенных под кнопкой «Выбор обрабатываемого блока» (рис. 4).
Замеры внутри блока принимаются для оценки погрешности, указанной в специальном окне. Признаком учета замера является красный
цвет его условного обозначения. При необходимости пользователь, устанавливая курсор мыши на замер и нажимая ее левую клавишу, может выбрать дополнительные, окружающие блок замеры, которые будут учтены при расчете погрешности (при этом цвет их знака меняется на красный).
Повторный выбор замера приводит к его исключению из расчетов (цвет условного знака сменится на черный). Целесообразность исключения замера, находящегося внутри контура оценочного четырехугольника, может быть вызвана необходимостью обеспечения
равномерности размещения замеров по его площади.
Включение и исключение каждого замера из обработки сопровождается выполнением расчета среднеквадратической погрешности.
После завершения выбора замеров «нажимается» кнопка «Расчет по блоку» - на экране возникает панель расчета (рис. 5). На панели указываются номер блока, наименование его вершин (скважин), значение критерия разведанности, имя обрабатываемого признака и
Расчет коэффициентов пересчета
Признак Гипсометрия
№ 2ее(*і»*в Змримм 4(«[шкм К{мтефмй Факт псгреиюсть-Г 0.7 0.9 1.0 1.2 1.4
1 9754 9750 9752 26 24 18 23 26 31 36 □
2 9754 3245 9750 27 23 19 24 27 32 38
3 3245 3380 9750 58 1.7 41 52 58 70 81
4 9750 3258 3380 01 20 01 01 01 01 01
5 3245 9750 1 9 20 1.3 1 7 1 9 23 27
І Б56 2072 3380 138 1.7 ¡317 124 138 166 19 3
7 1620 656 3380 168 3.1 11 8 151 168 202 235
6 3258 2072 3380 27 21 1 9 24 27 32 38 _ |
«Г с **>ел „ J 1ЛЛЛ Л 1 л *> * * с т Г «V 1 л А А 2І
Рис. 6. Панель подрежима «Расчет коэффициентов пересчета» при установлении величины коэффициентов пересчета критериев в погрешности
величина учтенной при расчетах систематической погрешности.
В находящейся на панели таблице указаны номера использованных замеров в горных выработках, значения признака в них и результаты промежуточных расчетов. В качестве указанной в последней строке таблицы погрешности горногеометрической модели принимается
среднеквадратическая разность значений признака по данным горных и геологоразведочных работ в замерах, отнесенных к оценочному блоку.
Приведенные в таблице результаты сопоставлений рекомендуется анализировать на предмет выявления ошибочных и нехарактерных результатов, появление которых связано с использованием в качестве замеров точек маркшейдерской сети не всегда однотипно расположенных по отношению к пласту: например, в кровли и на почве.
При необходимости удаление замера из обработки может быть осуществлено непосредственно из таблицы. Для этого ячейка, подлежащая удалению, выделяется (например, ячейка замера 291 на рис. 5), нажимается правая кнопка мыши - на экране появляется контекстное меню, содержащее предложение по удалению замера, после чего сопоставление может быть игнорировано (удалено).
Если пользователь удовлетворен результатами расчета, то нажимается кнопка «Принять для обработки», и результаты будут сохранены для дальнейшего использования. В противном случае «нажимается» кнопка «Не принимать в обработку», либо кнопка «Вернуться» (возврат к выбору замеров, относящихся к
рассматриваемому блоку). «Нажатие» любой из трех перечисленных кнопок закрывает панель обработки.
Нажатие кнопки «Выбор обрабатываемого блока» приведет к рассмотрению следующего оценочного четырехугольника. По завершению обработки группы блоков (в том числе и при перерыве в работе) результаты можно сохранить в файле.
После завершения обработки блоков нажатием кнопки «Обработка результатов» переходят собственно к расчету коэффициентов
К илиКр .
При старте подрежима появляется панель, в таблице которой приведены критерии
разведанности и фактические погрешности признака по оценочным блокам (рис. 6).
Для поиска оптимального значения
коэффициента пересчета в верхнюю строку таблицы необходимо ввести несколько возможных вариантов (от 1 до 12 вариантов) значений коэффициентов пересчета (например: 1,0; 1,2; 0,7; 0,9; 1,4 - рис. 6). После ввода любого количества вариантов нажимается кнопка «Расчет таблицы». При этом в графической части панели показывается график
зависимости вероятности от коэффициента Kg или Kp. Этот график отображается в виде ломаной прямой, соединяющей точки (рис. 6), или аппроксимирующей параболы.
Выбор варианта графика производится с помощью переключателя на варианты «Парабола» и «Ломаная». Для смены варианта построения графика необходимо изменить его с помощью переключателя и «нажать» «Расчет таблицы».
После построения графика с помощью переключателя устанавливается вероятность (0,68 или 0,80), для которой рассчитывается коэффициент.
Далее пользователь перемещает курсор в поле графика, нажимает и удерживает левую клавишу мыши (возникает изображение вертикальной стрелки). Не отпуская клавишу, стрелку перемещают по графику влево или вправо до тех пор, пока ее положение не совпадет с точкой пересечения графика и пунктирной линии вероятности (0,68 рис. 6). При этом перемещении соответствующее значение
коэффициентаотражается в окне «Коэффициент». Отпускание левой клавиши мыши фиксирует положение стрелки.
«Нажатие» клавиши «Записать коэффициент» приводит к передаче значения коэффициента в блок установок режима «Установки».
Полученные в итоге указанных действий значения коэффициентов пересчета Kg или Kp будут использоваться программой при выполнении категоризации запасов в режиме «Оценка геологоразведки».
Нажатие кнопки панели «Передать в Excel» в любом режиме приводит к автоматическому открытию редактора Excel, в который будет
экспортирована расчетная таблица (программа «MDZ» закрыта не будет). Кнопка «Рисунок в буфер» приводит к передаче изображения графика в стандартный буфер обмена, из которого его можно внедрить в документы текстовых, табличных и графических редакторов.
Помимо обеспечения автоматическими расчетами и графическими построениями основных этапов мониторинга достоверности запасов, программа «MDZ» содержит вспомогательные процедуры, позволяющие
выполнять:
- построение и анализ кривой разведанности,
используемой для оценки правомерности
интерполяции высотных отметок пласта;
- оценку правомерности интерполяции
признаков, модели которых отстраиваются путем линейной интерполяции методом многогранника, и расчет погрешности отстраиваемых изолиний;
- установление наличия и величины
систематической погрешности данных
геологоразведочных работ и ее исключения;
- расчеты элементов залегания пласта в точке замера, расположенной в замке складке;
- оценку геометрии формы оценочного блока;
- исключать из расчетов оценочные блоки с неоправданно низкими значениями критериев, либо заменять значения их критериев на заданное пороговое значение.
Представление результатов теоретических исследований в форме реализованных в программе алгоритмов потребовало от авторов значительных затрат труда (какие-либо сторонние специалисты к разработке программы не привлекались). Но авторы не сожалеют о затраченном на это времени, ибо эта реализация позволила выявить и устранить все возможные неопределенности действий и, в ряде случаев, потребовало их научного решения. Наличие компьютерного обеспечения переводит технологию мониторинга достоверности запасов из научно-исследовательской стадии в инженерную.
СПИСОК ЛИТЕРАТУРЫ
1. Классификация запасов и прогнозных ресурсов твердых полезных ископаемых: приказ МПР РФ от
11.12.2006 № 278. - М., 2006. - 6 с.
2. Рогова, Т. Б. Методические рекомендации по проведению количественной оценки степени соответствия геологических моделей месторождения угля его истинному состоянию /Т. Б. Рогова, О.П. Никифорова, С.В. Шаклеин, В.А. Коткин, В.Г. Малухин, В.О. Ярков.-Кемерово: ОЭРН, 2011. - 86 с.
□ Авторы статьи:
Шаклеин Сергей Васильевич, докт. техн. наук, вед. научн. сотр. филиала ИВТ Со РАН, проф. каф. маркшейдерского дела, кадастра и геодезии КузГТУ E-mail: [email protected].
Рогова Тамара Борисовна, канд. техн. наук, доцент каф. маркшейдерского дела, кадастра и геодезии КузГТУ. Тел. 906-928-71-35. E-mail: [email protected]....
