автореферат диссертации по безопасности жизнедеятельности человека, 05.26.02, диссертация на тему:Разработка методов оценки экосистемных рисков в зонах воздействия выбросов на объектах газовой промышленности

□□3057800

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

Демидова Ольга Анатольевна

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

05 26 02 — Безопасность в чрезвычайных ситуациях в нефтяной и газовой

промышленности

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

Москва-2007

003057800

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

Демидова Ольга Анатольевна

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

05 26 02 - Безопасность в чрезвычайных ситуациях в нефтяной и газовой

промышленности

АВТОРЕФЕРАТ

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

Москва - 2007

Работа выполнена в Обществе с ограниченной ответственностью «Научно-исследовательский институт природных газов и газовых технологий -ВНИИГАЗ» и в Некоммерческом партнерстве «Центр по экологической оценке «Эколайн»

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

доктор технических наук Казак А С

Официальные оппоненты

доктор технических наук, профессор Бухгалтер Э Б

доктор технических наук, профессор Акимов В А

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

Московский филиал ОАО «Гипроспецгаз»

Защита диссертации состоится 22 мая 2007 г в 1330 ч на заседании диссертационного совета Д 511 001 02 при ООО «ВНИИГАЗ» по адресу 142717, Московская область, Ленинский район, пос Развилка

С диссертацией можно ознакомиться в библиотеке ООО «ВНИИГАЗ»

Автореферат разослан « апреля 2007 г

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

Курганова И Н

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы исследования. Строительство и эксплуатация производственных объектов ведут к преобразованию природных комплексов (экосистем) в зонах воздействия Минимизация атмотехногенных воздействий на природную среду приобретает особое значение для газовой отрасли в связи с ростом объемов добычи, транспортировки, переработки и потребления газа а, следовательно, и объемов эмиссии поллютантов Так, в 2005 г валовые выбросы загрязняющих веществ (ЗВ) в атмосферу на объектах газовой промышленности России достигли 2351 тыс т Свыше 10% от этого объема составляют оксиды азота и серы, которые признаны значимыми химическими факторами опасности для экосистем

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

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

Цель работы - разработка методов количественной оценки экосистемных рисков в зонах воздействия выбросов ЗВ от объектов газовой промышленности для обеспечения экологической безопасности их эксплуатации и предотвращения ЧС техногенного характера Основные задачи исследования

1 Анализ современных принципов и методов оценки атмотехногенных воздействий на экосистемы, включая методы оценки экосистемных рисков, применительно к объектам газовой промышленности

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

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

4 Анализ результатов оценки рисков и разработка рекомендаций по управлению рисками для экосистем территорий исследования

В качестве объектов демонстрационных исследований выступали наземные экосистемы

• зоны атмотехногенного воздействия объектов газовой промышленности (ГП) Боливарианской Республики Венесуэла (ГП Венесуэлы) (региональный уровень),

• территории размещения объектов Средне-Тиманского бокситового рудника (СТЕР) в Республике Коми (локальный уровень)

Территория СТЕР испытывает комплексное воздействие выбросов ЗВ от объектов СТЕР и от компрессорных станций газопровода «Ямал - Центр»

Предмет исследования — методы количественной оценки негативных воздействий на экосистемы, связанных с выбросами ЗВ на технологических объектах

Научная новизна. Разработаны новые методические подходы к количественной оценке экосистемных рисков, связанных с выбросами соединений серы, азота и некоторых тяжелых металлов (ТМ), развивающие положения методологии критических нагрузок (КН) поллютантов, направленные на предотвращение развития ЧС техногенного характера Предложено теоретическое обоснование и разработан алгоритм расчета показателей экосистемного риска на основе вероятностного моделирования величин КН и их превышений по методу Монте-Карло Обоснована схема проведения оценки экологических рисков в рамках работ по оценке воздействий действующих или проектируемых производственных объектов на окружающую среду Предложен метод определения допустимых уровней атмотехногенной нагрузки приоритетных ЗВ для зон воздействия выбросов ЗВ от технологических объектов

Впервые выполнен расчет величин КН соединений серы и азота и их превышений для наземных экосистем Южной Америки (в границах зоны атмотехногенного воздействия объектов ГП Венесуэлы)

Разработанные методические подходы были впервые применены для количественной оценки атмотехногенных воздействий на экосистемы в рамках демонстрационных исследований для выбранных ЗВ и сценариев воздействия Результаты оценки были использованы для выработки рекомендаций по управлению экосистемными рисками с целью предотвращения развития ЧС техногенного характера

Защищаемые положения

1 Концепция количественной оценки экосистемных рисков, связанных с выбросами производственных объектов, на основе величин КН поллютантов

2 Алгоритм расчета показателей экосистемного риска на основе моделирования величин превышений КН по методу Монте-Карло

3 Обоснование схемы проведения оценки экосистемных рисков, связанных с выбросами соединений серы, азота и ТМ, при оценке воздействий на окружающую среду действующих и/или проектируемых производственных объектов газовой промышленности

4 Метод определения безопасных уровней атмотехногенной нагрузки приоритетных ЗВ для экосистем, затронутых выбросами объектов газовой промышленности

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

нием в окружающую среду соединений серы и азота, ТМ Результаты оценки экосистемных рисков предназначены для выработки оптимальных стратегий по снижению выбросов производственных объектов до экологически безопасных уровней с целью предотвращения ЧС техногенного характера

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

Результаты диссертационного исследования использованы при разработке Генеральной схемы развития газовой промышленности Боливарианской Республики Венесуэла и при экологическом обосновании расширения СТЕР

Апробация и публикация результатов исследований. Материалы исследований по теме диссертации докладывались и обсуждались на заседаниях кафедры геохимии ландшафта и географии почв Географического факультета МГУ им Ломоносова, на 24-й и 26-й конференциях Международной ассоциации оценки воздействия в Ванкувере (2004 г), в Ставанге (2006 г), на Шестой субрегиональной встрече по деятельности в области эффект-ориентированных методов оценки под эгидой Конвенции о трансграничном загрязнении воздуха на большие расстояния ЕЭК ООН (Москва-Пущино, 2004 г), на 8-й Всероссийской научно-практической конференции «Геоинформатика в нефтегазовой и горной отраслях» (Москва, 2007 г )

По материалам диссертации опубликовано 10 работ, в том числе 3 работы в журналах, включенных в «Перечень » ВАК Минобрнауки РФ

Структура и объем диссертации. Диссертация состоит из введения, четырех глав, выводов и списка литературы из 195 наименований Работа изложена на 172 страницах машинописного текста и содержит 38 рисунков, 15 таблиц и 4 приложения

Автор выражает глубокую благодарность д т н АС Казаку, д б н , проф ВН Башкину, к гн ИВ Припутиной, доценту Центрально-Европейского Университета О М Черпу за их ценные идеи и методическую поддержку Большую помощь в сборе и обработке исходных данных для исследования оказали начальник Департамента по экологии, охране труда и промышленной безопасности ЗАО «Коми Алюминий» Е А Юркин и главный специалист Департамента М В Котова, специалисты ЗАО «Группа компаний «ШАНЭКО», а также В М Марунин, Н В Павлов и А В Танканаг Автор выражает признательность коллективу НП «Центр по экологической оценке «Эколайн» и лично директору Центра к х н М В Хотулевой за содействие в подготовке и проведении диссертационного исследования

Оценка экосистемных рисков для района СТБР осуществлялась на основе данных, предоставленных ЗАО «Коми Алюминий» и ОАО «Боксит Тимана» Правом собственности на результаты данной части практического компонента диссертационного исследования обладает ЗАО «Коми Алюминий»

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

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

Глава 1 посвящена анализу принципов и методов оценки атмотехногенных воздействий производственных объектов на экосистемы применительно к объектам газовой промышленности

В работах российских и зарубежных ученых (В С Безеля, А А Быкова, Е Л Воробейчика, К Н Дьяконова, Т В Звонковой, 3 Тгеууеек, IС БтгсЬек, N3^0) указывается на приоритет экосистемного подхода к выбору и характеристике реципиентов техногенного воздействия Совершенствование методологии оценки воздействий на экосистемы идет, на наш взгляд, в следующих направлениях

• развитие количественных методов прогнозирования воздействий, включая математическое моделирование и геоинформационные системы (ГИС),

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

• уход от детерминистского подхода к прогнозированию воздействий и развитие методов оценки экосистемных рисков

Рассмотрена методология анализа экологических рисков (АЭР) и ее роль при подготовке экологического обоснования строительства и эксплуатации промышленных объектов, в том числе объектов газовой отрасли

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

Концептуальные основы количественной оценки экосистемных рисков, связанных с выбросами производственных объектов, на основе интегрального показателя — устойчивости экосистем к аэрогенному поступлению поллютантов - рассмотрены в работах А С Казака и В Н Башкина При этом роль локальных (экосистемных) нормативов воздействия выполняют величины КН поллютантов Величина критической нагрузки представляет собой такой объем выпадений поллютанта, который не приведет к необратимому изменению структуры и функций экосистем-реципиентов в течение длительного периода времени (50-100 лет)

Рассмотрено состояние современной методологии КН и опыт ее применения в различном масштабе, в том числе при оценке воздействия на экосистемы, расположенные в зонах воздействия выбросов магистральных газопроводов (региональный масштаб) Показано, что для предупреждения отложенных ЧС техногенной природы необходима дальнейшая разработка методов характеристики экосистемных рисков и определение общего порядка проведения оценки рисков в рамках данной методологии

В главе 2 изложены методические подходы к количественной оценке экосистемных рисков в зонах воздействия объектов газовой промышленности, на основе величин КН поллютантов

Следуя концепции экологического риска И. V КоПига, экосистемный риск определен как комплексный показатель, характеризующий вероятность развития негативных изменений в состоянии экосистем, и величину таких изменений

В данной работе количественную оценку экосистемных рисков предложено проводить на основе расчета и пространственного анализа величин превышений критических нагрузок приоритетных ЗВ Превышение критической нагрузки поллютанта X (Ех(Х)) рассчитывается согласно формуле (1)

Ех(Х) = Х1кр-СЦХ), (1)

где Ех(Х) — превышение критической нагрузки поллютанта X; данная

величина может принимать положительные, отрицательные значения или равняться нулю, Хл,р — аэрогенное поступление поллютанта X в экосистемы (выпадения), СЬ(Х) — величина КН поллютанта X

Таким образом, величина КН поллютанта играет роль референтной дозы — общепринятого норматива воздействия в исследованиях по оценке экологических рисков

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

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

Функция экосистемного риска (Лх) является функцией распределения Она описывает вероятностное распределение случайной величины рх, которая

характеризует площадь ареала превышений КН поллютанта относительно площади группы рецепторных участков (в %) (М(Ех(Х)>0))

Ях=Г(() = Р(рх<1) (2)

Функция экосистемного риска позволяет рассчитать

• вероятность того, что превышения КН будут отмечены на территории, меньшей или равной заданному значению М(Ех(Х)>0),

• вероятность того, что величина М(Ех(Х)>0) будет находиться в заданном интервале значений

Для вычисления случайной величины рх необходимо определить вероятности превышения КН поллютанта для каждого рецепторного участка (р,1) и площади этих рецепторных участков (5,)

Для расчета величин р,' предложено проводить моделирование превышений КН по методу Монте-Карло Традиционно расчет превышений КН осуществляется на основе оптимальных значений входных параметров В рамках данного подхода расчет величин Ех(Х) осуществляется многократно, для каждого модельного эксперимента по случайному закону формируется набор значений входных параметров на основе имеющихся массивов таких значений Вариабельными (случайно изменяемыми) могут быть выбраны некоторые входные параметры для расчета величин превышений КН

Согласно центральной предельной теореме, при большом количестве рецепторных участков небольшого размера функция распределения случайной величины рх как суммы независимых малых случайных величин хорошо аппроксимируется функцией нормального распределения Для проверки нормальности функции экосистемного риска предложено подготовить выборку значений М(Ех(Х)>0) и проверить гипотезу о нормальном распределении данного параметра с помощью критериев согласия (например, критерия Шапиро-Уилка) Если распределение не подчиняется нормальному закону, то для расчетов показателей экосистемного риска необходимо вычислять эмпирическую функцию распределения рх.

При проведении оценки экосистемных рисков на основе КН поллютантов предложено следовать формальной процедуре оценки экологических рисков, включающей четыре этапа

На этапе идентификации опасности следует установить источники и состав выбросов проектируемого предприятия, определить возможные сценарии воздействия Необходимо также охарактеризовать реципиенты воздействия (экосистемы в границах зоны атмотехногенного воздействия действующего или проектируемого объекта) и провести их ранжирование (выделить особо ценные и/или уязвимые экосистемы) и провести подразделение территории исследования на рецепторные участки На основе имеющейся информации о факторах опасности и реципиентах производится качественная характеристика экосистемных эффектов и определяется список

поллютантов, для которых целесообразно проведение детальной оценки риска (приоритетных ЗВ)

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

Характеристику экосистемных рисков предложено проводить в два этапа На первом этапе необходимо осуществить детерминистический расчет превышений КН на основе оптимальных значений входных параметров В случае выявления рецепторных участков с Ех(Х)>0, на втором этапе целесообразно провести имитационное моделирование величин превышений КН по методу Монте-Карло и расчет показателей экосистемного риска Рекомендуемое число модельных экспериментов - не менее 1000

Приемлемость атмотехногенного воздействия зависит от ценности затронутых экосистем Для групп рецепторных участков, которые имеют статус особо ценных и/или уязвимых, величины КН не должны быть превышены на 100% их территории Для территории санитарно-защитной зоны производственного объекта допустимым можно считать наличие превышений КН на 50% территории В остальных случаях предложено следовать принципу «95%-ой защищенности» экосистем

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

Разработанные методические подходы позволяют выполнять количественную оценку не только величины прогнозируемых изменений в состоянии экосистем, но и вероятности их наступления, учитывают природную вариабельность входных параметров Результаты оценки рисков для экосистем предложено использовать для расчета безопасных для экосистем уровней эмиссии приоритетных ЗВ (управления экосистемными рисками), а также для ранжирования отдельных проектных альтернатив (технологических, планировочных)

Сферу применения предлагаемых подходов ограничивает современное состояние методологии КН Спектр поллютантов, для которых разработаны алгоритмы расчета КН, ограничен в него входят соединения серы, азота, некоторые тяжелые металлы (РЬ, Сс1, Си, Ъа, Н§), биогеохимические модели ориентированы на наземные (в первую очередь, лесные) и пресноводные экосистемы

Врезка 1. Алгоритмы и входные параметры для расчета КН поллютантов для наземных экосистем Максимальная КН серы (СЬтах@)) С1т11х($) = ВС^ер + ВСК - С1Лер - Всир, - АЫСщсм), КН питательного азота (СЬпш(Ы)) СЬШ(Ы) =АГ,„ +АГ„Р/+ Л^асс/( 1 ~/<к)- ЛШ[е(СП1), Максимальная КН азота (СЬтах(Ы)) С1тах(Ы) = Ыт +Ныр[ + СЬтах(8)!{ 1 -/¡¡е), КН тяжелого металла (ТМ) (СЦМ)) СЬ (М) = Мир/ + Мщ^)

Условн обозн Наименование входного параметра Способ расчета/используемые значения

BCdep 2+ 2+ Поступление катионов Са , , К , N3 из атмосферы, экв /га/год Входными параметрами для расчета являются количество атмосферных осадков и концентрации химических элементов в атмосферных выпадениях

Cltfep Поступление анионов хлора из атмосферы, экв /га/год

BCW Внутрипочвенное выветривание 2+ 2+ + + катионов Са , , К , N3 , экв /га/год Зависит от класса текстуры почв и среднегодовой температуры

Bcupt Долговременное депонирование 2+ 2+ + Са , , К в древесной биомассе, экв /га/год Рассчитываются на основе показателя ежегодного прироста древесины (Ywood) и концентрации химических элементов в древесине ([Вс]др, [N]dp, [М]др)

Nup, Долговременное депонирование азота в древесной биомассе, экв /га/год

Mupt Долговременное депонирование ТМ в древесной биомассе, г/га/год

ANCle(cr,t) Критическое вымывание щелочности, экв /га/год Определяется величиной слоя инфильтрации осадков (Q[e) и критическим значением рН почв для разных видов растительности

Mie(cnt) Критическое вымывание ТМ из рассматриваемого почвенного слоя, г/га/год Рассчитывается на основе геохимических характеристик рН, содержания органического вещества в почве ([OMJs), содержания илистой фракции в почве ([clay]), концентрации растворенных органических веществ в почвенном растворе ([DOCJss)

N,m Долговременная иммобилизация азота в почве, экв /га/год Определяется режимом увлажнения и типом растительности

Nlefacc) Допустимое вымывание азота, экв /га/год Определяется величиной слоя инфильтрации осадков (Qie) и критической концентрацией азота в почвенном растворе ([NJacc) (зависит от вида растительности)

fde Коэффициент денитрификации азота в почве Зависит от режима увлажнения почв

В главе 3 представлены результаты апробации предлагаемых подходов к оценке экосистемных рисков на региональном уровне в рамках работ по экологическому обоснованию Генеральной схемы развития ГП Венесуэлы В работе использованы данные только открытых литературных источников

Идентификация опасности

Состав природного газа, добываемого на различных месторождениях Венесуэлы, характеризуется наличием таких основных примесей как N2, Н28, СгНб Технология добычи и транспортировки газа включает его сжигание на компрессорных станциях (КС), что приводит к выбросам оксидов азота (Ж)х), оксида углерода (СО) и диоксида серы (БОг) Кроме того, при эксплуатации объектов ГП в атмосферу поступает природный газ (его главным компонентом является метан) и меркаптаны

На сегодняшний день основными источниками выбросов ЗВ в ГП Венесуэлы являются объекты добычи и транспортировки газа, включая КС и дожимные компрессорные станции (ДКС) на промыслах В связи с планируемой модернизацией газотранспортной системы ее вклад в загрязнение атмосферы существенно возрастет, поэтому границы зоны атмотехногенного воздействия объектов ГП Венесуэлы определялись относительно трасс магистральных газопроводов Согласно экспертным оценкам, вследствие перекрывания зон влияния выбросов отдельных КС, воздействие их выбросов может прослеживаться на расстоянии до 50 км от трассы газопровода

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

БОг и >ГОХ были выбраны в качестве приоритетных ЗВ для детальной оценки атмотехногенных рисков для экосистем С ними связаны следующие негативные экосистемные эффекты повышение кислотности почв и изменение видового состава биоценозов вследствие эвтрофирования

Расчеты экосистемного риска проводились для сценария воздействия, соответствующего современному уровню атмотехногенной нагрузки В составе зоны атмотехногенного воздействия объектов ГП Венесуэлы выделены две группы рецепторных участков 1) зона воздействия «Запад» с высокими уровнями атмотехногенной нагрузки в районе оз Маракайбо в связи с деятельностью предприятий химической промышленности, 2) зона воздействия «Центр-Восток», для которой повышенная атмотехногенная нагрузка характерна только для столичного региона и отдельных районов штатов Монагас и Анако

Оценка экспозиции

Оценка экспозиции включала определение величин аэрогенного поступления (выпадений) серы {Sjep) и азота {А1',^,)', а также основных катионов (BCjnp) и анионов хлора (С/,^), по данным исследований химического состава атмосферных осадков. Расчетные величины выпадений представлены в табл. 1.

Таблица 1

Величины выпадений основных катионов (ВС&р), анионов хлора {С1^р), серы (Sdep), азота (Njt,p) для юны атмотехногенного воздействия объектов

( И Венесуэлы

Группа рецепторных участков ВСаер S/lcp

эк в,/г а/год

Зона воздействия «Запад» 370 (45-1480) п =12 138 (50-330) п =72 219 (169-800) а =8 443 (150-957) п =1

Зона воздействия «Центр-Восток» 71 (31-244) п=8 120 (69-238) п =8

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

Оценка '¡ффекитч

Оценка экосистем и ых рисков, связанных с выбросами этих пол л юта НТО в, потребовала расчета величин максимальной критической нагрузки серы (а„Н1Х(5)), которая определяет устойчивость к кислотности, критической нагрузки питательного азота (С/.;„,,(ЗД) и максимальной критической нагрузки азота (а„,1а(М)) (врезка 1),

Расчетные Значения величин максююяыюй КН серы для наземных экосистем зоны атмотехногенного воздействия объектов ГП Венесуэлы варьируют от 100 до 3000 экв./та/год. Для большей части территории исследования допустимый уровень кислотных выпадений составляет 2000 экв./кг в год и выше. Лесные экосистемы, включая заболоченные леса на востоке территории исследования, характеризуются достаточно высокой устойчивостью к подкислению (а„Н1Х(3) от 1000 до 2000 экв./кг в год). Самые

1 Величины аэрогенного поступления ТМ принято рассчитывать в г/га/год, серы и азота - в эквивалентах (окв.)/га/год. Величина нагрузки в I экв./га'год соответствует 0,016 кЙ>/га/год и 0,014 кгМ/га/тед,

низкие значения ЕЛтах(5) (менее 250 экз./га/год) установлены для участков на побережье Карибского моря, занятых как лесами, так и опустыненными степями и пустынями, и горных лесов на маломощных почвах. Доля таких экосистем не превышает 2% от площади территории исследования.

Величина КН питательного азота в наибольшей степени зависит от относительно долговременного депонирования данного элемента в древесной биомассе растительности Ареалы с низкими значениями СЬ,Ш,(М) (ниже

250 экв./га/год) составляют участки с травянистой растительностью, большей устойчивостью отличаются кустарниковые экосистемы; для лесных экосистем характерны высокие значения СЬт„(Ъ1) (1000 экв./га/год и выше). Пространственное распределение величин СЬ,Ш,(И) иллюстрирует рис. 1.

УСЛОВНЫЕ ОБОЗНАЧЕНИЯ

CLnul(N), эка,/га/год газовые промыслы

<250 —— система магистральных газопроводов

250 - 500 А газокомпрессорная станция

500 - 1000 - границы зоны атмотехногеннсго воздействия

^И 1000- 1500 ГГ7] группы рецепторных участков

^■>1500 i^m рецепторные участки, для которых не определен

^^ тип растительности

Рис. 1. Величины КН питательного азота (CLnul(N)) для наземных экосистем шны атмотехмогенного воздействия объектов ГП Венесуэлы

На 70% площади исследуемой территории величины максимальной КН азота составляют 2000 экв./га/год и выше, что сопоставимо с современным уровнем нагрузок в Европе. Наиболее высокие значения CLmiU(N) (3500 экв./га/год и выше) определены для тропических н мангровых лесов к районе дельты р. Ориноко и побережья Карибского моря. Минимальные значения данного показателя характерны для района опустыненных степей и пустынь и урбанизированных территорий, где рассчитаны низкие значения CLmax(S). Суммарная площадь этих экосистем составляет менее 1% площади зоны территории исследования.

Характеристика риска

В рамках характеристики рисков для наземных экосистем зоны атмотехногенного воздействия объектов ГП Венесуэлы были проведены: 1) расчет превышений КН соединений серы (Ех(5)% превышений КН соединений азота (Ех(Щ) и совместных превышений КН соединений серы и азота (Ех(3+№)) с использованием оптимальных значений входных параметров (детерминистический расчет) и 2) моделирование этих величин по методу Монте-Карло (вероятностный расчет).

В качестве вариабельных входных параметров при вероятностном расчете превышений КН были выбраны показатели содержания основных катионов и азота а древесной биомассе ([Вс],Ьр, /Л7ф) и величины выпадений основных катионов, анионов хлора, а также соединений серы и азота (ВС^р, С^р, Л^ер) ввиду отсутствия пространственно распределенных данных по выпадениям. Для всей территории исследования был принят единый критерий приемлемости воздействия - 95%-ая защищенность экосистем.

Результаты детерминистического расчета превышений КН показали, что величины современной нагрузки соединений серы и азота для отдельных рецеи горных участков превышают расчета О-допустимые уровни. Положительные значения Ех($') отмечены только для зоны воздействия «Запад», превышения КН по азоту установлены для экосистем в пределах обеих групп рецепторных участков. Пространственное распределение величин Ех(5+И) для территории исследования показано на рис. 2,

Марзкзиоо

УСЛОВНЫЕ ОБОЗНАЧЕНИЯ

Ех(5+М). зкв./га/год ^ газовые промыслы

О — система магистральных газопроводов

<100 Д газокомпрессорная станций

вц то . 200 - границы зоны атмотехногенного воздействия

■■ 200- 300 ГО! группы рецепторных участков

рецегтторные участки, для которых не определен тип растительности

Рис. 2. Величины совместных превышений КН соединений серы и азота (Ех(8+Ы)) для наземных экосистем зоны атмотехногениого воздействия объектов ГЦ Венесуэлы

Рецепторные участки с наибольшими значениями Ех('£+Л£)>0 относятся к бассейну оз Маракайбо (зона воздействия «Запад») и приурочены к природным комплексам опустыненных степей и пустынь и территориям с городской застройкой В зоне воздействия «Центр-Восток» выделено небольшое число рецепторных участков с минимальными значениями совместных превышений КН серы и азота (менее 100 экв /га/год) (рис 2)

Для рассматриваемой территории исследования площадь ареала совместных превышений КН соединений серы и азота (М(Ех(Б+Ы)>0)) составляет менее 5%, поэтому уровень современной нагрузки этих ЗВ может считаться приемлемым (табл 1) Однако для зоны воздействия «Запад» величины М(Ех(Ы)>0) и М(Ех(Б+И)>0) составили около 7,0%, превысив установленный критерий приемлемости атмотехногенного воздействия

Площади ареалов превышений КН, полученные в результате вероятностного расчета, в два и более раз выше значений, полученных на основе оптимальных значений входных параметров (табл 2) Для зоны воздействия «Центр-Восток», где детерминистический расчет показал отсутствие превышений КН соединений серы, также выявлены рецепторные участки с Ех(Б)>0

Таблица 2

Площадь ареалов превышений критических нагрузок (М(Ех(8)>0), М(Ех(Ы)>0), М(Ех(Б+Ы)>0), %) для наземных экосистем зоны атмотехногенного воздействия объектов ГП Венесуэлы

Детерминистический расчет Вероятностный расчет (р,'>0)

Факторы опасности Группа рецепт участков Б N 5+ЛГ Б N Я+ЛГ

Зона воздействия «Запад» 2,5 6,9 7,0 29,1 23,6 64,4

Зона воздействия «Центр-Восток» 0,0 4,0 4,0 5,2 9,8 14,7

Зона атмотехногенного воздействия объектов ГП Венесуэлы <0,1 4,1 4,1 5,5 10,0 15,3

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

нормального распределения На основе функций экосистемного риска были выведены 6 классов значений М(Ех(Х)>0) и проведен расчет вероятностей превышения КН для отдельных интервалов значений М(Ех(Х)>0) ((Р(Ех(Б)>0), Р(Ех(Ы)>0) и Р(Ех(Б+Ы)>0)) В табл 3 представлены значения параметра Р(Ех(8+Ы)>0) Только для зоны воздействия «Запад» результаты детерминистического и вероятностного расчетов существенно расходятся

С вероятностью 92,9% площадь ареала совместных превышений КН соединений серы и азота варьирует в пределах 15-20% от площади данной группы рецепторных участков, в то время как детерминистическое значение М(Ех(Б+И)>0) составляет 7,0% Для зоны воздействия «Центр-Восток» ареал совместных превышений КН соединений серы и азота занимает менее 5% от площади группы рецепторных участков с Р = 100%

Таким образом, приемлемость атмотехногенного воздействия для зоны воздействия «Центр-Восток» показывают результаты и детерминистического, и вероятностного расчетов Вместе с тем, для зоны воздействия «Запад» с вероятностью свыше 99% размер ареала превышений КН больше допустимого критерия приемлемости Наиболее вероятно, что величины М(Ех(Б)>0) и М(Ех(М)>0) составят от 5 до 10% площади данной группы выделов, величина М(Ех(8+^>0) - от 15 до 20%

Таблица 3

Вероятности совместных превышений КН соединений серы и азота (Р(Ех(Б+Ы)>0), %) для наземных экосистем зоны атмотехногенного воздействия объектов ГП Венесуэлы

Группа рецепторных участков М(Ех(Б+1Ч)>0), %

0-5 5-10 10-15 15-20 20-25 25-100

Зона воздействия «Запад» <0,1 <0,1 <0,1 92,1 7Д 0,0

Зона воздействия «Центр-Восток» 100 0,0 0,0 0,0 0,0 0,0

Зона атмотехногенного воздействия объектов ГП Венесуэлы 100 0,0 0,0 0,0 0,0 0,0

Примечание жирным шрифтом выделены значения Р(Ех(8+Щ>0), соответствующие детерминистическим значениям площадей ареалов превышений КН (табл 1)

Неопределенность выполненных оценок в первую очередь связана с пробелами в исходных данных Основными факторами неопределенности являются 1) отсутствие данных о пространственном распределении показателей выпадений химических элементов (по данным мониторинга или моделирования), 2) ограниченные данные о содержании азота и основных катионов в биомассе растительности и соотношении различных пулов биомассы в различных типах леса, 3) отсутствие данных о продуктивности травянистых экосистем Это снижает корректность оценок величин выпадений и депонирования основных катионов и азота в древесной биомассе

Моделирование величин КН и их превышений с помощью метода Монте-Карло позволило использовать все доступные данные по экосистемам территории исследования и экосистемам-аналогам и повысить, таким образом, достоверность расчетов (ошибка моделирования не превышает 3,4%)

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

Управление риском

На основе величин КН соединений ссры и азота был произведен расчет допустимых уровней поступления соединений серы и азота ( А^и'асф ^скр(жс)) в наземные экосистемы зоны атмотехногенного воз,действия объектов ГП Венесуэлы согласно формулам (3) и (4).

г(5) (3)

А -С!

N

(¡ер{асс)"

ат«(М),с1ПШ{М)>аП111(щ

(4)

Показатель А^ер<асс) определяет допустимый совместный объем выпадений кислотообразующих соединений.

Моделирование по методу Монте-Карло показало, что для экосистем зоны атмотехногенного воздействия объектов ГГ1 Венесуэлы с вероятностью 99% величина показателя А^ер(асс) не превысит 4068 экв./га/год, показателя ЫсЬ>р(асс) ~ 4598 экв./га/год. С помощью эмпирических функций распределения величин А<1ер(асс) и №с1ер(цсс) были рассчитаны вероятности их попадания в определенные интервалы значений. Результаты расчетов иллюстрируют диаграммы на рис. 3. Показано, что с наибольшей вероятностью величина А&р(асс) составит ог 2500 до 3000 экв./га/год, величина ^(^р(асс) ~ 50 ДО 500 экв./га/год.

?4000

ц 3600-4000 о

а 3000-3500

5 2500-3000 а

й 2000-2500

"8 1500-2000 и

"Е 5000-1500 а>

§ 500-1000 <500

17 14 13 13

а)

*45О0 В 2

ч 4000-4600 1 1

2 3600-4000 СП 8 а

¡9 3000-3500 Г-1 17

0 2500-3000 С3 8 2000-2500 'СЗ 11

1 1600-2000 I-18

1000-1500 г-112

г 500-1000 |—| я

<500 I

6)

10 20 30 40 50 60 70 80 80 100 Р,%

10 20 30 40 50 60 70 80 00 100

Рис, 3. Вероятностное распределение значений Лйер(асс) (а) и ^аер(асс) (6) для наземных экосистем зоны атмотехногенного воздействия объектов ГП Венесуэлы

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

В главе 4 изложены результаты апробации предлагаемых подходов к оценке экосистемных рисков на локальном уровне в рамках оценки воздействия на окружающую среду расширения Средне-Тиманского бокситового рудника (СТБР), расположенного в Княжпогостском районе Республики Коми

Строительство рудника ведется с 1998 г на запасах Вежаю-Ворыквинского месторождения бокситов Поскольку до строительства рудника и сопутствующей инфраструктуры район месторождения не был освоен в хозяйственном отношении, большая часть территории исследования представлена ненарушенными или слабонарушенными ландшафтами Выбросы вредных веществ в атмосферу, в особенности пылевые выбросы, рассматриваются в качестве значимого экологического аспекта деятельности данного предприятия Помимо этого, территория рудника входит в состав зоны потенциального воздействия выбросов магистрального газопровода «Ямал-Центр»

Идентификация опасности

К числу основных видов деятельности, определяющих выбросы СТБР и связанные с ними воздействия на качество атмосферного воздуха и биоту, относятся добычные работы в карьерах, транспортировка и дробление боксита В составе выбросов рудника преобладают пыль неорганическая с различным содержанием 8Ю2, СО, N0*, 802, сажа, летучие органические соединения

Демонстрационное исследование охватывало фрагмент зоны атмотехногенного воздействия СТБР площадью 8623 га В качестве реципиентов воздействия выступили наземные экосистемы Почвенный покров территории исследования формируют в основном низкогорные глееподзолистые почвы в сочетании с низкогорными торфянисто-подзолисто-глееватыми и низкогорными торфяно-подзолисто-глеевыми почвами Основными типами растительности являются еловые зеленомошные и долгомошные леса, а также верховые болота, чередующиеся с вторичными елово-березовыми лесами

Ландшафтную структуру территории исследования слагают 12 типов урочищ В ее границах имеются участки водоохранных зон (ВОЗ) двух водотоков (р Ворыквы, руч Черного) и санитарно-защитной зоны (СЗЗ) СТБР, которые составляют, соответственно, 25% и 22% от рассматриваемой площади Экосистемы в пределах ВОЗ были отнесены нами к особо ценным реципиентам воздействия

В рамках данного демонстрационного исследования оценка экосистемных рисков проводилась для двух сценариев воздействия Сценарий 1 соответствует уровню аэрогенного поступления ЗВ до начала разработки

месторождения Сценарий 2 характеризует современный (по состоянию на конец 2003 г) уровень атмотехногенной нагрузки на экосистемы

В силу низкой буферной емкости и низкого содержания основных катионов в почвах, а также низкой биологической продуктивности биоценозов, наземные экосистемы территории исследования чувствительны к нагрузке кислотообразующих и эатрофирующих соединений серы и азота (Ж)*, 802) Другими значимыми химическими факторами опасности для экосистем являются микроэлементы, отличающиеся высоким фоновым содержанием в почвах рассматриваемой территории и накапливающиеся в почвах, растительности и снежном покрове зоны аэрогенного загрязнения рудника По данным многолетнего мониторинга геологической среды, такими элементами являются Мп, N1, Си, А1, Ре, Т1, V, Хп, Со, Сг, РЬ

В качестве приоритетных ЗВ были выбраны соединения серы и азота, а также три приоритетных ТМ РЬ, Си и Zn Детальная оценка экосистемного риска проводилась для следующих групп рецепторных участков а) каждого из 12 типов урочищ, б) фрагмента СЗЗ СТЕР, в) территории ВОЗ

Оценка экспозиции

Расчет величин аэрогенного поступления химических элементов и соединений в экосистемы-реципиенты осуществлялся на основе данных об их накоплении в снежном покрове за 2000-2003 гг Для сценария 2 на основе расчетных величин выпадений в каждой точке опробования были получены поля выпадений Расчет проводится на основе данных о климатических параметрах, определяющих рассеивание примесей в атмосфере, и источниках выбросов с помощью УПРЗА «Гарант-Универсал», версия 3 0 Данные таблицы 4 иллюстрируют расчетные величины ВС^ер, С1с/ер, 8с[ер, РЪ^ер, Сщер, 2пс!ер

Таблица 4

Величины выпадений основных катионов (ВС&р), анионов хлора (С1с/ер), серы (5^), азота (Л^) и приоритетных тяжелых металлов для района

СТБР

всаер С1(1ер РЬйер Си,1ср 2п<1ер

экв./га/год г/га/год

Сценарий 1 506 и = 7 172 п = 7 128 и = 7 206 п = 7 2,3 п = 7 4,1 п = 7 45,1 п = 7

Сценарий 2 280-560 п = 13 65-320 п = 13 100-220 п = 13 40-340 п — 13 0,5-18,5 п =28 0,5-56,5 п=28 5,0-185,0 п = 28

Примечание для сценария 1 приведены средние значения параметров, которые использовались для детерминистического расчета величин КН и их превышений, для сценария 2 - максимальные и минимальные значения каждого параметра, п - число пунктов опробования

Для сценария 2 величины Sjep превысили фоновые значения на 99% территории исследования, Njep - на 84% территории Фоновые значения выпадений превышены на 62% территории исследования по РЬ, 80% территории - по Си, и 17% территории - по Zn

Оценка эффектов

Расчет максимальной критической нагрузки серы (CLmax(S')) и максимальной критической нагрузки азота (CLmax(N')) проводился для двух рассматриваемых сценариев воздействия, поскольку эти величины зависят от уровня аэрогенного поступления анионов хлора (Cljep) и основных катионов (ВС&р) Величины критической нагрузки «питательного» азота (CLnut(N)) и критических нагрузок TM (CL(Pb), CL(Cu), CL(Zn)) не зависят от величин выпадений

Проведенные расчеты показали некоторое снижение величин CLmax(S) и CLmax(N) при современном уровне атмотехногенной нагрузки (сценарий 2), что обусловлено ростом величин СЦер, который недостаточно компенсируется поступлением основных катионов из атмосферы

Рис 4 иллюстрирует пространственное распределение величин максимальной КН серы для наземных экосистем, рассчитанных для сценария 2 Значения CLmax(S) возрастают от пойменных травянистых и болотных экосистем (менее 400 экв /га/год) к экосистемам еловых лесов, занимающим возвышенные хорошо дренированные участки (более 1500 экв/га/год) Минимальные значения CLmax(S) (менее 250 экв /га/год) получены для пойменных экосистем с преобладанием древесной растительности (ивняки), где значительный пул основных катионов «расходуется» на создание биомассы, что снижает потенциал этих экосистем в отношении нейтрализации кислотных выпадений

Пространственная динамика величин критической нагрузки питательного азота (CLnut(N')) в значительной степени зависит от типа растительности Устойчивость рассматриваемых экосистем к эвтрофирующей составляющей выпадений азота возрастает от экосистем хвойных лесов на почвах легкого механического состава (менее 450 экв /га/год) к заболоченным лесам (500-550 экв /га/год) и болотам (более 550 экв /га/год) Ареал максимальных значений (более 750 экв /га/год) приурочен к участкам заболоченных ивняков

Наименее чувствительны к суммарной нагрузке азота леса на переувлажненных торфяно-подзолисто-глеевых почвах (CLmax(N) свыше 4500 экв /га/год) благодаря сочетанию высоких показателей денитрификации и депонирования азота в древесной биомассе Участки с наименьшими значениями CLmax(N) (менее 1000 экв/га/год) совпадают с ареалами максимальной чувствительности к кислотным выпадениям (см выше)

С1тах(5), экв./гаЛод Н '250

250 - 500

условные обозначения

500 -1000 1000 - 1250 1250- 1500 =•1500

границы типов урочищ

границы сзкитарно-защитной зоны СТБР

границы водоохранных зон

границы рудных тел (I. II, III, МЖБ)

Рис, 4. Величины максимальной критической нагрузки серы для

наземных экосистем района СТБР (сценарий 2)

Величины К11 ТМ определяются в основном почвенными характеристиками (рН, содержанием органического вещества и др.). Наиболее заметна зависимость величин КН от содержания органического вещества для свинца: прослеживается четкое подразделение территории исследования на участки с торфяно-подзолистф-глеевыМИ и болотными почвами (СЦРЬ) свыше 30 г/га/год) и другими разновидностями почв (9-20 г/га/год). Вторым важным фактором степени защищенности экосистем от поступления ТМ является интенсивность выноса последних с внутрипочвенным стоком. С этим связаны высокие значения КН меди и цинка для лссных экосистем дренированных местообитаний (СЬ(Си) - от 40 до 60 г/га/год, СЦ2п) - от 400 до 450 г/га/год). Отмечено снижение величин КН с увеличением рН почвенного раствора до 5 и выше в ряду болотных почв, наиболее выраженное для меди. Это обусловлено уменьшением концентрации растворенного органического вещества, от которой зависит концентрация свободных ионов ТМ в почвенном растворе, оказывающих токсичное действие на биоту.

Характеристика риска

Для наземных экосистем района СТБР были выбраны следующие критерии приемлемости воздействия:

• для территории ВОЗ - 100%-ая защищенность экосистем (М(Ех(Х)>0) = 0%);

• для каждого из 12 типов урочищ - 95%-ая защищенность экосистем (М(Ех(Х)>0) < 5%);

• для территории СЗЗ - 50%-ая защищенность экосистем (М(Ех(Х)>0) < 50%).

В качестве вариабельных параметров выступали величины Всир1, РЬир,, Сиир,, ¿Пцр,, которые моделировались как независимые величины на основе данных о составе древостоя лесов территории исследования, а также геохимические параметры: рН и содержание органического вещества в почве {[ОМ]$). Массивы значений вариабельных параметров сформированы на основе данных нолевых геоботанических и геоэкологических исследований.

Детерминистический расчет превышений КН не выявил участков с превышениями КН при фоновых величинах атмотехногенной нагрузки (сценарий 1).

По состоянию на конец 2003 г., некоторые рецепторные участки наземных экосистем испытывали нагрузки соединений серы и азота, свинца и меди, превышающие рас чет но-допустимые уровни. Пространственное распределение величин превышений КН соединений серы и азота (Ех(5+М)) иллюстрирует рис. 5.

Рис. 5. Величины совместных превышении КН соединений серы и азота (Ех(8+Щ) для наземных экосистем района СТБР (сценарий 2)

Площадь ареалов превышений КН приоритетных ЗВ невелика, она составляет для соединений серы и азота - менее 2% от общей территории исследования, для свинца - 4%, для меди - около 1%. Вместе с тем, превышения КН выявлены для территории ВОЗ, для которой необходимо обеспечить 100%-ную защищенность экосистем. Установленные критерии приемлемости атмотехногеиного воздействия превышены для следующих групп рецепторных участков:

• для типов урочищ 6 и 10, территории ВОЗ - по соединениям серы и азота (величины М(Ех($+М)>0) равны 5,1%, 18,4% и 3,5% соответственно);

• для типов урочищ 2 и 9 - по свинцу (величины М(Ех(РЬ)>0) равны 13,0 и 14,0% соответственно),

• для типа урочищ 1, территории ВОЗ - по меди (величины М(Ех(Си)>0) равны 19,0 и 1,2% соответственно)

Вероятностный расчет совместных превышений КН соединений серы и азота для наземных экосистем показал, что для всех рецепторных участков с превышениями КН расчетные вероятности Ех(Б+М)>0 (р') были равны 1, те результаты детерминистического и вероятностного расчетов данного показателя идентичны

В случае свинца и меди вероятностные значения площадей ареалов превышений КН для всех групп рецепторных участков оказались выше значений, полученных в результате детерминистического расчета Так, для типов урочищ 2 и 9 вероятностные величины М(Ех(РЬ)>0) составили 13,5% и 22,8% соответственно Вероятностные величины М(Ех(Си)>0) достигли 38,6% для типа урочища 1 и 4,3% для территории ВОЗ

Для большинства типов урочищ распределения величин М(Ех(РЬ)>0) и М(Ех(Си)>0) не подчиняются нормальному закону, для них были вычислены эмпирические функции экосистемного риска Для крупных массивов рецепторных участков (например, типов урочищ 2 и 4, территории ВОЗ) показана корректность использования функций нормального распределения для расчета вероятностных показателей экосистемного риска

По сравнению с результатами расчетов для зоны атмотехногенного воздействия объектов ГП Венесуэлы, разброс значений превышений КН, полученных в результате вероятностного расчета, гораздо выше было выделено 9 классов значений М(Ех(РЬ)>0) и М(Ех(Си)>0) Расчетные вероятности превышения КН свинца (Р(Ех(РЬ)>0)) для некоторых групп рецепторных участков, представлены в табл 5

Для групп рецепторных участков, которые испытывают атмотехногенную нагрузку ТМ выше установленного критерия приемлемости по данным детерминистического расчета, вероятности достижения детерминистических значений превышений КН свинца и меди составили

• Р(Ех(РЬ)>0) для типа урочища 2 - 22,7%, для типа урочища 9 - 35,2%,

• Р(Ех(Си)>0) для типа урочища 1 - 23,3%, для территории ВОЗ - 99,9% Вероятность того, что в пределах ВОЗ не будут отмечены превышения КН

данных элементов, равна нулю

Основными факторами неопределенности расчетов риска для экосистем района СТБР стали 1) оценка нагрузки приоритетных ЗВ на основе косвенных данных (данных атмогеохимических исследований) и достаточно высокая погрешность (до 30%) расчетных величин выпадений химических элементов, 2) ограниченные данные о геохимических характеристиках некоторых типов почв (например, болотных остаточно-низинных), 3) оценка депонирования основных катионов, азота и ТМ в древесной биомассе на основе данных литературных

источников Это может быть причиной заниженных величин КН, в частности КН соединений серы и азота, для некоторых рецепторных участков

Таблица 5

Вероятности превышений критических нагрузок свинца (Р(Ех(РЬ)>0), %) для наземных экосистем района СТЕР (сценарий 2)

М(Ех(РЪ)>0), %

Группа рецепт, участков 0-5 5-10 10-15 15-20 20-25 25-30 30-50 50100

Тип урочища 2 0,5 76,7 22,7 <0,1 <0,1 <0,1 0,0 0,0

Тип урочища 4 40,0 59,9 <0,1 <0,1 <0,1 <0,1 <0,1 0,0

Тип урочища 9 9,0 45,0 35,2 9,9 0,8 <0,1 <0,1 <0,1

Территория СЗЗ СТБР <0,1 30,4 68,9 0,6 <0,1 <0,1 <0,1 0,0

Территория ВОЗ 97,5 2,4 <0,1 <0,1 <0,1 0,0 0,0 0,0

Расчетная площадка в целом 64,0 35,9 <0,1 <0,1 <0,1 <0,1 0,0 0,0

Примечания / - жирным шрифтом выделены значения Р(Ех(РЬ)>0), соответствующие детерминистическим значениям площадей ареалов превышений КН, 2 — в таблице представлены результаты расчетов для типов урочищ, для которых детерминистические значения М(Ех(РЪ)>0)> 1%

Согласно результатам многофакторного анализа, влияние параметров, связанных с типом растительности, на значения Ех(Я+И) для наземных экосистем гораздо ниже, чем для почвенных характеристик Для приоритетных ТМ почвенные характеристики также оказались гораздо более значимыми, чем величины долговременного депонирования ТМ в древесной биомассе Следовательно, в данном случае последний фактор неопределенности не может считаться значимым В целом, выбранные модели нелинейной регрессии с высокой степенью достоверности описывают зависимости превышений КН приоритетных ЗВ от входных параметров (значения коэффициентов детерминации Я выше 0,9) Основной вклад в величину ошибки результата вносят погрешности входных параметров, ошибка метода - не более 2%

На основе результатов оценки экосистемных рисков предложены организационно-технические решения по снижению атмотехногенных воздействий для территории исследования, включая предложения по организации и проведению экологического мониторинга для СТЕР, а также рекомендации по проведению оценки экосистемных рисков в зоне влияния СТБР на стадии его расширения

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

1 Разработаны методические подходы к оценке рисков для экосистем, затронутых выбросами объектов газовой промышленности, на основе методологии критических нагрузок ЗВ Предложены

• алгоритм расчета показателей экосистемного риска на основе моделирования величин превышений КН по методу Монте-Карло,

• схема оценки экосистемных рисков, связанных с выбросами ЗВ на производственных объектах,

• метод определения допустимых уровней атмотехногенной нагрузки приоритетных ЗВ для зон воздействия выбросов от производственных объектов

2 С помощью разработанных методических подходов выполнена оценка рисков для наземных экосистем, расположенных в зоне атмотехногенного воздействия объектов ГП Венесуэлы и в районе Средне-Тиманского бокситового рудника (СТЕР) (Республика Коми) в рамках работ по экологическому обоснованию расширения этих производственных объектов

• охарактеризованы экосистемы-реципиенты, произведено их ранжирование, определены приоритетные ЗВ (соединения серы и азота для экосистем Венесуэлы, соединения серы и азота, тяжелые металлы (РЬ, Си, Хх\) для экосистем района СТБР),

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

• рассчитаны величины КН приоритетных ЗВ для реципиентов,

• проведены детерминистический и вероятностный расчеты превышений КН приоритетных ЗВ, выведены функции экосистемного риска для групп рецепторных участков наземных экосистем, проведен анализ вероятностного распределения площадей ареалов превышений КН (М(Ех(Х)>0))

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

4 Обоснованы допустимые уровни поступления соединений серы и азота <Ас1ер(асф ^ер(асс)) Для зоны атмотехногенного воздействия объектов ГП Венесуэлы, данные показатели предложены в качестве экологических ограничителей при планируемом развитии газовой отрасли

5 Подготовлены рекомендации по смягчению атмотехногенного воздействия на экосистемы района СТБР, включая предложения по проведению оценки экосистемных рисков для зоны влияния выбросов проектируемых объектов рудника

Основные публикации по теме диссертации:

1 Демидова О А Современное состояние системы экологического страхования в Российской Федерации // Экономика природопользования, ВИНИТИ -2005 -№5 -С 27-59

2 Демидова О А Перспективы использования методов анализа риска в рамках процедуры ОВОС // Проблемы безопасности и чрезвычайных ситуаций -2007 - № 1 - С 47-54

3 Демидова О А Оценка экосистемных рисков для производственных объектов на основе методологии критических нагрузок // Экология и промышленность России -2007 - №3 -С 50-52

4 Башкин В Н , Самсонов Р О , Казак А С , Припутина И В , Демидова О А , Танканаг А В Использование ГИС-технологий для оценки величин критических нагрузок оксидов азота в зонах влияния газотранспортной системы Венесуэлы // Материалы 8-й Всероссийской конференции «Геоинформатика в нефтегазовой и горной отраслях», 27 февраля — 3 марта 2007 г, Москва -ГИС-Ассоциация, 2007 - 1 электрон опт диск (CD-ROM)

5 Демидова О А , Башкин В Н , Юркин Е А , Котова М В , Припутина И В Оценка экосистемных рисков, связанных расширением Средне-Тиманского бокситового рудника // Проблемы анализа риска - 2007 - № 1 (принято в печать)

6 Derrndova О , Cherp A Use of Risk Assessment for Addressing Human Health Implications in EIAs of Selected UK Waste Incineration Projects // Proceedings of the 2nd International Conference on Ecological Chemistry, 11-12 October, 2002, Chisinau, Moldova - Chisinau MRDA, 2002 - P 282-284

7 Demidova О , Bashkin V The Application of Critical Load And Level Approach for Assessing Ecosystem Risks in EIA and EA Practitioner's Approach // Proceedings of the 6th Subregional Meeting on Effect-Oriented Activities in the Countries of Eastern and South-Eastern Europe Relevant to UN ECE Convention on LRTAP, Moscow-Puschino, Russia, 16-18 September, 2004 - Smolensk Magenta, 2004 -P 110-118

8 Demidova О, Cherp A Risk Assessment for Improved Treatment of Health Considerations in EIA // Environmental Impact Assessment Review - 2005 -Vol 25(4) -P 411-429

9 Demidova О, Bashkin V Critical Load and Level Approach for Assessment of Ecosystem Risks in EIA Promises and Challenges // Power, Poverty and Sustainability The role of Impact Assessment, the 26th annual conference of the International Association for Impact Assessment, 23-26 May 2006, Stavanger, Norway Abstracts - IAIA, 2006 - P 46

10 Bashkin VN, Demidova OA Biogeochemical Approaches to Environmental Risk Assessment // Encyclopedia of Ecology / Editor-in-chief S E Jorgensen -Copenhagen Elsevier, 2007 - MS number 768 (in press)

Подписано к печати 9 04 2007 Заказ № №78530086 Тираж 100 экз 1 уч - изд л , ф-т 60x84/16 Отпечатано в ООО «ВНИИГАЗ» по адресу 142717, Московская область, Ленинский р-н, п Развилка, ООО «ВНИИГАЗ»

Введение.

Глава 1. Принципы и методы оценки атмотехногенных воздействий объектов газовой промышлености на экосистемы.

1.1. Концептуальные основы изучения воздействий производственных объектов на окружающую среду.

1.2. Нормирование атмотехногенных воздействий объектов газовой промышленности на окружающую среду.

1.2.1. Объекты газовой промышленности как источники атмотехногенных воздействий.

1.2.2. Принципы и методы нормирования атмотехногенных воздействий

1.3. Методы оценки риска для экологического обоснования проектов промышленного строительства (на примере объектов газовой промышленности).

1.4. Современные методы оценки воздействия химических факторов экологической опасности на экосистемы.

1.4.1. Индикаторы состояния экосистем как реципиентов воздействия

1.4.2. Методы прогнозирования воздействия на экосистемы.

1.4.3. Оценка значимости воздействий на экосистемы, нормирование экосистемных эффектов.

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

1.5.1. Источники неопределенности в прогнозе воздействий на экосистемы.

1.5.2. Методы оценки рисков для экосистем, связанных с загрязнением окружающей среды.

1.5.3. Методология критических нагрузок поллютантов и оценка риска для экосистем.

Глава 2. Методика количественной оценки экосистемных рисков, связанных с выбросами загрязняющих веществ на объектах газовой промышленности.

2.1. Концепция и алгоритм оценки экосистемных рисков на основе величин критических нагрузок поллютантов.

2.2. Схема проведения оценки экосистемных рисков, связанных с выбросами производственного объекта, на основе величин критических нагрузок поллютантов.

2.3.1 Идентификация опасности.

2.3.2 Оценка экспозиции.

2.3.3 Оценка эффектов.

2.3.4 Характеристика риска.

2.3. Источники информации для оценки экосистемных рисков и требования к исходным данным.

Глава 3. Оценка экосистемных рисков для зоны атмотехногенного воздействия объектов газовой промышленности Венесуэлы.

3.1 Введение.

3.2 Объект и методы исследования.

3.3 Идентификация опасности.

3.3.1 Характеристика факторов экологической опасности.

3.3.2 Определение и характеристика реципиентов.

3.3.3 Характеристика воздействий.

3.4 Оценка экспозиции.

3.5 Оценка эффектов.

3.6 Характеристика риска.

3.7 Управление риском.

Глава 4. Оценка экосистемных рисков для фрагмента зоны атмотехногенного воздействия Средне-Тиманского бокситового рудника (СТЕР) (Республика Коми).

4.1 Введение.

4.2 Объект и методы исследования.

4.3 Идентификация опасности.

4.3.1 Характеристика факторов экологической опасности.

4.3.2 Определение и характеристика реципиентов.

4.3.3 Характеристика воздействий.

4.4 Оценка экспозиции.

4.5 Оценка эффектов.

4.5.1 Критические нагрузки соединений серы и азота для наземных экосистем.

4.5.2 Критические нагрузки приоритетных ТМ для наземных экосистем

4.6 Характеристика риска.

4.6.1 Детерминистический расчет превышений КН приоритетных ЗВ для наземных экосистем.

4.6.2 Вероятностный расчет превышений КН приоритетных ЗВ для наземных экосистем.

4.7 Рекомендации по смягчению атмотехногенных воздействий на наземные экосистемы.

Выводы.

Актуальность темы исследования. Строительство и эксплуатация производственных объектов ведут к преобразованию природных комплексов (экосистем) в зонах воздействия. Чрезмерная эксплуатация и деградация «экосистемных услуг», таких, как поддержание запасов питьевой воды, качества воздуха и воды, регулирование регионального климата (Ecosystems., 2005), создает угрозу для жизнеобеспечения общества.

Минимизация атмотехногенных воздействий на природную среду приобретает особое значение для газовой отрасли в связи с ростом объемов добычи, транспортировки, переработки и потребления газа а, следовательно, и объемов эмиссии поллютантов. Так, в 2005 г. валовые выбросы загрязняющих веществ (ЗВ) в атмосферу на объектах газовой промышленности России достигли 2351 тыс. тонн (Экологический., 2005). Свыше 10% от этого объема составляют оксиды азота и серы, которые признаны значимыми химическими факторами опасности для экосистем.

В силу структурной и функциональной сложности экосистем результаты прогноза изменений в их состоянии характеризуются высокой степенью неопределенности. Низкая достоверность прогноза воздействий может стать причиной развития негативных эффектов, масштаб и интенсивность которых значительно превысят ожидаемые показатели (эффект «бомб замедленного действия» (Stigliani, 1991)). Такие явления можно классифицировать как чрезвычайные ситуации (ЧС) техногенного характера.

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

Цель работы - разработка методов количественной оценки экосистемных рисков в зонах воздействия выбросов ЗВ от объектов газовой промышленности для обеспечения экологической безопасности их эксплуатации и предотвращения ЧС техногенного характера.

Основные задачи исследования:

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

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

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

4. Анализ результатов оценки рисков и разработка рекомендаций по управлению рисками для экосистем территорий исследования.

В качестве объектов демонстрационных исследований выступали наземные экосистемы:

• зоны атмотехногенного воздействия объектов газовой промышленности (ГП) Боливарианской Республики Венесуэла (ГП Венесуэлы) (региональный уровень);

• территории размещения объектов Средне-Тиманского бокситового рудника (СТБР) в Республике Коми (локальный уровень).

Территория СТБР испытывает комплексное воздействие выбросов объектами СТБР и компрессорных станций газопровода «Ямал - Центр».

Предмет исследования - методы количественной оценки негативных воздействий на экосистемы, связанных с выбросами ЗВ на технологических объектах.

Научная новизна. Разработаны новые методические подходы к количественной оценке экосистемных рисков, связанных с выбросами соединений серы, азота и некоторых тяжелых металлов (ТМ), развивающие положения методологии критических нагрузок (КН) поллютантов, направленные на предотвращение развития ЧС техногенного характера. Предложено теоретическое обоснование и разработан алгоритм расчета показателей экосистемного риска на основе вероятностного моделирования величин КН и их превышений по методу Монте-Карло. Обоснована схема проведения оценки экологических рисков в рамках работ по оценке воздействий действующих или проектируемых производственных объектов на окружающую среду. Предложен метод определения допустимых уровней атмотехногенной нагрузки приоритетных ЗВ для зон воздействия выбросов ЗВ от технологических объектов.

Впервые выполнен расчет величин КН соединений серы и азота и их превышений для наземных экосистем Южной Америки (в границах зоны атмотехногенного воздействия объектов ГП Венесуэлы).

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

Защищаемые положения:

1. Концепция количественной оценки экосистемных рисков, связанных с выбросами производственных объектов, на основе величин КН поллютантов.

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

3. Обоснование схемы проведения оценки экосистемных рисков, связанных с выбросами соединений серы, азота и ТМ, при оценке воздействий на окружающую среду действующих и/или проектируемых производственных объектов газовой промышленности.

4. Метод определения безопасных уровней атмотехногенной нагрузки приоритетных ЗВ для экосистем, затронутых выбросами объектов газовой промышленности.

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

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

Результаты диссертационного исследования использованы при разработке Генеральной схемы развития газовой промышленности Боливарианской Республики Венесуэла и при экологическом обосновании расширения СТБР.

Апробация и публикация результатов исследований. Материалы исследований по теме диссертации докладывались и обсуждались на заседаниях кафедры геохимии ландшафта и географии почв Географического факультета МГУ им. Ломоносова; на 24-й и 26-й конференциях Международной ассоциации оценки воздействия (International Association for Impact Assessment) в Ванкувере (Канада) в 2004 г. и в Ставанге (Норвегия) в 2006 г.; на Шестой субрегиональной встрече по деятельности в области эффект-ориентированных методов оценки под эгидой Конвенции о трансграничном загрязнении воздуха на большие расстояния ЕЭК ООН (Москва-Пущино, 2004 г.); на 8-й Всероссийской научно-практической конференции «Геоинформатика в нефтегазовой и горной отраслях» (Москва, 2007 г.). По материалам диссертации опубликовано 10 работ, в том числе 3 работы - в журналах, включенных в «Перечень.» ВАК Минобрнауки РФ.

Структура и объем диссертации. Диссертация состоит из введения, четырех глав, выводов и списка литературы из 195 наименований. Работа изложена на 172 страницах машинописного текста и содержит 38 рисунков, 12 таблиц и 4 приложения.

ВЫВОДЫ

1. Разработаны методические подходы к оценке рисков для экосистем, затронутых выбросами объектов газовой промышленности, на основе методологии критических нагрузок ЗВ. Предложены:

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

• схема оценки экосистемных рисков, связанных с выбросами ЗВ на производственных объектах;

• метод определения допустимых уровней атмотехногенной нагрузки приоритетных ЗВ для зон воздействия выбросов от производственных объектов.

2. С помощью разработанных методических подходов выполнена оценка рисков для наземных экосистем, расположенных в зоне атмотехногенного воздействия объектов ГП Венесуэлы и в районе Средне-Тиманского бокситового рудника (СТБР) (Республика Коми) в рамках работ по экологическому обоснованию расширения этих производственных объектов:

• охарактеризованы экосистемы-реципиенты, произведено их ранжирование, определены приоритетные ЗВ (соединения серы и азота для экосистем Венесуэлы, соединения серы и азота, тяжелые металлы (Pb, Си, Zn) для экосистем района СТБР);

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

• рассчитаны величины КН приоритетных ЗВ для реципиентов;

• проведены детерминистический и вероятностный расчеты превышений КН приоритетных ЗВ, выведены функции экосистемного риска для групп рецепторных участков наземных экосистем, проведен анализ вероятностного распределения площадей ареалов превышений КН (М(Ех(Х)>0)).

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

4. Обоснованы допустимые уровни поступления соединений серы и азота dep(acc)> ^dep(acc)) зоны атмотехногенного воздействия объектов ГП

Венесуэлы; данные показатели предложены в качестве экологических ограничителей при планируемом развитии газовой отрасли.

5. Подготовлены рекомендации по смягчению атмотехногенного воздействия на экосистемы района СТБР, включая предложения по проведению оценки экосистемных рисков для зоны влияния выбросов проектируемых объектов рудника.

1. Акимов В.А., Лесных В.В., Радаев Н.Н. Риски в природе, техносфере, обществе и экономике. - М.: Деловой экспресс, 2004. - 352 с.

2. Базилевич Н.В. Биологическая продуктивность экосистем Северной Евразии. М.: Наука, 1993.-292 с.

3. Байдаков С.А., Серов Г.П. Правовое обеспечение охраны окружающей среды и экологической безопасности. М.: Анкил, 2003. - 464 с.

4. Башкин В.Н. Биогеохимия. М.: Научный мир, 2004. - 584 с.

5. Башкин В.Н. Оценка степени риска при критических нагрузках загрязняющих веществ на экосистемы // География и природные ресурсы. 1999. -№ 1. - С. 35-39.

6. Башкин В.Н. Управление экологическим риском. М.: Научный мир, 2005. -368 с.

7. Башкин В.Н., Евстафьева Е.В., Снакин В.В. и др. Биогеохимические основы экологического нормирования. М.: Наука, 1993. - 312 с.

8. Башкин В.Н., Казак А.С., Припутана И.В., Горлов Д.В. Оценка экологического риска при модернизации газопроводной системы «Средняя Азия Центр» // Охрана окружающей среды в нефтегазовом комплексе. - 2006. - № 5. - С. 5-13.

9. Башкин В.Н., Курбатов А.А., Припутана И.В. Показатели критических нагрузок вместо ПДК // Экология и промышленность. 2005. - №8. - С. 25-29.

10. Башкин В.Н., Курбатова А.С., Савин Д.С. Методологические основы оценки критических нагрузок поллютантов на городские экосистемы. М.: НИиПИЭГ, 2004. -64 с.

11. Биопродукционный процесс в лесных экосистемах Севера/ К.С. Бобкова, Э.П. Галенко, И.В. Забоева и др. СПб.: Наука, 2001. - 277 с.

12. Боровиков В.П. Программа Statistica для студентов и инженеров: 2-е изд. -М.: КомпьютерПресс, 2001. 301 с.

13. Босняцкий Г.П., Гриценко, А.И., Седых А.Д. Проблемы экологического мониторинга в газовой промышленности. М.: АО «Ника-5», 1993. - 80с.

14. Букс И.И., Фомин С.А. Экологическая экспертиза и оценка воздействия на окружающую среду (ОВОС): учеб. пособие. Кн.1. - М.: Изд-во МНЭПУ, 1999 -127 с.

15. Быков А.А., Мурзин Н.В. Проблемы анализа безопасности человека, общества и природы. СПб.: Наука, 1997. - 247 с.

16. Вараксин А.Н. Статистические модели регрессионного типа в экологии и медицине. Екатеринбург, 2006. - 256 с.

17. Веретенников А.В. Эколого-биологические основы повышения продуктивности таёжных лесов Европейского Севера. JI.: Наука, 1981. - 232 с.

18. Возмещение экологического ущерба (правовые и экономические аспекты проблем прошлого экологического ущерба) / под ред. С.А. Боголюбова и И.Н. Сенчени. М.: Изд-во НУМЦ, 2001. - 144 с.

19. Воробейчик E.JI., Садыков О.Ф., Фарафонтов М.Г. Экологическое нормирование техногенных загрязнений наземных экосистем: (локальный уровень). -Екатеринбург: Наука, 1994. 280 с.

20. Втюрин Г.М. Структура почвенного покрова таежной зоны европейского Северо-Востока. JI. Наука 1991. - 148 с.

21. Говорушко С.М. Экологическое сопровождение хозяйственной деятельности. Владивосток: Дальнаука, 2003. 271 с.

22. Голубев Г.Н. Геоэкология: учеб. для студентов вузов. М.: Геос, 1999. - 338 с.

23. ГОСТ Р ИСО 14001-98. Системы управления окружающей средой требования и руководство по применению: утв. Постановлением Госстандарта РФ от 21.10.1998 № 378. Введ. 01.04.1999. -М.: ИПК Издательство стандартов, 1998.

24. Гриценко А.И., Акопова Г.С., Максимов В.М. Экология, нефть и газ. Наука, М., 1997.-598 с.

25. Данилов-Данильян В.И., Залиханов М.Ч., Лосев К.С. Экологическая безопасность. Общие принципы и российский аспект. М.: Изд-во МНЭПУ. - 2001. -328 с.

26. Дмитриев В.В., Фурмин Г.Т. Экологическое нормирование и устойчивость природных систем: учеб. пособие. СПб.: Наука: РГГМУ, 2004. - 294 с.

27. Добровольский В.В. Основы биогеохимии: учеб. пособие. М.: Высшая школа, 1998-412 с.

28. Дьяконов К.Н., Дончева А.В. Экологическое проектирование и экспертиза: учебник для вузов. М.: Аспект Пресс, 2002. - 384 с.

29. Дьяконов К.Н., Звонкова Т.В. Теоретические, методологические и организационные вопросы эколого-географических экспертиз // Основы эколого-географи-ческой экспертизы / Под ред. К.Н. Дьяконова, Т.В. Звонковой. М. Изд-во МГУ, 1992.-С. 5-26.

30. Забоева И.Р. Почвы и земельные ресурсы Коми АССР. Сыктывкар: Коми книж. изд-во, 1975. - 344 с.

31. Ильин В.Б. Тяжелые металлы в системе почва-растение / отв. ред. И.Л. Клевенская; АН СССР, Сиб. отд-ние, Ин-т почвоведения и агрохимии. -Новосибирск: Наука: Сиб. отд-ние, 1991. 150 с.

32. Исаченко А.Г. Ландшафтоведение и физико-географическое районирование: учеб. для геогр. спец. ун-тов. М.: Высшая школа, 1991. - 365 с.

33. Концепция экологического нормирования в Российской Федерации // Материалы выездного совещания МПР России / МПР России. М., 2001. - 30 с.

34. Коплан-Дикс В.А., Алехова М.В. К вопросу разработки экологических нормативов качества атмосферного воздуха // Труды НИИ Атмосфера 2002 г. / НИИ Атмосфера 2002. - Режим доступа: http://iaap.narod.ru/stKDiks.htm. - Загл. с экрана.

35. Копцик Г.Н., Макаров М.И., Киселев В.В. Принципы и методы оценки устойчивости почв к кислотным выпадениям. М.: Изд-во МГУ, 1998. - 96 с.

36. Кузьмин И.Н., Махутов Н.А., Хетагуров С.В. Безопасность и риск: эколого-экономические аспекты. СПб.: Изд-во СПбГУЭФ, 1997. - 164 с.

37. Ч. 2. Апатиты: КНЦ РАН, 1996. - 196 с.

38. Методика определения ущерба окружающей природной среде при авариях на магистральных нефтепроводах: утв. Минтопэнерго РФ 01.11.1995 // СПС «КонсультантПлюс».

39. Методические рекомендации по анализу и управлению риском воздействия на здоровье населения вредных факторов окружающей среды / А.А. Быков, Л.Г. Соле-нова, Г.М. Земляная, В.Д. Фурман. М.: Анкил, 1999. - 72 с.

40. Методология оценки состояния экосистем: учеб. пособие / О.М. Кожова, JI.P. Изместьева, Б.К. Павлов и др. Ростов-на-Дону: Изд-во ООО «ЦВВР», 2000. -128 с.

41. О промышленной безопасности опасных производственных объектов: Федеральный закон от 21 июля 1997 г. № 116-ФЗ: принят ГД ФС РФ 20 июня 1997 г.: ред. от 9 мая 2005 г. // Собрание законодательства РФ. 1997. - № 30. - Ст. 3588.

42. О рыбохозяйственных нормативах: приказ Госкомрыболовства РФ от 28 апреля 1999 г. № 96. М., Госкомрыболовство РФ, 1999.

43. Об охране окружающей среды: Федеральный закон от 10 января 2002 г. № 7-ФЗ: принят ГД ФС РФ 20 декабря 2001 г.: ред. от 18 декабря 2006 г. с изм. и доп., вступившими в силу с 1 января 2007 г. // Собрание законодательства РФ. 2002. -№2.-Ст. 133.

44. Об экологической экспертизе: Федеральный закон от 23 ноября 1996 г. № 174-ФЗ: принят ГД ФС РФ 19 июля 1995 г.: ред. от 18 декабря 2006 г. с изм. и доп., вступившими в силу с 1 января 2007 г. // Собрание законодательства РФ. 1995. -№48.-Ст. 4556.

45. Отчет по объекту «Организация и ведение мониторинга геологической среды на Вежаю-Ворыквинском месторождении бокситов в 2001-2002 гг.» / Министерство природных ресурсов и охраны окружающей среды Республики Коми, ОАО «Боксит

46. Тимана», ООО «Геолог-1»; исполн. В.И. Безуглов, Н.В. Куделина, О.В. Афанасенко и др.; № ГР 89-01-7/2. Ухта, 2003. - 196 с.

47. Охрана окружающей среды. Книга 15 // Средне-Тиманский бокситовый рудник I очередь строительства. Проект / Коми Научный центр УрО РАН, Институт биологии; рук. А.И. Таскаев. - Сыктывкар, 2001.

48. Оценка и регулирование качества окружающей природной среды / Под ред. А.Ф. Порядина, А.Д. Хованского. М.: Прибой, 1996. - 348 с.

49. Оценка риска как инструмент социально-гигиенического мониторинга / Б.А. Кацнельсон, Л.И. Привалова, С.В. Кузьмин и др. Екатеринбург: Издательство АМБ, 2001.-261 с.

50. Положение о государственном природном ихтиологическом заказнике республиканского значения «Вымский»: утв. Постановлением Совета Министров Республики Коми от 1 марта 1993 г. № 110 // СПС «КонсультантПлюс».

51. Почвенная карта района расположения Средне-Тиманского бокситового рудника (I очередь): рисунок 6.1 / Оценка воздействия на окружающую среду.

52. Пояснительная записка. Книга 12 // Средне-Тиманский бокситовый рудник -I очередь строительства. Проект / Коми Научный центр УрО РАН, Институт биологии; рук. А.И. Таскаев; отв. исполн.: Т.Е. Евдокимова. 1:25 ООО. - Сыктывкар, 2001.

53. Почвы Коми АССР. М.: АН СССР, Коми филиал, 1958. - 198 с.

54. Проведение мониторинга состояния наземных экосистем в зоне влияния Средне-Тиманского бокситового рудника: отчет о результатах НИР / Коми Научный центр УрО РАН, Институт биологии; отв. исполнитель С.В. Дегтева. Сыктывкар, 2003.- 110 с.

55. Проведение мониторинга состояния наземных экосистем в зоне влияния Средне-Тиманского бокситового рудника: отчет о результатах НИР / Коми Научный центр УрО РАН, Институт биологии; отв. исполнитель С.В. Дегтева. Сыктывкар, 2002.- 128 с.

56. Пузаченко Ю.Г. Методологические основы географического прогноза и охраны среды. М.: УРАО, 1998. - 212 с.

57. РД 03-418-01. Методические указания по проведению анализа риска опасных производственных объектов: утв. Постановлением Госгортехнадзора России от 10.07.2001 № 30 //СПС «КонсультантПлюс». Введ. 01.10.2001.

58. РД 153-34.0-002-01. Временная методика оценки ущерба, возможного вследствие аварии гидротехнического сооружения: утв. Приказом Минэнерго РФ от 26.04.2001 № 130 // СПС «КонсультантПлюс».

59. Росновский И.Н. Устойчивость почв в экосистемах как основа экологического нормирования / отв. ред. В.Н. Воробьев. Томск: Изд-во Ин-та оптики атмосферы СО РАН, 2001-251 с.

60. Рубцов Д.М. Глеево-подзолистые почвы Притиманья. Сыктывкар, 1968. --41 с.

61. Сафонов В.С, Одишария Г.Э., Швыряев А.А. Теория и практика анализа риска в газовой промышленности. М.: НУМЦ Минприроды, 1996. - 207 с.

62. СНиП 11-02-96. Инженерные изыскания для строительства. Основные положения: утв. Постановлением Минстроя России от 29.10.1996 № 18-77 // ИС «СтройКонсультант». Взамен СНиП 1.02.07-87; введ. 01.11.1996.

63. Соболь И.М. Метод Монте-Карло. М.: Наука, 1985. - 64 с.

64. СП 11-102-97. Инженерно-экологические изыскания для строительства: утв. Письмом Госстроя России от 10.07.1997 № 9-1-1/69, согл. письмом Госкомэкологии РФ от 25.02.1997 № 02-12/23-577 // ИС «СтройКонсультант». Введ. 15.08.1997.

65. Справочник по теории вероятности и математической статистике / B.C. Королюк, Н.И. Портенко, А.В. Скороход, А.Ф. Турбин. М.: Наука. Главная редакция физико-математической литературы, 1985. - 640 с.

66. Устойчивость экосистем к эмиссиям магистральны газопроводов: / В.Н. Башкин, А.С. Казак, В.В. Снакин, И.П. Припутана, В.Р. Хрисанов, Б.И. Кочуров. Москва-Смоленск: Универсум, 2002. - 232 с.

67. Швыряев И.А. Оценка воздействия продуктов превращений выбросов компрессорных станций на природные объекты Крайнего Севера // Вестник РУДН. Сер. «Экология и безопасность жизнедеятельности». -2006. -№1 (14)-Р. 14-22.

68. Экоаудиторское заключение по деятельности предприятия ОАО «Боксит Тимана» / ООО «Экосеверсервис». Ухта, 2002. - 42 с.

69. Экологическая безопасность автомобильной дороги: понятие и количественная оценка: утв. распоряжением Минтранса РФ от 31.12.2002 № OC-1181-p // ИС «СтройКонсультант». Введ. 31.12.2002.

70. Экологическая оценка и экологическая экспертиза / О.М. Черп, В.Н. Вини-ченко, М.В. Хотулева и др. 3-е изд., пер. и доп. - М.: Социально-экологический союз, 2001.-312 с.

71. Экологические проблемы газовой промышленности / А.И. Гриценко, Г.П. Босняцкий, Ю.С. Шилов, А.Д. Седых; Всерос. НИИ природ, газов и газовых технологий (ВНИИГАЗ). М.: ВНИИГАЗ, 1993. - 94 с.

72. Экологический отчет ОАО «Газпром» за 2005 г. Электронный ресурс. // Экология. Режим доступа: http://www.gazprom.ru/documents/ Gazprom EcologyRus2005.pdf. - Загл. с экрана.

73. Экология подземного хранения газа / Э.Б. Бухгалтер, Е.В. Дедиков, Л.Б. Бухгалтер и др. М.: МАИК «Наука/Интерпериодика», 2002. - 431 с.

74. A Guide to Risk Assessment and Risk Management for Environmental Protection / United Kingdom Department of the Environment. London: HMSO Books, 1995. - 92 p.

75. Adler, A., Verwijst, Т., Aronsson, P. Estimation and Relevance of Bark Proportion in a Willow Stand // Biomass and Bioenergy. 2005. - Vol. 29. - P. 102-113.

76. Ale, B.J.M. Dealing with Risk in Environmental Policy in the Netherlands // Risk Assessment. Proceedings of an HSE Conference on Risk Assessment, 6-9 October 1992, London. Bootle, Merseyside: Health and Safety Executive, UK, 1992. - P. 65-72.

77. Alriksson, A. and Eriksson, H.M. Distribution of Cd, Cu, Pb and Zn in Soil and Vegetation Compartments in Stands of Five Boreal Tree Species in N.E. Sweden // Water, Air, and Soil Pollution: Focus 2001. - Vol. 1. - P. 461-475.

78. Alriksson, A. and Eriksson, H.M. Variations in Mineral Nutrient and С Distribution in the Soil and Vegetation Compartments of Temperate Tree Species in N.E. Sweden // Forest Ecology and Management. 1998. - Vol. 108. - P. 261-273.

79. Andrews, R.N.L. Environmental Impact Assessment and Risk Assessment: Learning from Each Other // Environmental Impact Assessment: Theory and Practice/ ed. by P. Wathern. L.: Unwin Hyman, 1990. - P. 85-97.

80. Arquiaga, М.С., Canter, L., and Nelson, D.I. Risk Assessment Principles in Environmental Impact Studies // Environmental Professional. 1992. - Vol. 14 (3). -P. 201-219.

81. Belyazid, S., Sverdrup, H., Alveteg, M. The Biogeochemical Models Family of the Biogeochemistry Group at Lund University // Proceedings of the 6th Subregional Meeting On Effect-Oriented Activities in the Countries of Eastern and South-Eastern Europe

82. Relevant to UN ECE Convention of LRTAP, Moscow-Puschino, Russia, September 16-18, 2004. Moscow: Institute for Urban Ecology, 2004 - P. 46-55.

83. Burger, J. and Gochfeld, M. On Developing Bioindicators for Human and Ecological Health // Environmental Monitoring and Assessment. 2001. - Vol. 66. - P. 23-46.

84. Cairns, J., McCormick, P.V., and Niederlehner B.R. A Proposed Framework for Developing Indicators of Ecosystem Health // Hydrobiologia. 1993. - Vol. 236. - P. 1-44.

85. Calculation of Critical Loads for Cadmium, Lead and Mercury. Background document to a Mapping Manual on Critical Loads of Cadmium, Lead and Mercury: Alterra Report 1104 / W. de Vries, G. Schutze, S. Lofts et al. Wageningen: Alterra, 2004. - 143 p.

86. Canter, L.W. Environmental Impact Assessment. 2nd ed. - N.Y.: McGraw-Hill, 1996.-660 p.

87. Canter, L.W. Methods for Effective EIA Practice // Environmental Methods Review: Retooling Impact Assessment for the New Century / ed. by A.L. Porter and J.J. Fittipaldi. -Fargo, ND: The Press Club, 1998. P. 58-68.

88. Carpenter, R.A. Risk assessment // Environmental and Social Impact Assessment / ed. by F. Vanclay and D.A. Bronstein. Chichester: Wiley, 1996. - P. 193-219.

89. Chen, C.S. Ecological Risk Assessment for Aquatic Species Exposed to Contaminants in Keelung River, Taiwan // Chemosphere. 2005. - Vol. 61. -P. 1142-1158.

90. Clark, J.M., Chapman, P.J., Adamson, J. K. and Lane S.N. Influence of Drought-Induced Acidification on the Mobility of Dissolved Organic Carbon in Peat Soils // Global Change Biology. 2005. - Vol. 11. - P. 791-809.

91. Crane, M., Grosso, A., Whitehouse, P., and Forrow, D. Risk Characterization in Direct Toxicity Assessment of The River Esk and The Tees Estuary // Ecotoxicology. -2004.-Vol. 13(5).-P.463-474.

92. Dale, V.H. and Beyeler, S.C. Challenges in the Development and Use of Ecological Indicators. Ecological Indicators. - 2001. - Vol. 1. - P. 3-10.

93. Dalva, M. and Moore, T.R. Sources and Sinks of Dissolved Organic Carbon in a

94. Forested Swamp Catchment // Biogeochemistry. 1991. - Vol. 15(1). - P. 1-9.

95. Demidova, O. and Cherp, A. Risk Assessment for Improved Treatment of Health Considerations in EIA // Environmental Impact Assessment Review. 2005. - Vol. 25(4). -P. 411-429.

96. Derivation of Partition Relationships to Calculate Cd, Cu, Ni, Pb and Zn Solubility and Activity in Soil Solutions: Alterra Report 305 / P.F.A.M. Romkens, J.E. Groenenberg, L.T.C. Bonten et al. Wageningen: Alterra, 2004. - 72 p.

97. Dezzeo, N. and Chacon, N. Nutrient Fluxes in Incident Rainfall, and Stemflow in Adjacent Primary and Secondary Forests of the Gran Sabana, Southern Venezuela // Forest Ecology and Management. 2006. - P. 234:218-226.

98. Ecosystem Health New Goals for Environmental Management / ed. by R. Costanza, B. Norton and B. Haskell. - Washington, DC: Island Press, 1992. - 269 p.

99. Ecosystems and Human Well-being: Synthesis. A report of the Millennium Ecosystem Assessment / Millennium Ecosystem Assessment Board. -Washington DC: Island Press, 2005.-166 p.

100. Eduljee, G. Risk assessment // Handbook on Environmental Impact Assessment: 2 volumes / ed. by J. Petts. Oxford: Blackwell Science Ltd., 1999. - Vol. 1. - P. 374-404.

101. Emission Reductions in the Natural Gas Sector Through Project-Based Mechanisms:

102. A Information Paper Электронный ресурс. / International Energy Agency. 2003. -Режим доступа: http://www.iea.org/textbase/papers/2003/devbase.pdf. - Загл. с экрана.

103. Environmental Hazard and Risk Assessment // Environmental Assessment Sourcebook / The World Bank. Washington, DC: World Bank, 1997. - Update No. 21. - 10 p.

104. Environmental Impact Assessment Training Resource Manual / United Nations Environmental Programme. 2nd edition. - Geneva: UNEP-DTIE Economics and Trade Branch, 2002. - 573 p.

105. Erickson, P.A. A Practical Guide to Environmental Impact Assessment. San Diego: Academic Press Inc., 1994. - 266 p.

106. Fairman, R., Mead, C.D., and Williams, P.W. Environmental Risk Assessment: Approaches, Experiences and Information Sources. Copenhagen: European Environment Agency, 1999.-219 p.

107. Guidance on EIA Scoping / European Commission. - Luxembourg: Office for Official Publications of the European Communities, 2001 - 35 p.

108. Guidelines for Ecological Risk Assessment / US Environmental Protection Agency. EPA/630/R-95/002F. Washington, DC: U.S. EPA, 1998. - 124 p.

109. Health Aspects of Environmental Assessment // Environmental Assessment Sourcebook / The World Bank. Washington, DC: World Bank, 1997. — Update No. 18. - 10 p.

110. Henriksen, A., Dillon, P.J., and Aherne, J. Critical Loads of Acidity for Surface Waters in South-Central Ontario, Canada: Regional Application of the Steady-State Water

111. Chemistry (SSWC) Model // Canadian Journal of Fisheries and Aquatic Sciences. 2002. -Vol. 59.-P. 1287-1295.

112. Holl, K.D., Cairns, J., Jr. Landscape Indicators in Ecotoxicology // Handbook of Ecotoxicology, ed. by D.J. Hoffman, B.A. Rattner, G.A. Burton et al. Boca Raton, F.L.: Lewis Publishers, 1995.-P. 185-197.

113. Intensive Monitoring of Forest Ecosystems in Europe. Technical Report 2003 / W. de Vries, G.J. Reinds, M. Posch et al.; UN/ECE and EC, Forest Intensive Monitoring Coordinating Institute. Geneva and Brussels, 2003.

114. Kirchhoff, D. and Doberstein, B. Pipeline risk assessment and risk acceptance criteria in the State of Sao Paulo, Brazil // Impact Assessment and Project Appraisal. 2006. -Vol. 24(3).-P. 221-234.

115. Kolluru, R.V., Bartel, S., and Stricoff, S. Risk Assessment and Management Handbook for Environmental, Health and Safety Professionals. NY: McGraw-Hill, 1996. -556 p.

116. Kooistra, L., Leuven, R.S.E.W., Wehrens, R., Nienhuis, P.H., and Buydens, L.M.C. A Procedure for Incorporating Spatial Variability in Ecological Risk Assessment of Dutch River Floodplains // Environmental Management. 2001. - Vol. 28 (3). - P. 359-373.

117. Kuylenstierna, J.C., Cambridge, H., Cinderby, S., and Chadwick, M. Terrestrial Ecosystem Sensitivity to Acidic Deposition in Developing Countries // Water, Air and Soil Pollution. 1995. - Vol. 85. - P. 2319-2324.

118. Lackey, R.T. Ecological Risk Analysis // Fundamentals of Risk Analysis and Risk Management, ed. by V. Molak. Boca Raton, F.L.: Lewis Publishers, 1997. - P. 87-97.

119. Lee, N. and George, С. Environmental Assessment in Developing and Transitional Countries: Principles, Methods and Practice. Chichester: Wiley, 2000. - 314 p.

120. Leith, H. Modelling the Primary Productivity of the World // Primary Productivity of the Biosphere. Ecological Studies 14 / ed. by H. Leith and R.H. Whittaker. NY.: Springer-Verlag, 1975.-P. 237-263.

121. Lemly, D.A. and Richardson, C.J. Guidelines for Risk Assessment in Wetlands // Environmental Monitoring and Assessment. 1997. - Vol. 47. - P. 117-134.

122. Leuven, R.S.E.W. and Poudevigne, I. Riverine Landscape Dynamics and Ecological Risk Assessment // Freshwater Biology. 2002. - Vol. 47. - P. 845-865.

123. Lewis, W.M., Hamilton, S.K., Jones, S., and Runnels, D. Major Element Chemistry, Weathering and Element Yields for the Caura River Drainage, Venezuela // Biogeochemistry. 1987.-P. 4: 159-181.

124. Livingston, P.A., Aydin, K., Boldt, J., Ianelli, J., and Jurado-Molina, J. A Framework for Ecosystem Impacts Assessment Using an Indicator Approach // ICES Journal of Marine Science. 2005. - Vol. 62. - P. 592-597.

125. Lofts, S., Spurgeon, D.J., Svendsen, C., and Tipping, E. Deriving Soil Critical Limits for Cu, Zn, Cd and Pb: a Method Based on Free Ion Concentrations // Environ. Sci. and Technol. 2004. - Vol. 38 (13). - P. 3623-3631.

126. Lorenz, C.M., Gilbert, A.J., and Cofino, W.P. Indicators for Transboundary River Basin Management // Environmental Indices: System Analysis Approach; ed. by Y.A. Pykh, D.E. Hyatt, R.J.M. Lenz. Oxford, UK: EOLSS Publishers Co. Ltd., 1999. - P. 313-328.

127. Machado, A., Przybylski, J., Montiel, V. et al. Efectos de La Composition Del Gas Natural Comprimido en Los Componentes Regulados de Los Gases de Escape En un Motorde Encendido рог Chispa. Caracas: INCI, Die., 2001, - Vol. 26, No. 12. - P. 591-596.

128. Manahan, S.E. Environmental Chemistry: 7th edition. Boca Raton: Lewis Publishers, 2000. - 898 p.

129. Modelling and Mapping of Critical Thresholds in Europe: Status Report 2001 / ed. by M. Posch, P.A.M. de Smet, J.-P. Hettelingh, R.J. Downing. Bilthoven, The Netherlands: RIVM, 2001.- 188 p.

130. Models Электронный ресурс. // The Biogeochemistry Group 2002. - Режим доступа: http://www2.chemeng.lth.se/models/index.shtml. - Загл. с экрана.

131. Models and Methods Электронный ресурс. // The Coordination Center for Effects. 2006. - Режим доступа: http://mv4v.mnp.n1/cce/methmod/#tcm:42-16260. - Загл. с экрана.

132. Montes, R., San Jose J.J. Chemical Composition and Nutrient Loading by Precipitation in the Trachypogon Savannas of the Orinoco Llanos, Venezuela // Biogeochemistry. -1989. -№. 7. P. 214-256.

133. Morales, J., Bifano, C., Escalona, A. Atmospheric Deposition of SO4-S and (NH4 + N03)-N at Two Rural Sites in the Western Maracaibo Lake Basin, Venezuela // Atmospheric Environment. 1998. - Vol. 32. - № 17. - P. 3051-3058.

134. Morales, J., Sanchez, L., Velasquez, H., and De Borrego, B. Nutrient Loading by Precipitation in the Maracaibo Lake Basin, Venezuela // Water, Air and Soil Pollution. -2001.-Vol. 130.-P. 511-516.

135. Organisation, Committee on Challenges to Modern Society Pilot Study, co-ordinated by E.F. Verheyen, K. Nagels and M. Coenen. Antwerpen: University of Antwerpen, 1997. -P. 49-53.

136. Naito, W., Miyamoto, K., Nakanishi, J., Masunaga, Sh., and Bartell, S.M. Evaluation of an Ecosystem Model in Ecological Risk Assessment of Chemicals // Chemosphere. -2003.-Vol. 53.-P. 363-375.

137. PDVSA Business Plan for 2004-2009 Электронный ресурс. // Publications and Papers // International Energy Agency. 2004. - Режим доступа: http://www.iea.org/ Textbase/work/2004/opec/pdvsa.pdf. - Загл. с экрана.

138. Petts, J. and Eduljee, G. Environmental Impact Assessment for Waste Treatment and Disposal Facilities. Chichester: Wiley, 1994. - 485 p.

139. Petts, J. Risk Assessment and Management for Waste Treatment and Disposal // Handbook of Environmental Risk Assessment and Management / ed. by P. Calow. -Oxford: Blackwell Science, 1998. P. 417-452.

140. Planes Estrategicos de PDVSA: Agosto 2005 Электронный ресурс. // Plan Siembra Petrolera 2005-2030 // Sitio oficial de Petroleos de Venezuela S.A. 2005. -Режим доступа: http://www.pdvsa.com/. - Загл. с экрана.

141. Power, М., and McCarthy, L.S. A Comparative Analysis of Environmental Risk Assessment/Risk Management // Environmental Science and Technology. 1998. -Vol. 32.-P. 224-231.

142. Publications Электронный ресурс. // The Coordination Center for Effects. 2004. - Режим доступа: http://www.mnp.nl/cce/publ/. - Загл. с экрана.

143. Ramos, Т.В., Caeiro, S., de Melo, J.J. Environmental Indicator frameworks to Design and Assess Environmental Monitoring Programs // Impact Assessment and Project Appraisal. 2004. - Vol. 22(1). - P. 47-62.

144. Risk Management and ELA: Research Needs and Opportunities / A.P. Grima, P. Timmerman, C.D. Fowle, and P. Byer. Quebec: Canadian Environmental Assessment Research Council, 1986. - 18 p.

145. Risk: analysis, perception, management / Royal Society. L.: Royal Society of London, 1992.-200 p.

146. Sander, P., Bergback, В., and Oberg, T. Uncertain Numbers and Uncertainty in the Selection of Input Distributions: Consequences for a Probabilistic Risk Assessment of Contaminated Land // Risk Analysis. 2006. - Vol. 26(5). - P. 1363-1375.

147. Sanhueza, E., Arias, M.C., Donoso, L., Graterol, N., Hermoso, M. Marti, I., Romero, J., Rondon, A., and Santana, M. Chemical Composition of Acid Rains in the Venezuelan Savannah Region // Tellus. 1992. - 44B. - P. 54-62.

148. Schroeder, L.D., Sjoquist, D.L., and Stephan, P.E. Understanding Regression Analysis: An Introductory Guide. Beverly Hills, CA: Sage Publications, 1986. - 96 p.

149. Serveiss V.B. Applying Ecological Risk Principles to Watershed Assessment and Management // Environmental Management. 2002. - Vol. 29. - P. 145-154.

150. Smith, E.P. and Cairns, J., Jr. Extrapolation Methods for Setting Ecological Standards for Water Quality: Statistical and Ecological Concerns // Ecotoxicology. 1993. -Vol. 2.-P. 203-219.

151. Smrchek, J.C and Zeeman, M.G. Assessing Risks to Ecological Systems from Chemicals // Handbook of Environmental Risk Assessment and Management / ed. by P. Calow. Oxford: Blackwell Science, 1998. - P. 417-452.

152. Stigliani, W.M. Chemical Time Bombs: Definition, Concepts and Examples: Executive Report 16 / International Institute for Applied Systems Analysis, Laxenburg, Austria. Laxenburg: IASA, 1991. - 23 p.

153. The Digital Soil Map of the World Электронный ресурс. / Food and Agriculture Organization of the United Nations (FAO). Version 3.5. - FAO, 1995.

154. The Plan for Sowing the Oil Электронный ресурс. // The New PDVSA Contact. -2005. №3. - 20 p. - Режим доступа: http://www.pdvsa.com/interface.en/database/ fichero/publicacion/925/19.PDF. - Загл. с экрана.

155. Treweek, J. Ecological Impact Assessment. Oxford, UK: Blackwell Science, 1999. -351 p.