Ошибка
  • Delete failed: '881470919214d6e255e94fd9de18fc42.php_expire'
  • Delete failed: '881470919214d6e255e94fd9de18fc42.php'
  • Delete failed: 'abf04075ca2b776a5d25119289e77d0f.php_expire'
  • Delete failed: 'abf04075ca2b776a5d25119289e77d0f.php'

Создать PDF Рекомендовать Распечатать

Алгоритм генерации борта карьера по производственному плану горнодобывающего предприятия

  • Автор (авторы):
    Арсланов Дмитрий Мерзагитович
  • Дата публикации:
    24.12.12
  • ВУЗ ИЛИ ОРГАНИЗАЦИЯ:
    Южно-Уральский государственный университет

Алгоритм генерации борта карьера по производственному плану горнодобывающего предприятия

Algorithm for Pit Edge GenerationbyProduction Plan of Mining Enterprise

Арсланов Дмитрий Мерзагитович

аспирант кафедры «Автоматика и управление»

Южно-Уральский государственный университет

arslanstudio@hotmail.com

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

 Summary: In article algorithm forpit edge generation is described that uses voxel model of rock bodies and mining reserve blocks. Pit edge at a voxel layeris shifted by the horizontal distance. It’s achieved by discrete convolution of layerand special image processing matrix filter that allows for pit edge stability parameters for variable rock types. Algorithm is intended for pit design and optimal open pit mine planning.

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

 Keywords: mining enterprise; production; open pit mine; pit edge generation; deposit voxel model; image processing matrix filter; discrete convolution; pit edge stability parameters for variable rock types; pit design; optimal planning.

Введение

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

 Существенный вклад в изучение проблемы оптимального планирования производства открытых горных работ внесли Г. Лерчс и И. Ф. Гроссман; М. Т. Пана; Д. Виттл; С. Д. Коробов и многие другие. В трудах этих ученых предложены различные оптимизационные методы для установления обоснованных границ карьеров и последовательности разработки запасов месторождения[12].

 Однако в них не рассматриваются вопросы синтеза детальных открытых горных выработок с учетом параметров устойчивости борта карьера по различным типам пород, принятой технологии открытой разработки, схемвскрытия. В результате предложенные методы не могут быть в чистом виде применены для решения задачи оптимального планирования, так как строят только оценочные оболочки карьеров, которые имеют низкий уровень детализации. Как следствие, они дают низкую точность в расчетах и требуют большого объема ручной обработки для получения окончательного результата. В итогеданные методы не обеспечивают достаточное для оперативного решения задачи быстродействие[4].

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

 Воксельная модель месторождения полезных ископаемых

 В основу алгоритма положено моделирование объектов горнодобывающего производства в виде объемных изображений воксельной графики [8].

 Рассмотрим изображенный на рисунке 1 участок недр, где введем правую тройку локальных базисных векторов ED = (eD(i3D)), i3D = 1, 2, 3, с общим началом в точке x0.А также правую тройку ортонормированных глобальных базисных векторов EW = (eW(i3D)), eW(i3D, i3D) = 1 с общим началом в опорном геодезическом пункте O = (0, 0, 0) [10].

Рисунок 1 – Участок недр.

 Положение x0, размеры и ориентация векторов локального базиса ED в глобальном базисе EW задается при помощи мировой матрицы преобразований W = (w(iM)), iM = (iM(i2D)),iM(i2D) = 1, 2, 3, 4, i2D = 1, 2:

x0 = WO’,eW = WeD.

 Для выполнения матричных преобразований векторы должны быть представлены в однородных или гомогенных координатах, для обозначения которых в данной работе применяется штрих [10]. На примере базисного вектора eвектор его гомогенных координат e’ = (e’(i4D)), i4D = 1, 2, 3, 4 найдется как

e’(i3D) = e(i3D)e’(4),e’(4) = 1.

 На примере вектора x0 координаты вектора после матричных преобразований могут быть найдены из его гомогенных координат по отношению:

 Множество векторов координат всех интересующих точек рассматриваемого участка недр в локальном базисе ED

 образует единичный куб XD = (xD), xD = (xD(i3D)), 0≤ xD(i3D) ≤ 1.

 Проведем дискретизацию внутреннего локального пространства единичного куба XDрассматриваемого участка недр путем построения объемной прямоугольной регулярной решеткиG3D = (g3D), g3D = (g3D(i3D)), изображенной на рисунке 2, согласно заданным размерам SB = (SB(i3D)) ячейки вдоль соответствующих векторов локального базиса ED.

Рисунок 2 – Дискретизация пространства участка недр.

 В результате становится известно количество SD = (SD(i3D)) узлов g3D объемной решетки G3D вдоль соответствующих векторов локального базиса ED:

 Связь между координатами g3D(i3D) и номерами b(i3D) узла g3D вдоль соответствующих векторов локального базисаEDзадается следующими соотношениями:

g3D(i3D) = (b(i3D) – 1)SB(i3D),

b(i3D) = 1, 2, …, SD(i3D), 0≤ g3D(i3D) ≤ 1.

 Пусть на рассматриваемом участке недр залегает NRтипов горных пород и подлежит добыче NDтипов полезных ископаемых. В результате обработки данных геологоразведочных и маркшейдерских работ, проведенных на участке недр, в узлах g3D объемной решетки G3Dстановятся известны значения трехмерных двоичных дискретных табличных функций тел горных пород R = (r(iR)), r(iR) = (r(iR, b)), iR = 1, 2, …, NR,

 При этом узел относится только к одной породе: r(iR1, b) = 1, r(iR2, b) = 0, iR1iR2.

 А также становятся известны значения трехмерных вещественных дискретных табличных функций пространственного распределения содержания полезных ископаемых D = (d(iD)),d(iD) = (d(iD, b)), iD = 1, 2, …, ND, 0 ≤ d(iD, b) ≤ 1.

 Данные функции представляют собой пространственные стационарные дискретные скалярные поля, описывающие геологическое строение участка недр, а также залегание на нем полезных ископаемых. Они представляют собой объемные изображения или воксельную модель. На рисунке 1 изображены поверхности равного значения разных скалярных полей.

 Пусть горнодобывающее предприятие планирует разработать на участке недр блок полезных ископаемых, который описывается трехмерной двоичной дискретной табличной функцией плана производства открытых горных работ s = (s(b)),

 Будем описывать производство открытых горных работ трехмерной двоичной дискретной табличной функцией карьера p = (p(b)),

 При этом разрешающая способностьSD объемной решеткиG3D достаточна для того, чтобы отобразить все необходимые элементы его борта.

Передаточные функции воксельной модели

 Борт карьера, представляет собой систему уступов. Основными параметрами уступа, показанными на рисунке 3, являются высотная отметка горизонта разработки месторождения hP= (hP(iR)), высота hB = (hB(iR)) и угол aS = (aS(iR)) откоса, а также ширина бермы wB = (wB(iR)) [6, 7].

Рисунок 3 – Параметры устойчивости борта карьера.

 Значения данных параметров, обеспечивающие устойчивое положение борта карьера, заданы для каждого типа пород iR при помощи специальных табличных функций, которые называются передаточными функциями воксельной модели TF = (hP, hB, wB, aS). С учетом функций тел горных пород соответствующие параметры оказываются сопоставленными каждому узлу g3D объемной решетки G3D в зависимости от типа породыiR.

 Таким образом, проложение w(iR, b) участка борта карьера, приходящееся на высотуSB(3) ячейки объемной решетки G3D, изображенное на рисунке 3, будет зависеть от типа породы iR, в котором находится участок борта карьера, и номера b(3) горизонтального слоя объемной решетки G3D вдоль базисного вектора eD(3), на котором вычисляется проложение:

h(iR, b) = x0(3) + SB(3)b(3) – hP(iR).

Оператор генерации участка борта карьера

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

 На начальном шаге алгоритма проводится инициализация значенийфункции карьераpнулевыми значениями: p(b) = 0, b(i3D) = 1, 2, …, SD(i3D), i3D = 1, 2, 3. Следом проводится инициализация значений функции карьера p в узлах g3D первого горизонтального слоя объемной решетки G3D значениями функции производственного плана s в этих узлах: p(b) = s(b),b(3) = 1, b(i2D) = 1, 2, …, SD(i2D), i2D = 1, 2.

 На каждом следующем последовательном шаге алгоритм вычисляет новые значения функции карьера pв узлах g3D текущего горизонтального слоя объемной решеткиG3D по данным нижележащего слоя. Номер текущего слоя вдоль локального базисного вектора eD(3) принимает значения b(3) = 2, 3, …, SD(3). Вычисления идут параллельно для каждого узла g3D слояb(3).Так достигается быстродействие алгоритма. Номера узлов слоя вдоль соответствующих локальных базисных векторов eD(i2D) принимают значения b(i2D) = 1, 2, …, SD(i2D).

 Для каждого типа породы iRв ходе вычислений находится произведение значений функции тела горных пород r(iR) и функции карьера p в узлах нижележащего слоя b(3) – 1:

Int(iR, beW(3)) = r(iR, beW(3))p(beW(3)).

 Так определяется пересечение сечения функции тела горных пород r(iR) и сечения функции карьера pна нижележащем слое b(3) – 1. Для каждого найденного пересеченияInt(iR, b eW(3)) строится приращение на величину проложения w(iR, b eW(3)), постоянную для данного пересеченияInt(iR, b eW(3))в данной породе iR на нижележащем слое b(3) – 1. Для этого используется оператор генерации участка борта карьера, который представляет собой квадратный матричный фильтр f = (f(iF)), iF = (iF(i3D)), iF(i2D) = 1, 2, …, SF, iF(3) = 1 обработки горизонтального слоя объемного изображения [11]. Матрица оператора имеет нечетную размерностьSF, равную

SF = 2RF + 1,

 где RF – радиус оператора, равный

 Элементы матрицы оператора принимают значения

 Обработка пересечения I(iR, b eW(3)) оператором f заключается в вычислении дискретной свертки [11]данных функций в соответствующих узлах beW(3):

b(i2D) RF, (1)

b(i2D) ≥ SD(i2D) – RF(i2D). (2)

 

 где RF = (RF, RF, 1).

 Генерация борта карьера

 Сечение функции карьера p на текущем горизонтальном слое b(3) объемной решетки G3D найдется как бинаризация суммы найденных дискретных сверток для каждой породы iR и сечения функции производственного плана s на текущем слоеb(3):

 Бинаризация представляет собой оператор пороговой обработки вида Bin(σ) = (Bin(σ(b)) [11],

 Для обозначения суммы значений в узлах g3D объемной решетки G3D использован вектор σ= (σ(b)).

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

Применение алгоритма

 Выполнение описанных шагов алгоритма позволяет построить модель карьераp разработки соответствующего планаs производства открытых горных работ. На рисунке 4 показана поверхность карьера, сгенерированная данным алгоритмом для воксельной модели участка недр, изображенной на рисунке 1.Внешняя линия является пересечением с рельефом земной поверхности. Устойчивая конфигурация борта карьера точно прослеживает геологическую ситуацию на участке недр. Для построения использовалась модель с размерами 1024 x 1024 x 175 вокселей. Размеры вокселя составляли 2 x 2 x 2 метра. Среднее время генерации всего борта карьера составило 0,27 секунды на персональном компьютере со следующей конфигурацией: материнская плата ASRock X58 Extreme; процессор Core i7 920 (2,8 ГГц); память PC3-10666 Hynix (4 ГБ); видеокарта GTX 580 (GV-N580SO-15I); жесткий диск Corsair CSSD-F60GB2 (60 ГБ); операционная система Windows 7 Home Premium x64.

Рисунок 4 – Пример работы алгоритма генерации борта карьера.

 Объем горной породы, вскрываемой построенным карьером, найдется как

 Объем добываемых полученным карьером полезных ископаемых определится как

 Объем извлекаемой найденным карьером пустой породы или вскрыши может быть рассчитан как

W(R, D, p) = R(R, p) – D(D, p).

 Коэффициент вскрыши построенного карьера составит

 Рассмотрим применение описанной воксельной модели месторождения и алгоритма генерации борта карьера для постановки задачи оптимального планирования производства открытых горных работ.

 Пусть горнодобывающему предприятию для производства горных работ отведены согласно государственной лицензии определенные границы, заданные в узлах g3D объемной решетки G3D трехмерной двоичной дискретной табличной функцией B = (B(b)),

 Согласно пункту 5 раздела 3 «Рациональное использование и охрана недр» закона Российской Федерации «О недрах» от 21.02.1992 г. № 2395-1 [3] горнодобывающее предприятие должно обеспечивать наиболее полное извлечение из недр запасов основных и совместно с ними залегающих полезных ископаемых и попутных компонентов. Тогда постановка задачи оптимального планированияможет быть сформулирована так: горнодобывающее предприятие должно добывать максимальный объем полезных ископаемых на месторождении, при этом производство в целом не должно быть убыточным и должно находиться в пределах допустимых границ. Постановка задачипри этом записывается в виде

 Здесь F0 – это граничный коэффициент вскрыши, который определяется по экономическим показателям и показывает, в каком отношении должен находиться объем вскрываемой карьером пустой породы и объем добываемых карьером полезных ископаемых, чтобы производство не было убыточным.

 Метод решения поставленной задачи предложен в [1]. Метод обработки и визуализации воксельной модели на графических процессорах описан в [2].

Заключение

 Воксельная модель дополняет существующую классификацию моделей месторождений полезных ископаемых. Ее отличиями от регулярной блочной модели являются высокая разрешающая способность дискретного мульти-скалярного поля, применение матричных преобразований, передаточных функций и интерполяции промежуточных значений, а также методов параллельной обработки при помощи матричных фильтров, в первую очередь ориентированных на высокопроизводительные вычисления общего назначения на графических процессорах (GPGPU) [9]. Тем самым достигается высокая точность моделирования и быстродействие. Единая воксельная модель способна заменить ряд разнородных моделей, используемых в горном деле для решения различных задач.

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

Литература

1. Арсланов Д. М. Генетический алгоритм поиска экономически эффективных границ карьера: [Электронный ресурс] // Управление экономическими системами: электронный научный журнал. 2012. № (45) УЭкС, 9/2012. URL: http://uecs.ru/instrumentalnii-metody-ekonomiki/item/1546-2012-09-20-07-35-09. (Датаобращения: 17.12.2012).

2. Арсланов Д. М. Метод воксельной растеризации и обработки: [Электронный ресурс] // RSDN Magazine. 2011. № 3. URL: http://www.rsdn.ru/article/alg/03-12-voxel.xml. (Дата обращения: 09.07.2012).

3. Закон РФ «О недрах» от 21.02.1992 № 2395-1. Раздел III. Рациональное использование и охрана недр: [Электронный ресурс] // Официальный сайт компании«Консультант Плюс». Правовые ресурсы. Федеральные законы и кодексы Российской Федерации (РФ): полный текст документов. 2012.URL:http://www.consultant.ru/popular/nedr/66_3.html. (Дата обращения: 17.09.2012).

4. Залит А. В. Опыт использования GEMS. // II Конференция пользователей системы диспетчеризации ГТК «Карьер» и горно-геологической системы GEMS, 23–26 октября 2006.

5. Коробов Д. С. Разработка оптимизационных методов горно-геометрического анализа при освоении рудных месторождений открытым способом: автореферат дис. д-ра технических наук.– М.: Изд-во МГГУ, 1994. – 34 с.

6. Методические указания к определению углов наклона бортов, откосов уступов и отвалов строящихся и эксплуатируемых карьеров. – Ленинград: ВНИМИ, 1972. – 165 с.

7. Ржевский В. В. Открытые горные работы: учебник для вузов. В 2 ч. Ч. 2: Технология и комплексная механизация.– 4-е изд. – М.: Недра, 1985. – 549 с.

8. Engel, K. Real-Time Volume Graphics / Klaus Engel, Markus Hadwiger, Joe M. Kniss, ChristofRezk-Salama, Daniel Weiskopf. – Wellesley: A K Peters, Ltd., 2006. – 488 p.

9. Green, S. DirectCompute Programming Guide: [Электронныйресурс] / Simon Green // NVIDIA GPU Computing Documentation. 2010. URL: http://developer.download.nvidia.com/compute/DevZone/docs/

html/DirectCompute/doc/DirectCompute_Programming_Guide.pdf. (Дата обращения: 09.07.2012).

10.  Moller, T. Real-TimeRendering / TomasAkenine-Moller, EricHaines, NatyHoffman. – 3rded. – Wellesley: AKPeters, Ltd., 2008. – 1027 p.

11.  NI Vision Concepts Manual: [Электронныйресурс] // National Instruments. 2005. URL: http://www.ni.com/pdf/manuals/372916e.pdf.(Дата обращения: 09.07.2012).

12.  Shahriar, K. A Study On The Optimization Algorithms For Determining Open-Pit and Underground Mining Limits: [Электронныйресурс] / K.Shahriar, K.Oraee, E.Bakhtavar. // Professor KazemOraee publication website. Paper presented at the 7th International Scientific Conference, Sofia, Bulgaria. 2007. URL: http://www.oraee.net/shenase/Portals/0/limits.pdf. (Дата обращения: 17.12.2012).

  vakperechen

ОБНОВЛЕННЫЙ СПИСОК ВАК 2016 г.
ОТ 19.04.2016  >> ПРОСМОТРЕТЬ
tass
 
ПО ВОПРОСАМ ПУБЛИКАЦИИ СТАТЕЙ И СОТРУДНИЧЕСТВА ОБРАЩАЙТЕСЬ:
skype SKYPE: vak-uecs
e-mail
MAIL: info@uecs.ru
phone
+7 (928) 340 99 00
 

АРХИВ НОМЕРОВ

(01) УЭкС, 1/2005
(02) УЭкС, 2/2005
(03) УЭкС, 3/2005
(04) УЭкС, 4/2005
(05) УЭкС, 1/2006
(06) УЭкС, 2/2006
(07) УЭкС, 3/2006
(08) УЭкС, 4/2006
(09) УЭкС, 1/2007
(10) УЭкС, 2/2007
(11) УЭкС, 3/2007
(12) УЭкС, 4/2007
(13) УЭкС, 1/2008
(14) УЭкС, 2/2008
(15) УЭкС, 3/2008
(16) УЭкС, 4/2008
(17) УЭкС, 1/2009
(18) УЭкС, 2/2009
(19) УЭкС, 3/2009
(20) УЭкС, 4/2009
(21) УЭкС, 1/2010
(22) УЭкС, 2/2010
(23) УЭкС, 3/2010
(24) УЭкС, 4/2010
(25) УЭкС, 1/2011
(26) УЭкС, 2/2011
(27) УЭкС, 3/2011
(28) УЭкС, 4/2011
(29) УЭкС, 5/2011
(30) УЭкС, 6/2011
(31) УЭкС, 7/2011
(32) УЭкС, 8/2011
(33) УЭкС, 9/2011
(34) УЭкС, 10/2011
(35) УЭкС, 11/2011
(36) УЭкС, 12/2011
(37) УЭкС, 1/2012
(38) УЭкС, 2/2012
(39) УЭкС, 3/2012
(40) УЭкС, 4/2012
(41) УЭкС, 5/2012
(42) УЭкС, 6/2012
(43) УЭкС, 7/2012
(44) УЭкС, 8/2012
(45) УЭкС, 9/2012
(46) УЭкС, 10/2012
(47) УЭкС, 11/2012
(48) УЭкС, 12/2012
(49) УЭкС, 1/2013
(50) УЭкС, 2/2013
(51) УЭкС, 3/2013
(52) УЭкС, 4/2013
(53) УЭкС, 5/2013
(54) УЭкС, 6/2013
(55) УЭкС, 7/2013
(56) УЭкС, 8/2013
(57) УЭкС, 9/2013
(58) УЭкС, 10/2013
(59) УЭкС, 11/2013
(60) УЭкС, 12/2013
(61) УЭкС, 1/2014
(62) УЭкС, 2/2014
(63) УЭкС, 3/2014
(64) УЭкС, 4/2014
(65) УЭкС, 5/2014
(66) УЭкС, 6/2014
(67) УЭкС, 7/2014
(68) УЭкС, 8/2014
(69) УЭкС, 9/2014
(70) УЭкС, 10/2014
(71) УЭкС, 11/2014
(72) УЭкС, 12/2014
(73) УЭкС, 1/2015
(74) УЭкС, 2/2015
(75) УЭкС, 3/2015
(76) УЭкС, 4/2015
(77) УЭкС, 5/2015
(78) УЭкС, 6/2015
(79) УЭкС, 7/2015
(80) УЭкС, 8/2015
(81) УЭкС, 9/2015
(82) УЭкС, 10/2015
(83) УЭкС, 11/2015
(84) УЭкС, 11(2)/2015
(85) УЭкС,3/2016
(86) УЭкС, 4/2016
(87) УЭкС, 5/2016
(88) УЭкС, 6/2016
(89) УЭкС, 7/2016
(90) УЭкС, 8/2016
(91) УЭкС, 9/2016
(92) УЭкС, 10/2016
(93) УЭкС, 11/2016
(94) УЭкС, 12/2016
(95) УЭкС, 1/2017
(96) УЭкС, 2/2017
(97) УЭкС, 3/2017
(98) УЭкС, 4/2017
(99) УЭкС, 5/2017
(100) УЭкС, 6/2017
(101) УЭкС, 7/2017
(102) УЭкС, 8/2017
(103) УЭкС, 9/2017
(104) УЭкС, 10/2017
(105) УЭкС, 11/2017
(106) УЭкС, 12/2017
(107) УЭкС, 1/2018
(108) УЭкС, 2/2018
(109) УЭкС, 3/2018
(110) УЭкС, 4/2018
(111) УЭкС, 5/2018
(112) УЭкС, 6/2018
(113) УЭкС, 7/2018
(114) УЭкС, 8/2018
(115) УЭкС, 9/2018
(116) УЭкС, 10/2018
(117) УЭкС, 11/2018
(118) УЭкС, 12/2018
(119) УЭкС, 1/2019
(120) УЭкС, 2/2019
(121) УЭкС, 3/2019
(122) УЭкС, 4/2019

 Федеральная служба по надзору в сфере связи и массовых коммуникаций

№ регистрации СМИ ЭЛ №ФС77-35217 от 06.02.2009 г.       ISSN: 1999-4516