CC BY
26
УДК 519.24: 621.396.01 П.В. Васильев
ПОСТРОЕНИЕ ДИНАМИЧЕСКИХ ГОРНОГЕОЛОГИЧЕСКИХ МОДЕЛЕЙ МЕСТОРОЖДЕНИЙ ДЛЯ ОПТИМИЗАЦИИ УПРАВЛЕНИЯ ЗАПАСАМИ МИНЕРАЛЬНОГО СЫРЬЯ
Семинар № 1
В задаче построения комплексных геологических моделей месторождений широко используются многочисленные признаки вещественного состава, структурнотекстурные характеристики и физикомеханические свойства горных пород и руд, формирующие базу данных опробования вдоль скважин и разведочных пересечений. Все они в той или иной степени влияют на технологические показатели добычи и переработки руды, определяя качество исходного минерального сырья [2].
Среди множества качественных характеристик полезного ископаемого выделяют [3]: показатели назначения (состав, содержание вредных или полезных примесей и др. потребительские свойства); технологичности (дробимость, измельчаемость, обога-тимость, крепость, абразивность, фло-тируемость, растворимость и др.); сохранности (слеживаемость, окисляе-мость, самовозгораемость и др.). При этом важен не только уровень характеристик, но и стабильность качества во времени. При колебаниях качества нарушается оптимальный режим переработки сырья, ухудшаются техникоэкономические показатели, снижается извлечение полезного компонента в готовую продукцию, повышаются потери.
В связи с этим, определение наиболее информативных интегральных
показателей квалиметрической оценки руд в массиве имеет важное практическое значение. Особенно это является актуальным для решения задачи построения экологоэкономической блочной модели месторождения с целью оптимизации этапов отработки запасов и поиска наилучших контуров ведения открытых горных работ.
1. Прогнозирование качества руды в массиве
С позиций геометрического анализа рассматриваются два класса структур - монолитные и дезинтегрированные. В породном массиве монолитная руда имеет трещиноватость природного или искусственного происхождения. В свою очередь дробленая или измельченная масса при складировании образует навалы с достаточно плотной упаковкой, происходит слеживание и коагуляции тонких фракций. И в том и другом случае система хорошо описывается статистикой множеств вкраплений и частиц. Исходя из этого вопрос оценки параметров монолитных и дезинтегрированных структур рассматривается на основе методов стереологии как раздела стохастической геометрии, позволяющей выполнять реконструкцию геометрических структур 3Э тел по измерениям, выполненным в пространствах более низких размерно-
стей (вдоль линий в Ю и по плоскостям в 2Э).
В общем случае пробы или куски руды характеризуются некоторым набором трех типов признаков [ ]: геометрических ц (таких как размер минеральных зерен, объем, площадь поверхности, фактор формы, случайная хорда частиц, площадь случайного сечения и т.д.), вещественных в (содержание химических элементов или сочетаний компонентов, состав минералов) и физических % (плотность, магнитная восприимчивость, электропроводность, микротвердость и т.д.). Таким образом, совокупность проб или частиц может быть описана дифференциальной функцией состояния
у(г), в,%) - плотностью распределения объема выборки проб (реальных или имитационных фрагментов) по элементарным фракциям [ц,ц + йц] •
[,в + йв] [,% + й%] Однако в данной многомерной форме функцию практически трудно использовать на практике, как это делается в задачах классификации горных пород и расчленении разрезов [7]. В горнопромышленной геологии, в частности при добыче однокомпонентного минерального сырья, необходимо осуществлять снижение размерности пространства признаков на основе поиска наиболее информативных геологотехнологических показателей. Если для произвольной геологической породы при классификации стоимость является условной величиной, то для выделения промышленных типов руд именно стоимость должна служить базовым критерием распознавания принадлежности пробы к кондиционным участкам. Величина бортового содержания, особенно для тонковкрап-ленных бедных руд, не учитывает структурно-текстурных различий и
физико-механических свойств технологических сортов руд. Определим интегральную стоимость пробы разведочного пересечения, которая далее должна использоваться в процедурах построения блочной-каркасной модели месторождения.
В суммарных затратах по переработке руды основную долю составляют расходы энергии на дезинтеграцию и разрушение [6]. Согласно определению [1] под энергопоглощением горной породы понимается расход полезной энергии для превращения массива в горную массу требуемого состава по крупности, перемещение и укладку в отвал для пустых пород и переработку полезного ископаемого. При этом подчеркивается, что расходуется полезная, а не фактическая затрачиваемая энергия, выражая собой закон перехода количества энергии в качество горной массы и её перемещение. Таким образом, величина энергопоглощения определяется как зависящая от свойств горных пород и технологических параметров производственного процесса и определяется выражением [1]:
Е ■ I
Е = (1)
где Е - удельное технологическое энергопоглощение, дж/кг; Е - сопротивление, оказываемое горной породой, кг; 1 - длина пути преодоления сопротивления и перемещения горной породы, м; Р - количество продукции, произведенной в рабочий цикл, кг.
В предлагаемом подходе подготовка горно-геологической модели месторождения для оптимизации управлениями запасов состоит в следующем. Для разведочных проб месторождения задаются V, с, т — функции плотности, определяемые в каждой точке трехмерного простран-
ства; при этом ^х, у, ъ, 1) — себестоимость руды; с(х, у, ъ, 1) - стоимость разработки руды и пород. Текущее время определяет существующую на данный момент рыночную стоимость полезного компонента. Тогда имеем по аналогии со стационарным случаем [8] следующее выражение для дисконтированной прибыли от разработки данного элементарного объема: т(х, у, ъ, 1) = ф, у, ъ, 1) - с(х, у, ъ, 1) (2) При этом себестоимость руды напрямую зависит от доли полезного компонента в ней с учетом полного коэффициента извлечения, горного и обогатительного производства. Стоимость разработки определяется, главным образом энергозатратами
Переходя от энергозатрат к стоимости получим для удельной стоимости пробы единичного интервала следующее выражение:
Определим интегральную функцию извлечения кусков руды или частиц в зависимости от их физических свойств в технологическом оборудовании как е(%). Тогда в каждом диапазоне масштабов разработки месторождения можно определить ожидаемые показатели рудоподготовки (выход у , качество вк , извлечение ёк):
п
Г к
‘=1 (3)
п
вк = Тк
1=1
ёк =ГЛ а,
где а - среднее содержание полезного компонента в пробе.
Представленные формулы позволяют определить основные показатели всего технологического процесса рудоподготовки - от вскрытия месторождения до получения готового концентрата. Таким образом, важная задача прогнозирования качества за-
пасов состоит в том, чтобы определить потенциальную функцию у(ц, Р,^) «объёмно-качественно-физического» состояния руды по данным о вещественном составе, структурно-текстурным характеристикам и физико-химическим свойствам массива месторождения. Подобного класса задачи могут быть эффективно решены лишь методами стохастического моделирования, поскольку в технологический процесс добычи и переработки вовлечены значительные объёмы исходного минерального сырья, неоднородные по своей природе.
Кинетика разрушения кусков руды моделируется с помощью комбинирования механизмов раскалывания, откалывания и абразивного истирания. При этом, однако, кусок рассматривается как гомогенная масса, характеризуемая усредненными физикохимическими показателями. Для изучения эффекта раскрытия полимине-ральной кристаллической горной породы под действием разрушающих
нагрузок, необходимо принять во внимание некоторые дополнительные влияющие факторы. При разрушении полиминеральной системы сеть трещин может быть сформирована различными путями:
- случайным образом, когда точки центров разрушения выбираются в пространстве в соответствии с пуассоновским процессом заданной интенсивности (подобные условия возникают преимущественно в промышленных мельницах и дробилках с хаотичными импульсными нагрузками на частицы в рабочей зоне);
- по межфазным границам раздела, когда природная прочность
срастания минеральных фаз достаточно мала или искусственно уменьшенная за счет специальных физических воздействий;
КРУПНОСТЬ ЧАСТИЦ, мкм
Рис. 1. Раскрытие рудной фазы и породы блока магнетитового кварцита
- вдоль зон и линий сосредоточения искусственно созданных дефектов кристаллической структуры с использованием особых методов дезинтеграции;
- по интеркристаллитным границам, то есть по контактам зерен минералов;
- по плоскостям спайности кристаллов;
- по регулярной решетке или мозаике разбиения.
В технологических процессах взрывания, дробления, самоизмельче-ния или шарового измельчения руд проявляются преимущественно случайное и межфазное разрушение. Марковская модель процесса сокращения крупности и раскрытия минеральных фаз рассматривается в работе [10]. Данная схема расчета раскрытия минеральных фаз и прогнозиро-
вания извлечения полезного компонента реализована в программном обеспечении МтНЬег, являющегося подпрограммой интегрированной системы веоЫоск. Функция плотности распределения раскрытия для двухфазной системы магнетитового кварцита представлена на рис. 1.
В результате прогнозирования извлечения полезного компонента для каждой пробы корректируется величина бортового содержания, необходимого для получения готовой продукции заданного качества.
На базе предложенных алгоритмов стохастического моделирования и разрушения микроструктур разработан модуль МтНЬег, включенный в состав интегрированной программной системы. Перспективными объектами исследований являются руды твёрдых полезных ископаемых с тонковрап-
ленными ценными минералами для обоснования их технологической классификации по сортам.
Подготовка блочной модели для оптимизация контуров карьера
Наиболее широко применяемый в настоящее время в составе ряда зарубежных горно-геологических интегрированных пакетов метод оптимизации контуров открытых горных работ основан на алгоритме теории графов, который был предложен в 1964 году Ёерчсом и Гроссманом. В настоящее время установлено, что данный алгоритм теории графов также имеет ряд существенных недостатков для месторождений со сложной геологией и длительным сроком жизни карьеров. Это связано с тем, что реальное развитие карьерного поля ограничивается стохастическими горно-технологическими процессами и геологическими свойствами этих выработок, что приводит в целом к неоптимальным результатам. Изменения в базе данных также требуют полного перезапуска известных алгоритмов (вариантов, подвижного конуса, метода из теории графов), приводящих к длительной компьютерной обработке и не позволяющих выполнять достоверную оценку результатов недропользования.
Пусть функция а(х, у, 2) определяет угол в каждой точке пространства и Б является семейством поверхностей, наклон которых ни в одной точке по отношению к постоянной горизонтальной плоскости не превосходит а. Если обозначить через V семейство объемов, соответствующих семейству поверхности Б, то требуется найти среди всех возможных формирований объёмов тот, который максимизирует дисконтированную (приведенную на весь срок эксплуатации с учетом ог-
раничений) ценность включаемых в объём выемки рудных тел. Требуется найти такие границы открытых горных работ, в которых при разработке месторождения будет получена максимальная прибыль. При этом должна быть достигнута минимизация количества отходов в условиях экономических, геологических и технологических ограничений, обеспечен минимальный экологический ущерб окружающей среде.
Проблема оптимального планирования добычи в карьерах для блочных моделей месторождений определяется как отыскание последовательности выемки блоков в таком порядке, чтобы максимизировать чистый дисконтированный доход при следующих основных ограничениях: производительность фабрики обогащения, объём извлечения в заданный период времени, ограничения на смешивание типов руд, ограничения, связанные со складированием руды и отвалами пород, ограничения по логистике, экологические допуски и налоговые сборы.
Вводятся следующие параметры и ограничения: Т - число этапов планирования добычи; N - число блоков рудного тела; с' - прибыль (в NPV) от добычи блока 1 в период 1; О - подмножество рудных блоков; Ш подмножество блоков пустой породы; _ масса блока 1; 1т - масса обогащаемой руды в период времени 1; Б1 - множество блоков, которые необходимо удалить до отработки блока 1.
если блок 1 отработан в
1,
период времени между 1 и 1
0, иначе.
10 - нижняя граница массы руды, обогащенной в период времени 1; ы‘0 -верхняя граница массы руды, посту-
пающей на обогащение в период времени t; u‘w - верхняя граница массы породы, добытой в период времени t.
Тогда формулировка проблемы оптимизации добычи в карьере записывается следующим образом [12]. Необходимо максимизировать:
T N N
z = X X -С)xt-1+X cTX
t=2 i=1 i=1
при следующих ограничениях:
XtX = mi
ieO
Xt(xt - xt-1) = mt, t = 2,3,...,T.
ieO
XtX < uW
ieW
Xt.(x‘ - x'-1) < u‘, t = 2,3,...,T.
l ^ l l ' W* 111
ieW
xt-1 < xt, t = 2,3,...,T. x' < x‘j , t = 1,2,...,T, j eSi; i = 1,2,...,N;
10 < m < u0, t = 1,2,...,T;
x‘ = 0,1, для всех i, t.
Указанные ограничения гарантируют, что должна обеспечиваться заданная производительность рудника и фабрик обогащения, а также то, что масса удаляемой в отвал породы не должна превышать предписанные верхние границы. Параметры ограничения также указывают на то, что каждый блок удаляется в заданный период времени (этап отработки) при заданных углах откосов бортов карьера.
Сравнение результатов оптимизации контуров карьеров с использованием пакета Whittle FourX (Gemcom) с учетом прогнозирования извлечения в каждом блоке и в предположении о его постоянстве для всего месторож-
дения в целом приведено в работах [9, 11]. Установлено, что для месторождений с несколькими типами руд и технологическими сортами, использование извлечения как константы существенно искажает как предельный контур карьера так и оптимальные последовательности этапов отработки.
Оптимизация системы разработки облицовочного камня
При промышленной добыче облицовочного камня пильные "блоки" нарезаются оборудованием непосредственно в карьере или каменоломне. Некачественный блок бракуется и отправляется в отвал.
На месторождении облицовочного известняка Канфанар в Хорватии блоки имеют следующие размеры: длина не менее 2,0 м, ширина не менее 1,2 м и высота не менее 1,5 м. Блок должен быть однородной окраски без каких-либо трещин, микротрещинноватости, пустот, пятен и т.д. По этой причине когда оценивается качество, то "хорошим материалом" считается 100 % годный блок. Предельный фактор утилизации колеблется для разных каменоломен и зависит от тектоники на месторождении, технологии добычи, однородности текстуры, цвета и т.д. В настоящий момент, в больших карьерах по добыче камня факторы утилизации варьируют от 2,5 % до 70 % по объёму не принимая в расчёт перекрывающие породы.
Исходные условия отработки запасов камня.
Рабочие плоскости при эксплуатации месторождения облицовочного камня выбираются в соответствии со следующими критериями:
1. Согласно стратификации месторождения камня (если она имеет-
ся; мраморные или гранитные месторождения обычно не имеют напластований).
2. Согласно плоскостям разрывных нарушений путём выбора нарушений, которые перпендикулярны друг другу во всех трёх плоскостях.
3. Принимая во внимание, что оборудование может нормально ра-бо-тать только в том случае, если рабочая поверхность выработки отклоняется от горизонтальной плоскости не более чем на 15 градусов.
4. Объединив все эти требования можно получить наилучший выход качественной продукции и высокий фактор утилизации.
Подсчет запасов облицовочного камня на месторождении известняков Канфанар ведётся по методу параллельных геологических разрезов, в предположении, что зону подсчета слагают:
- Кондиционный известняк - на 100 % годится для экскавации и из полностью подсчитанного объёма данного типа материала получают примерно 35 %-й выход коммерчески годных блоков.
- Средний известняк - этот материал не очень хорошего качества, но он всё ещё пригоден для добычи.
- Голубой известняк - материал окрашен в голубой цвет, не пригодный для рынка Необходимо лишь знать как много такого материала в месторождении.
- Некондиционный известняк -материал на 100 % не годится для блоков
Разработанный в [10] алгоритм оптимизации системы разработки для обеспечения максимального выхода высококачественной продукции включает следующие операции:
Шаг 1. Для рассеянных точек топографии выполняется триангуляция
Делоне и создается набор данных триангуляции, для точек контактов литологических пластов и полигонов складчатости создаются сеточные триангуляционные модели.
Шаг 2. На основе полученных наборов триангуляций для литологических поверхностей создаются оболочки рудных тел и внутри них - в соответствии с рудными сортами:
а) выполняется экструзия (вдавливание) триангуляции кровли каждого пласта до дна полной модели;
б) выполняются операции с множествами, чтобы получить оболочки тел пластов с конкретными типами пород.
Шаг 3. Подсчитываются исходные геологические запасы по каркасной модели
Шаг 4. В пределах рамок модели генерируется регулярная решетка с заданными размерами блоков и приращением азимутальной ориентацией модели по 1 градусу.
Шаг 5.Решетка блоков конвертируется в модель оболочек блоков, то есть параллелепипедов стандартных панелей камней.
Шаг 6. Выполняется операция пересечения множеств оболочек блоков и каркасных оболочек типов пород. Результирующая модель пересечений включает лишь панели внутри каркасных оболочек пластов. На этом шаге выполняется геометрическая процедура расчета пересечения ребер и вершин прямоугольных плит известняка с каркасными оболочками рудных тел или пластов.
Шаг 7. Выполняется операция разности между полученной моделью кондиционных плит известняка и модели складчатости, представленной набором триангуляций, чтобы найти все ненарушенные трещиноватостью плиты. Итоговая модель пересечений
включает лишь блоки вне разломов и структурных нарушений. При пересечении любой грани блока известняка с любым треугольником плоскости нарушения данный блок отбрасывается и выполняется переход к следующему блоку.
Шаг 8. Выполняется операция пересечения точек центров блоков с оболочками типов или сортов руд, чтобы пометить каждый блок соответствующим типом породы или сортом руды.
Шаг 9. Вычисляется полный объём и запасы модели для различного качества облицовочных блоков известняка. Определяется полученный выход годных блоков от максимально возможного выхода запаса, полученного на шаге 4.
Шаг 10. Если заданного максимума выхода блоков руды не достигнуто, то осуществляется переход к шагу 5. Иначе конец вычислений и печать отчёта.
Данный алгоритм оптимизации системы разработки реализован в виде программного модуля интегрированной системы моделирования и оптимизации отработки запасов месторождений облицовочного известняка.
Заключение
Для минимизации ошибок построения цифровых моделей месторождений и повышения достоверности подсчета промышленных запасов предложено использовать следующую методику:
• прогнозирование извлечения полезного компонента рудной массы осуществляется на этапе компьютерной обработки проб детальной и эксплуатационной разведок и при оперативном опробовании скважин и выработок;
• классификация руд по технологическим сортам выполняется с учетом расчетной величины извлече-
ния и ожидаемой энергоемкости процесса получения продукта с заданными качеством и выходом концентрата;
• оконтуривание залежей и геологических тел проводится путём триангуляции точек контактов проб «руда - порода» с последующей тет-раэдрализацией внутренних областей каркасных оболочек кондиционной руды; дополнительно используются нормали к поверхности по границам раздела;
• все пробы приводятся к одинаковому размеру путём деления интервалов опробования вдоль скважин на сегменты равной длины в пространстве, поскольку в операциях интерполяции дистанционное взвешивание проб различной длины или их проекций на вертикальную плоскость совершенно недопустимо;
• каркасная модель месторождения преобразуется в блочную модель, обеспечивая требуемую дискретность пространства для планирования порционных операций горнодобычных процессов при заданной селективности выемки пород;
• интерполяция геостатистиче-скими методами, такими как кригинг и стохастическая имитация, позволяет найти оптимальные значения показателей в блоках внутренних областей оболочек рудных тел для межсква-жинного пространства и на значительном удалении от опробованных горных выработков;
• ускорение интерполяции по методу кригинга (универсального, индикаторного и зонального) и методу ближайших соседей достигается за счет быстрого поиска проб по сети тетраэдров, вершинами которых являются точки центров проб;
• после построения комплексной детальной модели месторождения реализуется операция удаления бло-
ков, не соответствующих кондициям по минимальной мощности руды и породы или по метропроцентам;
• рассчитывается ценность каждого кондиционного блока;
• выполняется оптимизация выемки рудных блоков на основе алгоритмов теории графов и эволюционных алгоритмов - перерасчитывается проектные контуры предельного карьера и вложенные контура промежуточных этапов отработки запасов при максимизации дисконтированной прибыли (чистого дисконтированного дохода);
• корректируется стратегия отработки запасов на весь срок жизни горнодобывающего предприятия;
• осуществляется переход от перспективного к текущему и оперативному планированию добычи.
Кроме того, имитационно-иерархический подход к квалиметрии запасов рудных месторождений включает следующую последовательность методики оценки объемно-качественных характеристик:
1. Построение компьютерной геолого-технологической модели месторождения. Определение потенциальной раскрываемости сортов руд и минеральных фаз. Анализ сложности строения залежи.
2. Построение компьютерной эксплуатационной модели месторождения. Определение экономически оптимальной конфигурации карьерного поля или шахты для выемки запасов в течение планируемого срока
1. Анистратов Ю.И. Технологические потоки на карьерах. Энергетическая теория открытых горных работ. - М.: Глобус, -2005, - 304 с.
2. Ломоносов Г. Г. Формирование ка-
чества руды при открытой добыче. - М.:
Недра. - 1975, - 224 с.
жизни предприятия. Анализ сложности системы разработки.
3. Оценка интегральной стоимости месторождения по совокупности показателей вещественного состава и параметров системы отработки запасов.
При оценке сложности строения залежи пространственная модель разбивается последовательно на квадратные фрагменты уменьшающейся крупности D, покрывая технологически важный диапазон размеров (1м - 104 мкм), а затем совокупность фрагментов классифицируется по качеству и объемной доле в них полезной минеральной фазы. В сериях изображений, асимптотически воспроизводящих стереологиче-ские характеристики структуры, определяется совместная функция плотности вероятности образования фрагментов различного состава и размеров.
В 3Э варианте методики структура месторождения (рудные тела, однородные области минерализации и поля минеральных фаз) моделируются на основе построения диаграммы полиэдров Вороного. Технологический процесс добычи и переработки руды моделируется путём стохастического покрытия пространства залежи минерального сырья масштабируемыми пробными блоками, соответствующими объёму селективности выемки оборудования.
------------- СПИСОК ЛИТЕРАТУРЫ
3. Ершов В. В. Геолого-маркшейдерс-кое обеспечение управления качеством руд. - М.: Недра.- 1986, - 260 с.
4. Ермолов В. А. Геолого-экологичес-кое обеспечение управления качеством руд в процессах рудоподготовки. Учебное пособие. - М.: МГИ. - 1997, - 84 с.
5. Новожилов М.Г., Ройзен Я.Ш., Эр-перт А.М. Качество рудного сырья черной металлургии. - М.: Недра, - 1977, - 420 с.
6. Ревнивцев В.И., Гапонов Г.В., За-рогатский Ё.П. и др. Селективное разрушение минералов. - М.: Недра. - 1988. -285 с.
7. Родионов Д. А. Статистические решения в геологии. - М.: Недра, 1981. -231 с.
8. Коробов С. Д. Разработка оптимизационных методов горно-геометрического анализа при освоении рудных месторождений открытым способом. //Автореферат на соискание ученой степени д.т.н. - М.: МГГУ. 30 с.
9. Васильев П.В., Уиттл Д. Блочное моделирование месторождений при оперативной оценке запасов и оптимизации добычи руд. //Материалы 5-го международного симпозиума "Освоение месторождений
минеральных ресурсов и подземное строительство в сложных гидрогеологических условиях". Белгород, Виогем. Том 2, 1999, -С. 351-362.
10. Васильев П.В. Некоторые стохастические модели и алгоритмы квалимет-рии запасов минерального сырья. //Научные ведомости. Белгородский Государственный Университет. Серия: Информатика, Прикладная математика и Управление. Изд. БелГУ. -2006. -т. 1. -№3. -С.45-53.
11. Vassiliev P. Optimizing Open Pit Limits without and with Ore Dressing Predictions. 3nd Bienial Conference "Stratigic Mine Planning", Perth, WA, Mar, 1999, pp. 197-207.
12. Caccetta L. and Hill S.P. An application of branch and cut to open pit mine scheduling. Journal of Global Optimization, No.27, 2003, pp.349-365. НШ
— Коротко об авторе
Вacильeв П.В. - кандидат технических наук, вед. научный сотрудник маркшейдерского отдела ФГУП ВИОГЕМ.
Доклад рекомендован к опубликованию семинаром № 1 симпозиума «Неделя горняка-2007». Рецензент д-р техн. наук, проф. A.M. Гaльпepин.
----------------------------------- ДИССЕРТАЦИИ
ТЕКУЩАЯ ИНФОРМАЦИЯ О ЗАЩИТАХ ДИССЕРТАЦИЙ ПО ГОРНОМУ ДЕЛУ И СМЕЖНЫМ ВОПРОСАМ
Автор Название работыі Специальность Ученая степень
ТУЛЬСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ
МИСНИКОВ Олег Степанович Разработка научных принципов утилизации промышленных отходов с комплексным использованием ресурсов торфяных месторождений 25.00.3б д.т.н.
