Необходимым и решающим условием для обеспечения рентабельности отработки руды при общем низком уровне извлечения металла является создание геолого-технологических блочных моделей (ГТБМ), учитывающих не только содержание металла в руде, но и технологические свойства руды, среди которых наиболее важными являются коэффициент извлечение полезного компонента из руды, а также показатели физико-механических свойств руд. Преимущество ГТБМ по сравнению с моделью, построенной на основании содержания полезного компонента, связано с тем, что классификация руда/порода основывается на применении к блоку модели условия рентабельности в виде значения безубыточного бортового содержания.
Блоки представляют собой объемы горной массы, в которых оценены содержания по данным блочной модели. При этом также определена технология горных работ и переработки руды, позволяющая оценить наиболее вероятный диапазон затрат на горные работы и переработку руды. В общих чертах стратегия планирования состоит в том, что начальная стадия будет иметь высокий борт, а каждая последующая стадия — более низкий борт благодаря использованию рудных складов с учетом плавающего бортового содержания для текущих экономических условий и текущего этапа планирования. При этом порог разделения между рудой и породой определяется по безубыточному бортовому содержанию (ББС), обеспечивающему окупаемость затрат на добычу.
Формулы для аналитического расчета безубыточного бортового содержания при различных условиях добычи полезного ископаемого и его переработки является предметом исследований и дискуссий (за рамками компетенции автора статьи), поэтому в данной статье рассмотрим (на рис. 1) наиболее простую и показательную формулу, в которой учитываются основные показатели: а) З — затраты на добычу и переработку; б) Ц — цена металла; в) И — извлечение при переработке.
Рис. 1. Формула для аналитического расчета безубыточного бортового содержания и пример разделения между балансовой (розовая штриховка) и забалансовой рудой на фрагменте блочной модели, с указанием в каждом блоке: оцененных содержаний золота по данным блочной модели (слева); справа — рассчитанного безубыточного бортового содержания, (1,0 г/т)
Важно заметить, что несмотря на существующий в России приоритет оконтуривания и оценки запасов руды по утвержденному экспертизой ГКЗ бортовому содержанию, передовые горнодобывающие предприятия приходят к необходимости оконтуривания рентабельных участков недр через применение механизма плавающего бортового содержания по блочным моделям с учетом того, что цена металла и затраты могут изменяться в широких пределах по сравнению с усредненными значениями.
Первым шагом к построению геолого-технологической модели, является разработка программы отбора проб ГТК, что связано с необходимостью классификации интервалов опробования, которые являются неоднородными по природным признакам. Результатом и непременным условием разработки программы отбора является выделение принципов классификации руд на природные типы с целью их последующего объединения в композитные пробы с таким расчетом, чтобы в последующем определить рабочие зависимости между вещественными критериями (например, показателями химического и минералогического состава) и экспериментальными данными технологических исследований проб. Группировка частных проб проводится методом кластеризации интервалов опробования по принципу максимальной однородности интервалов частных проб в составе композитной пробы, путем определение классификационных параметров руд, применимых к частным пробам. Материал для каждой композитной пробы следует отбирать отдельно по каждому из выделенных природных типов.
Анализ геологических материалов позволяет составить лишь предварительную картину природной типизации руд. В ходе анализа данных геологической документации и рядового опробования на примере одного из проектов ГТК установлено, что наиболее информативной характеристикой для выделения природных типов руд являются показатели содержаний Au, Cu и Fe, а также степень окисления и литологическая характеристика. На основании гипотезы о том, что окисленные руды развиты в пределах зон интенсивной гипергенной минерализации, не удалось произвести предварительную типизацию из-за того, что в ходе проведения рядового опробования в рамках поисковых и оценочных работ в недостаточной степени уделено внимание выделению интервалов распространения гипергенных и вторичных минералов, а оцифрованная информация, содержащаяся в базе первичной геологической документации скважин и горных выработок, является малоинформативной и недостаточной для выделения минеральных разновидностей и природных типов руд. Лучшим решением в данной ситуации является возможность в результате построения нейросетевых моделей по данным химического анализа максимально достоверно рассчитать недостающие прогнозные значения степени окисления (там, где отсутствуют лабораторные данные) для рудных интервалов, намеченных для технологического опробования, что является вещественным фактором для природной типизации руд. Таким образом в специализированном пакете статистической обработки ST Neural Networks (SANN) был произведен Кластерный анализ обучение нейронной модели по сети Кохонена для поиска кластеров по двум типа переменных: а) Числовые входные переменные (содержания Au, Cu и Fe, а также степень окисления) и б) Категориальная входная переменные (литологическая характеристика).
Алгоритм Кохонена путем размещения центров радиальных элементов пытается распознать однородные кластеры в множестве наблюдений. В результате итеративной процедуры обучения алгоритм стремится к тому, чтобы элементы сети, соответствующие кластерам, расположенным близко друг к другу в пространстве входов, находились бы близко друг от друга и в топологически упорядоченной выходной сети (рис. 2). После того, как кластеры выявлены, все анализируемые пробы помечаются метками принадлежности к разным кластерам, соответствующим предварительно выделенным четырем природным типам руд. Таким образом кластерный анализ позволил выполнить кластеризации интервалов опробования по принципу максимальной однородности и подтвердил рабочую гипотезу выделения четырех природных типов по двум критериям:
-по 1-му критерию: слабоизмененные и интенсивноизмененные руды;
-по 2-му критерию: первичные и окисленные руды.
Рис. 2. Топологическая карта диалогового окна Нейронные сети. Результаты
По результатам прямых данных лабораторных исследований, а также выявленных рабочих зависимостей между обогатимостью руды и вещественными критериями имеется возможность выполнить построение ГТБМ, в которой для каждого блока геологической модели будут рассчитаны показатели извлечения золота с учетом неоднородности распределения в трехмерном пространстве. При необходимости в состав ГТБМ могут быть включены показатели физико-механических свойств руд. Существует несколько подходов к построению ГТБМ месторождений, основанных как на прямом, так и косвенном способе задания технологических параметров в блочную модель.
Прямой способ базируется на интерполяции параметров в блочную модель, используя для интерполяции только экспериментальные данные лабораторных исследований руд в ходе проведения геолого-технологического картирования.
Косвенный способ основан на прогнозировании показателей переработки руды по их зависимостям от геолого-минералогических факторов (способ наиболее перспективный для использования при ГТК). При отборе проб ГТК рекомендуется все композитные пробы разделять на две выборки: обучающую и прогнозную. На первом этапе лабораторных исследований выполняется необходимый комплекс технологического тестирования и изучения вещественного состава технологических проб обучающей выборки. В результате построения и анализа нейросетевых моделей и регрессионных уравнений на обучающей выборке, состоящей из 30–35 композитных технологических проб, устанавливается зависимость между коэффициентом извлечения золота и, например, содержанием Au, степенью окисления или содержанием сульфидов, что позволяет в последующем рассчитать технологические показатели обогатимости уже для прогнозной выборки (30–40 композитных проб) через зависимости технологических показателей от вещественных критериев обогатимости, полученных более дешевыми методами исследований, которые будут определены в ходе анализа обучающей выборки. Перспективность аналитического подхода основана на высокой экспрессности и надежности оценок технологических характеристик руд в сочетании с более низкой трудоемкостью выполнения работ.
В последующем рекомендуется использовать результаты разработанной ГТБМ для дальнейшего геолого-технологического изучения месторождения и в процессе производственной деятельности предприятия для проектирования и планирования производства.
Выводы
1. Состав руды неоднороден. Так же, как и содержание полезного компонента, коэффициент извлечения — величина переменная, определенным образом распределенная в трехмерном пространстве
2. Одной из основных проблем технологического моделирования является постоянная ограниченность количества лабораторных определений коэффициента извлечения.
3. Выявленные зависимости между обогатимостью руды и вещественными критериями, а также разработанные методики моделирования позволяют выполнить построение геолого-технологической блочной модели, в которой для каждого блока трехмерной модели будут рассчитаны показатели извлечения золота.
4. Классификация руда/порода при использовании ГТБМ основывается на применении к блоку модели условия рентабельности в виде безубыточного бортового содержания, составной частью которого является извлечение полезного компонента и затраты на переработку.
Комментарии, отзывы, предложения
Мальцев, 12.12.20 09:30:02 — Николаю
1. Вы правы для нейросетевой модели для поиска природных кластеров использовались в приведенном примере -1070 проб (сотни проб заметно кстати и по рисунку)
2. В мое понимание «безубыточного» разумеется, так же входят и технологические обстоятельства (упорная руда с пониженным извлечением)
СНС, 13.12.20 13:57:51 — Автору
Как влияет на результаты моделирования фундаментальная погрешность опробования и самородковый эффект.
Мальцев, 13.12.20 17:35:15 — СНС
И то и другое влияет в зависимости от того насколько то и другое велико. Однако невозможно объять необъятное, как и подробно ответить на ваш вопрос в рамках обсуждаемой темы статьи. Извините
СНС, 14.12.20 06:15:32 — Мальцев, 13.12.20
Подробно не надо. Попадет в результатах опробования пиковая проба, что с ней будет?
Мальцев, 14.12.20 08:25:36 — СНС
ураганные пробы при любом моделировании (нейросетевом или блочном) ограничиваются путем замены
СНС, 14.12.20 09:51:52 — Мальцев, 14.12.20
В подсчетах для ГКЗ до сих пор, по-моему, используют метод Когана. Но в зарубежных программах должно быть что-то другое? Или вы тоже по методу Когана режете?
Мальцев, 14.12.20 17:53:40 — СНС
Методы ограничения разные используем
Для ГКЗ - чаще Коган и график персентиль-функции
СНС, 15.12.20 15:48:06 — Мальцев, 14.12.20
Спасибо. При высокой изменчивости содержаний безубыточность вероятно определяется с некоторой погрешностью. Погрешность исходных данных и итоговой модели как-нибудь оценивается?
Генералов В.И., 15.12.20 17:39:44 — Е.Н.Мальцеву
Ох, уж эти англицизмы (персентиль-функции) и наукообразие! Пишите и мыслите по русски! Метод Когана, это директивный, единообразный (статичный) и простой (народный) по расчету метод. Существует множество других математических динамических методов усреднения (сглаживания) наблюдений. И они имеют практический (геологический) смысл, если объяснить пределы допущения к их применению. Ошибка всех "моделистов" заключается в том, что они показывают красочные компьютерные картины, не объясняя алгоритмов их составления.
Предлагаю поразмыслить об устоявшемся принципе исключения "ураганных содержаний" как таковых. Все единогласно относят их к "вредным содержаниям", которые необходимо непременно исключить из стройных рядов рядовых концентраций. Но с геологической точки зрения "ураганные пробы" объективны и имеют определенную геологическую позицию (читайте комментарии к месторождению Федоровское-1). По моему мнению, их надо тоже учитывать в подсчете запасов, в пределах участков их развития (гнезда). Учет гнезд в подсчете отражает неравномерность распределения металла в контуре рудного тела. Такую природную неравномерность можно обнаружить только в ненарушенном керне. Приведу учебный пример. В пробе из скважины RC-бурения (шлам) диаметром 124 мм было получено 5г/т золота на длину пробы в 1,0 м. Если бы бурение было колонковое с диаметром керна 63,5 мм и выходом более 95%, то была бы обнаружена одиночная кварц-сульфидная жилка мощностью 3 мм. Содержание золота в ней 500 г/т. Одиночная жилка имеет вертикальное падение. Содержание в 500г/т никакого отношения к рудному телу не имеет, так как не соблюдены кондиции по мощности и способу разработки.
Магадан, 16.12.20 03:32:56 — Генералову
Зато как научно: «график персентиль-функции»... Я сразу проникся и понял, что без науки не обойтись
Mad Miner, 16.12.20 10:06:40
Выбрасывание ураганных проб приводит в итоге (не всегда, конечно) к повышенным потерям с хвостами.
Старый, 16.12.20 10:14:07 — Генералов В.И., 15.12.20
Жив что-ли метод Когана? Я думал, что столько не живут. Ничего другого с тех пор не придумали? Я думал, что все пиковые пробы в современных программах моделирования урезаются автоматически и думать о них уже не надо. Хотя, какой то алгоритм урезания все равно должен быть, хоть в программах, хоть вручную.
С удовольствием почитал бы о современном положении. Написал бы кто-нибудь. Может вы или Мальцев или вместе?
СНС, 18.12.20 16:25:25 — Мальцеву
Если вы не возражаете, я вам сброшу перевод статьи на редактирование
"Контроль содержания ценного компонента при подземной добыче с учетом теоретических основ опробования: примеры с золотодобывающих предприятий". В ней про погрешности опробования интересно написано. А у нас этот вопрос, по-моему никак не двигается
Мальцев, 24.12.20 09:05:10 — СНС
Конечно, любая полезная информация приветствуется, хотя погрешности опробования - это отдельная (от ГТБМ) тема
Евгений Мальцев - evg65@mail.ru
Николай, 11.12.20 21:30:27 — Автору
Вы пишете про обучающую и прогнозную выборки композитных проб, причём каждая из них исчисляется первыми десятками (согласно тексту статьи). Если для построения регрессионных моделей этого может хватить, то как быть с нейросетевой моделью, размер обучающей выборки для которой исчисляется сотнями случаев?
Термин "ББТ" видится несколько некорректным, раз уж добываемые и перерабатываемые объёмы руды определяются напрямую по рентабельности этого процесса. Ведь содержание в блоке/домене/выемочной единице может оказаться и выше безубыточного, но технологические обстоятельства (упорная руда с пониженным извлечением) заставят этот блок опустить в "забаланс".