МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДЛЯ ОПРЕДЕЛЕНИЯ МАССОВОГО И ВАЛОВОГО ВЫБРОСА АТМОСФЕРНЫХ ПОЛЛЮТАНТОВ ПРОМЫШЛЕННОГО ПРЕДПРИЯТИЯ
Аннотация и ключевые слова
Аннотация (русский):
В процессе функционирования технологического оборудования промышленного предприятия в атмосферу поступает несколько десятков поллютантов, оказывающих вредное воздействие на здоровье населения, окружающую среду, сельскохозяйственные посевы и т. д. Для определения концентрации данных загрязнителей на контролируемых объектах и территориях используются модели адвективного переноса и диффузии, вызванной турбулентным взаимодействием воздушных масс. В основу этих моделей положено допущение, что весовой выброс поллютантов в атмосферу на всем интервале решения задачи остается неизменным. Однако на практике весовой выброс существенно трансформируется в зависимости от изменений состава и производительности работающего оборудования, производственной программы и т. д. Математическая модель, разработанная для определения валового и массового выброса атмосферных поллютантов промышленного предприятия, позволяет более точно определить величину весового расхода загрязнителей и, благодаря этому, существенно повысить точность прогнозирования степени загрязнения контролируемых объектов и территорий атмосферными выбросами промышленных предприятий.

Ключевые слова:
имитационная модель, выброс загрязнителей, промышленное предприятие, автоматизированное управление
Текст
Введение В системах объектового управления Министерства Российской Федерации по делам гражданской обороны, чрезвычайным ситуациям и ликвидации последствий стихийных бедствий (МЧС России) и экологического мониторинга предприятий при математическом моделировании процесса загрязнения объектов и территорий атмосферными поллютантами широкое распространение получили формальные модели адвективного переноса и диффузии, вызванной турбулентным взаимодействием воздушных масс [1]. В основу этих моделей положено допущение, что величины массового и валового выброса атмосферных загрязнителей промышленного предприятия на интервале моделирования остаются постоянными, хотя в действительности они изменяются во времени в зависимости от состава и производительности технологического оборудования, его износа и интенсивности эксплуатации, особенностей производственной программы и т. д. Использование данного допущения обусловливают, на наш взгляд, сложность аналитического описания процесса формирования атмосферных выбросов технологическим оборудованием, случайный характер включения или отключения оборудования на интервале моделирования, различная продолжительность и эффективность работы воздушных фильтров и др. Один из способов повышения точности прогнозирования загрязненности контролируемых объектов и территорий атмосферными выбросами связан с использованием численных имитационных моделей, описывающих процесс формирования и переноса атмосферных поллютантов. Это обстоятельство обусловливает практическую направленность цели нашего исследования - разработка имитационной модели для определения величины массового и валового выброса загрязнителей атмосферы промышленного предприятия. Постановка задачи Целью нашего исследования стала разработка имитационной модели процесса формирования атмосферных поллютантов технологическим оборудованием промышленного предприятия для определения величины массового и валового выброса его загрязнителей атмосферы. Модель должна быть ориентирована на использование в составе информационных систем МЧС объектового уровня, а также в системах экологического мониторинга промышленного предприятия. Для достижения цели исследования необходимо было решить следующие задачи: - выполнить классификацию производственного оборудования по способу определения величины массового и валового выброса загрязнителей атмосферы, формируемых в процессе его функционирования; - разработать математическую модель для определения валового и массового выброса атмосферных поллютантов по каждой группе технологического оборудования; - предложить и обосновать алгоритм практического использования разработанной математической модели. Математическая модель Модель позволяет определить массовую и валовую величину атмосферных выбросов в зависимости от типа предприятия, количества и вида используемого оборудования, а также от продолжительности и интенсивности его работы. Входной информацией для модели являются принадлежность предприятия к той или иной отрасли народного хозяйства; список вредных веществ, выбрасываемых в атмосферу в процессе функционирования; производственное оборудование, при работе которого образуются вредные вещества; смены и дни работы производственного оборудования; законы распределения дискретных случайных величин, характеризующих моменты его включения и выключения. Выходными переменными математической модели служат величины массового и валового выброса загрязняющих веществ . Для разработки математической модели предварительно необходимо осуществить классификацию производственного оборудования предприятия по способу расчета величины атмосферных выбросов. Для этого все производственное оборудование, в процессе функционирования которого образуются атмосферные поллютанты, декомпозируется на 10 групп. Для каждой из этих групп при расчете выбросов используются: - одна из характеристик оборудования, например марка станка; - удельных эмиссий загрязняющих веществ, приходящихся на единицу массы, длины или объема используемого материала; - балансовый метод; - известная интенсивность испарения поллютантов с единицы поверхности, например гальванической ванны; - методика инвентаризации выбросов, образующихся при производстве радиоэлектронной аппаратуры; - методика инвентаризации выбросов, образующихся при ремонте резинотехнических изделий; - способ инвентаризации выбросов, образующихся при проведении медницких работ; - способ инвентаризации выбросов, образующихся при мойке деталей, узлов и агрегатов; - алгоритм определения выбросов пыли во время механической обработки древесины; - методика определения выбросов пыли во время перегрузки сыпучих материалов. Указанные методики были систематизированы в [2], где они использовались при решении задачи инвентаризации промышленных выбросов в атмосферу. Адаптация этих методик к требованиям разрабатываемой системы имитационного моделирования, в основу которой положен программный продукт Anylogic [3], осуществлена ниже. В частности, при расчете и сделано допущение, что поллютанты от источников загрязнений попадают в помещения, оснащенные газопылеуловителями, работающими без утечек воздуха. Удаление поллютантов производится с помощью местной и общеобменной вентиляции. Поясним основные этапы формирования имитационной модели на примере оборудования первой и второй группы. Определение массового и валового выброса для первой группы оборудования Расчет ведется по одной из характеристик оборудования. Данный метод расчета величины массового и валового выброса применяется, например, при механической обработке металлов. Определение выполняется с использованием справочных данных, приведенных в соответствующих отраслевых методиках [2]. При расчете и для первой группы оборудования используются следующие зависимости: ; ; ; где - массовое количество i-го поллютанта, попадающего в газопылеуловители местной вентиляционной системы от работающих источников загрязнений первой группы оборудования; - массовое количество i-го поллютанта, поступившего от работающего оборудования; - коэффициент эффективности местных вытяжных устройств; - количество поллютантов; - коэффициент производительности j-й единицы работающего оборудования первой группы; - степень очистки воздуха газопылеуловителями; - концентрация i-го поллютанта входе газопылеуловителей; - концентрация i-го поллютанта на выходе газопылеуловителей; - массовое количество i-го поллютанта, не уловленное местными вытяжными устройствами и удаляемое общеобменной вытяжной вентиляцией; l - количество интервалов времени Δtj, в течение которых работало оборудование первой группы за год; - количество загрязнителей, выбрасываемых в атмосферу j-й единицей оборудования первой группы; символ «1» в верхнем индексе означает принадлежность к оборудованию первой группы. В момент решения задачи работа j-й единицы оборудования первой группы определяется по функции распределения непрерывной случайной величины , характеризующей вероятность работы каждой единицы оборудования первой группы на временном интервале решения задачи. Данная зависимость определяется путем обработки экспериментальных данных на стадии адаптации имитационной системы к конкретному производству. Суммирование единиц работающего оборудования позволяет найти Определение массового и валового выброса для второй группы оборудования К этой группе относится оборудование сварочного и литейного производства, термообработки металлов, производства деталей способом порошковой металлургии и др. Расчет ведется по удельным эмиссиям загрязняющих веществ, приходящихся на единицу массы, длины или объема используемого материала. При этом для определения и используются следующие зависимости: ; ; ; ; ; где - массовое количество i-го поллютанта, попадающего в газопылеуловители местной вентиляционной системы; - массовое количество i-го поллютанта, поступившего от работающего оборудования второй группы; - коэффициент эффективности местных вытяжных устройств; - степень очистки воздуха газопылеуловителями; - концентрация i-го поллютанта на входе газопылеуловителей; - концентрация i-го поллютанта на выходе газопылеуловителей; - коэффициент производительности j-й единицы работающего оборудования второй группы; - массовое количество i-го поллютанта, не уловленное местными вытяжными устройствами и удаляемое общеобменной вытяжной вентиляцией; l - количество интервалов времени ; - величина удельного выделения атмосферного загрязнителя на единицу массы (длины, или объема) j-го используемого материала; - количество j-го материала, перерабатываемого в момент решения задачи t = tp; - масштабирующий коэффициент, используемый для согласования размерностей; - количество перерабатываемых материалов; символ «2» в верхнем индексе означает принадлежность к оборудованию второй группы. При t = tp величина определяется по функции распределения непрерывной случайной величины характеризующей вероятность переработки j-го материала на временном интервале решения задачи. Аналогично были построены блоки имитационной модели, определяющие величину массового и валового выброса для остальных восьми групп оборудования. Для подтверждения возможности практического использования разработанной модели на рисунке приведены результаты вычислительного эксперимента по определению веса поллютантов, поступающих в атмосферу при эксплуатации оборудования второй группы. г Результаты вычислительного эксперимента по определению веса поллютантов, поступающих в атмосферу в течение суток от участка ручной электродуговой сварки штучными электродами Участок оснащен шестью сварочными постами, являющимися источниками загрязнений: ИВ1, ИВ2, …, ИВ6 соответственно. При сварке использовались электроды ОЗС-4 и УОНИ-13/45. Часовой расход электродов первой марки составлял 1,8 кг, второй - 1,5 кг. В верхней части рисунка пунктиром показаны сварочные посты, отключаемые на соответствующем временном интервале, закрашенные прямоугольники обозначают включаемые посты. До начала моделирования сварка не производилась. Результаты вычислительного эксперимента, выполненного для оборудования второй группы, подтверждают наличие существенной неравномерности массового выброса атмосферных поллютантов промышленного предприятия на временном интервале одни сутки. Анализ полученных результатов позволяет рекомендовать разработанную модель к использованию в автоматизированных системах экологического контроля и мониторинга, т. к. в них при проведении расчетов часто используются усредненные по месяцам или годам значения величин массового и валового выброса, что снижает точность прогнозирования. Алгоритм моделирования Общий алгоритм использования разработанной имитационной модели при определении массового и валового выброса атмосферных поллютантов имеет следующий вид. Начало алгоритма. 1. Определение загрязняющих веществ , для которых необходимо определить величину массового и валового выброса. 2. Выбор интервала моделирования, на котором будут определяться величины массового и валового выброса загрязняющих веществ , а также шага моделирования. 3. Определение наименований и количества технологического оборудования, в процессе функционирования которого образуются загрязняющие вещества . 4. Отнесение каждой единицы используемого технологического оборудования к соответствующей классификационной группе. 5. Выбор соответствующей имитационной модели для расчета величины массового и валового выброса по каждой группе оборудования. 6. Определение функций распределения непрерывных случайных величин, характеризующих вероятность работы каждой единицы оборудования на временном интервале решения задачи. 7. Расчет величин массового и валового выброса загрязняющего вещества для i-й группы оборудования. 8. Определение величин массового и валового выброса загрязняющего вещества по всему предприятию на выбранном шаге моделирования и на всем интервале прогнозирования. 9. Выдача сообщения о величинах массового и валового выброса загрязняющих веществ на выбранном интервале моделирования. Конец алгоритма. Разработанная имитационная модель была использована для определения массового и валового выброса атмосферных загрязнителей предприятия ОАО «СЭПО-ЗЭМ» (г. Саратов), что позволило на интервалах времени 1 месяц и 1 год уменьшить максимальную величину выбросов поллютантов за счет изменения графика работы технологического оборудования. Результаты моделирования и сама модель использовались также при модернизации системы автоматизированного управления производственными процессами ОАО «СЭПО-ЗЭМ». Адаптация имитационной модели к работе в составе информационных систем данного предприятия выполнялась в соответствии с рекомендациями [4-8]. Заключение Предложенная имитационная модель процесса образования атмосферных загрязнителей позволяет, в зависимости от используемого технологического оборудования и производственной программы, в динамике определить величину массового и валового выброса различных видов поллютантов, образующихся при функционировании промышленного предприятия. Использование данной модели в составе автоматизированных систем экологического мониторинга дает возможность на различных временных интервалах повысить точность прогнозирования загрязнения контролируемых объектов и территорий атмосферными выбросами промышленных предприятий.
Список литературы

1. Марчук Г. И. Математическое моделирование в проблеме окружающей среды / Г. И. Марчук. М.: Наука, 1982. 320 с.

2. Квашнин И. М. Промышленные выбросы в атмосферу. Инженерные расчеты и инвентаризация / И. М. Квашнин. М.: АВОК-PRESS, 2005. 390 с.

3. Кушникова Е. В. Имитационное моделирование загрязнения территорий выбросами промышленных предприятий / Е. В. Кушникова, О. А. Торопова // Проблемы управления в социально-экономических и технических системах: cб науч. ст. по материалам Всерос. науч. конф. Саратов, 2012. С. 20-23.

4. Кушников В. А. Архитектура прикладного программного обеспечения для формального анализа свойств целей и синтеза критериев управления сложными социальными и экономическими системами / В. А. Кушников, Е. В. Кушникова // Вестн. Саратов. гос. техн. ун-та. 2009. № 4 (43), вып. 2. С. 199-201.

5. Аветисян Ю. А. Математические модели и алгоритмы оперативного управления процессами ликвидации чрезвычайных ситуаций / Ю. А. Аветисян, В. А. Кушников, А. Ф. Резчиков, В. А. Родичев // Мехатроника, автоматизация, управление. 2009. № 11. С. 43-47.

6. Соляник Н. А. Математическое моделирование процесса загрязнения атмосферного воздуха в зоне влияния промышленных предприятий / Н. А. Соляник, В. А. Кушников // Вестн. Саратов. гос. техн. ун-та. 2009. №. 1 (37). С. 104-108.

7. Кушников В. А. Управление в человеко-машинных системах с автоматизированной процедурой коррекции целей / В. А. Кушников, А. Ф. Резчиков, А. Д. Цвиркун // Автоматика и телемеханика. 1998. № 7. С. 168-175.

8. Пшеничников И. С. Анализ выполнимости планов мероприятий в системе автоматизированного управления мостостроительной организацией / И. С. Пшеничников, В. А. Кушников, Е. И. Шлычков, Д. Ф. Резчиков // Мехатроника, автоматизация, управление. 2006. № 11. С. 45-49.