Многокритериальная оптимизация показателей взрывных работ на карьере Текст научной статьи по специальности «Математика»

Организационная структура компании ГУП «Мосгаз»

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

Для поиска подходов к решению этой задачи было выполнено исследование существующей структуры и проанализированы информационные потоки подразделений ГУП «Мосгаз», на основании чего построены про-цессно-событийные модели, отражающие принципы

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

Техническое задание проекта информатизации выбранного предприятия включает в себя следующие требования:

• проектирование инфраструктуры корпоративной компьютерной сети;

• использование клиент-серверных технологий и стандартных серверов баз данных;

• проектирование структуры и разработка программного обеспечения хранилища данных для накопления, хранения и обработки информации, полученной из разнородных источников;

• обеспечение прозрачности данных (возможность стандартного доступа к данным хранилища из других модулей корпоративной информационной системы);

• обеспечение разграничения уровней доступа для пользователей.

На основе проведенного анализа и сформулированных в техническом задании на разработку информационной системы целей и условий функционирования создаваемой ИИАС предложена интегрированная архитектура системы на основе хранилища данных.

В рамках предложенной архитектуры, информационная система:

• охватывает главные аспекты функционирования организации, являясь отражением основных процессов и результатов деятельности;

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

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

МНОГОКРИТЕРИАЛЬНАЯ ОПТИМИЗАЦИЯ ПОКАЗАТЕЛЕЙ ВЗРЫВНЫХ РАБОТ НА КАРЬЕРЕ

Д.К. Потресов, С.И. Сапожников

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

альной оптимизации (МКО) является одним из наиболее сложных. Это связано с тем, что оптимальное решение х задачи МКО, в общем слу-

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

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

1) степень дробления зависит от длины заряда в скважине - l3ap, расстояния между скважинами в ряду - a и между рядами скважин - b;

2) ширина развала горной массы после взрыва зависит от длины заряда в скважине - l3ap, расстояния между скважинами в ряду - a и между рядами скважин - b, от величины интервала замедления - т;

3) расстояние безопасного удаления для объектов зависит от длины заряда в скважине - l3ap, расстояния между скважинами в ряду - a.

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

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

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

Ye=|e* min, (1)

где £ - подвергаемая оценке степень дробления горной породы; £ - оптимальная степень дробления горной породы.

Критерий оптимальности для оценки ширины развала горной массы YB после взрыва имеет следующую математическую запись:

YB = B ^ min, (2)

где B - подвергаемая оценке ширина развала горной массы; - коэффициент, определяемый по формуле:

Г1, B > Bmin =1 ■ , (3)

[«», B < Bmin

где Bmin - минимально допустимая ширина развала горной массы, определяемая используемым оборудованием.

(8)

Критерий оптимальности для оценки расстояний безопасного удаления YR представим следующим видом:

Yr = R ^ min, (4)

где R - подвергаемое оценке безопасное расстояние; - коэффициент, определяемый по формуле:

_ = 11, R < Rmax

= К R > Rmax , (5)

где Rmax - максимально допустимое безопасное расстояние.

Рассмотрим некоторое пространство решений W, определенное на четырехмерном линейном пространстве действительных чисел R4:

W = {Xt = Озар ,a,b, т )t е R4 },

(6)

1зар ,a,b, т е R,t = 1, 2, ...

Пусть на W наложены некоторые ограничения в виде областей ßa, ßb, П|зар и ßT допустимых значений a, b, 1зар и т:

I е ß, ,aе ß ,bе а ,те ß . (7)

зар |зар a b т

Тогда с учетом (6) и (7) множество допустимых решений X может быть записано в виде:

X = { Xq е W : I е ß. ,a еа ,

q зар .зар a

bе ßb,те ат},q = 1, 2, ...

Совокупность критериев Y, характеризующих качество решения Xq, является трехмерным вектором, определенным в трехмерном линейном пространстве вещественных чисел R3, и связь между решениями и значениями критериев устанавливается отображением f, действующим из пространства решений W в критериальное пространство R3, то есть

f : X ^ Y с R3. (9)

Таким образом, множество достижимости Y в критериальном пространстве R3 может быть задано следующим образом:

Y =_{ (y., yB,yR)r е r3 }, yj = fj(Xq)е R, ^

j = 13 , r = 1, 2, ...

Предпочтения ЛПР в критериальном пространстве задаются бинарным отношением строгого порядка:

y' h>y" « yj < yj',y' * y',jе J . (11)

Решаемая математическая задача обладает следующими особенностями:

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

- анализу подвергаются критерии, выраженные математическими функциями и формирую-

щие невыпуклые множества достижимых целей;

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

Выбор метода МКО технико -экономических показ ателей взрывных работ на карьере

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

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

Методы без участия ЛПР могут применяться в системах, время выработки управляющих воздействий в которых является критичным, например, в системах операторского контроля и управления (8САБА), но непригодны для использования на более высоких уровнях, так как эти методы оперируют только одним критерием и полностью исключают участие ЛПР в процессе принятия решения.

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

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

тер решает очередную вспомогательную задачу оптимизации.

Многокритериальную информацию все чаще представляют с использованием современных компьютерных средств. Развитием итеративных человеко-машинных методов стали итеративные методы, основанные на визуализации паретовой границы, ориентированные на итеративный процесс принятия решений в форме диалога с пользователем, организуемого с помощью средств мультимедиа. Такие методы наиболее рационально применять для решения многокритериальных задач МЕ8-уровня, в которых математический аппарат должен сочетаться с удобным и информативным графическим человеко-машинным интерфейсом. Представим сравнительную характеристику итеративных методов, основанных на визуализации паретовой границы (табл. 1).

С точки зрения решаемой задачи наиболее важными являются удобство использования метода при трех критериях и простота его применения. Из таблицы 1 видно, что в наибольшей степени удовлетворяющим принятым критериям оценки методов является метод достижимых целей (МДЦ).

Тем не менее этот метод обладает двумя недостатками: наиболее выигрышным случаем применения классического МДЦ является решения задач МКО с фиксированным количеством стратегий, кроме того, метод крайне субъективен.

Таблица 1

Методы Критерий оценки метода

Удобство использования при трех критериях Простота применения метода Ограничение субъективного участия ЛПР

Метод проектирования на паретову границу - - +

Метод Корхонена-Лааксо - + +

Метод Pareto Step + - -

МДЦ + + -

Проблему субъективности МДЦ можно решить путем разработки сценария принятия решения, основанной на нахождении оптимального решения в ограничениях, выявленных в диалоге с ЛПР с помощью МДЦ.

Касательно другого недостатка, наиболее выигрышным случаем применения классического МДЦ является решение задач МКО с фиксированным количеством стратегий (вариантов решений), в которых возможно или представление ЛПР всего множества достижимых целей, или, по крайней мере, множество достижимых целей может быть получено перебором всех решений. В рассматриваемой задаче анализу подвергаются критерии, выраженные математическими нелинейными мультимодальными функциями, и количество

возможных вариантов решений стремится к бесконечности.

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

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

Таблица 2

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

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

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

Содержательное описание модернизированного МДЦ

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

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

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

Количество классов также определяется сред-неквадратическим отклонением по данному критерию, то есть при коэффициенте функции сред-неквадратического отклонения, равном 1/15, данный критерий будет разбит на 15 классов.

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

Для уменьшения субъективизма МДЦ применяется сценарий ограничения субъективного участия ЛПР в процессе принятия решения.

Рассмотрим сценарий ограничения субъективного участия ЛПР в процессе принятия решения. Его схема представлена на рисунке. После того как карты решений МДЦ построены, ЛПР может решать задачу тремя способами.

Первый из них основывается на стандартных механизмах МДЦ: ЛПР с помощью визуального анализа по карте решений интерактивно выбирает рациональную критериальную точку.

Второй способ заключается в сведении задачи к однокритериальной. Для этого от ЛПР требуется задать основной критерий, после чего задача решается методом основного критерия.

Наконец, ЛПР может использовать механизм поддержки принятия решений, заключающийся в следующем. Задача решается методом уступок. На каждой итерации ЛПР в визуальной форме на карте решений задает уступки для текущего критерия, что позволяет избежать основного недостатка метода уступок - отсутствие наглядности. Другой его недостаток - отсутствие адаптации к анализу критериев, представленных мультимодальными нелинейными функциями, - устраняется за счет

Критерий оценки метода

Анализ Анализ Схож- Просто-

Методы мультимо- унимо- дение та при-

дальных дальных метода менения

функций функций

Методы слу-

чайного поис- + + - +

ка

Метод каса- +

тельных

Градиентные методы - + + +

Методы уско-

ренного поиска вдоль греб- - + + -

ня

Сценарий ограничения субъективного участия _ЛПР в процессе принятия решения_

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

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

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

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

По сравнению с классическим, модернизированный МДЦ обладает следующими преимуществами:

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

- при зондировании пространства решений за счет применения градиентного метода происходит более жесткий отбор интересующих ЛПР точек;

- наблюдается снижение количества экспериментальных точек во множестве достижимых целей;

- разработан сценарий ограничения субъективного участия ЛПР в процессе решения задачи.

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

ОТКРЫТЫЕ ПРОГРАММНЫЕ СИСТЕМЫ С ПРИМЕНЕНИЕМ ГЕОИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ В ГОРНОЙ ПРОМЫШЛЕННОСТИ

В.М. Шек, Е.А. Конкин

Геоинформационные технологии (ГИС-тех-нологии) используются во всех отраслях народного хозяйства. Уровень применения геоинформационных программ в горном деле за последние тридцать лет серьезно повысился. Это привело к изменению самого характера деятельности горных компаний, а также к закономерному росту их производительности. Сегодня на мировом рынке коммерческих компьютерных программ для горных предприятий работают десятки фирм, предлагающих более 1000 программных продуктов, предназначенных для автоматизации различных функций управления горным производством.

Первая волна новых ГИС-технологий была связана с созданием простых моделей месторождений для оценки тоннажа и содержаний. Автома-

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

В начале 1970-х годов горная промышленность получила трехмерное цифровое блочное моделирование и геостатистический анализ минеральных ресурсов. Геологи научились использовать эти модели для прогнозирования запасов месторождений. Результатом стало улучшение качества их работы и достоверности оценки минеральных ресурсов. Эта волна нововведений была довольно продолжительной, ее действие закончилось в конце 1980-х годов.