
книги из ГПНТБ / Суменков М.С. Математические методы планирования открытых горных работ
.pdfскими параметрами локальных задач, что объясняется следую щими основными причинами:
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 |
+Ш] |
|
|
||
|
|
SZJ—T\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
и считаем, что оптимальное решение получено.