1. ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
1.2. Связь работы с научными программами, планами, темами
2. АНАЛИЗ РАНЕЕ ПРОВЕДЕННЫХ ИССЛЕДОВАНИЙ
3.1. Методика геохимических исследований
3.2. Методика геофизических исследований
3.1. Результаты геохимических исследований
3.2. Результаты геофизических исследований
СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ
Тема работы актуальна в связи с тем, что зона сочленения Донбасса с Приазовьем представляет собой глубинный разлом, по которому идет сток глубинного вещества. Разрывные нарушения, большое разнообразие пород, многочисленные рудопроявления и небольшие мощности покровных отложений способствуют накоплению многих тяжелых металлов в почвах. Разломные зоны, перекрытые четвертичными отложениями, имеют различную активность и проницаемость и по-разному влияют на накопление токсичных веществ в почвах. На состав почв влият также состав коренных почвообразующих пород, многие из которых содержат элементы в опасных для человека концентрациях. В связи с этим актуальной проблемой является установление таких опасных участков комплексом новых геофизических (структурно-геодинамических) и традиционных геохимических методов.
1.2. Связь работы с научными программами, планами, темами
Работа связана с научно-исследовательской работой кафедры ПИ и ЭГ Н8-07 «Исследование геологических особенностей и строения зон экологического риска Донецкой области и разработка методики их прогнозирования», с темой «Сохранение окружающей среды и устойчивое развитие» основных научных направлений и важнейших проблем фундаментальных исследований в области природных, технических и гуманитарных наук на 2009-2013 г.
Целью исследований является установление геофизических и геохимических особенностей пород и грунтов зоны сочленения Донбасса с Приазовьем, которые позволяют выявлять экологически опасные участки земель сельскохозяйственного назначения.
1. Анализ современного состояния проблемы по литературным данным и фондовым материалам.
2. Проведение геофизических полевых исследований методом СГДК-А и установления особенностей проявления разломных зон в электромагнитном поле, оценка их проницаемости для токсичных элементов.
3. Установление характера влияния состава коренных пород на концентрацию токсичных элементов в почвах по геохимическим данным.
4. Построение карт геохимиеских аномалий элементов, превышающих ПДК и региональный фон.
5. Разработка рекомендаций по прогнозированию экологически опасных участков земель сельскохозяйственного назначения.
Породы и почвы Комсомольского, Камышевахского и Петровского участков зоны сочленения Донбасса с Приазовьем
Геофизические и геохимические особенности горных пород и почв зоны сочленения Донбасса с Приазовьем.
Геофизические методы (азимутальный способ структурно-геодинамического картирования); литогеохимические методы, методы математической статистики обработки данных.
Установлены особенности влияния геологических, тектонических и геодинамических факторов на формирование экологически опасных концентраций химических элементов в почвах зоны сочленения Донбасса с Приазовьем.
Проведенные исследования позволяют выделить экологически опасные участки, связанные с природными аномалиями элементов геологической и тектонической природы в почвах. Такой подход позволяет дифференцировать земли сельскохозяйственного назначения по их ценности и экологической безопасности.
Проведение полевых геофизических исследований и выявление геофизических аномалий, установление характера зависимости состава почв от типа подстилающей породы; выявление особенностей накопления химических элементов в почвах над зонами разломов; построение эколого-геохимических карт и установление участков загрязнения.
Результаты исследований докладывались на конференции «География, геоэкология, биология, геология: опыт научных исследований в контексте международного сотрудничества», 2010 г., ДНУ им. Олеся Гончара, Днепропетровск; на конференции «География, геоэкология, биология, геология: опыт научных исследований в контексте международного сотрудничества», 2011 г., ДНУ им. Олеся Гончара, Днепропетровск; на конференции «Современные тенденции научной парадигмы географического образования Украины», 2010 г., ДИСО, Донецк; на конференции «Охрана окружающей среды и рациональное использование природных ресурсов», 2010 г., ДонНТУ, Донецк. По теме работы имеется 4 публикации.
Комплекс методов структурно-геодинамического картирования (СГДК) был разработан в 70-е годы прошлого столетия для возможности оценки современной активности тектонических структур. Основателями данных методов являются Ю. С. Рябоштан и Е. П. Тахтамиров. Один из наиболее экспрессных методов этого комплекса – электромагнитный азимутальный [1, с. 97]. С целью изучения энергомассопереноса исследования данным методом были проведены В. И. Алехиным [2, 3]. Автором установлено, что вдоль тектонически активных структур происходит сток энергии и глубинного вещества, отмечается формирование месторождений и рудопроявлений эндогенного генезиса различного возраста.
В пределах зоны сочленения Донбасса с Приазовьем геофизические исследования методом СГДК-А проведены в недостаточном объеме. На исследуемых же участках зоны сочленения Донбасса с Приазовьем (Комсомольский, Камышевахский, Петровский) исследования данного рода ранее не проводились и будут проведены впервые.
Изучению геологии и рудоносности зоны сочленения Донбасса с Приазовьем (Волновахской зоны разломов) посвящены исследования И. С. Усенко, Л. Г. Бернадской, Г. А. Макухиной, Д. Е. Айзенберга, Ю. Б. Устиновского, Б. С. Панова и многих других [4]. Обобщенный итог многолетних исследований данной площади приведен в монографии С. В. Нечаева «Минерализация Волноваской зоны разломов». Исследователями этой зоны установлено, что большинство выявленных рудных комплексов тяготеют к разрывным структурам участка. К этим комплексам относятся жилы кобальтистского пирита, кальцитовые, доломитовые жилы, зоны окварцевания, точки ртутной и флюоритовой минерализации и других металлов.
Накоплению руд свинца на исследуемой площади посвящена работа Б. С. Панова и В. И. Алехина «Структурно-тектонические особенности и минеральный состав нового полиметаллического проявления зоны Южно-Донбасского глубинного разлома» [5]. Авторами установлено, что формирование полиметаллического оруденения приурочено к разрывным нарушениям.
В. А. Корчемагиным, В. А. Дудником, Б. С. Пановым, В. И. Алехиным установлена связь полей деформаций и напряжений с рудной минерализацией [6].
Проблеме влияния разломов на состав почв посвящена работа М. В. Пристинской, П. В. Койнаша, В. И. Алехина [7], статья В. И. Алехина, Б. С. Панова, П. В. Койнаша [8]. Авторами установлено, что некоторые геохимические аномалии в почвах контролируются разрывными структурами. При этом в почвах площади исследований обнаружены экологически опасные аномалии многих химических элементов: цинк, мышьяк, фосфор, кобальт, никель и др.
Результаты исследований почв в зоне влияния разрывных структур Приазовья представлены также в работе Е. М. Шеремета, К. В. Волковой, А. А.Омельченко. Проведенные исследования позволили выявить аномальные концентрации в почвах токсичных элементов, расположенных в зоне влияния Малоянисольской и Кальмиусской зоны разломов [9].
3.1 Методика геохимических исследований
В работе использованы данные геохимических съемок по вторичным ореолам рассеивания элементов в почвах, выполненные Приазовской геолого-съемочной партией. Анализы проб почв выполнены спектральным полуколичественным методом.
Методика геохимических исследований включает в себя, в первую очередь, подготовку и формирование выборок элементов по каждому участку и в зависимости от типа подстилающей породы. С помощью программного пакета SPSS производится расчет статистических характеристик по каждому участку и по каждому типу горных пород. Затем полученные характеристики сравниваются с целью установления особенностей накопления элементов в почвах над разными типами подстилающих пород. Нормирование содержаний элементов по ПДК и региональному фону производится по следующей формуле:
где n – коэффициент превышения регионального фона или ПДК; С – содержание элемента в конкретной пробе; Сф (ПДК) – значение регионального фона или ПДК для исследуемого элемента.
В случае если n=2 и больше, содержание элемента можно считать экологически опасным.
Картирование участков загрязнения осуществлялось с использованием компьютерной программы SURFER. Моноэлементные карты в изолиниях строились для элементов, концентрации которых на участке испытывают наибольшие колебания и превышают ПДК или 2-х кратный уровень регионального фона [10].
3.2 Методика геофизических исследований
Для проведения полевых геофизических исследований выбран метод структурно-геодинамического картирования азимутальным способом (СГДК-А). Выбор данного метода обоснован тем, что он предназначен именно для выявления и картирования геодинамических зон. Метод прошел проверку в различных странах, является легким в освоении и применении, оперативен и экономически выгоден, а также обладает высокой точностью и достоверностью.
Способ СГДК-А основан на явлении азимутальной неоднородности электропроводности поверхностного слоя покровных отложений в связи с геодинамическими процессами в недрах. Для реализации способа СГДК-А разработан прибор – электронный фиксатор аномалий (ЭФА). Прибор ЭФА позволяет оперативно выявлять геодинамические зоны по изменению анизотропии электропроводности в покровных отложениях до глубины 2,5-3 м [11].
На каждой точке (пикете) профиля установка ЭФА ориентировалась с помощью компаса на север. Затем в этом положении измерялась электропроводность, и брус перемещался в горизонтальной плоскости по часовой стрелке для снятия замеров в других направлениях. Измерения проводились с угловым шагом в 30 градусов. Каждому фиксированному положению приемника присваивался порядковый номер (код) от 0 до 12. Исходная позиция имела номер 0. Последний замер брался в позиции 12. Последний замер дублировал измерения в исходной позиции и являлся контрольным. В случае отличия первого и последнего замера более чем на 4 единицы, измерения на точке наблюдения повторялись. Такой прием обеспечивал высокое качество наблюдений [11].
Снятые замеры отражают электропроводность грунтов в различных направлениях. Их анализ позволяет установить направление с максимальной электропроводностью в пределах каждого из 4 квадрантов (секторов) круга на каждом полевом пикете съемки.
Для выявлений аномальной анизотропии электропроводности грунтов использовались три показателя: К1 – степень устойчивости ориентировок максимальной электропроводности по профилю; К2 – степень отличия ориентировок максимальной электропроводности на пикете от глобального фона; К3 – степень отличия ориентировок максимальной электропроводности на пикете от фона участка съемки [11]. По этим трем показателям строились графики. Аномальные участки интерпретировались как геодинамические зоны, связанные с разрывными нарушениями.
4.1 Результаты геохимических исследований
С целью определения зависимости формирования состава почв от типа подстилающей горной породы, были исследованы почвы Комсомольского, Камышевахского и Петровского участков зоны сочленения Донбасса с Приазовским кристаллическим массивом (рис 4.1).
Рисунок 4.1 – Схема расположения исследуемых участков на территории зоны сочленения Донбасса с Приазовьем.
1 – граниты; 2 – базальты; 3 – трахиты; 4 – андезиты; 5 – известняки доломитизированные; 6 – известняки; 7 – песчано-глинистые сланцы; 8 – ультрамафиты; 9 – разрывные нарушения; 10 – предполагаемые разрывные нарушения; 11 – Комсомольский участок; 12 – Камышевахский участок; 13 – Петровский участок.
Всего по Комсомольскому участку было изучено 218 проб почв, отобранным по 9 профилям (рис. 4.2); по Камышевахскому – 1120 проб, отобранным по 52 профилям; по Петровскому – 805 проб почв, отобранным по 26 профилям.
1 – граниты и мигматиты протерозоя; 2 – базальты девона; 3 – трахиты пермского возраста; 4 – точки отбора проб; 5 – номера профиля; 6 – разрывные нарушения.
В пробах были исследованы содержания следующих элементов – P, Pb, Ti, W, V, Mn, Ni, Cr, Ge, Co, Bi, Be, Ba, Nb, Mo, Sn, Li, Cu, Zr, Yb, Y, La, Zn, Ag, Sr. С использованием геологической карты пробы были разделены на группы по соответствующим типам подстилающих пород исследуемых участков (граниты, трахиты, базальты, андезиты). Для каждой группы проб с использованием программного пакета SPSS были рассчитаны описательные средние и максимальные значения.
Кроме того, содержания элементов в пробах были сопоставлены с ПДК и региональным фоном с целью оценки экологической безопасности почв исследуемой территории.
Результаты статистического анализа проб почв Комсомольского участка выявили явные изменения характеристик таких элементов как P, Ti, Mn, Ni, Cr, Sr в зависимости от типа подстилающей породы (табл 4.1). Так, среднее содержание P в почвах меньше над гранитами (60*10-3 %) в сравнении с трахитами (66*10-3 %) и базальтами (66*10-3 %).
Таблица 4.1. Статистические характеристики элементов в почвах участка Комсомольского над разными типами подстилающих горных пород (значения содержаний в n*10-3 %).
| Хим. эл-т | Граниты | Базальты | Трахиты | ||||||
| Средн. знач. | Мин. знач. | Макс. знач. | Средн. знач. | Мин. знач. | Макс. знач. | Средн. знач. | Мин. знач. | Макс. знач. | |
| P | 64,05 | 25 | 100 | 66 | 50 | 100 | 66 | 50 | 70 |
| Ti | 505,4 | 300 | 1000 | 540 | 300 | 700 | 540 | 300 | 700 |
| Mn | 69,9 | 30 | 150 | 70 | 70 | 70 | 58 | 30 | 70 |
| Ni | 5,3 | 1,5 | 20 | 6,5 | 5 | 15 | 4,6 | 3 | 5 |
| Cr | 8,3 | 3 | 30 | 8,2 | 7 | 10 | 10,4 | 7 | 15 |
| Sr | 7,7 | 7 | 15 | 13,6 | 7 | 50 | 8,2 | 7 | 10 |
Среднее содержание Ti также меньше в почвах на гранитах (500*10-3 %), чем на базальтах и трахитах. На базальтах и трахитах это значение составляет 550*10-3 %. Значения содержаний Mn меньше в почвах на трахитах (58*10-3 %), чем на базальтах и гранитах, где они составляют в среднем 70*10-3 %. Максимальные значения Ni приурочены к базальтам, где они в среднем равны 7*10-3 %. В трахитах его среднее значение равно 4,5*10-3 %, а в гранитах – 5*10-3 %.
Средние содержания Cr в почвах достигают больших значений на трахитах (10,5*10-3 %), базальты и граниты характеризуются меньшим средним содержанием данного элемента (8*10-3 %).
Среднее содержание Sr в почвах на базальтах практически вдвое превышает содержание данного элемента в почвах на трахитах и гранитах.
По результатам сопоставления содержаний элементов с региональным фоном и ПДК можно сделать следующие выводы. Содержания Ti незначительно, но повсеместно превышает региональный фон. Однако на 10 точке 18 профиля наблюдается двойное его превышение, что может быть связано с техногенным загрязнением почв. В отдельных точках наблюдается тройное превышение фона таким элементом как Ni. 8 и 11 точка 21 профиля характеризуются двойным превышением ПДК данным элементом. Двойное превышение фона в точках 19, 21 профиля 17 обусловлено влиянием Васильевского разлома, отделяющего базальты от гранитов. Содержание Со повсеместно превышает фоновые концентрации. В отдельных точках наблюдается тройное и двойное превышение ПДК таким элементом как Cr, что также может быть связано с разрывными нарушениями и базальтами. Превышения регионального фона отмечены и для Сu, Zn, Sr, Ba.
Результаты статистического анализа проб почв Камышевахского участка выявили явные изменения характеристик таких элементов как P, Ti, Mn, Sr в зависимости от типа подстилающей породы (табл. 4.2).
Таблица 4.2. Статистические характеристики элементов в почвах участка Камышевахского над разными типами подстилающих горных пород (значения содержаний в n*10-3 %).
| Хим. эл-т | Граниты | Базальты | Трахиты | Андезиты | ||||||||
| Средн. знач. | Мин. знач. | Макс. знач. | Средн. знач. | Мин. знач. | Макс. знач. | Средн. знач. | Мин. знач. | Макс. знач. | Средн. знач. | Мин. знач. | Макс. знач. | |
| P | 67 | 50 | 150 | 75,2 | 50 | 200 | 76,4 | 50 | 150 | 72,9 | 50 | 150 |
| Ti | 515 | 300 | 1000 | 618,9 | 300 | 1000 | 584,5 | 500 | 1000 | 654,4 | 500 | 1500 |
| Mn | 69,5 | 30 | 100 | 71,8 | 30 | 150 | 70,7 | 50 | 100 | 73,1 | 70 | 20 |
| Sr | 7,6 | 2,5 | 10 | 7,5 | 2,5 | 20 | 8,1 | 2,5 | 30 | 9,2 | 7 | 20 |
Минимальные значения P приурочены к почвам над гранитами, где его среднее содержание составляет 67*10-3 %, в то время как над андезитами, базальтами и трахитами оно соответственно равно – 72,9*10-3 %, 75,2*10-3 %, 74,6*10-3 %.
Максимальные значения Ti наблюдаются в почвах над андезитами (654*10-3 %), а минимальные над гранитами – 515*10-3%. Содержания Sr достигают своих максимальных значений в почвах, приуроченных к андезитам (9,2*10-3 %), а вот над гранитами и базальтами его средние значения содержаний падают до 7,6*10-3 % и 7,5*10-3 % соответственно. Значения содержаний Mn меньше в почвах на трахитах (70,7*10-3 %) и гранитах (69,5*10-3 %), чем на базальтах и андезитах, где они составляют 71,9*10-3 % и 73,1*10-3 %.
Сопоставление содержаний элементов с региональным фоном и ПДК показало следующие результаты. Содержание в почвах Ti превышают региональный фон на всем участке, в отдельных пробах наблюдается двойное превышение. 1-7 проба 2 профиля характеризуется двойным превышением ПДК таким элементом как V, что может быть обусловлено техногенным загрязнением почв. Двойное превышения ПДК указанным элементом наблюдается и в пробе №1 51 профиля, что может быть связано с переходом почв на андезиты. Кстати, в этой же пробе P превышает фон в 2 раза. Cu и Ba также превышают региональный фон незначительно, но повсеместно.
Результаты статистического анализа проб почв Петровского участка выявили явные изменения характеристик таких элементов как P, V, Mn, Sr в зависимости от типа подстилающей породы (табл. 4.3).
Таблица 4.3. Статистические характеристики элементов в почвах Петровского участка над разными типами подстилающих горных пород (значения содержаний в n*10-3 %).
| Хим. эл-т | Граниты | Базальты | Андезиты | ||||||
| Средн. знач. | Мин. знач. | Макс. знач. | Средн. знач. | Мин. знач. | Макс. знач. | Средн. знач. | Мин. знач. | Макс. знач. | |
| P | 71.4 | 50 | 150 | 76.1 | 50 | 300 | 71.1 | 70 | 100 |
| V | 7.7 | 1.5 | 15 | 8.2 | 1.5 | 20 | 10 | 10 | 10 |
| Mn | 69,8 | 20 | 100 | 71.8 | 50 | 100 | 73.1 | 70 | 70 |
| Sr | 8.6 | 7 | 15 | 8.75 | 2.5 | 30 | 9.5 | 7 | 20 |
Максимальные значения содержаний Р характерны для почв над базальтами (76,1 *10-3 %), в сравнении с гранитами и андезитами, где они составляют 71,4*10-3 % и 71,1*10-3 % соответственно. V, Mn и Sr достигают максимальных значений содержаний в почвах над андезитами (10*10-3 %, 73,1*10-3 % и 9,5*10-3 %).
Результаты сопоставления содержаний элементов с ПДК и региональным фоном следующие. Повсеместно наблюдается превышение фона для P. На точках 20-27 19-го профиля, точках 22-31 20-го профиля, точках 5-9 25-го профиля отмечены двойным превышением фона, что может быть обусловлено влиянием базальтов на формирование аномалий данного элемента. Ti превышает региональный фон в 2 раза практически на всей территории участка (фон для Ti 482*10-3 %, а содержание в пробах повсеместно – 1000*10-3 %). На большей части участка отмечается двойное превышение фона таким элементом как Ва.
Таким образом, для всех исследуемых участков зоны сочленения Донбасса с Приазовьем наблюдается значительное превышение регионального фона. Учитывая то, что данные элементы относятся к первому и ко второму классам опасности, и на данном участке осуществляется сельскохозяйственная деятельность, их концентрации представляют экологическую опасность.
По результатам анализа состава почв и его сопоставления с типом подстилающих пород можно сделать следующие выводы. Для таких пород как трахиты, характерны большие содержания Р, Ti, Cr. Для базальтов характерны повышенные концентрации в почвах Р, Ti, Mn, Sr, Ni, для андезитов – Ti, Sr, Mn, V.
4.1 Результаты геофизических исследований
Съемка методом СГДК-А была проведена на участке Комсомольском. Было пройдено два профиля: профиль № 1 с шагом наблюдения 5 м составил 50 пикетов; профиль № 2 с шагом наблюдения 10 м – 84 пикета.
Результаты обработки и интерпретации данных геофизической съемки выявили на профиле 1 явную аномалию в районе 25–32 пикетов, отражающую положение Васильевского разлома (рис. 4.3).
Рисунок 4.3 – Карта распределения Co в почвах на участке «Комсомольский».
1 – Васильевский разлом; 2 – геодинамические зоны по данным СГДК-А; 3 – профиль СГДК-А; 4 – изолинии концентраций кобальта; 5 – точки отбора проб почв; 6 – контур участка геохимической съемки
Анализ содержаний элементов в пределах зоны влияния выявленного разлома показал следующие результаты. Значительное превышение ПДК и регионального фона отмечено для Co (в 2 раза), Zn (в 3 раза), Ni (в 3 раза), Mo (в 2 раза), Be (в 1,5 раза). Графики изменения содержаний элементов в почвах в пределах 17 профиля геохимической съемки, где геофизическим методом выявлен разлом, представлены на рисунка 4.4, 4.5. Зона разлома (Васильевский взброс) проявлена на данных рисунках в пределах пикетов 18-20 геохимического профиля № 17. Этот участок соответствует геофизической аномалии СГДК-А на профиле № 1 (пикеты 25-32).
Рисунок 4.4 – Графики содержаний элементов по профилю № 17 участка Комсомольского.
1 – Co; 2 – Zn; 3 – Ni.
Рисунок 4.5 – Графики содержаний элементов по профилю № 17 участка Комсомольского.
1 – Mo; 2 – Be.
Положение геодинамической зоны в плане, которая сопровождается аномалиями многих элементов в почвах, показано на рисунке 4.6.
Рисунок 4.6 – Графики показателей СГДК-А по профилю № 1 участка Комсомольского.
1 – К1; 2 – К2; 3 – К3.
Таким образом, проведенные исследования показали тенденцию накопления токсических элементов в пределах зон влияния разломов и трещиноватости, а также влияние на состав почв особенностей состава горных пород. Наибольшие влияние на состав почв по большинству элементов оказывают базальты и трахиты. Эти факторы совместно с техногенными определяют возникновение экологически опасной ситуации на исследуемых участках. Эти же факторы следует учитывать при оценке стоимости земель сельскохозяйственного назначения.
1. Алехин В. И., Аноприенко А. Я., Анциферов А. В., Купенко В. И., Панов Б. С., Приходько С. Ю. Геодинамическое картирование: методы и апаратура. – Донецк: ДонНТУ, Технопарк ДонГТУ «Унитех», 2007. – с. 144.
2. Алехин В. И. Энергомассоперенос в геодинамических структурах локального уровня и его практическое значение // Проблемы общей и региональной геологии. – Новосибирск: Наука, 1971. – с. 33-36.
3. Алехин В. И. Проницаемость и неотектоническая активность разломов Приазовского блока Украинского щита в связи с оценкой их рудоносности // Наукові праці Донецького державного технічного університету. Серія: гірничо-геологічна. Випуск 45. – Донецьк, ДонДТУ, 2002. – с. 38-44.
4. Нечаев С.В. Минерализация Волновахской зоны разломов. – Киев: Наукова думка, 1970. – с. 172.
5. Алехин В. И., Панов Б. С. Структурно-тектонические особенности и минеральный состав нового полиметаллического проявления зоны Южно-Донбасского глубинного разлома // Науковий вісник національного гірничого університету. – № 6. – Днепропетровск, 2005. – с. 12-15.
6. Корчемагин В. А., Дудник В. А., Панов Б. С., Алехин В. И. О связи полей тектонических деформаций и напряжений с рудоносностью в Донбассе // Геофизический журнал. – № 1, т.27. – 2005. – с. 97-109.
7. Пристинская М. В., Койнаш П. В., Алехин В. И. Влияние разломов на формирование экологически опасных аномалий химических элементов в почвах // Охорона навколишнього середовища та раціональне використання природних ресурсів. Зб. доп. 2 міжнар. наук. конф. аспірантів і студентів. – Том 2. – Донецьк, 2003. – С. 117-118.
8. Алехин В. И., Панов Б. С., Койнаш П. В., Коренев В.В. Геохимические особенности почв на карбонатных породах в зоне влияния Южно-Донбасского глубинного разлома // Наук. Праці ДонНТУ. Серія гірничо-геологічна. Вип. 55. – 2003. – с. 120-125.
9. Шеремет Е. М., Волкова К. В., Омельченко А. А. Эколого-геохимические исследования Приазовья с учётом влияния геодинамически активных разломов // Проблеми гірського тиску. Випуск № 12. – Донецьк, ДонНТУ, 2004.–с. 43-53.
10. Методические указания к выполнению курсового проекта по дисциплине «Эколого-геохимическая съемка». Составил: В. И. Алехин. – Донецк, ДонНТУ, 2006. – с. 22.
11. Алехин В. И., Ребенко Я. В., Шелест К. В. Экологически опасные геодинамические зоны и геофизический метод их выявления // Сучасні тенденції наукової парадигми географічної освіти України/ Зб. наукових статей VI міжвуз. науково-практичн. конф. – Донецьк: ДІСО, 2010. – с. 94-97.