Оценка остроты экологической ситуации при крупномасштабном эколого-геохимическом картировании городов с применением математического моделирования тема диссертации и автореферата по ВАК РФ 04.00.02, кандидат геолого-минералогических наук Иванов, Юрий Владимирович

• Иванов, Юрий Владимирович
• кандидат геолого-минералогических наук
• 1999

Оглавление диссертации кандидат геолого-минералогических наук Иванов, Юрий Владимирович

Введение.

1. Методика.

2. Экологические исследования в г. Россошь и на прилегающей территории.

2.1. Общая характеристика.

2.2. Результаты опробования почвы.

2.3. Результаты опробования снега.

2.4. Результаты опробования верхнего подземного водоносного горизонта

2.5. Результаты медицинских исследований.

2.6. Результаты математического моделирования.

3. Экологические исследования в г. Ливны

3.1. Общая характеристика.

3.2. Результаты опробования почвы.

3.3. Результаты опробования снега.

3.4. Результаты медицинских исследований

3.5. Результаты математического моделирования.

4. Экологические исследования в г. Воронеж

4.1. Общая характеристика.

4.2. Результаты опробования почвы.

4.3. Результаты опробования снега.

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

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

Введение диссертации (часть автореферата) на тему «Оценка остроты экологической ситуации при крупномасштабном эколого-геохимическом картировании городов с применением математического моделирования»

АКТУАЛЬНОСТЬ ПРОБЛЕМЫ. Несмотря на значительный общий спад промышленного производства в последние годы уровень техногенной нагрузки на окружающую среду, особенно в городах, сохраняется на достаточно высоком уровне. Наиболее эффективным методом оценки величины техногенной нагрузки является геохимическая съёмка различных депонирующих сред. Так, для определения загрязненности атмосферы в качестве таковых обычно используются почва и снег. При этом, почва отражает в большей степени постоянное многолетнее загрязнение, а состав снеговой пыли - сезонное, современное. Оценка уровней загрязненности в этом случае производится по степени превышения содержания определяемых ингредиентов в почве или снеговой пыли над их предельно допустимыми (ПДК) или природными концентрациям.

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

В связи с появлением и распространением в последнее десятилетие ПЭВМ исследователи для оценки экологической ситуации получили возможность привлекать дополнительный фактический материал, обработка которого требует значительного 9 объема вычислений. Наиболее привлекательной в этом плане выглядит идея оценки опасности достигнутого загрязнения окружающей среды, особенно в городах, по уровню заболеваемости населения. Основным недостатком большинства работ, основанные на этом принципе яценки ОЭС, является их относительная незавершенность. Дело, как правило, кончается или расчетом парных коэффициентов корреляции меж^у частотой отдельных заболеваний и содержанием некоторых загрязняющих компонентов, или оценкой роли различных факторов в изменчивости того или иного заболевания. То есть, в конечном итоге, даже при значимых связях заболеваемости с загрязнением, исследования не доводились до логического конца - разбраковки конкретных геохимических аномалий по их экологической опасности.

ЦЕЛЬ И ЗАДАЧИ ИССЛЕДОВАНИЙ. Целью настоящей работы является демонстрация на большом фактическом материале эффективности оценки ОЭС с применением математического моделирования в предлагаемом варианте.

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

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

10 в соответствии с оптимизированными уравнениями регрессии перевести загрязненность в каждой точке опробования в уровень заболеваемости и построить карту рассчитанной заболеваемости ( для нескольких заболеваний -карту средней рассчитанной заболеваемости); провести анализ карт техногенной нагрузки и рассчитанной заболеваемости с оценкой ОЭС в пределах каждой аномалии техногенной нагрузки; провести анализ рассчитанных и фактических уровней заболеваемости по участкам; для случаев сложного функционального строения составить карты зонирования по ОЭС с учетом результатов опробования нескольких депонирующих сред, отражающих загрязненность атмосферы,

НАУЧНАЯ НОВИЗНА. В настоящей работе впервые проводится оценка ОЭС в гг. Россошь, Ливны и Воронеж методом математического моделирования по результатам геохимического картирования различных депонирующих сред с учетом фактического уровня заболеваемости детского населения. При этом, элементы новизны имеют не только полученные результаты, но и сам метод в предлагаемом варианте.

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

ПРАКТИЧЕСКАЯ ЗНАЧИМОСТЬ РЕЗУЛЬТАТОВ РАБОТЫ. В результате районирования территории указанных городов по остроте экологической ситуации с учетом величины и характера техногенной нагрузки и уровня заболеваемости детского населения получена достаточно объективная и достоверная информация, необходимая для принятия решений в различных областях, начиная от рационального землепользования и планов генеральной застройки, до страховой медицины. Нагдяд

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

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

АПРОБАЦИЯ РАБОТЫ. Применение математического моделирования для оценки ОЭС в предлагаемом варианте использовалось автором при экологических исследованиях за период с 1989 по 1997гг. так же в г Нововоронеж, р.ц. Грибанов-ский - Воронежской области; гг. Данков, Елец - Липецкой области; Красногорск -Московской области; Краснопресненском районе г. Москвы. В 1993г. данный метод, изложенный-в конкурсной работе по теме: «Разработка критериев и показателей для районирования территорий по степени экогеохимической нагрузки» был признан Г лав н ы м^нау чно-техническим управлением Минэкологии России эффективным при районировании городских и прилегающих территорий по степени остроты экологической ситуации (отзыв ГНТУ N 246 от 1.04.93., хгрил. ).

Всего, по теме диссертации было опубликовано работ, не считая производственных ^ научно - производственных отчетов. Результаты исследований были представлены на Международном симпозиуме: «Контроль и реабилитация окружающей среды» (^Гомск, 1998); на Международной конференции: «Геоэкологические проблемы устойчивого развития городской среды» (Воронеж, 1996); на Всероссийской конференции: «Геологическое картографирование» (Москва, 1998); на региональной конференции: «Проблемы интеграции экологической и хозяйственной политики в Черноземном Центре России» (Мичуринск, 1995); на городской научно - производственны^ конференции: «Высокие технологии в практике учреждений здравоохранения г. Воронежа»(Воронеж, 1995); на научно-практической конференции, посвященной 70-летию Воронежского биосферного государственного заповедника (Воронеж, 1997).

ФАКТИЧЕСКИЙ МАТЕРИАЛ. В основу работы положен материал, собранный за периоды работы автором с 1986 по 1993гг. главным геохимиком Комплексного

12 производственного предприятия (КПП) «Геокомплекс», а с 1993 по 1998гг. - научным сотрудником кафедры исторической геологии Воронежского госуниверситета.

КПП «Геокомплекс», бывшая Центральная опытно - методическая экспедиция (ЦОМЭ), проводило среднемасштабные исследования во многих городах центральной части России, включая л Воронеж, и, так называемого теперь «ближнего зарубежья». Воронежский университет с 1992г. проводит комплексные исследования масштаба 1: 200000 в пределах листов N-3 7-ХVI,XXII. В процессе этих исследований собран большой фактический материал, дополнительная математическая обработка косорото показала на .достаточную эффективность метода моделирования при оценке ОЭС городских территорий.

СТРУКТУРА И ОБЪЕМ РАБОТ. Диссертация состоит из 5 глав, введения, заключения, списка литературы. Общий объем работы состоит из 140 страниц машинописного текста, 4 рисунков, 24 таблиц, 28 приложений. Список использованной литератууы содержит 51 наименование.

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

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

В третьей главе приводятся результаты аналогичных исследований в г. Ливны. Выявлен^! особенности проявления техногенной нагрузки в загрязненности почвы, снега и заболеваемости населения, характерные для малых городов.

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

13

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

При написании настоящей работы на различных этапах решающая помощь автору была оказана: доктором геолого-минералогических наук профессором А. Д. Савко, создавшему необходимые условия для ее проведения; доктором геолого-минералогических наук профессором В. Л. Бочаровым, убедившему автора в необходимости этой работы и оказывающему ему постоянные консультации; сотрудником КПП «Геокомплекс» В. П. Лепендиным - в сборе информации.

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

14

1. МЕТОДИКА.

Метод оценки техногенной нагрузки по содержанию микроэлементов в депонирующих средах при площадных экологических исследованиях унаследован в несколько адаптированном виде от поисковой геохимии. А так как объекты и методы опробования, аналитическая база и принципы камеральной обработки результатов аналогичны таковым при поисковых работах, то методической основой проводимых экологических исследований являются инструкции и методические указания по геохимической съёмке [36] . Особенности экологических исследований отражены в отраслевых рекомендациях [33, 34].

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

Опробование почвы и снегового покрова проводилось по сети с переменной, как правило, плотностью.

Так, в окрестностях г. Россошь [4], на площади 20x20 км., исключая застроенную о

часть города, плотность опробования составляла 1 пр./км' - для почвы и 0.5 пр./км" для снега. Непосредственно на территории города площадью 60 км2 почвенные пробы

2 2 отбирались с плотностью 4 пр/км , а снеговые - 2 пр/км . Для оценки специализации загрязнения в районе отдельных предприятий дополнительно отбиралось по 5-10 почвенных и снеговых проб. Загрязненность грунтовых вод оценивалась по результатам анализа проб, отобранных из колодцев. Всего на территории города и прилегающих окрестностей в 1990-91гг. отобрано около 700 почвенных, 200 снеговых и 50 водных проб.

В г. Ливны и на непосредственно прилегающей к нему территории [5], на площади 40 км2, в 1992г. было отобрано около 500 почвенных и 55 снеговых проб с постоянной плотностью около 12 пр/км2 - для почвы и 1 пр/км2 - для снега.

В г. Воронеж [48, 49], на его правобережной части, на площади около 250 км , включая прилегающие районы, отобрано 1000 почвенных проб в 1985г., 300 снеговых

15 л проб в 1986г. и 600 - в 1987г. со средней плотностью 4 пр./ км - для почвы и около 2,5 пр./ км - для снега.

Отбор почвенных проб производился с поверхности на глубину до 5 см. методом конверта, когда каждая проба состояла из 5ти частных, отобранных по углам и из центра квадрата 10x10м. Вес объединенной пробы 450 - 500 г.

Снеговые пробы сечением 0,16-0,25 м отбирались по всей мощности снегового покрова в конце зимы в местах, где влияние вторичных аэрозолей сведено к минимуму. Вес снеговых проб в зависимости от мощности снегового покрова составлял 15-20кг.

Предварительная обработка почвенных проб сводилась к их просушиванию, просеиванию через сито с диаметром ячей 0,7-1,0мм.с последующим квартованием до получения 100 граммовых навесок, которые отправлялись в лабораторию.

Предварительная обработка снеговых проб включала в себя: растапливание снега, фильтрацию талой воды вначале через сито с диаметром ячей 0,1мм., а затем - через беззольный фильтр, просушивание, взвешивание, прокаливание вместе с фильтром отфильтрованного осадка при Т - 600 ° С, истирание прокаленного остатка и отправку его в лабораторию.

Анализ почвенных и снеговых проб проводился полуколичественным спектральным методом на 36 основных микроэлементов (ТМ): отобранных в г. Россошь - в лаборатории Белгородской экспедиции, в г. Дивны - в лаборатории ПГО «Воронежгеология», в г. Воронеж - для почвенных проб - в лаборатории Бронницкой экспедиции, для снеговых - в лаборатории Александровской экспедиции. Водные пробы из колодцев, отобранные в г. Россошь, анализировались на ТМ в лаборатории ПГО «Воронежгеология», а на макрокомпоненты - химическим методом в лаборатории Воронежского университета.

Основным критерием определения уровня содержания того или иного ингредиента являются его природные ( фоновые ) концентрации в аналогичных средах опробования, а основным базовым показателем - коэффициент концентрации (KJ, вычисляемый по формуле: = с/сф, о) где: СА - содержание элемента в конкретной пробе, Сф - фоновое (природное) содержание этого элемента (табл. 1).

16

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

Таблица 1

Фоновые содержания ТМ в почве и снеге,

Наименован. ТМ Единиц. изм.(%) Фоновые содержания тяжелых металлов.

Почва Снег

Россошь Ливны Воронеж Россошь Ливны Воронеж

Медь 1*10-3 1.5 2.18 3,25 21.8 21.8 41,0

Цинк -1-1- 3.25 5.79 6,62 88.1 88.1 94,2

Свинец -1-1- 0.83 2.41 1,5 78.6 78.6 33,4

Никель -1-1- 2.13 2.3 2,45 2.67 2.67 8,2

Кобальт -1-1- 0.65 0.65 0,7 0.22 0.22 1,4

Хром -1-1- 7 4.3 3,6 13.1 13.1 чэ о

Ванадий -1-1- 8.5 5.2 4,0 6.9 6.9 11,78

Молибден 1*10 ^ 9 9.5 8,8 62.2 62.2 99,8

Серебро -1-1- - 0.136 Г0 7.3 7.3 5,3

Марганец 1*10"3 34 46.9 99,2 46.2 46.2 55,3

Мышьяк -1-1- - - - - -

Сурьма -1-1- - - 0,3 3 3

Вольфрам -1-1- - - 0,32 0.3 0.3 1,07

Олово -1-1- 0.44 0.1-6 0,17 4 4 2,09

Висмут 0.175 0.121 - 0.27 0.27 0,46

Барий -1-1- 10 17.7 15,7 19 19 55,53

Скандий -1-1- - - 0,86 - - 1,1

Титан -1-1- 200 85 457 190 190 599

Литий -1-1- 2.63 1.16 - - -

Бериллий -1-1- 0.15 0.135 - - -

Ниобий -1-1- 0.3 0.45 - 0.82 0.82

Тантал -1-1- - - - - -

Иттрий -1-1- 0.4 0.9 1,4 2.91 2.91 1,73

Иттербий -1-1- - 0.13 0.33 0.33

Цирконий -1-1- 4 20 10,2 3.6 3.6 20,8

Гафний -1-1- - - - - -

Лантан -1-1- - - - - -

Церий -1-1- - 7.6 - - -

Кадмий -1-1- - - - - -

Индий -1-1- - - - - -

Фосфор -1-1- 42.5 40 64 70.1 70.1 116,7

Германий -1-1- 0.05 - - 0.13 0.13 0,32

Галлий -1-1- - 1.38 - 4.0 4.0 1,86

Бор -1-1- - - 2,7 - - 9,75

Таллий -1-1- - - - - -

Стронций -1-1- 1.5 3.7 2,38 - - 20

Ртуть 1x10"6 - - 1,9 - - 1,9

17

7= ¿^-(N-1) (2)

1=1 где: N -количество элементов, участвующих в расчете показателя 2; Кк - коэффициент концентрации:

Таблица 2.

Уровни загрязнения почвы по величине Z

Уровень загрязнения Значение Ъ

Минимальный Ъ< 8

Слабый 8<г<1б

Средний \Ь<Ъ<Ъ2

Высокий 32<г<64

Очень высокий 64<г<128

MaкeимaJ 1 ьный 7>\2%

Как видно из формул 1 и 2 значение показателя Z функционально связано с тремя исходными параметрами: Сх, Сф и N.

Содержание ТМ в почве - Сх зависит от её типа, особенно фракционного состава, способности к самоочищению, величины, характера и продолжительности техногенной нагрузки. Так, природные концентрации ТМ в «легких» разновидностях даже одного и того же типа почв обычно в 1.5-2.0 ниже, чем в глинистых. К тому же, песчанистые почвы обладают и повышенной способностью к самоочищению и до определенного уровня техногенная нагрузка может не проявляться в них повышением концентраций металлов. Тяжелые почвы, особенно обогащенные гумусом, наоборот, имеют повышенные сорбционные свойства и начинают накапливать ТМ даже при невысокой техногенной нагрузке. Кроме того, в результате перемешивания с подстилающими материнскими отложениями и за счет завоза различных почво-грунтов для рекультивации в городах формируется особый тип почв, которым иногда невозможно найти природный аналог. Значительное влияние на содержание ТМ в почве оказывает, кроме величины, характер техногенной нагрузки. В некоторых очагах с интенсивным загрязнением в результате изменения рН и ЕЬ почвы и, как следствие, миграционной способности элементов, на фоне общего увеличения содержания одних из них может отмечаться относительное уменьшение концентраций других. Немаловажное влияние на содержание ТМ в почве оказывает л продолжительность техногенной нагрузки.

18

Для приведения в соответствие загрязненности почвы с величиной техногенной нагрузки требуется не один год.

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

Практикалроведения этих, работ доказала, что наиболее эффективна и надёжна такая идентификация математическим методом. В качестве математической основы для данной цели выбран коэффициент пропорционального сходства «Имбри и Парди» - Яху, вычисляемый по формуле:

Я =-—--(3} П? * 1 ;

V г=1 ¿=1 где Х% и У; - содержание одноименных элементов в пробах X и V.

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

19

При оценке загрязнённости снеговой пыли часто складывается ситуация, когда фоновые концентрации ТМ в пробах, отобранных на удалении 150-200 км. от промышленных объектов, близки к их среднему содержанию по городу [9]. Данный факт объясняется в первую очередь различным фракционным составом пыли - более тонким на удалении от города и с относительно более высоким в этой пыли содержанием ТМ. В первую очередь это касается таких металлов, как РЬ, Си, фоновое содержание которых в снеге обычно на порядок выше такового в почве. Кроме естественной дифференциации металлов в атмосфере, существенное влияние в пределах городской черты на фракционный и минеральный состав пыли, оказывают вторичные аэрозоли, разубоживающие содержание ТМ. Кроме того, сами по себе аномалии, выделенные по содержанию ТМ в снеговом покрове в пределах города, являются обычно более контрастными, а интенсивность их для одних и тех же постоянных источников загрязнения, как правило в 3-4 раза выше почвенных. Поэтому аномалии, выделенные по загрязненности снега и почвы, по расположению в плане и, особенно, интенсивности имеют между собой часто слабую корреляцию.

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

Определённое в л ияниела-величину показателя оказывает и количество элементов - N, имеющих значимые содержания и участвующих в его расчете. Из 36-40 элементов, Iр которые проводится спектральный анализ, лишь 20-22 имеют в лочве устойчивые значимые содержания. Из остальных 8-10 имеют содержания на пороге чувствительности данного метода и обнаруживаются лишь при повышенных концентрациях. Поэтому при расчете показателя Ъ для проб, отобранных в пределах очагов загрязнения и вне их, участвует различное количество элементов, в связи с чем появление некоторых аномалий, -особеннох. невысоким уровнем загрязнения, в большей мере связано с увеличением количества элементов, а не содержания каждого из них.

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

20

Ещё более условным является оценка экологической опасности техногенной нагрузки по значению показателя Ъ в соответствии с ГОСТом 17.4.1.02-83 [13] по шкале, приведенной в таблице 3.

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

Таблица 3.

Категория загрязнения Величина Ъ Изменения показателя здоровья населения

Допустимая Ъ< 16 Наиболее низкий уровень заболеваемости населения

Умеренно опасная 16 < ъ < 32 Увеличение общей заболеваемости

Опасная 32 < Ъ < 128 Увеличение общей заболеваемости детей с хроническими заболеваниями, нарушение функционального состояния

Чрезвычайно опасная Ъ> 128 Увеличение общей заболеваемости детей, опасные нарушение репродуктивной функции женщин (увеличение токсикоза беременных, числа преждевременных родов, мертворожденных, гипотрофии новорожденных)

Более перспективным в этом плане возможно является расчет данных показателей с учетом поправочных да класс опасности коэффициентов. Таковыми могут быть или «комплексный индекс загрязнения» [27], при определении .которого каждый коэффициент концентрации^включенныйз суммарный до-казателъ, возводится вхтепень в зависимости от классаопасности, или « суммарный показатель опасности», рассчитываемый относительно нормативных ПДК, или, если они отсутствуют, относительно ПДК, определенных, к примеру, по методике, предложенной Г.А. Вострокнутовым [12] - «суммарный показатель токсичности».

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

21 ности могут отличаться. В настоящей работе для оценки опасности загрязнения грунтовых вод традиционным методом применялась оценочная шкала, предложенная легендой ЦРГЦ для целей мелкомасштабного картирования, фрагмент которой приведен в таблице 4.

Таблица 4.

Экологическое состояние по загрязнению подземных вод

Подземная гидросфера. Гидрохимическая обстановка

ТО V I 1 ч 11 о!V химические '^пацдтт'гт Т Л ^ТТЛИОТТ^Г т 1 ! I ! \ 1Л. 1у1 11! Iл 1 Ы. ну ^IX !

1 класс 2 класс 3-4 класс Оценка состояния приопасности опасности опасности родной среды 1.0 < 1.0 < 1.0 Допустимое

1 - 1.5 1.0-2.0 1.0 - 5.0 Умеренно-опасное

1.5-2.0 2.0-5.0 5.0- 10 Опасное

2.0-3.0 5.0 - 10 10- 15 Высоко опасное

3 > 10 > 15 Чрезвычайно опасное

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

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

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

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

Сбор информации о детской заболеваемости по городам проводился несколькими организациями, использующими при этом различный методический подход. Так, по г. Россошь кафедрой госпитальной педиатрии Воронежской медицинской академии были собраны данные о заболеваемости по 19 медицинским участкам за 1993 г. по принятой годовой форме отчетности. По г. Ливны городская детская поликлиника представила уже сгруппированные данные о среднем количестве стационарных больных по 13 участкам за 1991 - 1992 гг. По г. Воронеж Центральной опытно - методической экспедицией по 203 детским медицинским участкам были собраны данные о количестве хронических больных в 1998 г. по отдельным видам заболеваний и их группам.

Наиболее полные, пригодные для дальнейшей обработки, данные были получены для Россоши - по шести группам и двум отдельным заболеваниям: органов пищеварения (ОР), дыхания (OD), мочеполовой системы (MS), сердечно-сосудистой системы (KAV), вегетативно-сосудистой дистонии (VSO), нервной системы (NER), заболеванию щитовидной железы (SHT) и пиелонефриту (МРР). В г. Ливны - для 11 из 13 участков - по шести группам заболеваний: инфекционным и паразитарным (IP), органов пищеварения: (ОР), дыхания (OD), мочеполовой системы (MS), эндокринной системы (EYS ) и кроветворных органов (КО). В г. Воронеж, на его правобережной части, для 140 участков по десяти отдельным заболеваниям и их группам: органов пищеварения, всего (ОР), холециститу (ОРХ), органов дыхания, всего (OD), пневмонии (ODP), бронхиальной астме (ODA), астматическому бронхиту (ODB), органов мочеполовой системы, всего (MPS), пиелонефриту (МРР), дискенезии желчно-выводящих путей (DIS), всем заболеваниям, кроме ОРВИ (VSE).

Для дальнейших исследований фактическое количество больных по каждому участку и заболеванию - St приводилось к 1000 чел. по формуле: где: S{ - приведённый уровень заболевания «I», N" - количество детей, закреплённых за участком.

23

Для первичной оценки фактической заболеваемости по городам было проведено их районирование. При этом, по г. Воронеж учитывалось данные только о количестве хронических заболеваний ОР, а по гг. Ливны и Россошь - о всех заболеваниях. В последнем случае соизмеримость частоты нескольких заболеваний, имеющих разные абсолютные значения, достигалась ранговым методом. Суть данного метода заключается в присвоении каждому участку ранга в зависимости от его места в ранжированном для каждого заболевания ряду. Затем для каждого участка определяется сумма рангов всех заболеваний, по значениям которых в изолиниях составляются схемы общей заболеваемости. Так, третьему участку железнодорожной поликлиники, имеющему самую высокую (из девятнадцати) приведенную заболеваемость ОБ в г. Россошь (табл. 6) присваивается 19 ранг. Сумма же рангов десяти заболеваний для этого участка составляет 108, при том, как его среднее для всех участков значение равно 80 (10 х 8). То есть, по г. Россошь для участков, имеющих ранговое значение больше 80, общую заболеваемость следует считать повышенной относительно средней.

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

Расчёт среднего содержания каждого компонента - С, по участкам проводился по формуле: где: /V' - количество проб, участвующих в расчете среднего содержания по каждому участку.

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

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

5)

24

1 /

6)

М j=l где: Я - количество участков (наблюдений), по которым проводится корреляция.

Значение «г» изменяется от +1 до -1 в зависимости от степени синхронности колебания значений сопоставляемых параметров Л^иС,,, если считать, что связь между ними имеет линейный характер типа: Ьа Ск + си ( 7 )

Определение коэффициентов регрессии: Ъи и с у, проводилось «методом наименьших квадратов», по которому:

IV,- 1 1

Ъи =

1=1 N

Ес,г)2

8) N у=1 N

Си^^А-ЬиЪС^/Ы

1 7=1

9)

Вес изменчивости содержания отдельного металла / в изменчивости каждого заболевания I при линейном характере их связи определялся по величине:

N N N ]=\ ]=\ /к 2

N 1 ЦЫ 2 X N ¿=1 ---~2 -(¿С^/М >1

М j=l

1=1

10)

У=1 .7=1

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

И) р ас. ' ** * "

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

25 poro эти колебания асинхронны. В этом случае перегиб параболы, графически отражающей уравнение второй степени, будет находиться в районе значения заболеваемости, который можно принять за местный фоновый уровень, вычислив по формуле: el - 4а,с,,

Su. = и i¡

4 а,;

12)

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

Уровень связи между содержанием металла и заболеваемостью в этом случае определяется уже по значению множественного коэффициента корреляции -К, а вес изменчивости содержания в заболеваемости - по величине Я , которая рассчитывается в общем случае по формуле:

Л ( N а, 1

Ео s~i2 7=' 1

ЛГ !LSjI^C,

N \

Y. N f N

13) v ./=1

7=1 где: п - общее количество параметров, относительно которых проводится множественная регрессия; т и (п-т) - соответственно количество параметров второй и первой степени, относительно которых проводится множественная регрессия.

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

26

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

В настоящей работе при моделировании использовались в основном простейшие линейные уравнения множественной регрессии типа:

V, = ЬиСи + ЪпСп.+ ЬИСц л си ( 14 )

В полиэлементные модели включались только те металлы, которые по данным спектрального анализа имели устойчивые значимые содержания. Это, как правило, 2527 элементов из 36 определяемых - в снеговом покрове, 20-22 - в почве, 5-10 - в воде. После расчета уравнений множественной регрессии со всем комплексом металлов для каждого заболевания по формуле 13 определялся их общий вес - К:. Дальнейшие исследования проводились уже только с теми заболеваниями, изменчивость которых определялась загрязненностью металлов не менее, чем на 65%.

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

ЕГ-=Я2п-Я2,и ( 15)

Какие-либо жесткие критерии оценки значимости удельного веса каждого металла отсутствуют. Все определяется полнотой представленного материала и предъявляемыми к нему требованиями. Поэтому практически для этих целей вместо 1<1 использовался показатель Я'и, рассчитываемый по формуле 16, который кроме веса уравнения учитывает и его надежность - обеспеченность модели количеством опытов (в рассматриваемых случаях - количеством медицинских участков, по которым имеются данные о заболеваемости). я'п = -п2/№), (16)

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

27 ределяемых элементов при оптимизации по этому показателю обычно 8-12 исключаются с уменьшением общего веса всего на 2-4%.

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

17)

1=1

При значениях (}п > 1.8 расхождения результатов интерполяции модели с фактическими значениями, очевидно за счет приближений, допущенных на разных стадиях расчета, как правило, превышают ожидаемые, определяемые значением К2ц . В случаях, когда эти отклонения достигают значительных величин и не удовлетворяют предъявляемые требования, результаты моделирования следует признать отрицательными, а саму модель непригодной для дальнейшей интерполяции.

Для оценки устойчивости модели с учетом ее достоверности часто использовался синтетический показатель ¥ц, рассчитываемый по формуле:

Уп = Си1К'п (18)

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

Следующим шагом исследований является пересчет содержания ТМ в отдельных точках опробования по определенным оптимизированным моделям в частоту той или иной заболеваемости. Когда значимые связи с содержанием ТМ имеют несколько заболеваний, то для сопоставления их частоты, имеющей разные абсолютные значения, они приводятся к фоновым уровням, принцип расчета которых приведен выше. А для каждой точки опробования в таких случаях определяется средний уровень заболеваемости - Я/ по формуле:

28 к V 1 К

19) где: К - количество заболеваний, для которых рассчитывались модели.

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

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

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

29

Значительно более достоверными и простыми являются карты зонирования по ОЭС составленные по данным отражающим постоянное загрязнение. На картах этого типа выделяются только те очаги, которые подтверждается результатами почвенного и снегового опробования. Графически это достигается методом «пересечения» площадей с одним и тем же рассчитанным по загрязнению почвы и снега уровнем заболеваемости.

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

Все этих три варианта карт рассмотрены при зонировании по ОЭС территории г. Воронеж.

Иногда, при анализе природы очагов повышенной заболеваемости является оправданным составление карт, построенных по значению отклонения рассчитанной заболеваемости от фактической - И«*?/, определяемому по формуле:

20) р ас. *фщ>т.

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

30

Заключение диссертации по теме «Геохимия», Иванов, Юрий Владимирович

ЗАКЛЮЧЕНИЕ

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

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

99

В общем случае применение математического моделирования для оценки остроты экологической ситуации при крупномасштабном картировании значительно повышает эффективность информации, полученной в процессе экологических исследований. Этот метод позволяет максимально учитывать специфику техногенного загрязнения в каждом конкретном случае и его проявленность в заболеваемости населения. Благодаря данному методу возможно сопоставление условных уровней суммарного загрязнения разных сред опробования, часто для одних и тех же очагов не совпадающих по интенсивности. Поэтому наиболее эффективно применение математического моделирования при комплексной оценке остроты экологической ситуации по результатам опробования нескольких депонирующих сред. Это, в свою очередь, позволяет более объективно определять тип загрязнения, что немаловажно, как для выявления источников загрязнения, так и для выбора природоохранных мероприятий Резко повышает данный метод эффективность спектрального анализа, обычно применяемого при экологических исследованиях. Широкий спектр определяемых им элементов и информация о фактической заболеваемости позволяют при моделировании оценивать условную токсичность отдельных металлов без учета формы их нахождения. А, если учесть высокую «поточность» этого метода и низкую себестоимость одного эл.ан., то применение его в комплексе с моделированием предполагает особые перспективы.

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

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

100

Список литературы диссертационного исследования кандидат геолого-минералогических наук Иванов, Юрий Владимирович, 1999 год

1. Беспамятное Г. Б., Кротов Ю. Д. Предельно допустимые концентрации химических веществ в окружающей среде. -Л., 1985. -528 с.

2. Богацкий В. В., Витязь В. И., Коллеганов Ю. М. и др. Методические рекомендации по практическому осуществлению пространственно-статистических исследований. -Новосибирск, 1977. -48 с.

3. Бочаров В.Л., Иванов Ю.В. Автотранспорт и экология малых городов // Экологический вестник Черноземья. -Воронеж, 1996. Вып. 2. -С. 38.

4. Бочаров В.Л., Иванов Ю.В. Зонирование городских территорий по остроте экологической ситуации ( на примере г. Россошь Воронежской области) // Вестник Воронежского университета. Серия геологическая. -Воронеж, 1996. Вып. 2. -С. 108-112.

5. Бочаров В.Л., Иванов Ю.В. Эколого-геохимические методы оценки загрязненности малых городов ( на примере г. Ливны Орловской области ) // Вестник Воронежского университета. Серия геологическая. -Воронеж, 1997. Вып. 4. -С. 142-148.

6. Бочаров В.Л., Иванов Ю.В. Экологическое крупномасштабное картирование в городах с применением математического моделирования // Контроль и реабилитация окружающей среды. Материалы международного симпозиума. -Томск, 1998. -С. 122-125.

7. Бочаров В.Л., Иванов Ю.В. Математическое моделирование при оценке экологической опасности техногенного загрязнения // Экологи. Экологическое образование. Нелинейное мышление. -Воронеж, 1998. -С. 63-71.

8. Бочаров В.Л., Иванов Ю.В., Бурляев В. А. Некоторые приемы статистической обработки материалов при экогеохимических исследованиях101

9. Вестник Воронежского университета. Серия геологическая. -Воронеж, 1998. Вып. 6. -С. 168-173.

10. Виноградов А.П. Геохимия редких и рассеяных химических элементов в почвах. -М., 1950. -278 с.

11. Вострокнутов Г. А. Математическая модель предельно допустимых количеств химических элементов // Всероссийская научно-практическая конференция «Геологическое картографирование». Тезисы докладов. Ч. И. -М.: Изд-во Геоинформмарк., 1998. -С. 153-155.

12. ГОСТ 17.4.1.02-83. Охрана природы. Почвы. Классификация химических веществ для контроля загрязнения. Госстандарт. -М., 1983.

13. Гольцман Ф.М. Статистические модели интерпретации. -М., 1971. -327 с.

14. Григорян С. В. Первичные геохимические ореолы при поисках и разведке рудных месторождений. -М., 1987. -408 с.

15. Девис Дж. С. Статистический анализ данных в геологии. Кн. 1. -М., 1990.- 319 с.

16. Девис Дж. С. Статистический анализ данных в геологии. Кн. 2. -М., 1990. 427 с.102

17. Зверев В. П. Роль подземных вод в миграции химических элементов. -М., 1982. -184 с.

18. Иванов В.В. Экологическая геохимия элементов. Справочник. Кн. 1 .М., 1994. -303 с.

19. Иванов Ю.В. Зонирование городских территорий по остроте экологической ситуации // Проблемы интеграции экологической и хозяйственной политики в Черноземном Центре России: Тез. докл. 1-ой регион, научно-практ. Конф. 4.1. -Мичуринск, 1995. -С. 43-45.

20. Иванов Ю.В. Геохимический метод и математическое моделирование при оценке загрязненности атмосферы городов // Геоэкологические проблемы устойчивого развития городской среды. -Воронеж, 1996. -С. 190-192.

21. Канцерогенные вещества // Справочник. Материалы международного агентства по изучению рака. -М., 1987. -334 с.103

22. Ковалевский В. С. Исследования режима подземных вод в связи с их эксплуатацией. -М., 1986. -200с.

23. Копчак В.В. Разработка методики составления комплексной территориальной схемы охраны окружающей среды // НТО по теме А-35. -Л., 1985.

24. Крамбейн У., Грейбилл Ф. Статистические модели в геологии. -М.,1969. -397 с.

25. Ким Дж. -О., Мьюллер Ч. У., Клекка У. Р. и др. Факторный, дискри-минантный и кластерный анализ. -М.,1989. -215 с.

26. Кулиш С. А., Науменко К. Д., Сивых В. Б. и др. Математические методы при планировании и управлении горным производством. -М., 1978. -320 с.

27. Купер Дж., Макгиллем К. Вероятностные методы анализа сигналов и систем. -М., 1989. -376 с.

28. Математические методы в поисково-разведочной практике. -Иркутск,1970. -352 с.

29. Методические рекомендации по геохимической оценке загрязнения территории городов химическими элементами / Под ред. В. А. Ревича, Ю. Е.Саета. -М.: Изд-во ИМГРЭ, 1982. -112с.

30. Методические указания по оценке загрязнения почвы химическими веществами // Минздрав СССР. Главное санитарно-эпидемиологическое управление. -М., 1987.

31. Перельман А. И. Геохимия. -М., 1989. -528 с.

32. Питулько В. М., Резников И. Н., Ульянов Н. К. Литохимические метод съемки и поисков // Методическое пособие по геологической съемке масштаба 1:50000. Вып. 15. -Л., 1985. -199 с.

33. Применение математических методов при поисках и разведке месторождений твердых полезных ископаемых // Методическое руководство. -М., 1987. -262 с.104

34. Протасова Н.А. Редкие и рассеянные элементы в почвах Среднерусской возвышенности. -М., 1985.

35. Протасова Н. А., Щербаков А. П., Копаева М. Т. Редкие и рассеянные элементы в почвах Центрального Черноземья. -Воронеж, 1992. -166 с.

36. Смирнова А. Я., Умнякова Л. В., Гольдберг В. М. Грунтовые воды и их естественная защищенность от загрязнения на территории Воронежской области. -Воронеж, 1986. -108 с.

37. Принципы и методика геохимических исследований при прогнозировании и поисках рудных месторождений // Методические рекоменда105ции. Под редакцией А. А. Смыслова, В. А. Рудника, Н. М. Динкова и др. -Л., 1979. -348 с.

38. Справочник по геохимическим поискам полезных ископаемых // Под редакцией А. П. Соловова. -М., 1990. -336 с.

39. Тшценко Н. Ф. Охрана атмосферного воздуха. -М., 1991. -362 с.1. Фондовая

40. Комплексная геохимическая оценка г. Воронеж с прилегающей территорией: Произв. отчет / Центральная опытно-методической экспедиция.; Руководитель Ю. В. Иванов. Инв. № 139. -Воронеж, 1989. - 129 с. - (Фонды ЦОМЭ).

41. Литомониторинг на промплощадке Россошанского химзавода им. 60 -летня СССР: Отчет о НИР / Воронеж, гос. ун-т; Руководитель В. Ф. Черняев. Воронеж, 1990. -150с. -(Фонды Комитета экологии и природных ресурсов г. Россошь).

42. Экологические исследования г. Россошь с прилегающей территорией: Произв. отчет / КПП "Геокомплекс"; Руководитель В. А. Чернов. -Инв. № 190. -Воронеж, 1991. -153 с. (Фонды КПП "Геокомплекс).