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

004610337 На правах рукописи

ГЕРБЕР ЮЛИЯ ВАЛЕРЬЕВНА

МЕТОД КОНТРОЛЯ И ПРОГНОЗИРОВАНИЯ РАСПРОСТРАНЕНИЯ ПРОМЫШЛЕННЫХ ГАЗОВЫХ ВЫБРОСОВ В ГОРОДАХ НА ПРИМЕРЕ НОВОМОСКОВСКА

05.11.13 «Приборы и методы контроля природной среды, веществ, материалов и изделий»

Автореферат диссертации на соискание учёной степени кандидата технических наук

Москва-2010

1 4 ОКТ 20Ю

004610337

Работа выполнена на кафедре «Автоматизация производственных процессов» Новомосковского института Российского химико-технологического университета им. Д.И.Менделеева и на кафедре «Мониторинг и автоматизированные системы контроля» Московского государственного университета инженерной экологии.

Научный руководитель:

доктор технических наук, профессор Беляев Юрий Иванович

Официальные оппоненты: доктор технических наук, профессор

Бирюков Валентин Васильевич, кандидат технических наук, доцент Лопатин Александр Геннадьевич

Ведущая организация:

НПО «Химавтоматика», г. Москва

Защита состоится 28 октября 2010 г. в 1400 часов на заседании диссертационного совета Д 212.145.02 Московского государственного университета инженерной экологии по адресу: 105066, г. Москва, ул. Старая Басманная, д. 21/4, аудитория В-23.

С диссертацией можно ознакомиться в библиотеке университета.

Отзывы в двух экземплярах заверенные печатью, просим направлять по адресу: 105066 г. Москва, ул. Старая Басманная, д. 21/4.

Автореферат разослан ___ сентября 2010 г.

Учёный секретарь диссертационного совета к.т.н., доцент

Н.В. Мокрова

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность работы. Стремительное развитие научно-технического прогресса и образование городов привело к появлению новых проблем в развитии биосферы, непосредственно связанной с нарушением экологического баланса среды обитания человека.

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

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

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

Для достижения поставленной цели необходимо решить следующие задачи:

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

-разработка методики составления карты прогноза экологического состояния воздушной среды города;

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

- разработка математической модели распространения облака газа в приземном слое атмосферы;

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

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

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

- разработан алгоритм и программа прогнозирования распространения облака газа в виде карты экологической обстановки воздушной среды города;

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

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

- предложена методика расчёта распространения облака газа по территории города с учётом ландшафта местности;

- предложена методика составления картограмм линий тока движения воздуха, как на местности, так и в жилых районах по экспериментальным данным.

Практическая значимость. Результаты работы предлагается использовать при создании систем экологического мониторинга воздушной среды города для служб ропотребнадзора, МЧС, метеорологии.

Реализация результатов. Результаты теоретических и экспериментальных исследований использовались территориальным управлением роспотребнадзора Новомосковского района, г. Донского, г. Ки-мовска и Узловского района при подготовке ежегодного муниципального доклада о санитарной эпидемиологической обстановке в муниципальном образовании г. Новомосковска.

Апробация и публикации. По теме диссертационной работы опубликовано девять печатных работ, в том числе две в журналах, рекомендованных ВАК.

Результаты работы докладывались и обсуждались на конференциях: Города России: проблемы строительства, инженерного обеспечения, благоустройства и экологии. (Пенза, 2009 г.), Окружающая среда и здоровье. (Пенза, 2009 г.), РХТУ им. Д.И.Менделеева, Новомосковский институт, 2008 -2010 гг.

Структура и объём работы. Диссертационная работа состоит из введения, четырёх глав, заключения, списка использованных источников и пяти приложений. Работа изложена на 145 страницах машинописного текста, содержит 7 таблиц и 40 рисунков. Список использованных источников включает 130 наименований.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

Во второй главе излагается разработанный метод контроля и прогнозирования распространения промышленных газовых выбросов в городах, который заключается в том, что составляется математическая модель и на её основе рассчитывается характеристика экологической обстановки (рис. 1).

Топотра- Физические

фическая свойства Координаты

карта параметров источника

города выброса выброса

-i-.......Ч-

Входные данные

Реперные высоты местности Метеоусловия

i

Математическая модель распространения облака газа

т

Карта прогноза экологической обстановки воздушной среды города

Картограмма линий тока движения воздуха

Рис. 1. Метод контроля и прогнозирования воздушной среды города

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

- уравнение ландшафта местности задаётся в виде ряда тригонометрических полиномов:

т-1

2(х,у) = соэ( ж рхкх)соз(ж руку) (1)

ыо

где х=(0,... ,1) и >>=((),... ,1); т - число гармоник; Ак - коэффициент к-ой гармоники; рхК и рук - соответственно порядок гармоники по координатам х и у.

- система уравнений турбулентности приземного слоя атмосферы: [¿фг.у)

*"М (2)

<ъ>

где с(х, у) - концентрация газа; х, у - координаты точки местности; Ьх,Ьу - размеры участка местности; Ух(х,у),Уу(х,у)~ проекция скоростей на оси координат;

- уравнение конвективного переноса газа:

(3)

«г

где к - коэффициент; р^ш, /%озд - плотность газа и воздуха.

- уравнение диффузии газа:

е = (4)

где/) - коэффициент диффузии газа.

Начальное условие распространения облака газа:

во (хист - Лет)= производит ельность источника выброса , (5)

где (*ист,уИС1)~ координаты источника выброса.

Граничные условия остановки распространения облака газа: 14с, ПРИ К(х,У) = О

„при Уу(х,у) = О' (6)

где ДСост, уос1 - координаты остановки облака газа.

Время остановки распространения облака газа равно моменту когда с(х, у) < спдк, спдк - предельно допустимая концентрация.

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

В расчётных соотношениях для учета задаются вид газа коэффициент диффузйи и отношение плотности газа к плотности возд>"ха при данных метеоусловиях.

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

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

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

В вертикальном направлении на частицу действует сила конвекции, которая пропорциональна разности плотностей среды и примеси: В = к(р0-р) (7)

где р0 - плотность среды, кг / мъ; р - плотность примеси, кг / м3; к ~ коэффициент.

Моделирование процесса диффузии проводится по зависимостям, основанным на первом законе Фика. На рис. 2 представлена расчётная схема перемещения частицы газа в одной элементарной ячейке трёхмерной модели.

Механизм переноса газа в элементарной ячейке учитывается следующим образом: считается, что в ячейках п1, п2, п5, пб перенос осуществляется под воздействием турбулентного переноса и диффузии, действующей в горизонтальном направлении; в ячейки пЗ, п4 перенос осуществляется под воздействием конвективной силы и диффузии, действующей в вертикальном направлении.

Рис. 2. Расчётная схема перемещения частицы газа в одной элементарной ячейке трёхмерной модели

Расчётное соотношение изменения концентрации в текущем объёме запишется в виде:

An = С-и1 - п2) + (иЗ - и4) + («5 - иб) = D + D^^ + D^Jlfl +

Ах Ах Az ,ол

_ _ _ (о)

+Д +D £I_£L + D £s_£°.. Az Ay Ay

Для расчётной модели турбулентности воздушного потока исходные данные, полученные в реперных узлах в виде значений скоростей V(x, у) и направлений ветра, а(х, у) преобразуются в две ортогональные горизонтальные проекции векторов скоростей Ау(х, у) и Ах(х, у):

ГАх(х,у) = V(x,y)cos(a(x,y)) [ЛуО, У) = ^(x,y)sin(a(x,y))' которые задаются двумя интерполяционными уравнениями:

М-1

£3 , (Ю)

Af{x, у) = &! ксрк (ж, у), ^ е [0... Lv j

где у) - тригонометрические полиномы; ахк,а у- оценка амплитуды для к-й гармоники в соответствующей моделируемой проекции.

На основе расчётов модели прогноза распространения облака газа формируется карта прогноза экологической обстановки города. В качестве примера на рис. 3 изображён фрагмент карты прогноза экологической обстановки воздушной среды г. Новомосковска.

Карта прогноза экологической обстановки воздушной среды г. Новомосковска

Газ аммиак

Ветер северо-западный

Иш¡8!Щi!LШйШ¿¿¿&

Рис. 3. Фрагмент карты прогноза экологической обстановки воздушной среды г. Новомосковска

На ней указан источник выброса аммиака, изображены промышленная и жилая зоны. Приведены расчёты распространения облака аммиака при северо-западном направлении ветра при разных баллах шкалы Бофорта. Там же указано время остановки движения облака газа. Представленная и разработанная карта прогноза даёт достаточно полную информацию о степени опасности выброса газа.

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

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

Для предварительного анализа степени опасности рассматривались наиболее неблагоприятные ситуации по направлению попадания облака газа на жилые районы города. Эта ситуация отражена в документе, который назван картой прогноза экологической обстановки жилого района города. Расчёты распространения облака газа по жилому району ведутся до полной остановки движения облака газа. На рис.4 представлена карта прогноза распространения облака газа по одной из улицы города. При указанном направлении можно наблюдать распространение ССЬ из автомобильных выбросов на жилой район города.

.........

Рис.4 Карта прогноза распространения облака газа по улицы города

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

1) что имеется п детекторов, определённым способом размещённых на отрезке от [0,Ь] в узлах {х,}.

2) что ка'ндое наблюдение можно представить в виде

= Ж)+ #(*,), «=1(11)

где }!(х[) - действительное значение сигнала детектора поля в точке х. аддитивная погрешность измерения в точке xi со стати-

стическими характеристиками.

М4 (л,) = 0 У»с1,и

где М - оператор математического ожидания.

Необходимо на основе таких исходных данных построить интерполяционное уравнение вида:

у(х)=Хат(х) х £ [о,ь], (13)

где - искомые коэффициенты; щ (х) - п первых функций полной в Ь2 системе полиномов, которое даёт наилучшее в определённом смысле приближение к истинной функции у{х) на фоне воздействия шума в узлах измерения.

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

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

к, номер гармоники Ак, коэффициент гармоники рхК порядок гармоники по координате X РУк порядок гармоники по координате у

0 217,11 0 0

1 9,89 2 0

2 -12.04 0 1

3 -9,99 1 1

4 -5,44 0

5 -4,90 0 8

6 -6,05 3 3

7 594 5 8

-5,15 7 5

9 -3,64 6 0

10 -3,50 4 0

11 7,64 3 1

12 -2,73 0 5

13 5,39 5 4

14 4,99 5 3

15 -4,00 6 6

В результате проведенных расчётов была получена модель ландшафта исследуемой местности. Графическое изображение результатов моделирования представлено на рис. 5.

vws»

X

,х

Рис. 5. Изображение рельефа местности Новомосковского района

В этой же главе были рассмотрены влияния различных факторов на погрешность прогнозирования распространения облака газа. К ним

относятся: погрешность метеоданных (скорость и направление ветра); погрешность задания коэффициента диффузии; влияние неоднородности поверхности района.

1. Влияние погрешность метеоданных. Исследование погрешности расчёта обусловлено неточностью задания направления движения ветра. Для этого рассмотрим две ситуации, когда из одного источника облако газа распространяется в двух направлениях смещенных на угол А^(рис. 6).

Ксточнкк выбросы

Ь - расстояние от источника выброса до фронта облака; Аср - угол отклонения направления ветра;

Д51... ДУп - отклонение облака; Д5=_ДТ Аср) при малых углах

Аср.

Распространение облака газа в двух направлениях, смещённых на угол Аср. Отклонение облака газа можно рассчитать по формуле:

Рис.1

Д!>=£ —Дя>, V

(14)

скорость ветра;Д^ -

где к - коэффициент; Ь - длина облака газа; V -угол отклонения направления ветра.

2. Влияние погрешности задания коэффициента диффузии на определение геометрических размеров облака. Одним из главных параметров облака газа после выброса является его геометрические размеры, которые зависят от скорости ветра, коэффициента диффузии газа и времени от начала выброса.

Площадь облака рассчитывается по формуле :

Б = (15)

где - площадь облака при V = 0, в случае, когда распространение происходит за счёт диффузии

5\=Ю, (16)

где О - диффузия; к - коэффициент; V - скорость ветра; Т- время от начала выброса.

Погрешность расчета ширины облака рассчитывается по формуле:

АН = АтШ. (17)

где к - коэффициент; АО - погрешность отклонения задания коэффициента диффузии.

3. Влияние неоднородности поверхности района. На рис. 7 изображены эпюры распространения облака хлора по району города без учёта рельефа (рис. 7 а) и с учётом рельефа местности (рис. 76).

б

Рис. 7. Эпюры распространения облака хлора Длина и ширина облака получились в обоих случаях 5x1 км и 3,5x0,8 км. Погрешность определения составила по ширине 200 м, по длине соответственно 1500 м.

В результате сопоставительного анализа было выявлено, что основными источниками погрешности прогнозирования являются метеоданные (скорость и направление движения ветра).

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ РАБОТЫ

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

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

3. Предложена структура карты прогноза экологической обстановки воздушной среды города.

4. Предложена методика составления карты прогноза экологической обстановки воздушной среды жилого района города.

5.Составлены карты прогноза экологической обстановки воздушной среды города Новомосковска и его отдельных характерных жилых районов.

6.Разработана программа «Прогноз» для прогнозирования распространения облака газа по ландшафту местности.

7. Результаты диссертационной работы использованы территориальным управлением роспотребнадзора Новомосковского района, г. Донского, г. Кимовска и Узловского района при подготовке ежегодного муниципального доклада о санитарной эпидемиологической обстановке в муниципальном образовании г. Новомосковска.

Список работ, в которых отражены основные результаты:

1. Беляев Ю.И., Гербер Ю.В. Система мониторинга распространения газообразных веществ в атмосфере города /«Датчики и системы», №10(125), Москва, 2009,- С. 13 - 16.

2. Беляев Ю.И., Вент Д.П., Гербер Ю.В., Вепренцева О.Н., Латышенко К.П. Экологический мониторинг атмосферы города/ «Химическое и нефтегазовое машиностроение», №8, Москва, 2009 .-С. 32 - 34.

3.Вент Д.П., Беляев Ю.И., Гербер Ю.В., Вепренцева О.Н. Модели распространения загрязнений атмосферы города Информатика, Экология, Экономика. Вестник международной академии системных исследований (МАСИ). М.: 2007. С. 26-28.

4. Беляев Ю.И., Гербер Ю.В., Вепренцева О.Н., Методика оперативного прогнозирования экологического состояния атмосферы города. Информатика, Экология, Экономика. Вестник международной академии системных исследований (МАСИ). М.: 2007. С. 26 - 28.

5. Беляев Ю.И., Гербер Ю.В., Вепренцева О.Н.. Система мониторинга атмосферы города IX научно-техническая конференция молодых ученых, аспирантов, студентов; Тезисы докладов. / РХТУ им. Д.И. Менделеева, Новомосковский институт. Новомосковск, 2008. С. 182 - 183.

6. Беляев Ю.И., Гербер Ю.В., Бизина С.А. Оперативное прогнозирование состояния атмосферы города, XXVII научная конференция про-фессорско -преподавательского состава и сотрудников Тезисы докладов/НИ РХТУ им. Д.И. Менделеева. Ч. 1. Новомосковск, 2009. - С. 142-143.

7. Беляев Ю.И., Гербер Ю.В. Информационные и компьютерные технологии прогнозирования состояния атмосферы города. XI международная научно-практическая конференция. Города России: проблемы строительства, инженерного обеспечения, благоустройства и экологии. МК -14-9, Пенза, 2009. - С. 27 - 30.

8. Беляев Ю.И., Гербер Ю.В., Микаилов Р.И. Основные принципы моделирования атмосферы города как объекта экологического мониторинга. VI Международная научно-практическая конференция. Окружающая среда и здоровье. МК-/4-9, Пенза, 2009. - С. 12 - 15.

9.Беляев Ю.И., Гербер Ю.В., Моисеева Й.Д. Бизина СЛ. Оперативное прогнозирование состояния атмосферы города. Успехи в химии и химической технологии: сб. науч. Тр. Том XXII, № 7(101). - М. : РХТУ им. Д.И. Менделеева, 2009. - С. 122-123.

Заказ № 943/777_Объем 1.0 п. л._Тираж 100 экз.

ГОУ ВПО «РХТУ им. Д.И. Менделеева» Новомосковский институт (филиал) Издательский центр

ВВЕДЕНИЕ

1 АНАЛИЗ МЕТОДОВ ЭКОЛОГИЧЕСКОГО МОНИТОРИНГА

1.1 Общие сведения о методах

1.1.3 Характеристики источников выбросов

1.1.4 Методы оценки дисперсии

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

1.2 Классификация существующих методов и моделей

1.3 Штатная модель служб ГО и возможности её совершенствования

1.4 Модель Паскуилла-Гиффорда

1.5 Модель Института экспериментальной метеорологии

1.6 Трехмерные модели переноса и диффузии примеси и их упрощенные варианты 39 Выводы

2 МЕТОД КОНТРОЛЯ И ПРОГНОЗИРОВАНИЯ РАСПРОСТРАНЕНИЯ ПРОМЫШЛЕННЫХ ГАЗОВЫХ ВЫБРОСОВ

В ГОРОДАХ

2.1 Метод контроля и прогнозирования распространения промышленных газовых выбросов в городах

2.2 Составление математической модели распространения облака газа

2.3 Моделирование турбулентности приземного слоя атмосферы

2.4 Моделирование конвективного переноса газа

2.5 Термодинамические потоки

2.6 Описание явлений переноса в газах

2.7 Уравнение конвективного переноса газа

2.8 Моделирование процесса диффузии

2.9 Моделирование рельефа местности с помощью помехоустойчивой интерполяции

2.10 Построение картограмм линий тока движения воздуха

2.11 Карта прогноза экологической обстановки воздушной среды города

Выводы

3 МЕТОД КОНТРОЛЯ И ПРОГНОЗИРОВАНИЯ СОСТОЯНИЯ ВОЗДУШНОЙ СРЕДЫ ЖИЛОГО РАЙОНА ГОРОДА «

3.1 Методика составления карты прогноза экологической обстановки жилого района города

3.2 Экспериментальное исследование жилого района и построение картограмм линий тока движения воздуха

3.3 Рассмотрение неблагоприятных ситуаций по направлению попадания облака газа на жилой район

Выводы

4 ОСНОВНЫЕ СОСТАВЛЯЮЩИЕ МЕТОДА КОНТРОЛЯ И ПРОГНОЗИРОВАНИЯ РАСПРОСТРАНЕНИЯ ПРОМЫШЛЕННЫХ ГАЗОВЫХ ВЫБРОСОВ В ГОРОДАХ

4.1 Основные положения помехоустойчивой интерполяции

4.2 Описание интерфейса и настроек программы «Прогноз»

4.3 Построение математической модели рельефа местности

4.4 Построение картограмм линий тока движения воздуха

4.5 Построение картограмм распространения газового облака при заданном направлении ветра

4.6 Исследование возможностей помехоустойчивой интерполяции для построения картограмм линий тока движения воздуха по жилым районам

4.7 Исследование точности метода контроля прогнозирования распространения промышленных газовых выбросов в городах

4.7.1 Неточность метеоданных

4.7.2 Влияние точности задания коэффициента диффузии на определение геометрических размеров облака

4.7.3 Влияние рельефа местности

Выводы

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

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

Данная проблема актуальна и для нашего региона. Тульская область отнесена к 12 наиболее экологически неблагополучным регионам России. В результате усиленной хозяйственной деятельности наиболее критическая экологическая ситуация сложилась в Новомосковском промышленном районе, который включает Новомосковский, Узловский и Кимовский районы. Район выделяется высоким сосредоточением' промышленности и занимает третье место в стране по концентрации химических производств на единицу площади. К химии следует добавить угольную промышленность и интенсивное сельское хозяйство с избыточным применением химических удобрений. Все это является мощным ударом по экологической обстановке. Основными источниками загрязнения атмосферы города Новомосковска являются промышленные предприятия такие, как ГРЭС, HAK "Азот", АО "Гипс Кнауф", ООО "Проктер энд Гэмбл" и крупные котельные. Они выбрасывают в воздух кроме пыли вредные вещества, которые взаимодействуя между собой создают неблагоприятный экологический фон.

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

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

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

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

Для достижения поставленной цели необходимо решить следующие задачи:

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

- разработка методики составления карты прогноза экологического состояния воздушной среды города;

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

- разработка математической модели распространения облака газа в приземном слое атмосферы;

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

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

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

- разработанный алгоритм и программа прогнозирования распространения облака газа в виде карты экологической обстановки воздушной среды города;

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

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

- предложенная методика расчёта распространения облака газа по территории города с учётом ландшафта местности;

- предложенная методика составления картограмм линий тока, движения воздуха, как на местности, так и в жилых районах по экспериментальным данным.

Результаты теоретических и экспериментальных исследований позволили:

- решить задачу оценки экологического состояния воздуха города;

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

- предложить структуру карты прогноза экологической обстановки воздушной среды города;

- предложить методику составления карты прогноза экологической обстановки воздушной среды жилого района города;

- разработать программу «Прогноз» для прогнозирования распространения облака газа по ландшафту местности.

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

Результаты диссертационной работы были использованы при подготовке территориальным управлением роспотребнадзора Новомосковского района, г. Донского, г. Кимовска и Узловского района ежегодного муниципального доклада о санитарной эпидемиологической обстановке в муниципальном образовании г. Новомосковска.

По теме диссертационной работы опубликовано девять печатных работ, в том числе две в журналах, рекомендованных ВАК.

Результаты работы докладывались и обсуждались на конференциях: Города России: проблемы строительства, инженерного обеспечения, благоустройства и экологии. (Пенза, 2009 г.), Окружающая среда и здоровье. (Пенза, 2009 г.), РХТУ им. Д.И.Менделеева, Новомосковский институт, 2008 -2009 гг.

Диссертационная работа состоит из введения, четырёх глав, заключения, списка использованных источников и пяти приложений. Работа изложена на 145 страницах машинописного текста, содержит 7 таблиц и 40 рисунков. Список использованных источников включает 130 наименований

Выводы

1. Приведены основные положения помехоустойчивой интерполяции.

2. Разработана программа «Прогноз», предназначенная для расчёта экологического портрета воздушной среды города и его жилых районов.

3. Построена математическая модели рельефа местности города Новомосковска.

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

5. Приведены картограммы распространения газового облака при заданном направлении ветра.

6. Исследованы возможностей помехоустойчивой интерполяции для построения картограмм линий тока движения воздуха по жилым районам.

7. Исследована точность метода контроля и прогнозирования распространения газовых выбросов в городах и в результате сопоставительного анализа метода было выявлено, что основными источниками погрешности прогнозирования являются метеоданные (скорость и направление движения ветра).

ЗАКЛЮЧЕНИЕ

Можно сделать следующие выводы по работе.

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

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

3. Предложена структура карты прогноза'экологической обстановки воздушной среды города.

4. Предложена методика составления карты прогноза экологической обстановки воздушной среды жилого района города.

5. Составлены карты прогноза экологической обстановки воздушной среды города Новомосковска и его отдельных характерных жилых районов.

6. Разработана программа «Прогноз» для прогнозирования распространения облака газа по ландшафту местности.

7. Результаты диссертационной работы использованы территориальным управлением роспотребнадзора Новомосковского района, г. Донского, г. Ки-мовска и Узловского района при подготовке ежегодного муниципального доклада о санитарной эпидемиологической обстановке в муниципальном образовании г. Новомосковска (приложение Г).

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

Материалы диссертационной работы используются в учебных курсах и научно-исследовательской практике Новомосковского института РХТУ им. Д.И. Менделеева (приложение Д).

1. Руководство по организации контроля состояния природной среды в районе расположения АЭС/ Под ред. К.П. Махонько.-Л.: Гидрометеоиздат, 1990, 264 с.

2. Techniques and decision making in the assessment of off-site consequences of an accident in a nuclear facility / Safety series, N.86, International Atomic Energy Agency. Vienne. 1987. -185 p.

3. Методика прогнозирования масштабов заражения сильнодействующими ядовитыми веществами при авариях (разрушениях) на химически опасных объектах и транспорте. Руководящий документ РД 52.04.253-90.-Л.: Гидрометеоиздат, 1991.-23 с.

4. Учет дисперсионных параметров атмосферы при выборе площадок для атомных электростанций. Руководство по безопасности АЭС. Международное агентство по атомной энергии. Вена, 1980. -106 с.

5. Методика расчета концентраций в атмосферном воздухе вредных веществ, содержащихся в выбросах предприятий. ОНД-86,- Л.: Гидрометеоиздат, 1987. 93 с.

6. Вызова Н.Л., Гаргер Е.К., Иванов В.Н. Экспериментальные исследования атмосферной диффузии и расчет pacпpocтpaнeнияí примеси,-Л.: Гидрометеоиздат, -1991. 54 с.

7. Монин A.C., Яглом A.M. Статистическая гидромеханика: Механика турбулентности,- М.: Наука, 1965. 720 с.

8. Пененко В.В., Алоян А.Е. Модели и методы для задач охраны окружающей среды. Новосибирск: Наука. Сиб.отд-ние, 1985.- 256 с.

9. Пащенко С.Э., Сабельфельд К.К. " Атмосферный и техногенный аэрозоль (кинетические, электронно-зондовые и численные методы исследования). Новосибирск. Ч. 1. 1992 . - 190 с.

10. Бородулин А.И., Майстренко Г.М., Чалдин Б.М. Статистическое описание распространения аэрозолей в атмосфере: метод и приложения.-Новосибирск: Изд-во Новосибирского ун-та. 1992. 123 с.

11. Берлянд М.Е. Прогноз и регулирование загрязнения атмосферы. Л.: Гидрометеоиздат, 1985.- 272 с.

12. Динамическая метеорология. Л.: Гидрометеоиздат, 1967. - 607 с.

13. Атмосферная турбулентность и моделирование распространения примесей. /Под ред. Ф.Т.М. Ньюстадта и Х.Ван Дона. Л.: Гидрометеоиздат. 1985.-351 с.

14. Методика прогнозирования масштабов заражения сильнодействующими ядовитыми веществами при авариях (разрушениях) на химически опасных объектах и транспорте. Руководящий документ РД 52.04.253-90. Л.: Гидрометеоиздат. 1991. - 23 с.

15. Руководство по организации контроля состояния природной среды в районе расположения» АЭС/ Под ред. К.П. Махонько.-Л.: Гидрометеоиздат. 1990. 264 с.

16. Turner D.B. Addendum to TUPOS Incorporator of a Hesitant Plume Algorithm. 1986.

17. Методика расчета концентраций в атмосферном воздухе вредных веществ,содержащихся в выбросах промышленных предприятий (ОНД-86). -Л.: Гидрометиоиздат, 1987.

18. ЕРА-600/8-86/0.27. U.S. Environmental Protection Agency, Research Triange Park, NC (available only from* NTIS, Accession Number PB86-241 031/AS).

19. Briggs G.A. Analytical parameterization of diffusion: the convective boundary layer//J. Clim. Appl. Met. 1985.V. 24/Pp. 1167-1186.

20. Turner D.B., Bender L.W., Paumier J.O., Boone P.F. Evaluation of the TUPOS air quality dispersion model using data from EPRI KINCAID field study //Atmos.Env. 1991. V. 25A. N.10. Pp. 2187-2201.

21. Venkatram A. Dispersion from an elevated source in a convective boundary layer//Atmos.Env. 1980. V. 14.N.l.Pp. 1-10.

22. Руководство по организации контроля состояния природной среды в районе расположения /АЭС. Под ред. К.П. Махонько. Л.: Гидрометеоиздат, 1990.- 264 с.

23. Вызова H.J1., Гаргер Е.К., Иванов В.Н. Экспериментальные исследования атмосферной диффузии и расчет распространения, примеси. -Л.: Гидрометеоиздат, 1991. 120 с.

24. Simpson I.R., Clarkson T.S. Dry plume: a computer model for predicting the behaviour of plumes in the atmosphere. Scientific report 19. New Zealand Meteorological Service, Wellington, 1986. -79 p.

25. Перри С.Г. Модель диффузии EPA для сложного рельефа: структура и характеристики. В сб.: Международная конференция ВМО по моделированию загрязнения атмосферы и его применениям. Л.: Гидрометеоиздат,1986. - 120 с.

26. Techniques and decision making in the assessment of off-site consequences of an accident in a. nuclear facility /Safety series, N.86, International Atomic Energy Agency. -Vienne, 1987. 185 p.

27. Пененко B.B., Короткое М.Г. Численная модель для исследования изменений климата и качества атмосферы мезорегионального масштаба // Математические проблемы экологии. Новосибирск: Изд-во ИМ СО РАН, 1994.- 160 с.

28. Yoshida A. Two-dimensional numerical simulation of thermal structure of urban polluted, atmosphere (effects of aerosol characteristics) //Atmos. Env., 1991. V. 25B.N. l.Pp. 17-23.

29. Годунов C.K., Рябенький B.C. Разностные схемы. M.:, Наука, 1973. -400 с.

30. Bianconi R., Tamponi М. A mathematical model of diffusion from a steady source of short duration in a finite mixing layer //Atmos.Env., 1993. V. 27A. N. 5, Pp. 781-792.

31. Методические указания по проведению анализа риска опасных производственных объектов: РД 03-418-01. введ.01.10.2001. - М., 2001. -25 с.

32. Chrysikopoulos C.V., Hildmann L.M., Roberts P.V. A three-dimensional steady-state atmospheric dispersion-deposition model for emission from a ground-level area source //Atmos. Env., 1992. V. 26A. N.5. Pp. 747-757.

33. Динамическая метеорология. Теоретическая метеорология / Под ред. Д. Л.Лайхтмана. Л.: Гидрометеоиздат. 1976. - 200 с.

34. Kitabayashi К. Wind tunnel simulation of airflow and pollutant diffusion offer complex terrain//Atm.Env. V. 25A. 1991. N 7. Pp. 1155-1161.

35. Защита атмосферы от промышленных загрязнений Справочник. Изд.: В 2-х ч. 4.2 Пер с английского. /Под редакцией Калверта С., Инглунда Г.М. М.: Металлургия, 1988. 712 с.

36. Государственный комитет РФ по охране окружающей среды. Количественный химический анализ атмосферного воздуха и выбросов в атмосферу. ПНДФ 13.1.2.-97. М.: 1997.

37. Контроль химических и биохимических показателей окружающей среды. Сб. под ред. Исаева В.М. П.: 1998, - 896 с.

38. Аварии и катастрофы. Предупреждение и ликвидация последствий. Учебное пособие. Книга 3. /Под редакцией.: В.А. Котляревского и A.B. Забегаева, М.; Изд-во АСВ, 1998 416 с.

39. Берлянд М. Е. Современные проблемы атмосферной диффузии и загрязнения атмосферы. -Л.: Гидрометеоиздат, 1975. —448 с.

40. Бурков В.Н. Экологическая безопасность/В.Н.Бурков, A.B. Щепкин. -2003

41. К.В.Глаголев, А.К.Морозов Физическая термодинамика, электронное учебное пособие МГТУ им. Н.Э. Баумана http://fh.bmstu.ru/phys/bib/physbook/tom2/about.htm

42. Калыгин В.Г. Промышленная экология. М. : Издательский центр «Академкнига», 2007.- 432 с.

43. Беляев Ю.И. Помехоусточивые методы и средства контроля параметров полей в системах технологического и экологического мониторинга//диссертация док.т.н 05.11.13. М: 2001, -405 с.

44. Латышенко К.П., Попов A.A. Информационно-измерительные системы для экологического мониторинга М.: МГУИЭ,2010. 230 С.

45. Беляев Ю.И. Синтез линейных схем оценивания скалярного поля методом стохастической интерполяции. Киев Автоматика.-1987,- №4. -43 с.

46. Временная типовая методика определения экономической эффективности осуществления природоохранных мероприятий и оценки экономического ущерба, причиняемого народному хозяйству загрязнением окружающей Среды. М.: Стройиз- дат, 1963.

47. Гавриленко В.В. Экологическая минералогия и геохимия месторождений полезных ископаемых. СПб., 1993. 150 с.

48. Жукова Г . С ., Митрохин^ С И Дарсалия В . Ш Дифференциальные уравнения / РХТУ им . Д . И . Менделеева . М ., 1999. 366с .

49. Закон РФ «Об охране окружающей природной среды»,1993

50. Кандель А.,Байатг У.Дж. Нечеткие множества, нечеткая алгебра, нечеткая статистика// ТИИЭР. Т.66, № 12 1978.

51. Кафаров В.В. Принципы создания безотходных химических производств. М.: Химия, 1982.

52. Ковалевский B.C. Влияние изменений гидрогеологических условий на окружающую среду. М., 1994. 138 с.

53. Мильков Ф.Н. Физическая география: современное состояние, закономерности, проблемы. В.: ВГУ. 1981.

54. Беляев Ю.И., Кораблев И.В., Вент Д.П. Помехоустойчивый контроль параметров полей в системах экологического мониторинга // Приборы. 2003. №1(31). С.35-38

55. Беляев Ю.И., Гербер Ю.В. Система мониторинга распространения газообразных веществ- в атмосфере города /«Датчика и системы», №10(125), Москва, 2009.- С. 13 16.

56. Беляев Ю.И., Вент Д.П., Гербер Ю.В., Вепренцева О.Н., Латышенко К.П. Экологический мониторинг атмосферы города/ «Химическое и нефтегазовое машиностроение», №8, Москва, 2009 .-С. 32 34.

57. Овчаренко Т.Б., Чернышов Д.Н. Отечественные приборы экологического контроля. -М.: ЦНИИТЭнефтехим, 1995, №2-3, -60 с.

58. Берленд М.Е. Прогноз и регулирование загрязнения атмосферы. Л.: Гидрометеоиздат, 1989, - 272 с.

59. Берленд М.Е. Современные проблемы атмосферной диффузии и загрязнения атмосферы. Л.: Гидрометоиздат, 1979, - 448 с.

60. Плотников Н.И., Карцев A.A., Рогинец Н.И. Научно-методические основы экологической гидрогеологии. М., 1992. 62 с.

61. Полоцкий Л . M Лапшенков Г . И . Автоматизация химических производств . M .: Химия , 1982. 296 с .

62. Современные проблемы инженерной геологии и гидрогеологии территорий городов и городских агломераций/ Тез. докл. III Всес. семинара. М., 1987. 408 с.

63. Эдельштейн Ю.Д., Вент ДП. Концептуальный подход к построению АКСМ «Новомосковск». Вестник Академии: Информатика, Экология, Экономика, T.I 4.1- М.:РАДСИ,1977.

64. Защита атмосферы от промышленных загрязнений Справочник. Изд.: В 2-х ч. 4.2 Пер с английского. /Под редакцией Калверта С., Инглунда Г.М. М.: Металлургия, 1988. 712 с.

65. Защита атмосферы от промышленных загрязнений Справочник. Изд.: В 2-х ч. 4.2 Пер- с английского. /Под редакцией Калверта С., Инглунда Г.М. М1: Металлургия, 1988. 712 с.

66. Маршалл В. Основные опасности химических производств: Пер. с англ.// Под ред. Б. Б. Чайванова, A'. Н. Черноплекова. -М:: Мир, 1989. -672 с.

67. Беляев Ю.И. Синтез линейных схем оценивания скалярного поля методом стахостичекой интерполяции- Киев: Автоматика,1987№4. 32-34 с.

68. Васильев Ф.П. Численные методы решения экстремальных задач,- М.: Наука.1988 522 с.

69. Горчиев А. А., Рафиев Р. М. О прогнозе загрязнения воздуха над Апшеронским полуостровом. // Метеорология и гидрология,-1980. №2 -21 -26 с.

70. ГОСТ 17.0.01-76 Система стандартов в области охраны природы и улучшения использования природных ресурсов

71. ГОСТ Р 8.589-2001 ГСИ Контроль загрязнения окружающей природной среды. Метрологическое обеспечение. Основные положения

72. ГОСТ Р 8.563-96 ГСИ Методика выполнения измерений

73. ГОСТ Р ИСО 14001-98 Система управления окружающей средой. Требования и руководство по применению

74. ГОСТ Р ИСО 14050- 99 Управление окружающей средой. Словарь

75. Пигурнов Б. В. Интерактивная машинная графика в автоматизированных системах контроля загрязнения воздуха. Энергетика и электрификация, 1982 - №4 - С. 43-47.

76. Практикум по численным методам прогноза погоды / Под ред. Л. С. Гандина. Л.: Гидрометеоиздат, 1978. - 212 с.

77. Таран Т. А., Пигурнов Е. В., Хомяков А. Т. Анализ алгоритмов оперативного прогноза загрязнения воздуха, базирующихся на статистических моделях // Сб.: Проблемы контроля и защиты атмосферы от загрязнения.- Киев: На-укова думка, 1983 Вып. 9,- С. 34 - 37.

78. Литвин В.А. Многокритериальная автоматизированная региональная система моделирования эффективных стратегий. М.: Гидрометеоиздат, 1988, -188 с.

79. Бокс Дж„ Дженкинс Г.Д. Анализ временных рядов. Прогноз и управление. Вып.1. Пер. с англ. М., «Мир», 1974,-406 с. 95'. Венцель Е.С. Теория вероятностей: М.: 1962, -160 с.

80. Донченко В.К. Экологическая экспертиза/ В.К. Донченко. М.: Издательский центр «Академия», 2010.-528 с.

81. Берлянд М.Е., Генихович E.JI., Грачёва И.Г. Об усовершенствовании методов расчета загрязнения атмосферы // Сб. тр. ГТО им. А.И. Воейкова, 1987.Вып №1- 4-11с.

82. Вожодаев М.В., Ивлева Т.П. Прогноз высоких уровней загрязнения воздуха, создаваемого выбросами автортрансрпорта// Информ. Бюлл. ГГо им. А.И. Воейкова,1987.Вып 511- 4-11с.99.

83. Берлянд М.Е., Генихович E.JL, Оникул Р.и. Моделирование загрязнения* атмосферы выбросами из низких и холодных источнтков// Метеоролия и гидрология. 1990- №5.- М.: Гидрометиоиздат, 1990. 5-17с.

84. Другов Ю.С. Экологическая аналитическая химия/ Ю.С. Другов. Спбю: Изд-во Анатолия, 2000. 432 с.

85. Загорский В.П. Информационно-измерительные системы коллективного пользования/ В.П. Загорский,- Минск: Наука и техника, 1987 104 с.

86. Капиев Р.Э. Измерительно-вычислительные комплексы/ Р.Э. Капиев. -JI. :Энергоатомиздат,1988.- 176 с.

87. Кашкин В.Б. Дистанционное зондирование Земли из космоса/ В.Б. Кашкин, А.И. Сухинин. М.: Логос, 2001 - 264 с.

88. Вент Д.П., Беляев Ю.И., .Гербер Ю.В., Вепренцева О.Н. Модели распространения загрязнений атмосферы города Информатика, Экология, Экономика. Вестник международной академии системных исследований (МАСИ). М.: 2007. С. 26 28.

89. Беляев Ю.И., Гербер Ю.В., Микаилов Р.И. Основные принципы моделирования атмосферы города как объекта экологического мониторинга. VI Международная научно-практическая конференция. Окружающая среда и здоровье. МК-24-9, Пенза, 2009. С. 12 - 15.

90. Васильков Ю.В. Компьютерные технологии вычислений в математическом планировании/ Ю.В. Васильков, H.H. Василькова. М.: Финансы и статистика, 2002.- 256 с.

91. МИ 1730-87 Погрешности косвенных измерений характеристик процессов. Методика расчётов

92. Клюев В.В. Экологическая диагностика: Энциклопедия/ В.В. Клюев-М.: Машиностроение,2000.- 496 с.

93. Тартаковский Д.Ф: Метрология, стандартизация и технические средства измерений/ Д.Ф. Тартаковский, A.C. Ястребов.- М.: Высшая школа, 2001.- 346 с.

94. Сергеев А.Г. Метрология/ А.Г. Сергеев, В.В. Крохин,- М.: Логос,2001.-376 с.

95. Беляев Ю.И., Гербер Ю.В., Моисеева И.Д. Бизина С.А. Оперативное прогнозирование состояния атмосферы города. Успехи в химии и химической'технологии: сб. науч. Тр. Том XXII, № 7(101). М. : РХТУ им. Д.И. Менделеева, 2009. - С. 122-123.

96. Кафаров В.В. Принципы создания безотходных химических производств. М.: Химия, 1982.

97. Кафаров В.В., Дорохов H.H., Елисеев П.И., Вербато Е.Г. Интерактивные задачи экспертных систем управления// Доклады АН СССР, Т.305,№ 5,1989.

98. Кафаров В.В., Дорохов И.Н., Марков Е.П. Системный анализ процессов химической технологии. Применение метода нечетких процессов химической технологии. Применение метода нечетких множеств. М.: Наука, 1983.

99. Перельман А.И. Геохимия. М., 1979. 423 с.

100. Современные проблемы инженерной геологии и гидрогеологии территорий городов и городских агломераций/ Тез. докл. III Всес. семинара. М., 1987. 408 с.

101. Emmernegger С. Le programme NADUF dans le cadre la surveillance qualitative descours d'eau Suisses. uGas-wasser-Abwasser", 1985.

102. Эдельштейн Ю.Д., Котельников A.A., Наумов В.Ю., Мягкова Г.И. Автоматизированные системы экологического мониторинга. Ч.Ш. Учебное пособие/Под редакцией д.т.н. проф. Д.П.Вента, РХТУ им. Д.И.Менделеева, Новомосковский институт, -2000.

103. Гавриленко В.В. Экологическая минералогия и геохимия месторождений полезных ископаемых. СПб., 1993. 150 с.

104. Munn R.E. Global. Environmental monitoring system. Action plane for phase 1. SCOPE, 1983.

105. Вредные химические вещества: Справ.изд./ Под ред В.А.Филова и др. Л., 1988 -Î990. Т. 1 4.

106. Методика расчета распространения аварийных выбросов основанная на модели рассеивания тяжелого газа //Безопасность труда промышленности 2004. №9? С. 38-42.

107. Сборник временных отраслевых методик для определения загрязняющих веществ в атмосферном воздухе в. районах размещения предприятий нефтеперерабатывающей 1985, с. 24.

108. Овчаренко Т.Б., Чубукова Н.М., Чернышов Д.Н. Приборы и средства экологического контроля. М.: ЦНИИТЭнефтехим, 1996, 1, -42 с.

109. Горстко А.Б. Познакомьтесь с математическим моделированием. М.: Знание, 1991.

110. Виженский В.А., Петрухин В.А. Мониторинг фонового загрязнения природных сред,- Л.: Гидрометеоиздат, 1990. Вып. 6.