Атрибуты сейсмических волновых полей и их использование при решении задач инженерной геологии тема диссертации и автореферата по ВАК РФ 25.00.10, кандидат геолого-минералогических наук Крылаткова, Надежда Анатольевна

• Крылаткова, Надежда Анатольевна
• кандидат геолого-минералогических наук
• 2008

Оглавление диссертации кандидат геолого-минералогических наук Крылаткова, Надежда Анатольевна

Глава 1.

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

Состояние и пути совершенствования методов наземной 10 инженерной сейсморазведки

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

Многоволновая инженерная сейсморазведка - ключ к 16 успешному решению сложных инженерно-геологических задач Роль количественных характеристик сейсмической записи при 23 изучении разреза с участками пониженной прочности

Глава 2 Получение и использование сейсмических атрибутов при 32 изучении инженерно-геологической среды 2.1 Особенности сейсмических волновых полей для некоторых 32 типов инженерно-геологических моделей среды 2.2. Атрибуты сейсмических волновых полей и способы их 45 определения

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

Глава 3 Примеры применения атрибутных геосейсмических разрезов для 70 изучения ослабленных зон в геологических средах 3.1 Выявление закарстованных зон на трассе железной дороги 71

Свердловск-Богданович 3.2. Выявление подземных горных выработок на площадках 82 промышленного и гражданского строительства в г. Березовский 3.3 Решения задач мониторинга состояния грунтовых оснований в г. 87 Каменск-Уральском

Заключение 96

Список литературных источников 97

ВВЕДЕНИЕ

Актуальность проблемы

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

Сейсморазведочные работы для изучения геологического разреза при инженерно-геологических изысканиях, как правило, выполняются на тех участках, которые классифицируются как сложные, при этом перед ними ставятся задачи по локализации и детальному изучению неоднородностей различного генезиса. Чаще всего для строительных целей или мониторинга оснований уже существующих сооружений исследователей интересуют пространственное распределение прочностных показателей, и, в частности, положение участков пониженной прочности. Характерными чертами этих участков является повышенный процент содержания пустот в горной породе по отношению к вмещающим породам и анизотропия физических свойств. Снижению прочности горных пород способствуют различные инженерно-геологические процессы [26], такие как разуплотнение горных пород, выветривание, суффозия, карстообразование и т.д. Эти процессы во многом формируют сложный характер распределения упругих свойств в геологической среде.

Изучение сложно построенных сред с помощью инженерной сейсморазведки на основе традиционно применяемых кинематических методов интерпретации сталкивается с принципиальными трудностями, связанными с погрешностями определения кинематических характеристик волн, неоднозначностью решения обратных задач и т.д. При применении методов сейсморазведки, которые ориентированы на использование или однотипных по поляризации волн (продольных, поперечных, поверхностных релеевских волн и волн Лява), или волн с определенным характером их распространения (отражение, преломление, дифракция и рассеяние), получаемые результаты отображают структурно-вещественное строение изучаемой среды зачастую недостаточно полно и однозначно. Конечно, основным способом преодоления этих недостатков является комплексирование геофизических методов [53]. Однако важной особенностью сейсморазведки является то, что она сама может выступать как комплекс методов, характеризующих объект исследования посредством привлечения информации о динамике и кинематике прохождения через изучаемую среду различных типов сейсмических волн. Перспективы одновременного использования различных сейсмических волн и, дополнительно, их динамических параметров достаточно очевидны и были поняты геофизиками давно. Однако эффективная реализация этих возможностей стала возможна только при применении современных цифровых сейсмостанций с существенно более широким, чем ранее, динамическим и частотным диапазоном регистрируемых сейсмических волн, а также с использованием компьютерных технологий, в значительной степени облегчивших анализ волновой картины [1].

В последние десятилетия в сейсморазведке МОГТ, ориентированной главным образом на поиски и разведку месторождений углеводородов, помимо кинематических параметров (времен прихода, скоростей распространения сейсмических волн, амплитудных аномалий) при прогнозировании геологического разреза (ПГР) и фильтрационно-емкостных свойств (ФЕС) горных пород успешно используются дополнительные характеристики волнового поля [85]. Эти характеристики часто формируют как функции базовых динамических или кинематических параметров сейсмических записей, которые не являются в строгом смысле параметрами волнового поля. Они занимают среди обычно используемых характеристик некоторое особое положение и получили обобщенное название — «сейсмические атрибуты». Количество рассматриваемых при интерпретации данных МОГТ сейсмических атрибутов составляет сегодня несколько десятков и продолжает увеличиваться. Получаемые поля значений сейсмических атрибутов позволяют создавать новые виды изображений геологической среды, которые облегчают процесс решения структурных задач, задач ПГР и прогноза ФЕС. Широкое применение атрибутов в нефтяной сейсморазведке открыло совершенно уникальную возможность создания тонкослоистых детальных моделей геологического пространства. Применение аналогичного подхода в инженерной сейсморазведке методом ОГТ также открывает новые возможности видения и оценки строения и состояния геологической среды, хотя следует отметить, что сам термин атрибуты здесь обычно не рассматривается. Так в работах Санфирова И.А. и др. [72,73] при малоглубинных сейсмических исследованиях МОГТ с целью определения ослабленных зон в геологическом разрезе для расчета комплексного параметра-индикатора рассматриваются четыре величины: эффективные скорости, максимум амплитудно-частотного спектра, отношение сигнал/помеха, средние арифметические абсолютных значений амплитуд. К сожалению, в основном методе инженерной сейсморазведки - методе преломленных волн (МПВ), сейсмические атрибуты не нашли еще широкого применения, хотя отдельные элементы подобного анализа (изучение затухания амплитуд и частот [51,56], характера полей времен и скоростных законов [17]) известны.

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

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

1) анализ возможности применения сейсмических атрибутов в инженерной сейсморазведке и оценка эффективности использования этих атрибутов;

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

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

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

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

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

При написании работы применены следующие методы исследования:

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

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

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

- анализ информативности атрибутов сейсмических записей;

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

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

Положения, выносимые на защиту:

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

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

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

Научная новизна исследования:

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

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

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

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

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

Практическая значимость:

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

- способов анализа сейсмических записей при проведении традиционных инженерно-сейсмических исследований МПВ;

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

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

Личный вклад автора заключается в разработке и внедрении:

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

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

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

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

Апробация работы

Основные положения и результаты диссертационного исследования докладывались и обсуждались на следующих конференциях:

- Международной конференции "Экологическая геофизика и геохимия", Москва-Дубна, 1998 г.

- Международном симпозиуме "Карстоведение — XXI век: Теоретическое и практическое значение", Пермь, 25-30 мая 2004 г.

- Научно-практической конференции "Проблемы инженерных изысканий для строительства в Уральском регионе", г. Екатеринбург, 8 июля 2004 г.

- Международной научно-практической конференции "Инжгео - 2008", Геленджик, 25-30 апреля 2008 г.

Публикации

По результатам выполненных исследований опубликовано 7 работ, отражающих основные результаты исследований, в том числе статья, в журнале, рекомендованном ВАК - Крылаткова Н.А. Результаты инженерно-сейсмических исследований карста на Свердловской железной дороге / Н.А. Крылаткова // Горный журнал. - 2008. - № 8. С. 190192.

Структура, содержание и объем работы. Диссертация состоит из введения, трех глав, заключения, списка литературы (120 наименований). Основное содержание диссертационной работы изложено на 102 страницах машинописного текста, иллюстрированного 8 таблицами и 45 рисунками.

Заключение диссертации по теме «Геофизика, геофизические методы поисков полезных ископаемых», Крылаткова, Надежда Анатольевна

Выводы

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

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

Особенность атрибутов МПВ в отличие от атрибутов МОГТ заключается в том, что они получаются по сейсмограмме, а не по временному разрезу. Поэтому при практическом использовании следует отдать предпочтение более надежно определяемым оконным или комплексным атрибутам.

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

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

Получение некоторых атрибутов возможно на базе выполнения инверсии сейсмических записей.

3. ПРИМЕРЫ ПРИМЕНЕНИЯ АТРИБУТНЫХ ГЕОСЕЙСМИЧЕСКИХ РАЗРЕЗОВ ДЛЯ ИЗУЧЕНИЯ ОСЛАБЛЕННЫХ ЗОН В ГЕОЛОГИЧЕСКИХ СРЕДАХ

Возможности применения атрибутов записей МПВ изучались на нескольких участках и площадках инженерно-геологических исследований в течение ряда лет. Это, прежде всего, работы 2001-2008 гг. на участках Свердловской железной дороги: Богданович - Пышминская (1914 -1915 км), Кунгур Пермь (1516-1520 км), Свердловск-Каменск-Уральский (94-95 км), Богданович-Каменск-Уральский (90-91, 244-245, 258-259 км) с целью получения информации о положении карстоопасных зон. В период 2007-2008 гг. проводились исследования территорий под гражданское и промышленное строительство в г. Березовский (на 3 объектах), г. Ноябрьском (ХМАО), а также исследования устойчивости фундаментов эксплуатирующихся жилых домов в г. Каменск-Уральском (3 объекта). Перед сейсморазведочными работами стояли задачи определения участков пониженной прочности связанных с подземными пустотами и зонами трещиноватости.

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

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

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

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

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

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

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

3.1. Выявление закарстованных зон

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

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

Работы па каретоопасных участках Свердловской железной дороге

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

Участок 1913-1914 км перегона Богданович - Пышминская

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

Провалы, связанные с карстом наблюдались на этом участке с 70-х годов. Для защиты железнодорожного полотна и уменьшения нагрузки на грунтовое основание пути укреплены челноками длиной от 10,0 м до 200,0 м. Сами железнодорожные пути расположены на насыпи высотой 0,7-1,0 м. Карстовые формы наиболее ярко проявляются к югу от железнодорожных путей под газопроводом ПК 1914.0 — 1913.9. На участке ПК 19140+00 были обнаружены и затампонированы соединенные понором провалы, размеры которых составляют 9,0x7,0x7,0 м и 5,0x4,0x4,0 м.

К началу работы в 2001 году имелось следующее представление о геологическом строении изучаемой территории. Карстоопасный участок в районе железнодорожной станции Богданович имеет размеры 700 х 400 м и находится в пределах развития карбонатного карста. Верхняя часть геологического разреза сложена известняками нижнего карбона Ci и терригенными отложениями верхнего карбона Сз, имеющими западное простирание. По данным бурения, проводившегося летом 2001 года силами НПЦ «Карст», в западной части участка геологический разрез сверху до глубины 4.0-11.6 метров представлен рыхлыми грунтами, в основном суглинками, реже песками, супесями, а также дресвяными и щебенистыми грунтами. По результатам инженерно-геологического изучения (отчет Дубейковского С.Г. и Афанасиади Э.И.) граница сильно закарстованной зоны на участке исследований располагается в районе ПК 19138+50 - ПК 19140+60. Кровля известняков в середине участка располагается на глубине от 2 до 4 метров, западнее и восточнее она погружается до глубин 19-25 метров.

Толща известняков сложена круто (53-70°) падающими пластами различной прочности. Известняки средне- и сильнотрещиноватые, в верхней части выветрелые до песка, супеси и рухляков. Наиболее сильно разрушены пласты известняковых брекчий, менее разрушены органогенные детритовые известняки [20]. Мощность зоны выветривания 2-3 метра, по трещинам до 20 метров. Трещины круто падающие (до 79-805)— секут породу через 3-10 метров. Общая закарстованность участка около 10 %. Полости примерно в 40 % случаев заполнены глиной или суглинками с примесью дресвы и щебня, в 45 % - глиной текучей и мягкопластичной консистенции; в 15 % - полости не заполнены или заполнены водой. Незаполненные полости приурочены, как правило, к уровню грунтовых вод и располагаются в верхних частях древних заполненных карстовых полостей. Размер полостей колеблется от десятков сантиметров до нескольких метров. Размеры известных провалов на данной территории составляют от 2 до 7 метров, при глубине от 1.5 до 9 метров.

Участок 1516-1517 км железнодорожного перегона Пермь - Кунгур

Объект в значительной степени отличается от предыдущего литологией карстующихся пород. В геоморфологическом отношении участок представляет собой обширную платообразную поверхность, развитую на известняках артинского яруса, в которую врезаны системы каньонообразных долин. Характерная черта рельефа - широкое развитие карстовых форм. Участок находится у тылового шва первой надпойменной террасы р. Бабки. В коренных породах под аллювием, мощность которого составляет от 10 до 18 метров, находится зона интенсивной трещиноватости и закарстованности. Основное направление трещиноватости вертикальное и горизонтальное, заметно разбиение на блоки. Коренные породы пермского возраста представлепы мелко- и среднезернистыми песчаниками с карбонатным цементом, а также гипс-ангидритовыми и карбонатными пачками. Ниже по разрезу с глубины 20-25 метров залегают слабо трещиноватые доломиты. Характерная черта карстующихся пород - загипсованность. Она определяет интенсивность развития карста, строение полостей, наличие зон разуплотнения, в которых карбонатная мука находится в текучем состоянии. Углы падения пластов небольшие.

Карстопроявления на участке связаны с тремя горизонтами растворимых горных пород. Горизонт разрушенных карбонатных пород обладает благоприятными для образования карста условиями. За счет растворения кровли сульфатных пород образуются зоны разуплотнения, в которых карбонатная мука находится в текучем состоянии. Кроме того, при растворении гипсов развиваются полости в карманах под козырьками сульфатных пород. Гипсовый горизонт отличается интенсивным развитием карста. Выявлены карстовые полости значительной высоты (1.5 метра) и мощные разрушенные зоны. В некоторых случаях встречается этажное строение полостей, разделенных незначительными прослоями гипсов. Мощность разрушенных зон меняется от 0.1 до 1.0 метра. Рухляк представлен щебнем, дресвой гипсов. Ангидритовый горизонт характеризуется незначительным развитием карста. Ангидриты в результате процессов гидратации и выщелачивания деформированы и разрушены, что проявляется в увеличении степени их неоднородности. Обнаружены маломощные разрушенные зоны от 0.1-0.3 метра, представленные щебнем, пластинами ангидрита.

Методика полевых работ

На обоих участках основной объем полевых сейсмических работ выполнялся методом преломленных волн (МПВ). Основным видом наблюдений являлось продольное сейсмическое профилирование с получением систем встречных и нагоняющих годографов. Профили располагались в основании насыпи железнодорожных путей параллельно железной дороге. На каждой стоянке сейсмоприемников записывались сейсмограммы от нескольких пупктов возбуждения. При регистрации применялось вертикальное направление возбуждения колебаний и вертикальная ориентировка сейсмоприемников. Шаг сейсмоприемников составлял 1 или 2 м, а шаг пунктов возбуждения равнялся 10, 20 или 46 метров при работах МПВ и 2-4 метра при работах MOB. Длина годографа 46 м.

Для наблюдений были использованы сейсморазведочные станции с цифровой записью ("Урал-Мини", разработчик Сенин JI.H., Институт геофизики УРО РАН, и "Эхо-3"), обладающие возможностью накопления слабых сигналов. Это позволило обеспечить необходимую для решения поставленных задач глубину исследований при малой энергии источника - 15-20 м.

Интерпретация данных

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

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

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

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

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

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

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

Рнс.3.1. Годографы первых вступлений по профилям а-01.01 (к северу от полотна дороги) и 02.01 (к югу) по участку 1913-1914 км, ст. Богданович-Пышминская (по материалам Крылаткова С.М. Крьшатковой Н.А., 2002 г)

Третий шаг заключался в привязке скоростных построений к геологическому разрезу и выяснению характерных значений скорости для грунтов различного состава и состояния. Результатом такой операции явилась классификация скоростных разрезов по величине скорости н ее вертикального градиента в различных условиях. При этом анализировались отсутствие наличие полостей, глубина их залегания, положение уровня грунтовых вод, мощность и состав четвертичных отложений. Так для участка ст. Богданович в интервале глубин от 0 до 5 метров, соответствующих рыхлым отложениям, можно было выделить два типа разреза. Первый - нязкоскоростной (от 200 до 500 м/с) с линейной зависимостью скорости от глубины и малым градиентом скорости, второй -высокоскоростной (от 200 до 1000-1500 м/с) с экспоненциальной зависимостью и с высоким градиентом скорости. Первый тип скоростного разреза характеризует участки пониженной скорости, соответствующие разуплотнению грунтов над карстовыми полостями, расположенными в толще известняков (скрытые полости), либо областям трещиноватости. Второй тип скоростного разреза - нормальный разрез с резким изменением скорости вблизи уровня грунтовых вод или скальных грунтов. В обоих случаях нельзя дать однозначную геологическую интерпретацию скоростного разреза. Другой тип неоднозначности обнаружился при сопоставлении с бурением разрезов 2 типа, когда не исключены как возможность резкого перехода в разрезе к обводненным рыхлым грунтам, так и присутствие карстовой полости, заполненной водой вблизи УГВ. Это послужило толчком для последующего привлечения дополнительных данных. Глубже по разрезу в отсутствие карстовых полостей скорость продольных волн плавно увеличивается от 2000 м/с до 3000-3500 м/с, что соответствует постепенному увеличению прочности известняков с глубиной. В области скрытого карста наблюдается локальное уменьшение скорости, а в интервале от глубины залегания УГВ до максимальной глубины исследования значения скоростей в разрезе близки к 2000 м/с. Для разрезов в отсутствие известняков скорость возрастает по линейному закону от 1500-1600 м/с до 2000-2500 м/с.

Пример сопоставления скоростных разрезов с данными бурения по Кунгурскому участку приведен на рис. 3.2.

C47S3g

•в/

-8

-10

7*7" у.

B.S и

42

-10

328

О О а

С4778с

2 ВО

Z7

-6

-10 Ш

0,7 Z2

0 ■6.7

9,е о -2 -4 -в •8 •10

75- о о -----

0 S

0 I с о 1

1000 2000 3000

1000 2000 3000

Рис. 3.2 Пример сопоставления скоростных разрезов с данными бурения по Кунгурскому участку; а -участок с резким изменением скорости вблизи УГВ, б - участок с резким изменением скорости вблизи кровли разрушенных известняков н инверсией скорости, соответствующей разрушенной зоне

Анализ разрезов в соответствии с обучением на материалах бурения позволил классифицировать данные по наличию-отсутствию ослабленных участков. Кроме того, были определены характерные значения скоростей для изученного интервала геологического разреза (табл. 3.1). Сведения о скоростях учитывались при составлении геосейсмических разрезов.

Четвертый шаг - изучение латеральной изменчивости скорости в сопоставлении со скважшшыми данными выполнялся по вертикальным картам скоростей при анализе по профилям, а для участка Богдановича дополнительно по горизонтальным картам скоростей. На обоих объектах по данным вертикальных карт было отмечено, что скважины с полостями попадают на низкоскоростные участки, окаймленные резкими горизонтальными градиентами скорости (рис. 3.4). Для неизмененных участков, а также участков с некарстующимися породами характерен плавный горизонтальный градиент скорости. Горизонтальные карты были составлены для глубин 2,4,10 и 15 м. Глубина 2 м соответствовала подошве четвертичных суглинков, 5м — близка к средней отметке УГВ, 10 м - отражала среднюю глубину расположения менее выветрелых грунтов, 15 м -соответствовала минимальной глубине исследований.

ЗАКЛЮЧЕНИЕ

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

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

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

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

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

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

Список литературы диссертационного исследования кандидат геолого-минералогических наук Крылаткова, Надежда Анатольевна, 2008 год

1. Авербух А.Г., 1975, Интерпретация материалов сейсморазведки преломленными волнами. -М.: Недра, 223 с.

2. Аки К., Ричарде П., 1983, Количественая сейсмология: Теория и методы. Т. 1.- М.:Мир.-520 е., Т.2.-М.:Мир .-360 с.

3. Альтяшева Е.П., Николаев И.В., Потапов О.А., 1991, Использование рефрагированных волн в задачах сейсморазведки. Разведочная геофизика: Обзор / М.: «Геоинформмарк», 68 с.

4. Ампилов Ю.П., 2004. Сейсмическая интерпретация: опыт и проблемы. М.: Геоинформмарк.-286 с.

5. Афанасиади Э.И., Бодин В.В., Грязное О.Н., Дубейковский С.Г., Нещсткин О.Б., 2000, Изучение, оценка и прогноз закарстованности карбонатных массивов железнодорожных трасс. Изв. УГГА, Серия: Геология и геофизика, 10. С.229-233.

6. Берзон И.С., Епинатьева A.M., Парийская Г.Н., Стародубовская С.П., 1962, Динамические характеристики сейсмических волн в реальных средах, М.: Изд-во АН СССР, 511с.

7. Боганик Г.Н., Номоконов В.П., 1994, Выявление карстовых образований высокоразрешающей сейсморазведкой MOB. Геофизика, 2.- С.52-53.

8. Бондарев В.И., 1974. Рекомендации по применению сейсмической разведки для изучения физико-механических свойств рыхлых грунтов в естественном залегании для строительных целей. М.: Стройиздат. — 142 с.

9. Бондарев В.И., 1983, Кинематические особенности годографов первых волн, связанных с полупространством, перекрытым однородным или градиентным слоем // Изв. вузов. Геол. и разведка, 3, 65-74.

10. Бондарев В.И., 1997. Сейсмический метод определения физико-механических свойств нескальных грунтов. Екатеринбург: Изд-во УПТА.- 220 с.

11. Бондарев В.И., Крылатков С.М., 2006, Получение сейсмических изображений геологической среды. Учебное пособие для вузов. Екатеринбург: Изд-во УГГУ 304 с.

12. Бондарев В.И., Крылатков С.М., Пущина Ж.В., Крылаткова Н.А., 1995. Интерпретация годографов первых волн в инженерной сейсморазведке. Екатеринбург: Изд-во УПТА. - 84 с.

13. Бондарик Г.К., Пендин В.В., Ярг Л.А., 2007, Инженерная геодинамика, М.; КДУ. - 440 с.

14. Воронков O.K., 1965, Методика и результаты изучения малых глубин методом КМПВ в условиях многолетней мерзлоты. Геология и геофизика, СО АН СССР, 4. С. 875-883.

15. Воронков O.K., Акатов Ю.Е., 1967, Сейсморазведочные работы на карсте. Геология и геофизика, 6, с.98-105

16. Гальперин Е.И., 1977. Поляризационный метод сейсмических исследований. М.: Недра, 319 с.

17. Горбова С.В., 2006, Карбонатный карст Сухоложско-Каменского района (восточный склон Среднего Урала), Автореферат на соискание ученой степени кандидата геолого-минералогических наук, Екатеринбург, Изд-во УГГУ, 24 с.

18. Горяинов Н.Н., Ляховицкий Ф.М., 1979, Сейсмические методы в инженерной геологии. — М.: Недра. 143 с.

19. Давыдова Е.А., Ворошилова М.С., Самапин М.А., Холмянская Н.А., 2002, Методика выявления малоамплитудных разломов по спектрально- временным параметрам. Геофизика, 4, Москва, С. 7-10.

20. Зинченко B.C., 2005, Петрофизические основы гидрогеологической и инженерно-геологической интерпретации геофизических данных: Учебное пособие для студентов вузов. -М.-Тверь: Изд-во АИС. — 392 с.

21. Инженерно-геологические изыскания: Справочное пособие, 1989 / Н.Ф. Архипов, Е.С. Карпышев, J1.A. Молоков, В.А. Парфинович. -М: Недра-288 с.

22. Инженерно-геологические изыскания: Справочное пособие/ Н.Ф Арипов, Е.С. Карпышев, JT.A. Молоков, В.А. Парфиянович. — М.: Недра, 1989. — 288 е.: ил.

23. Интерпретация данных сейсморазведки: Справочник, 1990 / Под редакцией О.А. Потапова. -М.: Недра.- 448 с.

24. Каляшин С.В., 2008, Особенности спектрального состава поверхностной сейсмической волны при возбуждении и распространении. Геофизика, 4. - С. 22-27.

25. Канарейкин Б.А., Прихода А.Г., Сагайдачная О.М., Щербаков В.В., 2004 Сейсмотомографические исследования грунтов при инженерно-геологических изысканиях: Геофизика, 5, С. 27-30.

26. Карасик В.М. , 1993. Изучение скоростей сейсмических волн комплексом методов. — М.: Недра- 220 с.

27. Кашик А.С., 1998, Изучение многомерных многопараметровых пространств на ЭВМ. Их формирование и представление методами динамической визуализации. (Философия и идеология): Геофизика, 1, 84-95 с.

28. Кожевников В. Н., 1985 Изучение и оценка гидрогеологических условий в районах развития современных карстово-суффозионных процессов. М.: ЦП Hi 10. -165с.

29. Комплексные инженерно-геофизические исследования при строительстве гидротехнических сооружений, 1990 / А.И. Савич, Б.Д. Куюнджич, В.И. Коптев и др.; Под ред. А.И. Савича, Б.Д. Куюнджича. -М.: Недра. — 462 с.

30. Крылаткова Н.А., 1998 Изучение двумерно-неоднородных сред по годографам первых волн / В.И. Бондарев, Р.В. Кузнецов, С.М. Крылатков, Н.А. Крылаткова // Известия УГГГА. -Вып. 8. Серия: Геология и геофизика. С. 137-142.

31. Крылаткова Н.А. 2003 Исследование карстоопасных участков железных дорог с помощью сейсморазведки /С.М. Крылатков, Н.А. Крылаткова, О.Б. Нещеткин // Известия УГГГА. -Вып. 18. Серия: Геология и геофизика. С. 137-142.

32. Крылаткова Н.А., 2008 Результаты инженерно-сейсмических исследований карста на Свердловской железной дороге // Известия вузов. Горный журнал. № 8 .- С.190-192.

33. Кутепов В.М., Кожевникова В.Н., 1989, Устойчивость закарстованных территорий, М.: Наука. — 151 с.

34. Лаврова Л.Д. Интерпретация данных наземного сейсмического профилирования в одной из моделей неоднородного скального основания. Сборник научных трудов Гидропроекта, М.: Гидропроект, Вып.89. С.45-51

35. Леоненко М.В., 2004, Опыт использования геофизических методов при оценке карстоопасности в районах покрытого карста. // Материалы международного симпозиума "Карстоведение XXI век: теоретическое и практическое значение". - Пермь. С. 228-230

36. Лисин В.П., Щаранский А.И., Коган И.И. , 2006 Геофизический мониторинг состояния грунтов берегового склона нефтяного терминала// Тезисы докладов Второй международной научно-практической конференции «Инженерная геофизика 2006». — Геленджик. — С. 38-40.

37. Ломтадзе В.Д., 1977, Инженерная геология, Л.: Недра. 479 с.

38. Ляховицкий Ф.М, Хмелевской В.К., Ященко З.Г., 1989, Инженерная геофизика. — М.: Недра, 252 с.

39. Метод преломленных волн, 1990/ A.M. Епинатьева, Г.М. Голошубин, А.Л. Литвин и др.: Под ред. A.M. Епинатьевой. М.: Недра. - 297 е.: ил.

40. Михальцев А.В., Мушин И.А., Погожев В.М., 1990. Обработка динамических параметров в сейсморазведке. М.: Недра.- 189 с.

41. Назарный С. А., 1974, Интерпретация данных малоглубинной сейсморазведки с помощью полей скоростей. Л.: Недра. 81 с.

42. Никитин А.А., Хмелевской В.К. 2004, Комплексирование геофизических методов: учебник для вузов. — Тверь: ООО « Издательство ГЕРС»,. 294 с

43. Никитин В.Н., 1981. Основы инженерной сейсмики. М.: Изд-во МГУ. - 177 с.

44. Огильви А.А. Основы инженерной геофизики: Учеб. для вузов/ Под ред. В.А. Богословского. -М.: Недра, 1990.-501 е.: ил.

45. Палагин В.В., Попов А.Я., Дик П.И., 1989, Сейсморазведка малых глубин. -М.: Недра. 210 с.

46. Пийп В.Б. Ефимова Е.А. 1985, Определение скоростных полей по материалам инженерной сейсморазведки. Вестник МГУ. Сер. Геология, 3, С. 51-56

47. Подьяпольский Г.С., Васильев Ю.И., 1960, О волне релеевского типа на несвободной поверхности. Известия академии наук СССР, серия геофизическая, 9, стр 1289-1308

48. Поперечные и обменные волны в сейсморазведке, 1967/ Под ред. Н.Н. Пузырева. — М.: Недра. -288 с.

49. Применение сейсмоакустических методов в гидрогеологии и инженерной геологии, 1992.:Мин-во геологии СССР, Всесоюз. науч.- исслед. ин-т гидрогеол. и инж. геол.- Под ред. Н.Н. Горяинова. М.: Недра. 264 с.

50. Птецов С.Н., 1989. Анализ волновых полей для прогнозирования геологического разреза. -М.: Недра. 135с

51. Пузырев Н.Н., 1997. Методы и объекты сейсмических исследований. Новосибирск: Издательство СО РАН. - 301 с.

52. Пузырев Н.Н., Тригубов А.В., Бродов Л.Ю. и др., 1985. Сейсмическая разведка методом поперечных и обменных волн. — М.: Недра. 277 с.

53. Рекомендации по изучению карста геофизическими методами 1986 / ПНИИИС — М.: Стройиздат. 112 с.

54. Ризниченко Ю.В., 1985, Сейсморазведка слоистых сред. М.: Недра, 184 с.

55. Рослов Ю.В., Ефимова Н.Н. , 2007, Построение скоростных моделей среды в системе сейсмической томографии ХТОМО//Тезисы докладов Третьей международной научно-практической конференции «Инженерная и рудная геофизика 2007». —Геленджик. — С. 198-200.

56. Саваренский И.А., Миронов Н.А. , 1995 Руководство по инженерно-геологическим изысканиям в районах развития карста. М.: ПНИИС Минстроя России. -165 с.

57. Савич А.И., Коптев В.И., Никитин В.Н., Ященко З.Г. , 1969 Сейсмоакустические методы изучения массивов скальных пород. — М.: Недра. 240 с.

58. Савич А.И., Куинджич Б.Д., Коптев В.И. и др., 1990. Комплексные инженерно — геофизические исследования при строительстве гидротехнических сооружений. М.: Недра. -449 с.

59. Савич А.И., Ященко З.Г. Исследование упругих и деформационных свойств горных пород сейсмоакустическими методами. М.: Недра, 1979. - 214 с.

60. Савич А.И., Ященко З.Г., 1979. Исследование упругих и деформационных свойств горных пород сейсмоакустическими методами. -М.: Недра. 214 с.

61. Санфиров И.А., Ярославцев А.Г., 2004, Опыт применения сейсморазведки для решения инженерно-геологических задач ОГТ. Геофизика, 4, С. 27-30.

62. Сейсмическая разведка методом поперечеых и обменных волн/ Пузырев Н.Н., Тригубов А.В., Бродов Л.Ю. и др. М.: Недра, 1985. - 277 с.

63. Сейсмическая томография. С приложениями в глобальной сейсмологии и разведочной геофизике: Пер. с англ., 1990 / Под ред. Г. Нолета.-М.: Мир.— 416 с.

64. Спасский Б.А., 1986. Основы цифровой обработки данных сейсморазведки. Пермь: Изд-во ПГУ.-96 с.

65. Справочник геофизика. Т.4. Сейсморазведка. М.: Недра, 1981.

66. Старобинец А.Е. Выделение и интерпретация дифрагированных и квазидифрагированных волн. -М.: Недра, 1988. 199 е.: ил.

67. Телегин А.Н., 2004. Сейсморазведка методом преломленных волн. — СПб.: Издательство С.-Петербургского университета, 187 с.

68. Толмачев В.В., Ройтер Ф., 1990, Инженерное карстоведение, М.: Наука. - 151 с.

69. Федотов С.А., 2006 Высокоразрешающие технологии сейсморазведки и электроразведки для инженерных исследований// Тезисы докладов Второй международной научно-практической конференции «Инженерная геофизика 2006». — Геленджик. — С. 36-38.

70. Хоменко В.П., 1986, Карстово-суффозионныс процессы и их прогноз М.: Наука, 97 с.

71. Хромова И.Ю.:, 2007 Технология построения цифровой сейсмогеологической модели на примере программного комплекса Landmark. М:. 1 часть — 315 с.

72. Циммер В.А., 1991, Метод рефрагированных волн при исследовании геологических сред. Алма-Ата: Изд-во КазВИРГ, 87 с.

73. Шерифф Р., Грегори А., 1982, Сейсмическая стратиграфия.- М.: Мир, Т.1.-373 е., Т.2.-460 с.

74. Ямщиков B.C., 1984, Волновые процессы в массиве горных пород.: Учебник для вузов. — М.: Недра, 271 с.

75. Beaty K.S., Schmitt D.R., Sacchi М., 2002, Simulated annealing inversion of multimode Rayleigh wave dispersion curves for geological structure: Geophys. June 2002, 151, pp 622-631.

76. Bodet L., Van Wijk K., Bitri A., Abraham O., Cote P., 2005, Surface-wave inversion limitations from laser-doppler physical modeling: Journal of Environmental and Engineering Geophysics, June 2005, Volume 10, Issue 2, pp 151-162.

77. Brown A R., 2004. Interpretation of Three-Dimensional Seismic Data.- Tulsa, AAPG Memoir 42.542 p.

78. Chen C., Liu J., Xia J., Li Zh., 2006, Integrated geophysical techniques in detecting hidden dangers in river embankments: Journal of Environmental and Engineering Geophysics, June 2006, Volume 11, Issue 2, pp 83-94.

79. De Bremaecker J.C., 1958, Transmission and reflection of Rayleigh waves at corners: Geophysics, April 1958, Volume 23, Issue 2, pp 253-266.

80. De Nil D., 2005, Characteristics of surface waves in media with significant vertical variations in elasto-dynamic properties: Journal of Environmental and Engineering Geophysics, September 2005, Volume 10, Issue 3, pp 263-274.

81. Ewing W.M., Jardetzky W.S., Press F., 1957, Elastic waves in layered media: International ed., McGraw-Hill Book Company inc.

82. Gelis C., Leparoux D., Virieux J., Bitri A., Operto S., Grandjean G., 2005, Numerical modeling of surface waves over shallow cavities: Journal of Environmental and Engineering Gcophysics, June 2005, Volume 10, Issue 2, pp 111-121.

83. Gu H., Cai Ch., Wang Y., 2006, Investigation of fractures using seismic computerized crossholc tomography: Journal of Environmental and Engineering Geophysics, June 2006, Volume 11, Issue 2, pp 143-150.

84. Ivanov J., Miller R.D., Xia J., Steeples D., Park C.B., 2006, Joint analysis of refractions withsurface waves: an inverse solution to the refraction-traveltime problem: Geophysics, November-December 2006, Volume 71, Issue 6, pp 131 -13 8.

85. Ivanov J., Park C.B., Miller R.D., Xia J., 2005, Analyzing and filtering surface-wave energy by muting shot gathers: Journal of Environmental and Engineering Geophysics, September 2005, Volume 10, Issue 3, pp 307-322.

86. Lai C.G., Foti S., Rix G.J., 2005, Propagation of data uncertainty in surface wave inversion: Journal of Environmental and Engineering Geophysics, June 2005, Volume 10, Issue 2, pp 219-228.

87. Miller C.R., Allen A.L., Speece M.A., El-Werr A-K., Link C.A., 2005, Land streamer aided geophysical studies at Saqqara, Egypt: Journal of Environmental and Engineering Geophysics, December 2005, Volume 10, Issue 4, pp 371-380.

88. Neducza В., 2007, Stacking of surface waves: Geophysics, March-April 2007, Volume 72, Issue 2, pp V51-V58.

89. O'Neill A., Matsuoka Т., 2005, Dominant higher surface-wave modes and possible inversion pitfalls: Journal of Environmental and Engineering Geophysics, June 2005, Volume 10, Issue 2, pp

90. Park C.B., Miller R.D., Xia J., 1999, Multichannel analysis of surface waves: Geophysics, May-June 1999, Volume 64, Issue 3, pp 800-808.

91. Reservoir Geophysics, 1992,- Investigation in Geophysics, № 7. Edited by Robert E Sheriff. -Tulsa, Society of Exploration Geophysicists. 400 p.

92. Rydcn N., Park C.B., 2006, Fast simulated annealing inversion of surface waves on pavement using phase-velocity spectra: Geophysics, July-August 2006, Volume 71, Issue 4, pp R49-R58.

93. Sheehan J.R., Doll W.E., Mandell W.A., 2005, An evaluation of methods and available software for seismic refraction tomography analysis: Journal of Environmental and Engineering Geophysics, March 2005, Volume 10, Issue 1, pp 21-34.

94. Steeples D.W., Schmeissner C.M., Macy B.K., 1995, The evolution of shallow seismic exploration methods: Journal of Environmental and Engineering Geophysics, July 1995, Volume 0, Issue 1, pp 15

95. Taner M.T., Koehler F., Sheriff R.E., 1979. Complex seismic trace analysis: Geophysics, Volume 44, Issue 6, 1041-1063 pp.

96. Van Ovcrmeeren R.A., 2001. Hagedoorn's plus-minus method: the beauty of simplicity. Geophysical Prospecting, 2001, 49, 687-696.

97. Wisen R., Christiansen A.V., 2005, Laterally and mutually constrained inversion of surface wave seismic data and resistivity data: Journal of Environmental and Engineering Geophysics, September 2005, Volume 10, Issue 3, pp 251-262.

98. Xia J., Chen C., Tian G., Miller R.D., Ivanov J., 2005, Resolution of high-frequency Rayleigh-waved data: Journal of Environmental and Engineering Geophysics, June 2005, Volume 10, Issue 2, pp 99-110.

99. Xia J., Miller R.D.,Park C.B., 1999, Estimation of near-surface shear-wave velocity by inversion of Rayleigh waves: Geophysics, May-June 1999, Volume 64, Issue 3, pp 691-700.

100. Yamanaka H., 2005, Comparison of performance of heuristic search methods for phase velocity inversion in shallow surface wave method: Journal of Environmental and Engineering Geophysics, June 2005, Volume 10, Issue 2, pp 163-173.185.201.24.