You are currently viewing Программное обеспечение мониторинга достоверности запасов угледобывающих предприятий

Программное обеспечение мониторинга достоверности запасов угледобывающих предприятий

Содержание

Действующая отечественная классификация запасов [1], а также разрабатываемый «Российский Кодекс публичной отчетности о результатах геологоразведочных работ, ресурсах и запасах твердых полезных ископаемых», относящийся к кодексам семейства CRIRSCO, предполагают необходимость использования при оценке запасов специальных количественных методов.

В связи с этим авторами в содружестве с группой работников ФГУ «Государственная комиссия по запасам полезных ископаемых» (ГКЗ) были разработаны «Методические рекомендации по проведению количественной оценки степени соответствия геологических моделей месторождения угля его истинному состоянию» [2].

Данные рекомендации прошли соответствующее рассмотрение и даны к практическому применению экспертно-техническим советом ФГУ «ГКЗ» (протокол от 22.05.2007) и признаны в качестве официального нормативно-методического документа, рекомендованного «Обществом экспертов России по недропользованию» (протокол экспертно-технического совета ОЭРН от 12.05. 2011).

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

Алгоритмы расчета коэффициентов и исходные данные

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

Блок схема алгоритма расчета коэффициентов пропорциональности рисунок 1:

Рисунок 1 – Общая схема порядка ведения мониторинга достоверности запасов
Рисунок 1 – Общая схема порядка ведения мониторинга достоверности запасов

Основные этапы мониторинга согласно алгоритмам расчетов:

  1. Подготовка геологоразведочной информации по планам подсчета запасов и геологическим разрезам. Исходные данные для расчетов:
    • Номер скважины.
    • Плоские геодезические координаты точек пластоподсечения (X, Y, м).
    • Высотная отметка пласта в точке пластоподсечения (Z, м).
    • Дирекционный угол направления линии падения пласта в точке замера (А, градус).
    • Угол падения пласта в точке замера (δ, °).
    • Данные о значениях признаков, характеризующих форму, условия залегания и качественные свойства пласта (pn, например, мощность угольных пачек my, зольность пачек Ay, и т. п.).
  2. Квадриангулирование разведочной сети (разбиение на оценочные четырехугольные блоки с вершинами — разведочными скважинами) и расчет критериев разведанности гипсометрии (λ) и принятых к оценке геометризируемых признаков (Δi).
  3. Накопление по мере отработки пласта информации по данным горных работ. Исходные данные для расчетов:
    • Номер замера в горной выработке.
    • Плоские геодезические координаты точек измерения (X , Г, м).
    • Высотная отметка пласта в точке замера (Z, м).
    • Данные о значениях признаков, по которым построены горно-геометрические модели.
  4. Расчет фактических погрешностей (RF) моделей гипсометрии и признаков путем сравнения данных геологоразведочных и горных работ в оценочных блоках (при отработке участка пласта площадью в 20 и более оценочных блоков).
  5. Корректировка коэффициентов перехода (Kg, Кp) от критериев разведанности (λ, Δ) к погрешностям моделей (R).
  6. Построение картограмм разведанности горно-геометрических моделей признаков.

Для автоматизации работ по мониторингу достоверности запасов разработано специальное программное обеспечение, реализующее разработанные авторами алгоритмы, — программа «MDZ», правомочность применения которой подтверждена протоколами экспертно-технических советов ФГУ «ГКЗ» и ОЭРН.

Данная программа обеспечивает выполнение работ по оценке достоверности запасов как при отсутствии, так и при наличии горных работ (в режиме мониторинга) и имеет дружественный интерфейс, ориентированный на непрофессионального пользователя.

При старте программы пользователю предлагается основное меню, состоящее из 7 режимов:

  • «Данные».
  • «Расчет критериев».
  • «Адаптация».
  • «Установки».
  • «Сервис».
  • «Информация».
  • «Выход».

Работа программы начинается с режима «Установки», где указываются наименования рассматриваемых признаков обязательный признак «гипсометрия» и 8 дополнительных признаков:

  • Горно-геометрические модели для оценки.
  • Число отображаемых знаков в их значениях.
  • Предельные значения погрешностей для различных категорий запасов.
  • Коэффициенты перехода от критериев разведанности к погрешностям.

Для всех дополнительных признаков, достоверность моделей которых оценивается с помощью дельта-критериев разведанности, указывается вид используемого для классификации критерия:

  • Относительный.
  • Абсолютный.

Здесь же выбирается уровень надежности оценок (0,68 пли 0,80).

Режим «Данные» имеет два подрежима:

  • «Данные геологоразведки».
  • «Данные горных работ».

Эти подрежимы обеспечивают, соответственно, ввод данных по геологоразведочным скважинам и по данным замеров в горных выработках. Данные горных работ используются только при выполнении оценки достоверности в режиме мониторинга. Ввод данных соответствует общим правилам и настройкам Windows.

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

  • «Выбор признака».
  • «Оценка геологоразведки».

После выбора признака запускается подрежим «Оценка геологоразведки».

Работа с расчетно-графической средой программы

При старте подрежима на панели в масштабе показывается положение точек пластоподсечений. Пользователь имеет возможность подписать номера, увеличить или уменьшить изображение.

Разделение сети на систему выпуклых четырехугольников с вершинами — точками пластоподсечений выполняется интерактивно в режиме «Квадриангулирование».

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

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

Рисунок 2 – Фрагмент панели подрежима «Оценка геологоразведки» при выполнении квадриангулирования
Рисунок 2 – Фрагмент панели подрежима «Оценка геологоразведки» при выполнении квадриангулирования

Управление режимом «Квадриангулирование» осуществляется кнопками панели:

  • Кнопка «Просмотр блоков» обеспечивает возможность удаления любого из ранее выделенных четырехугольников.
  • Нажатие кнопки «Удаление блоков» приводит к удалению всех ранее выделенных блоков.

После завершения квадриангулирования осуществляется расчет критериев разведанности:

  • Для этого нажимается кнопка «Расчет критериев».
  • После на экране появляется панель результата расчета критериев (при оценке гипсометрии используются ламбда-критерии, при оценке остальных признаков — дельт а-критерии) в табличной форме.

Таблица расчета содержит колонки с порядковым номером оценочного четырехугольного блока в частности:

  • Номера скважин — его вершин.
  • Значения критериев разведанности (для гипсометрии ламбда-критериев, в вертикальном и в нормальном пласту направлении, для остальных признаков — абсолютный и относительный дельта-критерии).
  • Ожидаемая погрешность модели.
  • Площадь проекции оценочного блока в тысячах квадратных метрах.
  • Рекомендуемая категория запасов.

Погрешность модели гипсометрии выражается только в метрах, моделей остальных признаков в абсолютных или в относительных единицах (в зависимости от размерности, указанной пользователем в режиме).

При «нажатии» кнопки «Категоризация» выполняется построение картограммы разведанности, квадриангулированная часть сети скважин окрашивается соответствующими различным категориям запасов цветами в зависимости от величины ожидаемых погрешностей (рисунок 3):

Рисунок 3 – Панель построения картограммы разведанности
Рисунок 3 – Панель построения картограммы разведанности

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

  • Для реализации предусмотренной мониторингом процедуры корректировки коэффициентов пересчета критериев разведанности в погрешности признаков используется режим «Адаптация».
  • Собственно установление коэффициентов пересчета критериев разведанности в погрешности признаков (ранее упомянутых коэффициентов Kg и Кp) осуществляется в подрежиме «Коэффициенты перехода».
  • При старте режима появляется панель, на которой показывается положение пластоподсечений (в виде окружностей) и замеров в горных выработках (в форме «косых» крестов). В верхней части панели указывается наименование признака.

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

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

  • «Все последовательно»
  • «Необработанные — последовательно»
  • «Произвольный порядок»

Контур выбранного блока показывается на графической части панели, его номер и значение критерия разведанности указываются в окнах, расположенных под кнопкой «Выбор обрабатываемого блока» (рисунок 4):

Рисунок 4 – Фрагмент панели подрежима «Коэффициенты перехода» в режиме выбора относящихся к блоку замеров в горных выработках
Рисунок 4 – Фрагмент панели подрежима «Коэффициенты перехода» в режиме выбора относящихся к блоку замеров в горных выработках

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

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

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

Включение и исключение каждого замера из обработки сопровождается выполнением расчета среднеквадратической погрешности.

После завершения выбора замеров «нажимается» кнопка «Расчет по блоку» — на экране возникает панель расчета (рисунок 5):

Рисунок 5 – Фрагмент панели подрежима «Коэффициенты перехода» е режиме выбора расчета погрешности по блоку
Рисунок 5 – Фрагмент панели подрежима «Коэффициенты перехода» е режиме выбора расчета погрешности по блоку

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

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

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

  • В кровли.
  • На почве.

При необходимости удаление замера из обработки может быть осуществлено непосредственно из таблицы. Дня этого ячейка, подлежащая удалению, выделяется (например, ячейка замера 291 на рисунок 5), нажимается правая кнопка мыши — на экране появляется контекстное меню, содержащее предложение по удалению замера, после чего сопоставление может быть игнорировано (удалено).

Если пользователь удовлетворен результатами расчета, то нажимается кнопка «Принять для обработки», и результаты будут сохранены для дальнейшего использования. В противном случае «нажимается» кнопка «Не принимать в обработку», либо кнопка «Вернуться» (возврат к выбору замеров, относящихся к рассматриваемому блоку). «Нажатие» любой из трех перечисленных кнопок закрывает панель обработки.

Нажатие кнопки «Выбор обрабатываемого блока» приведет к рассмотрению следующего оценочного четырехугольника. По завершению обработки группы блоков (в том числе и при перерыве в работе) результаты можно сохранить в файле.

Обработка результатов расчета

После завершения обработки блоков нажатием кнопки «Обработка результатов» переходят собственно к расчету коэффициентов Kg или Кp.

При старте подрежима появляется панель, в таблице которой приведены критерии разведанности и фактические погрешности признака по оценочным блокам (рисунок 6):

Рисунок 6 – Панель подрежима «Расчет коэффициентов пересчета» при установлении величины коэффициентов пересчета критериев в погрешности
Рисунок 6 – Панель подрежима «Расчет коэффициентов пересчета» при установлении величины коэффициентов пересчета критериев в погрешности

Для поиска оптимального значения коэффициента пересчета в верхнюю строку таблицы необходимо:

  • Ввести несколько возможных вариантов (от 1 до 12 вариантов) значений коэффициентов пересчета (например: 1,0; 1,2; 0,7; 0,9; 1,4 — рисунок 6).
  • После ввода любого количества вариантов нажимается кнопка «Расчет таблицы». При этом в графической части панели показывается график зависимости вероятности от коэффициента Kg или Кp. Этот график отображается в виде ломаной прямой, соединяющей точки (рис. 6), или аппроксимирующей параболы.
  • Выбор варианта графика производится с помощью переключателя на варианты «Парабола» и «Ломаная». Для смены варианта построения графика необходимо изменить его с помощью переключателя и «нажать» «Расчет таблицы».
  • После построения графика с помощью переключателя устанавливается вероятность (0,68 или 0,80), для которой рассчитывается коэффициент.
  • Далее пользователь перемещает курсор в поле графика, нажимает и удерживает левую клавишу мыши (возникает изображение вертикальной стрелки). Не отпуская клавишу, стрелку перемещают по графику влево или вправо до тех пор, пока ее положение не совпадет с точкой пересечении графика и пунктирной линии вероятности (0,68 рисунок 6).
  • При этом перемещении соответствующее значение коэффициента отражается в окне «Коэффициент». Отпускание левой клавиши мыши фиксирует положение стрелки. «Нажатие» клавиши «Записать коэффициент» приводил к передаче значения коэффициента в блок установок режима «Установки».

Полученные в итоге указанных действий значения коэффициентов пересчета Kg или Кp будут использоваться программой при выполнении категоризации запасов в режиме «Оценка геологоразведки».

Нажатие кнопки панели «Передать в Excel» в любом режиме приводит к автоматическому открытию редактора Excel, в который будет экспортирована расчетная таблица (программа «MDZ» закрыта не будет). Кнопка «Рисунок в буфер» приводит к передаче изображения графика в стандартный буфер обмена, из которого его можно внедрить в документы текстовых, табличных и графических редакторов.

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

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

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

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

Список литературы

  1. Классификация запасов и прогнозных ресурсов твердых полезных ископаемых: приказ МПР РФ от 11.12.2006 № 278. — М., 2006. — 6 с.
  2. Методические рекомендации по проведению количественной оценки степени соответствия геологических моделей месторождения угля его истинному состоянию /Т.Б. Рогова, О.П. Никифорова, С.В. Шаклеин, В.А. Коткин, В.Г. Малухин, В.О. Ярков.-Кемерово: ОЭРН, 2011. — 86 с.

Источник: Программное обеспечение мониторинга достоверности запасов угледобывающих предприятий / Т.Б. Рогова, С.В. Шаклеин // Вестник КузГТУ. — 2012. — №1. — C. 20-26.

Статья в формате docx

Добавить комментарий

Gekoms LLC

Коллектив экспертов большая часть опыта и знаний которых востребованы в области промышленной автоматизации, разработке технически сложного оборудования, программировании АСУТП, управлении электроприводом. Телефон: +7(812) 317-00-87 Email: info@gekoms.com Сайт: https://gekoms.org