Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Суменков М.С. Математические методы планирования открытых горных работ

.pdf
Скачиваний:
10
Добавлен:
23.10.2023
Размер:
8.51 Mб
Скачать

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

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

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

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

Неувязку локальных задач можно показать на примере оп­ ределения граничного коэффициента вскрыши .по формуле

 

 

 

Кг=

 

с " ~ с ° ,

 

 

 

 

 

(П.1)

 

 

 

 

 

Св

 

 

 

 

 

 

 

где Сп и С 0

— себестоимость

добычи

полезного

ископаемого

 

 

подземным

и

открытым

(без

учета

погашения

 

 

вскрыши)

способами, руб/т;

 

 

 

 

 

 

С в

— себестоимость

вскрышных работ, ,руб/м3.

 

При определении граничного коэффициента

вскрыши по фор­

муле

(II.1)

проектные

технико-экономические

 

показатели

по

данному карьеру пока еще не определены, поэтому значения

С п

С 0 принимаются как средние

(нормативные) показатели по карь­

ерам,

имеющим в какой-то

мере аналогичные

горно-геологиче­

ские

условия залегания

месторождений.

Полученные

в

резуль­

тате

выполнения проекта

показатели

Сп ,

С 0 и Св ,

как

правило,

отличаются от первоначально принятых из-за специфики нового карьера. Если их подставить в формулу (II.1), то полученное уточненное значение граничного коэффициента вскрыши будет отличаться от первоначально вычисленного. Это отклонение до­ стигает иногда 50%- Следовательно, изменится и конечная глу­ бина карьера. Возникает информационная неувязка, например,

между задачами определения

глубины карьера (она

уточнена)

и вскрытия месторождения,

требующая уточнения

способа

вскрытия.

 

 

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

1)метод декомпозиции;

2)метод агрегации;

3)комбинированный метод.

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

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

декомпозиции был

впервые

применен В. В. Данцигом л Р. Вул-

фом в 1960

г.

к

задачам

линейного программирования т г

И. В. Розеном

в 1963 г. к задачам нелинейного программиро­

вания.

 

 

 

 

М е т о д

а г р е г а ц и и основан на изучении простых локаль­

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

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

Автоматизированная система проектирования открытых гор­ ных работ (АСЛО) должна обеспечивать создание автоматизи­ рованной системы управления горным предприятием, которая бы функционировала в оптимальном режиме в соответствии с за­ данным критерием оптимальности. Кроме того, АСПО должна предусматривать:

1.Широкое использование эвристического проектирования, позволяющего формализовать опыт 'проектировщиков.

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

3.Использование ЭВМ для решения различных проектных

задач.

4.Использование информационной системы нормативных технико-экономических показателей.

Большое количество и сложность задач, решаемых при про­ ектировании карьера, не позволяют создать единую модель все­ го горного предприятия. Поэтому АСПО в структурном отноше­ нии должна состоять из ряда подсистем, соответствующих ло­ кальным задачам оптимизации. Локальные подсистемы согласо­

вываются

как по входам

(исходная информация),

так

и по

выходам

(результаты вычисления).

 

 

 

Различная технико-экономическая основа локальных проект­

ных задач

вызывает необходимость и формализации

их различ­

ными математическими языками. Одни задачи

целесообразно

описывать

в терминах линейного программирования,

другие —

в терминах

динамического

программирования,

третьи — в

тер-

минах метода статистических испытаний и т. д. Следовательно, каждая подсистема ЛСПО функционирует самостоятельно в со­

ответствии

с алгоритмом

решения задачи и локальным

критери­

ем оптимальности, а отдельные подсистемы оказывают

взаимное

влияние через'ВХОДЫ

I I выходы.

 

Для дальнейших

выкладок определим граф локальности за­

дач, под

которым будем

понимать графическое расположение

последовательности выполнения проектных локальных задач, отображающее их взаимосвязь. Каждой локальной задаче по­ ставлена в соответствие вершина графа, а взаимосвязями задач является ребро графа. Следовательно, каждое ребро графа ло­

кальности показывает, что из найденных параметров

в

результа­

те решения t-й задачи хотя бы один используется

в

качестве

исходной информации для решения у'-й задачи.

 

 

 

 

 

Информационной неувязкой і-н локальной,

задачи

по пара­

метру k назовем абсолютное значение

разности

 

 

 

 

 

 

А(* =

К * а *

I .

 

 

 

(П-2)

где a n h — проектное значение

k-vo параметра на выходе

послед­

ней локальной задачи графа;

 

 

 

 

 

о,/£ —численное значение /г-го параметра

на входе

г'-й

ло­

кальной

задачи.

 

 

 

 

 

 

 

Информационная неувязка локальных задач может быть

уменьшена или сведена практически к нулю за

счет

следующих

мероприятий:

 

 

 

 

 

 

 

 

1. Формирование рационального графа локальности

задач,

при котором наибольшее количество входящих

параметров

оп­

ределялось бы на

выходе предшествующих локальных

 

задач.

2. Применение

процедуры

последовательных

приближений.

3. Изменение топологии графа локальности

па каждом

цик­

ле последовательных приближений.

 

 

 

 

 

 

4. Применение ЭВМ для обоснования нормативных парамет­ ров, используемых в качестве исходной информации при расче­ тах и не являющихся выходом какой-либо локальной задачи.

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

Вычислительная процедура последовательных приближений состоит в последовательном выполнении циклов решения задач графа локальности. Начиная со второго цикла вычислений в ав­ томатизированной системе проектирования карьеров начинает функционировать обратная связь, позволяющая использовать технико-экономические параметры проекта z-ro цикла вычисле­ ний для решения задач графа локальности на 2 + 1 - м цикле вы-

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

 

 

Aik<eik,

(П.З)

где ЄЇ/І допустимая

информационная неувязка k-ro параметра

на входе t-й локальной

задачи.

 

Численное значение

параметра

ЄІЛ может быть определено из

условия необходимой

его точности

вычислений.

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

стигается

применением комбинированного 'метода

построения

АСПО. Изменение топологии графа осуществляется

изменением

в каждой

локальной задаче количества охватываемых технико-

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

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

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

Ниже рассмотрим математические модели некоторых локаль­ ных задач автоматизированной системы проектирования карьст ров, для которых предлагаются свои методы решения, опробо­ ванные на ЭВМ.

§. 2. ОБОСНОВАНИЕ ПРОМЕЖУТОЧНЫХ И КОНЕЧНЫХ КОНТУРОВ К А Р Ь Е Р О В

Одной из основных и наиболее сложных технико-экономнче- ских задач открытой разработки месторождений является уста­ новление оптимальных промежуточных и конечных контуров карьеров [10]. Промежуточными контурами выделяются очереди (этапы) разработки месторождения, а конечными — экономиче­ ски целесообразная зона открытых горных работ. Сложность за­ дачи оптимизации контуров карьеров вызывается в основном следующими причинами:

1.Необходимостью количественного учета динамики эконо­ мических параметров.

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

3.Потребностью в постоянном уточнении элементов залега­ ния месторождения и качественной характеристики полезного ископаемого в процессе эксплуатационной разведки месторож­ дения.

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

1.Пригодные для разработки открытым и подземным спосо­

бами.

2.Разрабатываемые только открытым способом (асбест, гра­ нит, глина, тальк, бедные железные руды и т. д . ) .

Условию наиболее эффективной разработки месторождения отвечает принятие в качестве критерия оптимальности макси­ мума приведенной прибыли при обработке всего месторож­ дения.

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

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

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

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

1. Вычерчиваются погоризонтные геологические планы ме­ сторождения. Причем типы руд по различным признакам выде­ ляются на различных планах (например, ца геологических пла­ нах одного комплекта выделяются типы руд по содержанию же­ леза, другого—но содержанию ванадия и т. Д:).

24

2.Наносятся технически целесообразные варианты конечны\ контуров карьера, внутри каждого из них — варианты промежу­ точных контуров.

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

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

которая позволяет формировать

граф

многовариантности плана

и присваивать каждому ребру

графа

экономическую оценку —

величину приведенной прибыли

ГОКа.

 

Каждой вершине графа соответствует определенное положе­ ние горных работ в пространстве, а каждому ребру графа — оп­ ределенный варианто-год развития горных работ. Следователь­ но, каждый конечный контур карьера будет характеризоваться отдельным пучком графа многовариантности.

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

6. Проводится повторное моделирование горных работ на ЭВМ для нанесения промежуточных положений линии фронта разработки на логоризонтные планы. Необходимость .повторного моделирования вызвана тем, что память Э В М недостаточна для хранения всех промежуточных результатов вычислений.

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

Изложенная методика оптимизации конечных и промежуточ­ ных контуров карьера имеет следующие основные преимущества:

1.Возможность выполнения трудоемкой работы по геомет­ рическому анализу карьерного поля на ЭВМ, особенно для ме­ сторождений неправильной формы.

2.Учет взаимосвязи промежуточных и конечных контуров карьера.

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

4.Учет изменения эффективности вариантов вскрытия и раз-

вирня горных работ при изменении промежуточных контуров карьеров.

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

Упрощенный пример построения графа многовариантности промежуточных контуров карьера для пояснения сущности мето­ дики показан на рис. 1.

Выделенные варианты контуров карьера по поверхности при-

Рис. 1. Схема построения графа многовариантности промежуточ­ ных контуров карьера. Шифры:

О—исходное

положение

горных работ;

7 — конечные

контуры карьера;

/ , . . .

б — контуры

этапов (очередей)

разработки

месторождения.

ведены на рис. 1,о, а на рис. 1,6 показан граф многовариантно­

сти формирования этапов разработки

месторождения.

 

 

 

Каждый путь графа от начального события до конечного со­

ответствует определенному

варианту

выделения

этапов

разра­

ботки

месторождения. Так,

варианты

{0, 1, 3,

7)

и (0,

2,

3,

7)

различаются

только контурами первого этапа,

а варианты

(0,

2,

3, 7)

и (0, 2,

4, 7) —только

контурами

второго

этапа.

 

 

 

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

§ 3. ВЫБОР МЕСТА З А Л О Ж Е Н И Я ОБОГАТИТЕЛЬНОЙ ФАБРИКИ

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

Экономическая постановка задачи следующая. Пусть имеет­ ся п месторождений и задана очередность их швода в эксплуата­ цию. Кроме того, известна продолжительность строительства

26

каждого рудника, динамика изменения их производительноег.ч по горным условиям и производительность обогатительной фаб­ рики по сырой руде. Требуется выбрать точку расположения обо­ гатительной фабрики таким образом, чтобы максимизировать величину прибыли ГОКа за период отработки всех месторожде­ ний при условии транспортирования всей руды на данную фаб­ рику и таким образом учесть влияние выбранного местораспо­ ложения обогатительной, фабрики на технико-экономические по­

казатели

работы ГОКа.

 

 

Для построения

экономико-математической модели

сформу­

лированной задачи

поступаем следующим образом.

Критерий

оптимальности

(приведенная

прибыль за период отработки всех

месторождений)

выражается

следующей функцией

 

 

т

 

 

 

 

П -- Z J

 

(14-Е v

LjiZuЕ

V

л1

 

 

 

 

 

(1 + Е»"У

 

 

(II.4)

 

 

 

 

 

 

 

 

где П — общая прибыль ГОКа

за время

Т отработки

п место­

рождений,

приведенная

к началу

планируемого

перио­

да, руб.;

 

 

 

 

 

 

 

 

 

 

 

Д ' доход ГОКа от реализации

продукции

по

отпускным

ценам в t-u году, полученной из

руд 1-й залежи;

 

3'.— эксплуатационные

(за

вычетом

 

затрат

на реновацию)

и капитальные затраты

ГОКа

в

г"-м году, связанные с

добычей

и переработкой руды г'-й залежи, руб.;

 

б[ (т)—коэффициент

прироста

прибыли в т-м году по отноше­

нию к т—1-му году за счет действия технического

прог­

ресса;

 

 

 

 

 

 

 

 

 

 

 

"г(т) —коэффициент

снижения

себестоимости

1 т/км в т-м го­

ду по отношению к (т—1)-му году за счет действия тех­

нического

прогресса;

 

 

 

 

 

 

 

Ет— нормативный коэффициент

приведения

затрат;

 

/; — расстояние

транспортирования от

/-й залежи

до

обога­

тительной

фабрики,

км;

 

 

 

 

 

 

С.— себестоимость одного т/км в t-м

году перевозки руды от

г'-й залежи до фабрики (без учета отчислений на рено­

вацию),

руб.;

 

 

 

 

 

 

 

 

QJ— объем перевозок руды

в t-м

году от г'-й залежи до фаб­

рики, т;

 

 

 

 

 

 

 

 

 

 

 

/„+) — расстояние

от обогатительной

фабрики до

железнодо­

рожной

станции МПС,

км;

 

 

 

 

 

 

С'п,, — себестоимость

одного т/км в ^-м году перевозки

готовой

продукции с обогатительной фабрики до железнодо­

рожной станции МПС, руб.;

 

 

 

 

у'—коэффициент

выхода готовой продукции

из руды і-й за­

лежи

в if-м году.

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

6 x 0 0 =

?~3\

1.

6 » ( T ) = - f r - I .

(И.5)

 

 

д х - і _ 3 х - і

 

 

 

сх-\

 

 

Для упрощения

дальнейших

математических

выкладок вве­

дем следующие

обозначения:

 

 

 

 

 

 

 

 

 

 

 

т

 

 

t

 

 

 

 

 

 

 

 

[Д\-3\)

 

П [ І +в,(т)]

 

 

^-ZAZJ

 

 

^

 

 

.

(її.-»

 

 

 

 

 

г

 

 

 

 

 

 

 

 

%^

 

C{Q/

П

[1

+Ш]

 

 

 

 

SZJT\f

 

 

v

'

 

<IIJ>

 

 

 

 

л

Г

 

 

 

 

 

 

 

 

 

 

 

c i + i Q ' v )

п

а (т)і

 

 

 

 

 

 

— < r f b — •

<"••>

Тогда функция

(11.4)

примет вид

 

 

 

 

 

 

 

 

П = К - S V I - W I ' » + I -

 

(" - 9)

 

 

 

 

 

<=i

 

 

 

 

 

 

Представим

рудники точками

плоскости

х, у. Пусть t-му руд­

нику соответствует

точка

с координатами

а,-, 6,, і1 , " . Тогда

расстояния от каждой t'-й залежи

до обогатительной

фабрики

будут определяться как

 

 

 

 

 

 

 

 

 

 

і ^ у ^ - х у +

 

іь-уу

 

 

и расстояние

от обогатительной

фабрики до железнодорожной

станции МПС

 

 

 

 

 

 

 

 

 

 

 

 

 

1п+1

=

п+ і—х)г + {bn+ і — # .

 

 

Требуется

определить

такие

координаты

(х, у) обогатитель­

ной фабрики, чтобы целевая функция

 

 

 

 

п = ^ - 2 1

К Vfa-xF

 

+

fa-y)*

(п. ю)

 

 

 

 

i=i

 

 

 

 

 

 

 

достигала максимальной

величины.

 

 

 

 

 

28

В точке х, у, где достигается оптимальное значение функции (11.10), принимает минимальное значение функция

С = "І'я,- V{a-xf + (bi-yf, (П. 11)

і—і

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

Для нахождения

оптимального

значения

функции

(11.11)

предлагается следующий итерационный процесс. Строится

по­

следовательность точек Xh, Ук, сходящаяся к точке,

дающей

ре­

шение задачи, k—\, 2, . . .

 

 

точки х0 , у^

Последовательность начинается с произвольной

в качестве которой может быть выбрана точка

начала координат

(0, 0) . Если известна

точка {Xk, Ук),

то следующую

точку

нахо­

дят по правилу

 

 

 

 

 

 

 

 

 

 

 

(11.12)

(11.13)

 

V(at

-

xk? + (Ь£ -

ykf,

 

если

Vi^-x.f

+

(6, - у*? > е

' * Н

-

 

,

 

 

 

 

 

 

 

I Е ' - Е С Л И

Viai-x^+Xb-yJ*

 

 

<

г

 

 

•и є — достаточно малое положительное

число.

 

Итерационный процесс

получения

точек

последовательности

продолжается

 

до тех пор, пока

не будет выполнено следующее

неравенство:

 

 

 

 

 

 

 

 

 

 

 

У ( х к + 1 - х к у

+

( у к +

1 - у к у < г ,

(И.14)

т. е. расстояние между соседними точками

последовательности

не должно превосходить достаточно малой

величины.

В

случае,

если неравенство

(П.14)

выполнено,

полагаем

х = х к + 1 , у — у к + 1

и считаем, что оптимальное решение получено.

Соседние файлы в папке книги из ГПНТБ