В работе предлагается модификация численного метода фейерверков однокритериальной оптимизации для решения задач многокритериальной оптимизации. Метод относится к метаэвристическим алгоритмам, он не гарантирует нахождения точного решения, но может найти достаточно хорошее приближенное решение. Рассматриваются многокритериальные задачи оптимизации с числовыми критериями, имеющими одинаковую важность. Допустимое решение задачи представляется вектором из действительных чисел, значение каждой компоненты которого принадлежит определенному отрезку. Под оптимальным решением понимается решение, оптимальное по Парето. Так как точных решений, оптимальных по Парето, может быть бесконечно много, рассматривается способ нахождения приближения, состоящего из конечного числа решений, оптимальных по Парето. Модификация основана на процедуре недоминируемой сортировки, которая является основной процедурой для управления процессом поиска приближенного решения. Недоминируемая сортировка – это ранжирование решений на основе значений компонент числового вектора, полученных с помощью вычисления критериев. Каждая компонента соответствует определенному критерию, а множество решений разбивается на непересекающиеся подмножества. Первое подмножество – это решения, оптимальные по Парето, второе подмножество – это решения, оптимальные по Парето, если не учитывать первое подмножество, последнее подмножество – это решения, оптимальные по Парето, если не учитывать все предыдущие подмножества. После такого разбиения принимается решение о генерировании новых допустимых решений. Работа метода протестирована на общеизвестных задачах многокритериальной оптимизации с двумя критериями: ZDT2, LZ01. Задачи отличаются структурой расположения решений, оптимальных по Парето. Так LZ01 имеет достаточно сложную структуру решений, оптимальных по Парето. В заключении обсуждаются вопросы о дальнейшем направлении исследований и о возможности модификации метода для задач многокритериальной оптимизации с произвольными, а не параллелепипедными ограничениями.