MATHEMATICAL MODEL FOR DETERMINATION OF MASS AND GROSS EMISSIONS OF ATMOSPHERIC POLLLUTANTS OF THE INDUSTRIAL ENTERPRISES
Journal: VESTNIK OF ASTRAKHAN STATE TECHNICAL UNIVERSITY. SERIES: MANAGEMENT, COMPUTER SCIENCE AND INFORMATICS
Volume 2015 № 4
, 2015
Rubrics: MATHEMATICAL MODELING
Authors:
1.
Institute of Precision Mechanics and Control, Russian Academy of Sciences
2.
Institute of Precision Mechanics and Control, Russian Academy of Sciences
Type:
Article
Pages:
from 134
to 140
Status:
Published
Received:
05.10.2015
Accepted:
25.10.2015
Published:
25.10.2015
Subject area:
VAC 05.12.2013 Системы, сети и устройства телекоммуникаций
VAC 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
VAC 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
VAC 05.13.10 Управление в социальных и экономических системах
VAC 05.13.18 Математическое моделирование, численные методы и комплексы программ
VAC 05.13.19 Методы и системы защиты информации, информационная безопасность
UDK 574:517.977
GRNTI 20.01 Общие вопросы информатики
GRNTI 28.01 Общие вопросы кибернетики
GRNTI 49.01 Общие вопросы связи
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 82.01 Общие вопросы организации и управления
VAC 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
VAC 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
VAC 05.13.10 Управление в социальных и экономических системах
VAC 05.13.18 Математическое моделирование, численные методы и комплексы программ
VAC 05.13.19 Методы и системы защиты информации, информационная безопасность
UDK 574:517.977
GRNTI 20.01 Общие вопросы информатики
GRNTI 28.01 Общие вопросы кибернетики
GRNTI 49.01 Общие вопросы связи
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 82.01 Общие вопросы организации и управления
Language:
Russian
Keywords:
simulation model, emission of pollutants, industrial enterprise, automated control
Abstract (English):
In the course of functioning of the processing equipment of the industrial enterprise some tens pollutants causing harmful effects on health of the population, environment, agricultural crops, etc. come to the atmosphere. For minimization of the damage caused by atmospheric pollutants, the issue of the automated control of load distribution between the processing equipment of the enterprise is solved so that the size of damage do not exceed the preset value. For determination of concentration of pollutants on the controlled objects and territories, the models of advective transfer and diffusion caused by turbulent interaction of air masses are used. These models are based on the assumption that weight emission of pollutants in the atmosphere on all interval of the task solution remains invariable. However, in practice it significantly changes, depending on the changes of the structure and productivity of the working equipment, the production program, etc. The mathematical model, developed for determination of gross and mass emission of atmospheric pollutants of the industrial enterprise, helps determine more precisely the volume of the weight expense of pollutants and, thus, increase the accuracy of the solution of the problem of minimization of damage from atmospheric emissions of the industrial enterprise.
In the course of functioning of the processing equipment of the industrial enterprise some tens pollutants causing harmful effects on health of the population, environment, agricultural crops, etc. come to the atmosphere. For minimization of the damage caused by atmospheric pollutants, the issue of the automated control of load distribution between the processing equipment of the enterprise is solved so that the size of damage do not exceed the preset value. For determination of concentration of pollutants on the controlled objects and territories, the models of advective transfer and diffusion caused by turbulent interaction of air masses are used. These models are based on the assumption that weight emission of pollutants in the atmosphere on all interval of the task solution remains invariable. However, in practice it significantly changes, depending on the changes of the structure and productivity of the working equipment, the production program, etc. The mathematical model, developed for determination of gross and mass emission of atmospheric pollutants of the industrial enterprise, helps determine more precisely the volume of the weight expense of pollutants and, thus, increase the accuracy of the solution of the problem of minimization of damage from atmospheric emissions of the industrial enterprise.
Keywords:
simulation model, emission of pollutants, industrial enterprise, automated control
simulation model, emission of pollutants, industrial enterprise, automated control
Введение В системах объектового управления Министерства Российской Федерации по делам гражданской обороны, чрезвычайным ситуациям и ликвидации последствий стихийных бедствий (МЧС России) и экологического мониторинга предприятий при математическом моделировании процесса загрязнения объектов и территорий атмосферными поллютантами широкое распространение получили формальные модели адвективного переноса и диффузии, вызванной турбулентным взаимодействием воздушных масс [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. Marchuk G. I. Matematicheskoe modelirovanie v probleme okruzhayuschey sredy / G. I. Marchuk. M.: Nauka, 1982. 320 s.
2. Kvashnin I. M. Promyshlennye vybrosy v atmosferu. Inzhenernye raschety i inventarizaciya / I. M. Kvashnin. M.: AVOK-PRESS, 2005. 390 s.
3. Kushnikova E. V. Imitacionnoe modelirovanie zagryazneniya territoriy vybrosami promyshlennyh predpriyatiy / E. V. Kushnikova, O. A. Toropova // Problemy upravleniya v social'no-ekonomicheskih i tehnicheskih sistemah: cb nauch. st. po materialam Vseros. nauch. konf. Saratov, 2012. S. 20-23.
4. Kushnikov V. A. Arhitektura prikladnogo programmnogo obespecheniya dlya formal'nogo analiza svoystv celey i sinteza kriteriev upravleniya slozhnymi social'nymi i ekonomicheskimi sistemami / V. A. Kushnikov, E. V. Kushnikova // Vestn. Saratov. gos. tehn. un-ta. 2009. № 4 (43), vyp. 2. S. 199-201.
5. Avetisyan Yu. A. Matematicheskie modeli i algoritmy operativnogo upravleniya processami likvidacii chrezvychaynyh situaciy / Yu. A. Avetisyan, V. A. Kushnikov, A. F. Rezchikov, V. A. Rodichev // Mehatronika, avtomatizaciya, upravlenie. 2009. № 11. S. 43-47.
6. Solyanik N. A. Matematicheskoe modelirovanie processa zagryazneniya atmosfernogo vozduha v zone vliyaniya promyshlennyh predpriyatiy / N. A. Solyanik, V. A. Kushnikov // Vestn. Saratov. gos. tehn. un-ta. 2009. №. 1 (37). S. 104-108.
7. Kushnikov V. A. Upravlenie v cheloveko-mashinnyh sistemah s avtomatizirovannoy proceduroy korrekcii celey / V. A. Kushnikov, A. F. Rezchikov, A. D. Cvirkun // Avtomatika i telemehanika. 1998. № 7. S. 168-175.
8. Pshenichnikov I. S. Analiz vypolnimosti planov meropriyatiy v sisteme avtomatizirovannogo upravleniya mostostroitel'noy organizaciey / I. S. Pshenichnikov, V. A. Kushnikov, E. I. Shlychkov, D. F. Rezchikov // Mehatronika, avtomatizaciya, upravlenie. 2006. № 11. S. 45-49.
